期刊
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
卷 17, 期 8, 页码 5312-5321出版社
AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.1c00185
关键词
-
资金
- HPCI System Research project [hp190097, hp200129, hp200135]
- MEXT
- MEXT/KAKENHI [19H05645]
Recent years have seen an advancement in molecular dynamics (MD) simulations with larger time steps using the hydrogen-mass-repartitioning (HMR) scheme. By optimizing HMR scaling factors and incorporating accurate temperature/pressure evaluations, stable and reliable MD trajectories with larger time steps have been achieved, providing more efficiency in comparison to conventional methods.
In recent years, molecular dynamics (MD) simulations with larger time steps have been performed with the hydrogen-mass-repartitioning (HMR) scheme, where the mass of each hydrogen atom is increased to reduce high-frequency motion while the mass of a non-hydrogen atom bonded to a hydrogen atom is decreased to keep the total molecular mass unchanged. Here, we optimize the scaling factors in HMR and combine them with previously developed accurate temperature/pressure evaluations. The heterogeneous HMR scaling factors are useful to avoid the structural instability of amino acid residues having a five- or six-membered ring in MD simulations with larger time steps. It also reproduces kinetic properties, namely translational and rotational diffusions, if the HMR scaling factors are applied to only solute molecules. The integration scheme is tested for biological systems that include soluble/membrane proteins and lipid bilayers for about 200 mu s MD simulations in total and give consistent results in MD simulations with both a small time step of 2.0 fs and a large, multiple time step integration with time steps of 3.5 fs (for fast motions) and 7.0 fs (for slower motions). We also confirm that the multiple time step integration scheme used in this study provides more accurate energy conservations than the RESPA/C1 and is comparable to the RESPA/C2 in NAMD. In summary, the current integration scheme combining the optimized HMR with accurate temperature/pressure evaluations can provide stable and reliable MD trajectories with a larger time step, which are computationally more than 2-fold efficient compared to the conventional methods.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据