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

A Molecular View of Methane Activation on Ni(111) through Enhanced Sampling and Machine Learning

A combination of machine learned interatomic potentials (MLIPs) and enhanced sampling simulations is used to investigate the activation of methane on a Ni(111) surface. The work entails the development and iterative refinement of MLIPs, initially trained on a dataset constructed via ab initio molecular dynamics (AIMD) simulations, supplemented by adaptive biasing forces, to enrich the sampling of catalytically relevant configurations. Our results reveal that by incorporating collective variables that capture the behavior of the reactant molecule, as well as additional frames that describe the dynamic response of the catalytic surface, it is possible to enhance considerably the accuracy of predicted energies and forces. By employing enhanced sampling schemes in the refinement of the MLIP, we systematically explore the potential energy surface, leading to a refined MLIP capable of predicting DFT-level energies and forces and replicating key geometric characteristics of the catalytic system. The resulting free energy landscapes at several temperatures provide a detailed view of the thermodynamics and dynamics of methane activation. Specifically, as methane approaches and dissociates on the catalytic surface, the process involves the dynamic interplay of CH 4 and the Ni catalyst that includes both enthalpic and entropic contributions. The progression towards the transition state involves an CH 4 moiety that is increasingly restrained in its ability to rotate or translate, while the stage following the transition state is characterized by a notable rise of the Ni atom that interacts with the cleaved C–H bond. Furthermore, this leads to an increase in the mobility of the adsorbed species, a feature that becomes more pronounced at higher temperatures.

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

Capturing Coupled Structural and Electronic Motions During Excited-State Intramolecular Proton Transfer via Computational Multiedge Resonant Inelastic X-ray Scattering

Proton-transfer processes build the foundation of many chemical processes. In Excited State Intramolecular Proton Transfer (ESIPT) processes the proton transfer process is impulsively started through light. Here, in this study, we explore the changes in coupled atomic and electronic motions during and following ESIPT through computational time-resolved Resonant Inelastic X-ray Scattering (trRIXS). Excited-state Ab Initio Molecular Dynamics (AIMD) simulations combined with Time-Dependent Density Functional Theory (TDDFT) calculations were performed for 10-hydroxybenzo[h]quinoline to obtain trRIXS signatures. RIXS at both the nitrogen and oxygen K-edges were computed to resolve the dynamics electronic structure from both the proton donor and acceptor perspective. The results show how RIXS can reveal the local electronic structure, the coupling between different electronic states and how electronic structure and coupling change during the proton transfer process. Additionally, we observe a strong correlation between spectral changes and structural changes during ESIPT.

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

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

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

Breaking the Brønsted–Evans–Polanyi Relation with Dual-Metal Sites

Linear scaling relationships impose inherent limitations on catalyst activity; the Brønsted−Evans−Polanyi (BEP) relation, which correlates activation and reaction energies, is a prominent example. Here we report a dual-metal site catalyst (DMSC) on ceria that breaks the BEP relation for C−C coupling of methyl intermediates an elementary step in methane coupling to form ethane. The DMSC structure on CeO 2 (111) was discovered by density-functional theory (DFT) structural exploration and confirmed to be stable via ab initio thermodynamics and ab initio molecular dynamics. Homonuclear and heteronuclear DMSCs of Ni, Pd, Pt, Fe, Ru, Os, Co, Rh, and Ir (45 pairs in total) were examined for methyl affinity and methyl−methyl coupling activation energy. We found that many heteronuclear DMSCs break the BEP linear scaling due to a mixed low-affinity/high-affinity coadsorption of the two methyl groups, decoupling the step responsible for the activation energy (E a ) at the low-affinity site from the overall reaction energy (ΔE) determined by both sites. This mechanism of breaking the BEP relationship via the DMSCs offers a catalyst design principle for C−C coupling reactions.

