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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within {approx}30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H{sub 2}O{sup +}+ H{sub 2}O {yields} OH + H{sub 3}O{sup +}. The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM/MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with {approx}40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H{sub 2}O{sup +}{sub (aq)} species.

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

Evolution of local atomic structure during solidification of Al2Au liquid: An ab initio study

The local atomic structure evolution in Al2Au alloy during solidification from 2000 K to 400 K was studied by ab initio molecular dynamics simulations and analyzed using the structure factor, pair correlation functions, bond angle distributions, the Honeycutt-Anderson (HA) index and Voronoi tessellation methods. It was found that the icosahedral-like clusters are negligible in the Al2Au stable liquid and supercooled liquid states, and the most abundant clusters are those having HA indices of 131 and 120 or Voronoi indices of < 0,4,4,0 >, < 0,3, 6,0 > and < 0,4,4,2 > with coordination numbers of 8, 9 and 10, respectively. These clusters are similar to the local atomic structures in the CaF2-type Al2Au crystal, revealing the existence of structure heredity between liquid and crystalline phase in Al2Au alloy. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Ab Initio Study of Carbon Impurities in Cu$_2$ ZnSnS$_4$

Cu 2 ZnSnS 4 (CZTS) is a highly promising absorber material for efficient and low-cost solar cells. In a cost-effective method of production of high-quality CZTS, nanocrystals of precursor compounds, solved in an amine solution, are deposited on a substrate and sulfurized. Here, this can generally lead to considerable residual carbon in the final samples, in the form of separate amorphous layers as observed experimentally, but also incorporated in the CZTS matrix. In this study, preferred configurations, bonding characteristics, formation energies, and ionization levels of carbon-related defects in CZTS are calculated, including interstitials, single- and double-carbon antisites, and complexes of those with intrinsic defects. For different growth conditions of carbon, the set of dominant defects is determined, and the resulting total carbon concentrations are discussed with respect to experimental values.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2014

Ab Initio Thermodynamic Model for Magnesium Carbonates and Hydrates

An ab initio thermodynamic framework for predicting properties of hydrated magnesium carbonate minerals has been developed using density-functional theory linked to macroscopic thermodynamics through the experimental chemical potentials for MgO, water, and CO2. Including semiempirical dispersion via the Grimme method and small corrections to the generalized gradient approximation of Perdew, Burke, and Ernzerhof for the heat of formation yields a model with quantitative agreement for the benchmark minerals brucite, magnesite, nesquehonite, and hydromagnesite. The model shows how small differences in experimental conditions determine whether nesquehonite, hydromagnesite, or magnesite is the result of laboratory synthesis from carbonation of brucite, and what transformations are expected to occur on geological time scales. Because of the reliance on parameter-free first principles methods, the model is reliably extensible to experimental conditions not readily accessible to experiment and to any mineral composition for which the structure is known or can be hypothesized, including structures containing defects, substitutions, or transitional structures during solid state transformations induced by temperature changes or processes such as water, CO2, or O2 diffusion. Demonstrated applications of the ab initio thermodynamic framework include an independent means to evaluate differences in thermodynamic data for lansfordite, predicting the properties of Mg analogs of Ca-based hydrated carbonates monohydrocalcite and ikaite which have not been observed in nature, and an estimation of the thermodynamics of barringtonite from the stoichiometry and a single experimental observation.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2019

Adsorption Behavior of Organic Molecules: A Study of Benzotriazole on Cu(111) with Spectroscopic and Theoretical Methods

The adsorption of organic molecules on solid substrates is important to applications in fields such as catalysis, photovoltaics, corrosion inhibition, adhesion, and sensors. The molecular level description of the surface–molecule interaction and of the adsorption structures in these complex systems is crucial to understand their properties and function. Here, we present an investigation of one such system, benzotriazole (BTAH) on single-crystal Cu(111) in vacuum conditions. BTAH is the most widely used corrosion inhibitor for copper and thus a molecule of great industrial relevance. We show that the co-application of a wide range of spectroscopic techniques with theoretical methods provides unique insight in the description of the atomistic details of the adsorbed structures. Specifically, spectroscopic photoemission, absorption, and standing wave experiments combined with ab initio computational modeling allowed us to identify that benzotriazole forms overlayers of intact BTAH when deposited at low temperature, and it dissociates into BTA and H at room temperature and above. The dissociated molecule then forms complex structures of mixed chains and dimers of BTA bound to copper adatoms. Our work also reveals that copper adatoms at low concentrations, such as the theoretically predicted superstructures, cannot be resolved by means of current X-ray photoelectron spectroscopy as the modeled Cu 2p spectra are practically indistinguishable from those for a Cu surface without adatoms. Overall this study significantly deepens understanding of BTAH on Cu, a system studied for more than 50 years, and it highlights the benefits of combining spectroscopic and computational methods to obtain a complete picture of a complex adsorption system.

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

