Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive force field”

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

Highly Altered State of Proton Transport in Acid Pools in Charged Reverse Micelles

Transport mechanisms of solvated protons of 1 M HCl acid pools, confined within reverse micelles (RMs) containing the negatively charged surfactant sodium bis(2-ethylhexyl) sulfosuccinate (NaAOT) or the positively charged cetyltrimethylammonium bromide (CTABr), are analyzed with reactive force field simulations to interpret dynamical signatures from TeraHertz absorption and dielectric relaxation spectroscopy. We find that the forward proton hopping events for NaAOT are further suppressed compared to a nonionic RM, while the Grotthuss mechanism ceases altogether for CTABr. We attribute the sluggish proton dynamics for both charged RMs as due to headgroup and counterion charges that expel hydronium and chloride ions from the interface and into the bulk interior, thereby increasing the pH of the acid pools relative to the nonionic RM. For charged NaAOT and CTABr RMs, the localization of hydronium near a counterion or conjugate base reduces the Eigen and Zundel configurations that enable forward hopping. Thus, localized oscillatory hopping dominates, an effect that is most extreme for CTABr in which the proton residence time increases dramatically such that even oscillatory hopping is slow.

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

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2020

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗
DOE OSTI2025

Atomistic characterization of hydration-dependent fuel cell ionomer nanostructure: validation by vibrational spectroscopy

The development of Nafion alternatives for fuel cells and electrolyzers requires a fundamental understanding of hydration-dependent ion-exchange site acid/base chemistry. We present here reactive force field (ReaxFF) molecular dynamics (MD) simulations of Nafion at varying molar water/ion-exchange-site ratios (λ), which we correlate to our experimental and density functional theory-based vibrational spectra. ReaxFF describes the formation and breaking of covalent bonds, enabling simulations of proton exchange between sulfonic acid/sulfonate groups and water/hydronium. Our MD simulations determine the λ-dependent equilibrium proportions of protonated and deprotonated sites. We find that protonated sites persist across all λ from 0 to 20, challenging the widely accepted notion that all sites are ionized above a threshold λ value (e.g., 3 or 4). Our simulations generate hundreds of realistic exchange site environments, the characterization of which are based on 6 Å radii sulfur-centered ‘inner-spheres’. These inner- and outer-sphere regions elucidate an interplay of stereoelectronic factors that influence protonation states, including the number of inner-sphere waters (Λ). Our simulations produce broad λ-dependent distributions of Λ values, representing non-uniform exchange site hydration. Moreover, we demonstrate that these distributions are specific to both protonated and deprotonated sites, with significant overlap between the two distributions for all non-zero λ. These distributions underpin the IR spectra of hydrated membranes, with each exchange site contributing an IR spectrum characteristic of its protonation state and Λ. We expect that these nanostructural characterizations of Nafion exchange sites will contribute to the development of new ionomers.

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

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2022

Thermal conductivity tensor of γ and ε -hexanitrohexaazaisowurtzitane as a function of pressure and temperature

Using reverse non-equilibrium molecular dynamics simulations, we have determined the dependences on temperature and pressure of the thermal conductivity tensors for the monoclinic γ and ε polymorphs of hexanitrohexaazaisowurtzitane (HNIW or CL20). In this work, a recently developed non-reactive force field [X. Bidault and S. Chaudhuri, RSC Adv. 9, 39649–39661 (2019)], designed to study polymorphism and phase transitions in CL20, is employed. The effects of temperature and pressure are investigated between 200 and 500 K and up to 0.5 GPa for γ-CL20 and 2 GPa for ε-CL20. In order to obtain the full thermal conductivity tensor, $κ_{ij}$, for the monoclinic crystals, four distinct heat propagation directions are used. We find that $κ_{ij}$ for both polymorphs is more isotropic than for other energetic molecular crystals, including α- and γ-RDX, β-HMX, and PETN, with a maximum difference of 9.8% between orientations observed at 300 K and 0 GPa for γ-CL20 and a maximum difference of 4.8% for ε-CL20. The average thermal conductivity, $\bar{κ}$, of ε-CL20 is 6.4% larger than that of γ-CL20 at 300 K and 0 GPa. Analytic linear functions of the inverse temperature and the pressure are provided, which fit the data well and can be used to predict the thermal conductivity of both polymorphs for any orientation, pressure, and temperature in and around the fitting range. Our predictions agree reasonably well with the limited available experimental data, for which the polymorph type is unknown.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determining the hydronium pKα at platinum surfaces and the effect on pH-dependent hydrogen evolution reaction kinetics

