Heat capacity estimators for random series path-integral methods by finite-difference schemes

Abstract

Previous heat capacity estimators used in path integral simulations either have large variances that grow to infinity with the number of path variables or require the evaluation of first and second order derivatives of the potential. In the present paper, we show that the evaluation of the total energy by the T-method estimator and of the heat capacity by the TT-method estimator can be implemented by a finite difference scheme in a stable fashion. As such, the variances of the resulting estimators are finite and the evaluation of the estimators requires the potential function only. By comparison with the task of computing the partition function, the evaluation of the estimators requires k + 1 times more calls to the potential, where k is the order of the difference scheme employed. Quantum Monte Carlo simulations for the Ne_13 cluster demonstrate that a second order central-difference scheme should suffice for most applications.

11 pages, 4 figures

References in corpus (11)

• Simulated Tempering: A New Monte Carlo Scheme
• Parallel Tempering Algorithm for Conformational Studies of Biological Molecules
• A Biased Monte Carlo Scheme for Zeolite Structure Solution
• Phase changes in 38 atom Lennard-Jones clusters. I: A parallel tempering study in the canonical ensemble
• The Approach to Ergodicity in Monte Carlo Simulations
• Optimal series representations for numerical path integral simulations
• Optimal Energy Estimation in Path-Integral Monte Carlo Simulations
• Numerical implementation of some reweighted path integral methods
• Energy estimators for random series path-integral methods
• Upon the existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods
• Local Variational Principle

Cited by in corpus (17)

• Structure, superfluidity, and quantum melting of hydrogen clusters
• Thermodynamics and equilibrium structure of Ne_38 cluster: Quantum Mechanics versus Classical
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• The inefficiency of re-weighted sampling and the curse of system size in high order path integration
• Phase changes in selected Lennard-Jones X_{13-n}Y_n clusters
• Direct evaluation of the temperature dependence of the rate constant based on the quantum instanton approximation
• Upon the existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods
• Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene
• Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo
• Imaginary time Gaussian dynamics of the Ar_3 cluster
• Path-integral virial estimator for reaction rate calculation based on the quantum instanton approximation
• First-order corrections to semiclassical Gaussian partition functions for clusters of atoms
• Convergence Characteristics of the Cumulant Expansion for Fourier Path Integrals
• On the efficient Monte Carlo implementation of path integrals
• Geometric Phase Effect in Thermodynamic Properties and in the Imaginary-Time Multi-Electronic-State Path Integral Formulation
• Probing the partition function for temperature-dependent potentials with nested sampling
• Structural information about the ${\text{Ar}}_6$ cluster with the frozen Gaussian imaginary time propagator