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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A global potential energy surface for ArH2

We describe a simple analytic representation of the ArH2 potential energy surface which well reproduces the results of extensive ab initio electronic structure calculations. The analytic representation smoothly interpolates between the dissociated H2 and strong bonding limits. In the fitting process, emphasis is made on accurately reproducing regions of the potential expected to be important for high temperature (ca. 3000 K) collision processes. Overall, the anisotropy and H2 bond length dependence of the analytic representation well reproduce the input data.

Schwenke, David W.↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2017

An experimental and theoretical study of the thermal decomposition of C 4 H 6 isomers

The chemistry of small unsaturated hydrocarbons, such as 1,3–butadiene (1,3–C 4 H 6 ), 1,2–butadiene (1,2–C 4 H 6 ), 2–butyne (2–C 4 H 6 ) and 1–butyne (1–C 4 H 6 ), is of central importance to the modeling of combustion systems. These species are important intermediates in combustion processes, and yet their high-temperature chemistry remains poorly understood, with various dissociation and isomerization pathways proposed in the literature. Here we investigate the thermal decompositions of 1,3–C 4 H 6 , 1,2–C 4 H 6 , 2–C 4 H 6 and 1–C 4 H 6 inside a diaphragmless shock tube, at post shock total pressures of 26–261 Torr and temperatures ranging from 1428–2354 K, using laser schlieren densitometry. The experimental work has been complemented by high-level ab initio calculations, which collectively provide strong evidence that formally direct dissociation is the major channel for pyrolysis of 1,3–C 4 H 6 and 2–C 4 H 6 ; these paths have not been previously reported but are critical to reconciling the current work and disparate literature reports. The reaction mechanism presented here simulates the current experiments and experimental data from the literature very well. As a result, pressure and temperature dependent rate coefficients are given for the isomerization, formally direct and direct dissociation paths.

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

Model-free Rayleigh weight from x-ray Thomson scattering measurements

X-ray Thomson scattering (XRTS) has emerged as a powerful tool for the diagnostics of matter under extreme conditions. In principle, it gives one access to important system parameters such as the temperature, density, and ionization state, but the interpretation of the measured XRTS intensity usually relies on theoretical models and approximations. In this context, a key property is given by the Rayleigh weight that describes the electronic localization around the ions. Here, we show that it is possible to extract the Rayleigh weight directly from the experimental data without the need for any model calculations or simulations. As a practical application, we consider an experimental measurement of strongly compressed Be at the National Ignition Facility [Döppner et al., Nature 618, 270–275 (2023)]. We demonstrate that experimental results for the Rayleigh weight open up new avenues for the interpretation of XRTS experiments by matching the measurement with ab initio simulations such as density functional theory or path integral Monte Carlo. Interestingly, this new procedure leads to significantly lower density compared to previously used chemical models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as “Thiel’s set.” Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted “tier2+aug2” is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

Liu, Chi↗
DOE OSTI2016

Ab initio study of vacancy formation in cubic LaMnO{sub 3} and SmCoO{sub 3} as cathode materials in solid oxide fuel cells

Doped LaMnO{sub 3} and SmCoO{sub 3} are important solid oxide fuel cell cathode materials. The main difference between these two perovskites is that SmCoO{sub 3} has proven to be a more efficient cathode material than LaMnO{sub 3} at lower temperatures. In order to explain the difference in efficiency, we need to gain insight into the materials’ properties at the atomic level. However, while LaMnO{sub 3} has been widely studied, ab initio studies on SmCoO{sub 3} are rare. Hence, in this paper, we perform a comparative DFT + U study of the structural, electronic, and magnetic properties of these two perovskites. To that end, we first determined a suitable Hubbard parameter for the Co d–electrons to obtain a proper description of SmCoO{sub 3} that fully agrees with the available experimental data. We next evaluated the impact of oxygen and cation vacancies on the geometry, electronic, and magnetic properties. Oxygen vacancies strongly alter the electronic and magnetic structures of SmCoO{sub 3}, but barely affect LaMnO{sub 3}. However, due to their high formation energy, their concentrations in the material are very low and need to be induced by doping. Studying the cation vacancy concentration showed that the formation of cation vacancies is less energetically favorable than oxygen vacancies and would thus not markedly influence the performance of the cathode.

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

Experimental evaluation of shear modulus scaling of dynamic strength at extreme pressures

