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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Lattice stability and high-pressure melting mechanism of dense hydrogen up to 1.5 TPa

Lattice stability and metastability, as well as melting, are vital features of the physics and chemistry of dense hydrogen. Using ab initio molecular dynamics (AIMD), the classical superheating limit and melting line of metallic hydrogen are investigated up to 1.5 TPa. The computations show that the classical superheating degree is about 100 K, and the classical melting curve becomes flat at a level of 350 K when beyond 500 GPa. This data allows us to estimate the well depth and the potential barriers that must be overcome when the crystal melts. Inclusion of nuclear quantum effects (NQE) using path integral molecular dynamics (PIMD) predicts that both superheating limit and melting temperature are lowered to below room temperature, but the latter never reaches absolute zero. Detailed analysis indicates that the melting is thermally activated, rather than driven by pure zero-point motion (ZPM). This argument was further supported by extensive PIMD simulations, demonstrating the stability of Fddd structure against liquefaction at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 ·nH 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Finally, water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

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

Pressure-dependent kinetics of cyclopentene consumption

Recent experiments have uncovered deficiencies in the ability of current kinetic models for cyclopentane combustion to predict species profiles as a function of O2 concentration. One likely source of these discrepancies flows from the limited availability of relevant kinetic data for key functionalized intermediates, requiring the use of distant analogies to estimate their rates. This work fills this gap for cyclopentene, the most prominent C5 intermediate in low-temperature cyclopentane combustion. Ab initio transition-state-theory based master equation methods are used to calculate pressure-dependent rate constants across four potential energy surfaces (C5H8 + OH, C5H8 + HO2, C5H7, and C5H7 + O2). Potential energy surfaces are characterized at the CCSD(T)-F12/cc-pVDZ-F12//revDSD-PBEP86-D3BJ/def2-TZVPP level of theory. Pressure-dependent rate constants and branching fractions are determined via 1D master equation calculations. The results are used to analyze the competition for radicals between cyclopentene and cyclopentane in the intermediate temperature region, and to identify the dominant pathways flowing from these initiation channels. On the OH initiation surface, we find good agreement with recently published experimental rate constants at high temperature, and we elucidate the complex temperature dependence below 800 K, where 𝜋𝜋 addition becomes prominent. On the HO2 initiation surface, we find an interesting alternative route to bicyclic ether formation via well-skipping over the HO2 adduct well. On the cyclopentenyl radical surface, the major product is cyclopentadiene plus H. On the cyclopentenyl + O2 surface, we find that resonance stabilization largely prevents the allylic radical from undergoing oxidation, whereas the alkylic radical has a somewhat elevated propensity for oxidation pathways relative to cyclopentyl.

ab initio calculations↗
DOE OSTI2015

Knockout beyond the dripline

In this contribution, we will describe neutron and proton removal from {sup 9}C and {sup 7}Be which are two particularly interesting nuclei entering the nucleo-synthesis pp-chain [1, 2]. Neutron and proton removal reactions have been used in the past twenty years to probe the single-particle structure of exotic nuclei. The core parallel-momentum distribution can give information on the angular momentum and spin of the nucleon initial state while the total removal cross section is sensitive to the asymptotic part of the initial wave function and also to the reaction mechanism. Because knockout is a peripheral reaction from which the Asymptotic Normalization Constant (ANC) of the single-particle wave function can be extracted, it has been used as an indirect method to obtain the rate of reactions like {sup 8}B(p,γ){sup 9}C or {sup 7}Be(p,γ){sup 8}B. Nucleon removal has recently been applied by the HiRA collaboration [3] to situations in which the remaining “core” is beyond the drip line, such as {sup 8}C and {sup 6}Be, unbound by one or more protons, and whose excitation-energy spectrum can be obtained by the invariant-mass method. By gating on the ground-state peak, “core” parallel-momentum distributions and total knockout cross sections have been obtained similar to previous studies with well-bound “cores”. In addition for each projectile, knock out to final bound states has also been obtained in several cases. We will report on the theoretical description and comparison to this experimental data for a few cases for which advances in the accuracy of the transfer-to-the continuum model [4, 5] have been made [6]. These include the use, when available, of “ab-initio” overlaps for the initial state [7] and in particular their ANC values [8]. Also, the construction of a nucleus-target folding potential for the treatment of the core-target S-matrix [9] using for the cores “ab-initio” densities [10] and state-of-the-art n−{sup 9}Be optical potentials [11]. Preliminary results and open problems will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Accurate nuclear radii and binding energies from a chiral interaction

