Reactions of singly-reduced ethylene carbonate in lithium battery electrolytes: a molecular dynamics simulation study using the ReaxFF

We have conducted quantum chemistry calculations and gas- and solution-phase reactive molecular dynamics simulation studies of reactions involving the ethylene carbonate (EC) radical anion EC(-) using the reactive force field ReaxFF. Our studies reveal that the substantial barrier for transition from the closed (cyclic) form, denoted c-EC(-), of the radical anion to the linear (open) form, denoted o-EC(-), results in a relatively long lifetime of the c-EC(-) allowing this compound to react with other singly reduced alkyl carbonates. Using ReaxFF, we systematically investigate the fate of both c-EC(-) and o-EC(-) in the gas phase and EC solution. In the gas phase and EC solutions with a relatively low concentration of Li(+)/x-EC(-) (where x = o or c), radical termination reactions between radical pairs to form either dilithium butylene dicarbonate (CH(2)CH(2)OCO(2)Li)(2) (by reacting two Li(+)/o-EC(-)) or ester-carbonate compound (by reacting Li(+)/o-EC(-) with Li(+)/c-EC(-)) are observed. At higher concentrations of Li(+)/x-EC(-) in solution, we observe the formation of diradicals which subsequently lead to formation of longer alkyl carbonates oligomers through reaction with other radicals or, in some cases, formation of (CH(2)OCO(2)Li)(2) through elimination of C(2)H(4). We conclude that the local ionic concentration is important in determining the fate of x-EC(-) and that the reaction of c-EC(-) with o-EC(-) may compete with the formation of various alkyl carbonates from o-EC(-)/o-EC(-) reactions.     

ReaxFF(MgH) reactive force field for magnesium hydride systems

We have developed a reactive force field (ReaxFF(MgH)) for magnesium and magnesium hydride systems. The parameters for this force field were derived from fitting to quantum chemical (QM) data on magnesium clusters and on the equations of states for condensed phases of magnesium metal and magnesium hydride crystal. The force field reproduces the QM-derived cell parameters, density, and the equations of state for various pure Mg and MgH(2) crystal phases as well as and bond dissociation, angle bending, charge distribution, and reaction energy data for small magnesium hydride clusters. To demonstrate one application of ReaxFF(MgH), we have carried out MD simulations on the hydrogen absorption/desorption process in magnesium hydrides, focusing particularly on the size effect of MgH(2) nanoparticles on H(2) desorption kinetics. Our results show a clear relationship between grain size and heat of formation of MgH(2); as the particle size decreases, the heat of formation increases. Between 0.6 and 2.0 nm, the heat of formation ranges from -16 to -19 kcal/Mg and diverges toward that of the bulk value (-20.00 kcal/Mg) as the particle diameter increases beyond 2 nm. Therefore, it is not surprising to find that Mg nanoparticles formed by ball milling (20-100 nm) do not exhibit any significant change in thermochemical properties.     

Molecular dynamics simulations of laser-induced incandescence of soot using an extended ReaxFF reactive force field

Laser-induced incandescence (LII) of soot has developed into a popular method for making in situ measurements of soot volume fraction and primary particle sizes. However, there is still a lack of understanding regarding the generation and interpretation of the cooling signals. To model heat transfer from the heated soot particles to the surrounding gas, knowledge of the collision-based cooling as well as reactive events, including oxidation (exothermic) and evaporation (endothermic) is essential. We have simulated LII of soot using the ReaxFF reactive force field for hydrocarbon combustion. Soot was modeled as a stack of four graphene sheets linked together using sp(3) hybridized carbon atoms. To calculate the thermal accommodation coefficient of various gases with soot, graphene sheets of diameter 40 ? were used to create a soot particle containing 2691 atoms, and these simulations were carried out using the ReaxFF version incorporated into the Amsterdam Density Functional program. The reactive force field enables us to simulate the effects of conduction, evaporation, and oxidation of the soot particle on the cooling signal. Simulations were carried out for both reactive and nonreactive gas species at various pressures, and the subsequent cooling signals of soot were compared and analyzed. To correctly model N(2)-soot interactions, optimization of N-N and N-C-H force field parameters against DFT and experimental values was performed and is described in this paper. Subsequently, simulations were performed in order to find the thermal accommodation coefficients of soot with various monatomic and polyatomic gas molecules like He, Ne, Ar, N(2), CO(2), and CH(4). For all these species we find good agreement between our ReaxFF results and previously published accommodation coefficients. We thus believe that Molecular Dynamics using the ReaxFF reactive force field is a promising approach to simulate the physical and chemical aspects of soot LII.     

