We consider the problem of the extrapolation of path-integral Monte Carlo (PIMC) data to infinite Trotter number P. Finite-P data, being even functions of P, have high-P dependence that is generally well described by a quadratic fit, a0+a1P-2, where a0 is the exact quantum value. However, in order to get convergence it is often necessary to run PIMC codes with rather high P values, which implies long computer times and larger statistical errors of the data. It is well known that also for harmonic systems the finite-P data are not exact; nevertheless, they can be easily calculated by Gaussian quadrature. Starting from this observation, we suggest an easy way to correct PIMC data for anharmonic systems in order to take into account the harmonic part exactly, with strong improvement of the extrapolation to P=. Lower Trotter numbers are thus required, with the advantages of computer-time saving and much better accuracy of the extrapolated values, without any change in the PIMC code. In order to demonstrate the effectiveness of the approach, we report finite-P data processing for a single anharmonic particle, whose finite-P data are obtained by the matrix-squaring method, and for a chain of atoms with Morse interaction.

EXTRAPOLATION TO INFINITE TROTTER NUMBER IN PATH-INTEGRAL MONTE-CARLO SIMULATIONS OF SOLID-STATE SYSTEMS / A. CUCCOLI; A. MACCHI; G. PEDROLLI; V. TOGNETTI; R. VAIA. - In: PHYSICAL REVIEW. B, CONDENSED MATTER. - ISSN 0163-1829. - STAMPA. - 51:(1995), pp. 12369-12379. [10.1103/PhysRevB.51.12369]

EXTRAPOLATION TO INFINITE TROTTER NUMBER IN PATH-INTEGRAL MONTE-CARLO SIMULATIONS OF SOLID-STATE SYSTEMS

CUCCOLI, ALESSANDRO;TOGNETTI, VALERIO;
1995

Abstract

We consider the problem of the extrapolation of path-integral Monte Carlo (PIMC) data to infinite Trotter number P. Finite-P data, being even functions of P, have high-P dependence that is generally well described by a quadratic fit, a0+a1P-2, where a0 is the exact quantum value. However, in order to get convergence it is often necessary to run PIMC codes with rather high P values, which implies long computer times and larger statistical errors of the data. It is well known that also for harmonic systems the finite-P data are not exact; nevertheless, they can be easily calculated by Gaussian quadrature. Starting from this observation, we suggest an easy way to correct PIMC data for anharmonic systems in order to take into account the harmonic part exactly, with strong improvement of the extrapolation to P=. Lower Trotter numbers are thus required, with the advantages of computer-time saving and much better accuracy of the extrapolated values, without any change in the PIMC code. In order to demonstrate the effectiveness of the approach, we report finite-P data processing for a single anharmonic particle, whose finite-P data are obtained by the matrix-squaring method, and for a chain of atoms with Morse interaction.
1995
51
12369
12379
A. CUCCOLI; A. MACCHI; G. PEDROLLI; V. TOGNETTI; R. VAIA
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificatore per citare o creare un link a questa risorsa: https://hdl.handle.net/2158/355061
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 17
  • ???jsp.display-item.citation.isi??? 16
social impact