Results for “Importance Sampling” — Page 49

NASA NTRS1988

Performance analysis of the Alliant FX/8 multiprocessor using statistical clustering

Results for two distinct, real, scientific workloads executed on an Alliant FX/8 are discussed. A combination of user concurrency and system overhead measurements was taken for both workloads. Preliminary analysis shows that the first sampled workload is comprised of consistently high user concurrency, low system overhead, and little paging. The second sample has much less user concurrency, but significant paging and system overhead. Statistical cluster analysis is used to extract a state transition model to jointly characterize user concurrency and system overhead. A skewness factor is introduced and used to bring out the effects of unbalanced clustering when determining states with important transitions. The results from the models show that during the collection of the first sample, the system was operating in states of high user concurrency approximately 75 percent of the time. The second workload sample shows the system in high user concurrency states only 26 percent of the time. In addition, it is ascertained that high system overhead is usually accompanied by low user concurrency. The analysis also shows a high predictability of system behavior for both workloads.

Dimpsey, Robert Tod↗

DOE OSTI2023

Action Effectiveness Monitoring and Research of Dredged Material Placement at Woodland Islands (Final Report)

Benthic invertebrates play critical ecosystem roles including the breakdown of organic matter, sediment mixing, and nutrient cycling. In the lower Columbia River and estuary, benthic invertebrates provide a foraging resource for threatened and endangered juvenile salmon. Driven by the goal to create low velocity, shallow water, and riparian shrub habitats to benefit juvenile salmon, the USACE placed 237,000 CY of dredged material, resulting in the deposition of 13.5 acres of sand on the off-channel margins at Woodland Islands. The Pacific Northwest National Laboratory (PNNL) designed and implemented an action effectiveness research study to understand how dredged material placement at Woodland Islands affected sediment conditions and benthic invertebrate assemblages. A Before-After-Control-Impact, or BACI, study design was used to evaluate the response of environmental conditions and benthic invertebrates to dredged material placement. The design included one impact site and two control sites. Sampling occurred for two years prior to dredged material placement and two years after. Our study found that spatial variation was a significant factor for both environmental and biological response variables which suggests local conditions are important considerations for mechanisms affecting benthic assemblages. We found that off-channel habitats across all locations sampled were producing benthic invertebrates, many of which are common prey items for juvenile salmon and steelhead—e.g., insects, chironomids, crustaceans, and corophium. There was a significant BACI effect (i.e., an effect at the dredged material placement site, relative to conditions at the control sites) for concentration of carbon and ammonium in sediment, but not for phosphorous concentrations or for percent sand. The estimated abundance for the three invertebrate response variables—total abundance, total chironomid abundance, and total corophium abundance—was significantly lower at the impact site after dredged material placement, compared to the control sites. At Woodland Island, the estimated mean abundance for all invertebrates combined decreased 28% after dredged material placement. Estimated mean abundances of chironomid and corophium decreased by 8% and 88%, respectively. While invertebrate abundances were lower after dredged placement, the composition of benthic invertebrates was similar before and after placement suggesting that as the new habitat feature evolves recolonization will likely follow. These findings provide a foundation for understanding potential benefits and consequences of repurposing dredged material for habitat creation in the LCRE.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1986

Predatory prokaryotes: predation and primary consumption evolved in bacteria

Two kinds of predatory bacteria have been observed and characterized by light and electron microscopy in samples from freshwater sulfurous lakes in northeastern Spain. The first bacterium, named Vampirococcus, is Gram-negative and ovoidal (0.6 micrometer wide). An anaerobic epibiont, it adheres to the surface of phototrophic bacteria (Chromatium spp.) by specific attachment structures and, as it grows and divides by fission, destroys its prey. An important in situ predatory role can be inferred for Vampirococcus from direct counts in natural samples. The second bacterium, named Daptobacter, is a Gram-negative, facultatively anaerobic straight rod (0.5 x 1.5 micrometers) with a single polar flagellum, which collides, penetrates, and grows inside the cytoplasm of its prey (several genera of Chromatiaceae). Considering also the well-known case of Bdellovibrio, a Gram-negative, aerobic curved rod that penetrates and divides in the periplasmic space of many chemotrophic Gram-negative bacteria, there are three types of predatory prokaryotes presently known (epibiotic, cytoplasmic, and periplasmic). Thus, we conclude that antagonistic relationships such as primary consumption, predation, and scavenging had already evolved in microbial ecosystems prior to the appearance of eukaryotes. Furthermore, because they represent methods by which prokaryotes can penetrate other prokaryotes in the absence of phagocytosis, these associations can be considered preadaptation for the origin of intracellular organelles.

