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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ab Initio Structures and Energetics of Hydrated Flat and Terrace-Step Surfaces of Forsterite (Mg 2 SiO 4 )

Forsterite (Mg 2 SiO 4 ), a model divalent metal silicate mineral, has been extensively studied in the context of mineral carbonation. Although dissolution is a key step in this process, the mechanisms by which forsterite dissolves under high CO 2 conditions remain poorly understood. Atomistic simulations could aid in exploring these mechanisms, but it is essential first to understand the structures and energetics of the relevant forsterite surfaces. We present an ab initio study of the structure and surface energy at 0 K of the flat $(010), (110), (001), (111), (021), (101)$ and $(120)$ faces of forsterite using the density functional PBE Hamiltonian and a plane-wave basis set. Dry surfaces became stabilized upon hydration through the formation of bonds between surface Mg and O from water, as well as by the formation of hydrogen bonds. According to surface energy values, the stability order of the hydrated forsterite faces was found to be $(120) < (101) < (021) < (111) < (001) < (110) < (010)$. We also investigated the energetics of the terrace-step $(0\bar{41})$ surface as a model site for forsterite dissolution. Among all the facets, the $(0\bar{41})$ surface is the least stable termination in water. Hydration of Mg atoms on the $(0\bar{41})$ surface increases their susceptibility to dissolution. The presence of a step and its hydration destabilizes the terraces, making step retreat more likely than a dissolution front advancing along the [010] direction. This research will support future simulations to investigate forsterite dissolution in water under CO 2 -rich conditions.

PBE Hamiltonian↗
DOE OSTI2023

Determination of the vacancy distribution over Al cation sites in γ-Al 2 O 3

Although gamma-alumina (γ-Al 2 O 3 ) is an extensively used material with wide-ranging applications due to its inherently high surface area and acidity, its atomic structure is still not fully understood. γ-Al 2 O 3 is described as having a spinel-like structure, where the O sublattice has a face-centered cubic (FCC) arrangement and Al cations are placed in the spinel tetrahedral and octahedral interstitial sites. Achieving the correct stoichiometry of Al 2 O 3 , however, requires the introduction of Al vacancies into some of the interstitial sites. Despite the importance of accurately describing the structure of γ-Al 2 O 3 , the distribution of vacancies between tetrahedral and octahedral sites remains unclear, in part because of the usually poor crystalline quality of γ-Al 2 O 3 that has often been used in previous studies. To determine the actual cation distribution in γ-Al 2 O 3 , single-crystalline γ-Al 2 O 3 was investigated using a correlative approach of experimental and simulated selected-area electron diffraction (SAED) and high-resolution electron energy-loss spectroscopy (EELS). Comparison of the reflection intensities in single-crystal SAED to simulated SAED from models with varied vacancy distributions revealed that vacancies exist primarily on tetrahedral sites, contrary to the placement of vacancies on octahedral sites proposed in several common models. Comparison of EELS spectra—acquired with the highest energy resolution reported so far for γ-Al 2 O 3 —with ab initio multiple scattering EELS simulations confirmed the distribution of vacancies on tetrahedral sites. These results enable more accurate modeling of γ-Al 2 O 3 to better predict its properties in existing and future applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluation of Density-Functional Tight-Binding Methods for Simulation of Protic Molecular Ion Pairs

In this work, we benchmark the accuracy of the density-functional tight-binding (DFTB) method, namely the long-range corrected second-order (LC-DFTB2) and third-order (DFTB3) models, for predicting energetics of imidazolium-based ionic liquid (IL) ion pairs. We compare the DFTB models against popular density functionals such as LC-ωPBE and B3LYP, using ab initio domain-based local pair-natural orbital coupled cluster (DLPNO-CC) energies as reference. Calculations were carried out in the gas phase, as well as in aqueous solution using implicit solvent methods. We find that the LC-DFTB2 model shows excellent performance in the gas phase and agrees well with reference energies in implicit solvent, often outperforming DFTB3 predictions for complexation energetics. Our study identifies a range of opportunities for use of the LC-DFTB method and quantifies its sensitivity to protonation states and the types of chemical interactions between ion pairs.

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

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio benchmarks of neutrinoless double- β decay in light nuclei with a chiral Hamiltonian

Here, we report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-β (0νββ) decays in light nuclei with mass numbers ranging from A = 6 to A = 22. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two plus three nucleon interaction from chiral effective field theory is employed, and both isospin-conserving (ΔT = 0) and isospin-changing (ΔT = 2) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of ΔT = 0 transitions are in good agreement among all calculations, at the level of 10%. For ΔT = 2, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps toward improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the 0νββ decays in heavier candidate nuclei, where quasiexact benchmarks are not available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2024

Anomalous Interlayer Exciton Diffusion in WS 2 /WSe 2 Moiré Heterostructure

