Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

High-pressure phase transition of alkali metal–transition metal deuteride Li 2 PdD 2

A combined theoretical and experimental study of lithium palladium deuteride (Li 2 PdD 2 ) subjected to pressures up to 50 GPa reveals one structural phase transition near 10 GPa, detected by synchrotron powder x-ray diffraction, and metadynamics simulations. We report the ambient-pressure tetragonal phase of Li 2 PdD 2 transforms into a monoclinic $C2/m$ phase that is distinct from all known structures of alkali metal–transition metal hydrides/deuterides. The structure of the high-pressure phase was characterized using ab initio computational techniques and from refinement of the powder x-ray diffraction data. In the high-pressure phase, the PdD 2 complexes lose molecular integrity and are fused to extended [PdD 2 ] ∞ chains. The discovered phase transition and new structure are relevant to the possible hydrogen storage application of Li 2 PdD 2 and alkali metal–transition metal hydrides in general.

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

Localized basis set for plutonium

The implementation of optimal strictly localized atomic orbitals basis for plutonium (Pu) using norm-conserving pseudopotential density-functional theory (DFT) is presented. The basis set was applied to the α, β, γ, δ, δ', and ε phases of Pu, δ-Pu surface, and δ-PuGa alloys. The computed properties of the Pu phases and δ-Pu surface were in good agreement with both available experimental data and prior DFT calculations based on plane-wave methodologies. Results for the δ-PuGa alloys were also in good agreement with experimental data. Finally, the reliability of the basis set was further demonstrated by using ab initio molecular dynamics to model the diffusion coefficient and activation barrier for atomic diffusion in a δ-PuGa alloy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

In-Silico Analysis of High Refractive Index Materials Through Principles of Materials Design

The intent of the paper is to use specific principles of Materials Design that were developed and applied in the electronics industry for enabling understanding and design of improved high refractive index materials. Further, by combining first-principle based ab-initio, semiempirical interatomic potential methods, and machine learning approaches in conjunction with experimental data, we identified specific determinants of high refractive index materials, which can be critically applied for informing materials design and accelerating discovery. Specifically, it was demonstrated that chalcogenides and perovskites as bulk materials can exhibit higher refractive indices with appropriate engineering of specific aspects of the materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Reaction pathways of cyclopentanone oxidation intermediates

Despite the promising role of cyclopentanone as a bio-derived fuel, thermodynamic and kinetic data are lacking for low-temperature oxidation regimes. This talk will present ab initio calculations of the subsequent reactivity that results from O2-addition to 2- and 3-oxo cyclopentyl radicals, including expected reaction classes such as intra-H migration, HO2-elimination, cyclic ether formation, and β-scission along with their thermodynamic parameters. Some of the rates are similar to the analogous reactions of cyclopentane, but some other reactions of cyclopentanone are very different. Comparison of the reactivity of α- and β- substituted ketone will also be discussed. These findings will be valuable for evaluating the potential of cyclopentanone as a biofuel.

09 BIOMASS FUELS↗
DOE OSTI2017

Ab initio calculations of reactions of light nuclei

An ab initio (i.e., from first principles) theoretical framework capable of providing a unified description of the structure and low-energy reaction properties of light nuclei is desirable as a support tool for accurate evaluations of crucial reaction data for nuclear astrophysics, fusion-energy research, and other applications. We present an efficient many-body approach to nuclear bound and scattering states alike, known as the ab initio no-core shell model with continuum. In this approach, square-integrable energy eigenstates of the A-nucleon system are coupled to (A-A)+A target-plus-projectile wave functions in the spirit of the resonating group method to obtain an efficient description of the many-body nuclear dynamics both at short and medium distances and at long ranges. We show that predictive results for nucleon and deuterium scattering on 4 He nuclei can be obtained from the direct solution of the Schröedinger equation with modern nuclear potentials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Magnetic Transitions in Iron Porphyrin Halides by Inelastic Neutron Scattering and Ab-initio Studies of Zero-Field Splittings