Development of a ReaxFF reactive force field for ettringite and study of its mechanical failure modes from reactive dynamics simulations

Ettringite is a hexacalcium aluminate trisulfate hydrate mineral that forms during Portland cement hydration. Its presence plays an important role in controlling the setting rate of the highly reactive aluminate phases in cement paste and has also been associated with severe cracking in cured hardened cement. To understand how it forms and how its properties influence those of hardened cement and concrete, we have developed a first-principles-based ReaxFF reactive force field for Ca/Al/H/O/S. Here, we report on the development of this ReaxFF force field and on its validation and application using reactive molecular dynamics (RMD) simulations to characterize and understand the elastic, plastic, and failure response of ettringite at the atomic scale. The ReaxFF force field was validated by comparing the lattice parameters, pairwise distribution functions, and elastic constants of an ettringite crystal model obtained from RMD simulations with those from experiments. The predicted results are in close agreement with published experimental data. To characterize the atomistic failure modes of ettringite, we performed stress-strain simulations to find that Ca-O bonds are responsible for failure of the calcium sulfate and tricalcium aluminate (C3A) column in ettringite during uniaxial compression and tension and that hydrogen bond re-formation during compression induces an increase in plastic strain beyond the material's stress-strain proportionality limit. These results provide essential insight into understanding the mechanistic role of this mineral in cement and concrete degradation, and the ReaxFF potential developed in this work serves as a fundamental tool to further study the kinetics of hydration in cement and concrete.     

ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation

To investigate the initial chemical events associated with high-temperature gas-phase oxidation of hydrocarbons, we have expanded the ReaxFF reactive force field training set to include additional transition states and chemical reactivity of systems relevant to these reactions and optimized the force field parameters against a quantum mechanics (QM)-based training set. To validate the ReaxFF potential obtained after parameter optimization, we performed a range of NVT-MD simulations on various hydrocarbon/O2 systems. From simulations on methane/O2, o-xylene/O2, propene/O2, and benzene/O2 mixtures, we found that ReaxFF obtains the correct reactivity trend (propene > o-xylene > methane > benzene), following the trend in the C-H bond strength in these hydrocarbons. We also tracked in detail the reactions during a complete oxidation of isolated methane, propene, and o-xylene to a CO/CO2/H2O mixture and found that the pathways predicted by ReaxFF are in agreement with chemical intuition and our QM results. We observed that the predominant initiation reaction for oxidation of methane, propene, and o-xylene under fuel lean conditions involved hydrogen abstraction of the methyl hydrogen by molecular oxygen forming hydroperoxyl and hydrocarbon radical species. While under fuel rich conditions with a mixture of these hydrocarbons, we observed different chemistry compared with the oxidation of isolated hydrocarbons including a change in the type of initiation reactions, which involved both decomposition of the hydrocarbon or attack by other radicals in the system. Since ReaxFF is capable of simulating complicated reaction pathways without any preconditioning, we believe that atomistic modeling with ReaxFF provides a useful method for determining the initial events of oxidation of hydrocarbons under extreme conditions and can enhance existing combustion models.     

Lithium salt electride with an excess electron pair--a class of nonlinear optical molecules for extraordinary first hyperpolarizability

