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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Multiple time step integrators in ab initio molecular dynamics

Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we present two schemes that enable efficient time-scale separation in ab initio calculations: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.

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

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Unified ab initio description of Fröhlich electron-phonon interactions in two-dimensional and three-dimensional materials

Ab initio calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation, so as to eliminate artificial interactions between periodic images of the 2D layer. While these techniques are well established, the transition of the Fröhlich coupling from three to two dimensions has not been investigated. Furthermore, it remains unclear what error one makes when describing 2D systems using the standard bulk formalism in a periodic supercell geometry. In this work, we generalize previous work on the ab initio Fröhlich electron-phonon matrix element in bulk materials by investigating the electrostatic potential of atomic dipoles in a periodic supercell consisting of a 2D material and a continuum dielectric slab. We obtain a unified expression for the matrix element, which reduces to the existing formulas for three-dimensional and 2D systems when the interlayer separation tends to zero or infinity, respectively. This expression enables an accurate description of the Fröhlich matrix element in 2D systems without resorting to Coulomb truncation. We validate our approach by direct ab initio density-functional perturbation theory calculations for monolayer BN and MoS 2 , and we provide a simple expression for the 2D Fröhlich matrix element that can be used in model Hamiltonian approaches. The formalism outlined in this work may find applications in calculations of polarons, quasiparticle renormalization, transport coefficients, and superconductivity, in 2D and quasi-2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Natural orbitals for the ab initio no-core configuration interaction approach

Ab initio no-core configuration interaction (NCCI) calculations for the nuclear many-body problem have traditionally relied upon an antisymmetrized product (Slater determinant) basis built from harmonic oscillator orbitals. The accuracy of such calculations is limited by the finite dimensions which are computationally feasible for the truncated many-body space. We therefore seek to improve the accuracy obtained for a given basis size by optimizing the choice of single-particle orbitals. Natural orbitals, which diagonalize the one-body density matrix, provide a basis which maximizes the occupation of low-lying orbitals, thus accelerating convergence in a configuration-interaction basis, while also possibly providing physical insight into the single-particle structure of the many-body wave function. We describe the implementation of natural orbitals in the NCCI framework and examine the nature of the natural orbitals thus obtained, the properties of the resulting many-body wave functions, and the convergence of observables. After taking 3 He as an illustrative testbed, we explore aspects of NCCI calculations with natural orbitals for the ground state of the p-shell neutron halo nucleus 6 He .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2002

Quantal Study of the Exchange Reaction for N + N2 using an ab initio Potential Energy Surface

The N + N2 exchange rate is calculated using a time-dependent quantum dynamics method on a newly determined ab initio potential energy surface (PES) for the ground A" state. This ab initio PES shows a double barrier feature in the interaction region with the barrier height at 47.2 kcal/mol, and a shallow well between these two barriers, with the minimum at 43.7 kcal/mol. A quantum dynamics wave packet calculation has been carried out using the fitted PES to compute the cumulative reaction probability for the exchange reaction of N + N2(J=O). The J - K shift method is then employed to obtain the rate constant for this reaction. The calculated rate constant is compared with experimental data and a recent quasi-classical calculation using a LEPS PES. Significant differences are found between the present and quasiclassical results. The present rate calculation is the first accurate 3D quantal dynamics study for N + N2 reaction system and the ab initio PES reported here is the first such surface for N3.

Wang, Dunyou↗
DOE OSTI2021

Importance of bulk excitations and coherent electron-photon-phonon scattering in photoemission from PbTe(111): Ab initio theory with experimental comparisons

