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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Theoretical Studies of the Quasi-Bound Vibrational Levels in the B State of CH3

Consider the calculation of the oscillator strength of th B yielded by X electronic transition in CH3. The X state is strongly bound, but in the B state, CH3 is pre-dissociated by H atom tunneling through a potential energy barrier to form CH2+H. If we turn off pre-dissociation by making the potential outside of the barrier repulsive, the calculated 0-0 oscillator strength for CD3 from high quality ab initio potential and transition moment data along with accurate ro-vibrational wavefunctions is 0.0212, which is a factor of two greater than experiment. This discrepancy is largely due to the neglect of pre-dissociation. In this poster we explore techniques for including pre-dissociation in the calculations.

Schwenke, David W.↗
DOE OSTI2024

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Molecular dynamics simulation of Fe-Si alloys using a neural network machine learning potential

Interatomic potential development using machine learning (ML) approaches has attracted a lot of attention in recent years because these potentials can effectively describe the structural and dynamical properties of complex materials at the atomistic level. Here, in this work, we present the development of a neural network (NN) deep ML interatomic potential for Fe-Si alloys, and we demonstrate the effectiveness of the NN-ML potential in predicting the structures and energies of liquid and crystalline phases of Fe-Si alloys in comparison with the results from ab initio molecular dynamics simulations or experimental data. The developed NN-ML potential is also used to perform molecular dynamics simulations to study the structures of Fe-Si alloys with various compositions under rapid solidification conditions. The short-ranged orders in the rapidly solidified Fe-Si alloys are also analyzed by a cluster alignment method.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Potential Energy Curves and Transport Properties for the Interaction of He with Other Ground-state Atoms

The interactions of a He atom with a heavier atom are examined for 26 different elements, which are consecutive members selected from three rows (Li - Ne, Na - Ar, and K,Ca, Ga - Kr) and column 12 (Zn,Cd) of the periodic table. Interaction energies are determined wing high-quality ab initio calculations for the states of the molecule that would be formed from each pair of atoms in their ground states. Potential energies are tabulated for a broad range of Interatomic separation distances. The results show, for example, that the energy of an alkali interaction at small separations is nearly the same as that of a rare-gas interaction with the same electron configuration for the dosed shells. Furthermore, the repulsive-range parameter for this region is very short compared to its length for the repulsion dominated by the alkali-valence electron at large separations (beyond about 3-4 a(sub 0)). The potential energies in the region of the van der Waals minimum agree well with the most accurate results available. The ab initio energies are applied to calculate scattering cross sections and obtain the collision integrals that are needed to determine transport properties to second order. The theoretical values of Li-He total scattering cross sections and the rare-gas atom-He transport properties agree well (to within about 1%) with the corresponding measured data. Effective potential energies are constructed from the ab initio energies; the results have been shown to reproduce known transport data and can be readily applied to predict unknown transport properties for like-atom interactions.

Partridge, Harry↗
DOE OSTI2024

Second-harmonic generation tensors from high-throughput density-functional perturbation theory

Optical materials play a key role in enabling modern optoelectronic technologies in a wide variety of domains such as the medical or the energy sector. Among them, nonlinear optical crystals are of primary importance to achieve a broader range of electromagnetic waves in the devices. However, numerous and contradicting requirements significantly limit the discovery of new potential candidates, which, in turn, hinders the technological development. In the present work, the static nonlinear susceptibility and dielectric tensor are computed via density-functional perturbation theory for a set of 579 inorganic semiconductors. The computational methodology is discussed and the provided database is described with respect to both its data distribution and its format. Several comparisons with both experimental and ab initio results from literature allow to confirm the reliability of our data. The aim of this work is to provide a relevant dataset to foster the identification of promising nonlinear optical crystals in order to motivate their subsequent experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

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

Toward Realistic Models of the Electrocatalytic Oxygen Evolution Reaction