Electrocatalytic hydrogen evolution reaction (HER) is critical for green hydrogen generation and exhibits distinct pH-dependent kinetics that have been elusive to understand. A molecular-level understanding of the electrochemical interfaces is essential for developing more efficient electrochemical processes. Here we exploit an exclusively surface-specific electrical transport spectroscopy (ETS) approach to probe the Pt-surface water protonation status and experimentally determine the surface hydronium pK a = 4.3. Quantum mechanics (QM) and reactive dynamics using a reactive force field (ReaxFF) molecular dynamics (RMD) calculations confirm the enrichment of hydroniums (H 3 O + * ) near Pt surface and predict a surface hydronium pK a of 2.5 to 4.4, corroborating the experimental results. Importantly, the observed Pt-surface hydronium pK a correlates well with the pH-dependent HER kinetics, with the protonated surface state at lower pH favoring fast Tafel kinetics with a Tafel slope of 30 mV per decade and the deprotonated surface state at higher pH following Volmer-step limited kinetics with a much higher Tafel slope of 120 mV per decade, offering a robust and precise interpretation of the pH-dependent HER kinetics. These insights may help design improved electrocatalysts for renewable energy conversion.

Zhong, Guangyan↗
DOE OSTI2025

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry↗
DOE OSTI2024

An examination of the performance of molecular dynamics force fields: Silicon and silicon dioxide reactive ion etching

Classical molecular dynamics (MD) simulations of plasma-surface interactions were performed of physical sputtering and reactive ion etching (RIE), with predictions based on several force fields. In this paper, we focus mainly on SiO 2 but include some results for Si substrates as well. Here, we compare predictions from these MD simulations to experimental studies of SiO 2 physical sputtering (by Ar + ions), RIE of Si, and RIE of SiO 2 (both using F atoms and Ar + ions). MD results using different published force fields are compared to reported yields from published vacuum beam experiments. The near-surface depth profiles predicted using different force fields are compared. One motivation for the present study is to document the nature and magnitude of differences in the predictions for selected systems and conditions of practical interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On the Network Topology of Cross-Linked Acrylate Photopolymers: A Molecular Dynamics Case Study

A reactive molecular dynamics approach is used to simulate cross-linking of acrylate polymer networks. By employing the same force field and reactive scheme and studying three representative multifunctional acrylate monomers, we isolate the importance of the nonreactive moieties within these model monomers. Analyses of reactive trajectories benchmark the estimated gel points, cyclomatic character, and spatially resolved cross-linking tendencies of the acrylates as a function of conversion. The insights into the similarities and differences of the polymerization and resulting networks presented here suggest molecular mechanics as a useful tool in the rational design of photopolymerization resins.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations

In this investigation, force field-based molecular dynamics (MD) simulations have been employed to generate detailed structural representations for a range of amorphous quaternary CaO-MgO-Al{sub 2}O{sub 3}-SiO{sub 2} (CMAS) and ternary CaO-Al{sub 2}O{sub 3}-SiO{sub 2} (CAS) glasses. Comparison of the simulation results with select experimental X-ray and neutron total scattering and literature data reveals that the MD-generated structures have captured the key structural features of these CMAS and CAS glasses. Based on the MD-generated structural representations, we have developed two structural descriptors, specifically (i) average metal oxide dissociation energy (AMODE) and (ii) average self-diffusion coefficient (ASDC) of all the atoms at melting. Both structural descriptors are seen to more accurately predict the relative glass reactivity than the commonly used degree of depolymerization parameter, especially for the eight synthetic CAS glasses that span a wide compositional range. Hence these descriptors hold great promise for predicting CMAS and CAS glass reactivity in alkaline environments from compositional information.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ReaxFF-based molecular dynamics study of bio-derived polycyclic alkanes as potential alternative jet fuels

