Search NASA⌕ Search

SEARCH · Search NASA

Results for “Extrapolation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗
DOE OSTI2025

New proton emitter 188 At implies an interaction unprecedented in heavy nuclei

We report the discovery of a new atomic nucleus 188 At, which is the heaviest proton-emitting isotope known to date. The new activity was observed through the 107 Ag( 84 Sr, 3n) 188 At fusion-evaporation reaction using the focal-plane spectrometer of the gas-filled recoil separator in the Accelerator Laboratory of the University of Jyväskylä, Finland. To fully interpret the experimental data, we have expanded the non-adiabatic quasiparticle model to treat nuclei in the beyond-lead region. The description reproduced the measured decay rate and pointed towards emission from an extremely prolate-deformed state with a dominant s 1/2 proton component in the wave function. The Thomas-Ehrman shift can be enhanced in low angular momentum states, but such effects have not been observed in heavy nuclei. The single-proton separation energy of 188 At deviates from that extrapolated from the systematics, which can be interpreted as the first evidence of this effect in heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗
DOE OSTI2025

Data-driven organic solubility prediction at the limit of aleatoric uncertainty

Abstract Small molecule solubility is a critically important property which affects the efficiency, environmental impact, and phase behavior of synthetic processes. Experimental determination of solubility is a time- and resource-intensive process and existing methods for in silico estimation of solubility are limited by their generality, speed, and accuracy. This work presents two models derived from the FASTPROP and CHEMPROP architectures and trained on BigSolDB which are capable of predicting solubility at arbitrary temperatures for a wide range of small molecules in organic solvent. Both extrapolate to unseen solutes 2–3 times more accurately than the current state-of-the-art model and we demonstrate that they are approaching the aleatoric limit (0.5–1$$\log S$$ log S ) of available test data, suggesting that further improvements in prediction accuracy require more accurate datasets. The FASTPROP-derived model (called FASTSOLV) and the CHEMPROP-based model are open source, freely accessible via a Python package and web interface, highly reproducible, and up to 2 orders of magnitude faster than current alternatives.

Science & Technology - Other Topics↗
DOE OSTI2025

Learning plasma dynamics and robust rampdown trajectories with predict-first experiments at TCV

The rampdown phase of a tokamak pulse is difficult to simulate and often exacerbates multiple plasma instabilities. To reduce the risk of disrupting operations, we leverage advances in Scientific Machine Learning (SciML) to combine physics with data-driven models, developing a neural state-space model (NSSM) that predicts plasma dynamics during Tokamak à Configuration Variable (TCV) rampdowns. The NSSM efficiently learns dynamics from a modest dataset of 311 pulses with only five pulses in a reactor-relevant high-performance regime. The NSSM is parallelized across uncertainties, and reinforcement learning (RL) is applied to design trajectories that avoid instability limits. High-performance experiments at TCV show statistically significant improvements in relevant metrics. A predict-first experiment, increasing plasma current by 20% from baseline, demonstrates the NSSM’s ability to make small extrapolations. The developed approach paves the way for designing tokamak controls with robustness to considerable uncertainty and demonstrates the relevance of SciML for fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Trivalent titanium in high-titanium lunar ilmenite

Lunar mare basalts are often rich in titanium, hosted predominantly within the mineral ilmenite (Fe 2+ Ti 4+ O 3 ). Here, we examine ilmenite in a ~3.8 billion-year-old mare basalt (Apollo rock 75035) using high-resolution electron microscopy and electron energy loss spectroscopy. A key finding is that 75035 ilmenite is itself enriched in Ti, beyond the end member of the conventional solid solution series. Using energy loss near-edge spectroscopy, we determine that the excess Ti is trivalent, with Ti 3+ accounting for 13% of the total Ti content. This discovery confirms the presence of trivalent Ti in lunar ilmenite, long hypothesized based on the Moon’s reducing environment. Accounting for the change in implied stoichiometry, a review of literature data suggests that Ti 3+ may be present in ilmenite across a wide range of lunar samples. We extrapolate known relationships from the literature to connect Ti 3+ to redox conditions, estimating the oxygen fugacity during crystallization of 75035 to be below the iron-wüstite buffer, ΔIW≤ − 1.6. Further quantifying the relationship between Ti valence state and oxygen fugacity would allow Ti 3+ -bearing ilmenite to serve as an oxybarometer able to access the reducing conditions found on many planetary bodies.

