Generic parameters of trajectory-extending kinetic Monte Carlo for calculating diffusion coefficients

Wu Jhao Tien, Chi-cheng Chiu

研究成果: Article

2 引文 (Scopus)

摘要

One of the challenging applications of molecular dynamics (MD) simulations is to determine the dynamic properties such as the diffusion coefficient of the molecule of interest, particularly slow dynamic systems such as hydrogels and polymer melts. Recently, Neyertz et al. proposed a trajectory- extending kinetic Monte Carlo (TEKMC) algorithm combining both MD and kinetic Monte Carlo to probe the penetrant diffusion within the glassy polymer systems (S. Neyertz and D. Brown, Macromolecules 43, 9210, 2010). Yet, the original TEKMC relies on the manual adjustments of the key parameters of the sampling interval τ and the discretizing grid size r grid , which limits its applicability to systems with unknown kinetic properties. Here, we reviewed the theoretical background of kinetic Monte Carlo to establish the generic criteria for selecting TEKMC parameters. Also, we modified and expanded the TEKMC algorithm for bulk fluid systems. The modified TEKMC algorithm were applied to systems with various kinetic properties, including Lennard Jones liquid, bulk water, Li + liquid electrolyte, and Li + polymer electrolyte. The diffusion coefficients obtained from the modified TEKMC and the generic parameter selections were promising and robust compared with the conventional MD results. With the proposed TEKMC approach, one can extend the MD trajectories to unambiguously characterize the diffusion behavior in the long-time diffusive regime.

原文English
文章編號065311
期刊AIP Advances
8
發行號6
DOIs
出版狀態Published - 2018 六月 1

指紋

diffusion coefficient
trajectories
kinetics
molecular dynamics
polymers
grids
electrolytes
penetrants
liquids
macromolecules
dynamic characteristics
adjusting
sampling
intervals
probes
fluids
water
molecules

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)

引用此文

@article{a4a70d79cc4e43e2afe73d79b6030995,
title = "Generic parameters of trajectory-extending kinetic Monte Carlo for calculating diffusion coefficients",
abstract = "One of the challenging applications of molecular dynamics (MD) simulations is to determine the dynamic properties such as the diffusion coefficient of the molecule of interest, particularly slow dynamic systems such as hydrogels and polymer melts. Recently, Neyertz et al. proposed a trajectory- extending kinetic Monte Carlo (TEKMC) algorithm combining both MD and kinetic Monte Carlo to probe the penetrant diffusion within the glassy polymer systems (S. Neyertz and D. Brown, Macromolecules 43, 9210, 2010). Yet, the original TEKMC relies on the manual adjustments of the key parameters of the sampling interval τ and the discretizing grid size r grid , which limits its applicability to systems with unknown kinetic properties. Here, we reviewed the theoretical background of kinetic Monte Carlo to establish the generic criteria for selecting TEKMC parameters. Also, we modified and expanded the TEKMC algorithm for bulk fluid systems. The modified TEKMC algorithm were applied to systems with various kinetic properties, including Lennard Jones liquid, bulk water, Li + liquid electrolyte, and Li + polymer electrolyte. The diffusion coefficients obtained from the modified TEKMC and the generic parameter selections were promising and robust compared with the conventional MD results. With the proposed TEKMC approach, one can extend the MD trajectories to unambiguously characterize the diffusion behavior in the long-time diffusive regime.",
author = "Tien, {Wu Jhao} and Chi-cheng Chiu",
year = "2018",
month = "6",
day = "1",
doi = "10.1063/1.5035553",
language = "English",
volume = "8",
journal = "AIP Advances",
issn = "2158-3226",
publisher = "American Institute of Physics Publising LLC",
number = "6",

}

TY - JOUR

T1 - Generic parameters of trajectory-extending kinetic Monte Carlo for calculating diffusion coefficients

AU - Tien, Wu Jhao

AU - Chiu, Chi-cheng

PY - 2018/6/1

Y1 - 2018/6/1

N2 - One of the challenging applications of molecular dynamics (MD) simulations is to determine the dynamic properties such as the diffusion coefficient of the molecule of interest, particularly slow dynamic systems such as hydrogels and polymer melts. Recently, Neyertz et al. proposed a trajectory- extending kinetic Monte Carlo (TEKMC) algorithm combining both MD and kinetic Monte Carlo to probe the penetrant diffusion within the glassy polymer systems (S. Neyertz and D. Brown, Macromolecules 43, 9210, 2010). Yet, the original TEKMC relies on the manual adjustments of the key parameters of the sampling interval τ and the discretizing grid size r grid , which limits its applicability to systems with unknown kinetic properties. Here, we reviewed the theoretical background of kinetic Monte Carlo to establish the generic criteria for selecting TEKMC parameters. Also, we modified and expanded the TEKMC algorithm for bulk fluid systems. The modified TEKMC algorithm were applied to systems with various kinetic properties, including Lennard Jones liquid, bulk water, Li + liquid electrolyte, and Li + polymer electrolyte. The diffusion coefficients obtained from the modified TEKMC and the generic parameter selections were promising and robust compared with the conventional MD results. With the proposed TEKMC approach, one can extend the MD trajectories to unambiguously characterize the diffusion behavior in the long-time diffusive regime.

AB - One of the challenging applications of molecular dynamics (MD) simulations is to determine the dynamic properties such as the diffusion coefficient of the molecule of interest, particularly slow dynamic systems such as hydrogels and polymer melts. Recently, Neyertz et al. proposed a trajectory- extending kinetic Monte Carlo (TEKMC) algorithm combining both MD and kinetic Monte Carlo to probe the penetrant diffusion within the glassy polymer systems (S. Neyertz and D. Brown, Macromolecules 43, 9210, 2010). Yet, the original TEKMC relies on the manual adjustments of the key parameters of the sampling interval τ and the discretizing grid size r grid , which limits its applicability to systems with unknown kinetic properties. Here, we reviewed the theoretical background of kinetic Monte Carlo to establish the generic criteria for selecting TEKMC parameters. Also, we modified and expanded the TEKMC algorithm for bulk fluid systems. The modified TEKMC algorithm were applied to systems with various kinetic properties, including Lennard Jones liquid, bulk water, Li + liquid electrolyte, and Li + polymer electrolyte. The diffusion coefficients obtained from the modified TEKMC and the generic parameter selections were promising and robust compared with the conventional MD results. With the proposed TEKMC approach, one can extend the MD trajectories to unambiguously characterize the diffusion behavior in the long-time diffusive regime.

UR - http://www.scopus.com/inward/record.url?scp=85048618857&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85048618857&partnerID=8YFLogxK

U2 - 10.1063/1.5035553

DO - 10.1063/1.5035553

M3 - Article

AN - SCOPUS:85048618857

VL - 8

JO - AIP Advances

JF - AIP Advances

SN - 2158-3226

IS - 6

M1 - 065311

ER -