Molecular dynamics (MD) simulation of complex chemistry typically involves thousands of atoms propagating over millions of time steps, generating a wealth of data. Traditionally these data are used to calculate some aggregate properties of the system and then discarded, but we propose that these data can be reused to study related chemical systems. Using approximate chemical kinetic models and methods from statistical learning, we study hydrocarbon chemistries under extreme thermodynamic conditions. We discover that a single MD simulation can contain sufficient information about reactions and rates to predict the dynamics of related yet different chemical systems using kinetic Monte Carlo (KMC) simulation. Our learned KMC models identify thousands of reactions and run 4 orders of magnitude faster than MD. The transferability of these models suggests that we can viably reuse data from existing MD simulations to accelerate future simulation studies and reduce the number of new MD simulations required.

Download full-text PDF

Source
http://dx.doi.org/10.1021/acs.jpca.8b09947DOI Listing

Publication Analysis

Top Keywords

kinetic monte
8
monte carlo
8
thousands reactions
8
molecular dynamics
8
chemical systems
8
transferable kinetic
4
models
4
carlo models
4
models thousands
4
reactions learned
4

Similar Publications

Physical effects of crowdant size and concentration on collective microtubule polymerization.

Biophys J

January 2025

Div. of Biology, IISER Pune, Dr. Homi Bhabha Road, Pashan, Pune 411008, India. Electronic address:

The polymerization of cytoskeletal filaments is regulated by both biochemical pathways, as well as physical factors such as crowding. The effect of crowding in vivo emerges from the density of intracellular components. Due to the complexity of the intracellular environment, most studies are based on either in vitro reconstitution or theory.

View Article and Find Full Text PDF

Gas-phase and water-mediated mechanisms for the OCS + OH reaction.

Phys Chem Chem Phys

January 2025

Departamento de Físico-Química, Instituto de Química - Universidade Federal da Bahia, Rua Barão de Jeremoabo, 147, Salvador, Bahia, 40170-115, Brazil.

We report a computational study of the gas-phase and water-mediated mechanisms for the oxidation of carbonyl sulfide (OCS) by the hydroxyl radical. To achieve reliable results, we employ a dual-level strategy within interpolated single-point energies (VTST-ISPE) at the CCSD(T)/aug-cc-pVTZ//M06-2X/aug-cc-pVTZ level of theory. In the gas-phase mechanism, we have determined the rate constants by kinetic Monte Carlo simulation in the interval of temperatures of 250-550 K.

View Article and Find Full Text PDF

Mechanistic models of humoral kinetics following COVID-19 vaccination.

J R Soc Interface

January 2025

Population Health Sciences, Bristol Medical School, University of Bristol, Oakfield Grove, Bristol, BS8 2BN, UK.

COVID-19 vaccine programmes must account for variable immune responses and waning protection. Existing descriptions of antibody responses to COVID-19 vaccination convey limited information about the mechanisms of antibody production and maintenance. We describe antibody dynamics after COVID-19 vaccination with two biologically motivated mathematical models.

View Article and Find Full Text PDF

Context: Natural fluorapatite (FAP) has been investigated as an adsorbent for the removal of dyes such as methylene blue (MB) and crystal violet (CV) from aqueous solutions. Effective dye removal is crucial for water treatment, particularly for industrial wastewater containing toxic dyes. FAP, a naturally abundant material, was characterized using XRD, FTIR, and SEM analysis.

View Article and Find Full Text PDF

Variable relative biological effectiveness (RBE) of carbon radiotherapy may be calculated using several models, including the microdosimetric kinetic model (MKM), stochastic MKM (SMKM), repair-misrepair-fixation (RMF) model, and local effect model I (LEM), which have not been thoroughly compared. In this work, we compared how these four models handle carbon beam fragmentation, providing insight into where model differences arise. Monoenergetic and spread-out Bragg peak carbon beams incident on a water phantom were simulated using Monte Carlo.

View Article and Find Full Text PDF

Want AI Summaries of new PubMed Abstracts delivered to your In-box?

Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!