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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

On the applicability of one- and many-electron quantum chemistry models for hydrated electron clusters

We evaluate the applicability of a hierarchy of quantum models in characterizing the binding energy of excess electrons to water clusters. In particular, we calculate the vertical detachment energy of an excess electron from water cluster anions with methods that include one-electron pseudopotential calculations, density functional theory (DFT) based calculations, and ab initio quantum chemistry using MP2 and eom-EA-CCSD levels of theory. The examined clusters range from the smallest cluster size (n = 2) up to nearly nanosize clusters with n = 1000 molecules. The examined cluster configurations are extracted from mixed quantum-classical molecular dynamics trajectories of cluster anions with n = 1000 water molecules using two different one-electron pseudopotenial models. We find that while MP2 calculations with large diffuse basis set provide a reasonable description for the hydrated electron system, DFT methods should be used with precaution and only after careful benchmarking. Strictly tested one-electron psudopotentials can still be considered as reasonable alternatives to DFT methods, especially in large systems. The results of quantum chemistry calculations performed on configurations, that represent possible excess electron binding motifs in the clusters, appear to be consistent with the results using a cavity structure preferring one-electron pseudopotential for the hydrated electron, while they are in sharp disagreement with the structural predictions of a non-cavity model.

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

Homochiral and racemic MicroED structures of a peptide repeat from the ice-nucleation protein InaZ

The ice-nucleation protein InaZ from Pseudomonas syringae contains a large number of degenerate repeats that span more than a quarter of its sequence and include the segment GSTSTA. Ab initio structures of this repeat segment, resolved to 1.1 Å by microfocus X-ray crystallography and to 0.9 Å by the cryo-EM method MicroED, were determined from both racemic and homochiral crystals. The benefits of racemic protein crystals for structure determination by MicroED were evaluated and it was confirmed that the phase restriction introduced by crystal centrosymmetry increases the number of successful trials during the ab initio phasing of the electron diffraction data. Both homochiral and racemic GSTSTA form amyloid-like protofibrils with labile, corrugated antiparallel β-sheets that mate face to back. The racemic GSTSTA protofibril represents a new class of amyloid assembly in which all-left-handed sheets mate with their all-right-handed counterparts. This determination of racemic amyloid assemblies by MicroED reveals complex amyloid architectures and illustrates the racemic advantage in macromolecular crystallography, now with submicrometre-sized crystals.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Evolution of short- and medium-range order in the melt-quenching amorphization of Ge 2 Sb 2 Te 5

Phase-change memory takes advantage of the fast phase transition between amorphous and crystalline phases of phase-change materials (e.g.,Ge 2 Sb 2 Te 5 or GST). To date, while the “SET” process (crystallization of GST glass) has been intensively studied, studies on the “RESET” process (melt-quenching amorphization of GST) are still limited. Here, we explore the structural changes of GST upon rapid cooling by ab initio molecular dynamics simulations and atomistic cluster alignment (ACA) analysis. Different from other methods which only focus on the nearest bonding atoms, the ACA method can study both the short- and medium-range order clusters containing atoms beyond the first-neighboring shell and enables us to explore the changes of cluster structures in a larger scale. The results reveal that low-coordinated octahedral clusters tend to become high-coordinated ones, and Ge-centered octahedral structures change to tetrahedrons whereas Sb-centered tetrahedrons transform to octahedral structures during the amorphization process. Interestingly, tetrahedrons show aggregation in liquid and supercooled liquid in contrast to 6-fold octahedrons which present notable aggregation in amorphous GST. Furthermore, our study showed that wrong bonds (Ge–Ge, Sb–Sb, Ge–Sb and Te–Te bonds) can promote the formation of large rings, and irreducible rings tend to separate into smaller and larger rings as the temperature is decreased. Our conclusions provide useful insights into the formation process and the structure of amorphous GST, which is valuable for facilitating the application of phase change materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Robust determination of the chemical potential in the pole expansion and selected inversion method for solving Kohn-Sham density functional theory

Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving Kohn-Sham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method. The main idea of the new method is not to find the exact chemical potential at each self-consistent-field (SCF) iteration but to dynamically and rigorously update the upper and lower bounds for the true chemical potential, so that the chemical potential reaches its convergence along the SCF iteration. Instead of evaluating the Fermi operator for multiple times sequentially, our method uses a two-level strategy that evaluates the Fermi operators in parallel. In the regime of full parallelization, the wall clock time of each SCF iteration is always close to the time for one single evaluation of the Fermi operator, even when the initial guess is far away from the converged solution. We demonstrate the effectiveness of the new method using examples with metallic and insulating characters, as well as results from ab initio molecular dynamics.

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

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. Further, the violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optically controlled single-valley exciton doublet states with tunable internal spin structures and spin magnetization generation

