The objective of pure-sampling quantum Monte Carlo is to calculate physical properties that are independent of the importance sampling function being employed in the calculation, save for the mismatch of its nodal hypersurface with that of the exact wave function. To achieve this objective, we describe a pure-sampling algorithm that combines features of forward-walking methods of pure-sampling and reptation quantum Monte Carlo. The importance sampling is performed by using a single-determinant basis set composed of Slater-type orbitals. We implement our algorithm by systematically increasing an algorithmic parameter until the properties sampled from the electron distributions converge to statistically equivalent values, extrapolated in the limit of zero time-step. In doing so, we are able to unambiguously determine the values for the ground-state fixed-node energies and one-electron properties of various molecules. These quantities are free from importance sampling bias, population control bias, time-step bias, extrapolation-model bias, and the finite-field approximation. We applied our algorithm to the ground-states of lithium hydride, water and ethylene molecules, and found excellent agreement with the accepted literature values for the energy and a variety of other properties for those systems. Some of our one-electron properties of ethylene had not been calculated before at any level of theory. In a detailed comparison, we found reptation quantum Monte Carlo, our closest competitor, to be less efficient by at least a factor of two. It requires different sets of time-steps to accurately determine the ground-state energy and one-electron properties, whereas our algorithm can achieve the same objective by using a single set of time-step values.
The objective of pure-sampling quantum Monte Carlo is to calculate physical properties that are independent of the importance sampling function being employed in the calculation, save for the mismatch of its nodal hypersurface with that of the exact wave function. To achieve this objective, we describe a pure-sampling algorithm that combines features of forward-walking methods of pure-sampling and reptation quantum Monte Carlo. The importance sampling is performed by using a single-determinant basis set composed of Slater-type orbitals. We implement our algorithm by systematically increasing an algorithmic parameter until the properties sampled from the electron distributions converge to statistically equivalent values, extrapolated in the limit of zero time-step. In doing so, we are able to unambiguously determine the values for the ground-state fixed-node energies and one-electron properties of various molecules. These quantities are free from importance sampling bias, population control bias, time-step bias, extrapolation-model bias, and the finite-field approximation. We applied our algorithm to the ground-states of lithium hydride, water and ethylene molecules, and found excellent agreement with the accepted literature values for the energy and a variety of other properties for those systems. Some of our one-electron properties of ethylene had not been calculated before at any level of theory. In a detailed comparison, we found reptation quantum Monte Carlo, our closest competitor, to be less efficient by at least a factor of two. It requires different sets of time-steps to accurately determine the ground-state energy and one-electron properties, whereas our algorithm can achieve the same objective by using a single set of time-step values.
The objective of pure-sampling quantum Monte Carlo is to calculate physical properties that are independent of the importance sampling function being employed in the calculation, save for the mismatch of its nodal hypersurface with that of the exact wave function. To achieve this objective, we report a pure-sampling algorithm that combines features of forward walking methods of pure-sampling and reptation quantum Monte Carlo (RQMC). The new algorithm accurately samples properties from the mixed and pure distributions simultaneously in runs performed at a single set of time-steps, over which extrapolation to zero time-step is performed. In a detailed comparison, we found RQMC to be less efficient. It requires different sets of time-steps to accurately determine the energy and other properties, such as the dipole moment. We implement our algorithm by systematically increasing an algorithmic parameter until the properties converge to statistically equivalent values. As a proof in principle, we calculated the fixed-node energy, static α polarizability, and other one-electron expectation values for the ground-states of LiH and water molecules. These quantities are free from importance sampling bias, population control bias, time-step bias, extrapolation-model bias, and the finite-field approximation. We found excellent agreement with the accepted values for the energy and a variety of other properties for those systems.
Abstract This paper introduces a successful new method for extending diffusion QMC with importance sampling to allow the direct calculation of pure ‘¢ 2 distributions in place of mixed ‘¢ 1/Jr distributions. The method removes any bias introduced with a trial function (except that due to node locations) in the calculation of observable quantities which do not commute with the Hamiltonian. It has the advantage, relative to the descendent weighting or forward walking technique,a of a controlled walker population which reduces statistical errors. The authors have named it the reptation quantum Monte Carlo (RQMC) method. Its basic object is a path or sequence of configurations called a “reptile” which may be altered by an action called “reptation,” terms adopted from related work in the polymer area. In conventional diffusion QMC with importance sampling walkers are subjected to diffusion, drift, and a local energy-based multiplication to produce a 1/J’I/Jr distribution. Without multiplication a ‘l/Jr’l/Jr distribution is produced. In RQMC multiplication is replaced by a Metropolis rejection procedure, also based on local energies, which produces a 1/J’!/J distribution of configurations in a reptile. The generation of a new reptile, by removing the tail and adding a new head, is based on diffusion and drift for proposed moves coupled with an acceptance test. Overall, the implementation is quite similar to that for variational QMC. As described, the method was tested first for the hydrogen atom, followed by the successful demonstration of a calculation for a helium fluid, 64 4He atoms in a cubic box with periodic boundary conditions. The advantages of RQMC relative to branching diffusion QMC were pointed out for treatment of clusters, films, and superfluids as well as for estimation of electronic forces.
This work investigates mathematical details and computational aspects of Metropolis-Hastings reptation quantum Monte Carlo and its variants, in addition to the Bounce method and its variants. The issues that concern us include the sensitivity of these algorithms' target densities to the position of the trial electron density along the reptile, time-reversal symmetry of the propagators, and the length of the reptile. We calculate the ground-state energy and one-electron properties of LiH at its equilibrium geometry for all these algorithms. The importance sampling is performed with a single-determinant large Slater-type orbitals (STO) basis set. The computer codes were written to exploit the efficiencies engineered into modern, high-performance computing software. Using the Bounce method in the calculation of non-energy-related properties, those represented by operators that do not commute with the Hamiltonian, is a novel work. We found that the unmodified Bounce gives good ground state energy and very good one-electron properties. We attribute this to its favourable time-reversal symmetry in its target density's Green's functions. Breaking this symmetry gives poorer results. Use of a short reptile in the Bounce method does not alter the quality of the results. This suggests that in future applications one can use a shorter reptile to cut down the computational time dramatically.
Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.
… as suggested by Fermi, the time-independent Schrödinger equation … can be interpreted as describing the behavior of a system of particles each of which performs a random walk, i.e., diffuses isotropically and at the same time is subject to multiplication, which is determined by the value of the point function V. N. Metropolis and S. Ulam, 1949 Summary In the projector quantum Monte Carlo method, one uses a function of the hamiltonian to sample a distribution proportional to the exact ground-state wavefunction, and thereby computes exact matrix elements of it. An importance sampling transformation makes the algorithm much more efficient. In this chapter we introduce and develop the diffusion Monte Carlo method, which involves drifting, branching random walks. For any excited state, including any system with more than two electrons, one encounters the sign problem, limiting the direct application of these algorithms for most fermion systems. Instead, by using approximate fixed-node or fixed-phase boundary conditions, one can achieve efficiency similar to variational Monte Carlo. We also discuss the application of the projector method in a basis of Slater determinants. In this chapter, we discuss a different quantum Monte Carlo method, projector Monte Carlo (PMC). This general method was first suggested by Fermi [1049]; see the quote at the start of this chapter by two of the inventors of the Monte Carlo method. An implementation of PMC was tried out in the early days of computing [1050]. Advances in methodology, in particular importance sampling, resulted in a significant large-scale application: the exact calculation of the ground-state properties of 256 hard-sphere bosons by Kalos, Levesque, and Verlet [1051] in 1974. Calculations for electronic systems and the fixednode approximation were introduced by Anderson [1052, 1053]. One of the most important projector MC algorithms, the diffusion Monte Carlo algorithm with importance sampling for fermions, was used to compute the correlation energy of homogeneous electron gas by Ceperley and Alder [109] in 1980; the resulting HEG correlation energy was crucial in the development of density functional calculations. Types and properties of projectors In this method, a many-body projector G ( R , R ) = Ĝ is repeatedly applied to filter out the exact many-body ground state from an initial state; the operation of the projector is carried out with a random walk, hence the name of this class of methods.
Fermion quantum Monte Carlo (QMC) methods that work in a general basis space may be more effective for some problems than the traditional diffusion QMC method. Projection of the ground state energy is achieved by random walks in the space of Slater determinants whose one‐particle orbitals are expressed in terms of the chosen basis set. A complication is that the determinants will in general acquire complex phases. This a consequence of ground state Monte Carlo projection using the Hubbard‐Stratonovich transformation of the two‐body interaction. To control the resulting “sign” decay, we describe a method we have recently introduced for the propagation of phaseless determinants. The approximation relies on importance sampling with a trial wave function. The approximation has features in common with diffusion MC fixed node and lattice‐model constrained path methods. Using a planewave basis and non‐local pseudopotentials, we apply the method to Si atom, dimer, and 2, 16, 54 atom (216 electrons) bulk supercells. Single Slater determinant wave functions from density functional theory calculations were used as |ΨT〉 with no additional optimization. The calculated binding energy of Si2 and cohesive energy of bulk Si are in excellent agreement with experiments and are comparable to the best existing theoretical results.