Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initio molecular dynamics”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First-Principles Treatment of Vibrational Broadening in X-Ray Excited Valence Band Spectra for n-SrTiO 3 (001)

The valence band maximum (VBM) is an important quantity for semiconductors as it locates the Fermi level relative to the band edge. Accurate measurement of this quantity in near-surface regions of semiconductors by photoemission is a first step toward determining the electronic properties of heterostructures involving these materials. While extrapolating the leading edge of the valence band to the energy axis in photoemission spectra is a widely used way to find the VBM, this method can be ambiguous if the leading edge exhibits multiple slopes. Another way to determine the VBM is to fit the leading edge to an appropriately broadened, cross-section modulated theoretical density of states (DOS). Three kinds of broadening that should be included for maximum accuracy are those due to: (1) finite instrumental resolution, (2) valence hole lifetime, and (3) vibrational excitations. While steps (1) and (2) are straightforward to implement, (3) is more difficult because the appropriate amount of broadening is not known a priori. Here, we demonstrate that explicit inclusion of vibrational broadening using ab initio molecular dynamics facilitates accurate VBM determination for n-SrTiO 3 (001). The total DOS is constructed by summing time-averaged projections at elevated temperature onto s - , p - , and d orbitals for the constituent atoms and modulating with the associated photoemission cross sections. Subsequent convolutions of the total DOS, first with a Gaussian of width equal to the experimental energy resolution and second with a Lorentzian to simulate valence hole lifetime effects, yield line shapes that reproduce the experimental leading edges rather well. The VBM is then given by the energy at which the vibrationally broadened total DOS (prior to the convolutions) goes to zero. The VBMs generated by this method quantitatively agree with those resulting from extrapolating from the middle of the measured leading edge for SrTiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stabilization of Hydroxide Ions at the Interface of a Hydrophobic Monolayer on Water via Reduced Proton Transfer

We report a joint study using surface-specific sum-frequency vibrational spectroscopy and ab initio molecular dynamics simulations, respectively, on a pristine hydrophobic (sub)monolayer hexane-water interface, namely, the hexane/water interface with varied vapor pressures of hexane and different pHs in water. We show clear evidence that hexane on water revises the interfacial water structure in a way that stabilizes the hypercoordinated solvation structure and slows down the migration of hydroxide ion (OH – ) relative to that in bulk water. Furthermore, this mechanism effectively attracts the OH – to the water-hydrophobic interface with respect to its counterion. The result illustrates the striking difference of proton transfer of hydrated OH – at the interface and in the bulk, which is responsible for the intrinsic charging effect at the hydrophobic interface.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Understanding the lattice thermal conductivity of SrTiO 3 from an ab initio perspective

Here, we present a detailed analysis of the structure dependence of the lattice thermal conductivity of SrTiO 3 . We have used both ab initio molecular dynamic simulations and density functional theory calculations to decouple the effect of different structural distortions on the thermal conductivity. We have identified two main mechanisms for tuning the thermal conductivity when a distortion is applied. First, the modification of the acoustic-modes energy dispersion when a change in the lattice parameters is imposed and second, the low energy polar modes. In particular and counterintuitively, we have found that an increase in the angle of the oxygen octahedral rotations increases the thermal conductivity due to its coupling to these polar modes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Probing ferroelectricity by x-ray absorption spectroscopy in molecular crystals

We carry out x-ray absorption spectroscopy experiment at the oxygen $K$ edge in croconic acid (C 5 H 2 O 5 ) crystal as a prototype of ferroelectric organic molecular solid, whose electric polarization is generated by proton transfer. The experimental spectrum is well reproduced by the electron-hole excitation theory simulations from configuration generated by ab initio molecular dynamics simulation. When inversion symmetry is broken in the ferroelectric state, the hydrogen bonding environment on the two bonded molecules become inequivalent. Such a difference is sensitively probed by the bound excitation in the pre-edge, which is strongly localized on the excited molecules. Our analysis shows that a satellite peak in the pre-edge will emerge at higher excitation energy, which serves as a clear signature of ferroelectricity in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio free energies of liquid metal alloys: Application to the phase diagrams of Li-Na and K-Na

Comparison of free energies between different phases and different compositions underlies the prediction of alloy phase diagrams. To allow direct comparison, consistent reference points for the energies or enthalpies are required, and the entropy must be placed on an absolute scale, yielding absolute free energies. Here we derive absolute free energies of liquids from ab-initio molecular dynamics (AIMD) by combining the directly simulated enthalpies with an entropy derived from simulated densities and pair correlation functions. Additionally, as an example of the power of this method we calculate the phase diagrams of two binary alkali metal alloys, Li-Na and K-Na, revealing a critical point and liquid-liquid phase separation in the former case, and a deep eutectic in the latter. Good agreement with experimental data demonstrates the power of this simple method.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Alkali metal halide–coated perovskite redox catalysts for anaerobic oxidative dehydrogenation of n -butane