mineralogy↗
DOE OSTI2024

Dynamic mode decomposition of nonequilibrium electron-phonon dynamics: accelerating the first-principles real-time Boltzmann equation

Abstract Nonequilibrium dynamics governed by electron–phonon ( e -ph) interactions plays a key role in electronic devices and spectroscopies and is central to understanding electronic excitations in materials. The real-time Boltzmann transport equation (rt-BTE) with collision processes computed from first principles can describe the coupled dynamics of electrons and atomic vibrations (phonons). Yet, a bottleneck of these simulations is the calculation of e –ph scattering integrals on dense momentum grids at each time step. Here we show a data-driven approach based on dynamic mode decomposition (DMD) that can accelerate the time propagation of the rt-BTE and identify dominant electronic processes. We apply this approach to two case studies, high-field charge transport and ultrafast excited electron relaxation. In both cases, simulating only a short time window of ~10% of the dynamics suffices to predict the dynamics from initial excitation to steady state using DMD extrapolation. Analysis of the momentum-space modes extracted from DMD sheds light on the microscopic mechanisms governing electron relaxation to a steady state or equilibrium. The combination of accuracy and efficiency makes our DMD-based method a valuable tool for investigating ultrafast dynamics in a wide range of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
DOE OSTI2026

Quantifying the phase diagram and Hamiltonian of S = 1/2 kagome antiferromagnets: bridging theory and experiment

Spin-1/2 kagome antiferromagnets are leading candidates for realizing quantum spin liquid (QSL) ground states. While QSL ground states are predicted for the pure Heisenberg model, understanding the robustness of the QSL to additional interactions that may be present in real materials is a forefront question in the field. Here we employ large-scale density-matrix renormalization group simulations to investigate the effects of next-nearest neighbor exchange couplings J 2 and Dzyaloshinskii-Moriya interactions D, which are relevant to understanding the prototypical kagome materials herbertsmithite and Zn-barlowite. By utilizing clusters as large as XC12 and extrapolating the results to the thermodynamic limit, we precisely delineate the scope of the QSL phase, which remains robust across an expanded parameter range of J 2 and D. Direct comparison of the simulated static and dynamic spin structure factors with inelastic neutron scattering reveals the parameter space of the Hamiltonians for herbertsmithite and Zn-barlowite, and, importantly, provides compelling evidence that both materials exist within the QSL phase. These results establish a powerful convergence of theory and experiment in this most elusive state of matter.

