Reactive Force Research Articles

Development of a ReaxFF reactive force field for ternary phosphate-based bioactive glasses.

Phosphate-based glasses (PBGs) in the CaO-Na2O-P2O5 system have diverse applications as biomaterials due to their unique dissolution properties. A reactive force field (ReaxFF) has been developed to simulate these materials at the atomic level. The ReaxFF parameters of PBGs, including the interaction between phosphorus and calcium atoms, have been developed using a published code based on genetic algorithms. The training data, including the atomic charges, atomic forces, bond and angle parameters, and different differential energies, are chosen and measured from static quantum-mechanical calculations and abinitio molecular dynamics of PBGs. We did different short- and medium-range structural analyses on the bulk simulated PBGs with different compositions to validate the developed potential. Radial and angular distribution functions and coordination numbers of network formers and modifiers, as well as the network connectivity of the glass, are in agreement with experimental and previous simulations using both shell-model classical force fields and abinitio simulated data; for example, the coordination number of phosphorus is 4.0. This successful development of ReaxFF parameters being able to describe the bulk PBGs enables us to work on the dissolution behavior of the glasses, including the interaction of water molecules with PBGs, in future works.

Combustion characteristics and nitrogen conversion mechanism in ammonia/coal Co-firing process

The coupled combustion of coal with the carbon-free fuel NH3 is an effective technology for reducing CO2 emissions. Based on elemental analysis, proximate analysis, Fourier Transform infrared spectroscopy (FTIR), Carbon-13 Nuclear Magnetic Resonance (13C NMR), and X-ray photoelectron spectroscopy (XPS) characterizations, a coal macromolecular model was constructed. The combustion characteristics of coal/NH3 co-firing were investigated by reactive force field molecular dynamics simulations (ReaxFF MD). The transformation mechanism of elemental N was elucidated. The influence of temperature, oxygen concentration, and NH3 co-firing ratio (CR) on products was investigated. The results showed that during the non-isothermal combustion phase, the coal macromolecules decomposed more rapidly when the heating rate was low. In the isothermal phase, high temperatures promoted secondary reactions. The production of inorganic gases increased with increasing temperature. As the molar fraction of NH3 increased, the production of CO2 and CO was effectively reduced. There were three main pathways for NO formation. As NH3 decomposed to form NH2, NH2 was subsequently converted to NH, and NH reacted with O2 to form NO. The NH2 could also react with O radicals to form HNO, which was then converted to NO. Besides, both thermal NOx and prompt NOx were generated since the N2 was involved in the reaction. In addition, the production of NOx increased significantly at oxygen concentrations greater than 1.4.

Effect of H2O on macroscopic flame behaviors and combustion reaction mechanism of 1,1-difluoroethane (R152a)

1,1-difluoroethane (R152a) is one of the most prospective low GWP refrigerants but may trigger fires and explosions once it leaks. Understanding the flame propagation process and combustion mechanisms at different relative humidity conditions is essential to use flammable refrigerants safely. In this study, a high-speed camera has recorded the flame propagation of R152a combustion near the flammability limit. Combining with macroscopic experimental phenomena, we revealed the effect of H2O on the microscopic mechanism of R152a combustion by using reactive force field molecular dynamics simulations with a reliable force field. The promotion effects of H2O on the oxidation of R152a were revealed from an atomistic perspective. The differences of oxidation intermediates and products in different environments were analyzed for the first time. The H2O/O2 environment was the most potent promoter of R152a decomposition, followed by the O2 and the pyrolysis environment. The O2/H2O environment reduced the apparent activation energy of R152a, significantly enhanced the consumption rate of R152a, and enriched the number of species. O2 reacts with H2O or H to form OH to accelerate the reaction process. The H2O could provide more OH active fragments for the reaction. Moreover, it promotes the formation of HF and H2, H2O + F→HF + OH, H2O + H→H2 + HO. This study aims to provide a guiding theory for the safe application and disaster prevention of flammable refrigerants.

Multi-methodological analysis of hydrogen desorption from graphene