Manipulating quantum states through light–matter interactions has been actively pursued in two-dimensional materials research. Significant progress has been made toward the optical control of the valley degrees of freedom in semiconducting monolayer transition-metal dichalcogenides, based on doubly degenerate excitons from their two distinct valleys in reciprocal space. Here, we introduce a type of optically controllable doubly degenerate exciton states that come from a single valley, dubbed as single-valley exciton doublet (SVXD) states. They are unique in that their constituent holes originate from the same valence band, making possible the direct optical control of the spin structure of the excited constituent electrons. Combining ab initio GW plus Bethe–Salpeter equation ( GW -BSE) calculations and a theoretical analysis method, we demonstrate such SVXD in substrate-supported monolayer bismuthene—which has been successfully grown using molecular beam epitaxy. In each of the two distinct valleys in the Brillouin zone, strong spin–orbit coupling and C 3 v symmetry lead to a pair of degenerate 1s exciton states (the SVXD states) with opposite spin configurations. Any coherent linear combinations of the SVXD in a single valley can be excited by light with a specific polarization, enabling full manipulation of their internal spin configurations. In particular, a controllable net spin magnetization can be generated through light excitation. Our findings open routes to control quantum degrees of freedom, paving the way for applications in spintronics and quantum information science.

2D materials↗
DOE OSTI2025

General features of the stellar matter equation of state from microscopic theory, new maximum-mass constraints, and causality

The profile of a neutron star probes a very large range of densities, from the density of iron up to several times the density of saturated nuclear matter, and thus no theory of hadrons can be considered reliable if extended to those regions. We emphasize the importance of taking contemporary ab initio theories of nuclear and neutron matter as the baseline for any extension method, which will unavoidably involve some degree of phenomenology. We discuss how microscopic theory, on the one end, with causality and maximum-mass constraints, on the other, set strong boundaries to the high-density equation of state. We present our latest neutron star predictions where we combine polytropic extensions and parametrizations guided by speed of sound considerations. The predictions we show include our baseline neutron star cooling curves.

chiral effective field theory↗
NASA NTRSDate not supplied

Overview of NASA Research Activities in Shock Layer Kinetics and Radiation

Shock layer radiation is an important heating mechanism for entry probes to most planetary destinations and re-entry to Earth from beyond low Earth orbit. Our understanding of shock layer radiation phenomena has improved tremendously over the last decade thanks to NASA investment in fundamental radiation research through Entry Systems Modeling’s (ESM) Shock Layer Kinetics and Radiation (SLKR) task. This talk will overview the recent activities within SLKR, including validation of models through ground testing, flight instrumentation and remote observations; development of first principles ab initio calculations for reaction mechanisms and spectroscopic databases and advancing numerical and computational methods for prediction of radiation in reacting hypersonic flows. 1 Sr. Research Scientist, Aerothermodynamics Branch, and AIAA Associate Fellow.

Brett A. Cruden↗
DOE OSTI2019

BaZrSe 3 : Ab initio study of anion substitution for bandgap tuning in a chalcogenide material

Recently, transition metal perovskite chalcogenide materials have been introduced as possible candidates for solar cell applications. Here, we provide accurate theoretical calculations for BaZrS 3 and two phases of SrZrS 3 , which have been recently synthesized and their optical properties elaborated. In this study, we consider the substitution of S in BaZrS 3 with Se to form BaZrSe 3 . Evolutionary methods are used to find the optimal structure of this compound, and accurate calculations of its optoelectronic properties are presented. Using phonon frequency calculations and ab initio molecular dynamics, we assess the stability of this compound. We find that BaZrSe 3 is likely to be stable under typical conditions, with a low bandgap and high optical absorption coefficients. This indicates that BaZrSe 3 could be useful for solar cell applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Ab initio calculation of the cross sections for electron impact vibrational excitation of CO via the {sup 2}Π shape resonance

The stabilization method is used to calculate the complex potential energy curve of the {sup 2}Π state of CO{sup −} as a function of bond length, with the refinement that separate potentials are determined for p-wave and d-wave attachment and detachment of the excess electron. Using the resulting complex potentials, absolute vibrational excitation cross sections are calculated as a function of electron energy and scattering angle. The calculated cross sections agree well with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

An Accurate Quartic Force Field and Vibrational Frequencies for HNO and DNO

An accurate ab initio quartic force field for HNO has been determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with the correlation consistent polarized valence triple zeta (cc-pVTZ) basis set. Improved harmonic frequencies were determined with the cc-pVQZ basis set. Fundamental vibrational frequencies were determined using a second-order perturbation theory analysis and also using variational calculations. The N-0 stretch and bending fundamentals are determined well from both vibrational analyses. The H-N stretch, however, is shown to have an unusually large anharmonic correction, and is not well determined using second-order perturbation theory. The H-N fundamental is well determined from the variational calculations, demonstrating the quality of the ab initio quartic force field. The zero-point energy of HNO that should be used in isodesmic reactions is also discussed.

Dateo, Christopher E.↗
DOE OSTI2025

Quantifying uncertainties in α -nucleus reaction dynamics informed from first principles

