Quantum Monte Carlo (QMC) methods are a class of stochastic techniques that can be used to compute the properties of electronic systems accurately from first principles. This thesis is mainly concerned with the development of trial wave functions for QMC. An extension of the backflow transformation to inhomogeneous electronic systems is presented and applied to atoms, molecules and extended systems. The backflow transformation I have developed typically retrieves an additional 50% of the remaining correlation energy at the variational Monte Carlo level, and 30% at the diffusion Monte Carlo level; the number of parameters required to achieve a given fraction of the correlation energy does not appear to increase with system size. The expense incurred by the use of backflow transformations is investigated, and it is found to scale favourably with system size. Additionally, I propose a single wave function form for studying the electron-hole system which includes pairing effects and is capable of describing all of the relevant phases of this system. The effectiveness of this general wave function is demonstrated by applying it to a particular transition between two phases of the symmetric electron-hole bilayer, and it is found that using a single wave function form gives a more accurate physical description of the system than using a different wave function to describe each phase. Both of these developments are new, and they provide a powerful set of tools for designing accurate wave functions. Backflow transformations are particularly important for systems with repulsive interactions, while pairing wave functions are important for attractive interactions. It is possible to combine backflow and pairing to further increase the accuracy of the wave function. The wave function technology that I have developed should therefore be useful across a very wide range of problems.
Quantum Monte Carlo (QMC) methods are a class of stochastic techniques that can be used to compute the properties of electronic systems accurately from first principles. This thesis is mainly concerned with the development of trial wave functions for QMC. An extension of the backflow transformation to inhomogeneous electronic systems is presented and applied to atoms, molecules and extended systems. The backflow transformation I have developed typically retrieves an additional 50% of the remaining correlation energy at the variational Monte Carlo level, and 30% at the diffusion Monte Carlo level; the number of parameters required to achieve a given fraction of the correlation energy does not appear to increase with system size. The expense incurred by the use of backflow transformations is investigated, and it is found to scale favourably with system size. Additionally, I propose a single wave function form for studying the electron-hole system which includes pairing effects and is capable of describing all of the relevant phases of this system. The effectiveness of this general wave function is demonstrated by applying it to a particular transition between two phases of the symmetric electron-hole bilayer, and it is found that using a single wave function form gives a more accurate physical description of the system than using a different wave function to describe each phase. Both of these developments are new, and they provide a powerful set of tools for designing accurate wave functions. Backflow transformations are particularly important for systems with repulsive interactions, while pairing wave functions are important for attractive interactions. It is possible to combine backflow and pairing to further increase the accuracy of the wave function. The wave function technology that I have developed should therefore be useful across a very wide range of problems.
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.
Quantum Monte Carlo (QMC) is one of the most promising methods for solving quantum many-body problems. QMC combines known analytical properties of wave functions, results from other methods, and stochastic techniques, into a powerful tool for investigation of the electronic structure of real materials. Furthermore, the method has a favorable scaling ( O[N3] in the number of electrons), perfect scalability on parallel architectures, and can be applied to a wide range of systems. Our work has focused both on expanding the capabilities of QMC for calculation of quantities beyond energies, such as interatomic forces, and on application of the method to challenging problems in nanostructure materials research. For the calculation of interatomic forces we have implemented a finite difference correlated sampling method. The correlated sampling enables us to avoid the problem of statistical noise in evaluating the energy differences and allows us to obtain forces with an average 1% error. We have analyzed the effect of the method's main approximations on the estimate of the forces, both by theoretical arguments and by performing a variety of tests on small systems. Furthermore, we have used this newly developed method for estimating equilibrium geometries in a problem where other methods have failed. QMC enabled us to treat the electron correlation effects with high accuracy and to provide accurate predictions in a number of interesting applications. In particular, for carbon rings of C4N+2 stoichiometry we have resolved a long-standing problem regarding the competition between the aromaticity and second order Jahn-Teller effects. We have found that the dominant mechanism at small and intermediate sizes, which stabilizes the bond angle and bond length alternated geometries, is the secondorder Jahn-Teller effect while the aromatic isomer is always found to be a transition state. These high accuracy calculations which have involved multi-reference wave functions, extensive variational optimizations, and the use of massively parallel platforms demonstrate clearly the unique capabilities of the QMC approach.
The fixed-node diffusion quantum Monte Carlo (DMC) method is the most accurate method known for calculating the energies of large many-particle quantum systems. The key element of the method is the development of accurate trial many-body wavefunctions which control the statistical efficiency of the calculations and the accuracy obtained. Accurate wavefunctions can be obtained by building correlation effects on top of mean field descriptions such as density functional theory. The wavefunctions can be improved by introducing multi-determinants, pairing functions, and backflow transformations. The calculations are expensive, but the method scales well with system size and calculations on 1000 particles are possible. Some recent applications of the DMC method to atoms, molecules and solids will be presented.
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
Trial wave function based quantum Monte Carlo is a promising family of methods for the solution of quantum mechanical Hamiltonians. These fundamentally non-perturbative methods can be used to treat bosons or fermions in weakly or strong interacting systems and in any phase. Properties computed for particles with bosonic statistics can be converged to their exact value at relatively moderate computational cost using diffusion quantum Monte Carlo. For fermions the problem is slightly more complicated. Exact methods, such as released-node diffusion quantum Monte Carlo, exist for computing unbiased fermion properties, however their computational expense increases rapidly with the Bose-Fermi energy gap and system size. It is possible to use a fixed-node version of diffusion quantum Monte Carlo, but it introduces a bias due to the nodes of the trial wave function. In this thesis we work towards reducing trial wave function bias in fixed-node calculations and then perform several benchmark studies. We begin with a pedagogical overview of the three methods used throughout the thesis, variational, diffusion, and reptation quantum Monte Carlo. Then we discuss the problem, trial wave function bias in the energy and other observables. Next we outline the Hamiltonians we are typically interested in for electronic systems and trial wave function forms used to solve them. The algorithm used to optimize the trial wave functions is presented along with special considerations for some particular cases. Next we present three studies of the pressure of the electron gas using improved estimators in variational, diffusion, and reptation quantum Monte Carlo. Benchmark results for a set of molecules are then presented for the massive multi-determinant expansion and optimization algorithm tailored to suit it. We conclude with two studies of the interacting Fermi liquid, one on the momentum properties of the electron gas in three dimensions and another on the entanglement properties of interacting Fermi liquids in two dimensions.
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely: the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudo potential, and the basis set for QMC calculations. We also introduce a new strategy for the definition of the atomic orbitals involved in the Jastrow - Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.