Oxidative dehydrogenation (ODH) of n-butane has the potential to efficiently produce butadiene without equilibrium limitation or coke formation. Despite extensive research efforts, single-pass butadiene yields are limited to <23% in conventional catalytic ODH with gaseous O 2 . This article reports molten LiBr as an effective promoter to modify a redox-active perovskite oxide, i.e., La 0.8 Sr 0.2 FeO 3 (LSF), for chemical looping–oxidative dehydrogenation of n-butane (CL-ODHB). Under the working state, the redox catalyst is composed of a molten LiBr layer covering the solid LSF substrate. Characterizations and ab initio molecular dynamics (AIMD) simulations indicate that peroxide species formed on LSF react with molten LiBr to form active atomic Br, which act as reaction intermediates for C–H bond activation. Meanwhile, molten LiBr layer inhibits unselective CO 2 formation, leading to 42.5% butadiene yield. The redox catalyst design strategy can be extended to CL-ODH of other light alkanes such as iso-butane conversion to iso-butylene, providing a generalized approach for olefin production.

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

Altered morphology and diffusivity of water confined in MXenes: Machine learning–accelerated computations combined with experiments

Nanoconfined water exhibits unique properties compared to bulk water due to limited quantities, frustrated hydrogen bonding, and surface interactions, which are fundamental for energy storage and transport applications. We integrate machine learning–accelerated ab initio molecular dynamics with x-ray diffraction (XRD) and inelastic neutron scattering (INS) to systematically analyze the thermodynamic and dynamic behavior of water confined between functionalized (-F, -O, and -OH) two-dimensional (2D) Ti 3 C 2 T x MXene layers. As water intercalates between layers, the interlayer spacing exhibits layer-dependent staging characteristics. The water polarization can be flipped by the count and morphology of intercalated molecules interacting with MXene surface groups, resulting in varying electrostatic potential profiles. On the basis of interfacial electrostatic potential, hydrogen bond lifetime, and molecular orientation, we establish a linear combination of exponential model describing water diffusivity. These computational insights align well with experimental x-ray and neutron measurements, suggesting strategies for tuning water morphology and transport by tailoring MXene surface chemistry and water content for electrochemical energy storage and nanofluidic applications.

