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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Molecular-scale insights into the electrical double layer at oxide-electrolyte interfaces

The electrical double layer (EDL) at metal oxide-electrolyte interfaces critically affects fundamental processes in water splitting, batteries, and corrosion. However, limitations in the microscopic-level understanding of the EDL have been a major bottleneck in controlling these interfacial processes. Herein, we use ab initio-based machine learning potential simulations incorporating long-range electrostatics to unravel the molecular-scale picture of the EDL at the prototypical anatase TiO 2 -electrolyte interface under various pH conditions. Our large-scale simulations, capable of capturing interfacial water dissociation/recombination reactions and electrolytic proton transport, provide unprecedented insights into the detailed structure of the EDL. Moreover, the larger capacitance of the EDL under basic relative to acidic conditions, originating from the higher affinity of the cations for the oxide surface, is found to give rise to distinct charging mechanisms on negative and positive surfaces. Our results are validated by the agreement between the computed EDL capacitance and experimental data.

Chemical physics↗

DOE OSTI2014

Accurate first-principles calculations for {sup 12}CH{sub 3}D infrared spectra from isotopic and symmetry transformations

Accurate variational high-resolution spectra calculations in the range 0-8000 cm{sup −1} are reported for the first time for the monodeutered methane ({sup 12}CH{sub 3}D). Global calculations were performed by using recent ab initio surfaces for line positions and line intensities derived from the main isotopologue {sup 12}CH{sub 4}. Calculation of excited vibrational levels and high-J rovibrational states is described by using the normal mode Eckart-Watson Hamiltonian combined with irreducible tensor formalism and appropriate numerical procedures for solving the quantum nuclear motion problem. The isotopic H→D substitution is studied in details by means of symmetry and nonlinear normal mode coordinate transformations. Theoretical spectra predictions are given up to J = 25 and compared with the HITRAN 2012 database representing a compilation of line lists derived from analyses of experimental spectra. The results are in very good agreement with available empirical data suggesting that a large number of yet unassigned lines in observed spectra could be identified and modeled using the present approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Transport coefficients of warm dense matter from Kohn-Sham density functional theory

We present a comprehensive study of transport coefficients including DC electrical conductivity and related optical properties, electrical contribution to the thermal conductivity, and the shear viscosity via ab initio molecular dynamics and density functional theory calculations on the “priority 1” cases from the “Second Charged-Particle Transport Coefficient Workshop” [Stanek et al., Phys. Plasmas (to be published 2024)]. The purpose of this work is to carefully document the entire workflow used to generate our reported transport coefficients, up to and including our definitions of finite size and statistical convergence, extrapolation techniques, and choice of thermodynamic ensembles. In pursuit of accurate optical properties, we also present a novel, simple, and highly accurate algorithm for evaluating the Kramers–Kronig relations. These heuristics are often not discussed in the literature, and it is hoped that this work will facilitate the reproducibility of our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Shock compression of strongly correlated oxides: A liquid-regime equation of state for cerium(IV) oxide

Here, the shock Hugoniot for full-density and porous CeO 2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ = 2.5 to 20 g/cm 3 . The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO 2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO 2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO 2 . In addition, a revised multiphase SESAME-type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2018

Material Properties for the Interiors of Massive Giant Planets and Brown Dwarfs

Here, we present thermodynamic material and transport properties for the extreme conditions prevalent in the interiors of massive giant planets and brown dwarfs. They are obtained from extensive ab initio simulations of hydrogen–helium mixtures along the isentropes of three representative objects. In particular, we determine the heat capacities, the thermal expansion coefficient, the isothermal compressibility, and the sound velocity. Important transport properties such as the electrical and thermal conductivity, opacity, and shear viscosity are also calculated. Further results for associated quantities, including magnetic and thermal diffusivity, kinematic shear viscosity, as well as the static Love number k 2 and the equidistance, are presented. In comparison to Jupiter-mass planets, the behavior inside massive giant planets and brown dwarfs is stronger dominated by degenerate matter. We discuss the implications on possible dynamics and magnetic fields of those massive objects. The consistent data set compiled here may serve as a starting point to obtain material and transport properties for other substellar H–He objects with masses above one Jovian mass and finally may be used as input for dynamo simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Experimental and kinetic modeling study of tetralin: A naphtheno-aromatic fuel for gasoline, jet and diesel surrogates