Chemisorption of hydrogen on graphene has been extensively investigated in last decades, but the reported values of over several eV, depending on different factors, such as the local structural environment, substrates or previously bound hydrogen.In this work, we combine Temperature-Programmed Desorption experiments with simulations using empirical reactive force fields and Density Functional Theory calculations. We observe that desorption occurs through processes with different activation energies, which we are able to identify by analyzing simulation data and post-processing experimental ones through a procedure based on Arrhenius fit with feedback. The result is the assignment of experimentally observed desorption peaks to different desorption processes, namely desorption from isolated sp3 H, or couples of H in conformation cis or trans, or edge sp3 or sp2 H.The protocol we outline is generally applicable and might be particularly useful in more complex cases with larger multiplicity of desorption peaks.

ЧИСЕЛЬНЕ МОДЕЛЮВАННЯ ПРОЦЕСУ ВПЛИВУ ПРОДУКТІВ ДЕТОНАЦІЇ НА НАДЗВУКОВИЙ ПОТІК ПОВІТРЯ

Currently, there is a large number of control bodies used in modern rocket technol-ogy. Gasdynamic systems are usually used in rocket engines to change the thrust vector, less often - to create asymmetric air flow around the case. But the application of the det-onation process for such scheme can have better energy and dynamic characteristics compared to existing systems. The driving force in such scheme is created not only by the reactive mass force of detonation products ejected from the gas generator, but also by the effect of an intense shock wave on the pattern of supersonic flow around the rocket. Therefore, it is necessary to evaluate the performance and efficiency of the experiment. The purpose of the study is a priori numerical study of the detonation wave effect on the supersonic flow in the nozzle. Modeling was carried out in the Solid Works application software package. The geometric parameters of the model coincide with the experimental model. It is a flat nozzle, one of the walls is longer than the other and imitates the surface of the rocket. Air, accelerating in the nozzle to supersonic speeds with Mach number M=1.3...2, flows around this wall. Setting devices affecting the flow in it, it is possible to investigate their influence and the possibility of creating a lateral (controlling) force. When detonation is used to disturb the flow, an intense shock wave is generated in it, which changes the pres-sure distribution on the wall surface. Some holes for pressure sensors are made in the wall to fix this disturbance over time. A priori numerical modeling was carried out in order to estimate the parameters of the interaction of detonation products with the supersonic air flow. The pressure range and duration of the process have been determined, which allows selection of equipment and planning of the experiment. Patterns of velocity and pressure distribution in the model over time were obtained.

Study on the evolutionary mechanism of the cavity inside aluminum nanoparticles while reacting with oxygen

Aluminum nanoparticles are often considered attractive choices for fuels in energetic materials. In this study, molecular dynamics simulation with the reactive force field (ReaxFF) was used to study the oxidation process of three different particle sizes (4 nm, 5 nm, 6 nm) of ANPs (aluminum nano-particles). The cavity formation mechanism was obtained during the ANPs oxidation process. The cavity size, the oxide layer thickness, and the active aluminum content of ANPs with time were obtained. The results indicated that the formation of cavities is highly dependent on particle size. Large particle size prevents the formation of cavities. In 0–100 ps, the average growth rate of the oxide layer is 0.33, 0.25, and 0.20 nm/ps, respectively. Within 100–250 ps, the average growth rate of the oxide layer is −0.017, 0.038, and 0.117 nm/ps, respectively. Our results emphasize the effect of particle size on cavity formation, which provides insights into the atomic-scale oxidation mechanism of metal nanoparticles.

Ab initio informed machine learning potential for tribochemistry and mechanochemistry: Application for eco–friendly gallate lubricant additive

Mechanochemistry and tribochemistry processes involve multiple physical/chemical interactions induced by extreme conditions including molecular confinement, high temperatures and mechanical stress applied. Simulating these processes by molecular dynamics is very challenging. While force fields fall short reproducing the enhanced reactivity arising by quantum effects, ab initio molecular dynamics is severely limited by the complexity of the systems of interest, their sizes, and the long-time scale on which relevant events take place. In this work using an active learning approach, a landmark deep neural network potential has been developed which reproduces the accuracy of ab initio interactions at the classical molecular dynamics computational cost and permits to successfully simulate the tribochemical processes occurring at the interface between butyl gallate molecules and iron substrates under tribological conditions. The simulations of the dynamics of the Fe-gallate system when sliding under an imposed external load reveal the key atomistic mechanisms underlying the formation of the friction reducing lubricant tribofilm and permit to characterize the tribological properties of the explored systems, clearly exposing the shortcoming of reactive force field based approaches. The successful development of neural network potentials, making it possible to push the limits of molecular dynamics marrying accuracy with system sizes and long time scales, paves the route toward a new area in computational tribochemistry.

