Estimation of error in the diffusion coefficient from molecular dynamics simulations

被引:70
作者
Chitra, R
Yashonath, S
机构
[1] INDIAN INST SCI,SOLID STATE & STRUCT CHEM UNIT,BANGALORE 560012,KARNATAKA,INDIA
[2] INDIAN INST SCI,SUPERCOMP EDUC & RES CTR,BANGALORE 560012,KARNATAKA,INDIA
来源
JOURNAL OF PHYSICAL CHEMISTRY B | 1997年 / 101卷 / 27期
关键词
D O I
10.1021/jp9703059
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
Long molecular dynamics simulations (120 ns) of guests confined in zeolite NaCaA as well as that of pure fluid have been carried out in the microcanonical ensemble to obtain an estimate of the error associated with the diffusion coefficient (D). It is found that the error is about 50% for a 1 ns run for argon parameters in NaCaA zeolite. It is found that shorter runs lead to an overestimation of the value of D. It is also found that a linear fit to a region of the msd (mean squared displacement) below similar to 200 ps leads to an overestimation of D. For runs that are sufficiently long, the initial configuration has no effect on D. Further, the values obtained from the msd and the velocity autocorrelation function converge for the runs reported here. It is shown that the error in D depends on the nature of the system and its underlying potential energy surface. The calculated statistical inefficiency suggests that the averages of properties such as the total energy will be uncorrelated over blocks longer than 150 ps. Dependence of the error in D on the simulation length and size is reported.
引用
收藏
页码:5437 / 5445
页数:9
相关论文
共 20 条
[1]   STUDIES IN MOLECULAR DYNAMICS .1. GENERAL METHOD [J].
ALDER, BJ ;
WAINWRIGHT, TE .
JOURNAL OF CHEMICAL PHYSICS, 1959, 31 (02) :459-466
[2]  
Allen M. P., 1987, J COMPUTER SIMULATIO, DOI DOI 10.2307/2938686
[3]   DIFFUSION ANOMALY IN SILICALITE AND VPI-5 FROM MOLECULAR-DYNAMICS SIMULATIONS [J].
BANDYOPADHYAY, S ;
YASHONATH, S .
JOURNAL OF PHYSICAL CHEMISTRY, 1995, 99 (12) :4286-4292
[4]  
ENGLESFIELD EJ, 1982, COMPUTER SIMULATION
[5]  
HENSON NJ, 1993, J COMPUT AIDED MATER, V41, P1
[6]   OPTIMIZED INTERMOLECULAR POTENTIAL FUNCTIONS FOR LIQUID HYDROCARBONS [J].
JORGENSEN, WL ;
MADURA, JD ;
SWENSON, CJ .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1984, 106 (22) :6638-6646
[7]   PREDICTION OF LOW OCCUPANCY SORPTION OF ALKANES IN SILICALITE [J].
JUNE, RL ;
BELL, AT ;
THEODOROU, DN .
JOURNAL OF PHYSICAL CHEMISTRY, 1990, 94 (04) :1508-1516
[8]   MOLECULAR STATISTICAL CALCULATION OF THE THERMODYNAMIC ADSORPTION CHARACTERISTICS OF ZEOLITES USING THE ATOM-ATOM APPROXIMATION .2. ADSORPTION OF NON-POLAR AND POLAR INORGANIC MOLECULES BY ZEOLITES OF TYPE-X AND TYPE-Y [J].
KISELEV, AV ;
DU, PQ .
JOURNAL OF THE CHEMICAL SOCIETY-FARADAY TRANSACTIONS II, 1981, 77 :1-15
[9]   EQUATION OF STATE CALCULATIONS BY FAST COMPUTING MACHINES [J].
METROPOLIS, N ;
ROSENBLUTH, AW ;
ROSENBLUTH, MN ;
TELLER, AH ;
TELLER, E .
JOURNAL OF CHEMICAL PHYSICS, 1953, 21 (06) :1087-1092
[10]   Cage-to-cage jumps of xenon in zeolite NaY: Transmission coefficients from molecular dynamics simulations [J].
Mosell, T ;
Schrimpf, G ;
Hahn, C ;
Brickmann, J .
JOURNAL OF PHYSICAL CHEMISTRY, 1996, 100 (11) :4571-4581