Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized sampling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

Phase evolution during conventional and reactive flash sintering of (Mg,Ni,Co,Cu,Zn)O via in situ X-ray diffraction

Reactive flash sintering (RFS) enables the simultaneous synthesis and sintering of ceramics and has been shown to affect the reaction pathway of different materials. Herein, in situ synchrotron X-ray diffraction (XRD) is used to investigate the (Mg,Ni,Co,Cu,Zn)O entropy-stabilized oxide formation during: (i) conventional heating and (ii) RFS under current rate-controlled mode. The same reaction pathway is verified in both instances: the starting rock-salt (RS), spinel (Co 3 O 4 ), tenorite (CuO), and wurtzite (ZnO) phases transform into a single RS phase with a (1 1 1) to (2 0 0) intensity ratio of 0.67, consistent with a random distribution of the cations into the structure. Pt lattice peak shift from the XRD patterns is used as standard to monitor the sample surface temperature, revealing a strong endothermic reaction during the RS single-phase formation (Pt peaks shift toward higher angles while increasing sample temperature/current density). Furthermore, in RFS, the single-phase RS structure is formed in just 60 s at a furnace temperature of 600°C and a current rate of 220 mA mm -2 /min. Therefore, RFS greatly accelerates the synthesis of (Mg,Ni,Co,Cu,Zn)O, however, it does not play a role in the reaction pathway for this material formation.

36 MATERIALS SCIENCE↗

SATYA: Defending Against Adversarial Attacks Using Statistical HypothesisTesting

The paper presents a new defense against adversarial attacks for deep neural networks. We demonstrate the effectiveness of our approach against the popular adversarial image generation method DeepFool. Our approach uses Wald’s Sequential Probability Ratio Test to sufficiently sample a carefully chosen neighborhood around an input image to determine the correct label of the image. On a benchmark of 50,000 randomly chosen adversarial images generated by DeepFool we demonstrate that our method SATYA is able to recover the correct labels for 95.76% of the images for CaffeNet and 97.43% of the correct label for GoogLeNet.

Raj, Sunny↗

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗

Density Dependence of the Phases of the v = 1 Integer Quantum Hall Plateau in Low Disorder Electron Gases

Recent magnetotransport measurements in low-disorder electron systems confined to GaAs/AlGaAs samples reveal that the v = 1 integer quantum Hall plateau is broken into three distinct regions. These three regions are associated with two phases with different types of bulk localization: the Anderson insulator is due to random quasiparticle localization, and the integer quantum Hall Wigner solid is due to pinning of a stiff quasiparticle lattice. Universal properties of the v = 1 plateau are highlighted: the structure of the stability diagram, the nonmonotonic dependence of the activation energy on the filling factor, and the alignment of features of the activation energy with features of the stability regions of the different phases are found to be similar in three samples spanning a wide range of electron densities. Quantitative differences between the samples are also discussed, such as the dependence of the onset temperature and the activation energy of the integer quantum Hall Wigner solid on the electron density. The findings provide insights into the localization behavior along the v = 1 integer quantum Hall plateau in the low disorder regime.

Anderson insulator↗

Estimation of hydraulic conductivity in a watershed using sparse multi-source data via Gaussian process regression and Bayesian experimental design

Enhanced water management systems depend on accurate estimation of subsurface hydraulic properties. However, geologic formations can vary significantly, so information from a single source (e.g., widely spaced boreholes) is insufficient in characterizing subsurface aquifer properties. Therefore, multiple sources of information are needed to complement the hydrogeology understanding of a region. Here, this study presents a numerical framework in which information from different measurement sources is combined to characterize the 3D random field in a multi-fidelity prediction model. Coupled with the model, a Bayesian experimental design was used to determine the best future sampling locations. The Upper Sangamon watershed in east-central Illinois was selected as the case study site, where the multi-fidelity Gaussian process model was used to estimate the hydraulic conductivity in the region of interest. Multi-source observation data were obtained from electrical resistivity and borehole pumping tests. The accuracy of the model prediction is dependent on the locations and the distribution of both high- and low-fidelity data. Furthermore, the multi-fidelity model was compared with the single-fidelity model. The uncertainties and confidence in the measurements and parameter estimates were quantified and used to design future cycles of data collection to further improve the confidence intervals.

