Development of a Mg/O ReaxFF Potential to describe the Passivation Processes in Magnesium-Ion Batteries

Reaction Dynamics of CO2 Hydrogenation on Iron Catalysts Using ReaxFF Molecular Dynamics Simulation

The conversion of CO2 to hydrocarbons using catalysts is a promising route to utilize CO2 and produce more valuable chemicals in a sustainable manner. Recent studies have shown that iron-based catalysts perform well for the hydrogenation of CO2. While the hydrogenation reaction mechanism in the gas phase is straightforward, when catalyzed by iron it has been demonstrated to involve various chemical transformations, and the selectivity and conversion are strongly dependent on the particle size. To further investigate the dependence of the reactivity of iron catalysts on cluster size, we performed reactive molecular dynamics simulations using the ReaxFF force field (ReaxFF-MD) for iron nanoclusters of various sizes in a CO2 and H2-rich environment. We demonstrated that the homogeneous hydrogenation of CO2 was correctly described by this ReaxFF model. The dissociation mechanism of CO2 on the Fe4, Fe16 clusters, and the bcc(100) Fe slab agrees with previous DFT results. The ReaxFF-MD simulations suggest a strong dependence of reactivity on the cluster size, with the Fe4 cluster having the highest reactivity. We show that ReaxFF-MD provides a route to understand reaction mechanisms in these nonequilibrium reactive processes where fast processes and local minima are important.

Surface Transfer Doping in MoO3-x/Hydrogenated Diamond Heterostructure

Surface transfer doping is proposed to be a potential solution for doping diamond, which is hard to dope for applications in high-power electronics. While MoO3 is found to be an effective surface electron acceptor for hydrogen-terminated diamond with a negative electron affinity, the effects of commonly existing oxygen vacancies remain elusive. We have performed reactive molecular dynamics simulations to study the deposition of MoO3-x on a hydrogenated diamond (111) surface and used first-principles calculations based on density functional theory to investigate the electronic structures and charge transfer mechanisms. We find that MoO3-x is an effective surface electron acceptor and the spatial extent of doped holes in hydrogenated diamond is extended, promoting excellent transport properties. Charge transfer is found to monotonically decrease with the level of oxygen vacancy, providing guidance for engineering of the surface transfer doping process.

Enhancing ReaxFF for molecular dynamics simulations of lithium-ion batteries: an interactive reparameterization protocol

Lithium-ion batteries (LIBs) have become an essential technology for the green economy transition, as they are widely used in portable electronics, electric vehicles, and renewable energy systems. The solid-electrolyte interphase (SEI) is a key component for the correct operation, performance, and safety of LIBs. The SEI arises from the initial thermal metastability of the anode-electrolyte interface, and the resulting electrolyte reduction products stabilize the interface by forming an electrochemical buffer window. This article aims to make a first-but important-step towards enhancing the parametrization of a widely-used reactive force field (ReaxFF) for accurate molecular dynamics (MD) simulations of SEI components in LIBs. To this end, we focus on Lithium Fluoride (LiF), an inorganic salt of great interest due to its beneficial properties in the passivation layer. The protocol relies heavily on various Python libraries designed to work with atomistic simulations allowing robust automation of all the reparameterization steps. The proposed set of configurations, and the resulting dataset, allow the new ReaxFF to recover the solid nature of the inorganic salt and improve the mass transport properties prediction from MD simulation. The optimized ReaxFF surpasses the previously available force field by accurately adjusting the diffusivity of lithium in the solid lattice, resulting in a two-order-of-magnitude improvement in its prediction at room temperature. However, our comprehensive investigation of the simulation shows the strong sensitivity of the ReaxFF to the training set, making its ability to interpolate the potential energy surface challenging. Consequently, the current formulation of ReaxFF can be effectively employed to model specific and well-defined phenomena by utilizing the proposed interactive reparameterization protocol to construct the dataset. Overall, this work represents a significant initial step towards refining ReaxFF for precise reactive MD simulations, shedding light on the challenges and limitations of ReaxFF force field parametrization. The demonstrated limitations emphasize the potential for developing more versatile and advanced force fields to upscale ab initio simulation through our interactive reparameterization protocol, enabling more accurate and comprehensive MD simulations in the future.