Zero-field splitting (ZFS) parameters of nondeuterated metalloporphyrins [Fe(TPP)X] (X = F, Br, I; H 2 TPP = tetraphenylporphyrin) are determined by inelastic neutron scattering (INS). The ZFS values are D = 4.49(9) cm –1 for tetragonal polycrystalline [Fe(TPP)F], and D = 8.8(2) cm –1 , E = 0.1(2) cm –1 and D = 13.4(6) cm –1 , E = 0.3(6) cm –1 for monoclinic polycrystalline [Fe(TPP)Br] and [Fe(TPP)I], respectively. Along with our recent report of the ZFS value of D = 6.33(8) cm –1 for tetragonal polycrystalline [Fe(TPP)Cl], these data provide a rare, complete determination of ZFS parameters in a metalloporphyrin halide series. The electronic structure of [Fe(TPP)X] (X = F, Cl, Br, I) has been studied by multireference ab initio methods: the complete active space self-consistent field (CASSCF) and the N-electron valence perturbation theory (NEVPT2) with the aim of exploring the origin of the large and positive zero-field splitting D of the 6 A 1 ground state. D was calculated from wave functions of the electronic multiplets spanned by the d 5 configuration of Fe(III) along with spin–orbit coupling accounted for by quasi degenerate perturbation theory. Results reproduce trends of D from inelastic neutron scattering data increasing in the order from F, Cl, Br, to I. A mapping of energy eigenvalues and eigenfunctions of the S = 3/2 excited states on ligand field theory was used to characterize the σ- and π-antibonding effects decreasing from F to I. This is in agreement with similar results deduced from ab initio calculations on CrX 6 3- complexes and also with the spectrochemical series showing a decrease of the ligand field in the same directions. A correlation is found between the increase of D and decrease of the π- and σ-antibonding energies e λ X (λ = σ, π) in the series from X = F to I. Analysis of this correlation using second-order perturbation theory expressions in terms of angular overlap parameters rationalizes the experimentally deduced trend. Furthermore, D parameters from CASSCF and NEVPT2 results have been calibrated against those from the INS data, yielding a predictive power of these approaches. Methods to improve the quantitative agreement between ab initio calculated and experimental D and spectroscopic transitions for high-spin Fe(III) complexes are proposed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Quasielastic electromagnetic scattering cross sections and world data comparisons in the GENIE Monte Carlo event generator

The usage of Monte Carlo neutrino event generators (MC$\textit{ν}$EGs) is a norm within the high-energy $\textit{ν}$ scattering community. The relevance of quasielastic (QE) energy regimes to $\textit{ν}$ oscillation experiments implies that accurate calculations of $\textit{ν}$A cross sections in this regime will be a key contributor to reducing the systematic uncertainties affecting the extraction of oscillation parameters. In spite of this, many MC$\textit{ν}$EGs utilize highly phenomenological, parametrized models of QE scattering cross sections. Moreover, a culture of validation of MC$\textit{ν}$EGs against prolific electron ($\textit{e}$) scattering data has been historically lacking. In this work, we implement new eA cross sections obtained from nuclear ab initio approaches in GENIE, the primary MC$\textit{ν}$EG utilized by the FNAL community. In particular, we utilize results from quantum MC methods which solve the many-body nuclear problem in the short-time approximation (STA), allowing consistent retention of two-nucleon dynamics which are crucial to explain available nuclear electromagnetic (electroweak) data over a wide range of energy and momentum transfers. This new implementation in GENIE is fully tested against the world QE electromagnetic data, finding agreement with available data below ~2 GeV of beam energy with the aid of a scaling function formalism. The STA is currently limited to study $\textit{A}$ ≤ 12 nuclei, however, its semi-inclusive multibody identity components are exportable to other many-body computational techniques such as auxiliary field diffusion MC which can reach $\textit{A}$ ≤ 40 systems while continuing to realize the factorization contained within the STA’s multinucleon dynamics. Together, these developments promise to make future experiments such as DUNE more accurate in their assessment of MC$\textit{ν}$EG systematics, $\textit{ν}$ properties, and potentially empower the discovery of physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Ab Initio Modeling on The Thermodynamic and Temperature-Dependent Elastic Properties of Subsystems of The FCC FeNiCoCr Medium Entropy Alloys (MEAs)

The stability, phonon spectra, thermodynamic, and temperature-dependent elasticity of subsystems of the FCC FeNiCoCr MEAs are systematically studied by the ab initio approach. Especially, the quasi-harmonic approximation (QHA) and the innovative Zentropy theory were utilized to predict the thermodynamic properties and elastic properties of FeNi, NiCo, FeNiCo, and FeNiCoCr MEAs with the consideration of magnetic transition. With the ensemble of the partition function based on the multiplicity of each magnetic configuration, the current work successfully reproduced the Curie temperature and the Schottky anomaly of heat capacity of these four MEAs purely based on the ab initio input. Meanwhile, the elastic properties of these alloys at finite temperatures are also successfully predicted with the consideration of magnetic transition. The overall results agree well with the available experimental data and CALPHAD prediction.

36 MATERIALS SCIENCE↗
NASA NTRS1982

The rotational spectra of HCNH/+/ and COH/+/ from quantum mechanical calculations