54 ENVIRONMENTAL SCIENCES↗

Strontium isoscape of sub-Saharan Africa allows tracing origins of victims of the transatlantic slave trade

Abstract Strontium isotope ( 87 Sr/ 86 Sr) analysis with reference to strontium isotope landscapes (Sr isoscapes) allows reconstructing mobility and migration in archaeology, ecology, and forensics. However, despite the vast potential of research involving 87 Sr/ 86 Sr analysis particularly in Africa, Sr isoscapes remain unavailable for the largest parts of the continent. Here, we measure the 87 Sr/ 86 Sr ratios in 778 environmental samples from 24 African countries and combine this data with published data to model a bioavailable Sr isoscape for sub-Saharan Africa using random forest regression. We demonstrate the efficacy of this Sr isoscape, in combination with other lines of evidence, to trace the African roots of individuals from historic slavery contexts, particularly those with highly radiogenic 87 Sr/ 86 Sr ratios uncommon in the African Diaspora. Our study provides an extensive African 87 Sr/ 86 Sr dataset which includes scientifically marginalized regions of Africa, with significant implications for the archaeology of the transatlantic slave trade, wildlife ecology, conservation, and forensics.

Science & Technology - Other Topics↗

The ASAS-SN catalogue of variable stars – IV. Periodic variables in the APOGEE survey

ABSTRACT We explore the synergy between photometric and spectroscopic surveys by searching for periodic variable stars among the targets observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 1924 periodic variables among more than $258\, 000$ APOGEE targets; 465 are new discoveries. We homogeneously classified 430 eclipsing and ellipsoidal binaries, 139 classical pulsators (Cepheids, RR Lyrae, and δ Scuti), 719 long-period variables (pulsating red giants), and 636 rotational variables. The search was performed using both visual inspection and machine learning techniques. The light curves were also modelled with the damped random walk stochastic process. We find that the median [Fe/H] of variable objects is lower by 0.3 dex than that of the overall APOGEE sample. Eclipsing binaries and ellipsoidal variables are shifted to a lower median [Fe/H] by 0.2 dex. Eclipsing binaries and rotational variables exhibit significantly broader spectral lines than the rest of the sample. We make ASAS-SN light curves for all the APOGEE stars publicly available and provide parameters for the variable objects.

Pawlak, Michał↗

SU-F-T-192: Study of Robustness Analysis Method of Multiple Field Optimized IMPT Plans for Head & Neck Patients

Purpose: Proton therapy is more sensitive to uncertainties than photon treatments due to protons’ finite range depending on the tissue density. Worst case scenario (WCS) method originally proposed by Lomax has been adopted in our institute for robustness analysis of IMPT plans. This work demonstrates that WCS method is sufficient enough to take into account of the uncertainties which could be encountered during daily clinical treatment. Methods: A fast and approximate dose calculation method is developed to calculate the dose for the IMPT plan under different setup and range uncertainties. Effects of two factors, inversed square factor and range uncertainty, are explored. WCS robustness analysis method was evaluated using this fast dose calculation method. The worst-case dose distribution was generated by shifting isocenter by 3 mm along x,y and z directions and modifying stopping power ratios by ±3.5%. 1000 randomly perturbed cases in proton range and x, yz directions were created and the corresponding dose distributions were calculated using this approximated method. DVH and dosimetric indexes of all 1000 perturbed cases were calculated and compared with the result using worst case scenario method. Results: The distributions of dosimetric indexes of 1000 perturbed cases were generated and compared with the results using worst case scenario. For D95 of CTVs, at least 97% of 1000 perturbed cases show higher values than the one of worst case scenario. For D5 of CTVs, at least 98% of perturbed cases have lower values than worst case scenario. Conclusion: By extensively calculating the dose distributions under random uncertainties, WCS method was verified to be reliable in evaluating the robustness level of MFO IMPT plans of H&N patients. The extensively sampling approach using fast approximated method could be used in evaluating the effects of different factors on the robustness level of IMPT plans in the future.

