Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

Ab Initio Prediction of Fluorescence Lifetimes Involving Solvent Environments by Means of COSMO and Vibrational Broadening

The fluorescence lifetime is a key property of fluorophores that can be utilized for microenvironment probing, analyte sensing, and multiplexing as well as barcoding applications. For the rational design of lifetime probes and barcodes, theoretical methods have been developed to enable the ab initio prediction of this parameter, which depends strongly on interactions with solvent molecules and other chemical species in the emitters immediate environment. In this work, we investigate how a conductor-like screening model (COSMO) can account for variations in fluorescence lifetimes that are caused by such fluorophore–solvent interactions. Therefore, we calculate vibrationally broadened fluorescence spectra using the nuclear ensemble method to obtain distorted molecular geometries to sample the electronic transitions with time-dependent density functional theory (TDDFT). The influence of the solvent on fluorescence lifetimes is accounted for with COSMO. For example, for 4-hydroxythiazole fluorophore containing different heteroatoms and acidic and basic moieties in aprotic and protic solvents of varying polarity, this approach was compared to experimentally determined lifetimes in the same solvents. Furthermore, our results demonstrate a good correlation between theoretically predicted and experimentally measured fluorescence lifetimes except for the polar solvents ethanol and acetonitrile that can specifically interact with the heteroatoms and the carboxylic acid of the thiazole derivative.

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

Analytic model of the γ -surface deviation and influence on the stacking fault width between partial dislocations

The stacking fault width (R$_{eq}$) between partial dislocations within an FCC crystalline lattice characterizes the onset of numerous plastic flow mechanisms, as well as the relationship between material strength and grain size. Continuum models traditionally consider a complete unit of dislocation slip ( b ) along the $\langle{110}$$\rangle$ direction distributed between two discrete partial dislocations, each with a fixed partial Burgers vector ( b p ), which bound a stacking fault. Across the stacking fault, the vectorial slip is assumed to be constant, yielding a constant intrinsic stacking fault energy density, $\gamma$$_{isf}$ Here, we demonstrate that the vectorial displacement path taken in accomplishing a complete unit of slip (b) deviates from the expected displacement path containing the local minima, $\gamma$$_{isf}$ leading to a correction in the nominally derived stacking fault width. The magnitude of the correction depends on both the net orientation of the dislocation within the lattice, and also the degree of relaxation of each partial Burgers vector along the $\langle{\bar{1}12}$$\rangle$ direction. We derive a simple analytic model for the corrected stacking fault width, explicitly accounting for the deviation, by introducing a variable ξ (0 < ξ < 1), which governs the magnitude of each partial dislocation component along the $\langle{\bar{1}12}$$\rangle$ direction. Significantly, our model predicts a correction to R$_{eq}$ of ~$\mathscr{O}$ (|b|) for the nominally screw dislocation, and has no noticeable influence on nominally edge dislocations. We apply our model towards computing the stacking fault width of several FCC metals, and demonstrate excellent agreement both with our own numerical data as well as that obtained from ab initio and Molecular Statics (MS) methods within the literature. Lastly, the results from this study demonstrate that, upon judicious application, discrete linear elastic models are successful in reproducing elastic interactions as computed from higher fidelity models on the spatial scale of metallic dislocation cores.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

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

Comprehensive End-To-End Design of Novel High Energy Density Materials: III. Fused Heterocyclic Energetic Compounds

In contrast to synthesis and experimental characterization that are usually expensive and time-consuming, a coupled combination of ab initio and group additive theoretical methods represents a powerful alternative to computationally design and characterize new materials with targeted properties. Heterocyclic energetic materials are considered attractive candidates to replace conventional high explosives, such as PETN, RDX, HMX, and TNT. Heterocycles enable an intriguing opportunity for combinatorial design of new classes of energetic materials with targeted performance and sensitivity properties and offer ways for an accurate analysis of structure–property correlations of the materials. By using a computational strategy that links methods of quantum chemistry and semiempirical statistical analysis, we designed, modeled, and characterized a series of novel linear and fused explosive compounds. Proposed materials with best combinations of performance and sensitivity parameters can be now considered as potential candidates for synthesis and experimental characterization.

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

Inelastic Triatom-Atom Quantum Close-Coupling Dynamics in Full Dimensionality: All Rovibrational Mode Quenching of Water Due to the H Impact on a Six-Dimensional Potential Energy Surface

The rovibrational level populations, and subsequent emission in various astrophysical environments, are driven by inelastic collision processes. The available rovibrational rate coefficients for water have been calculated using a number of approximations. We present a numerically exact calculation for the rovibrational quenching for all water vibrational modes due to collisions with atomic hydrogen. The scattering theory implements a quantum close-coupling (CC) method on a high level ab initio six-dimensional (6D) potential energy surface (PES). Total rovibrational quenching cross sections for excited bending levels were compared with earlier results on a 4D PES with the rigid-bender closecoupling (RBCC) approximation. General agreement between 6D-CC and 4DRBCC calculations are found, but differences are evident including the energy and amplitude of low-energy orbiting resonances. Quenching cross sections from the symmetric and asymmetric stretch modes are provided for the first time. The current 6D-CC calculation provides accurate inelastic data needed for astrophysical modeling.

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