This paper presents a fully ab initio many-body photoemission framework that includes coherent three-body electron-photon-phonon scattering to predict the transverse momentum distributions and the mean transverse energies (MTEs) of bulk photoelectrons from single-crystal photocathodes. The need to develop such a theory stems from the lack of studies that provide complete understanding of the underlying fundamental processes governing the transverse momentum distribution of photoelectrons emitted from single crystals. For example, initial predictions based on density-functional theory calculations of effective electron masses suggested that the (111) surface of PbTe would produce very small MTEs (≤ 15 meV), whereas our experiments yielded MTEs ten to twenty times larger than these predictions and also exhibited a lower photoemission threshold than predicted. The ab initio framework presented in this paper correctly reproduces both the magnitude of the MTEs from our measurements in PbTe(111) and the observed photoemission below the predicted threshold. Our results show that both photoexcitations into states that propagate in the bulk of the material and coherent many-body electron-photon-phonon scattering processes, which initial predictions ignored, play surprisingly important roles in photoemission from PbTe(111). Lastly, from the lessons learned, we recommend a procedure for rapid computational screening of potential single-crystal photocathodes for applications in next-generation ultrafast electron diffraction and X-ray free-electron lasers, which will enable significant advances in condensed matter research.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab initio translationally invariant nucleon-nucleus optical potentials

We combine the ab initio symmetry-adapted no-core shell model (SA-NCSM) with the single-particle Green's function approach to construct optical potentials rooted in first principles. Specifically, we show that total cross sections and phase shifts for neutron elastic scattering from a 4 He target with projectile energies between 0.5 and 10 MeV closely reproduce the experiment. In addition, we discuss an important new development that resolves a long-standing issue with spurious center-of-mass motion in the Green's function formalism for many-body approaches. Furthermore, the new development opens a path for first-principle predictions of cross sections for elastic scattering of single-nucleon projectiles, nucleon capture, and deuteron breakup reactions, feasible for a broad range of open-shell spherical and deformed nuclei in the SA-NCSM approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Deuteron-induced nucleon transfer reactions within an ab initio framework: First application to p -shell nuclei

Low-energy transfer reactions in which a proton is stripped from a deuteron projectile and dropped into a target play a crucial role in the formation of nuclei in both primordial and stellar nucleosynthesis, as well as in the study of exotic nuclei using radioactive beam facilities and inverse kinematics. Here, ab initio approaches have been successfully applied to describe the 3 H(d,n) 4 He and 3 He(d,p) 4 He fusion processes. An ab initio treatment of transfer reactions would also be desirable for heavier targets. In this work, we extend the ab initio description of (d,p) reactions to processes with light p-shell nuclei. As a first application, we study the elastic scattering of deuterium on 7 Li and the 7 Li(d,p) 8 Li transfer reaction based on a two-body Hamiltonian. We use the no-core shell model to compute the wave functions of the nuclei involved in the reaction, and describe the dynamics between targets and projectiles with the help of microscopic-cluster states in the spirit of the resonating group method. The shapes of the excitation functions for deuterons impinging on 7 Li are qualitatively reproduced up to the deuteron breakup energy. The interplay between d– 7 Li and p– 8 Li particle-decay channels determines some features of the 9 Be spectrum above the d+ 7 Li threshold. Our prediction for the parity of the 17.298 MeV resonance is at odds with the experimental assignment. Deuteron stripping reactions with p-shell targets can now be computed ab initio , but calculations are very demanding. Finally, a quantitative description of the 7 Li(d,p) 8 Li reaction will require further work to include the effect of three-nucleon forces and additional decay channels and to improve the convergence rate of our calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Ab Initio Simulations of Martensitic Phase Transformations in NiTi-based High Temperature Ternary Shape Memory Alloys: NiTiHf and NiTiZr

Ab initio simulations of phase stability and martensitic phase transitions are performed for NiTi-based ternary shape memory alloys (SMAs). Specifically, we considered NiTiHf and NiTiZr, which are highly studied for high temperature SMA applications. Previously, we performed investigations of ordered NiTi and related binaries [1,2]. However, similar approaches for chemically disordered compounds present additional difficulties. In this work, special quasi-random structures (SQS) were generated for various compositions, x∈[0,0.5], of Ni0.5Ti(0.5-x)Hfx and Ni0.5Ti(0.5-x)Zrx to capture chemical disorder of off-stoichiometric compounds. Phase stability was evaluated through analysis of finite temperature phonon spectra using temperature dependent effective potential (TDEP) method. Free energies for the cubic B2 phase of NiTiHf and NiTiZr were computed using ab initio thermodynamic integration (AITI) developed previously [1,2]. Free energies for monoclinic B19’ and orthorhombic B33 phases were evaluated via quasi harmonic approximations (QHA). Our results show a critical composition (xc) where the three phases of B2, B19’ and B33 meet, i.e. there is a tri-critical point. For x xc, the transition is between B33 and B2, i.e. it is not a shape memory transition. The approach presented here opens the door to ab initio based predictions of MTT for arbitrary ternary SMAs.