60 APPLIED LIFE SCIENCES↗

Conservative Estimation of Tail Probabilities from Limited Sample Data

Several sparse-sample uncertainty quantification (UQ) methods are compared for conservative but not overly conservative estimation of small tail probabilities involving responses that lay beyond specified thresholds in the tails of probability distributions. Sixteen very differently shaped distributions (or probability density functions, PDFs) and tail probability magnitudes ranging from 10 -5 to 10 -1 are considered in order for the study to be relevant to a wide range of risk analysis and quantification of margins and uncertainty (QMU) problems. The emphasis of the study is on limited data regimes ranging from N = 2 to 20 samples, reflective of most experimental and some expensive computational situations. Relatively simple sparse-sample UQ methods tested for this regime involve statistical tolerance interval "Equivalent Normal and related "Ensemble of Normals" and "Superdistribution (SD) approaches. (The independently derived SD is effectively equivalent to the Bayesian posterior predictive distribution given the assumptions of the derivation.) The performance of the methods was generally improved for N ≥ 5 samples with a generalized Jackknife resampling technique, which determines a tail probability estimate by averaging estimates from smaller sub-samples. Several quantitative metrics for method conservatism and accuracy of tail probability estimation are used to assess and rank the methods' performance over many random trials for each test PDF and probability magnitude. A variant of Bootstrap resampling was also tried, but did not significantly improve tail probability estimates in most cases. Detailed results are presented from over 100-million tests over the above factors that provide useful granular information on which methods or combination of methods perform best in various areas of the factor space.

97 MATHEMATICS AND COMPUTING↗

Deep Learning for In-Situ Layer Quality Monitoring during Laser-Based Directed Energy Deposition (LB-DED) Additive Manufacturing Process

Defects are a leading issue for the rejection of parts manufactured through the Directed Energy Deposition (DED) Additive Manufacturing (AM) process. In an attempt to illuminate and advance in situ quality monitoring and control of workpieces, we present an innovative data-driven method that synchronously collects sensing data and AM process parameters with a low sampling rate during the DED process. The proposed data-driven technique determines the important influences that individual printing parameters and sensing features have on prediction at the inter-layer qualification to perform feature selection. Three Machine Learning (ML) algorithms including Random Forest (RF), Support Vector Machine (SVM), and Convolutional Neural Network (CNN) are used. During post-production, a threshold is applied to detect low-density occurrences such as porosity sizes and quantities from CT scans that render individual layers acceptable or unacceptable. This information is fed to the ML models for training. Training/testing are completed offline on samples deemed “high-quality” and “low-quality”, utilizing only features recorded from the build process. CNN results show that the classification of acceptable/unacceptable layers can reach between 90% accuracy while training/testing on a “high-quality” sample and dip to 65% accuracy when trained/tested on “low-quality”/“high-quality” (respectively), indicating over-fitting but showing CNN as a promising inter-layer classifier.

36 MATERIALS SCIENCE↗

Random field optimization

Herein we present a new modeling paradigm for optimization that we call random field optimization. Random fields are a powerful modeling abstraction that aims to capture the behavior of random variables that live on infinite-dimensional spaces (e.g., space and time) such as stochastic processes (e.g., time series, Gaussian processes, and Markov processes), random matrices, and random spatial fields. This paradigm involves sophisticated mathematical objects (e.g., stochastic differential equations and space-time kernel functions) and has been widely used in neuroscience, geoscience, physics, civil engineering, and computer graphics. Despite of this, however, random fields have seen limited use in optimization; specifically, existing optimization paradigms that involve uncertainty (e.g., stochastic programming and robust optimization) mostly focus on the use of finite random variables. This trend is rapidly changing with the advent of statistical optimization (e.g., Bayesian optimization) and multi-scale optimization (e.g., integration of molecular sciences and process engineering). Our work extends a recently-proposed abstraction for infinite-dimensional optimization problems by capturing more general uncertainty representations. Moreover, we discuss solution paradigms for this new class of problems based on finite transformations and sampling, and identify open questions and challenges.

