Combination of the pair density approximation and the Takahashi–Imada approximation for path integral Monte Carlo simulations
We construct an accurate imaginary time propagator for path integral Monte Carlo simulations for heterogeneous systems consisting of a mixture of atoms and molecules. We combine the pair density approximation, which is highly accurate but feasible only for the isotropic interactions between atoms, with the Takahashi–Imada approximation for general interactions. We present finite temperature simulations results for energy and structure of molecules–helium clusters X{sup 4}He{sub 20} (X=HCCH and LiH) which show a marked improvement over the Trotter approximation which has a 2nd-order time step bias. We show that the 4th-order corrections of the Takahashi–Imada approximation can also be applied perturbatively to a 2nd-order simulation.
- OSTI ID:
- 22570195
- Journal Information:
- Journal of Computational Physics, Vol. 301; Other Information: Copyright (c) 2015 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
High order Chin actions in path integral Monte Carlo
On-the-Fly Generation of Differential Resonance Scattering Probability Distribution Functions for Monte Carlo Codes