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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm –3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.

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

Origins of Singlet Fission in Solid Pentacene from an ab initio Green’s Function Approach

We develop a new first-principles approach to predict and understand rates of singlet fission with an ab initio Green's-function formalism based on many-body perturbation theory. Starting with singlet and triplet excitons computed from a GW plus Bethe-Salpeter equation approach, we calculate the exciton--bi-exciton coupling to lowest order in the Coulomb interaction, assuming a final state consisting of two non-interacting spin-correlated triplets with finite center-of-mass momentum. For crystalline pentacene, symmetries dictate that the only purely Coulombic fission decay from a bright singlet state requires a final state consisting of two inequivalent nearly degenerate triplets of nonzero, equal and opposite, center-of-mass momenta. For such a process, we predict a singlet lifetime of 40 to 150 fs, in very good agreement with experimental data, indicating that this process can dominate singlet fission in crystalline pentacene. As a result, our approach is general and provides a framework for predicting and understanding multiexciton interactions in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Rovibrational spectroscopic constants and fundamental vibrational frequencies for isotopologues of cyclic and bent singlet HC{sub 2}N isomers

Through established, highly accurate ab initio quartic force fields, a complete set of fundamental vibrational frequencies, rotational constants, and rovibrational coupling and centrifugal distortion constants have been determined for both the cyclic 1 {sup 1} A' and bent 2 {sup 1} A' DCCN, H{sup 13}CCN, HC{sup 13}CN, and HCC{sup 15}N isotopologues of HCCN. Spectroscopic constants are computed for all isotopologues using second-order vibrational perturbation theory (VPT2), and the fundamental vibrational frequencies are computed with VPT2 and vibrational configuration interaction (VCI) theory. Agreement between VPT2 and VCI results is quite good, with the fundamental vibrational frequencies of the bent isomer isotopologues in accord to within a 0.1-3.2 cm{sup –1} range. Similar accuracies are present for the cyclic isomer isotopologues. The data generated here serve as a reference for astronomical observations of these closed-shell, highly dipolar molecules using new, high-resolution telescopes and as reference for laboratory studies where isotopic labeling may lead to elucidation of the formation mechanism for the known interstellar molecule: X {sup 3} A' HCCN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Towards the extraction of the crystal cell parameters from pair distribution function profiles

The approach based on atomic pair distribution function (PDF) has revolutionized structural investigations by X-ray/electron diffraction of nano or quasi-amorphous materials, opening up the possibility of exploring short-range order. However, the ab initio crystal structural solution by the PDF is far from being achieved due to the difficulty in determining the crystallographic properties of the unit cell. A method for estimating the crystal cell parameters directly from a PDF profile is presented, which is composed of two steps: first, the type of crystal cell is inferred using machine-learning approaches applied to the PDF profile; second, the crystal cell parameters are extracted by means of multivariate analysis combined with vector superposition techniques. The procedure has been validated on a large number of PDF profiles calculated from known crystal structures and on a small number of measured PDF profiles. The lattice determination step has been benchmarked by a comprehensive exploration of different classifiers and different input data. The highest performance is obtained using the k -nearest neighbours classifier applied to whole PDF profiles. Descriptors calculated from the PDF profiles by recurrence quantitative analysis produce results that can be interpreted in terms of PDF properties, and the significance of each descriptor in determining the prediction is evaluated. The cell parameter extraction step depends on the cell metric rather than its type. Monometric, dimetric and trimetric cells have top-1 estimates that are correct 40, 20 and 5% of the time, respectively. Promising results were obtained when analysing real nanocrystals, where unit cells close to the true ones are found within the top-1 ranked solution in the case of monometric cells and within the top-6 ranked solutions in the case of dimetric cells, even in the presence of a crystalline impurity with a weight fraction up to 40%.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural and electronic characteristics of amorphous Ge 8 Sb 2 Te 11