A new lithium salt electride with an excess electron pair is designed, for the first time, by means of doping two sodium atoms into the lithium salt of pyridazine. For this series of electride molecules, the structures with all real frequencies and the static first hyperpolarizability (beta 0) are obtained at the second-order M?ller-Plesset theory (MP2). Pyridazine H 4C 4N 2 becomes the lithium salt of pyridazine Li-H 3C 4N 2 as one H atom is substituted by Li. The lithium salt effect on hyperpolarizability is observed as the beta 0 value is increased by about 170 times from 5 to 859 au. For the electride effect, an electride H 4C 4N 2...Na 2 formed by doping two Na atoms into pyridazine, the beta 0 value is increased by about 3000 times from 5 to 1.5 x 10 (4) au. Furthermore, combining these two effects, that is, lithium salt effect and electride effect, more significant increase in beta 0 is expected. A new lithium salt electride Li-H 3C 4N 2...Na 2 is thus designed by doping two Na atoms into Li-H 3C 4N 2. It is found that the new lithium salt electride, Li-H 3C 4N 2...Na 2, has a very large beta 0 value (1.412 x 10 (6) au). The beta 0 value is 2.8 x 10 (5) times larger than that of H 4C 4N 2, 1644 times larger than that of Li-H 3C 4N 2, and still 93 times larger than that of the electride H 4C 4N 2...Na 2. This extraordinary beta 0 value is a new record and comes from its small transition energy and large difference in the dipole moments between the ground state and the excited state. The frequency-dependent beta is also obtained, and it shows almost the same trends as H 4C 4N 2 < Li-H 3C 4N 2 < H 4C 4N 2...Na 2 < Li-H 3C 4N 2...Na 2. This work proposes a new idea to design potential candidate molecules with high-performance NLO properties.     

Effects of fuel additives on the thermal cracking of n-decane from reactive molecular dynamics

Thermal cracking of n-decane and n-decane in the presence of several fuel additives are studied in order to improve the rate of thermal cracking by using reactive molecular dynamics (MD) simulations employing the ReaxFF reactive force field. From MD simulations, we find the initiation mechanisms of pyrolysis of n-decane are mainly through two pathways: (1) the cleavage of a C-C bond to form smaller hydrocarbon radicals, and (2) the dehydrogenation reaction to form an H radical and the corresponding decyl radical. Another pathway is the H-abstraction reactions by small radicals including H, CH(3), and C(2)H(5). The basic reaction mechanisms are in good agreement with existing chemical kinetic models of thermal decomposition of n-decane. Quantum mechanical calculations of reaction enthalpies demonstrate that the H-abstraction channel is easier compared with the direct C-C or C-H bond-breaking in n-decane. The thermal cracking of n-decane with several additives is further investigated. ReaxFF MD simulations lead to reasonable Arrhenius parameters compared with experimental results based on first-order kinetic analysis. The different chemical structures of the fuel additives greatly affect the apparent activation energy and pre-exponential factors. The presence of diethyl ether (DEE), methyl tert-butyl ether (MTBE), 1-nitropropane (NP), 3,6,9-triethyl-3,6,9-trimethyl-1,2,4,5,7,8-hexaoxonane (TEMPO), triethylamine (TEA), and diacetonediperodixe (DADP) exhibit remarkable promoting effect on the thermal cracking rates, compared with that of pure n-decane, in the following order: NP > TEMPO > DADP > DEE (～MTBE) > TEA, which coincides with experimental results. These results demonstrate that reactive MD simulations can be used to screen for fuel additives and provide useful information for more comprehensive chemical kinetic model studies at the molecular level.     

Simulations on the thermal decomposition of a poly(dimethylsiloxane) polymer using the ReaxFF reactive force field

