Search NASA⌕ Search

SEARCH · Search NASA

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 829 records · Page 46

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↗

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↗

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↗

Atomistic Modeling of the Electrode–Electrolyte Interface in Li-Ion Energy Storage Systems: Electrolyte Structuring

The solid electrolyte interface (SEI) forms as a result of side reactions between the electrolyte and electrode surfaces in Li-ion batteries and can adversely impact performance by impeding Li-ion transport and diminishing the storage capacity of the battery. To gain a detailed understanding of the impact of the SEI on electrolyte structure, atomistic molecular dynamics simulations of the electrode/electrolyte interface were performed in the presence and absence of the SEI under applied voltages. The composition of the SEI was guided by a wealth of data from experiments and allowed to vary across the simulations. A novel computational approach was implemented that showed significant computational speedup compared to fully polarizable electrode simulations, yet, retained the correct qualitative physics for the electrolyte. A force-matching algorithm was used to construct a new force field for the pure electrolyte, LiPF6 in ethylene carbonate, which was developed from ab initio molecular dynamics simulations. The electrode/electrolyte interface was included using a simple, physically motivated model, which includes the polarization of the conducting graphitic electrode by the electrolyte and the application of an external voltage. Changes in the structure of the electrolyte at the interface as a function of applied voltage, the thickness of the SEI layer, and composition of the SEI provide molecular level insight into the species present at these interfaces and potential clues to the effect of the SEI on transport. It is noted that, with increasing SEI thickness and LiF content, lithium ions are drawn closer to the SEI surface, which implies that these interfaces favor desolvation and promote more rapid lithium transport.

Jorn, Ryan↗

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↗

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↗

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↗

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.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Lattice dynamics and thermodynamics for δ -plutonium from density functional theory

Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Band and Impurity States in Dimond with the (MV){sup –} (M = Si, Ge, Sn) Centers Based on ab Initio Calculations

Ab initio calculations were used to determine the structure, electronic zone dispersion, partial (atomic and orbital) densities of states, and levels of impurity–vacancy cluster (MV){sup –} (M = Si, Ge, Sn) in diamond. Wavefunctions were determined in the density functional theory approximation using hybrid exchange-correlation functional, localized basis set of Gaussians, and a supercell consisting of 192 atoms. Results are compared with data on undoped diamond and with experimental results. The dependence of the direct and indirect energy gap on the cluster type is discussed. Calculations are carried out taking into account the spin–orbit interaction to identify impurity levels in the valence band. The effect of electron–electron interaction on the position of impurity levels is discussed. The dominant contribution of the d orbitals to the density of states at the impurity level in all clusters is detected.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