Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio”

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

The ab-initio density matrix renormalization group in practice

The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.

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

Comparison of alpha and beta tin for lithium, sodium, and magnesium storage: An ab initio study including phonon contributions

We present a comparative ab initio study of Li, Na, and Mg storage in tin, including phononic effects and phase competition between α and β Sn. Mg doping at low concentration is found to stabilize the β phase. On the contrary, Li and Na doping is shown to reverse the stability of the phases at room temperature: Li/Na-doped α-Sn is more stable than Li/Na-doped β-Sn up to a temperature of around 380/400 K. This may rationalize the formation of α-Sn upon lithiation and delithiation of β-Sn anodes reported in experimental studies. The changes in phase stability with Li/Na/Mg doping are directly related to the intercalation energies of Li/Na/Mg in one phase versus the other: at 300 K, Li/Na is easier intercalated in α-Sn (−0.37/−0.08 eV) than in β-Sn (0.06/0.49 eV), while Mg intercalation energy is, although positive (i.e., unfavored intercalation), lower in β-Sn (0.53 eV) than in α-Sn (0.66 eV). The temperature effect is found to affect significantly the intercalation energy, by up to 0.13 eV at 300 K. Analysis of diffusion barriers shows that Li, Na, and Mg diffusion in β-Sn is anisotropic with migration barriers along the (001) direction (respectively, 0.01, 0.22, and 0.07 eV) significantly lower than those in α-Sn (respectively, 0.20, 0.52, and 0.40 eV)

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

Muon Capture in Nuclei: an Ab Initio Approach Based on Quantum Monte Carlo Methods

An ab initio quantum Monte Carlo method is introduced for calculating total rates of muon weak capture in light nuclei with mass number $A \leq 12$. As a first application of the method, we perform a calculation of the rate in $^4$He in a dynamical framework based on realistic two- and three-nucleon interactions and realistic nuclear charge-changing weak currents. The currents include one- and two-body terms induced by $\pi$- and $\rho$-meson exchange, and $N$-to-$\Delta$ excitation, and are constrained to reproduce the empirical value of the Gamow-Teller matrix element in tritium. We investigate the sensitivity of theoretical predictions to current parametrizations of the nucleon axial and induced pseudoscalar form factors as well as to two-body contributions in the weak currents. The large uncertainties in the measured values obtained from bubble-chamber experiments (carried out over 50 years ago) prevent us from drawing any definite conclusions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Electronic and thermoelectric properties of InN studied using ab initio density functional theory and Boltzmann transport calculations

The thermoelectric properties of indium nitride in the most stable wurtzite phase (w-InN) as a function of electron and hole concentrations and temperature were studied by solving the semiclassical Boltzmann transport equations in conjunction with ab initio electronic structure calculations, within Density Functional Theory. Based on maximally localized Wannier function basis set and the ab initio band energies, results for the Seebeck coefficient are presented and compared with available experimental data for n-type as well as p-type systems. Also, theoretical results for electric conductivity and power factor are presented. Most cases showed good agreement between the calculated properties and experimental data for w-InN unintentionally and p-type doped with magnesium. Our predictions for temperature and concentration dependences of electrical conductivity and power factor revealed a promising use of InN for intermediate and high temperature thermoelectric applications. The rigid band approach and constant scattering time approximation were utilized in the calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A Guided Tour of ab initio Nuclear Many-Body Theory

Over the last decade, new developments in Similarity Renormalization Group techniques and nuclear many-body methods have dramatically increased the capabilities of ab initio nuclear structure and reaction theory. Ground and excited-state properties can be computed up to the tin region, and from the proton to the presumptive neutron drip lines, providing unprecedented opportunities to confront two- plus three-nucleon interactions from chiral Effective Field Theory with experimental data. In this contribution, I will give a broad survey of the current status of nuclear many-body approaches, and I will use selected results to discuss both achievements and open issues that need to be addressed in the coming decade.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ab initio quantum scattering calculations for the CO–O 2 system and a new CO–O 2 potential energy surface: O 2 and air broadening of the R(0) line in CO

We present ab initio calculations of the collisional broadening of the R(0) pure rotational line in CO (at 115 GHz) perturbed by O 2 . Our calculations are done in a fully quantum way by solving close-coupling quantum-scattering equations without any approximations. We also report a new, highly accurate CO–O 2 potential energy surface on which we did the quantum-scattering calculations. The calculated collisional broadening agrees with the available experimental data in a wide temperature range. The calculated collisional shift is negligible compared to the broadening, which is also consistent with the experimental data. We combine this result with our previous calculations for the same line in CO perturbed by N 2 [Józwiak et al., J. Chem. Phys. 154, 054314 (2021)]; the obtained air-perturbed broadening of the R(0) pure rotational line in CO and its temperature dependence perfectly agree with the HITRAN database. This result constitutes an important step toward developing a methodology for providing accurate ab initio reference data on spectroscopic collisional line-shape parameters for molecular systems relevant to the Earth’s atmosphere and for populating spectroscopic line-by-line databases.

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

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Ab initio phonon point defect scattering and thermal transport in graphene