A description is provided of ab initio molecular orbital calculations designed to provide accurate predictions for the J = 1 to 0 rotational line of the candidate interstellar molecules HCNH(+) and COH(+). The former is believed to be important in the formation of both HCN and HNC in the interstellar medium. The latter, a metastable isomer of HCO(+), was first proposed as an interstellar molecule by Herbst et al. (1976). Attention is given to thermochemical arguments that this molecule can be formed in the same reactions which are proposed to form HCO(+), taking into account theoretical data which establish its stability to intramolecular rearrangement. Rotational constants are derived by applying an empirical correction to the ab initio rotational constants.

Defrees, D. J.↗
DOE OSTI2024

The ab initio non-crystalline structure database: empowering machine learning to decode diffusivity

Non-crystalline materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven exploration and design of non-crystalline materials is hampered by the absence of a comprehensive database covering a broad chemical space. In this work, we present the largest computed non-crystalline structure database to date, generated from systematic and accurate ab initio molecular dynamics (AIMD) calculations. We also show how the database can be used in simple machine-learning models to connect properties to composition and structure, here specifically targeting ionic conductivity. These models predict the Li-ion diffusivity with speed and accuracy, offering a cost-effective alternative to expensive density functional theory (DFT) calculations. Furthermore, the process of computational quenching non-crystalline structures provides a unique sampling of out-of-equilibrium structures, energies, and force landscape, and we anticipate that the corresponding trajectories will inform future work in universal machine learning potentials, impacting design beyond that of non-crystalline materials. In addition, combining diffusion trajectories from our dataset with models that predict liquidus viscosity and melting temperature could be utilized to develop models for predicting glass-forming ability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗
DOE OSTI2025

Toward a Generalizable Prediction Model of Molten Salt Mixture Density with Chemistry‐Informed Transfer Learning

Optimally designing applications of molten salts requires knowledge of their thermophysical properties over a wide range of temperatures and compositions. There exist significant gaps in existing databases and this data can be challenging to experimentally measure due to high temperatures, salt corrosivity, and salt hygroscopicity. Existing databases have been used to create Redlich–Kister (RK) models for mixture density showing improved accuracy with respect to ideal mixing assumptions, but these models require subcomponent data measurements for each new system, therefore lacking generality. In order to address generalizability and data sparsity, a transfer learning procedure is proposed to train deep neural networks (DNNs) using a combination of semi‐empirical relationships (RK), data from the thermophysical arm of the molten salt thermal properties database and universal ab initio properties of component mixtures taken from the joint automated repository for various integrated simulations (JARVIS) classical force‐field inspired descriptors database to predict density in molten salts. Herein, it is shown that DNNs predict molten salt density with an r 2 over 0.99 and a mean absolute percentage error under 1%, outperforming alternative methods.

inorganic materials↗
DOE OSTI2025

Investigating the weak charge of 48 Ca using a dispersive optical model

A new nonlocal dispersive-optical-model analysis has been carried out for neutrons and protons in 48 Ca that reproduces the weak-form-factor measurement of CREX. In addition to elastic-scattering angular distributions, total and reaction cross sections, single-particle energies, the neutron and proton numbers, and the charge distribution, the CREX-measured weak form factor has been fit to extract the neutron and proton self-energies both above and below the Fermi energy. The resulting single-particle propagators yield a weak form factor of F w = 0.125 ± 0.05 and a neutron skin of R skin = 0.152 ± 0.05 fm, in good agreement with CREX. The rearrangement of the neutron distribution to accommodate such a thin neutron skin results in the high-momentum content of the neutrons exceeding that of the protons, in contrast to what is expected from high-energy two-nucleon knockout measurements by the CLAS collaboration and ab initio asymmetric matter calculations. The present analysis also emphasizes the importance of neutron experimental data in constraining weak charge observables necessary for a precise description of neutron densities. Notably, the neutron reaction cross section and further parity-violating experiments weak form factor measurements are essential to generate a unique way to determine the 48 Ca neutron distribution in this framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

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

Quantifying Uncertainties of Microscopic Nuclear Theories