Recent progress in the development of dynamic strength experimental platforms is allowing for unprecedented insight into the assumptions used to construct constitutive models operating in extreme conditions. In this work, we make a quantitative assessment of how tantalum strength scales with its shear modulus to pressures of hundreds of gigapascals through a cross-platform examination of three dynamic strength experiments. Specifically, we make use of Split–Hopkinson pressure bar and Richtmyer–Meshkov instability experiments to assess the low-pressure strain and strain rate dependence. Concurrent examination of magnetically driven ramp-release experiments up to pressures of 350GPa allows us to examine the pressure dependence. In this work, using a modern description of the shear modulus, validated against both ab initio theory and experimental measurements, we then assess how the experimentally measured pressure dependence scales with shear modulus. Furthermore, we find that the common assumption of scaling strength linearly with the shear modulus is too soft at high pressures and offer discussion as to how descriptions of slip mediated plasticity could result in an alternative scaling that is consistent with the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Deuterium Hugoniot: Pitfalls of thermodynamic sampling beyond density functional theory

Outstanding problems in the high-pressure phase diagram of hydrogen have demonstrated the need for more accurate ab initio methods for thermodynamic sampling. One promising method that has been deployed extensively above 100 GPa is coupled electron-ion Monte Carlo (CEIMC), which treats the electronic structure with quantum Monte Carlo (QMC). However, CEIMC predictions of the deuterium principal Hugoniot disagree significantly with experiment, overshooting the experimentally determined peak compression density by 7% and lower temperature gas-gun data by well over 20%. By deriving an equation relating the predicted Hugoniot density to underlying equation of state errors, we show that QMC and many-body methods can easily spoil the error cancellation properties inherent in the Rankine-Hugoniot relation, and very likely suffer from error addition. By cross validating QMC based on systematically improvable trial functions against post-Hartree-Fock many-body methods, we find that these methods introduce errors of the right sign and magnitude to account for much of the observed discrepancy between CEIMC and experiment. In conclusion, we stress that this is not just a CEIMC problem, but that thermodynamic sampling based on other many-body methods is likely to experience similar difficulties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1995

Nucleosynthesis in AGB stars: Observation of Mg-25 and Mg-26 in IRC+10216 and possible detection of Al-26