Non-NASA Center↗

NASA NTRS2022

The Impact of Sampling Strategy on the Cloud Droplet Number Concentration Estimated From Satellite Data

Cloud droplet number concentration (Nd) is of central importance to observation-based estimates of aerosol indirect effects, being used to quantify both the cloud sensitivity to aerosol and the base state of the cloud. However, the derivation of Nd from satellite data depends on a number of assumptions about the cloud and the accuracy of the retrievals of the cloud properties from which it is derived, making it prone to systematic biases. A number of sampling strategies have been proposed to address these biases by selecting the most accurate Nd retrievals in the satellite data. This work compares the impact of these strategies on the accuracy of the satellite retrieved Nd, using a selection of in situ measurements. In stratocumulus regions, the MODIS Nd retrieval is able to achieve a high precision (r2 of 0.5–0.8). This is lower in other cloud regimes but can be increased by appropriate sampling choices. Although the Nd sampling can have significant effects on the Nd climatology, it produces only a 20 % variation in the implied radiative forcing from aerosol–cloud interactions, with the choice of aerosol proxy driving the overall uncertainty. The results are summarised into recommendations for using MODIS Nd products and appropriate sampling.

Edward Gryspeerdt↗

NASA NTRS2018

Applying Modern Analytical Techniques to the Apollo Samples: A Potential Model for Future Sample Return Missions

From 1969-1972 the Apollo missions collected 382 kg of lunar samples from six distinct locations on the Moon. Studies of the Apollo sample suite have shaped our understanding of the formation and early evolution of the Earth-Moon system, and have had important implications for studies of the other terrestrial planets (e.g., through the calibration of the crater counting record). Despite nearly 50 years of research on Apollo samples, scientists are still developing new theories about the origin and evolution of the Moon. In order to resolve these questions, scientists need access to new lunar samples, particularly new plutonic samples. Although no new large plutonic samples (i.e., hand-samples) remain to be discovered in the Apollo sample collection, there are many large polymict breccias in the Apollo collection containing relatively large (1 cm or larger) previously identified plutonic clasts, as well as a large number of unclassified lithic clasts. In addition, new, previously unidentified plutonic clasts are potentially discoverable within these breccias. The question becomes how to non-destructively locate and identify new lithic clasts of interest while minimizing the contamination and physical degradation of the samples.

Zeigler, Ryan A.↗

NASA NTRS2004

Analyses of IR-Stealthy and Coated Surface Materials: A Comparison of LIBS and Reflectance Spectra and Their Application to Mars Surface Exploration

Identification of non-silicate samples on Mars, such as carbonates, sulfates, nitrates, or evaporites in general, is important because of their association with aqueous processes and their potential as exobiology sites. Infrared (IR) and thermal emission (TE) spectroscopy have been considered the primary tools for remote identification of these minerals. This includes current and future orbital assets such as TES on MGS, THEMIS on Mars Odyssey, OMEGA on Mars Express, CRISM on MRO, and now the Mini-TES on the MER rovers. While reflectance and emission spectroscopy have clearly been the method of choice for these missions, the technique is not always successful in mineral identifications due to dust, surface weathering chemistry, coatings, or surface texture. Here we describe and show IR spectra of several such samples, and then report on the relative success of LIBS analyses in determining the rock type.

Wiens, R. C.↗

DOE OSTI2024

Alpha-Quartz Plastic Strength Investigation Via Diffraction Experiments on Novaculite Using A D-Dia and Elastic Plastic Self- Consistent Interpretation [Thesis}