Here, we study the scattering of phonons from point defects and their effect on lattice thermal conductivity κ using a parameter-free ab initio Green's function methodology. Specifically, we focus on the scattering of phonons by boron (B), nitrogen (N), and phosphorus substitutions as well as single- and double-carbon vacancies in graphene. We show that changes of the atomic structure and harmonic interatomic force constants locally near defects govern the strength and frequency trends of the scattering of out-of-plane acoustic (ZA) phonons, the dominant heat carriers in graphene. ZA scattering rates due to N substitutions are nearly an order of magnitude smaller than those for B defects despite having similar mass perturbations. Furthermore, ZA phonon scattering rates from N defects decrease with increasing frequency in the lower-frequency spectrum in stark contrast to expected trends from simple models. ZA phonon-vacancy scattering rates are found to have a significantly softer frequency dependence (~ω 0 ) in graphene than typically employed in phenomenological models. The rigorous Green's function calculations demonstrate that typical mass-defect models do not adequately describe ZA phonon-defect scattering rates. Our ab initio calculations capture well the trend of κ vs vacancy density from experiments, though not the magnitudes. In conclusion, this work elucidates important insights into phonon-defect scattering and thermal transport in graphene, and demonstrates the applicability of first-principles methods toward describing these properties in imperfect materials.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Liquid Be, Ca and Ba. An orbital-free ab-initio molecular dynamics study

Several static and dynamic properties of liquid beryllium (l-Be), liquid calcium (l-Ca) and liquid barium (l-Ba) near their triple point have been evaluated by the orbital-free ab initio molecular dynamics method (OF-AIMD), where the interaction between valence electrons and ions is described by means of local pseudopotentials. These local pseudopotentials used were constructed through a force-matching process with those obtained from a Kohn-Sham ab initio molecular dynamics study (KS-AIMD) of a reduced system with non-local pseudopotentials. The calculated static structures show good agreement with the available experimental data, including an asymmetric second peak in the structure factor which has been linked to the existence of a marked icosahedral short-range order in the liquid. As for the dynamic properties, we obtain collective density excitations whose associated dispersion relations exhibit a positive dispersion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Atomic Diffusivities of Yttrium, Titanium and Oxygen Calculated by Ab Initio Molecular Dynamics in Molten 316L Oxide-Dispersion-Strengthened Steel Fabricated via Additive Manufacturing

Oxide-dispersion-strengthened (ODS) steels have long been viewed as a prime solution for harsh environments. However, conventional manufacturing of ODS steels limits the final product geometry, is difficult to scale up to large components, and is expensive due to multiple highly involved, solid-state processing steps required. Additive manufacturing (AM) can directly incorporate dispersion elements (e.g., Y, Ti and O) during component fabrication, thus bypassing the need for an ODS steel supply chain, the scale-up challenges of powder processing routes, the buoyancy challenges associated with casting ODS steels, and the joining issues for net-shape component fabrication. In the AM process, the diffusion of the dispersion elements in the molten steel plays a key role in the precipitation of the oxide particles, thereby influencing the microstructure, thermal stability and high-temperature mechanical properties of the resulting ODS steels. In this work, the atomic diffusivities of Y, Ti, and O in molten 316L stainless steel (SS) as functions of temperature are determined by ab initio molecular dynamics simulations. The latest Vienna Ab initio Simulation Package (VASP) package that incorporates an on-the-fly machine learning force field for accelerated computation is used. At a constant temperature, the time-dependent coordinates of the target atoms in the molten 316L SS were analyzed in the form of mean square displacement in order to obtain diffusivity. The values of the diffusivity at multiple temperatures are then fitted to the Arrhenius form to determine the activation energy and the pre-exponential factor. Given the challenges in experimental measurement of atomic diffusivity at such high temperatures and correspondingly the lack of experimental data, this study provides important physical parameters for future modeling of the oxide precipitation kinetics during AM process.

Chemistry↗
DOE OSTI2015

Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

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

Efficient ab initio calculations of electron-defect scattering and defect-limited carrier mobility

