A Caltech Library Service

An optimized initialization algorithm to ensure accuracy in quantum Monte Carlo calculations

Fisher, Daniel R. and Kent, David R., IV and Feldmann, Michael T. and Goddard, William A., III (2008) An optimized initialization algorithm to ensure accuracy in quantum Monte Carlo calculations. Journal of Computational Chemistry, 29 (14). pp. 2335-2343. ISSN 0192-8651.

[img] PDF
Restricted to Caltech community only
See Usage Policy.

[img] Image (GIF)
Restricted to Caltech community only
See Usage Policy.


Use this Persistent URL to link to this item:


Quantum Monte Carlo (QMC) calculations require the generation of random electronic configurations with respect to a desired probability density, usually the square of the magnitude of the wavefunction. In most cases, the Metropolis algorithm is used to generate a sequence of configurations in a Markov chain. This method has an inherent equilibration phase, during which the configurations are not representative of the desired density and must be discarded. If statistics are gathered before the walkers have equilibrated, contamination by nonequilibrated configurations can greatly reduce the accuracy of the results. Because separate Markov chains must be equilibrated for the walkers on each processor, the use of a long equilibration phase has a profoundly detrimental effect on the efficiency of large parallel calculations. The stratified atomic walker initialization (STRAW) shortens the equilibration phase of QMC calculations by generating statistically independent electronic configurations in regions of high probability density. This ensures the accuracy of calculations by avoiding contamination by nonequilibrated configurations. Shortening the length of the equilibration phase also results in significant improvements in the efficiency of parallel calculations, which reduces the total computational run time. For example, using STRAW rather than a standard initialization method in 512 processor calculations reduces the amount of time needed to calculate the energy expectation value of a trial function for a molecule of the energetic material RDX to within 0.01 au by 33%.

Item Type:Article
Related URLs:
URLURL TypeDescription
Goddard, William A., III0000-0003-0097-5716
Additional Information:© 2008 Wiley. Received: 2 November 2007; Revised: 11 January 2008; Accepted: 8 February 2008. Published online 12 May 2008. Contract/grant sponsor: Department of Energy; contract/grant number: DE-FG02-97ER25308. Contract/grant sponsor: DOE ASC Project; contract/grant number: B523297. Contract/grant sponsor: ARO-DURIP; contract/grant number: W911NF-07-1-0226. Contract/grant sponsor: ONR-HE; contract/grant number: N00014-05-1-0778. Contract/grant sponsor: ARO-MURI; contract/grant number: W911NF-05-1-0345. Contract/grant sponsor: DARPA-PROM; contract/grant number: N00014-06-1-0938. Contract/grant sponsor: Fannie and John Hertz Foundation.
Funding AgencyGrant Number
Department of EnergyDE-FG02-97ER25308
Department of EnergyB523297
Army Research OfficeW911NF-07-1-0226
Office of Naval ResearchN00014-05-1-0778
Army Research OfficeW911NF-05-1-0345
Defense Advanced Research Projects AgencyN00014-06-1-0938
Fannie and John Hertz FoundationUNSPECIFIED
Subject Keywords:quantum Monte Carlo; walker initialization; parallel computing; parallel efficiency
Issue or Number:14
Record Number:CaltechAUTHORS:FISjcc08
Persistent URL:
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:12218
Deposited By: Archive Administrator
Deposited On:29 Oct 2008 21:13
Last Modified:26 Nov 2019 11:15

Repository Staff Only: item control page