Back to Search
Start Over
Multi-level Monte Carlo path integral molecular dynamics for thermal average calculation in the nonadiabatic regime
- Publication Year :
- 2020
-
Abstract
- With the path integral approach, the thermal average in a multi-electronic-state quantum systems can be approximated by the ring polymer representation on an extended configuration space, where the additional degrees of freedom are associated with the surface index of each bead. The primary goal of this work is to propose a more efficient sampling algorithm for the calculation of such thermal averages. We reformulate the extended ring polymer approximation according to the configurations of the surface indexes, and by introducing a proper reference measure, the reformulation is recast as a ratio of two expectations of function expansions. By quantitatively estimating the sub-estimators, and minimizing the total variance of the sampled average, we propose a multi-level Monte Carlo path integral molecular dynamics method (MLMC-PIMD) to achieve an optimal balance of computational cost and accuracy.
- Subjects :
- Mathematics - Numerical Analysis
65C05, 82C31, 82B10
Subjects
Details
- Database :
- arXiv
- Publication Type :
- Report
- Accession number :
- edsarx.2006.09138
- Document Type :
- Working Paper