We report the detection in the circumstellar envelope IRC+10216 of millimeter lines of the rare isotopomers (25)MgNC and (26)MgNC, as well as of a line at 234433 MHz, which could be the J= 7-6 transition of (26)AlF (an alternate, although less likely identified would be the J= 9-8 transition of NaF). The derived Mg-24:Mg-25:Mg-26 isotopic abundance ratios (78 : 11+/- 1 : 11 +/-1) are consistent with the solar system values (79.0:10.0:11.0), following Anders & Grevesse 1989). According to new calculations of evolutionary models of 3 solar mass and 5 solar mass asymptotic giant branch (AGB) stars, these ratios and the previously measured N, O and Si isotopic ratios imply that the central star had an initial mass 3 solar mass (less than or equal to M(sub *, ini) less than 5 solar mass and has already experienced many 3rd dredge-up events. From this, it can be predicted that the Al-26/Al-27 isotopics ratio lies between 0.01 and 0.08; in fact, the value derived in the case that U234433 arises from (26)AlF is Al-26/Al-27 = 0.04. The identification of the (25)MgNC and (26)MgNC lines was made possible by ab-initio quantum mechanical calculations of the molecule geometrical structure. It was confirmed through millimeter-wave laboratory measurements. The quantum mechanical calculations are briefly described and the laboratory results presented in some detail. The rotation constants B, D, H and the spin-rotation constant gamma of (25)MgNC and (26)MgNC are determined from a fit of laboratory and astronomical data.

Guelin, M.↗
DOE OSTI2020

Theoretical study of the electronic, thermodynamic, and thermo-conductive properties of y-LiAlO2 with 6Li isotope substitutions for tritium production

Lithium aluminate has attracted researchers’ interests due to its wide applications. The structural, electronic, optical, lattice phonon thermodynamic, and thermo-conductive properties of ?-LiAlO2 with 6Li-isotope substitution are investigated by density functional theory and lattice phonon dynamics. The calculated results show that ?-LiAlO2 possesses a wide band-gap of 4.63 eV. The valence band is mainly formed by p orbitals of Li, O, and Al. The calculated lattice thermal conductivity of ?-LiAlO2 is found to be in good agreement over a wide range of temperatures with the available experimental data. The calculated dielectric matrix shows anisotropic behavior along z and x (or y) axes. Due to the wide band-gap, the optical conductivity of ?-LiAlO2 is mainly contributed from the electron hopping between valence-conduction bands. Substitution of 7Li with 6Li isotope in ?-LiAlO2 leads to observable differences in the lattice phonon frequencies at higher frequencies (>8THz) and to slight changes in corresponding thermal conductivity and the infrared and Raman spectra which are in good agreement with the measured data. The relevance of these properties for tritium production through 6Li absorbing neutron in nuclear reactors is discussed.

ab initio thermodynamics, ?-LiAlO2 and 6Li-isotope↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2023

Strong Axiality in a Dysprosium(III) Bis(borolide) Complex Leads to Magnetic Blocking at 65 K

Substituted dysprosocenium complexes of the type [Dy(Cp R ) 2 ] + exhibit slow magnetic relaxation at cryogenic temperatures and have emerged as top-performing single-molecule magnets. The remarkable properties of these compounds derive in part from the strong axial ligand field afforded by the cyclopentadiene anions, and the design of analogous compounds with even stronger ligand fields is one promising route toward identifying new single-molecule magnets that retain a magnetic memory at even higher temperatures. Here, we report the synthesis and characterization of a dysprosium bis(borolide) compound, [K(18-crown-6)][Dy(BC 4 Ph 5 ) 2 ] (1), featuring the dysprosocenate anion [Dy(BC 4 Ph 5 ) 2 ] - with a pseudoaxial coordination environment afforded by two dianionic pentaphenyl borolide ligands. Variable-field magnetization data reveal open magnetic hysteresis up to 66 K, establishing 1 as a top-performing single-molecule magnet among its dysprosocenium analogues. Ac magnetic susceptibility data indicate that 1 relaxes via an Orbach mechanism above ~80 K with U eff = 1500(100) cm -1 and τ 0 = 10 -12.0(9) s, whereas Raman relaxation and quantum tunneling of the magnetization dominate at lower temperatures. Compound 1 exhibits a 100 s blocking temperature of 65 K, among the highest reported for dysprosium-based single-molecule magnets. Ab initio spin dynamics calculations support the experimental U eff and τ 0 values and enable a quantitative comparison of the relaxation dynamics of 1 and two representative dysprosocenium cations, yielding additional insights into the impact of the crystal field splitting and vibronic coupling on the observed relaxation behavior. Importantly, compound 1 represents a step toward the development of alternatives to substituted dysprosocenium single-molecule magnets with increased axiality.

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

Ground state wave functions for single-band Hubbard models from the Gutzwiller conjugate gradient minimisation theory

The Gutzwiller conjugate gradient minimisation (GCGM) theory is an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. Previous applications of GCGM provides satisfying accuracy of ground-state energy of Hubbard models. In the current work, we address the problem of whether the obtained wave function is a good approximation for the true ground state by comparing the correlation functions with the benchmark data. Additionally, our results confirms the accuracy of the reproduced ground state of the regular Hubbard model, but with some exception for the frustrated Hubbard model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Hyperspherical coupled channel calculations of energy and structure of {sup 4}He–{sup 4}He–Li{sup +} and its isotopic combinations

The ground state vibrational energy and spatial features of {sup 4}He–{sup 4}He–Li{sup +} and its triatomic isotopic complexes are studied using the slow variable discretization (SVD) method in the hyperspherical coordinates for the zero total angular momentum. Our results show that the dominant structure of the system is an isosceles triangle with the shorter side associated with the two Li{sup +}–He distances using the sum-of-potential approximation. Corrections caused by the induced dipole-induced dipole interactions on the He atoms are also investigated. The effects are seen to be small and have a minor influence on the binding energy and the structure of present system. The results are also compared with the full ab initio calculations including all the three-body interactions and information of three-body corrections is obtained.

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

Theory and modeling of relevance to prompt-NO formation at high pressure

An improved understanding of NO x formation at high pressures would be of considerable utility to efforts to develop advanced combustion devices. Here a combination of theoretical and modeling studies are implemented in an effort to improve the accuracy of models for the prompt NO process, which is the dominant source of NO under many conditions, and to improve our understanding of the role of this process at high pressures. The theoretical effort implements state-of-the-art treatments of NCN thermochemistry, the interrelated CH + N 2 and NCN + H kinetics, and the kinetics of the NCN + OH reaction. For both reaction systems, we implement high level ab initio transition state theory based master equation simulations paying particular attention to the role of stabilization processes. For the NCN + H kinetics we include a treatment of inter-system crossing. The modeling effort focuses on exploring the role of pressure and prompt NO for premixed laminar flames at pressures ranging from 1 to 15 atm, via a comparison with the available experimental data. Additional simulations at higher pressures further explore the mechanistic changes at the pressures of relevance to applied combustion devices (e.g., 100 atm).

03 NATURAL GAS↗
DOE OSTI2021

Non-Adiabatic Excited State Molecular Dynamics Methodologies: comparison and convergence

Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum–classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.

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

Ab initio potential energy curves of the valence, Rydberg, and ion-pair states of iodine monochloride, ICl

We present for the first time a coherent ab initio study of 39 states of valence, Rydberg, and ion-pair character of the diatomic interhalogen ICl species through large scale multireference variational methods including spin-orbit effects coupled with quantitative basis sets. Various avoided crossings are responsible for a non-adiabatic behaviour creating a wonderful vista for its theoretical description. Our molecular constants are compared with all available experimental data with the aim to assist experimentalists especially in the high energy regime of up to ∼95 000 cm{sup −1}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