A Coupled Molecular Dynamics/Kinetic Monte Carlo Approach for Protonation Dynamics in Extended Systems

2014 ◽  
Vol 10 (10) ◽  
pp. 4221-4228 ◽  
Author(s):  
Gabriel Kabbe ◽  
Christoph Wehmeyer ◽  
Daniel Sebastiani
2015 ◽  
Vol 163 (3) ◽  
pp. A329-A337 ◽  
Author(s):  
Guillaume Blanquer ◽  
Yinghui Yin ◽  
Matias A. Quiroga ◽  
Alejandro A. Franco

2016 ◽  
Vol 18 (18) ◽  
pp. 13052-13065 ◽  
Author(s):  
Emanuel K. Peter ◽  
Joan-Emma Shea ◽  
Igor V. Pivkin

In this paper, we present a coarse replica exchange molecular dynamics (REMD) approach, based on kinetic Monte Carlo (kMC).


MRS Advances ◽  
2016 ◽  
Vol 1 (24) ◽  
pp. 1767-1772 ◽  
Author(s):  
Qian Yang ◽  
Carlos A. Sing-Long ◽  
Evan J. Reed

ABSTRACTKinetic Monte Carlo (KMC) methods have been a successful technique for accelerating time scales and increasing system sizes beyond those achievable with fully atomistic simulations. However, a requirement for its success is a priori knowledge of all relevant reaction pathways and their rate coefficients. This can be difficult for systems with complex chemistry, such as shock-compressed materials at high temperatures and pressures or phenolic spacecraft heat shields undergoing pyrolysis, which can consist of hundreds of molecular species and thousands of distinct reactions. In this work, we develop a method for first estimating a KMC model composed of elementary reactions and rate coefficients by using large datasets derived from a few molecular dynamics (MD) simulations of shock compressed liquid methane, and then using L1 regularization to reduce the estimated chemical reaction network. We find that the full network of 2613 reactions can be reduced by 89% while incurring approximately 9% error in the dominant species (CH4) population. We find that the degree of sparsity achievable decreases when similar accuracy is required for additional populations of species.


2003 ◽  
Vol 532-535 ◽  
pp. 531-535
Author(s):  
S. Baud ◽  
F. Picaud ◽  
C. Ramseyer

2001 ◽  
Vol 672 ◽  
Author(s):  
Sweta Somasi ◽  
Bamin Khomami ◽  
Ronald Lovett

ABSTRACTThe length and time scales of an atomistic simulation are often too small for any direct comparison with experimental observations. In order to study the coverage of pits (COPs) found on the Si (100) surface by epitaxial deposition, we first calculate rate of individual steps using molecular dynamics and then define a sequence of Monte-Carlo steps to study the effect of various factors on effective coverage of COPs.


Sign in / Sign up

Export Citation Format

Share Document