Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio theory”

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 OSTI2020

Ultrafast photoinduced band splitting and carrier dynamics in chiral tellurium nanosheets

Trigonal tellurium (Te) is a chiral semiconductor that lacks both mirror and inversion symmetries, resulting in complex band structures with Weyl crossings and unique spin textures. Detailed time-resolved polarized reflectance spectroscopy is used to investigate its band structure and carrier dynamics. The polarized transient spectra reveal optical transitions between the uppermost spin-split H4 and H5 and the degenerate H6 valence bands (VB) and the lowest degenerate H6 conduction band (CB) as well as a higher energy transition at the Lpoint. Surprisingly, the degeneracy of the H6 CB (a proposed Weyl node) is lifted and the spin-split VB gap is reduced upon photoexcitation before relaxing to equilibrium as the carriers decay. Using ab initio density functional theory (DFT) calculations, we conclude that the dynamic band structure is caused by a photoinduced shear strain in the Te film that breaks the screw symmetry of the crystal. The band-edge anisotropy is also reflected in the hot carrier decay rate, which is a factor of two slower along the c-axis than perpendicular to it. The majority of photoexcited carriers near the band-edge are seen to recombine within 30 ps while higher lying transitions observed near 1.2 eV appear to have substantially longer lifetimes, potentially due to contributions of intervalley processes in the recombination rate. These new findings shed light on the strong correlation between photoinduced carriers and electronic structure in anisotropic crystals, which opens a potential pathway for designing novel Te-based devices that take advantage of the topological structures as well as strong spin-related properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

One dimensional wormhole corrosion in metals

Corrosion is a ubiquitous failure mode of materials. Often, the progression of localized corrosion is accompanied by the evolution of porosity in materials previously reported to be either three-dimensional or two-dimensional. However, using new tools and analysis techniques, we have realized that a more localized form of corrosion, which we call 1D wormhole corrosion, has previously been miscategorized in some situations. Using electron tomography, we show multiple examples of this 1D and percolating morphology. To understand the origin of this mechanism in a Ni-Cr alloy corroded by molten salt, we combined energy-filtered four-dimensional scanning transmission electron microscopy and ab initio density functional theory calculations to develop a vacancy mapping method with nanometer-resolution, identifying a remarkably high vacancy concentration in the diffusion-induced grain boundary migration zone, up to 100 times the equilibrium value at the melting point. Deciphering the origins of 1D corrosion is an important step towards designing structural materials with enhanced corrosion resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pd-Ru pair on Pt surface for promoting hydrogen oxidation and evolution in alkaline media

Hydrogen oxidation reaction in alkaline media is critical for alkaline fuel cells and electrochemical ammonia compressors. The slow hydrogen oxidation reaction in alkaline electrolytes requires large amounts of scarce and expensive platinum catalysts. While transition metal decoration can enhance Pt catalysts’ activity, it often reduces the electrochemical active surface area, limiting the improvement in Pt mass activity. Here, we enhance Pt catalysts’ activity without losing surface-active sites by using a Pd-Ru pair. Utilizing a mildly catalytic thermal pyrolysis approach, Pd-Ru pairs are decorated on Pt, confirmed by extended X-ray absorption fine structure and high-angle annular dark-field scanning transmission electron microscopy. Density functional theory and ab-initio molecular dynamics simulations indicate preferred Pd and Ru dopant adsorption. The Pd-Ru decorated Pt catalyst exhibits a mass-based exchange current density of 1557 ± 85 A g –1 metal for hydrogen oxidation reaction, demonstrating superior performance in an ammonia compressor.

25 ENERGY STORAGE↗
DOE OSTI2021

Efficient diffusion of superdense lithium via atomic channels for dendrite-free lithium–metal batteries

