Abstract:A new idea is brought in this paper that molecular dynamics with a force field integrated a brand new reactive potential form which considered bond breaking and forming can be used to simulate chemical reaction. As an example, we report here a force field CRACK developed for studying pyrolysis process of n-alkane. CRACK splines together a harmonic well and a shallow Lennard-Jones well: the harmonic well models C-C covalent bond while the Lennard-Jones well models long-range van der Waals forces. These two wells are separated by an energy barrier which controls dissociation and recombination. To validate this new reactive potential form, CRACK has been used to simulate pyrolysis process of n-decane and n-octane by molecular dynamics (MD). The dissociation probabilities as well as product distributions were analyzed and the time history of pyrolysis was illustrated. The research on dissociation probabilities shows that the defined pyrolytic temperature could be altered by changing the displacement constant de which controls the height of dissociation barrier of C-C bond and it was also found that the pyrolytic temperature of n-octane was higher than of n-decane at the same de, which is consistent with the pyrolysis fundamental. The time history of decane pyrolysis revealed an inspirer result that CRACK force field certainly reflected the recombination between free radicals. Another exciting result is that, by analysis of the distribution of pyrolytic products, the main products distribution was observed to transfer from mid-short to mid-long chains and chain distribution tended to shorter or longer with the elevated temperature which is well consistent with the regularity of thermal cracking. We announce that molecular dynamics can be certainly used to simulate the process of chemical reaction as long as reactive potential is proper.

Accurate p K a Calculations in Proteins with Reactive Molecular Dynamics Provide Physical Insight Into the Electrostatic Origins of Their Values

A brand new reactive potential made molecular dynamics on chemical reaction possible

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Electrostatic Interactions in Protein Molecules: A Novel Model for the Calculation of the Intrinsic Pk a of Lysozyme

Continuum Electrostatics Approaches to Calculating p$K_a$s and $E_m$s in Proteins

An equilibrium constant pH molecular dynamics method for accurate prediction of pH-dependence in protein systems: Theory and application

Modeling chemical reactions in classical molecular dynamics simulations

Reactive Martini: Chemical Reactions in Coarse-Grained Molecular Dynamics Simulations

Accurate biomolecular simulations account for electronic polarization

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Connecting the quantum and classical mechanics simulation world: Applications of reactive step molecular dynamics simulations.

Determination of Amino Acids' p K a : Importance of Cavity Scaling within Implicit Solvation Models and Choice of DFT Functionals

Learning reduced kinetic Monte Carlo models of complex chemistry from molecular dynamics

Recalibrating the calcium trap in amino acid carboxyl groups via classical molecular dynamics simulations

AI2BMD: efficient characterization of protein dynamics with ab initio accuracy

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Accurately Predicting Protein p Ka Values Using Nonequilibrium Alchemy

Application of classical molecular dynamics for evaluation of proton transfer mechanism on a protein

Constant-pH Molecular Dynamics Simulations of Closed and Open States of a Proton-gated Ion Channel

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures