Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point

The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

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

Vacancy-Induced Tunable Kondo Effect in Twisted Bilayer Graphene

In single sheets of graphene, vacancy-induced states have been shown to host an effective spin-1/2 hole that can be Kondo screened at low temperatures. Here, we show how these vacancy-induced impurity states survive in twisted bilayer graphene (TBG), which thus provides a tunable system to probe the critical destruction of the Kondo effect in pseudogap hosts. Ab initio calculations and atomic-scale modeling are used to determine the nature of the vacancy states in the vicinity of the magic angle in TBG, demonstrating that the vacancy can be treated as a quantum impurity. Utilizing this insight, we construct an Anderson impurity model with a TBG host that we solve using the numerical renormalization group combined with the kernel polynomial method. We determine the phase diagram of the model and show how there is a strict dichotomy between vacancies in the AA/BB versus AB/BA tunneling regions. In AB/BA vacancies, the Kondo temperature at the magic angle develops a broad distribution with a tail to vanishing temperatures due to multifractal wave functions at the magic angle. Finally, we argue that scanning tunneling microscopy in the vicinity of the vacancy can act as a probe of both the critical single-particle states and the underlying many-body ground state in magic-angle TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The roles of three-nucleon force and continuum coupling in mirror symmetry breaking of oxygen mass region

With both three-nucleon force and continuum coupling included, we have developed a self-consistent ab initio Gamow shell model within the Gamow Hartree-Fock (GHF) basis obtained by the realistic interaction itself. With the chiral two-nucleon N 3 LO and three-nucleon N 2 LO interactions, the Gamow shell model has been applied to the mirror systems of Z = 8 neutron-rich isotopes and N = 8 proton-rich isotones, giving good agreements with data in binding energies, dripline positions and excitation spectra. The GHF calculated that the 0 ds/2 , 1 S1/2 and 1 p3/2 orbitals are resonances. The resonance states and their interplay with nonresonant continua play a crucial role in the descriptions of nuclei around driplines. Excitation spectra and Thomas-Ehrman shifts observed can be better described when both three-nucleon force and continuum coupling are considered in calculations. The three-nucleon force and continuum coupling produce a combined effect on the Thomas-Ehrman shift, e.g., for the 1/2 + resonance level of 19Na. The calculations help the understandings of related nuclear astrophysical processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Confinement-Controlled Rearrangements in Dioxolane Upgrading on H-ZSM-5 Revealed by Periodic DFT

Zeolitic Brønsted acid sites catalyze carbocation rearrangements central to upgrading biomass-derived oxygenates. Here we elucidate the mechanism of dioxolane conversion to methyl ethyl ketone and isobutanal on H-ZSM-5 using periodic density functional theory on the MFI model, complemented by ab initio molecular dynamics to probe confinement effects. Dioxolane adsorption at the Brønsted site is followed by protonation-assisted ring opening to form an oxocarbenium intermediate stabilized by the deprotonated framework. From this common intermediate, selectivity is governed by two competing rearrangements, namely, the 1,2-hydride shift with a free-energy barrier of 18.05 kcal mol –1 at 498 K leading toward MEK, and the 1,2-methyl shift with a higher barrier of 25.40 kcal mol –1 leading toward isobutanal. The hydride-shift channel is kinetically preferred over the methyl-shift channel, lowering the isobutanal/MEK ratio below the 3:1 limit expected for equal branching. Adsorption thermodynamics further indicate stronger stabilization of MEK than isobutanal within ZSM-5 channels, suggesting that confinement-controlled binding can bias product distributions in addition to intrinsic rearrangement barriers. These results highlight how Brønsted acidity and pore confinement jointly shape the rearrangement landscape in MFI zeolites.

adsorption↗
DOE OSTI2022

Atomically Resolved Electrically Active Intragrain Interfaces in Perovskite Semiconductors

Deciphering the atomic and electronic structures of interfaces is key to developing state-of-the-art perovskite semiconductors. However, conventional characterization techniques have limited previous studies mainly to grain-boundary interfaces, whereas the intragrain-interface microstructures and their electronic properties have been much less revealed. Herein using scanning transmission electron microscopy, we resolved the atomic-scale structural information on three prototypical intragrain interfaces, unraveling intriguing features clearly different from those from previous observations based on standalone films or nanomaterial samples. These intragrain interfaces include composition boundaries formed by heterogeneous ion distribution, stacking faults resulted from wrongly stacked crystal planes, and symmetrical twinning boundaries. The atomic-scale imaging of these intragrain interfaces enables us to build unequivocal models for the ab initio calculation of electronic properties. Our results suggest that these structure interfaces are generally electronically benign, whereas their dynamic interaction with point defects can still evoke detrimental effects. This work paves the way toward a more complete fundamental understanding of the microscopic structure-property-performance relationship in metal halide perovskites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2025

