Engineering Papers⌕ Search

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.

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

At least 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Local position-space two-nucleon potentials from leading to fourth order of chiral effective field theory

We present local, position-space chiral NN potentials through four orders of chiral effective field theory ranging from leading order (LO) to next-to-next-to-next-to-leading order (N 3 LO, fourth order) of the Δ-less version of the theory. The long-range parts of these potentials are fixed by the very accurate πN LECs as determined in the Roy-Steiner equations analysis. At the highest order (N 3 LO), the NN data below 190 MeV laboratory energy are reproduced with the respectable χ 2 /datum of 1.45. A comparison of the N 3 LO potential with the phenomenological Argonne v 18 (AV18) potential reveals substantial agreement between the two potentials in the intermediate range ruled by chiral symmetry, thus, providing a chiral underpinning for the phenomenological AV18 potential. Our chiral NN potentials may serve as a solid basis for systematic ab initio calculations of nuclear structure and reactions that allow for a comprehensive error analysis. In particular, the order by order development of the potentials will make possible a reliable determination of the truncation error at each order. Our new family of local position-space potentials differs from existing potentials of this kind by a weaker tensor force as reflected in relatively low D-state probabilities of the deuteron (P D ≲ 4.0 % for our N 3 LO potentials) and predictions for the triton binding energy above 8.00 MeV (from two-body forces alone). As a consequence, our potentials may lead to different predictions when applied to light and intermediate-mass nuclei in ab initio calculations and, potentially, help solve some of the outstanding problems in microscopic nuclear structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Utilizing Experimental-based Line Positions for Semi-empirical IR Line Lists

Molecular IR line lists computed from semi-empirically refined ab initio potential energy surface and high quality ab initio dipole moment surface may have line position accuracy typically within 0.01 - 0.05 cm−1. For high resolution spectroscopy databases, it is necessary to integrate the computed theoretical IR intensities with the line positions accurately determined from experiments and (or) reliably derived from Effective Hamiltonian models. However, it is not a trivial task, and the choices in reality are heavily contingent upon the coverage, consistency, and accuracy of the data available for a specific molecule or isotopologue. We will present several approaches applied to recent Ames IR line lists of CO 2 1, N 2 O 2 , and OCS to demonstrate how challenging such integration may become, and what we have learned From these projects. Discussions will focus on the relation between each different scenario and the corresponding choice or solution, including their advantages and limits, how the semi-empirical IR line lists can help, and what else may be needed for future improvements. We will emphasize that the “Theory+Experiment” synergy may still play significant role in the determination of the best line positions.

Effective Hamiltonian models↗
DOE OSTI2020

Ab initio leading order effective potentials for elastic nucleon-nucleus scattering

Calculating microscopic effective interactions (optical potentials) for elastic nucleon-nucleus scattering has already in the past led to a large body of work. For first-order calculations a nucleon-nucleon ( NN ) interaction and a one-body density of the nucleus were taken as input to rigorous calculations of microscopic full-folding calculations. Based on the spectator expansion of the multiple scattering series we employ a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction on the same footing in the structure as well as in the reaction calculation to obtai n an in leading-order consistent effective potential for nucleon-nucleus elastic scattering, which includes the spin of the struck target nucleon. The first order effective folding potential is computed by first deriving a nonlocal scalar density as well as a spin-projected momentum distribution. Those are then integrated with the off-shell Wolfenstein amplitudes $A$, $C$, and $M $. The resulting nonlocal potential serves as input to a momentum-space Lippmann-Schwinger equation, whose solutions are summed to obtain the nucleon-nucleus scattering observables. As a result, we calculate elastic scattering observables for $^4$He, $^6$He, $^8$He, $^{12}$C, and $^{16}$O in the energy r egime between 100 and 200~MeV projectile kinetic energy, and compare to available data. We also explore the extension down to about 70 MeV, and study the effect of ignoring the spin of the struck nucleon in the In our calculations we contrast elastic scattering off closed-shell and open-shell nuclei. We find that fo r closed-shell nuclei the approximation of ignoring the spin of the struck target nucleon is excellent. We only see effects of the spin of the struck target nucleon when considering $^6$He and $^8$He, which are nuclei with a $N/Z$ ratio larger than 1.

6 ≤ A ≤ 19↗
DOE OSTI2015

Ab initio study of the structural, electronic and optical properties of ZnTe compound