Distillate fuels contain significant proportions of naphtheno-aromatic components and tetralin is a suitable surrogate component to represent this molecular moiety. The presence of aromatic and naphthyl rings makes kinetic modeling of tetralin very challenging. Primary radicals formed during the oxidation of tetralin can be aryl, benzylic or paraffinic in nature. Using available information on reaction paths and rate constants of naphthenes and alkyl-aromatics, in this study a kinetic model of tetralin has been developed with emphasis on low-temperature chemistry and high-pressure conditions. Due to the lack of high-level quantum chemical calculations on reaction pathways of tetralin, analogous rates from ab-initio studies on benzylic and paraffinic radicals have been adopted here. Some modifications to the reaction rate rules are incorporated to account for the unique characteristics of tetralin's molecular structure. Important reaction channels have been identified using reaction path and brute force sensitivity analyses. In order to investigate the model performance at low temperatures, new experiments are carried out in a rapid compression machine on blends of tetralin and 3-methylpentane. Blending of low-reactivity tetralin with a high-reactivity alkane allowed the investigation of tetralin ignition at very low temperatures (665 – 856 K). The kinetic model developed in the current study is found to predict the current experiments and literature data adequately. The new model will aid in high-fidelity surrogate predictions at engine-relevant conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Electron–Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Quantum scattering of HC 5 N and para -H 2 on a new potential energy surface

In the interstellar medium (ISM), non-local thermodynamic equilibrium situations are common due to low density, and one needs to consider the effect of molecular collisions in order to interpret the observations. Among the species detected in the ISM, cyanopolyynes, with the general molecular formula HC 2n+1 N (n = 1, 2, …), are characterized by large dipole moments and small rotational constants and constitute an indispensable class of candidates for the sensitive tracers of local density and temperature. We present a study of the collisional (de-) excitation of HC 5 N by para -H 2 (p-H 2 ) in its ground rotational state, namely HC 5 N ( j 1 ) + H 2 ( j 2 = 0) → HC 5 N (j$_1^′$) + H2 (j$_2^′$ = 0), where j 1 (or j$_1^′$) and j 2 (or j$_2^′$) denote the initial (or final) rotational quantum numbers of HC 5 N and H 2 , respectively. We performed the quantum scattering calculations at low collision energy using a new four-dimensional ab initio potential energy surface. In the regime where p-H 2 remains in its rotational ground state, converged cross sections did not require including excited rotational states of p-H 2 in the rotational basis. State-to-state cross sections were computed by means of the quantum-mechanical close-coupling (CC) method and the coupled states (CS) approximation, and rate coefficients for the first 61 levels of HC 5 N were computed for the first time up to 20 K with the CC approach and up to 50 K with the CS method. CC and CS results were found to agree well at temperatures up to 20 K. Finally, these data should allow a more accurate derivation of the HC 5 N abundance in molecular clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

The muon g-2 experiment at Fermilab

The E989 collaboration has recently published the most precise measurement of the muon anomalous magnetic moment $a_\mu$ with an uncertainty of $0.46\,\mathrm{ppm}$. The new experimental world average of $a_\mu$ deviates by 4.2 standard deviations from the Standard Model prediction provided by the Muon g-2 Theory Initiative. Emerging results from ab-initio lattice QCD calculations for hadron-physics contributions allow to scrutinize this tantalizing hint for physics beyond the Standard Model for the first time in a three way comparison. To extract the value of $a_\mu$ a frequency comparison experiment is performed with spin-polarized muons confined in a superbly controlled electric and magnetic field environment. A measurement of the difference between the Larmor and the cyclotron frequency, the so-called anomalous spin precession frequency, is combined with a high-precision measurement of the magnetic field environment which is performed using nuclear magnetic resonance techniques. I discussed the most recent result of the Fermilab Muon g-2 experiment's first physics data run performed in 2018. I reported on the experimental improvements implemented for the subsequent science runs that will allow the collaboration to achieve the ultimate goal of $0.14\,\mathrm{ppm}$ uncertainty on $a_\mu$.

Fertl, Martin↗

DOE OSTI2016

Electronic Coupling Calculations for Bridge-Mediated Charge Transfer Using Constrained Density Functional Theory (CDFT) and Effective Hamiltonian Approaches at the Density Functional Theory (DFT) and Fragment-Orbital Density Functional Tight Binding (FODFTB) Level

In this paper, four methods to calculate charge transfer integrals in the context of bridge-mediated electron transfer are tested. These methods are based on density functional theory (DFT). We consider two perturbative Green's function effective Hamiltonian methods (first, at the DFT level of theory, using localized molecular orbitals; second, applying a tight-binding DFT approach, using fragment orbitals) and two constrained DFT implementations with either plane-wave or local basis sets. To assess the performance of the methods for through-bond (TB)-dominated or through-space (TS)-dominated transfer, different sets of molecules are considered. For through-bond electron transfer (ET), several molecules that were originally synthesized by Paddon-Row and co-workers for the deduction of electronic coupling values from photoemission and electron transmission spectroscopies, are analyzed. The tested methodologies prove to be successful in reproducing experimental data, the exponential distance decay constant and the superbridge effects arising from interference among ET pathways. For through-space ET, dedicated p-stacked systems with heterocyclopentadiene molecules were created and analyzed on the basis of electronic coupling dependence on donor-acceptor distance, structure of the bridge, and ET barrier height. The inexpensive fragment-orbital density functional tight binding (FODFTB) method gives similar results to constrained density functional theory (CDFT) and both reproduce the expected exponential decay of the coupling with donor-acceptor distances and the number of bridging units. Finally, these four approaches appear to give reliable results for both TB and TS ET and present a good alternative to expensive ab initio methodologies for large systems involving long-range charge transfers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17 O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and 17 O NMR quadrupolar relaxation rates are investigated in the series of complexes UO 2 2+ , UO 2 (CO 3 ) 3 4– , and (UO 2 ) 3 (CO 3 ) 6 6– . Car–Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO 2 2+ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl 17 O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO 2 (CO 3 ) 3 4– and (UO 2 ) 3 (CO 3 ) 6 6– complexes show a solvent dependence of the terminal carbonate 17 O and 13 C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin–orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique 13 C chemical shifts is correctly reproduced for (UO 2 ) 3 (CO 3 ) 6 6– . NMR relaxation rate data support large 17 O peak widths, but remain below those noted in the experimental literature. Furthermore, comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Water dimer equilibrium constant calculation: A quantum formulation including metastable states

