Reactive Molecular Dynamics of the Initial Oxidation Stages of Ni(111) in Pure Water: Effect of an Applied Electric Field

Coupling Process between Droplet and Iron Investigated by Reactive Molecular Dynamics Simulations

The droplet-to-iron electrochemical reaction is common in nature and industrial production, and it causes damage to the economy, safety, and the environment. The electrochemical reaction of droplet-to-iron is a coupling process of wetting and corrosion. Presently, investigations into electrochemical reactions mainly focus on the corrosions caused by a solution, and wetting is rarely considered. However, for the droplet-to-iron electrochemical reaction, the mechanism of charge transfer in the process is still unclear. In this paper, a reactive molecular dynamics simulation model for the droplet-to-iron electrochemical reaction is developed for the first time. The electrochemical reaction of droplet-to-iron is studied, and the interaction between droplet wetting and corrosion on iron is investigated. The effects of temperature, electric field strength, and salt concentration on the electrochemical reaction are explored. Results show that droplet wetting on the iron surface and the formation of a single-molecular-layer ordered structure are prerequisites for corrosion. The hydroxyl radicals that penetrate the ordered structure acquire electrons from iron atoms on the substrate surface under the action of Coulomb forces and form iron-containing oxides with these iron atoms. The corrosion products and craters lead to a reduced droplet height, which promotes droplet wetting on iron and further intensifies the droplet-to-iron electrochemical reaction.

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

The Role of Grain Boundaries in the Corrosion Process of Fe Surface: Insights from ReaxFF Molecular Dynamic Simulations

Intergranular corrosion is the most common corrosion phenomenon in Fe-based alloys. To better understand the mechanism of intergranular corrosion, the influence of grain boundaries on Fe-H2O interfacial corrosion was studied using molecular dynamics simulation based on a new Fe-H2O reaction force field potential. It is found that the corrosion rate at the polycrystalline grain boundary is significantly faster than that of twin crystals and single crystals. By the analysis of stress, it can be found that the stress at the polycrystalline grain boundary and the sigma5 twin grain boundary decreases sharply during the corrosion process. We believe that the extreme stress released at the grain boundary will promote the dissolution of Fe atoms. The formation of vacancies on the Fe matrix surface will accelerate the diffusion of oxygen atoms. This leads to the occurrence of intergranular corrosion.

Water in an External Electric Field: Comparing Charge Distribution Methods Using ReaxFF Simulations

The growing interest in the effects of external electric fields on reactive processes requires predictive methods that can reach longer length and time scales than quantum mechanical simulations. Recently, many studies have included electric fields in ReaxFF, a widely used reactive molecular dynamics method. In the case of modeling an external electric field, the charge distribution method used in ReaxFF is critical. The most common charge distribution method used in previous studies of electric fields is the charge equilibration (QEq) method, which assumes that the system is a contiguous conductor and that charge transfer can occur across any distance. In contrast, many systems of interest are insulators or semiconductors, and long-distance charge transfer should not occur in response to a small difference in potential. This study focuses on the limitations of the QEq method in the context of water in an external electric field. We demonstrate that QEq can predict unphysical charge distributions and exhibits properties that do not converge as a function of system size. Furthermore, we show that electric fields within the recently developed atom-condensed Kohn-Sham density functional theory (DFT) approximated to the second-order (ACKS2) approach address the major limitations of electric fields in QEq. With ACKS2, we observe more physical charge distributions and properties that converge as a function of system size. We do not suggest that ACKS2 is perfect in all circumstances but rather show specific cases where it addresses the major shortcomings of QEq in the context of an external electric field.

A theoretical study on the binding and electrolytic splitting of hydrogen fluoride on Ni(111) and Ni(211)

Anhydrous hydrogen fluoride (HF) is electrochemically activated on Ni electrodes in the Simons process. We study the electrode/electrolyte interface using DFT methods and find the splitting of HF to be exothermic even at low cell potentials.

Effect of electric field on polarization and decomposition of RDX molecular crystals: a ReaxFF molecular dynamics study

The polarization and decomposition of RDX crystal in external electric fields are simulated using molecular dynamics with ReaxFF force field. The results show that the molecular conformation of RDX is transformed from AAE to AAI in the electric field of 0.45 V/Å due to the polarization. In the process of the decomposition of RDX in the electric field, the energy of the system and the exothermic energy of the reaction increase with the increasing of the electric field intensity, and the overall characteristic time of the initial reaction is shorter than that without electric field, which indicates that the external electric field can greatly accelerate the decomposition reaction of RDX.

Role of surface oxidation on the size dependent mechanical properties of nickel nanowires: a ReaxFF molecular dynamics study

Highly reactive metallic nickel (Ni) is readily oxidized by oxygen (O₂) molecules even at low temperatures. The presence of the naturally resulting pre-oxide shell layer on metallic Ni nano materials such as Ni nanowires (NW) is responsible for degrading the deformation mechanisms and related mechanical properties. However, the role of the pre-oxide shell layer on the metallic Ni NW coupled with the complicated mechanical deformation mechanism and related properties have not yet been fully and independently understood. For this reason, the ReaxFF reactive force field for Ni/O interactions was used to investigate the effect of surface oxide layers and the size-dependent mechanical properties of Ni NWs under precisely controlled tensile loading conditions. To directly quantify the size dependent surface oxidation effect on the tensile mechanical deformation behaviour and related properties for Ni NWs, first, ReaxFF-molecular dynamics (MD) simulations were carried out to study the oxidation kinetics on the free surface of Ni NWs in a molecular O₂ environment as a function of various diameters (D = 5.0, 6.5, and 8.0 nm) of the NWs, but at the same length. Single crystalline, pure metallic Ni NWs were also studied as a reference. The results of the oxidation simulations indicate that a surface oxide shell layer with limiting thickness of ∼1.0 nm was formed on the free surface of the bare Ni NW, typically via dissociation of the O-O bonds and the subsequent formation of Ni-O bonds. Furthermore, we investigated the evolution of the size-dependent intrinsic mechanical elastic properties of the core-oxide shell (Ni/Ni_xO_y) NWs by comparing them with their un-oxidized counterparts under constant uniaxial tensile loading. We found that the oxide shell layer significantly decreases the mechanical properties of metallic Ni NW as well as facilitates the initiation of plastic deformation as a function of decreasing diameter. The disordered oxide shell layer on the Ni NW's surface remarkably reduces the yield stress and Young's modulus, due to the increased softening effects with the decreasing NW diameter, compared to un-oxidized counterparts. Moreover, the onset of plastic deformation occurs at a relatively low yielding strain and stress level for the smaller diameter of oxide-coated Ni NWs in comparison to their pure counterparts. Furthermore, for pure Ni NWs, Young's modulus, the yielding stress and strain slightly decrease with the decrease in the diameter size of Ni NWs.