Design Principles for Trap-Free CsPbX 3 Nanocrystals: Enumerating and Eliminating Surface Halide Vacancies with Softer Lewis Bases

We introduce a general surface passivation mechanism for cesium lead halide perovskite materials (CsPbX 3 , X = Cl, Br, I) that is supported by a combined experimental and theoretical study of the nanocrystal surface chemistry. A variety of spectroscopic methods are employed together with ab initio calculations to identify surface halide vacancies as the predominant source of charge trapping. The number of surface traps per nanocrystal is quantified by 1 H NMR spectroscopy, and that number is consistent with a simple trapping model in which surface halide vacancies create deleterious under-coordinated lead atoms. These halide vacancies exhibit trapping behavior that differs among CsPbCl 3 , CsPbBr 3 , and CsPbI 3 . Ab initio calculations suggest that introduction of anionic X-type ligands can produce trap-free band gaps by altering the energetics of lead-based defect levels. General rules for selecting effective passivating ligand pairs are introduced by considering established principles of coordination chemistry. Introducing softer, anionic, X-type Lewis bases that target under-coordinated lead atoms results in absolute quantum yields approaching unity and monoexponential luminescence decay kinetics, thereby indicating full trap passivation. This work provides a systematic framework for preparing highly luminescent CsPbX 3 nanocrystals with variable compositions and dimensionalities, thereby improving the fundamental understanding of these materials and informing future synthetic and post-synthetic efforts toward trap-free CsPbX 3 nanocrystals.

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

Theoretical study of the CO 2 –O 2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO 2 –O 2 ) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l –3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO 2 –O 2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ~ –243 cm –1 . Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ~ –218 cm –1 . The transition structure connecting the two minima is found at V ~ –211 cm –1 . We also performed comparisons with some CO 2 –X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO 2 and O 2 . Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO 2 –O 2 mixtures.

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

Electron-impact vibrational excitation of cyclopropane

We report a very detailed test of the ab initio discrete momentum representation (DMR) method of calculating vibrational excitation of polyatomic molecules by electron impact, by comparison of its results with an extensive set of experimental data, covering the entire range of scattering angles from 10{sup ∘} to 180{sup ∘} and electron energies from 0.4 to 20 eV. The DMR calculations were carried out by solving the two-channel Lippmann-Schwinger equation in the momentum space, and the interaction between the scattered electron and the target molecule was described by exact static-exchange potential corrected by a density functional theory (DFT) correlation-polarization interaction that models target’s response to the field of incoming electron. The theory is found to quantitatively reproduce the measured spectra for all normal modes, even at the difficult conditions of extreme angles and at low energies, and thus provides full understanding of the excitation mechanism. It is shown that the overlap of individual vibrational bands caused by limited experimental resolution and rotational excitation must be properly taken into account for correct comparison of experiment and theory. By doing so, an apparent discrepancy between published experimental data could be reconciled. A substantial cross section is found for excitation of the non-symmetric HCH twisting mode ν{sub 4} of A{sub 1}{sup ″} symmetry by the 5.5 eV A{sub 2}{sup ′} resonance, surprisingly because the currently accepted selection rules predict this process to be forbidden. The DMR theory shows that the excitation is caused by an incoming electron in an f-wave of A{sub 2}{sup ′} symmetry which causes excitation of the non-symmetric HCH twisting mode ν{sub 4} of the A{sub 1}{sup ″} symmetry and departs in p- and f-waves of A{sub 2}{sup ″} symmetry.

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

Thermal relaxation of molecular oxygen in collisions with nitrogen atoms

Investigation of O{sub 2}–N collisions is performed by means of the quasi-classical trajectory method on the two lowest ab initio potential energy surfaces at temperatures relevant to hypersonic flows. A complete set of bound–bound and bound–free transition rates is obtained for each precollisional rovibrational state. Special attention is paid to the vibrational and rotational relaxations of oxygen as a result of chemically non-reactive interaction with nitrogen atoms. The vibrational relaxation of oxygen partially occurs via the formation of an intermediate NO{sub 2} complex. The efficient energy randomization results in rapid vibrational relaxation at low temperatures, compared to other molecular systems with a purely repulsive potential. The vibrational relaxation time, computed by means of master equation studies, is nearly an order of magnitude lower than the relaxation time in N{sub 2}–O collisions. The rotational nonequilibrium starts to play a significant effect at translational temperatures above 8000 K. The present work provides convenient relations for the vibrational and rotational relaxation times as well as for the quasi-steady dissociation rate coefficient and thus fills a gap in data due to a lack of experimental measurements for this system.

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

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Hidden Mn magnetic-moment disorder and its influence on the physical properties of medium-entropy NiCoMn solid solution alloys