The electrocatalytic oxygen evolution reaction (OER) supplies the protons and electrons needed to transform renewable electricity into chemicals and fuels. However, the OER is kinetically sluggish; it operates at significant rates only when the applied potential far exceeds the reversible voltage. The origin of this overpotential is hidden in a complex mechanism involving multiple electron transfers and chemical bond making/breaking steps. Our desire to improve catalytic performance has then made mechanistic studies of the OER an area of major scientific inquiry, though the complexity of the reaction has made understanding difficult. While historically, mechanistic studies have relied solely on experiment and phenomenological models, over the past twenty years ab initio simulation has been playing an increasingly important role in developing our understanding of the electrocatalytic OER and its reaction mechanisms. In this Review we cover advances in our mechanistic understanding of the OER, organized by increasing complexity in the way through which the OER is modeled. We begin with phenomenological models built using experimental data before reviewing early efforts to incorporate ab initio methods into mechanistic studies. We go on to cover how the assumptions in these early ab initio simulations-no electric field, electrolyte, or explicit kinetics-have been relaxed. Through comparison with experimental literature, we explore the veracity of these different assumptions. We summarize by discussing the most critical open challenges in developing models to understand the mechanisms of the OER.

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

Directed gas phase preparation of ethynylallene (H 2 CCCHCCH; X 1 A') via the crossed molecular beam reaction of the methylidyne radical (CH; X 2 Π) with vinylacetylene (H 2 CCHCCH; X 1 A')

The gas-phase bimolecular reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A') was conducted at a collision energy of 20.3 kJ mol –1 under single collision conditions exploiting the crossed molecular beam experimental results merged with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations. The laboratory data reveal that the bimolecular reaction proceeds barrierlessly via indirect scattering dynamics through long-lived C 5 H 5 reaction intermediate(s) ultimately dissociating to C 5 H 4 isomers along with atomic hydrogen with the latter predominantly originating from the vinylacetylene reactant as confirmed by the isotopic substitution experiments in the D1-methylidyne–vinylacetylene reaction. Combined with ab initio calculations of the potential energy surface (PES) and statistical Rice–Ramsperger–Kassel–Marcus (RRKM) calculations, the experimental determined reaction energy of –146 ± 26 kJ mol –1 along with the distribution minimum of T(θ) at 90° and isotopic substitution experiments suggest ethynylallene (p1; ΔrG = –230 ± 4 kJ mol –1 ) as the dominant product. The ethynylallene (p1) may be formed with extensive rovibrational excitation, which would result in a lower maximum translational energy. Further, AIMD simulations reveal that the reaction dynamics leads to p1 (ethynylallene, 75%) plus atomic hydrogen with the dominant initial complex being i1 formed by methylidyne radical addition to the double C[double bond, length as m-dash]C bond in vinylacetylene. Altogether, combining the crossed molecular beam experimental results with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations, ethynylallene (p1) is expected to represent the dominant product in the reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A').

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

Ab Initio Computation of the Longitudinal Response Function in 40 Ca

In this work, we present a consistent ab initio computation of the longitudinal response function RL in 40 Ca using the coupled-cluster and Lorentz integral transform methods starting from chiral nucleon-nucleon and three-nucleon interactions. We validate our approach by comparing our results for R L in 4 He and the Coulomb sum rule in 40 Ca against experimental data and other calculations. For R L in 40 Ca we obtain a very good agreement with experiment in the quasielastic peak up to intermediate momentum transfers, and we find that final state interactions are essential for an accurate description of the data. This work presents a milestone towards ab initio computations of neutrino-nucleus cross sections relevant for experimental long-baseline neutrino programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗
DOE OSTI2017

Bayesian Estimation of Thermonuclear Reaction Rates for Deuterium+Deuterium Reactions

The study of d+d reactions is of major interest since their reaction rates affect the predicted abundances of D, 3 He, and 7 Li. In particular, recent measurements of primordial D/H ratios call for reduced uncertainties in the theoretical abundances predicted by Big Bang nucleosynthesis (BBN). Different authors have studied reactions involved in BBN by incorporating new experimental data and a careful treatment of systematic and probabilistic uncertainties. To analyze the experimental data, Coc et al. used results of ab initio models for the theoretical calculation of the energy dependence of S -factors in conjunction with traditional statistical methods based on χ 2 minimization. Bayesian methods have now spread to many scientific fields and provide numerous advantages in data analysis. Astrophysical S -factors and reaction rates using Bayesian statistics were calculated by Iliadis et al. In this paper we present a similar analysis for two d+d reactions, d(d, n) 3 He and d(d, p) 3 H, that has been translated into a total decrease of the predicted D/H value by 0.16%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

