Parametrization of a reactive force field for aluminum hydride

A reactive force field, REAXFF, for aluminum hydride has been developed based on density functional theory (DFT) derived data. REAXFF_(AlH_3) is used to study the dynamics governing hydrogen desorption in AlH_3. During the abstraction process of surface molecular hydrogen charge transfer is found to be well described by REAXFF_(AlH_3). Results on heat of desorption versus cluster size show that there is a strong dependence of the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. In the gas phase, it was observed that small alane clusters agglomerated into a bigger cluster. After agglomeration molecular hydrogen was desorbed from the structure. This thermodynamically driven spontaneous agglomeration followed by desorption of molecular hydrogen provides a mechanism on how mobile alane clusters can facilitate the mass transport of aluminum atoms during the thermal decomposition of NaAlH_4

The ReaxFF reactive force-field : development, applications and future directions

The reactive force-field (ReaxFF) interatomic potential is a powerful computational tool for exploring, developing and optimizing material properties. Methods based on the principles of quantum mechanics (QM), while offering valuable theoretical guidance at the electronic level, are often too computationally intense for simulations that consider the full dynamic evolution of a system. Alternatively, empirical interatomic potentials that are based on classical principles require significantly fewer computational resources, which enables simulations to better describe dynamic processes over longer timeframes and on larger scales. Such methods, however, typically require a predefined connectivity between atoms, precluding simulations that involve reactive events. The ReaxFF method was developed to help bridge this gap. Approaching the gap from the classical side, ReaxFF casts the empirical interatomic potential within a bond-order formalism, thus implicitly describing chemical bonding without expensive QM calculations. This article provides an overview of the development, application, and future directions of the ReaxFF method

Modeling the sorption dynamics of NaH using a reactive force field

We have parametrized a reactive force field for NaH, ReaxFFNaH, against a training set of ab initio derived data. To ascertain that ReaxFFNaH is properly parametrized, a comparison between ab initio heats of formation of small representative NaH clusters with ReaxFFNaH was done. The results and trend of ReaxFFNaH are found to be consistent with ab initio values. Further validation includes comparing the equations of state of condensed phases of Na and NaH as calculated from ab initio and ReaxFFNaH. There is a good match between the two results, showing that ReaxFFNaH is correctly parametrized by the ab initio training set. ReaxFFNaH has been used to study the dynamics of hydrogen desorption in NaH particles. We find that ReaxFFNaH properly describes the surface molecular hydrogen charge transfer during the abstraction process. Results on heat of desorption versus cluster size shows that there is a strong dependence on the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. To gain more insight into the structural transformations of NaH during thermal decomposition, we performed a heating run in a molecular dynamics simulation. These runs exhibit a series of drops in potential energy, associated with cluster fragmentation and desorption of molecular hydrogen. This is consistent with experimental evidence that NaH dissociates at its melting point into smaller fragments

Modeling of Hydrogen Storage Materials: A Reactive Force Field for NaH

Parameterization of a reactive force field for NaH is done using ab initio derived data. The parameterized force field(ReaxFFNaH) is used to study the dynamics governing hydrogen desorption in NaH. During the abstraction process of surface molecular hydrogen charge transfer is found to be well described by the parameterized force field. To gain more insight into the mechanism governing structural transformation of NaH during thermal decomposition a heating run in a molecular dynamics simulation is done. The result shows that a clear signature of hydrogen desorption is the fall in potential energy surface during heating

Thermal properties of fluorinated graphene

Large scale atomistic simulations using the reactive force field approach (ReaxFF) are implemented to investigate the thermomechanical properties of fluorinated graphene (FG). A new set of parameters for the reactive force field potential (ReaxFF) optimized to reproduce key quantum mechanical properties of relevant carbon-fluor cluster systems are presented. Molecular dynamics (MD) simulations are used to investigate the thermal rippling behavior of FG and its mechanical properties and compare them with graphene (GE), graphane (GA) and a sheet of BN. The mean square value of the height fluctuations $$ and the height-height correlation function $H(q)$ for different system sizes and temperatures show that FG is an un-rippled system in contrast to the thermal rippling behavior of graphene (GE). The effective Young's modulus of a flake of fluorinated graphene is obtained to be 273 N/m and 250 N/m for a flake of FG under uniaxial strain along arm-chair and zig-zag direction, respectively.Comment: To appear in Phys. Rev.

Thermal Decomposition of Hydrazines from Reactive Dynamics Using the ReaxFF Reactive Force Field