X-Ray Tomography. Porosity is commonly measured using mercury injection (MI) or water immersion porosimeter (WIP). Both MI and WIP utilize pressure to fill open voids with mercury or water respectively. A measure of the volume change of the sample is then used to estimate the percentage of open voids. However, for the purpose of rheological study it is important to get an accurate measure of open and closed voids. X-ray tomography has been selected as it offers a three-dimensional view into the sample that can quantify all pores limited only by the voxel size which is in the micron range for this study. Radiographs for tomography were collected at the Material Science and Technology Division at Los Alamos National Lab using the Carl Zeiss Xradio 520 instrument and Scout-and-Scan version 16.1 operating software. Two samples were imaged, the starting material and the deformed sample SiO2_65. 3001 radiographs were taken of the starting material with a 6 second exposure time using a 4x objective lens. The x-ray beam was set to 60 kilovoltage peak (kVp) and 5 watts. 1901 radiographs were taken of SiO2_65 with a 25 second exposure time using a 10x objective lens. The x-ray beam was set to 80 kVp and 7 watts. Radiograph files were analyzed by Brian Patterson using Avizo. Void and inclusion volumes were output by voxel sized (1.03 μm) slices, used to calculate a total percentage volume for the starting material.

36 MATERIALS SCIENCE↗

NASA NTRS2011

Principles, Techniques, and Applications of Tissue Microfluidics

The principle of tissue microfluidics and its resultant techniques has been applied to cell analysis. Building microfluidics to suit a particular tissue sample would allow the rapid, reliable, inexpensive, highly parallelized, selective extraction of chosen regions of tissue for purposes of further biochemical analysis. Furthermore, the applicability of the techniques ranges beyond the described pathology application. For example, they would also allow the posing and successful answering of new sets of questions in many areas of fundamental research. The proposed integration of microfluidic techniques and tissue slice samples is called tissue microfluidics because it molds the microfluidic architectures in accordance with each particular structure of each specific tissue sample. Thus, microfluidics can be built around the tissues, following the tissue structure, or alternatively, the microfluidics can be adapted to the specific geometry of particular tissues. By contrast, the traditional approach is that microfluidic devices are structured in accordance with engineering considerations, while the biological components in applied devices are forced to comply with these engineering presets. The proposed principles represent a paradigm shift in microfluidic technology in three important ways: Microfluidic devices are to be directly integrated with, onto, or around tissue samples, in contrast to the conventional method of off-chip sample extraction followed by sample insertion in microfluidic devices. Architectural and operational principles of microfluidic devices are to be subordinated to suit specific tissue structure and needs, in contrast to the conventional method of building devices according to fluidic function alone and without regard to tissue structure. Sample acquisition from tissue is to be performed on-chip and is to be integrated with the diagnostic measurement within the same device, in contrast to the conventional method of off-chip sample prep and subsequent insertion into a diagnostic device. A more advanced form of tissue integration with microfluidics is tissue encapsulation, wherein the sample is completely encapsulated within a microfluidic device, to allow for full surface access. The immediate applications of these approaches lie with diagnostics of tissue slices and biopsy samples e.g. for cancer but the approaches would also be very useful in comparative genomics and other areas of fundamental research involving heterogeneous tissue samples.

Wade, Lawrence A.↗

NASA NTRS1992

Apollo 15 impact melts, the age of Imbrium, and the Earth-Moon impact cataclysm

The early impact history of the lunar surface is of critical importance in understanding the evolution of both the primitive Moon and the Earth, as well as the corresponding populations of planetesimals in Earth-crossing orbits. Two endmember hypotheses call for greatly dissimilar impact dynamics. One is a heavy continuous (declining) bombardment from about 4.5 Ga to 3.85 Ga. The other is that an intense but brief bombardment at about 3.85 +/- Ga was responsible for producing the visible lunar landforms and for the common 3.8-3.9 Ga ages of highland rocks. The Apennine Front, the main topographic ring of the Imbrium Basin, was sampled on the Apollo 15 mission. The Apollo 15 impact melts show a diversity of chemical compositions, indicating their origin in at least several different impact events. The few attempts at dating them have generally not produced convincing ages, despite their importance. Thus, we chose to investigate the ages of melt rock samples from the Apennine Front, because of their stratigraphic importance yet lack of previous age definition.