Reactive Force Field Molecular Dynamics Investigation of NH3 Generation Mechanism during Protein Pyrolysis Process.

The mechanism of ammonia formation during the pyrolysis of proteins in biomass is currently unclear. To further investigate this issue, this study employed the AMS 2023.104 software to select proteins (actual proteins) as the model compounds and the amino acids contained within them (assembled amino acids) as the comparative models. ReaxFF molecular dynamics simulations were conducted to explore the nitrogen transformation and NH3 generation mechanisms in three-phase products (char, tar, and gas) during protein pyrolysis. The research results revealed several key findings. Regardless of whether the model compounds are actual proteins or assembled amino acids, NH3 is the primary nitrogen-containing product during pyrolysis. However, as the temperature rises to higher levels, such as 2000 K and 2500 K, the amount of NH3 decreases significantly in the later stages of pyrolysis, indicating that it is being converted into other nitrogen-bearing species, such as HCN and N2. Simultaneously, we also observed significant differences between the pyrolysis processes of actual proteins and assembled amino acids. Notably, at 2000 K, the amount of NH3 generated from the pyrolysis of assembled amino acids was twice that of actual proteins. This discrepancy mainly stems from the inherent structural differences between proteins and amino acids. In proteins, nitrogen is predominantly present in a network-like structure (NH-N), which shields it from direct external exposure, thus requiring more energy for nitrogen to participate in pyrolysis reactions, making it more difficult for NH3 to form. Conversely, assembled amino acids can release NH3 through a simpler deamination process, leading to a significant increase in NH3 production during their pyrolysis.

Resilient Si@C architecture fabricated by mechanochemical conjugation and thermal reduction for lithium storage: Simulation on interfacial state evolution

Synchronal comminution via high-energy bead milling not only efficiently pulverizes Si particles and exfoliates graphene oxide (GO) sheets but also facilitates Si surface functionalization and intensify Si-C conjugation. After a thermal reduction, a well-wrapped Si@C anode material can be fabricated for lithium storage. The discharge specific capacities of this Si@C composite sustain at 1110.3 mAh g−1 after 300 cycles at 1 A g−1 and 316.4 mAh g−1 for 800 cycles at 5 A g−1, which are respectively 1.35 and 2.07 times that of its random blended Si/C counterpart. By employing molecular dynamics (MD) simulations with ReaxFF reactive force field to investigate the mechanochemical reactions between Si and GO amidst high-energy milling, an atomic scale portrayal of Si crystal fragmentation and Si-C bonding can be depicted. Analyses focusing on the differential forces exerted on various Si crystal planes, coupled with potential energy (PE) trends, reveal a higher propensity for cleavage in (111) plane, with (311) plane followed. Radial distribution function (RDF) outcomes indicate that ball milling leads to substantial disruption of crystalline order in Si particles, generating a plethora of Si dangling bonds. This disruption stimulates the formation of SiOOC, SiOC and SiC bonds at the Si/C interface, thus enhancing the integration of GO with Si substrate. Subsequent thermal reduction processes not only restore the structural order of Si and GO but also consolidate Si/C bonding.

Revealing the nitrogen migration mechanism during pyrolysis and steam gasification of biomass: A combined ReaxFF MD and DFT study