The non-uniform aggregation of fast-diffused Li on an anode surface would aggravate its tip-effect-induced nucleation/growth, leading to the notorious dendrite growth in Li metal batteries (LMBs). Tuning the Li diffusion on the anode surface has been regarded previously as a mainstream method to induce its uniform deposition, while the diffusion of Li in the anode bulk is usually ignored. In this study, conceptually different from the classic surface modification, we propose a molecular tunnelling strategy to construct atomic channels in graphite bulk, which enables the fast diffusion of superdense Li. Density functional theory calculations and ab initio molecular dynamics simulations prove that the bulk diffusion through atomic channels could become a new and dominating path. Its reversible and efficient diffusion has been further visualized by in situ transmission electron microscopy. As a result, when coupled with high-loading LiFePO 4 cathodes (20 mg cm –2 ), a high areal capacity and 100% capacity retention are achieved over 370 cycles. Through this work a new strategy is developed based on the bulk-diffusion of superdense Li for dendrite-free LMBs, which can be pervasive in other high-performance energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2024

Exploring Cr and molten salt interfacial interactions for molten salt applications

Molten salts play an important role in various energy-related applications such as high-temperature heat transfer fluids and reaction media. However, the extreme molten salt environment causes the degradation of materials, raising safety and sustainability challenges. A fundamental understanding of material–molten salt interfacial evolution is needed. This work studies the transformation of metallic Cr in molten 50/50 mol% KCl–MgCl 2 via multi-modal in situ synchrotron X-ray nano-tomography, diffraction and spectroscopy combined with density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. Notably, in addition to the dissolution of Cr in the molten salt to form porous structures, a δ-A15 Cr phase was found to gradually form as a result of the metal–salt interaction. This phase change of Cr is associated with a change in the coordination environment of Cr at the interface. DFT and AIMD simulations provide a basis for understanding the enhanced stability of δ-A15 Cr vs. bcc Cr, by revealing their competitive phase thermodynamics at elevated temperatures and probing the interfacial behavior of the molten salt at relevant facets. This study provides critical insights into the morphological and chemical evolution of metal–molten salt interfaces. Finally, the combination of multimodal synchrotron analysis and atomic simulation also offers an opportunity to explore a broader range of systems critical to energy applications.

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

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
DOE OSTI2025

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electronic and magnetic properties of quasi-one-dimensional osmium halide OsCl 4

Using ab initio density functional theory, we study the electronic and magnetic properties of the van der Waals chain material OsCl 4 . In the nonmagnetic state, a strongly anisotropic band structure was observed, in agreement with its anticipated one-dimensional crystal geometry. Based on Wannier functions, we found that the four electrons of the 5d Os atom form a low-spin S = 1 state, with a large crystal field between the d xz/yz and d xy orbitals, corresponding to a strong Jahn–Teller distortion (Q 3 < 0). As a consequence, the magnetic properties are mainly contributed by the d xz/yz states. Furthermore, when a Mott gap develops after the introduction of the Hubbard U and Hund coupling J, we found that the staggered spin order is the most likely magnetic state, namely, spins arranged as (↑-↓-↑-↓) with π wavevector along the chain. In addition, the energy differences between various spin states are small, suggesting a weak magnetic exchange coupling along the chain. Overall, our results provide guidance to experimentalists and theorists working on quasi-one-dimensional osmium halides chain materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reaction processes at step edges on S-decorated Cu(111) and Ag(111) surfaces: MD analysis utilizing machine learning derived potentials

We report a variety of complexation, reconstruction, and sulfide formation processes can occur at step edges on the {111} surfaces of coinage metals (M) in the presence of adsorbed S under ultra-high vacuum conditions. Given the cooperative many-atom nature of these reaction processes, Molecular Dynamics (MD) simulation of the associated dynamics is instructive. However, only quite restricted Density Functional Theory (DFT)-level ab initio MD is viable. Thus, for M = Ag and Cu, we instead utilize the DeePMD framework to develop machine-learning derived potentials, retaining near-DFT accuracy for the M–S systems, which should have broad applicability. These potentials are validated by comparison with DFT predictions for various key quantities related to the energetics of S on M(111) surfaces. The potentials are then utilized to perform extensive MD simulations elucidating the above diverse restructuring and reaction processes at step edges. Key observations from MD simulations include the formation of small metal–sulfur complexes, especially MS 2 ; development of a local reconstruction at A-steps featuring an S-decorated {100} motif; and 3D sulfide formation. Additional analysis yields further information on the kinetics for metal–sulfur complex formation, where these complexes can strongly enhance surface mass transport, and on the propensity for sulfide formation.

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

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Chemistry↗
DOE OSTI2023

