Search NASA⌕ Search

SEARCH · Search NASA

Results for “ReaxFF 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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↗
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 OSTI2025

Diffusion of Atoms in Glassy Mixtures of Deuterium and Lithium

The diffusion coefficients of D and Li were calculated in amorphous, glassy Li+D mixtures for various concentrations of D in lithium over the temperature range 100–1000 K. The densities of the mixtures as a function of temperature were also determined. The diffusion coefficients were obtained by the analysis of the mean-squared displacement using molecular dynamics with ReaxFF, a reactive force field. In conclusion, the diffusion and density data obtained for Li+D mixtures were compared with the available experimental and calculated data reported in the literature.

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

Computational Advances in Ionic Liquid Applications for Green Chemistry: A Critical Review of Lignin Processing and Machine Learning Approaches

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

09 BIOMASS FUELS↗
DOE OSTI2020

Optimization of the Reax force field for the lithium–oxygen system using a high fidelity charge model

Studies using molecular dynamics (MD) have long struggled to simulate the failure modes of materials, predicting unrealistically high ductility and failing to capture brittle fracture. The primary cause of this shortcoming is an inadequate description of bond breaking. While reactive force fields such as ReaxFF show improvements compared to traditional force fields, the charge models used yield unphysical partial charges, especially during dissociation of ionic bonds. This flaw may be remedied by using the atom-condensed Kohn–Sham density functional theory (DFT) approximated to a second order (ACKS2) charge model for determining partial charges. In this work, we present a new ACKS2-enabled Reax force field for fracture simulations of lithium oxide systems, which was obtained by training against an extensive set of DFT, multireference configuration interaction (MRCI), and MRCI+Q reference data using genetic optimization techniques. This new force field significantly improves the bond breaking behavior, but still cannot fully capture the brittle fracture in MD simulations, suggesting more research is needed to improve simulation of brittle fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanistic study of pH effect on organic solvent nanofiltration using carboxylated covalent organic framework as a modeling and experimental platform

Covalent organic framework (COF) membranes for organic solvent nanofiltration (OSN) have attracted growing interest due to their high stability, well-defined pore size, and controllable pore functionality. In aqueous nanofiltration processes, solution pH has a significant effect on membrane permeability, selectivity, and fouling. However, our understanding of pH effects on OSN processes is still limited. In this work, we present experimental and reactive force field molecular dynamics (ReaxFF MD) simulation results on the effect of pH on OSN using a well-characterized carboxylated COF (C-COF) membrane. Filtration experimental and modeling results both demonstrate that pH changes can affect solvent permeance and solute rejection dramatically. Specifically, methanol permeance decreases with the increase of NaOH and HCl concentration, respectively. When Alcian Blue dye was used as the solute in methanol, the selectivity increased from 23% to 98% when changing the feed pH from 2.2 to 10.1, respectively. While ReaxFF MD results show that C-COF solvated pore size is only marginally affected by the pH of solution, methanol diffusion is affected by the size of methanol-solvated Na + and Cl - ions, which increase with increasing Na + and Cl - concentrations. Solution pH change also affected solute solvation layer thicknesses and C-COF pore charges, contributing to the significant AB dye selectivity increase with pH increase. To summarize, we have shown that dramatic changes in C-COF selectivity are not anomalous in OSN, but are instead controllable by varying feed solution pH to modify solvated solute radii and surface charge.

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

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
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 OSTI2020

Thickness and strain dependence of piezoelectric coefficient in BaTiO 3 thin films

We explore the thickness dependence of the converse piezoelectric coefficient (d 33 ) in epitaxial thin films of BaTiO 3 (BTO) grown on (001) SrTiO 3 substrates. Piezoresponse force microscope was performed using an atomic force microscope equipped with an interferometric displacement sensor allowing direct quantification of electromechanical coupling coefficients in BTO free from unwanted background contributions. We find that 80-nm-thick films exhibit a d 33 of ~20.5pm/V, but as the thickness is reduced, the d 33 reduces to less than 2 pm/V for a 10 nm film. To explain the atomistic origin of the effect, we performed molecular dynamics simulations with a recently developed ab initio-derived reactive force field, constructed using the ReaxFF framework. Simulations predict that under applied electric fields thin films of BaTiO 3 show an increasing thickness, with compressive strain, of the region screening the depolarization-field. This study confirms quantitatively the drop in piezoelectric performance in BTO ultrathin films and again highlights the importance of the screening mechanisms when films approach the ultrathin limits in dictating the functional behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D

Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. Here, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

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