HEPOM: Using Graph Neural Networks for the Accelerated Predictions of Hydrolysis Free Energies in Different pH Conditions

Hydrolysis is a fundamental family of chemical reactions where water facilitates the cleavage of bonds. The process is ubiquitous in biological and chemical systems, owing to water’s remarkable versatility as a solvent. However, accurately predicting the feasibility of hydrolysis through computational techniques is a difficult task, as subtle changes in reactant structure like heteroatom substitutions or neighboring functional groups can influence the reaction outcome. Furthermore, hydrolysis is sensitive to the pH of the aqueous medium, and the same reaction can have different reaction properties at different pH conditions. In this work, we have combined reaction templates and high-throughput ab initio calculations to construct a diverse data set of hydrolysis free energies. The developed framework automatically identifies reaction centers, generates hydrolysis products, and utilizes a trained graph neural network (GNN) model to predict ΔG values for all potential hydrolysis reactions in a given molecule. The long-term goal of the work is to develop a data-driven, computational tool for high-throughput screening of pH-specific hydrolytic stability and the rapid prediction of reaction products, which can then be applied in a wide array of applications including chemical recycling of polymers and ion-conducting membranes for clean energy generation and storage.

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

Chemical crystallography by serial femtosecond X-ray diffraction

Abstract Inorganic–organic hybrid materials represent a large share of newly reported structures, owing to their simple synthetic routes and customizable properties 1 . This proliferation has led to a characterization bottleneck: many hybrid materials are obligate microcrystals with low symmetry and severe radiation sensitivity, interfering with the standard techniques of single-crystal X-ray diffraction 2,3 and electron microdiffraction 4–11 . Here we demonstrate small-molecule serial femtosecond X-ray crystallography (smSFX) for the determination of material crystal structures from microcrystals. We subjected microcrystalline suspensions to X-ray free-electron laser radiation 12,13 and obtained thousands of randomly oriented diffraction patterns. We determined unit cells by aggregating spot-finding results into high-resolution powder diffractograms. After indexing the sparse serial patterns by a graph theory approach 14 , the resulting datasets can be solved and refined using standard tools for single-crystal diffraction data 15–17 . We describe the ab initio structure solutions of mithrene (AgSePh) 18–20 , thiorene (AgSPh) and tethrene (AgTePh), of which the latter two were previously unknown structures. In thiorene, we identify a geometric change in the silver–silver bonding network that is linked to its divergent optoelectronic properties 20 . We demonstrate that smSFX can be applied as a general technique for structure determination of beam-sensitive microcrystalline materials at near-ambient temperature and pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Homochiral and racemic MicroED structures of a peptide repeat from the ice-nucleation protein InaZ

The ice-nucleation protein InaZ from Pseudomonas syringae contains a large number of degenerate repeats that span more than a quarter of its sequence and include the segment GSTSTA. Ab initio structures of this repeat segment, resolved to 1.1 Å by microfocus X-ray crystallography and to 0.9 Å by the cryo-EM method MicroED, were determined from both racemic and homochiral crystals. The benefits of racemic protein crystals for structure determination by MicroED were evaluated and it was confirmed that the phase restriction introduced by crystal centrosymmetry increases the number of successful trials during the ab initio phasing of the electron diffraction data. Both homochiral and racemic GSTSTA form amyloid-like protofibrils with labile, corrugated antiparallel β-sheets that mate face to back. The racemic GSTSTA protofibril represents a new class of amyloid assembly in which all-left-handed sheets mate with their all-right-handed counterparts. This determination of racemic amyloid assemblies by MicroED reveals complex amyloid architectures and illustrates the racemic advantage in macromolecular crystallography, now with submicrometre-sized crystals.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Is interstellar detection of higher members of the linear radicals CnCH and CnN feasible?

Rotational constants and dipole moments for linear-chain radicals CnCH and CnN are estimated using a combinatiaon of ab initio molecular orbital calculations and observed data on the starting members of the series. CnCH with n = 0-5 have been observed by radioastronomy in carbon-rich interstellar clouds; higher members of the series have 2Pi ground states with large dipole moments and are strong candidates for observation. CN and C3N have also been observed by radioastronomy; higher members of the series, with the possible exception of C5N, have 2Pi ground states with near-zero dipole moments making their interstellar detection hopeless under present observational conditions. C5N can be a strong candidate only if it has a 2Sigma ground state, and best computations so far indicate that this is not the case.