Universal Effective Charges in the 𝑠⁢𝑑 and 𝑓⁡𝑝 Shells

The 247-keV state in 54 Sc, populated in the 𝛽 decay of 54 Ca, is reported here as a nanosecond isomer with a half-life of 26.0(22) ns. The state is interpreted as the 1 + member of the 𝜋⁢𝑓 7/2 ⊗ 𝜈⁢𝑓 5/2 spin-coupled multiplet, which decays to the 3 + ,𝜋⁢𝑓 7/2 ⊗ 𝜈⁢𝑝 1/2 ground state. The new half-life corresponds to a pure 𝐸⁢2 transition with a strength of 1.93(16) W.u., providing the most precise, unambiguous 𝐵⁡(𝐸⁢2) value in the neutron-rich 𝑓⁡𝑝 region to date for a nucleus with valence protons above 𝑍 = 20. Notably, it is roughly 4 times larger than the 𝐵⁡(𝐸⁢2;1/2 − → 5/2 − ) value in 55 Ca. The results, as compared to semiempirical and ab initio shell-model calculations, indicate (1) a weak 𝑁 = 34 subshell gap relative to 𝑁 = 32, (2) a large 𝐸⁢2 enhancement in Sc as compared to Ca due to 1⁢𝑝 − 1⁢ℎ proton excitations across 𝑍 = 28, and (3) empirical effective proton and neutron charges 𝑒 𝜋 = 1.30⁢(8)⁢𝑒 and 𝑒 𝜈 = 0.452⁢(7)⁢𝑒, respectively, that are in contrast to reports of 𝑒 𝜋 ≈1.1–1.15⁢𝑒 and 𝑒 𝜈 ≈ 0.6–0.8⁢𝑒 for 𝑓⁡𝑝-shell nuclei near 𝑁 = 𝑍. Here, we demonstrate that these reports are erroneous and that, in fact, a universal set of effective charges can be used across the 𝑠⁢𝑑 and 𝑓⁡𝑝 shells.

Beta decay↗
DOE OSTI2021

Two-dimensional multifunctional materials from endohedral fullerenes

A multifunctional two-dimensional (2D) material is theoretically predicted based on systematic ab initio calculations and model simulations for the honeycomb lattice of endohedral fullerene W@C 28 molecules. It has structural bistability, ferroelectricity, multiple magnetic phases, and excellent valley characters and can be easily functionalized by the proximity effect with magnetic isolators such as MnTiO 3 . Furthermore, we may also manipulate the valley Hall and spin transport properties by selectively switching a few W@C 28 molecules to the metastable phase. These conclusions contribute to integrating different functions in a single 2D material for technological innovations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Application of an efficient generator-coordinate subspace-selection algorithm to neutrinoless double- β decay

The generator coordinate method begins with the variational construction of a set of nonorthogonal mean-field states that span a subspace of the full many-body Hilbert space. These states are then often projected onto states with good quantum numbers to restore symmetries, leading to a set with members that can be similar to one another, and it is sometimes possible to reduce this set without greatly affecting results. Here, we propose a greedy algorithm that we call the energy-transition-orthogonality procedure (ENTROP) to select subsets of important states. As applied here, the approach selects on the basis of diagonal energy, orthogonality, and contribution to the matrix element that governs neutrinoless double-β decay. We present both shell-model and preliminary ab initio calculations of this matrix element for the decay of 76Ge, with quadrupole deformation parameters and the isoscalar pairing strength as generator coordinates. ENTROP converges quickly, reducing significantly the number of basis states needed for an accurate calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Liquid-Liquid Transition in Water from First Principles

A long-standing question in water research is the possibility that supercooled liquid water can undergo a liquid-liquid phase transition (LLT) into high- and low-density liquids. Here, we used several complementary molecular simulation techniques to evaluate the possibility of an LLT in an ab initio neural network model of water trained on density functional theory calculations with the SCAN exchange correlation functional. We conclusively show the existence of a first-order LLT and an associated critical point in the SCAN description of water, representing the first definitive computational evidence for an LLT in water from first principles.

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