To investigate the failure of the poly(dimethylsiloxane) polymer (PDMS) at high temperatures and pressures and in the presence of various additives, we have expanded the ReaxFF reactive force field to describe carbon-silicon systems. From molecular dynamics (MD) simulations using ReaxFF we find initial thermal decomposition products of PDMS to be CH(3) radical and the associated polymer radical, indicating that decomposition and subsequent cross-linking of the polymer is initiated by Si-C bond cleavage, in agreement with experimental observations. Secondary reactions involving these CH(3) radicals lead primarily to formation of methane. We studied temperature and pressure dependence of PDMS decomposition by following the rate of production of methane in the ReaxFF MD simulations. We tracked the temperature dependency of the methane production to extract Arrhenius parameters for the failure modes of PDMS. Furthermore, we found that at increased pressures the rate of PDMS decomposition drops considerably, leading to the formation of fewer CH(3) radicals and methane molecules. Finally, we studied the influence of various additives on PDMS stability. We found that the addition of water or a SiO(2) slab has no direct effect on the short-term stability of PDMS, but addition of reactive species such as ozone leads to significantly lower PDMS decomposition temperature. The addition of nitrogen monoxide does not significantly alter the degradation temperature but does retard the initial production of methane and C(2) hydrocarbons until the nitrogen monoxide is depleted. These results, and their good agreement with available experimental data, demonstrate that ReaxFF provides a useful computational tool for studying the chemical stability of polymers.     

Parameterization of a reactive force field using a Monte Carlo algorithm

Parameterization of a molecular dynamics force field is essential in realistically modeling the physicochemical processes involved in a molecular system. This step is often challenging when the equations involved in describing the force field are complicated as well as when the parameters are mostly empirical. ReaxFF is one such reactive force field which uses hundreds of parameters to describe the interactions between atoms. The optimization of the parameters in ReaxFF is done such that the properties predicted by ReaxFF matches with a set of quantum chemical or experimental data. Usually, the optimization of the parameters is done by an inefficient single‐parameter parabolic‐search algorithm. In this study, we use a robust metropolis Monte‐Carlo algorithm with simulated annealing to search for the optimum parameters for the ReaxFF force field in a high‐dimensional parameter space. The optimization is done against a set of quantum chemical data for MgSO₄ hydrates. The optimized force field reproduced the chemical structures, the equations of state, and the water binding curves of MgSO₄ hydrates. The transferability test of the ReaxFF force field shows the extend of transferability for a particular molecular system. This study points out that the ReaxFF force field is not indefinitely transferable. © 2013 Wiley Periodicals, Inc.     

First principle and ReaxFF molecular dynamics investigations of formaldehyde dissociation on Fe(100) surface

Detailed formaldehyde adsorption and dissociation reactions on Fe(100) surface were studied using first principle calculations and molecular dynamics (MD) simulations, and results were compared with available experimental data. The study includes formaldehyde, formyl radical (HCO), and CO adsorption and dissociation energy calculations on the surface, adsorbate vibrational frequency calculations, density of states analysis of clean and adsorbed surfaces, complete potential energy diagram construction from formaldehyde to atomic carbon (C), hydrogen (H), and oxygen (O), simulation of formaldehyde adsorption and dissociation reaction on the surface using reactive force field, ReaxFF MD, and reaction rate calculations of adsorbates using transition state theory (TST). Formaldehyde and HCO were adsorbed most strongly at the hollow (fourfold) site. Adsorption energies ranged from ?22.9 to ?33.9 kcal/mol for formaldehyde, and from ?44.3 to ?66.3 kcal/mol for HCO, depending on adsorption sites and molecular direction. The dissociation energies were investigated for the dissociation paths: formaldehyde → HCO + H, HCO → H + CO, and CO → C + O, and the calculated energies were 11.0, 4.1, and 26.3 kcal/mol, respectively. ReaxFF MD simulation results were compared with experimental surface analysis using high resolution electron energy loss spectrometry (HREELS) and TST based reaction rates. ReaxFF simulation showed less reactivity than HREELS observation at 310 and 523 K. ReaxFF simulation showed more reactivity than the TST based rate for formaldehyde dissociation and less reactivity than TST based rate for HCO dissociation at 523 K. TST‐based rates are consistent with HREELS observation. © 2013 Wiley Periodicals, Inc.     

Effect of formic acid addition on water cluster stability and structure