Electron-defect ($e - \mathrm{d}$) interactions govern charge carrier dynamics at low temperature, where they limit the carrier mobility and give rise to phenomena of broad relevance in condensed matter physics. Ab initio calculations of $e - \mathrm{d}$ interactions are still in their infancy, mainly because they require large supercells and computationally expensive workflows. Here we develop an efficient ab initio approach for computing elastic $e - \mathrm{d}$ interactions, their associated $e - \mathrm{d}$ relaxation times (RTs), and the low-temperature defect-limited carrier mobility. The method is applied to silicon with simple neutral defects, such as vacancies and interstitials. Contrary to conventional wisdom, the computed $e - \mathrm{d}$ RTs depend strongly on carrier energy and defect type, and the defect-limited mobility is temperature dependent. These results highlight the shortcomings of widely employed heuristic models of $e - \mathrm{d}$ interactions in materials. Finally, our method opens avenues for studying $e - \mathrm{d}$ scattering and low-temperature charge transport from first principles.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Specific interactions between DNA and regulatory protein controlled by ligand-binding: Ab initio molecular simulation

The catabolite activator protein (CAP) is one of the regulatory proteins controlling the transcription mechanism of gene. Biochemical experiments elucidated that the complex of CAP with cyclic AMP (cAMP) is indispensable for controlling the mechanism, while previous molecular simulations for the monomer of CAP+cAMP complex revealed the specific interactions between CAP and cAMP. However, the effect of cAMP-binding to CAP on the specific interactions between CAP and DNA is not elucidated at atomic and electronic levels. We here considered the ternary complex of CAP, cAMP and DNA in solvating water molecules and investigated the specific interactions between them at atomic and electronic levels using ab initio molecular simulations based on classical molecular dynamics and ab initio fragment molecular orbital methods. The results highlight the important amino acid residues of CAP for the interactions between CAP and cAMP and between CAP and DNA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Comparing (stochastic-selection) ab initio multiple spawning with trajectory surface hopping for the photodynamics of cyclopropanone, fulvene, and dithiane

Ab Initio Multiple Spawning (AIMS) simulates the excited-state dynamics of molecular systems by representing nuclear wavepackets in a basis of coupled traveling Gaussian functions, called trajectory basis functions (TBFs). New TBFs are spawned when nuclear wavepackets enter regions of strong nonadiabaticity, permitting the description of non-Born–Oppenheimer processes. The spawning algorithm is simultaneously the blessing and the curse of the AIMS method: it allows for an accurate description of the transfer of nuclear amplitude between different electronic states, but it also dramatically increases the computational cost of the AIMS dynamics as all TBFs are coupled. Recently, a strategy coined stochastic-selection AIMS (SSAIMS) was devised to limit the ever-growing number of TBFs and tested on simple molecules. In this work, we use the photodynamics of three different molecules—cyclopropanone, fulvene, and 1,2-dithiane—to investigate (i) the potential of SSAIMS to reproduce reference AIMS results for challenging nonadiabatic dynamics, (ii) the compromise achieved by SSAIMS in obtaining accurate results while using the smallest average number of TBFs as possible, and (iii) the performance of SSAIMS in comparison to the mixed quantum/classical method trajectory surface hopping (TSH)—both in terms of its accuracy and computational cost. We show that SSAIMS can accurately reproduce the AIMS results for the three molecules considered at a much cheaper computational cost, often close to that of TSH. We deduce from these tests that an overlap-based criterion for the stochastic-selection process leads to the best agreement with the reference AIMS dynamics for the smallest average number of TBFs.

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

Ab initio calculations of the carbon and oxygen isotopes: Energies, correlations, and superfluid pairing

We perform ab initio nuclear lattice calculations of the neutron-rich carbon and oxygen isotopes using high-fidelity chiral interactions. We find good agreement with the observed binding energies and compute correlations associated with each two-nucleon interaction channel. For the isospin T = 1 channels, we show that the dependence on T z provides a measure of the correlations among the extra neutrons in the neutron-rich nuclei. For the spin-singlet S-wave channel, we observe that any paired neutron interacts with the nuclear core as well as its neutron pair partner, while any unpaired neutron interacts primarily with only the nuclear core. For the other partial waves, the correlations among the extra neutrons grow more slowly and smoothly with the number of neutrons. These general patterns are observed in both the carbon and oxygen isotopes and may be universal features that appear in many neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

A high level Ab initio study of the anionic hydrogen-bonded complexes FH-CN-, FH-NC-, H2O-CN- and H2O-NC-

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗
NASA NTRS1989

A high-level ab initio study of the anionic hydrogen-bonded complexes FH-CN(-), FH-NC(-), H2O-CN(-), and H2O-NC(-)

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