Ab initio methods for polariton chemistry

Polariton chemistry exploits the strong interaction between quantized excitations in molecules and quantized photon states in optical cavities to affect chemical reactivity. Molecular polaritons have been experimentally realized by the coupling of electronic, vibrational, and rovibrational transitions to photon modes, which has spurred a tremendous theoretical effort to model and explain how polariton formation can influence chemistry. This tutorial review focuses on computational approaches for the electronic strong coupling problem through the combination of familiar techniques from ab initio electronic structure theory and cavity quantum electrodynamics, toward the goal of supplying predictive theories for polariton chemistry. Our aim is to emphasize the relevant theoretical details with enough clarity for newcomers to the field to follow, and to present simple and practical code examples to catalyze further development work.

Chemistry↗
DOE OSTI2023

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application (“Fragme∩t”), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Chemistry↗
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 OSTI2024

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. Further, the potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: $T$ m = $1969.23+19.15P-0.012P$ 2 . In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi-task graph neural networks for simultaneous prediction of global and atomic properties in ferromagnetic systems *

Abstract We introduce a multi-tasking graph convolutional neural network, HydraGNN, to simultaneously predict both global and atomic physical properties and demonstrate with ferromagnetic materials. We train HydraGNN on an open-source ab initio density functional theory (DFT) dataset for iron-platinum with a fixed body centered tetragonal lattice structure and fixed volume to simultaneously predict the mixing enthalpy (a global feature of the system), the atomic charge transfer, and the atomic magnetic moment across configurations that span the entire compositional range. By taking advantage of underlying physical correlations between material properties, multi-task learning (MTL) with HydraGNN provides effective training even with modest amounts of data. Moreover, this is achieved with just one architecture instead of three, as required by single-task learning (STL). The first convolutional layers of the HydraGNN architecture are shared by all learning tasks and extract features common to all material properties. The following layers discriminate the features of the different properties, the results of which are fed to the separate heads of the final layer to produce predictions. Numerical results show that HydraGNN effectively captures the relation between the configurational entropy and the material properties over the entire compositional range. Overall, the accuracy of simultaneous MTL predictions is comparable to the accuracy of the STL predictions. In addition, the computational cost of training HydraGNN for MTL is much lower than the original DFT calculations and also lower than training separate STL models for each property.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First-principles study of the low-temperature charge density wave phase in the quasi-one-dimensional Weyl chiral compound ( TaSe 4 ) 2 I

Using ab initio density functional theory, we study herein the lattice phase transition of quasi-one-dimensional (TaSe 4 ) 2 I. In the undistorted state, the strongly anisotropic semimetal band structure presents two nonequivalent Weyl points. In previous efforts, two possible Ta-tetramerization patterns were proposed to be associated with the low-temperature structure. Our phonon calculations indicate that the orthorhombic F222 CDW-I phase is the most likely ground state for this quasi-one-dimensional system. In addition, the monoclinic C2 CDW-II phase may also be stable according to the phonon dispersion spectrum. Since these two phases have very similar energies in our density functional theory calculations, both these Ta-tetramerization distortions likely compete or coexist at low temperatures. The semimetal-to-insulator transition is induced by a Fermi-surface-driven instability that supports the Peierls scenario, which affects the Weyl physics developed above T CDW . Furthermore, the spin-orbit coupling generates Rashba-like band splittings in the insulating charge density wave phases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reduction of electron-phonon coupling in warm dense iron

The electron-ion relaxation dynamics in warm dense iron were investigated by time-resolved x-ray absorption near-edge spectroscopy (XANES). A novel analysis, combining ab initio density functional theory (DFT) and two temperature model (TTM) simulations, was developed to calculate the x-ray absorption spectra as a function of delay time. Here we present experimental evidence of changes at the XAS L3 edge of iron that are consistent with the reduction of the electron-phonon coupling under warm dense matter conditions. The experimental results are in agreement with the model that takes both the electron ($T_e$) and the ion temperature ($T_i$) dependence of the thermophysical properties into consideration, while models where either constant electron-phonon coupling factor (G) or only $T_e$- dependent G are used do not agree with the observed relaxation dynamics of iron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