We report reactive dynamics (RD) studies on: the decomposition of bulk hydrazine (N_2H_4); the decomposition of bulk monomethyl-hydrazine (CH_3N_2H_3), hereafter referred to simply as methyl-hydrazine; the decomposition of hydrazine in the presence of hydrogen peroxide (H_2O_2); and decomposition hydrazine on catalytic surfaces Pt[100] and Pt[111] under various conditions. These studies use the ReaxFF reactive force field to describe the multitude of chemical reactions in these systems for a variety of reaction conditions in order to show that this approach leads to realistic decomposition mechanisms and rates. In particular, we determined how the decomposition of hydrazine is affected by temperature, pressure, and heating rate. We analyzed chemical reaction mechanism of the decomposition of hydrazine at the studied conditions and found that at lower temperatures the initial product from hydrazine decomposition is NH_3, whereas at higher temperatures H_2 and N_2 are the dominant early products. Prominent intermediates observed during these decompositions include N_2H_3, N_2H_2, and NH_2, in agreement with quantum mechanical studies (7.3 ps at 3000 K). As the heating rate is decreased, the onset for hydrazine decomposition shifts to lower temperatures. Using a constant heating rate, we found that higher pressure (increased density) favors formation of NH_3 over N_2 and H_2. In studies of the catalytic decomposition of hydrazine on surfaces Pt[100] and Pt[111], we found that the presence of a Pt-catalyst reduces the initial decomposition temperature of hydrazine by about 50%. We found that the Pt[100]-surface is 20 times more active for hydrazine decomposition than the Pt[111]-surface, in qualitative agreement with experiments. These studies indicate how ReaxFF RD can be useful in understanding the chemical processes involved in bulk and catalytic decomposition and in oxidation of reactive species under various reaction conditions

Multiparadigm modeling of dynamical crack propagation in silicon using a reactive force field

We report a study of dynamic cracking in a silicon single crystal in which the ReaxFF reactive force field is used for several thousand atoms near the crack tip, while more than 100 000 atoms are described with a nonreactive force field. ReaxFF is completely derived from quantum mechanical calculations of simple silicon systems without any empirical parameters. Our results reproduce experimental observations of fracture in silicon including changes in crack dynamics for different crack orientations

Valence energy correction for electron reactive force field

Reactive force fields (ReaxFF) are a classical method to describe material properties based on a bond-order formalism, that allows bond dissociation and consequently investigations of reactive systems. Semiclassical treatment of electrons was introduced within ReaxFF simulations, better known as electron reactive force fields (eReaxFF), to explicitly treat electrons as spherical Gaussian waves. In the original version of eReaxFF, the electrons and electron–holes can lead to changes in both the bond energy and the Coulomb energy of the system. In the present study, the method was modified to allow an electron to modify the valence energy, therefore, permitting that the electron\u27s presence modifies the three-body interactions, affecting the angle among three atoms. When a reaction path involving electron transfer is more sensitive to the geometric configuration of the molecules, corrections in the angular structure in the presence of electrons become more relevant; in this case, bond dissociation may not be enough to describe a reaction path. Consequently, the application of the extended eReaxFF method developed in this work should provide an improved description of a reaction path. As a first demonstration this semiclassical force field was parametrized for hydrogen and oxygen interactions, including water and water\u27s ions. With the modified methodology both the overall accuracy of the force field but also the description of the angles within the molecules in presence of electrons could be improved

Development and Validation of ReaxFF Reactive Force Field for Hydrocarbon Chemistry Catalyzed by Nickel

To enable the study of hydrocarbon reactions catalyzed by nickel surfaces and particles using reactive molecular dynamics on thousands of atoms as a function of temperature and pressure, we have developed the ReaxFF reactive force field to describe adsorption, decomposition, reformation and desorption of hydrocarbons as they interact with the nickel surface. The ReaxFF parameters were determined by fitting to the geometries and energy surfaces from quantum mechanics (QM) calculations for a large number of reaction pathways for hydrocarbon molecules chemisorbed onto nickel (111), (100) and (110) surfaces, supplemented with QM equations of state for nickel and nickel carbides. We demonstrate the validity and accuracy of ReaxFF by applying it to study the reaction dynamics of hydrocarbons as catalyzed by nickel particles and surfaces. For the dissociation of methyl on the (111), (100), and stepped (111) surfaces of nickel, we observe the formation of chemisorbed CH plus subsurface carbide. We observe that the (111) surface is the least reactive, the (100) surface has the fastest reaction rates, and the stepped (111) surface has an intermediate reaction rate. The importance of surface defects in accelerating reaction rates is highlighted by these results