GeTe-rich phase-change materials can be utilized in rewriteable optical memory due to the large contrast in reflectivity between amorphous and crystalline phases. Here we explored the structure and electronic properties of amorphous Ge 8 Sb 2 Te 11 using ab initio molecular dynamics simulations. The results indicate that amorphous Ge 8 Sb 2 Te 11 is dominantly composed of Ge-, Sb- and Te-centered octahedrons with distortions, while 30.4 % of Ge-centered clusters are in the form of tetrahedrons which are randomly distributed. The 5-fold rings possess a large proportion, and the Gesingle bondTe and Sbsingle bondTe bonds present larger formation energies than other bonds, leading to the ABAB bonding sequence (A: Ge and Sb, B: Te). The lone pair electrons locating at the opposite direction of bonds possess a large fraction of 14.8 %, which may enhance the distortions of local clusters. Importantly, these structural properties lead to the high stability of amorphous Ge 8 Sb 2 Te 11 and thus long data retention in the information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Communication: A benchmark-quality, full-dimensional ab initio potential energy surface for Ar-HOCO

A full-dimensional, global ab initio potential energy surface (PES) for the Ar-HOCO system is presented. The PES consists of a previous intramolecular ab initio PES for HOCO [J. Li, C. Xie, J. Ma, Y. Wang, R. Dawes, D. Xie, J. M. Bowman, and H. Guo, J. Phys. Chem. A 116, 5057 (2012)], plus a new permutationally invariant interaction potential based on fitting 12 432 UCCSD(T)-F12a/aVDZ counterpoise-corrected energies. The latter has a total rms fitting error of about 25 cm{sup −1} for fitted interaction energies up to roughly 12 000 cm{sup −1}. Two additional fits are presented. One is a novel very compact permutational invariant representation, which contains terms only involving the Ar-atom distances. The rms fitting error for this fit is 193 cm{sup −1}. The other fit is the widely used pairwise one. The pairwise fit to the entire data set has an rms fitting error of 427 cm{sup −1}. All of these potentials are used in preliminary classical trajectory calculations of energy transfer with a focus on comparisons with the results using the benchmark potential.

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

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2015

Unified description of 6 Li structure and deuterium- 4 He dynamics with chiral two- and three-nucleon forces

Here, we provide a unified ab initio description of the 6 Li ground state and elastic scattering of deuterium (d) on 4 He (α) using two- and three-nucleon forces from chiral effective field theory. We analyze the influence of the three-nucleon force and reveal the role of continuum degrees of freedom in shaping the low-lying spectrum of 6 Li. The calculation reproduces the empirical binding energy of 6 Li, yielding an asymptotic D- to S-state ratio of the 6 Li wave function in the d+α configuration of –0.027, in agreement with a determination from 6 Li– 4 He elastic scattering, but overestimates the excitation energy of the 3 + state by 350 keV. The bulk of the computed differential cross section is in good agreement with data. These results endorse the application of the present approach to the evaluation of the 2 H(α,γ) 6 Li radiative capture, responsible for the big-bang nucleosynthesis of 6 Li.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Energy of Phonons and Zero-Point Vibrations in Compressed Rare-Gas Crystals

A dynamic matrix of rare-gas crystals is constructed within the model of deformable and polarizable atoms based on the nonempirical short-range repulsion potential taking into account the three-body interaction and deformation of electron shells of dipole-type atoms within the pair and three-body approximations. Ab initio calculations of the phonon energy for compressed rare-gas crystals are performed at two and ten mean-value points of the Chadi–Kohen method in a wide pressure range. It is shown that the contribution of three-body forces to the phonon frequencies due to overlap of electron shells of neighboring atoms is small as compared to the pair interaction even at a high pressure and most pronounced for Xe. The contributions of the deformation of electron shells within the pair and three-body approximations differ for different mean-value points and increase with an increase in pressure. The zero-point energies calculated by the Chadi–Kohen method for crystals of the Ne–Xe series are compared with the known experimental data at p = 0 and calculation results obtained by other researchers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

First-principles chemical kinetic modeling of methyl trans -3-hexenoate epoxidation by HO 2

The design of innovative combustion processes relies on a comprehensive understanding of biodiesel oxidation kinetics. The present study aims at unraveling the reaction mechanism involved in the epoxidation of a realistic biodiesel surrogate, methyl trans -3-hexenoate, by hydroperoxy radicals using a bottom-up theoretical kinetics methodology. The obtained rate constants are in good agreement with experimental data for alkene epoxidation by HO 2 . The impact of temperature and pressure on epoxidation pathways involving H-bonded and non-H-bonded conformers was assessed. As a result, the obtained rate constant was finally implemented into a state-of-the-art detailed combustion mechanism, resulting in fairly good agreement with engine experiments.

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

