Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.
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.
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.