TY - JOUR
T1 - Connectivity-based parallel replica dynamics for chemically reactive systems
T2 - From femtoseconds to microseconds
AU - Joshi, Kaushik L.
AU - Raman, Sumathy
AU - Van Duin, Adri C.T.
PY - 2013/11/7
Y1 - 2013/11/7
N2 - Reactive force field methods such as AIREBO, ReaxFF and COMB, are extremely useful for studying physical and chemical interactions between molecules and materials. However, many chemical reactions have relatively high activation energies, putting them beyond the times-scale of conventional molecular dynamics (MD) at modest temperatures. To capture the low-temperature long-lived radical chemistry in atomistic simulations, we have developed a new transition detection scheme for performing Reactive Parallel Replica Dynamics (RPRD) simulation enabling an extended MD time-scales, essentially up to a microsecond using ReaxFF. In the newly implemented event detection scheme, the transition events are identified whenever there is a change in connectivity of any atom. 1-Heptene pyrolysis is chosen as a model system, and RPRD simulations are performed at temperatures as low as 1350K for up to 1 μs for a system consisting of 40 heptene molecules. The chemical mechanism and the product distribution that were obtained from the RPRD simulation are in reasonable agreement with shock tube experiments.
AB - Reactive force field methods such as AIREBO, ReaxFF and COMB, are extremely useful for studying physical and chemical interactions between molecules and materials. However, many chemical reactions have relatively high activation energies, putting them beyond the times-scale of conventional molecular dynamics (MD) at modest temperatures. To capture the low-temperature long-lived radical chemistry in atomistic simulations, we have developed a new transition detection scheme for performing Reactive Parallel Replica Dynamics (RPRD) simulation enabling an extended MD time-scales, essentially up to a microsecond using ReaxFF. In the newly implemented event detection scheme, the transition events are identified whenever there is a change in connectivity of any atom. 1-Heptene pyrolysis is chosen as a model system, and RPRD simulations are performed at temperatures as low as 1350K for up to 1 μs for a system consisting of 40 heptene molecules. The chemical mechanism and the product distribution that were obtained from the RPRD simulation are in reasonable agreement with shock tube experiments.
UR - http://www.scopus.com/inward/record.url?scp=84887778173&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84887778173&partnerID=8YFLogxK
U2 - 10.1021/jz4019223
DO - 10.1021/jz4019223
M3 - Article
AN - SCOPUS:84887778173
SN - 1948-7185
VL - 4
SP - 3792
EP - 3797
JO - Journal of Physical Chemistry Letters
JF - Journal of Physical Chemistry Letters
IS - 21
ER -