97 MATHEMATICS AND COMPUTING↗

Excitation spectrum and low-temperature magnetism in the disordered defect-fluorite Ho2⁢Zr2⁢O7

In this work, we report on the thermomagnetic characterization and crystalline-electric field (CEF) energy scheme of the disordered defect-fluorite Ho2⁢Zr2⁢O7. This structural phase is distinguished by the coexistence of magnetic frustration and extensive disorder, with Ho3+ and Zr4+ sharing randomly the same 4⁢𝑎 site with even 50% occupancy, and an average 1/8 oxygen vacancy per unit cell. AC magnetic susceptibility measurements performed on powder samples down to 0.5 K revealed signs of slowing spin dynamics without glassy behavior, including a frequency dependent peak at ∼1K. Yet, no evidence for long-range magnetic order is found down to 200 mK in the specific heat. Inelastic neutron scattering measurements show a weak, low-lying CEF excitation around 2 meV, accompanied by a broad level centered at 60 meV. To fit our observations, we propose an approach to account for structural disorder in the crystal-field splitting of the non-Kramers Ho3+. Our model provides an explanation to the broadening of the high-energy, single-ion excitations and suggests that the zirconate ground-state wavefunction has zero magnetic moment. Through the breaking of local symmetry, structural disorder enables the magnetic response observed in Ho2⁢Zr2⁢O7, allowing the mixing of low-lying states at finite temperatures. Finally, we show that this scenario is in good agreement with the bulk properties reported in this work.

Gardner, Jason [ORNL] (ORCID:0000000278234072)↗

A software package for modeling and simulating fault graphs

This report describes a novel fault graph modeling language and a simulation tool for executing models specified in the language. The modeling language has three primary features that distinguish it from similar reliability analysis tools. These are (1) a random variable modeling several distinct outcomes of a single fault; (2) chains of faults in which one fault triggers another; and (3) time to fail sampled from probability distributions including positive normal, exponential, Weibull with a minimum, or immediate. These features are motivated by their use in a historical analysis of centrifuge reliability.

97 MATHEMATICS AND COMPUTING↗

Genome collection processing for “Conserved upper thermal limits and small safety margins in soil copiotrophic bacteria”

We extracted the genomic DNA of 400 randomly selected isolates using a Quick-DNA Microprep Kit (Zymo Research D3020) according to the manufacturer’s protocol. We then submitted the extracted gDNA samples for short-read Illumina sequencing (200 Mbp) at SeqCoast Genomics (Portsmouth, NH, USA). After preprocessing the sequences using Trimmommatic (Bolger et al. 2014), we assembled the genomes using SPADES (Bankevich et al. 2012) and checked the quality of each assembly using QUAST (Gurevich et al. 2013). We processed the genome assemblies using a KBase (v1.4.0) pipeline (Allen et al. 2017; Arkin et al. 2018). Briefly, we used DRAM (v0.1.2) with default settings to annotate the genome assemblies. We then evaluated genome quality and possible contamination levels using CheckM (v1.0.18) (Parks et al. 2015) and retained genomes with completeness above 98% and contamination below 5% (n = 354), following the authors' guidelines. We then obtained taxonomic assignments for all remaining isolates using the Genome Taxonomy Database tool GTDB-Tk (v2.3.2, database version r214) (Chaumeil et al. 2019). We constructed a phylogenetic tree using the tool SpeciesTree (v2.2.0). We then trimmed the tree (using Trim SpeciesTree to GenomeSet- v1.4.0), retaining only tips within our collection with measured thermal performance.

