Path integrals with higher order actions: Application to realistic chemical systems

Analytical evaluations of the path integral Monte Carlo thermodynamic and Hamiltonian energies for the harmonic oscillator

By using the recently derived universal discrete imaginary-time propagator of the harmonic oscillator, both thermodynamic and Hamiltonian energies can be given analytically and evaluated numerically at each imaginary time step for any short-time propagator. This work shows that, using only currently known short-time propagators, the Hamiltonian energy can be optimized to the twelfth-order, converging to the ground state energy of the harmonic oscillator in as few as three beads. This study makes it absolutely clear that the widely used second-order primitive approximation propagator, when used in computing thermodynamic energy, converges extremely slowly with an increasing number of beads.

Anatomy of path integral Monte Carlo: Algebraic derivation of the harmonic oscillator's universal discrete imaginary-time propagator and its sequential optimization

The direct integration of the harmonic oscillator path integral obscures the fundamental structure of its discrete, imaginary time propagator (density matrix). This work, by first proving an operator identity for contracting two free propagators into one in the presence of interaction, derives the discrete propagator by simple algebra without doing any integration. This discrete propagator is universal, having the same two hyperbolic coefficient functions for all short-time propagators. Individual short-time propagator only modifies the coefficient function's argument, its portal parameter, whose convergent order is the same as the thermodynamic energy. Moreover, the thermodynamic energy can be given in a closed form for any short-time propagator. Since the portal parameter can be systematically optimized by matching the expansion of the product of the two coefficients, any short-time propagator can be optimized sequentially, order by order, by matching the product coefficient's expansion alone, without computing the energy. Previous empirical findings on the convergence of fourth and sixth-order propagators can now be understood analytically. An eight-order convergent short-time propagator is also derived.

Generalized fourth-order decompositions of imaginary time path integral: Implications of the harmonic oscillator

The imaginary time path integral formalism offers a powerful numerical tool for simulating thermodynamic properties of realistic systems. We show that, when second-order and fourth-order decompositions are employed, they share a remarkable unified analytic form for the partition function of the harmonic oscillator. We are then able to obtain the expression of the thermodynamic property and the leading error terms as well. In order to obtain reasonably optimal values of the free parameters in the generalized symmetric fourth-order decomposition scheme, we eliminate the leading error terms to achieve the accuracy of desired order for the thermodynamic property of the harmonic system. Such a strategy leads to an efficient fourth-order decomposition that produces third-order accurate thermodynamic properties for general systems.