Ryder, Graham↗

DOE OSTI2023

pH dependence of brown-carbon optical properties in cloud water

Abstract. Light-absorbing organic species present in aerosols, collectively called brown carbon (BrC), have important but highly uncertain effects on climate. Clouds likely represent a significant medium for secondary BrC production and for bleaching reactions, though the relative importance of the formation and loss processes in clouds is unknown at present. The acidity (or pH) of atmospheric particles and clouds affects the optical properties of BrC and bleaching rates. Given the wide variability of pH in the atmosphere (pH in particles and clouds ranges from −1 to 8), the optical properties of BrC and its bleaching behavior are expected to vary significantly, and the link between pH and BrC is yet another uncertainty in attempts to constrain its climate forcing effects. In this work, we characterize the pH dependence of BrC optical properties – including light absorption at 365 nm (Abs365), the mass absorption coefficient (MAC365), and the absorption Ångström exponent (AAE) – in bulk cloud water sampled from the summit of Whiteface Mountain, NY. In all samples (n=17), Abs365 and MAC365 increased linearly with increasing pH, highlighting the importance of reporting pH in studies of BrC in aqueous media. There was strong variability in the sensitivity of Abs365 to pH, with normalized slopes that ranged from 5.1 % to 17.2 % per pH unit. The normalized slope decreased strongly with increasing cloud water [K+], suggesting that the non-biomass-burning BrC has optical properties that are more sensitive to pH than BrC associated with biomass burning. AAE also showed a distinct pH dependence as it was relatively flat between pH 1.5–5 and then decreased significantly above pH 5. The cloud water composition was used to inform thermodynamic predictions of aerosol pH upwind and/or downwind of Whiteface Mountain and the subsequent changes in BrC optical properties. Overall, these results show that, in addition to secondary BrC production, photobleaching, and the altitudinal distribution, the climate forcing of BrC is quite strongly affected by its pH-dependent absorption.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Testing and Characterization to Develop a Mechanistic Explanation for Unsaturated Drift of Fiber Optic Sensors during High-Dose Irradiation

The primary limitation for any optical fiber-based sensor for nuclear reactor applications is radiation-induced attenuation (RIA) of the transmitted and/or reflected signals. Based on several recent studies, RIA is tolerable for some fused silica optical fibers with the proper choice of sensing wavelength and fiber dopants. For extreme temperature applications (> 1000°C), sapphire optical fibers have been proposed; however, recent optical transmission measurements performed on bulk sapphire samples showed prohibitively large RIA. For some sensors, radiation-induced dimensional changes in the fiber materials can also cause significant drift. Moreover, the drift that was observed in numerous experiments performed in the High Flux Isotope Reactor (HFIR), the Advanced Test Reactor, the Massachusetts Institute of Technology Reactor, and other international facilities far exceeded what would be expected based on compaction of fused silica glass. Clearly, additional work is needed to better understand the origins of both RIA and radiation-induced drift in both silica and sapphire optical fiber-based sensors before these sensors can be reliably deployed for nuclear applications. This work evaluated the underlying mechanisms that may be responsible for RIA and drift in silica and sapphire materials. First, detailed characterization was performed on bulk fused silica glass samples that were previously irradiated to different neutron fluences at different temperatures to better understand the structural changes that drive radiation-induced drift in the absence of coating effects that are discussed later. Results show that the non-monotonic compaction that occurs with increasing neutron fluence continues up to fast neutron fluences approaching 10 22 n/cm 2 , which has important implications for physics-based models that may be used to compensate for the sensor drift. Initial Raman spectroscopy and synchrotron x-ray diffraction provide insights into the nature of the structural changes. Next, detailed characterizations were performed on silica fibers with various coatings that were subjected to several different thermal treatments. The hypothesis is that the coatings convert to carbon-rich materials that compact under irradiation, putting a large compressive strain on the fiber. Out-of-pile testing confirms that both polyimide and acrylate fiber coatings convert to glassy carbon (GC) materials when heated under inert conditions, and the degree of order (i.e., graphitization) increases with increasing temperature. The results provide increasingly strong evidence that the combination of polymeric coatings and inert (or vacuum) conditions render fiber optic sensors susceptible to significant radiation-induced drift that would not otherwise exist in uncoated fibers. Finally, transmission electron microscopy was performed on bulk sapphire samples that were irradiated to two neutron fluences at different temperatures to gain insights into the potential mechanisms driving the prohibitive RIA at higher neutron fluences and temperatures. Contrary to previous hypotheses, results show that scattering from radiation-induced voids cannot explain the observed RIA. Similarly, models for scattering losses from dislocation loops also do not agree with the experimental results. Instead, fitting to the experimental data shows that increased absorption from aluminum vacancy centers is the most likely explanation for the the prohibitively large RIA that was observed at high irradiation temperature and dose. In addition, the voids that formed in these single-crystal samples were found to align along the basal plane (a-axis) as opposed to that seen in previous observations of c-axis alignment in polycrystalline samples, which could have important implications for anisotropic swelling and other phenomena that could affect sensor performance at high neutron fluence.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Testing and Characterization to Develop a Mechanistic Explanation for Unsaturated Drift of Fiber Optic Sensors during High-Dose Irradiation