Stacking van der Waals crystals allows for the on-demand creation of a periodic potential landscape to tailor the transport of quasiparticle excitations. We investigate the diffusion of photoexcited electron-hole pairs, or excitons, at the interface of WS 2 /WSe 2 van der Waals heterostructure over a wide range of temperatures. We observe the appearance of distinct interlayer excitons for parallel and antiparallel stacking and track their diffusion through spatially and temporally resolved photoluminescence spectroscopy from 30 to 250 K. While the measured exciton diffusivity decreases with temperature, it surprisingly plateaus below 90 K. Our observations cannot be explained by classical models like hopping in the moiré potential. A combination of ab initio theory and molecular dynamics simulations suggests that low-energy phonons arising from the mismatched lattices of moiré heterostructures, also known as phasons, play a key role in describing and understanding this anomalous behavior of exciton diffusion. Our observations indicate that the moiré potential landscape is dynamic down to very low temperatures and that the phason modes can enable efficient transport of energy in the form of excitons.

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

Density of states prediction for materials discovery via contrastive learning from probabilistic embeddings

Abstract Machine learning for materials discovery has largely focused on predicting an individual scalar rather than multiple related properties, where spectral properties are an important example. Fundamental spectral properties include the phonon density of states (phDOS) and the electronic density of states (eDOS), which individually or collectively are the origins of a breadth of materials observables and functions. Building upon the success of graph attention networks for encoding crystalline materials, we introduce a probabilistic embedding generator specifically tailored to the prediction of spectral properties. Coupled with supervised contrastive learning, our materials-to-spectrum (Mat2Spec) model outperforms state-of-the-art methods for predicting ab initio phDOS and eDOS for crystalline materials. We demonstrate Mat2Spec’s ability to identify eDOS gaps below the Fermi energy, validating predictions with ab initio calculations and thereby discovering candidate thermoelectrics and transparent conductors. Mat2Spec is an exemplar framework for predicting spectral properties of materials via strategically incorporated machine learning techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The O K -2 V spectrum of CO: the influence of the second core-hole

Using synchrotron radiation in the tender X-ray regime, a photoelectron spectrum showing the formation of single site double-core-hole pre-edge states, involving the K shell of the O atom in CO, has been recorded by means of high-resolution electron spectroscopy. The experimentally observed structures have been simulated, interpreted and assigned, employing state-of-the-art ab initio quantum chemical calculations, on the basis of a theoretical model, accounting for their so-called direct or conjugate character. Features appearing above the double ionization threshold have been reproduced by taking into account the strong mixing between multi-excited and continuum states. The shift of the σ* resonance below the double ionization threshold, in combination with the non-negligible contributions of multi-excited configurations in the final states reached, gives rise to a series of avoided crossings between the different potential energy curves.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The effect of rare-earth substitution on the Debye temperature of inorganic phosphors

In the quest to predict photoluminescent efficiency in rare-earth substituted inorganic phosphors, research has shown that materials consisting of ordered, rigid crystal structures tend to possess the highest photoluminescent quantum yields. A compound's Debye temperature (Θ D ), which can be calculated using ab initio calculations, is an ideal proxy for quantitatively comparing structural rigidity among different inorganic compounds, allowing potentially efficient phosphors to be selected from large crystal structure databases. However, the high computational cost of these calculations limits estimating Θ D for unsubstituted host crystal structures only. It is assumed that the low substitution concentration of the rare-earth luminescent center does not significantly influence a material's Debye temperature. This work evaluates the validity of this approximation by examining the effect of luminescent center substitution on a host structure's Θ D . Two well-known phosphors, (Y 1– x Ce x ) 3 Al 5 O 12 ( x = 0 – 0.05) and Ba 1-δ Eu δ MgAl 10 O 17 (δ = 0 – 0.15), were synthesized with varying rare-earth concentrations, while Θ D was computationally estimated and then determined by ultrasonic pulse-echo speed-of-sound and low-temperature heat capacity measurements. Overall, the ensuing results provide key implications for using Θ D as a proxy for structural rigidity in substituted inorganic compounds.

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

Efficient simulation of optical spectra via machine learning and physical decomposition of environmental effects

Simulations of optical spectra can provide key insights to aid experimental interpretation of electronic excitation phenomena. For chromophores in the condensed phase, these spectra, which incorporate the coupling between electronic excitation and molecular and solvent nuclear motions, can be simulated using excitation energies obtained from molecular dynamics simulations of the chromophore and solvent. Here, we present a hybrid scheme that exploits machine learning and physically informed spectral densities to show that as few as 25 ground and excited state energetic gradient calculations can be used to construct models that accurately predict environment-influenced vibronic coupling in optical spectra. We demonstrate our approach for the green fluorescent protein chromophore in water and the cresyl violet chromophore in methanol. We show that our hybrid approach, employing a machine learning model for the high-frequency spectral density and an ab initio parameterized Debye spectral density for the low-frequency, results in a systematic improvement of the optical absorption lineshape, leading to a simple machine learning scheme that can be used for the simulation of spectral densities and optical spectra.