Computational chemistry simulations were performed to determine the effect that the addition of a single formic acid molecule has on the structure and stability of protonated water clusters. Previous experimental studies showed that addition of formic acid to protonated pure water results in higher intensities of large-sized clusters when compared to pure water and methanol-water mixed clusters. For larger, protonated clusters, molecular dynamics simulations were performed on H(+)(H(2)O)(n), H(+)(H(2)O)(n)CH(3)OH, and H(+)(H(2)O)(n)CHOOH clusters, 19-28 molecules in size, using a reactive force field (ReaxFF). Based on these computations, formic acid-water clusters were found to have significantly higher binding energies per molecule. Addition of formic acid to a water cluster was found to alter the structure of the hydrogen-bonding network, creating selective sites within the cluster, enabling the formation of new hydrogen bonds, and increasing both the stability of the cluster and its rate of growth.     

A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part I: Force Field Development

This work presents a novel parametrization for the ReaxFF formalism as a means to investigate reaction processes of chlorinated organic compounds. Force field parameters cover the chemical elements C, H, O, Cl and were obtained using a novel optimization approach involving relaxed potential energy surface scans as training targets. The resulting ReaxFF parametrization shows good transferability, as demonstrated on two independent ab initio validation sets. While this first part of our two-paper series focuses on force field parametrization, we apply our parameters to the simulation of chlorinated dibenzofuran formation and decomposition processes in Part II.     

Atomistic-scale simulations of the initial chemical events in the thermal initiation of triacetonetriperoxide

To study the initial chemical events related to the detonation of triacetonetriperoxide (TATP), we have performed a series of molecular dynamics (MD) simulations. In these simulations we used the ReaxFF reactive force field, which we have extended to reproduce the quantum mechanics (QM)-derived relative energies of the reactants, products, intermediates, and transition states related to the TATP unimolecular decomposition. We find excellent agreement between the QM-predicted reaction products and those observed from 100 independent ReaxFF unimolecular MD cookoff simulations. Furthermore, the primary reaction products and average initiation temperature observed in these 100 independent unimolecular cookoff simulations match closely with those observed from a TATP condensed-phase cookoff simulation, indicating that unimolecular decomposition dominates the thermal initiation of the TATP condensed phase. Our simulations demonstrate that thermal initiation of condensed-phase TATP is entropy-driven (rather than enthalpy-driven), since the initial reaction (which mainly leads to the formation of acetone, O(2), and several unstable C(3)H(6)O(2) isomers) is almost energy-neutral. The O(2) generated in the initiation steps is subsequently utilized in exothermic secondary reactions, leading finally to formation of water and a wide range of small hydrocarbons, acids, aldehydes, ketones, ethers, and alcohols.     

亚烷基卡宾与亚烷基锂氟类卡宾化学键结构的量子拓扑研究

用ab initio(3-21G)方法对亚烷基卡宾H_2C—C的单线态及三线态结构进行了电子密度拓扑分析,说明了它们的亲电、亲核反应方向,讨论了亚烷基锂氟类卡宾H_2C—CLiF的4种构型,论证了该分子中不存在四元环结构、Li—C键以静电作用为主的特性,并预测了加成反应机理。     

Development and application of a ReaxFF reactive force field for hydrogen combustion

To investigate the reaction kinetics of hydrogen combustion at high-pressure and high-temperature conditions, we constructed a ReaxFF training set to include reaction energies and transition states relevant to hydrogen combustion and optimized the ReaxFF force field parameters against training data obtained from quantum mechanical calculations and experimental values. The optimized ReaxFF potential functions were used to run NVT MD (i.e., molecular dynamics simulation with fixed number of atoms, volume, and temperature) simulations for various H(2)/O(2) mixtures. We observed that the hydroperoxyl (HO(2)) radical plays a key role in the reaction kinetics at our input conditions (T ≥ 3000 K, P > 400 atm). The reaction mechanism observed is in good agreement with predictions of existing continuum-scale kinetic models for hydrogen combustion, and a transition of reaction mechanism is observed as we move from high pressure, low temperature to low pressure, high temperature. Since ReaxFF derives its parameters from quantum mechanical data and can simulate reaction pathways without any preconditioning, we believe that atomistic simulations through ReaxFF could be a useful tool in enhancing existing continuum-scale kinetic models for prediction of hydrogen combustion kinetics at high-pressure and high-temperature conditions, which otherwise is difficult to attain through experiments.     