The primary limitation for any optical fiber-based sensor for nuclear reactor applications is radiation-induced attenuation (RIA) of the transmitted and/or reflected signals. Based on several recent studies, RIA is tolerable for some fused silica optical fibers with the proper choice of sensing wavelength and fiber dopants. For extreme temperature applications (> 1000°C), sapphire optical fibers have been proposed; however, recent optical transmission measurements performed on bulk sapphire samples showed prohibitively large RIA. For some sensors, radiation-induced dimensional changes in the fiber materials can also cause significant drift. Moreover, the drift that was observed in numerous experiments performed in the High Flux Isotope Reactor (HFIR), the Advanced Test Reactor, the Massachusetts Institute of Technology Reactor, and other international facilities far exceeded what would be expected based on compaction of fused silica glass. Clearly, additional work is needed to better understand the origins of both RIA and radiation-induced drift in both silica and sapphire optical fiber-based sensors before these sensors can be reliably deployed for nuclear applications. This work evaluated the underlying mechanisms that may be responsible for RIA and drift in silica and sapphire materials. First, detailed characterization was performed on bulk fused silica glass samples that were previously irradiated to different neutron fluences at different temperatures to better understand the structural changes that drive radiation-induced drift in the absence of coating effects that are discussed later. Results show that the non-monotonic compaction that occurs with increasing neutron fluence continues up to fast neutron fluences approaching 10 22 n/cm 2 , which has important implications for physics-based models that may be used to compensate for the sensor drift. Initial Raman spectroscopy and synchrotron x-ray diffraction provide insights into the nature of the structural changes. Next, detailed characterizations were performed on silica fibers with various coatings that were subjected to several different thermal treatments. The hypothesis is that the coatings convert to carbon-rich materials that compact under irradiation, putting a large compressive strain on the fiber. Out-of-pile testing confirms that both polyimide and acrylate fiber coatings convert to glassy carbon (GC) materials when heated under inert conditions, and the degree of order (i.e., graphitization) increases with increasing temperature. The results provide increasingly strong evidence that the combination of polymeric coatings and inert (or vacuum) conditions render fiber optic sensors susceptible to significant radiation-induced drift that would not otherwise exist in uncoated fibers. Finally, transmission electron microscopy was performed on bulk sapphire samples that were irradiated to two neutron fluences at different temperatures to gain insights into the potential mechanisms driving the prohibitive RIA at higher neutron fluences and temperatures. Contrary to previous hypotheses, results show that scattering from radiation-induced voids cannot explain the observed RIA. Similarly, models for scattering losses from dislocation loops also do not agree with the experimental results. Instead, fitting to the experimental data shows that increased absorption from aluminum vacancy centers is the most likely explanation for the the prohibitively large RIA that was observed at high irradiation temperature and dose. In addition, the voids that formed in these single-crystal samples were found to align along the basal plane (a-axis) as opposed to that seen in previous observations of c-axis alignment in polycrystalline samples, which could have important implications for anisotropic swelling and other phenomena that could affect sensor performance at high neutron fluence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Efficient Bayesian inference with latent Hamiltonian neural networks in No-U-Turn Sampling