This work investigates the initial stages of the pyrolysis of HtH-1 (C 18 H 32 ; 2,2,7,7,8a,8b-hexamethyl-dodecahydrobiphenylene) and HtH-2 (C 18 H 34 ; 1,1',3,3,3',3'-hexamethyl-1,1'-bi(cyclohexane)), which are bio-derived polycyclic alkanes and potential jet fuels, using ReaxFF force field based molecular dynamics (MD) simulations. Global Arrhenius parameters, such as activation energies and pre-exponential factors, are calculated and used to analyze the overall decomposition kinetics of the fuels. HtH-1 decomposes faster than HtH-2 at the same temperature and density conditions, and they have a faster decomposition rate compared to some existing jet-fuels, such as JP-10. A systematic reaction analysis framework developed in this work is applied to determine a temperature-dependent decomposition mechanism. At lower temperature, the central C-C bond connecting the two cyclohexane rings is dominantly broken in both HtH-1 and HtH-2. However, C-CH 3 bond breaking becomes dominant with increasing temperature due to the large increase in entropy during this reaction. Major products from HtH-1 are C 5 H 8 and C 4 H 8 , and those from HtH-2 are C 4 H 8 and C 2 H 4 . The major products predict that HtH-1 has a higher sooting tendency than HtH-2, which is consistent with measurements. The impact of HtH-2 on the pyrolysis of HtH-1 is also investigated in their binary mixtures. HtH-1 and HtH-2 decompose by unimolecular reactions, and they rarely interact with each other during the pyrolysis of the mixtures. Furthermore, this work demonstrates that ReaxFF can be used to investigate pyrolysis and combustion chemistry of existing or future fuels and to contribute to the development of their chemical kinetic models without any a priori input and chemical intuition.

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

Role of Capping Agents in the Synthesis of Salicylate-Capped Zinc Oxide Nanoparticles

Capping agents are often used for controlling the size, aggregation, and properties of nanoparticles. To guide the design of improved nanomaterials for targeted performance, one can use mechanistic insights into the interactions between capping agents and nanoparticles. Here, we employ density functional theory (DFT), reactive force-field molecular dynamics (ReaxFF MD) simulations, and optical spectroscopy to study the interactions between salicylate, as a model capping agent, and zinc oxide (ZnO) nanoparticles. In this work, we find that salicylate strongly interacts with the nanoparticle via the formation of O–Zn bonds in a distorted six-membered coordination ring structure. We describe the mechanisms of capping of ZnO nanoparticles by salicylate via three different binding modes. Simulations indicate that salicylate undergoes dissociative adsorption at the highly active surface Zn sites via a hydrogen-transfer process, thereby forming a tridentate configuration. The water-mediated interaction also facilitates the dissociative adsorption, leading to two salicylate O atoms coordinating with a surface Zn atom, while the other salicylate O atom bonds with another surface Zn atom. For molecular adsorption, binding free energies indicate that salicylate binds more strongly to ZnO, often in a bidentate configuration, than water does. The formation of the salicylate–ZnO complex is substantiated by UV–visible and Fourier transform infrared spectra. We find that the C=O stretching mode of salicylate becomes softened when it interacts with the nanoparticle, suggesting chemisorption of salicylate on ZnO. Although DFT predicts strong interaction between salicylate and ZnO, ReaxFF MD simulation indicates the moderate interaction between these two components in aqueous solution. Water molecules in close contact with the nanoparticle surface undergo dissociation, thus resulting in a surface hydroxyl, a reactive oxygen species that may influence the nanoparticle’s catalytic properties. Overall, the atomic-level information provided here can guide the selection process of salicylate as an appropriate agent for ZnO nanoparticle synthesis when strong interactions of particles with capping agents are required.

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

A Graph Dynamical neural network approach for decoding dynamical states in ferroelectrics.