Tang, Jiawei [Southeast Univ., Nanjing (China)] (O↗
DOE OSTI2023

Many-Body Effects in Aqueous Systems: Synergies Between Interaction Analysis Techniques and Force Field Development

Interaction analysis techniques, including the many-body expansion (MBE), symmetry-adapted perturbation theory, and energy decomposition analysis, allow for an intuitive understanding of complex molecular interactions. We review these methods by first providing a historical context for the study of many-body interactions and discussing how nonadditivities emerge from Hamiltonians containing strictly pairwise-additive interactions. We then elaborate on the synergy between these interaction analysis techniques and the development of advanced force fields aimed at accurately reproducing the Born–Oppenheimer potential energy surface. In particular, we focus on ab initio–based force fields that aim to explicitly reproduce many-body terms and are fitted to high-level electronic structure results. These force fields generally incorporate many-body effects through ( a) parameterization of distributed multipoles, ( b) explicit fitting of the MBE, ( c) inclusion of many-atom features in a neural network, and ( d) coarse-graining of many-body terms into an effective two-body term. We also discuss the emerging use of the MBE to improve the accuracy and speed of ab initio molecular dynamics.

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

Highly Concentrated Electrolytes: Electrochemical and Physicochemical Characteristics of LiPF 6 in Propylene Carbonate Solutions

Highly concentrated electrolytes (HCEs) based on LiPF 6 in propylene carbonate (PC) have been examined as lithium-ion battery electrolytes. These HCEs have lower ionic conductivities and higher viscosities than ethylene carbonate (EC) electrolytes with 1.2 M LiPF 6 , but they have higher Li + ion transference numbers. Electrochemical cycling behaviour of LiNi 0.8 Co 0.015 Al 0.05 O 2 //graphite cells with 3.2 M LiPF 6 in PC resembles that of cells with EC-based electrolytes; the HCE cells have higher impedance, which can be lowered by increasing test temperature. By employing Raman and infrared spectroscopy, combined with density functional theory and ab initio molecular dynamics simulations, we reveal that the Li + solvation structure and speciation are key factors that determine cell performance. Two distinct regimes are observed as a function of salt concentration—in the conventional regime, the solvation number (SN) is mostly constant, while in the HCE regime it decreases linearly. Graphite exfoliation is suppressed only at very high salt concentrations ( > 2.4 M), where [ PC ] free /[Li + ] < 1 and P F 6 − f r e e > P C f r e e . Results from the Advanced Electrolyte Model indicate that Li + desolvation improves at higher LiPF 6 concentrations, thereby mitigating PC co-intercalation into the graphite. However, Li + ion transport is hindered in the HCEs, which increases impedance at both the oxide-positive and graphite-negative electrodes.

25 ENERGY STORAGE↗
DOE OSTI2022

First Principles Determination of the Potential-of-Zero-Charge in an Alumina-Coated Aluminum/Water Interface Model for Corrosion Applications

The surfaces of most metals immersed in aqueous electrolytes have a several-nanometer-thick oxide/hydroxide surface layer. This gives rise to the existence of both metal|oxide and oxide|liquid electrotlyte interfaces, and makes it challenging to correlate atomic length-scale structures with electrochemical properties such the potential-of-zero-charge (PZC). The PZC has been shown to be correlated the onset potential for pitting corrosion. Here, we conduct large-scale Density Functional Theory and ab initio molecular dynamics to calculate the PZC of a Al(111)| γ -Al 2 O 3 (110)| water double-interface model within the context of aluminum corrosion. By partitioning the multiple interfaces involved into binary components with additive contributions to the overall work function and voltage, we predict the PZC to be -1.53 V vs SHE for this model. Furthermore, we calculate the orbital energy levels of defects like oxygen vacancies in the oxide, which are critical parameters in theories associated with pitting corrosion. We predict that the Fermi level at the PZC lies above the impurity defect levels of the oxygen vacancies, which are therefore uncharged at the PZC. From the PZC estimate, we predict the voltage needed to create oxygen vacancies with net positive charges within a flatband approximation.

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

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗
DOE OSTI2025

Mechanistic Insights into Oxygen Reduction Reaction on Metal Cathode Over BaZrO 3 -Based Electrolyte under Moist Atmosphere: Interfacial Interactions and Reaction Pathways

The oxygen reduction reaction (ORR) is a critical process that often governs the overall performance of electrochemical systems such as proton conducting solid oxide fuel cells. Despite its significance, the current understanding about the ORR, especially when it involves H 2 O is still limited. In this study, the proton involved ORR (PI-ORR) on metal cathode (e.g., Pt) over BaZrO 3 (BZO) electrolyte is investigated using density functional theory (DFT) under moist atmosphere. Two scenarios are considered: one at relatively high temperature (900 K) in the presence of oxygen vacancy in BZO substrate and one at relatively low temperature (700 K) without oxygen vacancy. We systematically explored oxygen vacancy effects on adsorption sites and proton transfer pathways, calculated energy profiles and Gibbs free energies for elementary steps along the O 2 dissociated and the O 2 concerted reaction pathways. We also evaluated H 2 O impacts on intermediate oxygen species (*O and *OH) through Ab Initio molecular dynamics. These analyses revealed pathway-dependent reaction mechanisms and hydration effects at triple-phase boundary, advancing fundamental understanding of PI-ORR in metal cathode/Ba-based electrolyte systems for energy applications.

08 - HYDROGEN↗
DOE OSTI2024

The ab initio amorphous materials database: Empowering machine learning to decode diffusivity

Amorphous materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven ex- ploration and design of amorphous materials is hampered by the absence of a com- prehensive database covering a broad chemical space. In this work, we present the largest computed amorphous materials database to date, generated from sys- tematic and accurate ab initio molecular dynamics (AIMD) calculations. We also show how the database can be used in simple machine-learning models to connect properties to composition and structure, here specifically targeting ionic conductiv- ity. These models predict the Li-ion diffusivity with speed and accuracy, offering a cost-effective alternative to expensive density functional theory (DFT) calculations. Furthermore, the process of computational quenching amorphous materials provides a unique sampling of out-of-equilibrium structures, energies, and force landscape, and we anticipate that the corresponding trajectories will inform future work in uni- versal machine learning potentials, impacting design beyond that of non-crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

0BGRaman: Graph Network based Simulator for Forecasting Molecular Polarizability

This report presents the work performed under the GRaman project, sponsored by the PCSD LDRD Seed program. The project aimed at accelerating ab initio molecular dynamics simulation using Graph Networks. The Graph Network framework is a ML framework that has been successfully employed to simulate the dynamics of several physical systems: including water splashing in a container and flags moving with the wind. In this effort, we performed a data collection campaign for 3 different molecules of interest. We have built tools for preprocessing the trajectories obtained by simulating Raman Spectroscopy with NWChem and translating them into a suitable format for training. We have developed a training algorithm to train the Graph Network based simulators based on our data and developed a simulator that produces trajectories in the same NWChem format. While the tool has improved with each iteration of development and subsequent experiments, the current state of the tool does not allow to directly incorporate the technology within the NWChem framework because the trajectories produced by the tool are not yet accurate enough. However, the technology has proved to have good potential and it is certainly worth further research and development.

97 MATHEMATICS AND COMPUTING↗