The biomass pyrolysis/gasification technology takes advantages of high efficiency and carbon neutrality, which is promising for high value-added products fabrication. However, biomass possesses relatively high amount of N especially for seaweeds and municipal sludge wastes, and the lack of N migration mechanism has resulted in emission of N-containing pollutants and limited its industrial application. In this work, a simulation investigation combining reactive force field molecular dynamic and density functional theory calculations was performed to study the conversion mechanisms of N-containing species in the pyrolysis/gasification of biomass processes. The glutamic acid and cellulose were selected as the model compounds of N-containing and typical compositions of biomass, respectively, and six models containing multiple mass ratios of glutamic acid to cellulose and steam to feedstock were constructed. The results show that the majority of nitrogen in the original model is converted towards NH3-N (62%) during glutamic acid pyrolysis by direct cleavage, hydrogen capture, deamination, and decarboxylation. The presence of cellulose effectively inhibits the decomposition of NH3 at high temperature range, thereby increasing the yield of NH3-N (up to 70%) and restraining the conversion of nitrogen into mesomolecules. The introduction of H2O provides a substantial amount of •H radicals, which facilitates the conversion of nitrogen into NH3-N. During DFT calculations, eight reaction pathways for the formation of N-containing products were proposed on the basis of distribution of pyrolysis products. The glutamic acid cleavage reaction achieves the lowest energy barrier, and the deamination of glutamic acid towards a ternary cyclic ether is found to be the most feasible pathway for the NH3 production. This work can provide theoretical support for nitrogen pollutants mitigation during pyrolysis and steam gasification of biomass processes.

Photoresponsive Azobenzene Nanocluster-Modified Liposomes: Mechanism Analysis Combining Experiments and Simulations.

Stimuli-responsive behaviors and controlled release in liposomes are pivotal in nanomedicine. To this end, we present an approach using a photoresponsive azobenzene nanocluster (AzDmpNC), prepared from azobenzene compounds through melting and aggregation. When integrated with liposomes, they form photoresponsive vesicles. The morphology and association with liposomes were investigated by using transmission electron microscopy. Liposomes loaded with calcein exhibited a 9.58% increased release after UV exposure. To gain insights into the underlying processes and elucidate the mechanisms involved. The molecular dynamic simulations based on the reactive force field and all-atom force field were employed to analyze the aggregation of isomers into nanoclusters and their impacts on phospholipid membranes, respectively. The results indicate that the nanoclusters primarily aggregate through π-π and T-stacking forces. The force density inside the cis-isomer of AzDmpNC formed after photoisomerization is lower, leading to its easier dispersion, rapid diffusion, and penetration into the membrane, disrupting the densification.

Development of a Zr-Nb-H-O reactive force field for molecular dynamics simulations of in-reactor corrosion

The corrosion of Zircaloy fuel cladding is one of the crucial issues during the operation of pressurized water reactors (PWR). In this work, a Zr-Nb-H-O reactive force field (ReaxFF) is constructed to study the corrosion mechanisms and irradiation effect of Zr-Nb alloys at the atomic landscape. This ReaxFF can describe the interactions between water dissociation products and Zr-Nb alloys, and the stability as well as diffusion properties of irradiation defects in Zr bulk. The molecular dynamics (MD) simulations based on the present ReaxFF show that Nb can thicken the suboxide layer during the corrosion process, as well as the promotion effect of irradiation on corrosion varies under different PKA (Primary Knock-on Atom) energies. The ReaxFF described here has the potential to be a new tool for understanding the in-reactor corrosion mechanism under the irradiation environment.

The Catalytic Effect of Pt on Lignin Pyrolysis: A Reactive Molecular Dynamics Study

Lignin is the second-largest renewable resource in nature, second only to cellulose. Lignin is one of the most significant components of biomass, and it determines the behaviour of biomass in many thermochemical processes. However, limited studies have focused on the influence of metal catalysts on lignin pyrolysis. This study aims to develop a sustainable lignin catalytic pyrolysis technology to improve biomass energy-conversion efficiency, reduce dependence on fossil fuels, and promote the development of clean energy. In this study, the impact of Pt catalyst on the pyrolysis process of hardwood lignin was simulated by using reactive force field (ReaxFF) molecular dynamics. Through the comparison of the system without catalysts, the catalyst exhibited evident attraction to lignin macromolecules, prompting their decomposition at lower temperatures. Additionally, the catalyst has the strongest adsorption capacity for H radical. The activation energy of the reaction was calculated by kinetic analysis. It was found that the addition of catalysts significantly reduced the activation energy of the reaction. By revealing the effect of Pt catalyst on the lignin pyrolysis process, it provides a theoretical basis for biomass pyrolysis and the utilization of metal catalysts in industry.