Hessam Malmir↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Ab initio calculations of the optical properties of crystalline and liquid InSb

Ab initio calculations of the electronic and optical properties of InSb were performed for both the crystalline and liquid states. Two sets of atomic structure models for liquid InSb at 900 K were obtained by ab initio molecular dynamics simulations. To reduce the effect of structural peculiarities in the liquid models, an averaging of the two sets of the calculated electronic and optical properties corresponding to the two liquid models was performed. The calculated results indicate that, owing to the phase transition from crystal to liquid, the density of states around the Fermi level increases. As a result, the energy band gap opening near the Fermi level disappears. Consequently, the optical properties change from semiconductor to metallic behavior. Namely, owing to the melting of InSb, the interband transition peaks disappear and a Drude-like dispersion is observed in the optical dielectric functions. The optical absorption at a photon energy of 3.06 eV, which is used in Blu-ray Disc systems, increases owing to the melting of InSb. This increase in optical absorption is proposed to result from the increased optical transitions below 2 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Robust ab initio predictions for dimensionless ratios of 𝐸⁢2 and radius observables. II. Estimation of 𝐸⁢2 transition strengths by calibration to the charge radius

Converged results for 𝐸⁢2 observables are notoriously challenging to obtain in ab initio no-core configuration interaction approaches. Matrix elements of the 𝐸⁢2 operator are sensitive to the large-distance tails of the nuclear wave function, which converge slowly in an oscillator basis expansion. Similar convergence challenges beset ab initio prediction of the nuclear charge radius. However, we exploit systematic correlations between the calculated 𝐸⁢2 and radius observables to yield meaningful predictions for relations among these observables. In particular, we examine ab initio predictions for dimensionless ratios of the form 𝐵⁡(𝐸⁢2)/(𝑒 2 ⁢𝑟 4 ) for nuclei throughout the 𝑝 shell. Finally, meaningful predictions for 𝐸⁢2 transition strengths may then be made by calibrating to the ground-state charge radius if experimentally known.

ab initio calculations↗
DOE OSTI2022

Modeling the Solvation and Acidity of Carboxylic Acids Using an Ab Initio Deep Neural Network Potential

Formic and acetic acid constitute the simplest of carboxylic acids, yet they exhibit fascinating chemistry in the condensed phase such as proton transfer and dimerization. The go-to method of choice for modeling these rare events have been accurate but expensive ab-initio molecular dynamics simulations. Here, we present a deep neural network potential trained using accurate ab-initio data that can be used in tandem with enhanced-sampling methods to perform an efficient exploration of the free-energy surface of aqueous solutions of weak carboxylic acids. In particular, we show that our model captures proton dissociation and provides a good estimate of the pK a , as well as the dimerization of formic and acetic acid. This provides a suitable starting point for applications in different research areas where computational efficiency coupled with the accuracy of ab-initio methods is required.

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

Ab initio Bogoliubov coupled cluster theory for open-shell nuclei

Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass A≈130 on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking. Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory. Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in m scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation. Results: Proof-of-principle calculations in an N max =6 spherical harmonic oscillator basis for 16,18 O and 18 Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while 20 O and 20 Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations. Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Nature of Hyperfine Interactions in TbPc 2 Single-Molecule Magnets: Multiconfigurational Ab Initio Study