Coulomb sum rule for 4 He and 16 O from coupled-cluster theory

We demonstrate the capability of coupled-cluster theory to compute the Coulomb sum rule for the 4 He and 16 O nuclei using interactions from chiral effective field theory. We perform several checks, including a few-body benchmark for 4He. We provide an analysis of the center-of-mass contaminations, which we are able to safely remove. We then compare with other theoretical results and experimental data available in the literature, obtaining a fair agreement. This is a first and necessary step towards initiating a program for computing neutrino-nucleus interactions from first principles in coupled-cluster theory and supporting the experimental long-baseline neutrino program with a state-of-the-art theory that can reach medium-mass nuclei.

6 ≤ A ≤ 19A ≤ 5↗
DOE OSTI2014

A new ab initio potential energy surface and infrared spectra for the Ar–CS{sub 2} complex

We report a new three-dimensional potential energy surface for Ar–CS{sub 2} involving the Q{sub 3} normal mode for the υ{sub 3} antisymmetric stretching vibration of the CS{sub 2} molecule. The potential energies were calculated using the supermolecular method at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples, using augmented correlation-consistent quadruple-zeta basis set plus midpoint bond functions. Two vibrationally averaged potentials with CS{sub 2} at both the ground (υ = 0) and the first excited (υ = 1)υ{sub 3} vibrational states were generated from the integration of the three-dimensional potential over the Q{sub 3} coordinate. Each potential was found to have a T-shaped global minimum and two equivalent linear local minima. The radial discrete variable representation /angular finite basis representation method and the Lanczos algorithm were applied to calculate the rovibrational energy levels. The calculated band origin shift of the complex (0.0622 cm{sup −1}) is very close to the observed one (0.0671 cm{sup −1}). The predicted infrared spectra and spectroscopic parameters based on the two averaged potentials are in excellent agreement with the available experimental data.

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

Theoretical Studies on the Photoelectron and Absorption Spectra of MnO 4 – and TcO 4 –

The tetraoxo pertechnetate anion (TcO4-) is of great interest for nuclear waste management and radiopharmceuticals. To elucidate its electronic structure and to compare with that of its lighter congener MnO4-, the photoelectron and electronic absorption spectra of MnO4 - and TcO4 - are investigated with density functional theory (DFT) and ab initio wave function theory (WFT). The vertical electron detachment energies (VDEs) of MnO4 - obtained with the CR-EOM-CCSD(T) method are in good agreement with the lowest two experimental VDEs; the differences are less than 0.1 eV, representing a significant improvement over the IP-EOMCCSD(T) result in the literature. Combining our CCSD(T) and CR-EOMCCSD( T) results, the first five VDEs of TcO4 - are estimated between 5 and 10 eV with an estimated accuracy of about ±0.2 eV. The vertical excitation energies are determined by using TD-DFT, CR-EOM-CCSD(T), and RASPT2 methods. The excitation energies and the assignments of the spectra are analyzed and partly improved. They are compared with reported SAC-CI results and available experimental data. Both dynamic and nondynamic electron correlations are important in the ground and excited states of MnO4 - and TcO4 -. Nondynamical correlations are particularly relevant in TcO4 - for reliable prediction of excitation energies. In TcO4 - one Rydberg state interlaces but does not mix with the valence excited states, and it disappears in the condensed phase.

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

A computational study of diffusion in a glass-forming metallic liquid

In this study, liquid phase diffusion plays a critical role in phase transformations (e.g. glass transformation and devitrification) observed in marginal glass forming systems such as Al-Sm. Controlling transformation pathways in such cases requires a comprehensive description of diffusivity, including the associated composition and temperature dependencies. In our computational study, we examine atomic diffusion in Al-Sm liquids using ab initio molecular dynamics (AIMD) and determine the diffusivities of Al and Sm for selected alloy compositions. Non-Arrhenius diffusion behavior is observed in the undercooled liquids with an enhanced local structural ordering. Through assessment of our AIMD result, we construct a general formulation for Al-Sm liquid, involving a diffusion mobility database that includes composition and temperature dependence. A Volmer-Fulcher-Tammann (VFT) equation is adopted for describing the non-Arrhenius behavior observed in the undercooled liquid. Furthermore, the composition dependence of diffusivity is found quite strong, even for the Al-rich region contrary to the sole previous report on this binary system. The model is used in combination with the available thermodynamic database to predict specific diffusivities and compares well with reported experimental data for 0.6 at.% and 5.6 at.% Sm in Al-Sm alloys.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Ab Initio Characterization of Triatomic Bromine Molecules of Potential Interest in Stratospheric Chemistry