Jiang, Shengtao [SLAC National Accelerator Laborat↗
DOE OSTI2024

Unveiling the effect of composition on nuclear waste immobilization glasses’ durability by nonparametric machine learning

Abstract Ensuring the long-term chemical durability of glasses is critical for nuclear waste immobilization operations. Durable glasses usually undergo qualification for disposal based on their response to standardized tests such as the product consistency test or the vapor hydration test (VHT). The VHT uses elevated temperature and water vapor to accelerate glass alteration and the formation of secondary phases. Understanding the relationship between glass composition and VHT response is of fundamental and practical interest. However, this relationship is complex, non-linear, and sometimes fairly variable, posing challenges in identifying the distinct effect of individual oxides on VHT response. Here, we leverage a dataset comprising 654 Hanford low-activity waste (LAW) glasses across a wide compositional envelope and employ various machine learning techniques to explore this relationship. We find that Gaussian process regression (GPR), a nonparametric regression method, yields the highest predictive accuracy. By utilizing the trained model, we discern the influence of each oxide on the glasses’ VHT response. Moreover, we discuss the trade-off between underfitting and overfitting for extrapolating the material performance in the context of sparse and heterogeneous datasets.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Contextual subspace variational quantum eigensolver calculation of the dissociation curve of molecular nitrogen on a superconducting quantum computer

Abstract We present an experimental demonstration of the Contextual Subspace Variational Quantum Eigensolver on superconducting hardware. Calculating the potential energy curve of molecular nitrogen proves challenging for many conventional quantum chemistry techniques, since static correlation dominates in the dissociation limit. Our quantum simulations retain good agreement with the Full Configuration Interaction energy, outperforming all benchmarked single-reference wavefunction techniques in capturing the bond-breaking appropriately. Moreover, our methodology is competitive with multiconfigurational approaches but at a saving of quantum resource, meaning larger active spaces can be treated for a fixed qubit allowance. To achieve this result, we deploy an error mitigation/suppression strategy comprised of Dynamical Decoupling, Measurement-Error Mitigation and Zero-Noise Extrapolation. Circuit parallelization also provides passive noise-averaging and improves the effective shot yield to reduce the measurement overhead. Furthermore, we introduce a modified adaptive ansatz construction algorithm that incorporates hardware awareness into our variational circuits, minimizing the transpilation cost for the target qubit topology.

Physics↗
DOE OSTI2025

Chemical versus physical pressure effects on the structure transition of bilayer nickelates

The observation of high-T c superconductivity (HTSC) in concomitant with pressure-induced orthorhombic-tetragonal structural transition in bilayer La 3 Ni 2 O 7 has sparked hopes of achieving HTSC by stabilizing the tetragonal phase at ambient pressure. Chemical pressure, introduced by replacing La 3+ with smaller rare-earth R 3+ has been considered as a potential route. However, our experimental and theoretical investigation reveals that such substitutions, despite causing lattice contraction, actually produce stronger orthorhombic distortions, requiring higher pressures for the structural transition. A linear extrapolation of P c versus the average size of A-site cations (<$r_A$>), yields a putative critical value of (<$r_A$>) ≈ 1.23 Å for P c ≈ 1 bar. The negative correlation between P c and (<$r_A$>) indicates that replacing La 3+ with smaller R 3+ ions is unlikely to reduce P c to ambient pressure. Instead, substituting La 3+ with larger cations like Sr 2+ or Ba 2+ might be a feasible approach. Our results provide guidance for realizing ambient-pressure HTSC in bilayer nickelates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Using magnetic dynamics to measure the spin gap in a candidate Kitaev material

Spin-orbit entangled materials have attracted widespread interest due to the novel magnetic phenomena arising from the interplay between spin-orbit coupling and electronic correlations. However, the intricate nature of spin interactions within Kiteav materials complicates the precise measurement of low-energy spin excitations. Using Na 2 Co 2 TeO 6 as an example, we study these low-energy spin excitations using the time-resolved resonant elastic x-ray scattering (tr-REXS). Our observations unveil remarkably slow spin dynamics at the magnetic peak, whose recovery timescale is several nanoseconds. This timescale aligns with the extrapolated spin gap of ~1 μeV, obtained by density matrix renormalization group (DMRG) simulations in the thermodynamic limit. The consistency demonstrates the efficacy of tr-REXS in discerning low-energy spin gaps inaccessible to conventional spectroscopic techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Impact of Sound Ordinances on the Land-Based Wind Technical Potential of the United States

The rapid rise of sound ordinances has affected wind energy siting across state, county and local jurisdictions. Quantifying this impact on wind energy potential is required for effective energy planning. However, complex turbine sound modelling combined with evolving policies have contributed to limited understanding of these policies' impact on national wind technical potential. Here we develop an approach to evaluate wind turbine sound profiles across the United States and calculate setback distances from residential structures to translate sound-related policies into spatial turbine siting criteria. Compared to a baseline scenario that ignores sound, we find an 8% reduction in wind energy capacity when accounting for existing sound ordinances. Further, extrapolating a 45-A-weighted decibel (dBA) limit nationwide shows a 29% decline, disproportionately taken from high-quality and low-cost wind capacity. This reveals sound ordinances as a major factor in determining wind technical potential that could challenge US decarbonization goals reliant on wind energy.

17 WIND ENERGY↗
DOE OSTI2024

In-beam spectroscopy reveals competing nuclear shapes in the rare isotope 62 Cr

In recent decades, rareisotope facilities have enabled the study of short-lived, neutron-rich nuclei. Their measured properties indicate that shell structure changes in the regime of unbalanced neutron-to-proton ratios compared to that of stable nuclei. In the so-called islands of inversion in the nuclear chart – around the neutron-rich nuclei 32 Mg, 42 Si and 64 Cr, for example – the textbook shell-model picture predicts the presence of spherical shapes due to the respective magic neutron numbers of 20, 28, and 40 of these nuclei. However, nuclei in these regions turn out to be deformed in their ground states. Another hallmark of these islands is shape coexistence, where a nucleus assumes different shapes with excitation energy. Furthermore, we present evidence for this phenomenon from the observation of an excited 0 + state in 62 Cr, two neutrons away from the heart of the island of inversion around neutron number N = 40. We use large-scale shell-model calculations to interpret the results and report extrapolations for the doubly-magic nucleus 60 Ca.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