When sampling for Bayesian inference, one popular approach in the computational field is to use Hamiltonian Monte Carlo (HMC) and specifically the No-U-Turn Sampler (NUTS), which automatically decides the end time of the Hamiltonian trajectory. However, HMC and NUTS can require numerous numerical gradients of the target density and can prove slow in practice when relying on computationally expensive forward models. We propose Latent Hamiltonian neural networks (L-HNNs) with HMC and NUTS for solving Bayesian inference problems. Once trained, L-HNNs do not require numerical gradients of the target density during sampling, and hence numerous evaluations of the forward computational model. Moreover, L-HNNs satisfy important properties such as perfect time reversibility and Hamiltonian conservation, making them well-suited for use within HMC and NUTS because stationarity can be shown. We also propose the integration of L-HNNs in an online error monitoring scheme, in which numerical gradients of the target density are used for a few samples whenever the L-HNNs prediction errors are large. This online error monitor scheme prevents sample degeneracy in regions of low probability density and ensures robust uncertainty quantification. We demonstrate L-HNNs in NUTS with online error monitoring on several analytical examples involving complex, heavy-tailed, and high-local-curvature probability densities. We then demonstrate the applicability of L-HNNs in NUTS to two computational case studies, namely the Allen-Cahn stochastic partial differential equation and an elliptic partial differential equation with 25 and 50 inference parameters, respectively. Overall, the L-HNNs in NUTS with online error monitoring satisfactorily inferred these probability densities. In conclusion, compared to traditional NUTS, L-HNNs in NUTS with online error monitoring required 1–2 orders of magnitude fewer numerical gradients of the target density and improved the effective sample size (ESS) per gradient (which is a measure of both the sampling quality and the computational expense) by an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Status Update on the High Precision Isotopic Measurements on High Burnup LWR Fuel in 2020

The US Department of Energy (DOE) Office of Nuclear Energy (NE) is currently investigating the feasibility of directly disposing dual-purpose (storage and transportation) canisters (DPCs) in a spent nuclear fuel (SNF) repository. Criticality during the repository performance period (10,000 years or more) is one of the major concerns related to direct disposal of DPCs, specifically as the system undergoes degradation in the repository environment and timeframe. Oak Ridge National Laboratory (ORNL) is developing an as-loaded criticality analysis methodology using full (actinides + fission products) burnup credit that exploits the inherent criticality margin associated with actual canister-specific loading configuration. Burnup credit criticality analysis requires validation of the depletion/decay codes used to generate the burned isotopic inventory of an assembly by comparing the code-predicted inventory with experimentally determined isotopic data. Currently, isotopic measurement data for boiling water reactors (BWR) SNF are limited, and additional measurements will be highly beneficial for BWR burnup credit analysis, which is essential to demonstrate disposability of BWR DPCs. Moreover, new pressurized water reactor (PWR) samples of isotopic measurements will expand the PWR sample population and consequently will greatly improve sample statistics. In turn, this will reduce uncertainty in the computational determination of the isotopic composition of commercial SNF and eliminate additional penalties currently used for lack of data. As part of the isotopic analysis task, eight diverse, high burnup (HBU) samples from several PWR rods are being dissolved and isotopically analyzed to provide high-quality measurement data for PWR SNF to reducing the uncertainties associated with PWR isotopic depletion validation. The isotopics of interest focus on nuclides important to burnup credit, shielding, and decay heat in PWR SNF. Additionally, eight high burnup BWR fuel samples from the Limerick nuclear power plant have been identified for measurement to support reducing uncertainties associated with BWR isotopic depletion validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1979

Geochemical data synthesis and analysis