Snider, Andrew [Univ. of California, Merced, CA (U↗
DOE OSTI2025

Exploiting 20 Ne Isotopes for Precision Characterizations of Collectivity in Small Systems

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174⁢(8) stat ⁢(31) syst for NLEFT and 1.139⁢(6) stat ⁢(39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Characterizing the nonequilibrium response of FeRh thin films using time-domain thermoreflectance

Time-domain thermoreflectance (TDTR) characterization of FeRh throughout its first-order antiferromagnetic (AF) to ferromagnetic (FM) transition shows that the transient reflectance ΔR⁡(t)/R strongly depends on the magnetic order of the sample. Using TDTR, which uses optical pulses to induce small temperature excursions, we have found that Δ⁢R⁡(t)/R of the AF phase exhibits a large negative response, while the response of the FM phase is positive. This magnetic phase sensitivity has allowed us to study the transient response of both the AF and FM phases to the pump-pulse excitation and the mixed phase of the material. These results are significant since the ultrafast properties of antiferromagnetic materials and mixed antiferromagnetic and ferromagnetic materials are difficult to detect using other conventional techniques. We have found that the AF phase exhibits a strong subpicosecond decaying signal not observed in the FM phase. The magnetic phase dependence of the sign of Δ⁢R⁡(t)/R is qualitatively explained using the results of ab initio density functional theory calculations. Using the two-temperature model, we found that the change in the thermalization time across the transition is caused by differences in both the electronic heat capacity and the electron-phonon coupling factor of the AF and FM phases. The electron-phonon coupling constant in the AF phase is also determined using the two-temperature model conducted using the ntmpy code package. For the FM phase, we provide boundaries for the magnitude of the electron-phonon coupling factor for the FM phase. Finally, these results indicate that TDTR can be used to study the transient properties of magnetic materials that are otherwise challenging to probe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Examination of the sensitivity of quasifree reactions to details of the bound-state overlap functions

It is often stated that heavy-ion nucleon knockout reactions are mostly sensitive to the tails of the bound-state wave functions. In contrast, (p,2p) and (p,pn) reactions are known to access information on the full overlap functions within the nucleus. We analyze the oxygen isotopic chain and explore the differences between single-particle wave functions generated with potential models, used in the experimental analysis of knockout reactions, and ab initio computations from self-consistent Green's function theory. Contrary to common belief, we find that not only the tail of the overlap functions, but also their internal part is assessed in both reaction mechanisms, which are crucial to yield accurately determined spectroscopic information.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Light-Induced Renormalization of the Dirac Quasiparticles in the Nodal-Line Semimetal ZrSiSe

In nodal-line semimetals, linearly dispersing states form Dirac loops in the reciprocal space with a high degree of electron-hole symmetry and a reduced density of states near the Fermi level. The result is reduced electronic screening and enhanced correlations between Dirac quasiparticles. In this study we investigate the electronic structure of ZrSiSe, by combining time- and angle-resolved photoelectron spectroscopy with ab initio density functional theory (DFT) complemented by an extended Hubbard model ( DFT + U + V ) and by time-dependent DFT + U + V . We show that electronic correlations are reduced on an ultrashort timescale by optical excitation of high-energy electrons-hole pairs, which transiently screen the Coulomb interaction. Our findings demonstrate an all-optical method for engineering the band structure of a quantum material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Autoionizing Polaritons in Attosecond Atomic Ionization

Light-induced states and Autler-Townes splitting of laser-coupled states are common features in the photoionization spectra of laser-dressed atoms. The entangled light-matter character of metastable Autler-Townes multiplets, which makes them autoionizing polaritons, however, is still largely unexplored. We employ attosecond transient-absorption spectroscopy in argon to study the formation of polariton multiplets between the 3s –1 4p and several light-induced states. We measure a controllable stabilization of the polaritons against ionization, in excellent agreement with ab initio theory. Using an extension of the Jaynes-Cummings model to autoionizing states, we show that this stabilization is due to the destructive interference between the Auger decay and the radiative ionization of the polaritonic components. Furthermore, these results give new insights into the optical control of electronic structure in the continuum and unlock the door to applications of radiative stabilization in metastable polyelectronic systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Ab initio quantum many-body description of superconducting trends in the cuprates

Using a systematic ab initio quantum many-body approach that goes beyond low-energy models, we directly compute the superconducting pairing order and estimate the pairing gap of several doped cuprate materials and structures within a purely electronic picture. We find that we can correctly capture two well-known trends: the pressure effect, where the pairing order and gap increase with intra-layer pressure, and the layer effect, where the pairing order and gap vary with the number of copper-oxygen layers. From these calculations, we observe that the strength of superexchange and the covalency at optimal doping are the best descriptors for these trends. Our microscopic analysis further identifies that strong short-range spin fluctuations and multi-orbital charge fluctuations drive the development of the pairing order. Our work illustrates the possibility of a material-specific ab initio understanding of unconventional high-temperature superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