Search NASA⌕ Search

SEARCH · Search NASA

Results for “model potential calculations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Shallow Geothermal Potential of the Snake River Plain

The Snake River Plain (SRP) terrestrial heat flow and subsurface thermal regime are not well understood but are important for assessing the local geothermal resource potential, both for conventional and for Enhanced Geothermal Systems (EGS) development in the region. Resource evaluation for the SRP is complicated by the disparate data density, along with the known lateral advection of heat in the Eastern Snake River Plain Aquifer and vertical heat transport by fluids in the bounding faults, primarily in the southwestern section. Fortunately, recent studies, e.g., the Snake River Plain Play Fairway Analysis, the Idaho FORGE site, and site-specific investigations, which included drilling within the Camas Prairie and on the Mountain Home Air Force Base, near Twin Falls, and in the Eastern Snake River Plain as part of the HOTSPOT Project, add both additional drilling and geophysical data. The SMU Geothermal Laboratory has conducted detailed studies of SRP tectonics and heat flow since the 1970's and used this knowledge as part of the EGS geothermal potential estimation for the conterminous United States in 2006 and again in 2011, calculating geothermal potential from 3.5 km to 10 km depth. Recent temperature modeling refined the calculation methodology to estimate shallow (1 km to 4 km) resource potential using an improved thermal conductivity model and incorporation of shallow groundwater flow. By incorporating the new SRP geology, geophysics, and 206 thermal data sites into the SMU thermal modeling methodology, this project updates the resource estimate for the SRP, and generates new temperature-at-depth maps for the shallow subsurface (1 km to 4 km). The project results highlight the EGS potential resource areas (=150°C) and areas with more exploration risks based on minimal and/or low-quality data. The newest temperature modeling results suggest EGS potential is near five times greater in the SRP than previously estimated.

EGS↗
DOE OSTI2021

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Embedding a critical point in a hadron to quark-gluon crossover equation of state

Lattice QCD simulations have shown unequivocally that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Many model calculations predict the existence of a critical point at a value of the chemical potential where current lattice simulations are unreliable. We show how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid-gas phase transition and the three-dimensional Ising model. There are only two independent critical exponents; the relations α + 2β + γ = 2 and β(δ-1) = γ arise automatically, as does a relation between the two critical amplitudes. The resulting equation of state has parameters that may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the BNL Relativistic Heavy Ion Collider or in experiments at other accelerators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Identification of new transitions and levels in 163 Gd from $β$-decay studies

Background: Neutron-rich nuclei in the mass region around A = 160 have been and will continue to be of interest for the study of nuclear structure because of the rapid onset of deformation between 88 and 90 neutrons. The observation of detailed changes in nuclear structures within this mass region has provided and will continue to provide insight into the nuclear force. Purpose: Investigations of γ rays emitted following Eu 163 β -decay to Gd 163 have been performed for evaluation of the nuclear structure of Gd 163 . Method: Data were collected at the LeRIBSS station of the Holifield Radioactive Ion Beam Facility at Oak Ridge National Laboratory with an array of four Clover HPGe detectors for γ -rays and two plastic scintillators for β detection. The γ rays were identified as belonging to Gd 163 via mass selection and γ – γ – β , x-ray- γ , or γ – γ coincidences. Results: In total 107 new γ -ray transitions were observed in Gd 163 from 53 newly identified levels. Conclusions: We report the structure of Gd 163 has been identified for the first time. This structure has been evaluated in comparison to projected shell model, and potential energy surface calculations with good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Machine Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

Metal organic frameworks (MOFs) are a large class of porous materials and have garnered significant interest due to their large surface areas and their tunable physical and chemical properties. Numerous prior studies have been performed to screen large databases of this material class for promising DAC sorbent materials. These studies have often relied on classical model potentials. While density functional theory (DFT) calculations have been shown to be very accurate for modeling the interaction of CO2 with MOFs, such calculations are too computationally demanding for statistically significant adsorption predictions. To overcome this barrier, we developed methods for training models to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption using machine learned force fields (MLFFs). These methods were parametrized based on DFT calculations of CO2 in a flexible MOF and used to predict MOF structural properties as well as CO2 adsorption in several MOFs.

Findley, John↗
DOE OSTI2021

Theoretical study of CH 3 Cl-N 2 line shapes in the ν 1 band. Line mixing effects in QR doublets and QQ k sub-branches