The goal of this project was to produce modernized nuclear data libraries of light ion fusion reactions and actinide fission product yields. Such libraries are essential input in both stockpile science but also fundamental science (nucleosynthesis, superheavy nuclei). When such libraries are used to model complex processes involving exotic, very short-lived nuclei that have not been measured experimentally, the underlying nuclear data come from theoretical models, the predictive power of which needs to be quantified, together with its uncertainties. By teaming up nuclear theorists and statisticians, we built comprehensive theoretical frameworks based on the most fundamental theories of fusion and fission that we coupled for the first time with a Bayesian framework with machine learning techniques for uncertainty quantification. Our results include: (1) statistical emulators for neutron-alpha scattering, and the fission product mass yields of neutron-induced fission of 239 Pu, which will enable developers to generate random realizations of the corresponding data for uncertainty quantification (UQ) in applications; (2) the first complete and most accurate ab initio prediction of the S-factor for proton capture on 7 Be, which – combined with UQ – may become the recommended value at solar energies in the next review of solar fusion cross sections; (3) statistical emulators for fission mass yields in the neutron-induced fission of 239 Pu; (4) the first quantification of how uncertainties in nuclear forces impact fission properties that are relevant in nucleosynthesis simulations to explain how heavy elements are formed in the Cosmos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

Highly Accurate Potential Energy Surface, Dipole Moment Surface, Rovibrational Energy Levels, and Infrared Line List for (32)S(16)O2 up to 8000 cm(exp -1)

A purely ab initio potential energy surface (PES) was refined with selected (32)S(16)O2 HITRAN data. Compared to HITRAN, the root-mean-squares error (RMS) error for all J=0-80 rovibrational energy levels computed on the refined PES (denoted Ames-1) is 0.013 cm(exp -1). Combined with a CCSD(T)/aug-cc-pV(Q+d)Z dipole moment surface (DMS), an infrared (IR) line list (denoted Ames-296K) has been computed at 296K and covers up to 8,000 cm(exp -1). Compared to the HITRAN and CDMS databases, the intensity agreement for most vibrational bands is better than 85-90%. Our predictions for (34)S(16)O2 band origins, higher energy (32)S(16)O2 band origins and missing (32)S(16)O2 IR bands have been verified by most recent experiments and available HITRAN data. We conclude that the Ames-1 PES is able to predict (32/34)S(16)O2 band origins below 5500 cm(exp -1) with 0.01-0.03 cm(exp -1) uncertainties, and the Ames-296K line list provides continuous, reliable and accurate IR simulations. The Ka-dependence of both line position and line intensity errors is discussed. The line list will greatly facilitate SO2 IR spectral experimental analysis, as well as elimination of SO2 lines in high-resolution astronomical observations.

Sulfur Dioxide (SO2)↗
DOE OSTI2021

Resonances of A = 4T = 1 isospin triplet states within the ab initio no-core Gamow shell model

Here, the A = 4 nuclei, i.e., 4 H, 4 He, and 4 Li , establish an interesting isospin T = 1 isobaric system. 4 H and 4 Li are unbound broad resonances, whereas 4 He is deeply bound in its ground state but unbound in all its excited states. The present situation is that experiments so far have not given consistent data on the resonances. Few-body calculations have well studied the scatterings of the 4N systems. In the present work, we provide many-body calculations of the broad resonance structures, in an ab initio framework with modern realistic interactions. It occurs that, indeed, 4 H, 4 Li, and excited 4 He are broad resonances, which is in accordance with experimental observations. The calculations also show that the first 1 - excited state almost degenerates with the 2 - ground state in the pair of mirror isobars of 4 H and 4 Li, which may suggest that the experimental data on energy and width are the mixture of the ground state and the first excited state. The T = 1 isospin triplet formed with an excited state of 4 He and ground states of 4 H and 4 Li is studied, focusing on the effect of isospin symmetry breaking.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Abstract]

Historically, the free gas approximation has been used to treat the thermal scattering of neutrons with energies below a few electron-volts (eV) in unevaluated materials. However, this method inadequately reproduces neutron scattering at these energies. Until recently, only a limited number of materials had available thermal scattering law (TSL) files/libraries in the ENDF nuclear data libraries in this energy range. With advancements in atomistic modeling techniques, such as molecular dynamics, ab-initio molecular dynamics, and density functional theory, TSL libraries have become available for many more materials. This is particularly relevant due to the rising interest in several advanced reactor systems that require novel moderator and reflector materials. While quasi-integral and integral benchmarks have been designed to validate historically important moderator materials (such as light water and polyethylene), there is currently a lack of standard validation methods for TSLs, especially when multiple conflicting TSLs exist. To address this issue, the Oak Ridge National Lab Nuclear Data group has been working on utilizing inelastic neutron scattering (INS) measurements combined with transmission (i.e., total cross section) measurements to evaluate and validate TSLs for different materials. We plan to demonstrate how this method has worked on materials such as polyethylene, lucite, and polystyrene. In addition, we will compare the newly created libraries to ENDF libraries for these materials and explain why integral benchmarks should not be used for validation when multiple TSLs exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