The ab initio symmetry-adapted no-core shell model is a microscopic many-body method which naturally describes challenging collective and clustering features of atomic nuclei. Wave functions and observables computed with realistic nucleon-nucleon forces in this framework are tied to first principles, and are hence well-suited for rigorous uncertainty quantification. We discuss α-deuteron and α- 12 C cluster potentials informed by symmetry-adapted calculations, and propagate uncertainties in the effective binary cluster method as well as those in the nuclear interaction to reaction observables, namely scattering phase shifts, cross sections, partial widths and resonance energies. Here, we find that the overall uncertainties are dominated by those originating in the underlying nuclear force, speaking to the need for tighter constraints on realistic nucleon-nucleon interactions.

Ab initio↗
DOE OSTI2015

An accurate potential energy surface for the F + H{sub 2} → HF + H reaction by the coupled-cluster method

A three dimensional potential energy surface for the F + H{sub 2} → HF + H reaction has been computed by the spin unrestricted coupled cluster method with singles, doubles, triples, and perturbative quadruples [UCCSDT(2){sub Q}] using the augmented correlation-consistent polarised valence quadruple zeta basis set for the fluorine atom and the correlation-consistent polarised valence quadruple zeta basis set for the hydrogen atom. All the calculations are based on the restricted open-shell Hartree-Fock orbitals, together with the frozen core approximations, and the UCCSD(T)/complete basis set (CBS) correction term was included. The global potential energy surface was calculated by fitting the sampled ab initio points without any scaling factor for the correlation energy part using a neutral network function method. Extensive dynamics calculations have been carried out on the potential energy surface. The reaction rate constants, integral cross sections, product rotational states distribution, and forward and backward scattering as a function of collision energy of the F + HD → HF + D, F + HD → DF + H, and F + H{sub 2} reaction, were calculated by the time-independent quantum dynamics scattering theory using the new surface. The satisfactory agreement with the reported experimental observations previously demonstrates the accuracy of the new potential energy surface.

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

Influence of near-surface oxide layers on TiFe hydrogenation: mechanistic insights and implications for hydrogen storage applications

The inevitable formation of passivating oxide films on the surface of the TiFe intermetallic compound limits its performance as a stationary hydrogen storage material. Extensive experimental efforts have been dedicated to the activation of TiFe, i.e. oxide layer removal prior to utilization for hydrogen storage. However, development of an efficient activation protocol necessitates a fundamental understanding of the composition and structure of the air-exposed surface and its interaction with hydrogen, which is currently absent. Therefore, in this study we explored the growth and nature of oxide films on the most exposed TiFe surface (110) in depth using static and dynamic first-principles methods. We identified the lowest energy structures for six oxygen coverages up to approximately 1.12 nm of thickness with a global optimization method and studied the temperature effects and structural evolution of the oxide phases in detail via ab initio molecular dynamics (AIMD). Based on structural similarity and coordination analysis, motifs for TiO 2 and TiFeO 3 as well as Ti(FeO 2 ) x (x = 2, 3 or 5) phases were identified. On evaluating the interaction of the oxidized surface with hydrogen, a minimal energy barrier of 0.172 eV was predicted for H 2 dissociation while H migration from the top of the oxidized surface to the bulk TiFe was limited by several high-lying energy barriers above 1.4 eV. Our mechanistic insights will prove themselves valuable for informed designs towards new activation methods of TiFe and related systems as hydrogen storage materials.

08 HYDROGEN↗
DOE OSTI2023

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Investigation of the Thermodynamic and Electronic Properties of Double Perovskite Ca{sub 2}CoNbO{sub 6}

In the present work, a self-consistent ab initio calculation using the full-potential linearized augmented plane wave (FP-LAPW) method within the framework of the spin-polarized density functional theory (DFT) was used to study the structural, electronic, and thermodynamic properties of Ca{sub 2}CoNbO{sub 6} double perovskite compound. The generalized gradient approximation (GGA) described by Perdew–Burke–Ernzerhof (PBE) and GGA + U were used. The results obtained for the electronic properties show a ferrimagnetic and half-metallic behavior of the compound. The novelty of our work is the study of the thermodynamic properties of Ca{sub 2}CoNbO{sub 6} double perovskite such as heat capacity and Debye temperature which showed an important effect of pressure compared to the temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Effect of strain along C-axis NbS{sub 2}

We have studied electronic properties of double layered hexagonal structure of the Niobium Di-Sulphide (2H-NbS{sub 2}) superconductor for various strains introduced along the c-axis using ab-initio calculations. The DFT calculations based on Full Potential Linearized Augmented Plane Wave (FPLAPW) method are performed using the ELK code. The total energy curve (E vs a), Density of States (DOS) and the Band structure calculations obtained in this work are matching with the earlier reports. The Pressure-Volume (P-V) diagram for 2H-NbS{sub 2} was obtained using the Equation of State(EOS) calculations, which provides the relationship between the pressure and strain applied along the c-axis. The band structures for various strains ranging from 0 percent to 10 percent along c-axis in steps of 2 percent are obtained. We note that there are increasing number of bands crossing over the Fermi energy level with increase in strain. Thus, we conclude that with increasing strain along c-axis, number of conduction bands crossing the E{sub F} increases, which gives rise to more conduction states and hence higher conductivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