Role of sampling in evaluating classical time autocorrelation functions
Abstract
We analyze how the choice of the sampling weight affects efficiency of the Monte Carlo evaluation of classical time autocorrelation functions. Assuming uncorrelated sampling or sampling with constant correlation length, we propose a sampling weight for which the number of trajectories needed for convergence is independent of the correlated quantity, dimensionality, dynamics, and phase-space density. By contrast, it is shown that the computational cost of the "standard" algorithm sampling from the phase-space density may scale exponentially with the number of degrees of freedom. Yet, for the stationary Gaussian distribution of harmonic systems and for the autocorrelation function of a linear function of phase-space coordinates, the computational cost of this standard algorithm is also independent of dimensionality.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- September 2013
- DOI:
- 10.1063/1.4820420
- arXiv:
- arXiv:1212.3463
- Bibcode:
- 2013JChPh.139j4105Z
- Keywords:
-
- Gaussian distribution;
- Monte Carlo methods;
- sampling methods;
- 31.15.xv;
- 02.70.Uu;
- Molecular dynamics and other numerical methods;
- Applications of Monte Carlo methods;
- Physics - Chemical Physics;
- Nonlinear Sciences - Chaotic Dynamics;
- Physics - Computational Physics
- E-Print:
- 5 pages, 3 figures, submitted to Phys. Rev. Lett