Reactive force fields for aqueous and interfacial magnesium carbonate formation

Here we develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral–water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg 2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite–water and brucite–water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg–CO 3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral–water–CO 2 interface in water-poor conditions.

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

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2023

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger–Weber potentials. In this study, the model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. In conclusion, the model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

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

Carbon fiber synthesis from pitch: Insights from ReaxFF based molecular dynamics simulations

Despite the potential to lowering fabrication costs of pitch-based carbon fiber (CF) compared to polyacrylonitrile-based CF, large scale utilization is currently limited by the complexity in manufacturing. This is mainly due to the lack of consistency between the pitch structure/chemistry and CF properties under different processing conditions. Here, we employed full atomistic simulations based on reactive force fields to investigate the conversion of pitch precursors to CFs. Additionally, it was found that at the carbonization stage, the consumption rates of precursors, gas formations, and solid yields (carbonized CFs) are all dependent on the composition and arrangement of molecular moieties in the precursor. Specifically, precursors with sp 3 carbons in the cata-condensed core area have limited decomposition pathways (i.e., limited reaction sites), which leads to the formation of carbonized CFs with similar compositions. Furthermore, carbonized CFs possess distinct structural characteristics, and finally lead to the formation of CFs with different moduli. By analyzing CF crystallinity and elemental composition, general rationales are proposed for the selection of optimal precursors and processing conditions, and thus shed lights on parameter screening during industrial CF synthesis. For instance, inclusion of sp 3 carbons in cata- condensed rings can be preferred over the presence of sp 3 carbons on side chains.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Effect of Domain Size, Boundary, and Loading Conditions on Mechanical Properties of Amorphous Silica: A Reactive Molecular Dynamics Study

Mechanical properties are very important when choosing a material for a specific application. They help to determine the range of usefulness of a material, establish the service life, and classify and identify materials. The size effect on mechanical properties has been well established numerically and experimentally. However, the role of the size effect combined with boundary and loading conditions on mechanical properties remains unknown. In this paper, by using molecular dynamics (MD) simulations with the state-of-the-art ReaxFF force field, we study mechanical properties of amorphous silica (e.g., Young’s modulus, Poisson’s ratio) as a function of domain size, full-/semi-periodic boundary condition, and tensile/compressive loading. We found that the domain-size effect on Young’s modulus and Poisson’s ratio is much more significant in semi-periodic domains compared to full-periodic domains. The results, for the first time, revealed the bimodular and anisotropic nature of amorphous silica at the atomic level. We also defined a “safe zone” regarding the domain size, where the bulk properties of amorphous silica can be reproducible, while the computational cost and accuracy are in balance.

bimodular materials↗
DOE OSTI2020

Simulations of the Biodegradation of Citrate Based Polymers for Artificial Scaffolds Using Accelerated Reactive Molecular Dynamics

Here, we explore the reactivity and mechanical properties of poly(1,6-hexanediol-co-Citric acid) via ReaxFF molecular dynamics simulations. We implement an accelerated scheme within the ReaxFF framework to study the hydrolysis reaction of the polymer which is provided with sufficient amount of energy known as the restrain energy after a suitable pre-transition state configuration is obtained to overcome the activation energy barrier and the desired product is obtained. The validity of the ReaxFF force field is established by comparing the ReaxFF energy barriers of ester and ether hydrolysis with benchmark DFT values in literature. We perform chemical and mechanical degradation of polymer chain bundles at 300 K. We find that ester hydrolyzes faster than ether due to lower activation energy barrier of the reaction. The selectivity of the bond-boost scheme has been demonstrated by lowering the boost parameters of the accelerated simulation which almost stops the ether hydrolysis. Mechanical degradation of pre-hydrolyzed and intermittent hydrolyzed polymer bundles is performed along the longitudinal direction at two different strain rates. We find that the tensile modulus of the polymers are increases with increase in strain rates which shows that polymers show a strain dependent behavior. The tensile modulus of the Polyester-ether is higher than Polyester yet reaches yield-stress faster than Polyester. This makes Polyester more ductile than Polyester-ether.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