skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method

Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Computer Physics Communications
Additional Journal Information:
Journal Volume: 188; Journal Issue: C; Related Information: CHORUS Timestamp: 2016-09-04 15:24:53; Journal ID: ISSN 0010-4655
Country of Publication:

Citation Formats

Gilbreth, C. N., and Alhassid, Y.. Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method. Netherlands: N. p., 2015. Web. doi:10.1016/j.cpc.2014.09.002.
Gilbreth, C. N., & Alhassid, Y.. Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method. Netherlands. doi:10.1016/j.cpc.2014.09.002.
Gilbreth, C. N., and Alhassid, Y.. 2015. "Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method". Netherlands. doi:10.1016/j.cpc.2014.09.002.
title = {Stabilizing canonical-ensemble calculations in the auxiliary-field Monte Carlo method},
author = {Gilbreth, C. N. and Alhassid, Y.},
abstractNote = {},
doi = {10.1016/j.cpc.2014.09.002},
journal = {Computer Physics Communications},
number = C,
volume = 188,
place = {Netherlands},
year = 2015,
month = 3

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.cpc.2014.09.002

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • A model of triangulated random surfaces which is the discrete analog of the Polyakov string is considered. An algorithm is proposed which enables one to study the model by the Monte Carlo method in the grand canonical ensemble. Preliminary results on the determination of the critical index ..gamma.. are presented.
  • In this article the Taylor-expansion method is introduced by which Monte Carlo (MC) simulations in the canonical ensemble can be speeded up significantly, Substantial gains in computational speed of 20-40% over conventional implementations of the MC technique are obtained over a wide range of densities in homogeneous bulk phases. The basic philosophy behind the Taylor-expansion method is a division of the neighborhood of each atom (or molecule) into three different spatial zones. Interactions between atoms belonging to each zone are treated at different levels of computational sophistication. For example, only interactions between atoms belonging to the primary zone immediately surroundingmore » an atom are treated explicitly before and after displacement. The change in the configurational energy contribution from secondary-zone interactions is obtained from the first-order term of a Taylor expansion of the configurational energy in terms of the displacement vector d. Interactions with atoms in the tertiary zone adjacent to the secondary zone are neglected throughout. The Taylor-expansion method is not restricted to the canonical ensemble but may be employed to enhance computational efficiency of MC simulations in other ensembles as well. This is demonstrated for grand canonical ensemble MC simulations of an inhomogeneous fluid which can be performed essentially on a modern personal computer.« less
  • We report on the first successful attempt to apply the auxiliary-field quantum Monte Carlo technique to the calculation of ground-state properties of systems of many electrons interacting via a Coulomb potential. We have been able to substantially reduce the huge statistical fluctuations arising from the repulsive, long-range character of the interactions, which had so far hindered the application of this method to [ital realistic] Hamiltonians for atoms, molecules, and solids. Our technique is demonstrated with calculations of ground-state properties of the simplest molecular and solid-state systems, i.e., the H[sub 2] molecule and the homogeneous electron gas.
  • We extend the recently introduced phaseless auxiliary-field quantum Monte Carlo (QMC) approach to any single-particle basis and apply it to molecular systems with Gaussian basis sets. QMC methods in general scale favorably with the system size as a low power. A QMC approach with auxiliary fields, in principle, allows an exact solution of the Schroedinger equation in the chosen basis. However, the well-known sign/phase problem causes the statistical noise to increase exponentially. The phaseless method controls this problem by constraining the paths in the auxiliary-field path integrals with an approximate phase condition that depends on a trial wave function. Inmore » the present calculations, the trial wave function is a single Slater determinant from a Hartree-Fock calculation. The calculated all-electron total energies show typical systematic errors of no more than a few millihartrees compared to exact results. At equilibrium geometries in the molecules we studied, this accuracy is roughly comparable to that of coupled cluster with single and double excitations and with noniterative triples [CCSD(T)]. For stretched bonds in H{sub 2}O, our method exhibits a better overall accuracy and a more uniform behavior than CCSD(T)« less
  • Cited by 1