Line coupling and line mixing effects of CH 3 Cl lines in the ν 1 band perturbed by N 2 have been investigated. We have taken into account the non-diagonality of the exp(–S 2 ) operator within specified line spaces as well as the k-degeneracy of the transitions (due to the double degeneracy of the j,k levels with k ≠ 0). These transitions should be considered as doublets coupled by the line mixing process. A new problem appears in the calculation when the atom-atom potential model is introduced. In order to overcome this difficulty, a pragmatic approach is proposed. Comparisons of theoretically calculated matrix elements of W with measurements of QR(j,k) doublets as well as some QQ k sub-branches, which are strongly affected by line mixing, have been made. Finally, the results show that the formalism improved in this way leads to rather accurate predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Machine Learning Design of Perovskite Catalytic Properties

Abstract Discovering new materials that efficiently catalyze the oxygen reduction and evolution reactions is critical for facilitating the widespread adoption of solid oxide fuel cell and electrolyzer (SOFC/SOEC) technologies. Here, machine learning (ML) models are developed to predict perovskite catalytic properties critical for SOFC/SOEC applications, including oxygen surface exchange, oxygen diffusivity, and area specific resistance (ASR). The models are based on trivial‐to‐calculate elemental features and are more accurate and dramatically faster than the best models based on ab initio‐derived features, potentially eliminating the need for ab initio calculations in descriptor‐based screening. The model of ASR enables temperature‐dependent predictions, has well calibrated uncertainty estimates and online accessibility. Use of temporal cross‐validation reveals the model to be effective at discovering new promising materials prior to their initial discovery, demonstrating the model can make meaningful predictions. Using the SHapley Additive ExPlanations (SHAP) approach, detailed discussion of different approaches of model featurization is provided for ML property prediction. Finally, the model is used to screen more than 19 million perovskites to develop a list of promising cheap, earth‐abundant, stable, and high performing materials, and find some top materials contain mixtures of less‐explored elements (e.g., K, Bi, Y, Ni, Cu) worth exploring in more detail.

25 ENERGY STORAGE↗
DOE OSTI2021

Total absorption spectroscopy of the β decay of 101,102 Zr and 109 Tc

The β decay of 101,102 Zr and 109 Tc was studied using the technique of total absorption spectroscopy. The experiment was performed at the National Superconducting Cyclotron Laboratory using the Summing NaI(Tl) (SuN) detector in the first-ever application of total absorption spectroscopy with a fast beam produced via projectile fragmentation. The β -decay feeding intensity and Gamow-Teller transition strength distributions were extracted for these three decays. The extracted distributions were compared to three different quasiparticle random-phase approximation (QRPA) models based on different mean-field potentials. A comparison with calculations from one of the QRPA models was performed to learn about the ground-state shape of the parent nucleus. For 101 Zr and 102 Zr, calculations assuming a pure shape configuration (oblate or prolate) were not able to reproduce the extracted distributions. These results may indicate that some type of mixture between oblate and prolate shapes is necessary to reproduce the extracted distributions. For 109 Tc, a comparison of the extracted distributions with QRPA calculations suggests a dominant oblate configuration. The other two QRPA models are commonly used to provide β-decay properties in r-process network calculations. This work shows the importance of making comparisons between the experimental and theoretical β-decay distributions, rather than just half-lives and β-delayed neutron emission probabilities, as close to the r-process path as possible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

First-principles thermodynamics of Al 10 ⁢V: An analytical treatment of localized anharmonic modes

Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 ⁢V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10⁢ V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 ⁢V phase at elevated temperatures.

anharmonic lattice dynamics↗
DOE OSTI2024

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Anomalous ellipticity dependence of the generation of near-threshold harmonics in noble gases

We measured the ellipticity dependence of the harmonic yield of He, Ne, Ar, Kr, and Xe. It is found that the harmonic yield in general decreases monotonically with increasing laser ellipticity but anomalous dependence does occur for each atom for some harmonic orders. Additionally, we found that the degree of anomaly is weakest in helium but is stronger for the heavier atoms. The anomaly also depends on laser intensity. To explain these features, we developed a quantum trajectory Monte Carlo (QTMC) model based on the quantum path integral theory to study high-order harmonic generation (HHG). The model includes the effect of Coulomb potential from the target ion, thus it is capable of calculating near-threshold harmonics quantitatively. This model reveals that the presence of Coulomb potential would generate orbiting trajectories which are responsible for the anomalous ellipticity dependence near the threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Level scheme of Nb 92 and observation of an oblate collective rotational band