Ab Initio Investigation of Structural, Electronic, and Magnetic Properties of Cr-Doped ZnS and ZnSe in Wurtzite Structure

Using the first-principles calculations with the full potential linear muffin–tin orbital (FP-LMTO) method, we investigated the structural, electronic, and magnetic properties of Zn{sub 1−x}Cr{sub x}(S,Se)-diluted magnetic semiconductors (DMSs) in wurtzite structures with varying concentrations (x = 0.0625, x = 0.125, and x = 0.25) of Cr. The electronic properties indicated that Zn{sub 1−x}Cr{sub x}(S,Se), in all concentrations, exhibited half-metallic ferromagnetic (HMF) behavior with spin polarization of 100%. The density of states showed a hybridization between the p(S, Se) and 3d(Cr) states, which created an antibonding state in the gap that stabilized the ferromagnetic ground state linked to the double-exchange mechanism. Therefore, these compounds are highly likely candidates for spintronic applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

DeePKS-kit: A package for developing machine learning-based chemically accurate energy and density functional models

We introduce DeePKS-kit, an open-source software package for developing machine learning based energy and density functional models. DeePKS-kit is interfaced with PyTorch, an open-source machine learning library, and PySCF, an ab initio computational chemistry program that provides simple and customized tools for developing quantum chemistry codes. It supports the DeePHF and DeePKS methods. In addition to explaining the details in the methodology and the software, we also provide an example of developing a chemically accurate model for water clusters.

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

Ab initio study of the linear and nonlinear optical properties of chalcopyrite CdGeAs{sub 2}

We present an ab initio theoretical study of the electronic, linear and nonlinear optical properties of CdGeAs{sub 2} using a pseudopotential plane-wave method. Specifically, we evaluate the band structure, density of states, charge density, the dielectric function {epsilon}({omega}) and the second harmonic generation response susceptibility {chi}{sub 312}{sup (2)}(-2{omega};{omega},{omega}) over a large frequency range. As LDA underestimates the band gap, we have applied the GW approximation method to calculate the quasiparticle band structure and obtain an energy band gap in agreement with experiment. In this case the opening of the gap due to the GW correction can be used as scissor shift to calculate the linear and nonlinear optical properties. The intra- and inter-band contributions to the imaginary part of {chi}{sub 312}{sup (2)}(-2{omega};{omega},{omega}) are presented over a broad energy range. It is found that the small energy gap semiconductor CdGeAs{sub 2} has larger values of {epsilon}{sub 1}(0) and {chi}{sub 312}{sup (2)}(0) than other chalcopyrite structures. - Graphical abstract: Calculated total imaginary part of {chi}{sub 312}{sup (2)}(-2{omega};{omega},{omega}) spectra along with the intra-(2{omega})/(1{omega}) and inter-(2{omega})/(1{omega}) band contributions. Highlights: Black-Right-Pointing-Pointer We find that the charge density around the cation Cd of CdGeAs{sub 2} is spherical. Black-Right-Pointing-Pointer We do not have to rely on experimental gap to determine the self-energy. Black-Right-Pointing-Pointer It is the first time we have calculated the SHG susceptibility using pseudopotential plane-wave. Black-Right-Pointing-Pointer The SHG susceptibility beyond zero frequency was calculated.

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

Theoretical Studies on the Photoelectron and Absorption Spectra of MnO 4 – and TcO 4 –

The tetraoxo pertechnetate anion (TcO4-) is of great interest for nuclear waste management and radiopharmceuticals. To elucidate its electronic structure and to compare with that of its lighter congener MnO4-, the photoelectron and electronic absorption spectra of MnO4 - and TcO4 - are investigated with density functional theory (DFT) and ab initio wave function theory (WFT). The vertical electron detachment energies (VDEs) of MnO4 - obtained with the CR-EOM-CCSD(T) method are in good agreement with the lowest two experimental VDEs; the differences are less than 0.1 eV, representing a significant improvement over the IP-EOMCCSD(T) result in the literature. Combining our CCSD(T) and CR-EOMCCSD( T) results, the first five VDEs of TcO4 - are estimated between 5 and 10 eV with an estimated accuracy of about ±0.2 eV. The vertical excitation energies are determined by using TD-DFT, CR-EOM-CCSD(T), and RASPT2 methods. The excitation energies and the assignments of the spectra are analyzed and partly improved. They are compared with reported SAC-CI results and available experimental data. Both dynamic and nondynamic electron correlations are important in the ground and excited states of MnO4 - and TcO4 -. Nondynamical correlations are particularly relevant in TcO4 - for reliable prediction of excitation energies. In TcO4 - one Rydberg state interlaces but does not mix with the valence excited states, and it disappears in the condensed phase.

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

Influence of c-axis orientation and scandium concentration on infrared active modes of magnetron sputtered Sc{sub x}Al{sub 1−x}N thin films

