Chemical Dynamics Simulations of Energy Transfer and Unimolecular Decomposition in Collision‐Induced Dissociation (CID) and Surface‐Induced Dissociation (SID) (original) (raw)
Related papers
Thermal Decomposition of the Solid Phase of Nitromethane: Ab Initio Molecular Dynamics Simulations
The Car-Parrinello molecular dynamics simulations were employed to investigate thermal decomposition of the solid nitromethane. It is found that it undergoes chemical decomposition at about 2200 K under ambient pressure. The initiation of reactions involves both proton transfer and commonly known C-N bond cleavage. About 75 species and 100 elementary reactions were observed with the final products being H 2 O, CO 2 , N 2 , and CNCNC. It represents the first complete simulation of solid-phase explosive reactions reported to date, which is of far-reaching implication for design and development of new energetic materials.
Density-Dependent Liquid Nitromethane Decomposition: Molecular Dynamics Simulations Based on ReaxFF
The Journal of Physical Chemistry A, 2011
The decomposition mechanism of hot liquid nitromethane at various compressions was studied using reactive force field (ReaxFF) molecular dynamics simulations. A competition between two different initial thermal decomposition schemes is observed, depending on compression. At low densities, unimolecular CÀN bond cleavage is the dominant route, producing CH 3 and NO 2 fragments. As density and pressure rise approaching the ChapmanÀJouget detonation conditions (∼30% compression, >2500 K) the dominant mechanism switches to the formation of the CH 3 NO fragment via H-transfer and/or NÀO bond rupture. The change in the decomposition mechanism of hot liquid NM leads to a different kinetic and energetic behavior, as well as products distribution. The calculated density dependence of the enthalpy change correlates with the change in initial decomposition reaction mechanism. It can be used as a convenient and useful global parameter for the detection of reaction dynamics. Atomic averaged local diffusion coefficients are shown to be sensitive to the reactions dynamics, and can be used to distinguish between time periods where chemical reactions occur and diffusiondominated, nonreactive time periods.
Molecular Dynamics Simulation of Dissociation Kinetics
Journal of Thermophysics and Heat Transfer, 2001
The vibrational energy distribution and the degree of dissociation within a system of hydrogen and oxygen molecules was modeled using molecular dynamics (MD). The first step in this process was to model the atomic and molecular interactions. Since hydrogen and
Modelling and Simulation in Materials Science and Engineering, 2015
Ionic liquids based on the dicyanamide anion (DCA) are of interest as replacements for current hypergolic fuels, which are highly toxic. To better understand the reaction dynamics of these ionic liquid fuels, this study reports the results of molecular dynamics simulations performed for two predicted intermediate compounds in DCA-based ionic liquids/nitric acid (HNO 3) combustion, i.e. protonated DCA (DCAH) and nitro-dicyanamide-carbonyl (NDC). Calculations were performed using a ReaxFF reactive force field. Single component simulations show that neat NDC undergo exothermic decomposition and ignition. Simulations with HNO 3 were performed at both a low (0.25 g/ml) and high (1.00 g/ml) densities, to investigate the reaction in a dense vapor and liquid phase, respectively. Both DCAH and NDC react hypergolically with HNO 3 , and increased density led to shorter times for the onset of thermal runaway. Contrary to a proposed mechanism for DCA combustion, neither DCAH nor NDC are converted to 1,5-Dinitrobiuret (DNB) before thermal runaway. Details of reaction pathways for these processes are discussed.
The Journal of Physical Chemistry B, 2011
We studied the thermal decomposition and subsequent reaction of the energetic material nitromethane (CH 3 NO 2) using molecular dynamics with ReaxFF, a first principles-based reactive force field. We characterize the chemistry of liquid and solid nitromethane at high temperatures (2000À3000 K) and density 1.97 g/cm 3 for times up to 200 ps. At T = 3000 K the first reaction in the decomposition of nitromethane is an intermolecular proton transfer leading to CH 3 NOOH and CH 2 NO 2. For lower temperatures (T = 2500 and 2000 K) the first reaction during decomposition is often an isomerization reaction involving the scission of the CÀN bond the formation of a CÀO bond to form methyl nitrate (CH 3 ONO). Also at very early times we observe intramolecular proton transfer events. The main product of these reactions is H 2 O which starts forming following those initiation steps. The appearance of H 2 O marks the beginning of the exothermic chemistry. Recent quantum-mechanics-based molecular dynamics simulations on the chemical reactions and time scales for decomposition of a crystalline sample heated to T = 3000 K for a few picoseconds are in excellent agreement with our results, providing an important, direct validation of ReaxFF.
Molecular Dynamics Simulation of the Low-Temperature Partial Oxidation of CH 4
The Journal of Physical Chemistry A, 2009
Low-temperature partial oxidation of methane was investigated using reactive molecular dynamics (MD) and quantum mechanical (QM) methods. In particular, the ReaxFF hydrocarbon force field et al. J. Phys. Chem. A 2008, 112, 1040 was employed to simulate a [20 CH 4 + 10 O 2 ] model system at 500°C. The chemical mechanism of the partial oxidation of methane was proposed on the basis of analysis of the computed trajectory of this model system. The partial oxidation of methane was observed to be initiated by the abstraction of hydrogen from CH 4 by O 2 and the atomization of CH 4 itself. Subsequent radical recombination between hydrogen atoms and the dehydrogenation of CH 4 were the primary pathways by which H 2 was formed. In agreement with current models of low-temperature combustion, radicals including H 3 C-OO and H 2 C-OO were also observed during the MD simulation. The observed reaction mechanism was subsequently analyzed using QM methods. For instance, structural features of prominent radical species observed during the MD simulation were analyzed using density functional theory (DFT) and coupled-cluster (CCSD(T)) methods. Enthalpies of reaction of all observed chemical processes were calculated using DFT and the W1 composite method. Where possible, comparisons with experimental data were made.
Collision induced processes at the gas–surface interface
Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2002
The collision of energetic gas phase particles with adsorbed species can induce a variety of processes. Events of this kind can play an important role in the mechanisms governing heterogeneous catalysis at high pressures and elevated temperatures. Two collision induced processes (CIP) are described in this article. The first process discussed is collision induced desorption (CID). The CID of N 2 from Ru(001) is considered at both low and high coverage ranges. The interpretation of the experimental data using molecular dynamics (MD) simulations leads to the introduction of a new desorption mechanism involving surface corrugation and adsorbate frustrated rotational motion. The second process is collision induced migration (CIM), an event that has never been considered before neither experimentally nor theoretically. It is demonstrated, using MD simulations, that following energetic CIM, very long distances of more than 100A , can be covered by the adsorbates at low coverages. At high coverages, on the otherhand, these displacements become considerably shorter due to surface collisions with neighbors.
Various bonds in molecules accumulate energies under rising temperature until the energies are sufficient to promote dissociation. The activation energy of a chemical reaction AB →A + B and the bond dissociation energy, D(A-B) are calculated on the basis of the A―B bond stretching vibration ( cm-1). The activation energies, Ea = 883, 407, 249.31, and 437.7 kJmol-1 are calculated for nitrogen, hydrogen chloromethane and water respectively. The dissociation energies D(N-N) = 945.07, D(H-H) = 435.5, while for C2H5Cl the D(C-H) = 409.22 and D(C-Cl) = 341.75 kJmol-1 are also calculated. In each case, the theory confirms the experimental findings.
Reactive molecular dynamic simulations of hydrocarbon dissociations on Ni(111) surfaces
Surface Science, 2012
Empirical potential parameters for H, C and Ni elements have been developed for the ReaxFF force field in order to study the decomposition of small hydrocarbon molecules on nickel using molecular dynamics simulations. These parameters were optimized using the geometrical and energetic information obtained from density functional (DFT) calculations on a subset of hydrogen and methane reactions with nickel (111) surfaces. The resulting force field was then used to obtain a molecular perspective of the dynamics of the methane dissociative adsorption on Ni(111) as well as two other small alkane molecules, ethane and n-butane. NVT simulations of dissociative adsorption of methane over a range of temperatures enabled the estimation of the sticking coefficient for the adsorption as well as the activation energy of the first C-H bond breaking. The rate constants of each elementary step (both forward and reverse) of CH x dissociation on Ni(111) were obtained by monitoring the surface species and a microkinetic model was constructed as a result. Qualitative analyses of the simulations of ethane and n-butane decompositions on Ni(111) demonstrate that such reactive MD technique can also be used to obtain useful information on complex reaction networks.
The mechanism of energy transfer in H2O–H2O collisions – a molecular dynamics simulation
Chemical Physics, 1998
Earlier work on the activation-deactivation mechanism of gas phase unimolecular reactions is extended to the study of the detailed energy transfer mechanism in collisions of water molecules. Molecular dynamics simulations of binary collisions between a reactant water molecule at high internal energy with medium molecules at various selected initial temperatures are Ž. Ž. compared with results from approximate statistical theory. Energy transfer is related to i interaction strength, ii hard Ž. Ž. Ž. atom-atom encounters, iii multiple minima in the center of mass separation, iv collision lifetime and v anharmonicity of the intramolecular potential function. The observed trends are interpreted within the framework of the partially ergodic Ž. multiple encounter theory PEMET of collisional energy transfer. By comparison with typical stable molecule collisions the water-water collisions are more efficient as a reflection of the strong hydrogen bonding interactions. A good agreement between PEMET and molecular dynamics simulations over a wide range of interaction strengths and initial reactant energies is shown, indicating the possibility of a priori use of the PEMET model.