We present a full quantum evaluation of the water second virial coefficient B(T) based on the Takahashi-Imada second order approximation. As the associated trace Tr[e{sup −βH{sub A}{sub B}}−e{sup −βH{sub A}{sub B}{sup o}}] is performed in the coordinate representation, it does also include contribution from the whole continuum, i.e., resonances and collision pairs of monomers. This approach is compared to a Path Integral Monte Carlo evaluation of this coefficient by Schenter [J. Chem. Phys. 117, 6573 (2002)] for the TIP4P potential and shown to give extremely close results in the low temperature range (250–450 K) reported. Using a recent ab initio flexible potential for the water dimer, this new formulation leads to very good agreement with experimental values over the whole range of temperatures available. The virial coefficient is then used in the well known relation K{sub p}(T) = −(B(T) − b{sub M})/RT where the excluded volume b{sub M} is assimilated to the second virial coefficient of pure water monomer vapor and approximated from the inner repulsive part of the interaction potential. This definition, which renders b{sub M} temperature dependent, allows us to retrieve the 38 cm{sup 3} mol{sup −1} value commonly used, at room temperature. The resulting values for K{sub p}(T) are in agreement with available experimental data obtained from infrared absorption spectra of water vapor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Shock-induced metallization of polystyrene along the principal Hugoniot investigated by advanced thermal density functionals

To date, none of the ab-initio molecular dynamics simulations of polystyrene, often used as an ablator material in inertial confinement fusion targets, with the standard ground-state exchange-correlation (XC) functional in density-functional theory can satisfactorily agree with experiments in terms of reflectivity measurements. Here, we use recently developed thermal strongly constrained and appropriately normed Laplacian dependent meta-generalized gradient approximation XC density functional (T-SCAN-L) and thermal hybrid XC density functional (KDT0) to show that the inclusion of thermal and inhomogeneity effects is crucial for accurate prediction of structural evolution and corresponding insulator-metal transition (IMT) during shock compression. Optical reflectivity calculated as an indicator of IMT is in perfect accord with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1977

Theory of inelastic ion-atom scattering at low and intermediate energies

Ab initio calculations are presented of inelastic energy loss and ionization phenomena associated with Ar(+)-Ar collisions at small distances of closest approach and for laboratory collision energies ranging from several keV to several hundred keV. Outer-shell excitations are handled statistically; inner-shell excitations are calculated from the viewpoint of quasidiabatic molecular orbital promotion. Auger electron yield, average state of ionization, and average inelastic energy loss are calculated per collision as a function of distance of closest approach of the collision partners for several laboratory collision energies. Average charge-state probabilities per collision partner are calculated as a function of the average inelastic energy loss per atom. It is shown that the structure in the data is due to the underlying structure in the inner-shell independent-electron quasimolecular promotion probabilities.

Schmid, G. B.↗

NASA NTRS1995

Transition probabilities for the UV0.01 multiplet in N III

We report on large-scale ab initio multiconfiguration Hartree-Fock calculations for the UV0.01 multiplet, 2s(sup 2)2p(sup 2)P(sub J) - 2s2p(sup 2 4)P(sub J prime), in N III. The resulting transition probabilities agree very well with recent semiempirical calculations, and the lifetimes for two of the three upper levels agree with experiments. The deviation for the third level is discussed. Comparisons made with the highest quality IUE echelle spectra available -- those of RR Tel and V1016 Cyg (both photoionized sources with electron densities below 10(exp 8)/cu cm) -- show that computed branching ratios of lines sharing a common upper level are in agreement with observations to within uncertainties of +/- 10%. High-quality solar limb data or stellar data from the Hubble Space Telescope (HST) could, in principle, be used to determine whether the theoretical or measured lifetimes for the discrepant level are in error. Unfortunately, stellar data for high-density plasmas (N(sub e) greater than 10(exp 11)/cu cm are needed) do not yet exist, and existing solar data lack the photometric precision to address this problem.

Brage, Tomas↗