Inhibition mechanism of CHF3 on hydrogen–oxygen combustion: Insights from reactive force field molecular dynamics simulations

To mitigate the risks linked to hydrogen and oxygen (H2–O2) combustion through CHF3 additives, Reactive Force Field Molecular Dynamics (ReaxFF MD) simulations are performed in this study. The primary objective is to investigate the inhibition mechanism of CHF3 on H2–O2 combustion from 2000 K to 2800 K. The simulation results demonstrate that the reaction pathways of hydrogen combustion are changed under the extended second explosion limit, and the main radicals involved in elementary reactions transform from H, OH, and O to H, OH, and HO2. CHF3 predominantly engages in reactions with H radicals to impede the continuation of chain reactions by forming stable HF molecules. OH radicals react with modest amounts of secondary fluorides such as CHF2OH, CH2F, and so on, while HO2 radicals are combined with even fewer intermediates like CHF2O and CHFOH to prevent chain propagations. Moreover, comprehensive and novel reaction pathways are proposed for the inhibition of H2–O2 combustion by CHF3. The parameters of the reaction kinetics indicate that the ignition delay is advanced and the activation energy for the combustion process is increased under the influence of CHF3. This study is expected to provide practical guidance on the inhibition of H2–O2 combustion by CHF3.

Direct Synthesis of Layer‐Tunable and Transfer‐Free Graphene on Device‐Compatible Substrates Using Ion Implantation Toward Versatile Applications

Direct synthesis of layer‐tunable and transfer‐free graphene on technologically important substrates is highly valued for various electronics and device applications. State of the art in the field is currently a two‐step process: a high‐quality graphene layer synthesis on metal substrate through chemical vapor deposition (CVD) followed by delicate layer transfer onto device‐relevant substrates. Here, we report a novel synthesis approach combining ion implantation for a precise graphene layer control and dual‐metal smart Janus substrate for a diffusion‐limiting graphene formation to directly synthesize large area, high quality, and layer‐tunable graphene films on arbitrary substrates without the post‐synthesis layer transfer process. Carbon (C) ion implantation was performed on Cu–Ni film deposited on a variety of device‐relevant substrates. A well‐controlled number of layers of graphene, primarily monolayer and bilayer, is precisely controlled by the equivalent fluence of the implanted C‐atoms (1 monolayer ~4 × 1015 C‐atoms/cm2). Upon thermal annealing to promote Cu‐Ni alloying, the pre‐implanted C‐atoms in the Ni layer are pushed toward the Ni/substrate interface by the top Cu layer due to the poor C‐solubility in Cu. As a result, the expelled C‐atoms precipitate into a graphene structure at the interface facilitated by the Cu‐like alloy catalysis. After removing the alloyed Cu‐like surface layer, the layer‐tunable graphene on the desired substrate is directly realized. The layer‐selectivity, high quality, and uniformity of the graphene films are not only confirmed with detailed characterizations using a suite of surface analysis techniques but more importantly are successfully demonstrated by the excellent properties and performance of several devices directly fabricated from these graphene films. Molecular dynamics (MD) simulations using the reactive force field (ReaxFF) were performed to elucidate the graphene formation mechanisms in this novel synthesis approach. With the wide use of ion implantation technology in the microelectronics industry, this novel graphene synthesis approach with precise layer‐tunability and transfer‐free processing has the promise to advance efficient graphene‐device manufacturing and expedite their versatile applications in many fields.

Development of the ReaxFF Reactive Force Field for Li/Mn/O Battery Technology with Application to Design a Self-Healing Cathode Electrolyte Interphase

Development of the ReaxFF Reactive Force Field for Li/Mn/O Battery Technology with Application to Design a Self-Healing Cathode Electrolyte Interphase

Investigation of lubrication mechanism of phosphate ionic liquid by ReaxFF molecular dynamics simulations

