Abstract:In many scientific fields, there is an interest in understanding the way in which chemical networks evolve. The chemical networks which researchers focus upon have become increasingly complex, and this has motivated the development of automated methods for exploring chemical reactivity or conformational change in a "black-box" manner, harnessing modern computing resources to automate mechanism discovery. In this work, we present a new approach to automated mechanism generation which couples molecular dynamics and statistical rate theory to automatically find kinetically important reactions and then solve the time evolution of the species in the evolving network. The key to this chemical network mapping through combined dynamics and ME simulation approach is the concept of "kinetic convergence", whereby the search for new reactions is constrained to those species which are kinetically favorable at the conditions of interest. We demonstrate the capability of the new approach for two systems, a well-studied combustion system and a multiple oxygen addition system relevant to atmospheric aerosol formation.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jctc.1c00335?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jctc.1c00335</a>.Dynamical roaming trajectory (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.1c00335/suppl_file/ct1c00335_si_001.mpg">MPG</a>)IRC for concerted 1 MHE reaction (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.1c00335/suppl_file/ct1c00335_si_002.mpg">MPG</a>)TSSCDS algorithm details; current interfaces of ChemDyME with electronic structure theory packages; comparison of BXDE products at different temperatures; and full MESMER input files for each system (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.1c00335/suppl_file/ct1c00335_si_003.pdf">PDF</a>)This article has not yet been cited by other publications.

Analyzing Complex Reaction Mechanisms Using Path Sampling

Integrated Reaction Path Processing from Sampled Structure Sequences

Improved Sampling of Adaptive Path Collective Variables by Stabilized Extended-System Dynamics

Discover, Sample, and Refine: Exploring Chemistry with Enhanced Sampling Techniques

Entropic Path Sampling: Computational Protocol to Evaluate Entropic Profile along a Reaction Path

Combining transition path sampling with data-driven collective variables through a reactivity-biased shooting algorithm

Reaction Sampling and Reactivity Prediction Using the Stochastic Surface Walking Method

From Enhanced Sampling to Reaction Profiles

Enhanced sampling reveals the main pathway of an organic cascade reaction

Calculation of reaction constants using Transition Path Sampling with a local Lyapunov bias

Rich Dynamics Underlying Solution Reactions Revealed by Sampling and Data Mining of Reactive Trajectories

Constrained Broyden Dimer Method with Bias Potential for Exploring Potential Energy Surface of Multistep Reaction Process.

Computing Surface Reaction Rates by Adaptive Multilevel Splitting Combined with Machine Learning and Ab Initio Molecular Dynamics

Integrating Path Sampling with Enhanced Sampling for Rare-event Kinetics

Stochastic Surface Walking Reaction Sampling For Resolving Heterogeneous Catalytic Reaction Network: A Revisit To The Mechanism Of Water-Gas Shift Reaction On Cu

Combined Molecular Dynamics and Coordinate Driving Method for Automatic Reaction Pathway Search of Reactions in Solution.

ChemDyME: Kinetically Steered, Automated Mechanism Generation through Combined Molecular Dynamics and Master Equation Calculations

Multi-level Protocol for Mechanistic Reaction Studies Using Semi-local Fitted Potential Energy Surfaces

Differentiable Simulations for Enhanced Sampling of Rare Events

Effective Data-Driven Collective Variables for Free Energy Calculations from Metadynamics of Paths