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 955 records · Page 53

Charge-transfer energy in closed-shell ion-atom interactions

The importance of charge-transfer energy in the interactions between closed-shell ions and atoms is investigated. Ab initio calculations on H(plus)-He and Li(plus)-He are used as a guide for the construction of approximate methods for the estimation of the charge-transfer energy for more complicated systems. For many alkali ion-rate gas systems the charge-transfer energy is comparable to the induction energy in the region of the potential minimum, although for doubly charged alkaline-earth ions in rare gases the induction energy always dominates. Surprisingly, an empirical combination of repulsion energy plus asymptotic induction energy plus asymptotic dispersion energy seems to give a fair representation of the total interaction, especially if the repulsion energy is parameterized, despite the omission of any explicit charge-transfer contribution. More refined interaction models should consider the charge-transfer energy contribution.

Alvarez-Rizzatti, M.↗

The Vaporization of B2O3(l) to B2O3(g) and B2O2(g)

The vaporization of B2O3 in a reducing environment leads to formation of both B2O3(g) and B2O2(g). While formation of B2O3(g) is well understood, many questions about the formation of B2O2(g) remain. Previous studies using B(s) + B2O3(l) have led to inconsistent thermodynamic data. In this study, it was found that after heating, B(s) and B2O3(l) appear to separate and variations in contact area likely led to the inconsistent vapor pressures of B2O2(g). To circumvent this problem, an activity of boron is fixed with a two-phase mixture of FeB and Fe2B. Both second and third law enthalpies of formation were measured for B2O2(g) and B2O3(g). From these the enthalpies of formation at 298.15 K are calculated to be -479.9 +/- 41.5 kJ/mol for B2O2(g) and -833.4 +/- 13.1 kJ/mol for B2O3(g). Ab initio calculations to determine the enthalpies of formation of B2O2(g) and B2O3(g) were conducted using the W1BD composite method and show good agreement with the experimental values.

Jacobson, Nathan S.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An ab initio study of multiple phonon scattering resonances in silicon germanium alloys

We have computed here phonon scattering rates and density of states in silicon germanium alloys using Green's function calculations and density functional theory. This method contrasts with the virtual crystal approximation (VCA) used in conjunction with Fermi's golden rule, which cannot capture resonance states occurring through the interaction of substitutional impurities with the host lattice. These resonances are demonstrated by density of states and scattering rate calculations in the dilute limit and show broadening as the concentration increases. Although these deviations become significant from the VCA at high frequencies, the relaxation times obtained for these phonon modes are small in both the full scattering theory and the VCA, resulting in their negligible contribution to thermal transport.

36 MATERIALS SCIENCE↗

The Renner-Teller effect in HCCCl{sup +}(X{sup ~2}Π) studied by zero-kinetic energy photoelectron spectroscopy and ab initio calculations

The spin-vibronic energy levels of the chloroacetylene cation up to 4000 cm{sup −1} above the ground state have been measured using the one-photon zero-kinetic energy photoelectron spectroscopic method. The spin-vibronic energy levels have also been calculated using a diabatic model, in which the potential energy surfaces are expressed by expansions of internal coordinates, and the Hamiltonian matrix equation is solved using a variational method with harmonic basis functions. The calculated spin-vibronic energy levels are in good agreement with the experimental data. The Renner-Teller (RT) parameters describing the vibronic coupling for the H—C≡C bending mode (ε{sub 4}), Cl—C≡C bending mode (ε{sub 5}), the cross-mode vibronic coupling (ε{sub 45}) of the two bending vibrations, and their vibrational frequencies (ω{sub 4} and ω{sub 5}) have also been determined using an effective Hamiltonian matrix treatment. In comparison with the spin-orbit interaction, the RT effect in the H—C≡C bending (ε{sub 4}) mode is strong, while the RT effect in the Cl—C≡C bending mode is weak. There is a strong cross-mode vibronic coupling of the two bending vibrations, which may be due to a vibronic resonance between the two bending vibrations. The spin-orbit energy splitting of the ground state has been determined for the first time and is found to be 209 ± 2 cm{sup −1}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio phonon thermal transport in monolayer InSe, GaSe, GaS, and alloys

We compare vibrational properties and phonon thermal conductivities (κ) of monolayer InSe, GaSe and GaS systems using density functional theory and Peierls-Boltzmann transport methods. In going from InSe to GaSe to GaS, system mass decreases giving both increasing acoustic phonon velocities and decreasing scattering of these heat-carrying modes with optic phonons, ultimately giving κInSe< κGaSe< κGaS. This behavior is demonstrated by correlating the scattering phase space limited by fundamental conservation conditions with mode scattering rates and phonon dispersions for each material. We also show that, unlike flat monolayer systems such as graphene, thermal transport is governed by in-plane vibrations in InSe, GaSe and GaS, similar to buckled monolayer materials such as silicene. Alloying of InSe, GaSe and GaS systems provides an effective method for modulating their κ through intrinsic vibrational modifications and phonon scattering from mass disorder giving reductions ~2-3.5 times. This disorder also suppresses phonon mean free paths in the alloy systems compared to those in their crystalline counterparts. This work provides fundamental insights of lattice thermal transport from basic vibrational properties for an interesting set of two-dimensional materials.

Pandey, Tribhuwan (ORCID:0000000203737716)↗

Ab initio electroweak corrections to superallowed β decays and their impact on V ud

Radiative corrections are essential for an accurate determination of V ud from superallowed β decays. In view of recent progress in the single-nucleon sector, the uncertainty is dominated by the theoretical description of nucleus-dependent effects, limiting the precision that can currently be achieved for V ud . In this work, we provide a detailed account of the electroweak corrections to superallowed β decays in effective field theory (EFT), including the power counting, potential and ultrasoft contributions, and factorization in the decay rate. Here, we present a first numerical evaluation of the dominant corrections in light nuclei based on quantum Monte Carlo methods, confirming the expectations from the EFT power counting. Finally, we discuss strategies how to extract from data the low-energy constants that parametrize short-distance contributions and whose values are not predicted by the EFT. Combined with advances in ab initio nuclear-structure calculations, this EFT framework allows one to systematically address the dominant uncertainty in V ud , as illustrated in detail for the 14 O → 14 N transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Bergman-type medium range order in amorphous Zr77Rh23 alloy studied by ab initio molecular dynamics simulations

In recent years, some arguments about the existence of medium-range order (MRO) in the Zr-Rh system have been put forward. However, research on the structural features of the Zr-Rh binary alloy at the atomic level is still lacking. This study uses ab initio molecular dynamics simulations to systematically study the local structures of Zr77Rh23 from the liquid to the glassy states. Pair correlation function (PCF), coordination number (CN), Honeycutt–Anderson(HA) index, bond-angle distribution functions, and the Voronoi tessellation method are used to reveal a clear icosahedral-like configuration in the amorphous Zr77Rh23 alloy. It is noteworthy that the splitting in the second peak of the partial PDF implies the existence of a medium range order (MRO) in the Zr77Rh23 system. We obtain the local order in three-dimensional atomic density distributions by using a new atomistic cluster alignment (ACA) method. Interestingly, a Bergman-type MRO is observed in the glassy Zr77Rh23. Furthermore, the spatial distribution and connections of icosahedral-like clusters are shown to further demonstrate the MRO network. Our findings shed light on the nature of atomic local structures of amorphous Zr77Rh23 alloy and have important implications to understanding the formation of various MROs in metallic glasses.

Guo, Y. R.↗