The frictional characteristics of the interface between ionic liquid (IL) and metal can be significantly improved by the formation of a tribo-film, while the lubrication mechanism remains unclear. This study utilized reactive force field molecular dynamics (ReaxFF MD) simulation to investigate the tribological properties, interfacial tribochemical reactions, and formation mechanism of the tributyltetradecyl phosphonium di(2-ethylhexyl) phosphate ([P44414] [DEHP]) tribo-film. The findings indicate that the friction coefficients of ILs decrease as the normal pressure increases. The combination of normal pressure and shear force promotes the dissociation of ionic liquid chains and enhances the binding of free radicals to the end faces. The friction force is primarily influenced by the number of Fe-C bonds at the interface, which demonstrates a positive correlation with the normal pressure.

Comparative study on the molecular evolution of changqing petroleum coke during air and pure oxygen combustion: A ReaxFF method analysis

Based on reactive force field (ReaxFF) simulation, the molecular evolution mechanism of air combustion simulation and pure oxygen combustion simulation of Changqing petroleum coke was studied. On the basis of the original petroleum coke molecular model, N2 and O2 molecules were added to construct the air combustion simulation model and the pure oxygen combustion simulation model respectively. The number of molecules, main products, and typical hydrocarbon gases in the two simulation results were analyzed by setting up the LAMMPS in. file. Finally, the molecular changes of pyrrole and pyridine were described. The results show that the main reaction interval of pure oxygen combustion is earlier than that of air combustion. In the main reaction temperature range, the reaction rate of pure oxygen combustion is higher than that of simulated air combustion. At 170 ps, the two simulations reaction speed is the fastest time point; Pyridine is more stable than pyrrole. The research on the distinction and comparison of two kinds of simulation has certain guiding significance for the clean utilization of petroleum coke combustion process in industry.

Reaction path identification and validation from molecular dynamics simulations of hydrocarbon pyrolysis

AbstractCreation of complex chemical mechanisms for hydrocarbon pyrolysis and combustion is challenging due to the large number of species and reactions involved. Reactive molecular dynamics (RMD) enables the simulation of thousands of reactions and the discovery of previously unknown components of the reaction network. However, due to the inherent imprecision of reactive force fields, it is necessary to verify RMD‐obtained reaction paths using more accurate methods such as Density Functional Theory (DFT). We demonstrate a method for identification and confirmation of reaction pathways from RMD that supplement an established mechanism, using the example of benzene formation from n‐heptane and iso‐octane pyrolysis. We establish a validation workflow to extract reaction geometries from RMD and optimize transition states using the Nudged‐Elastic‐Band method on semi‐empirical and quantum mechanical levels of theory. Our findings demonstrate that the widely recognized ReaxFF parameterization, CHO2016, can identify known pathways from a established soot formation mechanism while also indicating new ones. We also show that CHO2016 underestimates hydrogen migration barriers by up to as compared to DFT and can lower activation barriers significantly for spin‐forbidden reactions. This highlights the necessity for validation or potentially even reparametrization of CHO2016.

Atomic-Level Insights into Defect-Driven Nitrogen Doping of Reduced Graphene Oxide

Nitrogen-doped graphene has been increasingly utilized in a variety of energy-related applications, serving as a catalyst or support material for fuel cells, and as an anode material for lithium-ion batteries, among others. The thermal reduction of graphene oxide (GO) in nitrogenous sources to incorporate nitrogen, producing nitrogen-doped reduced graphene oxide (NRGO), is the most favored method. Controlling atomic configurations of nitrogen-doped sites is the key factor for tailoring the physico-chemical properties of NRGO, but major challenges remain in identifying detailed atomic arrangements at nitrogen binding sites on highly defective and chemically functionalized GO surfaces. In this paper, we present atomistic-scale modeling of the nitrogen doping process of GO with different types of vacancy defects. Molecular dynamics simulations using a reactive force field indicate that the edge carbon atoms on defect sites are the dominant initiation location for nitrogen doping. Further, first-principles calculations using density functional theory present energetically favorable chemical transition pathways for nitrogen doping. The significance of this work lies in providing important chemical insights for the effective control of the desired properties of NRGO by suggesting a detailed mechanism of the nitrogen doping process of GO.