Cumulant expansion and analytic continuation in Monte Carlo simulation of classical Lennard-Jones clusters.

The Monte Carlo (MC) estimates of thermal averages are usually functions of system control parameters λ, such as temperature, volume, and interaction couplings. Given the MC average at a set of prescribed control parameters λ{0}, the problem of analytic continuation of the MC data to λ values in the neighborhood of λ{0} is considered in both classic and quantum domains. The key result is the theorem that links the differential properties of thermal averages to the higher order cumulants. The theorem and analytic continuation formulas expressed via higher order cumulants are numerically tested on the classical Lennard-Jones cluster system of N=13, 55, and 147 neon particles.