Pauzat, F.↗
NASA NTRS2015

Excited Vibrational Level Rotational Constants for SiC2: A Sensitive Molecular Diagnostic for Astrophysical Conditions

Silacyclopropynylidene, SiC2, is a known and highly abundant circumstellar molecule. Its spectrum has been established as a major component of lines observed toward the carbon-rich star IRC +10216 (CW Leonis). It has been detected in its low-lying v(sub 3) = 1 and 2 vibrational states as well as in various isotopic compositions. Increasing sensitivity and spatial resolution will enable many more emission or absorption lines to be detected. In order to detect new molecular species, unassigned lines of known species must be identified. This work uses established ab initio quartic force fields to produce data necessary for this classification of lines related to SiC2. Agreement between the theoretical vibrational frequencies and known rotational and spectroscopic constants is quite good, as good as 5 cm(exp -1) and 3 MHz, respectively in some cases. Additionally, experimentally unknown vibrational frequencies and rotational constants are provided for the first overtones and combination bands in addition to 3(sub v3), the second overtone of the low-lying antisymmetric stretch/carbide rotation mode. Frequencies of v(sub 3) = 3 low-J rotational transitions of the main isotopic species are also estimated from published data for v(sub 3) ≤ 2. Further, we determine rotational and centrifugal distortion parameters for which in most cases vibrational effects due to the v(sub 3) mode were reduced to first, and in several cases also to second order. These values may approximate equilibrium values better than the ground state values. The data produced herein will aid in the experimental and observational characterization of this known astromolecule in order to identify some of the unassigned lines for a known entity.

Astrochemistry↗
DOE OSTI2022

Precision measurement of Compton scattering in silicon with a skipper CCD for dark matter detection

Experiments aiming to directly detect dark matter through particle recoils can achieve energy thresholds of $\mathcal{O}(1\,\mathrm{eV})$. In this regime, ionization signals from small-angle Compton scatters of environmental $\gamma$-rays constitute a significant background. Monte Carlo simulations used to build background models have not been experimentally validated at these low energies. We report a precision measurement of Compton scattering on silicon atomic shell electrons down to 23$\,$eV. A skipper charge-coupled device (CCD) with single-electron resolution, developed for the DAMIC-M experiment, was exposed to a $^{241}$Am $\gamma$-ray source over several months. Features associated with the silicon K, L$_{1}$, and L$_{2,3}$-shells are clearly identified, and scattering on valence electrons is detected for the first time below 100$\,$eV. We find that the relativistic impulse approximation for Compton scattering, which is implemented in Monte Carlo simulations commonly used by direct detection experiments, does not reproduce the measured spectrum below 0.5$\,$keV. The data are in better agreement with $ab$$initio$ calculations originally developed for X-ray absorption spectroscopy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Ab Initio Molecular Dynamics Simulation of Divalent Metal Cation Incorporation in Calcite: Implications for Interpreting X-ray Absorption Spectroscopy Data

Calcite (CaCO 3 ) is a ubiquitous mineral with the ability to accommodate a wide range of impurities. Determination of the coordination environment and incorporation modes of impurities in calcite has historically relied primarily on the interpretation of extended X-ray absorption fine structure (EXAFS) spectroscopy. However, lack of standards combined with the large number of degrees of freedom involved in shell-by-shell fits have made the interpretation of EXAFS spectra challenging. As such, in this work, ab initio molecular dynamics (AIMD) simulations were performed to investigate the incorporation of seven divalent metal cation impurities, namely, Ba 2+ , Pb 2+ , Sr 2+ , Cd 2+ , Mn 2+ , Co 2+ , and Zn 2+ , in calcite. These cations span a wide range of sizes: 30.7 to –22.8% change in ionic radius with respect to Ca 2+ . The atomic trajectories were then used to compute EXAFS spectra for direct comparison with published experimental spectra. The simulations confirmed that all seven metal cations incorporate in calcite via substitution in the calcium site. The AIMD–EXAFS approach allowed for overcoming limitations of the shell-by-shell fitting approach, such as the difficulties in resolving weak backscatterers beyond the first coordination shell. As a result, the AIMD–EXAFS approach was able to provide a detailed and comprehensive characterization of the structural relaxation around divalent metal cation impurities in calcite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