With the goal of developing predictive ab initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO sat , yield accurate binding energies and radii of nuclei up to 40 Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective J π =3 - states in 16 O and 40 Ca are described accurately, while spectra for selected p- and sd-shell nuclei are in reasonable agreement with experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Benchmarking the performance of uncertainty quantification methods for neural network-based interatomic potentials

Machine-learned interatomic potentials (ML-IAPs) continue to gain popularity as accurate, computationally efficient replacements for traditional, physics-based interatomic potentials and expensive ab initio methods. Uncertainty quantification (UQ) of ML-IAPs is a growing area of research as UQ is critical in many applications of IAPs, such as developing curated datasets, active learning-based data augmentation, self-improving models, and estimating the uncertainty of molecular dynamics simulations. In this paper, we construct and benchmark a series of different neural network potentials (NNPs) with varying network architectures to determine the performance of these models with respect to both the mean and uncertainty calibration error. Each NNP method is specifically designed to predict either epistemic or aleatoric uncertainty with particular focus on the differences in behavior between the epistemic and aleatoric uncertainty estimates. We benchmark these methods using multiple datasets common in the ML-IAP literature. The results show that the aleatoric uncertainty from single-shot model architectures is a competitive alternative to ensemble-based epistemic uncertainty predictions in regions of sufficient data-density. However, in regions where the representative data is sparse, aleatoric uncertainty models tend to overpredict and epistemic methods tend to underpredict the actual model error. We conclude that the type of UQ is crucial when discussing performance of probabilistic model results as different methods have different performance characteristics depending on the regime in which they are evaluated. Therefore, the type of UQ method should be carefully evaluated against both the data characteristics and requirements for the intended application.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Intermolecular potential for thermal H2O-He collisions

Theoretical potentials for rotational excitation of H2O by He were constructed via several methods, all of which start with a large basis set SCF interaction. The semiempirical Hartree-Fock with damped dispersion model adds a damped long-range attraction with parameters adjusted to fit experimental total differential cross sections. Purely ab initio potentials add correlation energies obtained via perturbation theory (MP2 and MP4) or a variational method (ICF1). Scattering calculations were performed on all surfaces to compare wih available beam scattering and pressure broadening data and to assess sensitivity of state-to-state rates to uncertainties in the potential. From comparison with the limited experimental data, the ICF1 surface appears to be marginally better than the MP4 surface. Thermal rates calculated from this surface should be accurate to better than 50 percent, at least for the larger, more important rates.

Palma, Amedeo↗
DOE OSTI2016

CO Oxidation and Subsequent CO 2 Chemisorption on Alkaline Zirconates: Li 2 ZrO 3 and Na 2 ZrO 3

Here, two different alkaline zirconates (Li 2 ZrO 3 and Na 2 ZrO 3 ) were studied as possible bifunctional catalytic-captor materials for CO oxidation and the subsequent CO 2 chemisorption process. Initially, CO oxidation reactions were analyzed in a catalytic reactor coupled to a gas chromatograph, using Li 2 ZrO 3 and Na 2 ZrO 3 , under different O 2 partial flows. We found results clearly showed that Na 2 ZrO 3 possesses much better catalytic properties than Li 2 ZrO 3 . After the CO-O 2 oxidation catalytic analysis, CO2 chemisorption process was analyzed by thermogravimetric analysis, only for the Na 2 ZrO 3 ceramic. The results confirmed that Na 2 ZrO 3 is able to work as a bifunctional material (CO oxidation and subsequent CO 2 chemisorption), although the kinetic CO 2 capture process was not the best one under the physicochemical condition used in this case. For Na 2 ZrO 3 , the best CO conversions were found between 445 and 580 °C (100%), while Li 2 ZrO 3 only showed a 35% of efficiency between 460 and 503 °C. However, in the Na 2 ZrO 3 case, at temperatures higher than 580 °C its catalytic activity gradually decreases as a result of CO 2 capture process. Finally, all these experiments were compared and supported with theoretical thermodynamic data.

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

Geometry and electronic structure of iridium adsorbed on graphene

We report investigation of the geometry and electronic structure of iridium atoms adsorbed onto graphene through a combined experimental and theoretical study. Ir atoms were deposited onto a flake of graphene on a Pt(111) surface and found to form clusters even at low temperatures. The areal density of the observed clusters on the graphene flake suggests the clusters are most likely pairs of Ir atoms. Theoretical ab initio density functional (DFT) calculations indicate that these Ir dimers are oriented horizontally, near neighboring "bridge" sites of the graphene lattice, as this configuration has the strongest adsorption energy of all high-symmetry configurations for the Ir dimer. A large peak in the local density of states (LDOS) at the Dirac point energy was measured via scanning tunneling spectroscopy, and this result is reproduced by a DFT calculation of the LDOS. The peak at the Dirac point energy is found to be from the Ir s and p states. Finally, the LDOS in the monomer case was also calculated, and is found to significantly differ from the experimentally determined data, further supporting the hypothesis of low-temperature clustering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1995

