We propose to accelerate convergence toward full configuration interaction (FCI) energetics by using the coupled-cluster approach, in which singly and doubly excited clusters, needed to determine the energy, are iterated in the presence of their three- and four-body counterparts extracted from FCI quantum Monte Carlo (FCIQMC) propagations. Preliminary calculations for the water molecule at the equilibrium and stretched geometries show that we can accurately extrapolate the FCI energetics based on the early stages of FCIQMC propagations.
We propose to accelerate convergence toward full configuration interaction (FCI) energetics by using the coupled-cluster approach, in which singly and doubly excited clusters, needed to determine the energy, are iterated in the presence of their three- and four-body counterparts extracted from FCI quantum Monte Carlo (FCIQMC) propagations. Preliminary calculations for the water molecule at the equilibrium and stretched geometries show that we can accurately extrapolate the FCI energetics based on the early stages of FCIQMC propagations.
The convergence of full configuration interaction quantum Monte Carlo (FCIQMC) is accelerated using a quasi-Newton propagation (QN) which can also be applied to coupled cluster Monte Carlo (CCMC). The computational scaling of this optimized propagation is, keeping the additional computational cost to a bare minimum. Its effects are investigated deterministically and stochastically on a model system, the uniform electron gas, with Hilbert space size up to 10 40, and shown to accelerate convergence of the instantaneous projected energy by over an order of magnitude in the FCIQMC test case. Its capabilities are then demonstrated with FCIQMC on an archetypical quantum chemistry problem, the chromium dimer, in an all-electron basis set with Hilbert space size of about 10 22 yielding highly accurate FCI energies.
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.
In the last decade, stochastic versions of quantum chemistry methods such as coupled cluster Monte Carlo (CCMC) or full configuration interaction quantum Monte Carlo (FCIQMC) have made highly accurate energy calculations possible that are not accessible to the corresponding deterministic methods (full configuration interaction and coupled cluster) at the same accuracy. CCMC and FCIQMC parallelize well and exploit the sparsity in the wavefunction which decreases memory costs and makes calculations in larger systems tractable. With CCMC it is straightforward to set up high order coupled cluster calculations, such as CCSDTQ5, which includes quintuple excitations explicitly. In this thesis, the convergence of the energy accuracy with the coupled cluster levels up to CCSDTQ5 was tested on the uniform electron gas, a model solid system, for various degrees of electron correlation. This gave information on what coupled cluster level is needed to reach sufficient accuracy when modelling a solid system. Before large solid systems can be modelled, the CCMC and FCIQMC algorithms need to be optimised. The efficiency in one of the crucial steps in these algorithms, the spawn step, was improved, keeping computational and memory costs as low as possible. Furthermore, the convergence of CCMC and FCIQMC was accelerated by employing a quasi-Newton propagation. Using the model system information of what coupled cluster level is needed and having made great progress towards accelerating these methods, the computation of highly accurate energies in solid or large molecular systems should be more feasible in the future.
We introduce a family of methods for the full configuration interaction problem in quantum chemistry, based on the fast randomized iteration (FRI) framework [Lim, L.-H.; Weare, J. SIAM Rev. 2017, 59, 547; DOI: 10.1137/15M1040827 ]. These methods, which we term “FCI-FRI”, stochastically impose sparsity during iterations of the power method and can be viewed as a generalization of full configuration interaction quantum Monte Carlo (FCIQMC) without walkers. In addition to the multinomial scheme commonly used to sample excitations in FCIQMC, we present a systematic scheme where excitations are not sampled independently. Performing ground-state calculations on five small molecules at fixed cost, we find that the systematic FCI-FRI scheme is 11–45 times more statistically efficient than the multinomial FCI-FRI scheme, which is in turn 1.4–178 times more statistically efficient than the original FCIQMC algorithm.
Two distinct types of Quantum Monte Carlo (QMC) calculations are applied to electronic structure problems such as calculating potential energy curves and producing benchmark values for reaction barriers. First, Variational and Diffusion Monte Carlo (VMC and DMC) methods using a trial wavefunction subject to the fixed node approximation were tested using the CASINO code.[1] Next, Full Configuration Interaction Quantum Monte Carlo (FCIQMC), along with its initiator extension (i-FCIQMC) were tested using the NECI code.[2] FCIQMC seeks the FCI energy for a specific basis set. At a reduced cost, the efficient i-FCIQMC method can be applied to systems in which the standard FCIQMC approach proves to be too costly. Since all of these methods are statistical approaches, uncertainties (error-bars) are introduced for each calculated energy. This study tests the performance of the methods relative to traditional quantum chemistry for some benchmark systems. References: [1] R. J. Needs et al., J. Phys.: Condensed Matter 22, 023201 (2010). [2] G. H. Booth et al., J. Chem. Phys. 131, 054106 (2009).
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.