We present a study of fixed- and partial-node approximations in Slater determinant basis sets, using full configuration interaction quantum Monte Carlo (FCIQMC) to perform sampling. Walker annihilation in the FCIQMC method allows partial-node simulations to be performed, relaxing the nodal constraint to converge to the FCI solution. This is applied to ab initio molecular systems, using symmetry-projected Jastrow mean-field wave functions for complete active space (CAS) problems. Convergence and the sign problem within the partial-node approximation are studied, which is shown to eventually be limited in its use due to the large walker populations required. However, the fixed-node approximation results in an accurate and practical method. We apply these approaches to various molecular systems and active spaces, including ferrocene and acenes. This also provides a test of symmetry-projected Jastrow mean-field wave functions in variational Monte Carlo for a new set of problems. For trans -polyacetylene molecules and acenes, we find that the time to perform a constant number of fixed-node FCIQMC iterations scales as and, respectively, resulting in an efficient method for CAS-based problems that can be applied accurately to large active spaces.
We present a study of fixed- and partial-node approximations in Slater determinant basis sets, using full configuration interaction quantum Monte Carlo (FCIQMC) to perform sampling. Walker annihilation in the FCIQMC method allows partial-node simulations to be performed, relaxing the nodal constraint to converge to the FCI solution. This is applied to ab initio molecular systems, using symmetry-projected Jastrow mean-field wave functions for complete active space (CAS) problems. Convergence and the sign problem within the partial-node approximation are studied, which is shown to eventually be limited in its use due to the large walker populations required. However, the fixed-node approximation results in an accurate and practical method. We apply these approaches to various molecular systems and active spaces, including ferrocene and acenes. This also provides a test of symmetry-projected Jastrow mean-field wave functions in variational Monte Carlo for a new set of problems. For trans -polyacetylene molecules and acenes, we find that the time to perform a constant number of fixed-node FCIQMC iterations scales as and, respectively, resulting in an efficient method for CAS-based problems that can be applied accurately to large active spaces.
We have developed a new quantum Monte Carlo method for the simulation of correlated many-electron systems in full configuration-interaction (Slater determinant) spaces. The new method is a population dynamics of a set of walkers, and is designed to simulate the underlying imaginary-time Schrödinger equation of the interacting Hamiltonian. The walkers (which carry a positive or negative sign) inhabit Slater determinant space, and evolve according to a simple set of rules which include spawning, death and annihilation processes. We show that this method is capable of converging onto the full configuration-interaction (FCI) energy and wave function of the problem, without any a priori information regarding the nodal structure of the wave function being provided. Walker annihilation is shown to play a key role. The pattern of walker growth exhibits a characteristic plateau once a critical (system-dependent) number of walkers has been reached. At this point, the correlation energy can be measured using two independent methods--a projection formula and a energy shift; agreement between these provides a strong measure of confidence in the accuracy of the computed correlation energies. We have verified the method by performing calculations on systems for which FCI calculations already exist. In addition, we report on a number of new systems, including CO, O(2), CH(4), and NaH--with FCI spaces ranging from 10(9) to 10(14), whose FCI energies we compute using modest computational resources.
In this thesis we investigate the recently developed Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method. This method, unlike the traditional methods in Quantum Monte Carlo (QMC), doesn’t require the use of the uncontrolled fixed node approximation and so potentially it can yield far more accurate results. This means that it can be thought of as a hybrid between the methods used by quantum chemists and QMC, and thus has spawned a new field of stochastic quantum chemistry. The work described in this thesis can be split into three distinct but interrelated parts. We begin with an investigation of the underlying FCIQMC stochastic process. We show that FCIQMC is an example of Markov Chain Monte Carlo. This means we can compute a stochastic matrix from which all details about the Monte-Carlo simulation can be obtained. Unfortunately the size of the space scales unfavourably as a function of system size meaning that we have only managed to compute the matrix for the smallest interesting two determinant system. We then use these results to quantify population control bias in FCIQMC for a two determinant system supplementing these analytical results with empirical results to investigate more realistic systems. We then attempt to quantify the efficiency of the FCIQMC algorithms, defining a measure of efficiency. After this we investigate the dependence of our measure on the system size. The error bar of the most efficient FCIQMC algorithm will decay fastest as a function of computer time. We then draw conclusions about the applicability of the FCIQMC method. Finally we describe an implementation of FCIQMC on a novel data flow computer architecture. In our implementation we made a modification to the FCIQMC algorithm to fit the data flow paradigm. We investigate the efficiency of the modified FCIQMC algorithm.
Finding the ground state of a fermionic Hamiltonian using quantum Monte Carlo (QMC) is a very difficult problem, due to the Fermi sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCIQMC) mitigates some of the exponential variance by allowing annihilation of noise---whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCIQMC has been quite successful for quantum chemistry problems, its application to problems in condensed systems has been limited. In this paper, we apply the FCIQMC algorithm to the Fermi polaron problem, which provides an ideal test bed for improving the algorithm, although since we restrict the number of many-body excitations, our algorithm is more precisely a QMC implementation of the configuration-interaction method (CIQMC). In its simplest form, CIQMC is unstable for even a fairly small system sizes. However, with a series of algorithmic improvements, we are able to significantly increase its effectiveness. We modify fixed node QMC to work in these systems, introduce a well chosen importance sampled trial wave function, a partial node approximation, and a variant of release node. Finally, we develop a way to perform CIQMC directly in the thermodynamic limit.
Full configuration interaction quantum Monte Carlo for accurate ab initio nuclear structure calculations: algorithms and calculation details
Full configuration interaction quantum Monte Carlo (FCIQMC) is a stochastic many-body solver that has been widely applied to electronic, molecular, and condensed-matter systems. In this work we apply FCIQMC to ab initio nuclear structure calculations using interactions derived from chiral effective field theory. We describe the algorithm in detail, including imaginary-time propagation, excitation generation, estimator choices, the initiator approximation with adaptive shift correction, and reduced-density-matrix (RDM) sampling. Benchmark calculations in small model spaces, where deterministic full configuration interaction (FCI) results are available, validate the stochastic calculation of energies, radii, and RDM-based pure estimators. For large model spaces, we analyze the residual finite-walker bias through systematic walker-number convergence and infinite-walker extrapolations. We also demonstrate that FCIQMC can be extended beyond ground-state calculations by computing the low-lying spectrum of 6Li.
We provide a very simple adaptation of our recently published quantum Monte Carlo algorithm in full configuration-interaction (Slater determinant) spaces which dramatically reduces the number of walkers required to achieve convergence. A survival criterion is imposed for newly spawned walkers. We define a set of initiator determinants such that progeny of walkers spawned from such determinants onto unoccupied determinants are able to survive, while the progeny of walkers not in this set can survive only if they are spawned onto determinants which are already occupied. The set of initiators is originally defined to be all determinants constructable from a subset of orbitals, in analogy with complete-active spaces. This set is dynamically updated so that if a noninitiator determinant reaches an occupation larger than a preset limit, it becomes an initiator. The new algorithm allows sign-coherent sampling of the FCI space to be achieved with relatively few walkers. Using the N(2) molecule as an illustration, we show that rather small initiator spaces and numbers of walkers can converge with submilliHartree accuracy to the known full configuration-interaction (FCI) energy (in the cc-pVDZ basis), in both the equilibrium geometry and the multiconfigurational stretched case. We use the same method to compute the energy with cc-pVTZ and cc-pVQZ basis sets, the latter having an FCI space of over 10(15) with very modest computational resources.
A new quantum Monte Carlo (QMC) method is used to calculate exact, full configuration-interaction (FCI) energies of the neutral and cationic elements from Li to Mg, in a family of commonly used basis sets. Annihilation processes between positive and negative walkers enable the exact N-electron wave function to emerge as a linear superposition of the (factorially large) space of Slater determinants, with individual determinants being stochastically sampled. As a result, extremely large spaces (exceeding 10(15) determinants) become accessible for FCI calculations. No fixed-node approximation is necessary, and the only remaining source of error is the one-electron basis set, which can be systematically reduced by enlargement of the basis set. We have investigated the family of aug-cc-pVXZ Dunning basis sets up to X=5. The resulting ionization potentials are--with one exception (Na)--consistently accurate to within chemical accuracy. The anomalous case of Na suggests that its basis set may be improvable. Extrapolation schemes are examined as a way of further improving the values obtained, and although an improvement is seen in the mean-absolute error, the results of extrapolation are not uniformly better than the largest basis set calculations reported. More generally, these results demonstrate the utility of the QMC method to provide FCI energies for realistic systems and basis sets.