Structural, electronic and optical properties of ZnTe compound were calculated using Density Functional Theory (DFT) based on the pseudopotentials and planewaves (PP-PW) method as implemented in the ABINIT computer code, where the exchange–correlation functional is approximated using the local density approximation (LDA) and the generalized gradient approximation (GGA). The obtained results from either LDA or GGa calculation for lattice parameter, energy band gap and optical parameters, such as the fundamental absorption edge, the peaks observed in the imaginary part of the dielectric function, the macroscopic dielectric constants and the optical dielectric constant, are compared with the available theoretical results and experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

SMATool: Strength of materials analysis toolkit

The study of the strength of materials is a cornerstone in material science and engineering, playing a critical role in shaping the progress and application of materials in diverse industrial sectors. The strength of a material is meticulously examined to understand the behavior of the material under different stress conditions and environments, thereby guiding material selection and structural design. Herein, we introduce the SMATool, a computational toolkit for the efficient calculation and analysis of material strength at both zero and finite temperatures for 3D, 2D, 1D, and tubular 2D-based nanostructures and nanotubes, as well as 1D nanoribbons. The toolkit is capable of calculating tensile, shear, ultimate, yield, and indentation (Vickers' hardness) strengths in various dimensions, as well as the energy storage capacity. We conducted several calculations both at zero and finite temperatures to validate the accuracy and reliability of the developed software. Here, the results show that the SMATool package provides accurate predictions that align with existing data on material strength. SMATool integrates seamlessly with widely used electronic structure codes like VASP and Quantum Espresso, providing a user-friendly interface catering to academic researchers and industry professionals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling of Peptides with Classical and Novel Machine Learning Force Fields: A Comparison

The replacement of classical force fields (FFs) with novel neural-network-based frameworks is an emergent topic in molecular dynamics (MD) simulations. In contrast to classical FFs, which have proven their capability to provide insights into complex soft matter systems at an atomistic resolution, the machine learning (ML) potentials have yet to demonstrate their applicability for soft materials. However, the underlying philosophy, which is learning the energy of an atom in its surrounding chemical environment, makes this approach a promising tool. In particular for the exploration of novel chemical compounds, which have not been considered in the original parametrization of classical FFs. In this article, we study the performance of the ANI-2x ML model and compare the results with those of two classical FFs, namely, CHARMM27 and the GROMOS96 43a1 FF. We explore the performance of these FFs for bulk water and two model peptides, trialanine and a 9-mer of the α-aminoisobutyric acid, in vacuum and water. The results for water describe a highly ordered water structure, with a structure similar to those using ab initio molecular dynamics simulations. The energy landscape of the peptides described by Ramachandran maps show secondary structure basins similar to those of the classical FFs but differ in the position and relative stability of the basins. Details of the sampled structures show a divergent performance of the different models, which can be related either to the short-ranged nature of the ML potentials or to shortcomings of the underlying data set used for training. These findings highlight the current state of the applicability of ANI-2x ML potential for MD simulations of soft matter systems. Simultaneously, they provide insights for future improvements of current ML potentials.

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

Proton Abstraction Mediates Interactions between the Super Photobase FR0-SB and Surrounding Alcohol Solvent

Here, we report on the motional and proton transfer dynamics of the super photobase FR0-SB in the series of normal alcohols C 1 (methanol) through C 8 (n-octanol) and ethylene glycol. Steady-state and time-resolved fluorescence data reveal that the proton abstraction dynamics of excited FR0-SB depend on the identity of the solvent and that the transfer of the proton from solvent to FR0-SB*, forming FR0-HSB+*, fundamentally alters the nature of interactions between the excited molecule and its surroundings. In its unprotonated state, solvent interactions with FR0-SB* are consistent with slip limit behavior, and in its protonated form, intermolecular interactions are consistent with a much stronger interaction of FR0-HSB + * with the deprotonated solvent RO - . Furthermore, we understand the excited state population dynamics in the context of a kinetic model involving a transition state wherein FR0-HSB + * is still bound to the negatively charged alkoxide, prior to solvation of the two charged species. Data acquired in ethylene glycol confirm the hypothesis that the rotational diffusion dynamics of FR0-SB* are largely mediated by solvent viscosity while proton transfer dynamics are mediated by the lifetime of the transition state. Taken collectively, our results demonstrate that FR0-SB* extracts solvent protons efficiently and in a predictable manner, consistent with a ca. three-fold increase in dipole moment upon photoexcitation as determined by ab initio calculations based on the equation-of-motion coupled-cluster theory.

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

Impact of Modular Total Absorption Spectrometer measurements of β decay of fission products on the decay heat and reactor ν ¯ e flux calculation

Here, we report the results of a β-decay study of fission products 86 Br, 89 Kr, 89 Rb, 90 gs Rb, 90m Rb, 90 Kr, 92 Rb, 139 Xe, and 142 Cs performed with the Modular Total Absorption Spectrometer (MTAS) and on-line mass-separated ion beams. These radioactivities were assessed by the Nuclear Energy Agency as having high priority for decay heat analysis during a nuclear fuel cycle. We observe a substantial increase in β feeding to high excited states in all daughter isotopes in comparison to earlier data. This increases the average γ-ray energy emitted by the decay of fission fragments during the first 10 000 s after fission of 235 U and 239 Pu by approximately 2% and 1%, respectively, improving agreement between results of calculations and direct observations. New MTAS results reduce the reference reactor ¯ν e flux used to analyze reactor ¯ν e interaction with detector matter. The reduction determined by the ab initio method for the four nuclear fuel components, 235 U, 238 U, 239 Pu, and 241 Pu, amounts to 0.976, 0.986, 0.983, and 0.984, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Ab Initio Study of the Polarization, Electronic, Magnetic, and Optical Properties of Perovskite SrMO{sub 3} (M = Fe, Mn) Crystals and Thin Films Containing Magnetic Ions

The magnetic, electronic, and polarization properties of the SrFeO{sub 3} and SrMnO{sub 3} compounds with a perovskite structure are calculated using the density functional theory in the bulk and thin-film states. A ferroelectric instability is found to be absent in the bulk state, and the polar mode is softened in the thin-film state of SrMnO{sub 3} in the presence of tensile stresses in the substrate. As a result, a polar phase with a polarization of 23 μC/cm{sup 2} appears, which agrees with experimental data. The study of the magnetic and electronic properties demonstrates the existence of G-type antiferromagnetic ordering in SrMnO{sub 3} and the appearance of a dielectric gap of about 1.5 eV in its thin film. A ferromagnetic phase with metallic conduction in both the bulk and thin-film states is detected in SrFeO{sub 3}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗
DOE OSTI2024

Analytical ab initio hessian from a deep learning potential for transition state optimization

Identifying transition states—saddle points on the potential energy surface connecting reactant and product minima—is central to predicting kinetic barriers and understanding chemical reaction mechanisms. In this work, we train a fully differentiable equivariant neural network potential, NewtonNet, on thousands of organic reactions and derive the analytical Hessians. By reducing the computational cost by several orders of magnitude relative to the density functional theory (DFT) ab initio source, we can afford to use the learned Hessians at every step for the saddle point optimizations. We show that the full machine learned (ML) Hessian robustly finds the transition states of 240 unseen organic reactions, even when the quality of the initial guess structures are degraded, while reducing the number of optimization steps to convergence by 2–3× compared to the quasi-Newton DFT and ML methods. All data generation, NewtonNet model, and ML transition state finding methods are available in an automated workflow.

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

Structure and tunneling dynamics in a model system of peptide co-solvents: Rotational spectroscopy of the 2,2,2-trifluoroethanol⋯water complex

The hydrogen-bonding topology and tunneling dynamics of the binary adduct, 2,2,2-trifluoroethanol (TFE)⋯water, were investigated using chirped pulse and cavity based Fourier transform microwave spectroscopy with the aid of high level ab initio calculations. Rotational spectra of the most stable binary TFE⋯water conformer and five of its deuterium isotopologues were assigned. A strong preference for the insertion binding topology where water is inserted into the existing intramolecular hydrogen-bonded ring of TFE was observed. Tunneling splittings were detected in all of the measured rotational transitions of TFE⋯water. Based on the relative intensity of the two tunneling components and additional isotopic data, the splitting can be unambiguously attributed to the tunneling motion of the water subunit, i.e., the interchange of the bonded and nonbonded hydrogen atoms of water. The absence of any other splitting in the rotational transitions of all isotopologues observed indicates that the tunneling between g+ and g− TFE is quenched in the TFE⋯H{sub 2}O complex.

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

First-principles calculations of thermal, electrical, and thermoelectric transport properties of semiconductors

The transport properties of semiconductors are key to the performance of many solid-state devices (transistors, data storage, thermoelectric cooling and power generation devices, etc). An understanding of the transport details can lead to material designs with better performances. In recent years simulation tools based on first-principles calculations have been greatly improved, being able to obtain the fundamental ground-state properties of materials (such as band structure and phonon dispersion) accurately. Accordingly, methods have been developed to calculate the transport properties based on an ab initio approach. We focus on the thermal, electrical, and thermoelectric transport properties of semiconductors, which represent the basic transport characteristics of the two degrees of freedom in solids—electronic and lattice degrees of freedom. Starting from the coupled electron-phonon Boltzmann transport equations, we illustrate different scattering mechanisms that change the transport features and review the first-principles approaches that solve the transport equations. We then present the first-principles results on the thermal and electrical transport properties of semiconductors. The discussions are grouped based on different scattering mechanisms including phonon-phonon scattering, phonon scattering by equilibrium electrons, carrier scattering by equilibrium phonons, carrier scattering by polar optical phonons, scatterings due to impurities, alloying and doping, and the phonon drag effect. We show how the first-principles methods allow one to investigate transport properties with unprecedented detail and also offer new insights into the electron and phonon transport. The current status of the simulation is mentioned when appropriate and some of the future directions are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

Ames S-32 O-16 O-18 Line List for High-Resolution Experimental IR Analysis

By comparing to the most recent experimental data and spectra of the SO2 628 ν1/ν3 bands (see Ulenikov et al., JQSRT 168 (2016) 29-39), this study illustrates the reliability and accuracy of the Ames-296K SO2 line list, which is accurate enough to facilitate such high-resolution spectroscopic analysis. The SO2 628 IR line list is computed on a recently improved potential energy surface (PES) refinement, denoted Ames-Pre2, and the published purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface. Progress has been made in both energy level convergence and rovibrational quantum number assignments agreeing with laboratory analysis models. The accuracy of the computed 628 energy levels and line list is similar to what has been achieved and reported for SO2 626 and 646, i.e. 0.01-0.03 cm(exp −1) for bands up to 5500 cm(exp −1). During the comparison, we found some discrepancies in addition to overall good agreements. The three-IR-list based feature-by-feature analysis in a 0.25 cm(exp −1) spectral window clearly demonstrates the power of the current Ames line lists with new assignments, correction of some errors, and intensity contributions from varied sources including other isotopologues. We are inclined to attribute part of detected discrepancies to an incomplete experimental analysis and missing intensity in the model. With complete line position, intensity, and rovibrational quantum numbers determined at 296 K, spectroscopic analysis is significantly facilitated especially for a spectral range exhibiting such an unusually high density of lines. The computed 628 rovibrational levels and line list are accurate enough to provide alternatives for the missing bands or suspicious assignments, as well as helpful to identify these isotopologues in various celestial environments. The next step will be to revisit the SO2 828 and 646 spectral analyses.

Sulfur dioxide↗
DOE OSTI2018

Short Intramolecular O···O Contact in Some o -Dialkoxybenzene Derivatives Generates Efficient Hydrogen Bonding Acceptor Area

Statistical analysis of data extracted from the Cambridge Structural Database (CSD) has been used to investigate the crystal structure properties of odialkoxybenzene derivatives, compounds containing two ether oxygen acceptors in ortho positions of benzene ring. It has been shown that in more than 90% of cases, the fragment has predictable geometrical characteristics where the two ether oxygens form short interatomic O center dot center dot center dot O contact (2.57 A in average), while O-substitutents take trans position, both approximately coplanar with the benzene ring. Such arrangement of oxygen acceptors produces a large and uniform area of the negative electrostatic potential suitable for multiple hydrogen bonding. The acceptor abilities of the O center dot center dot center dot O system have been investigated by the statistical CSD analysis. The ab initio estimation of the interaction energy in the dimer of o-dimethoxybenzene (DMB) and H2O, employed as a model system, is achieved via high-level electron correlation CCSD(T) calculation with the complete basis set extrapolation. The interaction energy is estimated to be 6.5 kcal/mol. The results indicate the existence of a very flat potential the region between methoxy oxygens and that DMB water is a highly flexible system. The structural role of the O center dot center dot center dot O acceptor system is particularly interesting considering its ability to form multiple hydrogen bonding.

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

Ab initio study of intrinsic point defects in PbTe: an insight into phase stability

The stability of intrinsic point defects in PbTe, one of the most widely studied and efficient thermoelectric material, is explored by means of Density Functional Theory (DFT). The origin of n- and p-type conductivity in PbTe is attributed to particular intrinsic charged defects by calculating their formation energies. Here, these DFT calculated defect formation energies are then used in the Gibbs free energy description of this phase as part of the Pb-Te thermodynamic model built using the CALPHAD method, and in the resulting phase diagram it is found that its solubility lines and non-stoichiometric range agree very well with experimental data. Such an approach of using DFT in conjunction with CALPHAD for compound semiconductor phases that exhibit very small ranges of non-stoichiometry does not only make the process of calculating phase diagrams for such systems more physical, but is necessary and critical for the assessment of unknown phase diagrams.

36 MATERIALS SCIENCE↗