59 BASIC BIOLOGICAL SCIENCES↗

Sparse sampling methods and probe systems for analytical instruments

Sparse sampling approaches and probe systems for analytical instruments are disclosed providing for effective sub-sampling of a specimen and inpainting to reconstruct representations of actual information. The sub-sampling involves serial acquisition of contiguous measured values lying at positions along a scan path extending in a line toward a first direction and having random perturbations in a second direction. The perturbations are limited within a predetermined distance from the line. Inpainting techniques are utilized among the measured values to reconstruct a representation of actual information regarding the specimen.

Kovarik, Libor↗

Fault Graph (fg)

This software package offers a novel fault graph modeling language and a simulation tool for executing models specified in the language. The modeling language has three primary features that distinguish it from similar reliability analysis tools. These are (1) a random variable modeling several distinct outcomes of a single fault; (2) chains of faults in which one fault triggers another; and (3) time to fail sampled from probability distributions including positive normal, exponential, Weibull with a minimum, or immediate.

Nutaro, James↗

Scheduling strategies for the ESPRESSO follow-up of TESS targets

ABSTRACT Radial-velocity follow-up of stars harbouring transiting planets detected by TESS is expected to require very large amounts of expensive telescope time in the next few years. Therefore, scheduling strategies should be implemented to maximize the amount of information gathered about the target planetary systems. We consider myopic and non-myopic versions of a novel uniform-in-phase scheduler, as well as a random scheduler, and compare these scheduling strategies with respect to the bias, accuracy and precision achieved in recovering the mass and orbital parameters of transiting and non-transiting planets. This comparison is carried out based on realistic simulations of radial-velocity follow-up with ESPRESSO of a sample of 50 TESS target stars, with simulated planetary systems containing at least one transiting planet with a radius below 4R⊕. Radial-velocity data sets were generated under reasonable assumptions about their noise component, including that resulting from stellar activity, and analysed using a fully Bayesian methodology. We find the random scheduler leads to a more biased, less accurate, and less precise, estimation of the mass of the transiting exoplanets. No significant differences are found between the results of the myopic and non-myopic implementations of the uniform-in-phase scheduler. With only about 22 radial velocity measurements per data set, our novel uniform-in-phase scheduler enables an unbiased (at the level of 1 per cent) measurement of the masses of the transiting planets, while keeping the average relative accuracy and precision around 16 per cent and 23 per cent, respectively. The number of non-transiting planets detected is similar for all the scheduling strategies considered, as well as the bias, accuracy and precision with which their masses and orbital parameters are recovered.

Cabona, L.↗

Modeling intrinsic galaxy alignment in the MICE simulation

The intrinsic alignment (IA) of galaxies is potentially a major limitation in deriving cosmological constraints from weak lensing surveys. In order to investigate this effect, we assign intrinsic shapes and orientations to galaxies in the light-cone output of the MICE simulation, spanning ~ 5000 deg 2 and reaching redshift Z = 1.4. This assignment is based on a semianalytic IA model that uses photometric properties of galaxies as well as the spin and shape of their host halos. Advancing on previous work, we include more realistic distributions of galaxy shapes and a luminosity-dependent galaxy-halo alignment. The IA model parameters are calibrated against COSMOS and BOSS LOWZ observations. The null detection of IA in observations of blue galaxies is accounted for by setting random orientations for these objects. We compare the two-point alignment statistics measured in the simulation against predictions from the analytical IA models NLA and TATT over a wide range of scales, redshifts, and luminosities for red and blue galaxies separately. We find that both models fit the measurements well at scales above 8h -1 Mpc, while TATT outperforms NLA at smaller scales. The IA parameters derived from our fits are in broad agreement with various observational constraints from red galaxies. Finally, we build a realistic source sample, mimicking DES Year 3 observations and use it to predict the IA contamination to the observed shear statistics. We find this prediction to be within the measurement uncertainty, which might be a consequence of the random alignment of blue galaxies in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