Chen, Yiming [Vanderbilt Univ., Nashville, TN (Uni↗
DOE OSTI2026

Transition Metal Dichalcogenide MoS 2 : Oxygen and Fluorine Functionalization for Selective Plasma Processing

Low-temperature plasma processing is a promising technique for tailoring transition metal dichalcogenides (TMDs). For chalcogen substitution processing, a key challenge is to identify the ion energy window that enables selective chalcogen removal while preserving the metal lattice. Using ab initio molecular dynamics (AIMD), we demonstrate that oxygen and fluorine functionalization widen the processing window by significantly lowering the sulfur sputtering energy threshold (E sputt,S ) of MoS 2 from ∼30 to ∼10 eV via formation of sputtering products such as SO 2 and SF n . Additionally, we show that experimentally relevant cryogenic temperatures strongly affect E sputt,S (T). The dependence is confirmed via AIMD and also predicted by a mechanistic parameter-free theory, suggesting that E sputt (T) generalizes to other TMDs, functionalizations, and surface impact conditions. Our results highlight oxygen/fluorine functionalization, ionic impact angle, and material temperature to be key control parameters for selective, damage-controlled chalcogen removal in TMD processing.

Polyachenko, Yury [Princeton Plasma Physics Labora↗
DOE OSTI2025

Nanoconfinement of High Hydrogen-to-Metal Ratio Lanthanum Hydrides in Functionalized Carbon Hosts

Metal hydrides with a high hydrogen content are important for materials-based hydrogen storage and high-temperature superconductivity. Nanoconfinement of metal hydrides in porous hosts is a promising strategy to tune the thermodynamic stability and control the hydrogen-to-metal ratio. However, lanthanum hydride (LaH x ) nanoconfinement in porous materials has been limited due to the challenges associated with isolating and stabilizing nanoparticles of La or La-hydrides. Here we successfully demonstrated the chemical reduction of La­(III) salts to La(0), and subsequent infiltration into pure CMK-3 and nitrogen-doped CMK-3 (NCMK-3) porous carbons. Transmission electron microscopy measurements revealed a uniform distribution of LaH x species within the carbon hosts, while X-ray absorption and photoelectron spectroscopy provided detailed information about the local chemical environment. Sieverts measurements indicate that LaH x @NCMK-3 could desorb up to 0.75 wt % hydrogen, which is higher than non-nitrogen-functionalized CMK-3 (0.43 wt % H). Density Functional Theory and ab initio molecular dynamics calculations indicate that host–guest interaction energies are favorable for porous carbon with nitrogen defects, which is supported by experimental evidence. Moreover, high-pressure synchrotron X-ray diffraction measurements were conducted using a diamond anvil cell up to 60 GPa and reveal that the nitrogen-functionalized nanoporous carbon host favors the formation of higher H:La ratios in the presence of ammonia borane compared to the pure CMK-3 host. This approach could serve as a suitable platform for developing nanoscale superconducting materials at lower pressures and temperatures compared to bulk.

Shivanna, Mohana↗
DOE OSTI2020

Theoretical Design of Lithium Chloride Superionic Conductors for All-Solid-State High-Voltage Lithium-Ion Batteries

The development of solid electrolytes (SEs) is a promising pathway to improve the energy density and safety of conventional Li ion batteries. Several lithium chloride SEs, Li 3 MCl 6 (M = Y, Er, In, and Sc), have gained popularity due to their high ionic conductivity, wide electrochemical window, and good chemical stability. This study systematically investigated 17 Li 3 MCl 6 SEs to identify novel and promising lithium chloride SEs. Calculation results revealed that 12 Li 3 MCl 6 (M = Bi, Dy, Er, Ho, In, Lu, Sc, Sm, Tb, TI, Tm, and Y) were stable phase with a wide electrochemical stability window and excellent chemical stability against cathode materials and moisture. Li-ion transport properties were examined using bond valence site energy (BVSE) and ab initio molecular dynamics (AIMD) calculation. Li 3 MCl 6 showed the lower migration energy barrier in monoclinic structures, while orthorhombic and trigonal structures exhibited higher energy barriers due to the sluggish diffusion along the two-dimensional path based on the BVSE model. AIMD results confirmed the slower ion migration along the 2D path, exhibiting lower ionic diffusivity and higher activation energy in orthorhombic and trigonal structures. For the further increase of ionic conductivity in monoclinic structures, Li-ion vacancy was formed by the substitution of M 3+ with Zr 4+ . Zr-substituted phase (Li 2.5 M 0.5 Zr 0.5 Cl 6 , M = In, Sc) exhibited up to a fourfold increase in ionic conductivity. This finding suggested that the optimization of Li vacancy in the Li 3 MCl 6 SEs could lead to superionic Li 3 MCl 6 SEs.

25 ENERGY STORAGE↗
DOE OSTI2022

Nonsacrificial Additive for Tuning the Cathode–Electrolyte Interphase of Lithium-Ion Batteries

Solid–electrolyte interphases is essential for stable cycling of rechargeable batteries. The traditional approach for interphase design follows the decomposition of additives prior to the host electrolyte, which, as governed by the thermodynamic rule, however, inherently limits the viable additives. Here we report an alternative approach of using a nonsacrificial additive. This is exemplified by the localized high-concentration electrolytes, where the fluoroethylene carbonate (FEC) plays a nonsacrificial role for modifying the chemistry, structure, and formation mechanism of the cathode–electrolyte interphase (CEI) layers toward enhanced cycling stability. On the basis of ab initio molecular dynamics simulations, we further reveal that the unexpected activation of the otherwise inert species in the interphase formation is due to the FEC–Li + coordinated environment that altered the electronic states of reactants. In conclusion, the nonsacrificial additive on CEI formation opens up alternative avenues for the interphase design through the use of the commonly overlooked, anodically stable compounds.

25 ENERGY STORAGE↗
DOE OSTI2022

Synthesis of Fluorine-Doped Lithium Argyrodite Solid Electrolytes for Solid-State Lithium Metal Batteries

Solid-state lithium metal batteries (SSLMBs) that utilize novel solid electrolytes (SEs) have garnered much attention because of their potential to yield safe and high-energy-density batteries. Sulfide-based argyrodite-class SEs are an attractive option because of their impressive ionic conductivity. Recent studies have shown that LiF at the interface between Li and SE enhances electrochemical stability. However, the synthesis of F-doped argyrodites has remained challenging because of the high temperatures used in the state-of-the-art solid-state synthesis methods. In this work, for the f irst time, we report F-doped Li 5+y PS 5 F y argyrodites with a tunable doping content and dual dopants (F-/Cl- and F-/Br-) that were synthesized through a solvent-based approach. Among all compositions, Li 6 PS 5 F 0.5 C l0.5 exhibits the highest Li-ion conductivity of 3.5 x 10 -4 S cm -1 at room temperature (RT). Furthermore, Li symmetric cells using Li 6 PS 5 F 0.5 Cl 0.5 show the best cycling performance among the tested cells. X-ray photoelectron spectroscopy and ab initio molecular dynamics simulations revealed that the enhanced interfacial stability of Li 6 PS 5 F 0.5 C l0.5 SE against Li metal can be attributed to the formation of a stable solid electrolyte interphase (SEI)-containing conductive species (Li 3 P), alongside LiCl and LiF. These findings open new opportunities to develop high-performance SSLMBs using a novel class of F-doped argyrodite electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2020

Tribochemical Conversion of Methane to Graphene and Other Carbon Nanostructures: Implications for Friction and Wear

Tribochemistry involves chemical reactions occurring at sliding contact interfaces in the presence of gaseous and/or liquid media. It often leads to the formation of a solid reaction film (also termed boundary film) which controls friction and wear and hence the efficiency and reliability of moving mechanical systems (such as engines). Here we demonstrate tribochemical conversion of methane to graphene, nano-onion, and disordered carbons on the sliding surfaces of Ni-, Cu-, and CuNi-containing VN coatings at atmospheric pressure and room temperature, providing 2-3 orders of magnitude reduction in wear and similar to 50% reduction in friction compared to those of the uncoated steels. Transmission electron microscopy confirms that graphene forms preferably on metal rich nanoclusters of the composite coatings, while the carbon nano-onions are scattered throughout the carbon tribofilm. Ab initio molecular dynamics simulations elucidate underlying mechanisms involved in the tribochemical conversion of methane to carbon- based nanostructures in support of microscopic observations. These scientific findings may lead to new materials technologies that can use methane as a source for continuous and in situ lubrication. For example, there is an urgent need to curtail the uses of lubricating oils in natural gas compressors and engines as they contaminate the natural gas being compressed or burnt.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimizing the Solvent Reorganization Free Energy by Metal Substitution for Nanocage Catalysis

The supramolecular capsule Ga 4 L 6 12– has been found to catalyze a number of important chemical reactions, but the close proximity of a poorly organized external solvent environment can serve as an impediment to the reaction chemistry. Here, using ab initio molecular dynamic calculations and electric field analyses, we find that metal substitution of indium for gallium lowers the total and electrostatic activation barriers by ~3–4 kcal/mol in water solvent. Using an energy decomposition analysis of the interaction of water solvent with the metal vertices, we find that Pauli repulsion between the metal and water decreases, allowing the water coordination with the metal to increase, upon substitution of In for Ga. We therefore determine that the stabilization of the transition state is due to a better arrangement of water molecules around the metal vertices of In 4 L 6 12– and provides a proof of concept of how to redesign the nanocage scaffold in order to reduce the solvent reorganization energy.

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 OSTI2026

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗
DOE OSTI2026

NH 3 -Mediated Reactive Capture and Conversion: Integrating CO 2 Absorption from Flue Gas with CO Production via NH 4 HCO 3 Electrolysis

Efficient carbon capture and utilization require strategies that minimize energy penalties of CO 2 regeneration and compression. Reactive capture and conversion (RCC) address this challenge by integrating capture with direct electrochemical conversion. Here, we show an NH3-mediated tandem RCC system that couples capture of CO 2 from simulated flue gas (10% v/v CO 2 in N 2 ) with electroreduction of NH 4 HCO 3 to CO over a Ni single-atom catalyst (Ni-SAC). Speciation modeling and capture experiments revealed that a deep CO 2 capture with C/N ratio of 0.65 was achieved using 2.5 M NH 3 from simulated flue gas. Electrolysis of the resulting NH 4 HCO 3 on the Ni- SAC delivered an 85% CO Faradaic efficiency at 100 mA/cm 2 with excellent tolerance to NH 3 /NH 4 + as confirmed by DFT calculations and ab initio molecular dynamics (AIMD) simulations. Further, the technoeconomic analysis established a levelized total cost of CO manufacturing of $25.43/kmol, gauging the practical viability. Overall, this study holds great potential to decarbonize the chemical manufacturing industry while reducing synthetic production costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Ab Initio Simulation of Medium-Range Ordering in Ionic Glass Electrolytes LiSiPON and LiNaSiPON

Ionic glasses can exhibit a unique combination of optical transparency, electronic resistivity, ionic conductivity, and mechanical ductility. The origin of the relatively high ionic conductivity and ductility is poorly understood. Recently, these ionic glasses were found to have medium-range ordering (MRO) similar to that of metallic glasses. This MRO significantly impacts the properties of metallic glasses and is also expected to significantly impact the properties of ionic glasses. Here, this work used ab initio molecular dynamics (AIMD) simulations to study the effect of ionic glass composition on the MRO. The AIMD models showed that the degree of MRO increased as the temperature of the LiSiPON ionic glass decreased. The modeling also showed that the MRO is suppressed when 50% of Li is substituted with Na. This work lays the foundation for using AIMD to identify clear structure–property relationships in ionic glasses, allowing their properties to be optimized.

Osetsky, Yuri N. [Oak Ridge National Laboratory (O↗
DOE OSTI2021

Correlating Li-Ion Solvation Structures and Electrode Potential Temperature Coefficients

Temperature coefficients (TCs) for either electrochemical cell voltages or potentials of individual electrodes have been widely utilized to study the thermal safety and cathode/anode phase changes of lithium (Li)-ion batteries. However, the fundamental significance of single electrode potential TCs is little known. Here, we discover that the Li-ion desolvation process during Li deposition/intercalation is accompanied by considerable entropy change, which significantly contributes to the measured Li/Li + electrode potential TCs. To explore this phenomenon, we compare the Li/Li + electrode potential TCs in a series of electrolyte formulations, where the interaction between Li-ion and solvent molecules occurs at varying strength as a function of both solvent and anion species as well as salt concentrations. As a result, we establish correlations between electrode potential TCs and Li-ion solvation structures and further verify them by ab initio molecular dynamics simulations. We show that measurements of Li/Li + electrode potential TCs provide valuable knowledge regarding the Li-ion solvation environments and could serve as a screening tool when designing future electrolytes for Li-ion/Li metal batteries.

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

Probing the Local Site Disorder and Distortion in Pyrochlore High-Entropy Oxides

High-entropy oxides (HEOs) have attracted great interest in diverse fields because of their inherent opportunities to tailor and combine materials functionalities. The control of local order/disorder in the class is by extension a grand challenge toward realizing their vast potential. In this paper, we report the first examples of pyrochlore HEOs with five M-site cations, for Nd 2 M 2 O 7 , in which the local structure has been investigated by neutron diffraction and pair distribution function (PDF) analysis. The average structure of the pyrochlores is found to be orthorhombic Imma, in agreement with radius-ratio rules governing the structural archetype. The computed PDFs from density functional theory relaxed special quasirandom structure models are compared with real space PDFs in this work to evaluate M-site order/disorder. Reverse Monte Carlo combined with ab initio molecular dynamics and Metropolis Monte Carlo simulations demonstrates that Nd 2 (Ta 0.2 Sc 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7 is synthesized with its M-site local to nanoscale order highly randomized/disordered, while Nd 2 (Ti 0.2 Nb 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7+x exhibits a strong distortion of the TiO 6 octahedron and small degree of Ti chemical short-range order (SRO) on the subnanometer scale. Calculations suggest that this may be intrinsic, energetically favored SRO rather than due to sample processing. These results offer an important demonstration that the engineered variation of participating ions in HEOs, even among those with very similar radii, provides richly diverse opportunities to control local order/disorder motifs—and therefore materials properties for future designs. This work also hints at the exquisite level of detail that may be needed in computational and experimental data analysis to guide structure–property tuning in the emerging HEO materials class.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