ESTIMATION OF ERROR IN THE DIFFUSION-COEFFICIENT FROM MOLECULAR-DYNAMICS SIMULATIONS

Long molecular dynamics simulations (120 ns) of guests confined in zeo lite 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 abou t 50% for a 1 ns run for argon parameters in NaCaA zeolite. It is foun d 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 dis placement) below similar to 200 ps leads to an overestimation of D. Fo r runs that are sufficiently long, the initial configuration has no ef fect on D. Further, the values obtained from the msd and the velocity autocorrelation function converge for the runs reported here. It is sh own that the error in D depends on the nature of the system and its un derlying potential energy surface. The calculated statistical ineffici ency suggests that the averages of properties such as the total energy will be uncorrelated over blocks longer than 150 ps. Dependence of th e error in D on the simulation length and size is reported.