A correlated ab initio study of linear carbon-chain radicals CnH (n = 2-7)

Linear carbon-chain radicals CnH for n = 2-7 have been studied with correlation consistent valence and core-valence basis sets and the coupled cluster method RCCSD(T). Equilibrium structures, rotational constants, and dipole moments are reported and compared with available experimental data. The ground state of the even-n series changes from 2 sigma+ to 2 pi as the chain is extended. For C4H, the 2 sigma+ state was found to lie only 72 cm-1 below the 2 pi state in the estimated complete basis set limit for valence correlation. The C2H- and C3H- anions have also been characterized.

NASA Discipline Exobiology↗
DOE OSTI2018

Simulating the density of organic species in the atmosphere of Titan with a coupled ion-neutral photochemical model

We introduce a one-dimensional coupled ion-neutral photochemical kinetics and diffusion model to study the atmospheric composition of Titan in light of new theoretical kinetics calculations and scientific findings from the Cassini–Huygens mission. The model extends from the surface to the exobase. The atmospheric background, boundary conditions, vertical transport and aerosol opacity are all constrained by the Cassini–Huygens observations. The chemical network includes reactions between hydrocarbons, nitrogen and oxygen bearing species. It takes into account neutrals and both positive and negative ions with masses extending up to 116 and 74 u, respectively. We incorporate high-resolution isotopic photoabsorption and photodissociation cross sections for N 2 as well as new photodissociation branching ratios for CH 4 and C 2 H 2 . Ab initio transition state theory calculations are performed in order to estimate the rate coefficients and products for critical reactions. Critical reactions of production and loss for neutrals and ions are quantitatively assessed and thoroughly discussed. The vertical distributions of neutrals and ions predicted by the model generally reproduce observational data, suggesting that for the small species most chemical processes in Titan’s atmosphere and ionosphere are adequately described and understood; some differences are highlighted. Notable remaining issues include (i) the total positive ion density (essentially HCNH + ) in the upper ionosphere, (ii) the low mass negative ion densities (CN - , C 3 N - /C4H - ) in the upper atmosphere, and (iii) the minor oxygen-bearing species (CO 2 , H 2 O) density in the stratosphere. Pathways towards complex molecules and the impact of aerosols (UV shielding, atomic and molecular hydrogen budget, nitriles heterogeneous chemistry and condensation) are evaluated in the model, along with lifetimes and solar cycle variations.

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

Density-functional approach to the three-body dispersion interaction based on the exchange dipole moment

We implement and compute the density functional nonadditive three-body dispersion interaction using a combination of Tang-Karplus formalism and the exchange-dipole moment model of Becke and Johnson. The computation of the C{sub 9} dispersion coefficients is done in a non-empirical fashion. The obtained C{sub 9} values of a series of noble atom triplets agree well with highly accurate values in the literature. We also calculate the C{sub 9} values for a series of benzene trimers and find a good agreement with high-level ab initio values reported recently in the literature. For the question of damping of the three-body dispersion at short distances, we propose two damping schemes and optimize them based on the benzene trimers data, and the fitted analytic potentials of He{sub 3} and Ar{sub 3} trimers fitted to the results of high-level wavefunction theories available from the literature. Both damping schemes respond well to the optimization of two parameters.

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

Low temperature rate coefficients of the H + CH{sup +} → C{sup +} + H{sub 2} reaction: New potential energy surface and time-independent quantum scattering

The observed abundances of the methylidyne cation, CH{sup +}, in diffuse molecular clouds can be two orders of magnitude higher than the prediction of the standard gas-phase models which, in turn, predict rather well the abundances of neutral CH. It is therefore necessary to investigate all the possible formation and destruction processes of CH{sup +} in the interstellar medium with the most abundant species H, H{sub 2}, and e{sup −}. In this work, we address the destruction process of CH{sup +} by hydrogen abstraction. We report a new calculation of the low temperature rate coefficients for the abstraction reaction, using accurate time-independent quantum scattering and a new high-level ab initio global potential energy surface including a realistic model of the long-range interaction between the reactants H and CH{sup +}. The calculated thermal rate coefficient is in good agreement with the experimental data in the range 50 K–800 K. However, at lower temperatures, the experimental rate coefficient takes exceedingly small values which are not reproduced by the calculated rate coefficient. Instead, the latter rate coefficient is close to the one given by the Langevin capture model, as expected for a reaction involving an ion and a neutral species. Several recent theoretical works have reported a seemingly good agreement with the experiment below 50 K, but an analysis of these works show that they are based on potential energy surfaces with incorrect long-range behavior. The experimental results were explained by a loss of reactivity of the lowest rotational states of the reactant; however, the quantum scattering calculations show the opposite, namely, a reactivity enhancement with rotational excitation.

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

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1996