ReaxFF molecular dynamics simulation of thermal stability of a Cu3(BTC)2 metal-organic framework

The thermal stability of a dehydrated Cu(3)(BTC)(2) (copper(II) benzene 1,3,5-tricarboxylate) metal-organic framework was studied by molecular dynamics simulation with a ReaxFF reactive force field. The results show that Cu(3)(BTC)(2) is thermally stable up to 565 K. When the temperature increases between 600 K and 700 K, the framework starts to partially collapse. The RDF analysis shows that the long range correlations between Cu dimers disappear, indicating the loss of the main channels of Cu(3)(BTC)(2). When the temperature is above 800 K, we find the decomposition of the Cu(3)(BTC)(2) framework. CO is the major product, and we also observe the release of CO(2), O(2), 1,3,5-benzenetricarboxylate (C(6)H(3)(CO(2))(3), BTC) and glassy carbon. The Cu dimer is stable up to 1100 K, but we find the formation of new copper oxide clusters at 1100 K. These results are consistent with experimental findings, and provide valuable information for future theoretical investigations of Cu(3)(BTC)(2) and its application in adsorption, separation and catalytic processes.     

Development of a new parameter optimization scheme for a reactive force field based on a machine learning approach

Reactive molecular dynamics (MD) simulation is performed using a reactive force field (ReaxFF). To this end, we developed a new method to optimize the ReaxFF parameters based on a machine learning approach. This approach combines the k-nearest neighbor and random forest regressor algorithm to efficiently locate several possible ReaxFF parameter sets. As a pilot test of the developed approach, the optimized ReaxFF parameter set was applied to perform chemical vapor deposition (CVD) of an α-Al₂O₃ crystal. The crystal structure of α-Al₂O₃ was reasonably reproduced even at a relatively high temperature (2000 K). The reactive MD simulation suggests that the (110) surface grows faster than the (0001) surface, indicating that the developed parameter optimization technique could be used for understanding the chemical reaction in the CVD process. © 2019 Wiley Periodicals, Inc.     

Mechanism and kinetics for the initial steps of pyrolysis and combustion of 1,6-dicyclopropane-2,4-hexyne from ReaxFF reactive dynamics

We report the kinetic analysis and mechanism for the initial steps of pyrolysis and combustion of a new fuel material, 1,6-dicyclopropane-2,4-hexyne, that has enormous heats of pyrolysis and combustion, making it a potential high-energy fuel or fuel additive. These studies employ the ReaxFF force field for reactive dynamics (RD) simulations of both pyrolysis and combustion processes for both unimolecular and multimolecular systems. We find that both pyrolysis and combustion initiate from unimolecular reactions, with entropy-driven reactions being most important in both processes. Pyrolysis initiates with extrusion of an ethylene molecule from the fuel molecule and is followed quickly by isomerization of the fuel molecule, which induces additional radicals that accelerate the pyrolysis process. In the combustion process, we find three distinct mechanisms for the O(2) attack on the fuel molecule: (1) attack on the cyclopropane, ring expanding to form the cyclic peroxide which then decomposes; (2) attack onto the central single bond of the diyne which then fissions to form two C(5)H(5)O radicals; (3) attack on the alkyne-cyclopropane moiety to form a seven-membered ring peroxide which then decomposes. Each of these unimolecular combustion processes releases energy that induces additional radicals to accelerate the combustion process. Here oxygen has major effects both as the radical acceptor and as the radical producer. We extract both the effective activation energy and the effective pre-exponential factor by kinetic analysis of pyrolysis and combustion from these ReaxFF simulations. The low value of the derived effective activation energy (26.18 kcal/mol for pyrolysis and 16.40 kcal/mol for combustion) reveals the high activity of this fuel molecule.