In this work, the ab initio Korringa-Kohn-Rostoker method combined with the coherent potential approximation (CPA) was employed to investigate the electronic, magnetic, and transport properties of medium-entropy face-centered-cubic (fcc) NiCoMn solid solution alloys. By comparing the CPA electronic structure with that from supercell calculations, we uncovered an unconventional CPA ground state, which correctly distinguishes two equally populated Mn CPA components—with large spin moments but opposite orientations. Using the spin spiral calculations, we further demonstrated that this ground state is most energetically favorable in the presence of spin noncollinearity, and no significant longitudinal spin fluctuation is observed, justifying the applicability of the Heisenberg model. The finite-temperature magnetism was further studied using different approximations based on the Heisenberg model, and we found the Mn moments to be fully disordered at low temperature due to a small net effective Weiss field on Mn. In addition, the magnetic effect on the electron scattering at finite temperatures was evaluated and compared with other scattering mechanisms. Since the magnetization-induced electron scattering is almost saturated in the ground state, (full) spin disorder only yields a small addition to the resistivity, whereas the thermal displacements increase it modestly. Finally, we elucidate the role of hydrostatic pressure on the magnetic and transport properties. In conclusion, these findings reflect the importance of the magnetic signatures on the physical properties of alloys, and they provide a window into magnetism-controlled electronic structure and energy dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Real-time GW -BSE investigations on spin-valley exciton dynamics in monolayer transition metal dichalcogenide

We develop an ab initio nonadiabatic molecular dynamics (NAMD) method based on GW plus real-time Bethe-Salpeter equation (GW + rtBSE-NAMD) for the spin-resolved exciton dynamics. From investigations on MoS 2 , we provide a comprehensive picture of spin-valley exciton dynamics where the electron-phonon (e-ph) scattering, spin-orbit interaction (SOI), and electron-hole (e-h) interactions come into play collectively. In particular, we provide a direct evidence that e-h exchange interaction plays a dominant role in the fast valley depolarization within a few picoseconds, which is in excellent agreement with experiments. Moreover, there are bright-to-dark exciton transitions induced by e-ph scattering and SOI. Our study proves that e-h many-body effects are essential to understand the spin-valley exciton dynamics in transition metal dichalcogenides and the newly developed GW + rtBSE-NAMD method provides a powerful tool for exciton dynamics in extended systems with time, space, momentum, energy, and spin resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Ab initio thermodynamic characteristics of the formation of oxygen vacancies, and boron, carbon, and nitrogen impurity centers in anatase

Using the ab initio projector augmented wave (PAW) method, calculations are performed for the electronic energy-band structure of titanium dioxide having the structure of anatase doped with boron, nitrogen, and carbon. Thermodynamic characteristics are determined for the formation of impurity centers, such as the preference energy for the interstitial position, the energy of impurity oxidation, and the energy of oxygen vacancy formation. It is shown that under the conditions of thermodynamic equilibrium the interstitial position of boron atoms is stable, whereas carbon atoms, depending on the oxygen pressure, can occupy both interstitial positions and substitutional positions of oxygen atoms, and nitrogen atoms replace oxygen atoms. It is shown that the presence of oxygen vacancies promotes the thermodynamic stability of carbon and nitrogen atoms. The obtained densities of electronic states correspond to ESR spectroscopy data, which indicates the presence of spin-polarized electrons in the states of the oxygen vacancy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Symmetry reduction of tensor networks in many-body theory: I. Automated symbolic evaluation of SU(2) algebra

Abstract The ongoing progress in (nuclear) many-body theory is accompanied by an ever-rising increase in complexity of the underlying formalisms used to solve the stationary Schrödinger equation. The associated working equations at play in state-of-the-art ab initio nuclear many-body methods can be analytically reduced with respect to angular-momentum, i.e. SU (2), quantum numbers whenever they are effectively employed in a symmetry-restricted context. The corresponding procedure constitutes a tedious and error-prone but yet an integral part of the implementation of those many-body frameworks. Indeed, this symmetry reduction is a key step to advance modern simulations to higher accuracy since the use of symmetry-adapted tensors can decrease the computational complexity by orders of magnitude. While attempts have been made in the past to automate the (anti-) commutation rules linked to Fermionic and Bosonic algebras at play in the derivation of the working equations, there is no systematic account to achieve the same goal for their symmetry reduction. In this work, the first version of an automated tool performing graph-theory-based angular-momentum reduction is presented. Taking the symmetry-unrestricted expressions of a generic tensor network as an input, the code provides their angular-momentum-reduced form in an error-safe way in a matter of seconds. Several state-of-the-art many-body methods serve as examples to demonstrate the generality of the approach and to highlight the potential impact on the many-body community.

Tichai, A.↗