Doping of wurtzite aluminium nitride (AlN) with scandium (Sc) significantly enhances the piezoelectric properties of AlN. Sc{sub x}Al{sub 1−x}N thin films with different Sc concentrations (x = 0 to 0.15) were deposited by DC reactive magnetron sputtering. Infrared (IR) absorbance spectroscopy was applied to investigate the Sc concentration dependent shift of the IR active modes E{sub 1}(TO) and A{sub 1}(TO). These results are compared to ab initio simulations, being in excellent agreement with the experimental findings. In addition, IR spectroscopy is established as an economical and fast method to distinguish between thin films with a high degree of c-axis orientation and those exhibiting mixed orientations.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Erratum: “A finite difference Davidson procedure to sidestep full ab initio hessian calculation: Application to characterization of stationary points and transition state searches”

Certain key contributions to the field were not presented in the Introduction section of our recent paper. In addition to the variational method described in Ref. 14 and iterative diagonalization methods in Refs. 15–17, we wish to direct the attention of the reader to other significant efforts that have utilized similar principles for various applications. Reiher and Neugebauer developed a mode-tracking algorithm using the Davidson method for calculating only the modes relevant for vibrational analysis. Similarly, a block-Davidson approach via finite differentiation of the gradient was implemented for calculating hessian-vector products for normal mode analysis. Deglmann and Furche developed a method for stationary point characterization with density functional theory where full solution of all coupled-perturbed Kohn-Sham equations, the bottleneck in hessian calculation, could be avoided through iterative calculation of hessian-vector products to determine only the lowest eigenvalues.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗
DOE OSTI2018

Deep Potential Molecular Dynamics: A Scalable Model with the Accuracy of Quantum Mechanics

We introduce a scheme for molecular simulations, the deep potential molecular dynamics (DPMD) method, based on a many-body potential and interatomic forces generated by a carefully crafted deep neural network trained with ab initio data. The neural network model preserves all the natural symmetries in the problem. It is first-principles based in the sense that there are no ad hoc components aside from the network model. We show that the proposed scheme provides an efficient and accurate protocol in a variety of systems, including bulk materials and molecules. In all these cases, DPMD gives results that are essentially indistinguishable from the original data, at a cost that scales linearly with system size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Probing New Bosons and Nuclear Structure with Ytterbium Isotope Shifts

In this Letter, we present mass-ratio measurements on highly charged Yb 42+ ions with a precision of 4 ×10 -12 and isotope-shift measurements on Yb + on the 2 S 1/2 → 2 D 5/2 and 2 S 1/2 → 2 F 7/2 transitions with a precision of 4 ×10 -9 for the isotopes 168,170,172,174,176 Yb. We present a new method that allows us to extract higher-order changes in the nuclear charge distribution along the Yb isotope chain, benchmarking ab initio nuclear structure calculations. Additionally, we perform a King plot analysis to set bounds on a fifth force in the keV/c 2 to MeV/c 2 range coupling to electrons and neutrons.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Optimally tuned starting point for single-shot GW calculations of solids

The dependence of ab initio many-body perturbation theory within the $\textit{GW}$ approximation on the eigensystem used in calculating quasiparticle corrections limits this method's predictive power. Here, we investigate the accuracy of the recently developed Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional as a generalized Kohn-Sham starting point for single-shot $\textit{GW}$ ($G_0W_0$) calculations for a range of semiconductors and insulators. Comparison to calculations based on well-established functionals, namely, PBE, PBE0, and HSE, as well as to self-consistent GW schemes and to experiment, shows that band gaps computed via $G_0W_0$@WOT-SRSH have a level of precision and accuracy that is comparable to that of more advanced methods such as quasiparticle self-consistent $\textit{GW}$ and eigenvalue self-consistent $\textit{GW}$. In this work, we also find that $G_0W_0$@WOT-SRSH improves the description of states deeper in the valence band manifold. Finally, we show that $G_0W_0$@WOT-SRSH significantly reduces the sensitivity of computed band gaps to ambiguities in the underlying WOT-SRSH tuning procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Application of the AMPLE cluster-and-truncate approach to NMR structures for molecular replacement

Processing of NMR structures for molecular replacement by AMPLE works well. AMPLE is a program developed for clustering and truncating ab initio protein structure predictions into search models for molecular replacement. Here, it is shown that its core cluster-and-truncate methods also work well for processing NMR ensembles into search models. Rosetta remodelling helps to extend success to NMR structures bearing low sequence identity or high structural divergence from the target protein. Potential future routes to improved performance are considered and practical, general guidelines on using AMPLE are provided.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

Metal atom clusters are studied based on the application of ab initio quantum mechanical approaches. Because these large 'molecular' systems pose special practical computational problems in the application of the quantum mechanical methods, there is a special need to find simplifying techniques that do not compromise the reliability of the calculations. Research is therefore directed towards various aspects of the implementation of the effective core potential technique for the removal of the metal atom core electrons from the calculations.

Kahn, L. R.↗