The equilibrium structures, harmonic vibrational frequencies, quadratic force fields, dipole moments, and IR intensities of several triatomic bromine compounds of known or potential importance in stratospheric ozone depletion chemistry have been determined using the CCSD(T) electron correlation method in conjunction with a basis set of triple zeta double polarized (TZ2P) quality. Specifically, the molecules included in the present study are HOBr, HBrO, FOBr, FBrO, BrNO, BrON, Br2O, BrBrO, BrCN, BrNC, ClOBr, ClBrO, and BrClO. Very accurate isomeric energy differences have also been determined at the CCSD(T) level with atomic natural orbital basis sets that include through g-type functions. In most cases, the isomer with a normal neutral Lewis dot structure is the lowest energy form, with the single exception that FBRO is predicted to be 11.1 kcal/mol (0 K) lower in energy than FOBr. In all cases, however, the hypervalent isomer is more stable relative to the isomer with a normal Lewis dot structure as compared to the chlorine analogs. Consistent with this observation, the energy of the last three molecules given above increases in the order ClOBr less than ClBrO less than BrClO. The CCSD(T)/TZ2P geometries and vibrational frequencies are in good agreement with the available experimental data. Heats of formation are determined for all species using a combination of theoretical isomeric, homodesmic, and isodesmic reaction energies. The accuracy of these quantities is ultimately dependent on the reliability of the experimental heat of formation of HOBr.

Lee. Timothy J.↗
NASA NTRS1995

Ab Initio Characterization of Triatomic Bromine Molecules of Potential Interest in Stratospheric Chemistry

The equilibrium structures, harmonic vibrational frequencies, quadratic force fields, dipole moments, and IR intensities of several triatomic bromine compounds of known or potential importance in stratospheric ozone depletion chemistry have been determined using the CCSD(T) electron correlation method in conjunction with a basis set of triple zeta double polarized (TZ2P) quality. Specifically, the molecules included in the present study are HOBr, HBrO, FOBr, FBrO, BrNO, BrON, Br2O, BrBrO, BrCN, BrNC, ClOBr, ClBrO, and BrClO. Very accurate isomeric energy differences have also been determined at the CCSD(T) level with atomic natural orbital basis sets that include through g-type functions. In most cases, the isomer with a normal neutral Lewis dot structure is the lowest energy form, with the single exception that FBrO is predicted to be 11.1 kcal/mol (0 K) lower in energy than FOBr. In all cases, however, the hypervalent isomer is more stable relative to the isomer with a normal Lewis dot structure as compared to the chlorine analogs. Consistent with this observation, the energy of the last three molecules given above increases in the order ClOBr less than ClBrO less than BrClO. The CCSD(T)/TZ2P geometries and vibrational frequencies are in good agreement with the available experimental data. Heats of formation are determined for all species using a combination of theoretical isomeric, homodesmic, and isodesmic reaction energies. The accuracy of these quantities is ultimately dependent on the reliability of the experimental heat of formation of HOBr.

Lee, Timothy J.↗
DOE OSTI2015

Pressure-induced localisation of the hydrogen-bond network in KOH-VI

Using a combination of ab initio crystal structure prediction and neutron diffraction techniques, we have solved the full structure of KOH-VI at 7 GPa. Rather than being orthorhombic and proton-ordered as had previously be proposed, we find that this high-pressure phase of potassium hydroxide is tetragonal (space group I4/mmm) and proton disordered. It has an unusual hydrogen bond topology, where the hydroxyl groups form isolated hydrogen-bonded square planar (OH){sub 4} units. This structure is stable above 6.5 GPa and, despite being macroscopically proton-disordered, local ice rules enforce microscopic order of the hydrogen bonds. We suggest the use of this novel type of structure to study concerted proton tunneling in the solid state, while the topology of the hydrogen bond network could conceivably be exploited in data storage applications based solely on the manipulations of hydrogen bonds. The unusual localisation of the hydrogen bond network under applied pressure is found to be favored by a more compact packing of the constituents in a distorted cesium chloride structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