The development of thermodynamically consistent and physics-informed equation-of-state model through machine learning

Ab initio molecular dynamics (AIMD) simulations have become an important tool used in the construction of equations of state (EOS) tables for warm dense matter. Due to computational costs, only a limited number of system state conditions can be simulated, and the remaining EOS surface must be interpolated for use in radiation-hydrodynamic simulations of experiments. In this work, we develop a thermodynamically consistent EOS model that utilizes a physics-informed machine learning approach to implicitly learn the underlying Helmholtz free-energy from AIMD generated energies and pressures. The model, referred to as PIML-EOS, was trained and tested on warm dense polystyrene producing a fit within a 1% relative error for both energy and pressure and is shown to satisfy both the Maxwell and Gibbs–Duhem relations. In addition, we provide a path toward obtaining thermodynamic quantities, such as the total entropy and chemical potential (containing both ionic and electronic contributions), which are not available from current AIMD simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Utilizing Amino Acid Composition and Entropy of Potential Open Reading Frames to Identify Protein-Coding Genes

One of the main steps in gene-finding in prokaryotes is determining which open reading frames encode for a protein, and which occur by chance alone. There are many different methods to differentiate the two; the most prevalent approach is using shared homology with a database of known genes. This method presents many pitfalls, most notably the catch that you only find genes that you have seen before. The four most popular prokaryotic gene-prediction programs (GeneMark, Glimmer, Prodigal, Phanotate) all use a protein-coding training model to predict protein-coding genes, with the latter three allowing for the training model to be created ab initio from the input genome. Different methods are available for creating the training model, and to increase the accuracy of such tools, we present here GOODORFS, a method for identifying protein-coding genes within a set of all possible open reading frames (ORFS). Our workflow begins with taking the amino acid frequencies of each ORF, calculating an entropy density profile (EDP), using KMeans to cluster the EDPs, and then selecting the cluster with the lowest variation as the coding ORFs. To test the efficacy of our method, we ran GOODORFS on 14,179 annotated phage genomes, and compared our results to the initial training-set creation step of four other similar methods (Glimmer, MED2, PHANOTATE, Prodigal). We found that GOODORFS was the most accurate (0.94) and had the best F1-score (0.85), while Glimmer had the highest precision (0.92) and PHANOTATE had the highest recall (0.96).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Nanoscale fluctuation of stacking fault energy strengthens multi-principal element alloys

Chemical randomness and the associated energy fluctuation are essential features of multi-principal element alloys (MPEAs). Due to these features, nanoscale stacking fault energy (SFE) fluctuation is a natural and independent contribution to strengthening MPEAs. However, existing models for conventional alloys (i.e., alloys with one principal element) cannot be applied to MPEAs. The extreme values of SFEs required by such models are unknown for MPEAs, which need to calculate the nanoscale volume relevant to the SFE fluctuation. In the present work, we developed an analytic model to evaluate the strengthening effect through the SFE fluctuation, profuse in MPEAs. The model has no adjustable parameters, and all parameters can be determined from experiments and ab initio calculations. Further, this model explains available experimental observations and provides insightful guidance for designing new MPEAs based on the SFE fluctuation. It generally applies to MPEAs in random states and with chemical short-range order.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Prediction of exotic magnetic states in the alkali-metal quasi-one-dimensional iron selenide compound Na 2 FeSe 2

The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali-metal iron selenide compound Na 2 FeSe 2 is presented. The novelty of this material is that the valence of iron is Fe 2+ , contrary to most other iron-chain compounds with valence Fe 3+ . Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the ab initio band structure near the Fermi level. Including Hubbard and Hund couplings and studying the model via the density-matrix renormalization group and Lanczos methods, we constructed the ground-state phase diagram. A robust region where the block state ↑↑↓↓↑↑↓↓ is stabilized was unveiled. The analog state in iron ladders, employing 2×2 ferromagnetic blocks, is by now well established, but in chains a block magnetic order has not been observed yet in real materials. The phase diagram also contains a large region of canonical staggered spin order ↑↓↑↓↑↓↑ at very large Hubbard repulsion. At the block-to-staggered transition region, an exotic phase is stabilized with a mixture of both states: an inhomogeneous orbital-selective charge density wave with the exotic spin configuration ↑↑↓↑↓↓↑↓. Our predictions for Na 2 FeSe 2 may guide crystal growers and neutron-scattering experimentalists towards the realization of block states in one-dimensional iron selenide chain materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