Excited states of 92 Nb were investigated using the reactions 82 Se( 14 N, 4$n$) and 89 Y( 6 Li, 1 $p$2$n$). Based on these experimental results, the level scheme was extended abundantly. The structures of both the positive and negative parity states have been compared with shell model calculations using the SNET interaction. Potential energy surface calculation using the configuration-constrained method was performed to study a rotational-like band in 92 Nb, and possible configuration of the band head has been discussed. In conclusion, the calculation indicates that it may be an oblate collective rotational band.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Isospin symmetry holds beyond the proton drip line: First observation of excited states in 55Cu

Isospin symmetry, combined with shell-model description, has long been successful in reproducing nuclear structure and predicting mirror energy differences. However, these descriptions have been primarily limited to mirror nuclei with well-bound ground states. As one approaches the proton drip line, coupling to the continuum becomes important. Despite this, little is known about medium-mass nuclei at the drip line. Here we report the measurement of the excitation energies of several unbound states in a nucleus beyond the proton drip line, 55Cu, performed via γ-ray spectroscopy at the Radioactive Isotope Beam Factory (RIKEN). Several low-lying yrast states are observed to decay predominantly by γ-ray emission. In contrast, the state–predicted to be strongly populated in the one-proton knockout reaction–is notably absent, suggesting that it decays primarily via prompt proton emission. We show that even for such an unbound system, shell-model calculations based on a harmonic oscillator potential reproduce the observed mirror energy differences remarkably well in all the observed states. Based on this result, we propose that a general pattern might emerge where, in medium-mass systems at the proton drip line, the stronger Coulomb barrier together with the centrifugal barrier effectively confines the wave functions, so that the shell model is able to describe low-lying states, enabling accurate predictions of mirror energy differences between isobaric analogue states.

Pigliapoco, S [Dipartimento di Fisica e Astronomia↗
DOE OSTI2025

Systematic improvement of redox potential calculation of Fe(III)/Fe(II) complexes using a three-layer micro-solvation model

Electrochemical transformations of metal ions in aqueous media are challenging to model accurately due to the dynamic solvation structure surrounding ions at different charge states. Predictive modeling at the atomistic scale is essential for understanding these solvation architectures but is often computationally prohibitive. In this contribution, we present a simple, fast, and accurate three-layer micro-solvation model to evaluate the redox potential of metal ions in aqueous solutions. Our model, developed and validated for Fe 3+ /Fe 2+ redox potentials, combines the DFT-based geometry optimizations of the octahedral Fe complex with two layers of explicit water molecules to capture solute–solvent interactions and an implicit solvation model to account for bulk solvent effects. This approach yields accurate predictions for Fe 3+ /Fe 2+ redox potentials in water, achieving errors of 0.02 V with ωB97X-V, 0.01 V with ωB97X-D3, 0.04 V with ωB97M-V, and 0.02 V with B3LYP-D3 functionals. We further demonstrate the generality of our model by applying it to additional metal complexes, including the challenging Fe(CN) 6 3−/4− system, where our model successfully achieves close agreement with experimental values, with an error of 0.07 V and an average error of 0.21 V for all five systems. In summary, the presented simple solvation model has broad applicability and potential for enhancing computational efficiency in redox potential predictions across various chemical and industrial processes of metal ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Photovoltaic (PV) System Levelized Cost of Energy (LCOE) Evaluation with Grid Support Function Valuation and Service Lifetime Estimation

Photovoltaic (PV) systems play a critical role in renewable energy resource grid integration, and levelized cost of energy (LCOE) is commonly used to evaluate PV system feasibility in modern power grids. In this work, a revised PV system LCOE calculation model is derived to quantify the potential of LCOE reduction. Particularly, the grid support functions are valuated to offset the investment and operation costs of PV systems, which thereby reduces the LCOE. Meanwhile, PV system service lifetime is also estimated with the derived PV inverter reliability model, considering the critical components (i.e., semiconductor devices and capacitors). The case studies with field datasets are conducted to validate the effectiveness of the developed LCOE calculation model.

Zhao, Shijia↗