Ferroelectric materials such as BaTiO 3 show tremendous potential for emerging advances in memory devices, particular neuromorphic type devices. High density of memory can be obtained by stabilising polar domain walls at the nanoscale, regions of discontinuity between the well-defined polarization order parameter, but little is known about what controls their structure and dynamics in real nanoscale materials. Indeed, chiral polar domain walls have been observed in heterogeneous ferroelectrics, such as oxygen-deficient BaTiO 3 , but very little is known about how such polar-domains walls interact with defects. Indeed, a critical understanding of how dynamics of domain-walls depend on point-defects is crucial to create engineered ferroelectric memory devices. For this work, we perform large-scale simulations of nansocale domain-wall dynamics in pristine and defective BaTiO 3 using reactive force-field developed by us earlier (Phys. Chem. Chem. Phys., 2019, 21, 18240–18249), and capture their dynamical dependence on point defects using a graph dynamical neural-network approach, which we adapted to interrogate solids with well-defined order-parameters, and implemented using Pytorch based libraries. Our machine learning (ML) approach goes beyond the traditional post-processing methods to capture both spatial and temporal heterogeneities of large-scale molecular dynamics simulations of complex defective ferroelectric oxide materials. We crucially find that isolated oxygen vacancies introduce very localized spatial regions (~1–2 unit-cell in length) that show slow dipole relaxation due to formation of defect-dipoles, and that these defect-dipoles in turn slow the intrinsic dynamics of domain walls. Further, the roughness of domain walls, also influenced by vacancies, introduce dynamic heterogeneity along the domain-wall. As such we find a novel mechanism by which quenched disorder due to defects introduce dynamic heterogeneity thereby influencing response to external fields (particularly time varying fields) in a ferroelectric. Our study also emphasizes the need for creating digital twins of dynamical quantities to achieve autonomous in operando control of nanoscale switching.

42 ENGINEERING↗
DOE OSTI2022

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

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

Single-Step Direct Laser Writing of Multimetal Oxygen Evolution Catalysts from Liquid Precursors

We investigate a laser direct-write method to synthesize and deposit metastable, mixed transition metal oxides and evaluate their performance as oxygen evolution reaction catalysts. This laser processing method enabled the rapid synthesis of diverse heterogeneous alloy and oxide catalysts directly from cost-effective solution precursors, including catalysts with a high density of nanocrystalline metal alloy inclusions within an amorphous oxide matrix. The nanoscale heterogenous structures of the synthesized catalysts were consistent with reactive force-field Monte Carlo calculations. By evaluating the impact of varying transition metal oxide composition ratios, we created a stable Fe 0.63 Co 0.19 Ni 0.18 O x /C catalyst with a Tafel slope of 38.23 mV dec -1 and overpotential of 247 mV, a performance similar to that of IrO 2 . Synthesized Fe 0.63 Co 0.19 Ni 0.18 O x /C and Fe 0.14 Co 0.46 Ni 0.40 O x /C catalysts were experimentally compared in terms of catalytic performance and structural characteristics to determine that higher iron content and a less crystalline structure in the secondary matrix decreases the charge transfer resistance and thus is beneficial for electrocatalytic activity. Furthermore, this conclusion is supported by density-functional theory calculations showing distorted active sites in ternary metal catalysts are key for lowering overpotentials for the oxygen evolution reaction.

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

Effects of interlayer confinement and hydration on capacitive charge storage in birnessite

Nanostructured birnessite exhibits high specific capacitance and nearly ideal capacitive behaviour in aqueous electrolytes, rendering it an important electrode material for low-cost, high-power energy storage devices. The mechanism of electrochemical capacitance in birnessite has been described as both Faradaic (involving redox) and non-Faradaic (involving only electrostatic interactions). To clarify the capacitive mechanism, we characterized birnessite’s response to applied potential using ex situ X-ray diffraction, electrochemical quartz crystal microbalance, in situ Raman spectroscopy and operando atomic force microscope dilatometry to provide a holistic understanding of its structural, gravimetric and mechanical responses. These observations are supported by atomic-scale simulations using density functional theory for the cation-intercalated structure of birnessite, ReaxFF reactive force field-based molecular dynamics and ReaxFF-based grand canonical Monte Carlo simulations on the dynamics at the birnessite–water–electrolyte interface. We show that capacitive charge storage in birnessite is governed by interlayer cation intercalation. We conclude that the intercalation appears capacitive due to the presence of nanoconfined interlayer structural water, which mediates the interaction between the intercalated cation and the birnessite host and leads to minimal structural changes.

36 MATERIALS SCIENCE↗