An Ab Initio Based Potential Energy Surface for Water

We report a new determination of the water potential energy surface. A high quality ab initio potential energy surface (PES) and dipole moment function of water have been computed. This PES is empirically adjusted to improve the agreement between the computed line positions and those from the HITRAN 92 data base. The adjustment is small, nonetheless including an estimate of core (oxygen 1s) electron correlation greatly improves the agreement with experiment. Of the 27,245 assigned transitions in the HITRAN 92 data base for H2(O-16), the overall root mean square (rms) deviation between the computed and observed line positions is 0.125/cm. However the deviations do not correspond to a normal distribution: 69% of the lines have errors less than 0.05/cm. Overall, the agreement between the line intensities computed in the present work and those contained in the data base is quite good, however there are a significant number of line strengths which differ greatly.

Partridge, Harry↗
DOE OSTI2015

Ab initio calculations of three-nucleon scattering

Results on three-nucleon (3N) elastic scattering below the pion production threshold are discussed with an emphasis on the need for a three-nucleon force (3NF). The large discrepancies found between a theory based on numerical solutions of 3N Faddeev equations with (semi)phenomenological NN potentials only and data point to the need for 3NF’s. This notion is supported by the fact that another possible reason for the discrepancies in elastic nucleon-deuteron (Nd) scattering, relativistic effects, turned out to be small. Results for new generation of chiral NN forces (up to N{sup 4}LO) together with theoretical truncation errors are shown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Modeling of molecular nitrogen collisions and dissociation processes for direct simulation Monte Carlo

The Direct Simulation Monte Carlo (DSMC) method typically used for simulating hypersonic Earth re-entry flows requires accurate total collision cross sections and reaction probabilities. However, total cross sections are often determined from extrapolations of relatively low-temperature viscosity data, so their reliability is unknown for the high temperatures observed in hypersonic flows. Existing DSMC reaction models accurately reproduce experimental equilibrium reaction rates, but the applicability of these rates to the strong thermal nonequilibrium observed in hypersonic shocks is unknown. For hypersonic flows, these modeling issues are particularly relevant for nitrogen, the dominant species of air. To rectify this deficiency, the Molecular Dynamics/Quasi-Classical Trajectories (MD/QCT) method is used to accurately compute collision and reaction cross sections for the N{sub 2}({sup 1}Σ{sub g}{sup +})-N{sub 2}({sup 1}Σ{sub g}{sup +}) collision pair for conditions expected in hypersonic shocks using a new potential energy surface developed using a ReaxFF fit to recent advanced ab initio calculations. The MD/QCT-computed reaction probabilities were found to exhibit better physical behavior and predict less dissociation than the baseline total collision energy reaction model for strong nonequilibrium conditions expected in a shock. The MD/QCT reaction model compared well with computed equilibrium reaction rates and shock-tube data. In addition, the MD/QCT-computed total cross sections were found to agree well with established variable hard sphere total cross sections.

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

Dissection of the multichannel reaction of acetylene with atomic oxygen: from the global potential energy surface to rate coefficients and branching dynamics

The O( 3 P) + C 2 H 2 reaction is the first step in acetylene oxidation. The accurate kinetic data and the understanding of the reaction dynamics is of great importance. To this end, in this study a full-dimensional global potential energy surface (PES) for the ground triplet state of the O( 3 P) + C 2 H 2 reaction is constructed based on approximately 85,000 ab initio points calculated at the level of explicitly correlated unrestricted coupled cluster single, double, and perturbative triple excitations with the explicitly correlated polarized valence triple zeta basis set (UCCSD(T)-F12b/VTZ-F12). The PES is fit using the permutation invariant polynomial-neural network (PIP-NN) approach with a total root mean square error of 0.21 kcal mol –1 . The key topographic features of the PES, including multiple potential wells and saddle points along different reaction pathways, are well represented by this fit PES. The kinetics and dynamics of the O( 3 P) + C 2 H 2 reaction are investigated using the quasi-classical trajectory (QCT) method. The calculated rate coefficients are in good agreement with experimental data over a wide temperature range, especially when the temperature is lower than 1500 K. The product branch ratio has also been determined, which indicates the H + HCCO channel as the dominant reaction pathway at 298–3000 K, accounting for 80–90% of the overall rate coefficient, in agreement with experimental observations. The dynamics of the reaction is analyzed in detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