Shifted contour (SC) auxiliary field Monte Carlo (AFMC) for electronic structure is described in considerable algorithmic detail. Several applications demonstrate the high accuracy achievable in molecular electronic structure, including excited states, using the techniques of correlated sampling and variational AFMC. SC-AFMC is also applicable to large Hubbard lattices of strongly correlated electrons. Applications are shown, up to lattices of size 16 ×16 with 202 electrons. Finally, we present a method for exact coarse graining of the AFMC integrand.
Shifted contour (SC) auxiliary field Monte Carlo (AFMC) for electronic structure is described in considerable algorithmic detail. Several applications demonstrate the high accuracy achievable in molecular electronic structure, including excited states, using the techniques of correlated sampling and variational AFMC. SC-AFMC is also applicable to large Hubbard lattices of strongly correlated electrons. Applications are shown, up to lattices of size 16 ×16 with 202 electrons. Finally, we present a method for exact coarse graining of the AFMC integrand.
The shifted-contour auxiliary-field Monte Carlo (SCAFMC) approach has been recently developed by Rom, Charutz and Neuhauser as an extension of the auxiliary-field Monte Carlo (AFMC) method. AFMC replaces the difficult fully interacting electrons problem by an ensemble of simpler problems where the electrons interact with a fluctuating electric field but not with each other. SCAFMC is based on shifting the auxiliary-field contour of integration to pass through the (imaginary) stationary point, leading to numerical stability at long propagation times. The new approach converges to the full CI energy in electronic structure calculations (both ground and low-lying excited states). Here we expand the application of SCAFMC from atomic to molecular problems. First, we calculate ground-state energies of a highly correlated transition-metal system (Cr2) with a moderate (12 orbitals) active space size, and demonstrate that SCAFMC is able to extract the energies accurately. In addition, we use SCAFMC to calculate a C–C bond-stretch energy in ethane with complete active spaces of up to 28 orbitals.
The auxiliary field Monte Carlo (AFMC) technique has advantages over other ab initio quantum Monte Carlo methods for fermions, as it does not seem to require approximations for alleviating the sign problem and is directly applicable to excited states. Yet, the method is severely limited by a numerical instability, a numerical sign problem, prohibiting application to realistic electronic structure systems. Recently, the shifted contour auxiliary field method (SC-AFMC) was proposed for overcoming this instability. Here we develop a theory for the AFMC stabilization, explaining the success of SC-AFMC. We show that the auxiliary fields can be shifted into the complex plane in a manner that considerably stabilizes the Monte Carlo integration using the exact one-electron density. Practical stabilization can be achieved when an approximate Hartree–Fock density is used, proposing that an overwhelming part of the sign problem is removed by taking proper account of the Fermion mean-field contribution. The theory is demonstrated by application to H2.
The auxiliary‐field quantum Monte Carlo (AFQMC) method provides a computational framework for solving the time‐independent Schrödinger equation in atoms, molecules, solids, and a variety of model systems. AFQMC has recently witnessed remarkable growth, especially as a tool for electronic structure computations in real materials. The method has demonstrated excellent accuracy across a variety of correlated electron systems. Taking the form of stochastic evolution in a manifold of nonorthogonal Slater determinants, the method resembles an ensemble of density‐functional theory (DFT) calculations in the presence of fluctuating external potentials. Its computational cost scales as a low‐power of system size, similar to the corresponding independent‐electron calculations. Highly efficient and intrinsically parallel, AFQMC is able to take full advantage of contemporary high‐performance computing platforms and numerical libraries. In this review, we provide a self‐contained introduction to the exact and constrained variants of AFQMC, with emphasis on its applications to the electronic structure of molecular systems. Representative results are presented, and theoretical foundations and implementation details of the method are discussed. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Structure and Mechanism > Computational Materials Science Computer and Information Science > Computer Algorithms and Programming
The auxiliary-field Monte Carlo (AFMC) method, also known as the shell model Monte Carlo (SMMC), enables us to calculate microscopically statistical and collective properties of nuclei in the presence of strong correlations. These calculations can be carried out in shell model spaces that are many orders of magnitude larger than spaces that can be treated with conventional diagonalization methods. A recent major development has been the extension of the AFMC approach to heavy nuclei. Such applications to heavy nuclei have been a major challenge. On the conceptual level, a crucial question is whether a truncated spherical shell model Hamiltonian can describe the proper collectivity observed in heavy nuclei and, in particular, the rotational character of strongly deformed nuclei. On the technical level, the low excitation energies make it necessary to perform calculations down to much lower temperatures. At such low temperatures the propagator becomes ill-conditioned and requires the introduction of stabilization methods.
Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem.
Abstract This paper describes a different Q11C method for solving the Schrödinger equation, one unlike variational, diffusion, and Green’s function methods. It can be applied to many-body systems to produce solutions for ground states. For systems of fermions, antisymmetrization can be imposed to overcome the problem of node location, but the instability of the method is a serious shortcoming, so that the method as described is limited to treating relatively simple systems. Nevertheless, this auxiliary field quantum 1,fonte Carlo (AFMC) method was demonstrated to be successful for test cases reported in this paper, and in time has been improved to allow treatment of realistic systemsa. The method makes use of the imaginary-time propagation of an initial wavefunction toward the ground-state function. It requires use of a trial wavefunction which is the product of single-particle orbitals for bosons or an antisymmetrized product for fermions. A Hubbard-Stratonovich transformation of the propagator changes the many-body problem to a problem of a collection of single particles interacting with a fluctuating (auxiliary) field which replaces the many-body interactions. The authors developed an integral expression for the energy which, after discretization, could be evaluated by a Metropolis Monte Carlo integration forward in time. The algorithm was tested for one-dimensional systems of particles: 6 to 20 bosons with varied interactions, 2 to 6 neutron-proton pairs, others with repulsive potentials. Importance sampling was found to be critical and computational requirements to be large. For spin systems on a lattice the AFMC method was considered feasible, but for general multidimensional systems the method appeared to be limited by its computation requirements.