The construction of importance sampled diffusion Monte Carlo (DMC) schemes accurate to second order in the time step is discussed. A central aspect in obtaining efficient second order schemes is the numerical solution of the stochastic differential equation (SDE) associated with the Fokker-Plank equation responsible for the importance sampling procedure. In this work, stochastic predictor-corrector schemes solving the SDE and consistent with Itô calculus are used in DMC simulations of helium clusters. These schemes are numerically compared with alternative algorithms obtained by splitting the Fokker-Plank operator, an approach that we analyze using the analytical tools provided by Itô calculus. The numerical results show that predictor-corrector methods are indeed accurate to second order in the time step and that they present a smaller time step bias and a better efficiency than second order split-operator derived schemes when computing ensemble averages for bosonic systems. The possible extension of the predictor-corrector methods to higher orders is also discussed.

Improved diffusion Monte Carlo propagators for bosonic systems using Ito calculus

MELLA, MASSIMO;BRESSANINI, DARIO;MOROSI, GABRIELE;
2006-01-01

Abstract

The construction of importance sampled diffusion Monte Carlo (DMC) schemes accurate to second order in the time step is discussed. A central aspect in obtaining efficient second order schemes is the numerical solution of the stochastic differential equation (SDE) associated with the Fokker-Plank equation responsible for the importance sampling procedure. In this work, stochastic predictor-corrector schemes solving the SDE and consistent with Itô calculus are used in DMC simulations of helium clusters. These schemes are numerically compared with alternative algorithms obtained by splitting the Fokker-Plank operator, an approach that we analyze using the analytical tools provided by Itô calculus. The numerical results show that predictor-corrector methods are indeed accurate to second order in the time step and that they present a smaller time step bias and a better efficiency than second order split-operator derived schemes when computing ensemble averages for bosonic systems. The possible extension of the predictor-corrector methods to higher orders is also discussed.
2006
http://scitation.aip.org/content/aip/journal/jcp/125/18/10.1063/1.2371077
STOCHASTIC DIFFERENTIAL-EQUATIONS, NUMERICAL-INTEGRATION, RANDOM-WALK, SIMULATIONS, ALGORITHMS, OPERATORS, CLUSTERS, HELIUM
Hakansson, P.; Mella, Massimo; Bressanini, Dario; Morosi, Gabriele; Patrone, M.
File in questo prodotto:
File Dimensione Formato  
hakansson-jcp2006.pdf

accesso aperto

Descrizione: PDF editoriale
Tipologia: Altro materiale allegato
Licenza: Creative commons
Dimensione 446.29 kB
Formato Adobe PDF
446.29 kB Adobe PDF Visualizza/Apri

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

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11383/1744692
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 17
  • ???jsp.display-item.citation.isi??? 17
social impact