Data obtained at the Goddard Flight Center were collected for the purpose of completing analyses started at Goddard in order to maximize the scientific yield of the geochemistry program which was terminated in 1977. The major analytical task undertaken was to complete Gd analyses on a large number of samples already analyzed by mass spectrometry for other rare earth element abundances at Goddard. Gd values are important for pinning down the central part of the geochemically significant rare earth abundance pattern and are especially useful in the high precision definition of the utilitarian Eu anomaly. Isotope-dilution Gd abundances were obtained for 39 samples. The data are for 27 partition-coefficient samples, six Apollo 15 and 16 breccia samples, four terrestrial impactities, and associated rock standards.

Philpotts, J. A.↗

DOE OSTI2024

Robust training of machine learning interatomic potentials with dimensionality reduction and stratified sampling

Abstract Machine learning interatomic potentials (MLIPs) enable accurate simulations of materials at scales beyond that accessible by ab initio methods and play an increasingly important role in the study and design of materials. However, MLIPs are only as accurate and robust as the data on which they are trained. Here, we present DImensionality-Reduced Encoded Clusters with sTratified (DIRECT) sampling as an approach to select a robust training set of structures from a large and complex configuration space. By applying DIRECT sampling on the Materials Project relaxation trajectories dataset with over one million structures and 89 elements, we develop an improved materials 3-body graph network (M3GNet) universal potential that extrapolates more reliably to unseen structures. We further show that molecular dynamics (MD) simulations with the M3GNet universal potential can be used instead of expensive ab initio MD to rapidly create a large configuration space for target systems. We combined this scheme with DIRECT sampling to develop a reliable moment tensor potential for titanium hydrides without the need for iterative augmentation of training structures. This work paves the way for robust high-throughput development of MLIPs across any compositional complexity.

Qi, Ji (ORCID:0000000158089931)↗

NASA NTRS2019

Was the Sun a Slow Rotator? Sodium and Potassium Constraints from the Lunar Regolith

While the Earth and Moon are generally similar in composition, a notable difference between the two is the apparent depletion in moderately volatile elements in lunar samples. This is often attributed to the formation process of the Moon, and it demonstrates the importance of these elements as evolutionary tracers. Here we show that paleo space weather may have driven the loss of a significant portion of moderate volatiles, such as sodium and potassium, from the surface of the Moon. The remaining sodium and potassium in the regolith is dependent on the primordial rotation state of the Sun. Notably, given the joint constraints shown in the observed degree of depletion of sodium and potassium in lunar samples and the evolution of activity of solar analogs over time, the Sun is highly likely to have been a slow rotator. Because the young Sun's activity was important in affecting the evolution of planetary surfaces, atmospheres, and habitability in the early Solar System, this is an important constraint on the solar activity environment at that time. Finally, as solar activity was strongest in the first billion years of the Solar System, when the Moon was most heavily bombarded by impactors, evolution of the Sun's activity may also be recorded in lunar crust and would be an important well-preserved and relatively accessible record of past Solar System processes.

Saxena, Prabal↗

NASA NTRS2019

Was the Sun a Slow Rotator? Sodium and Potassium Constraints from the Lunar Regolith

While the Earth and Moon are generally similar in composition, a notable difference between the two is the apparent depletion in moderately volatile elements in lunar samples. This is often attributed to the formation process of the Moon, and it demonstrates the importance of these elements as evolutionary tracers. Here we show that paleo space weather may have driven the loss of a significant portion of moderate volatiles, such as sodium and potassium, from the surface of the Moon. The remaining sodium and potassium in the regolith is dependent on the primordial rotation state of the Sun. Notably, given the joint constraints shown in the observed degree of depletion of sodium and potassium in lunar samples and the evolution of activity of solar analogs over time, the Sun is highly likely to have been a slow rotator. Because the young Sun's activity was important in affecting the evolution of planetary surfaces, atmospheres, and habitability in the early Solar System, this is an important constraint on the solar activity environment at that time. Finally, as solar activity was strongest in the first billion years of the Solar System, when the Moon was most heavily bombarded by impactors, evolution of the Sun's activity may also be recorded in lunar crust and would be an important well-preserved and relatively accessible record of past Solar System processes.

Prabal Saxena↗