Lanthanide-based single-ion magnetic molecules can have large magnetic hyperfine interactions as well as large magnetic anisotropy. Recent experimental studies reported tunability of these properties by changes of chemical environments or by application of external stimuli for device applications. In order to provide insight onto the origin and mechanism of such tunability, here we investigate the magnetic hyperfine and nuclear quadrupole interactions for a 159 Tb nucleus in TbPc 2 (Pc = phthalocyanine) single-molecule magnets using multiconfigurational ab initio methods including spin–orbit interaction. Since the electronic ground and first-excited (quasi)doublets are well separated in energy, the microscopic Hamiltonian can be mapped onto an effective Hamiltonian with an electronic pseudospin S = 1/2. From the ab initio calculated parameters, we find that the magnetic hyperfine coupling is dominated by the interaction of the Tb nuclear spin with electronic orbital angular momentum. The asymmetric 4f-like electronic charge distribution leads to a strong nuclear quadrupole interaction with significant transverse terms for the molecule with low symmetry. The ab initio calculated electronic–nuclear spectrum including the magnetic hyperfine and quadrupole interactions is in excellent agreement with the experiment. We further find that the transverse quadrupole interactions significantly influence the avoided level crossings in magnetization dynamics and that the molecular distortions affect mostly the Fermi contact terms as well as the transverse quadrupole interactions. Synopsis Hyperfine interactions for 159 Tb nucleus in TbPc 2 single-molecule magnets are investigated from first-principles using multiconfigurational calculations. Strong nuclear quadrupole and magnetic hyperfine coupling are found, with the latter being dominated by the paramagnetic spin-orbital mechanism. Furthermore, we construct an ab initio pseudospin Hamiltonian and obtain the electronic-nuclear spectrum that is in excellent agreement with experiment. The Zeeman diagram is calculated, and the magnetization dynamics is discussed. The effects of molecular distortions are researched.

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

Efficacy of the SU(3) scheme for ab initio large-scale calculations beyond the lightest nuclei

We report on the computational characteristics of ab initio nuclear structure calculations in a symmetry-adapted no-core shell model (SA-NCSM) framework. We examine the computational complexity of the current implementation of the SA-NCSM approach, dubbed LSU3shell, by analyzing ab initio results for 6 Li and 12 C in large harmonic oscillator model spaces and SU3-selected subspaces. We demonstrate LSU3shell’s strong-scaling properties achieved with highly-parallel methods for computing the many-body matrix elements. Results compare favorably with complete model space calculations and significant memory savings are achieved in physically important applications. In particular, a well-chosen symmetry-adapted basis affords memory savings in calculations of states with a fixed total angular momentum in large model spaces while exactly preserving translational invariance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Neural network representation of electronic structure from ab initio molecular dynamics

Despite their rich information content, electronic structure data amassed at high volumes in ab initio molecular dynamics simulations are generally under-utilized. Here we introduce a transferable high-fidelity neural network representation of such data in the form of tight-binding Hamiltonians for crystalline materials. This predictive representation of ab initio electronic structure, combined with machine-learning boosted molecular dynamics, enables efficient and accurate electronic evolution and sampling. When it is applied to a one-dimension charge-density wave material, carbyne, we are able to compute the spectral function and optical conductivity in the canonical ensemble. The spectral functions evaluated during soliton-antisoliton pair annihilation process reveal significant renormalization of low-energy edge modes due to retarded electron-lattice coupling beyond the Born–Oppenheimer limit. The availability of an efficient and reusable surrogate model for the electronic structure dynamical system will enable calculating many interesting physical properties, paving the way to previously inaccessible or challenging avenues in materials modeling.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1992

A global ab initio potential for HCN/HNC, exact vibrational energies, and comparison to experiment

An ab initio (i.e., from first principles) calculation of vibrational energies of HCN and HNC is reported. The vibrational calculations were done with a new potential derived from a fit to 1124 ab initio electronic energies which were calculated using the highly accurate CCSD(T) coupled-cluster method in conjunction with a large atomic natural orbital basis set. The properties of this potential are presented, and the vibrational calculations are compared to experiment for 54 vibrational transitions, 39 of which are for zero total angular momentum, J = 0, and 15 of which are for J = 1. The level of agreement with experiment is unprecedented for a triatomic with two nonhydrogen atoms, and demonstrates the capability of the latest computational methods to give reliable predictions on a strongly bound triatomic molecule at very high levels of vibrational excitation.

Bentley, Joseph A.↗