Rigorous wave function embedding with dynamical fluctuations

The dynamical fluctuations in approaches such as dynamical mean-field theory (DMFT) allow for the self-consistent optimization of a local fragment, hybridized with a true correlated environment. We show that these correlated environmental fluctuations can instead be efficiently captured in a wave function perspective in a computationally cheap, frequency-independent, zero-temperature approach. This allows for a systematically improvable, short-time wave function analogue to DMFT, which entails a number of computational and numerical benefits. We demonstrate this approach to solve the correlated dynamics of the paradigmatic Bethe lattice Hubbard model, as well as detailing cluster extensions in the one-dimensional Hubbard chain where we clearly show the benefits of this rapidly convergent description of correlated environmental fluctuations

Spectral functions of strongly correlated extended systems via an exact quantum embedding

Density matrix embedding theory (DMET) [Phys. Rev. Lett., 109, 186404 (2012)], introduced a new approach to quantum cluster embedding methods, whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost, but was inherently limited by the construction of a bath designed to reproduce ground state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the 1D and 2D Hubbard model, where we obtain zero temperature, thermodynamic limit spectral functions, and show the trivial extension to two-particle Green functions. This advance therefore extends the scope and applicability of DMET in condensed matter problems as a computationally tractable route to correlated spectral functions of extended systems, and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.Comment: 6 pages, 6 figure

Linear-scaling and parallelizable algorithms for stochastic quantum chemistry

For many decades, quantum chemical method development has been dominated by algorithms which involve increasingly complex series of tensor contractions over one-electron orbital spaces. Procedures for their derivation and implementation have evolved to require the minimum amount of logic and rely heavily on computationally efficient library-based matrix algebra and optimized paging schemes. In this regard, the recent development of exact stochastic quantum chemical algorithms to reduce computational scaling and memory overhead requires a contrasting algorithmic philosophy, but one which when implemented efficiently can often achieve higher accuracy/cost ratios with small random errors. Additionally, they can exploit the continuing trend for massive parallelization which hinders the progress of deterministic high-level quantum chemical algorithms. In the Quantum Monte Carlo community, stochastic algorithms are ubiquitous but the discrete Fock space of quantum chemical methods is often unfamiliar, and the methods introduce new concepts required for algorithmic efficiency. In this paper, we explore these concepts and detail an algorithm used for Full Configuration Interaction Quantum Monte Carlo (FCIQMC), which is implemented and available in MOLPRO and as a standalone code, and is designed for high-level parallelism and linear-scaling with walker number. Many of the algorithms are also in use in, or can be transferred to, other stochastic quantum chemical methods and implementations. We apply these algorithms to the strongly correlated Chromium dimer, to demonstrate their efficiency and parallelism.Comment: 16 pages, 8 figure

Non-linear biases, stochastically-sampled effective Hamiltonians and spectral functions in quantum Monte Carlo methods

In this article we study examples of systematic biases that can occur in quantum Monte Carlo methods due to the accumulation of non-linear expectation values, and approaches by which these errors can be corrected. We begin with a study of the Krylov-projected FCIQMC (KP-FCIQMC) approach, which was recently introduced to allow efficient, stochastic calculation of dynamical properties. This requires the solution of a sampled effective Hamiltonian, resulting in a non-linear operation on these stochastic variables. We investigate the probability distribution of this eigenvalue problem to study both stochastic errors and systematic biases in the approach, and demonstrate that such errors can be significantly corrected by moving to a more appropriate basis. This is lastly expanded to include consideration of the correlation function QMC approach of Ceperley and Bernu, showing how such an approach can be taken in the FCIQMC framework.Comment: 12 pages, 7 figure

The intermediate and spin-liquid phase of the half-filled honeycomb Hubbard model

We obtain the phase-diagram of the half-filled honeycomb Hubbard model with density matrix embedding theory, to address recent controversy at intermediate couplings. We use clusters from 2-12 sites and lattices at the thermodynamic limit. We identify a paramagnetic insulating state, with possible hexagonal cluster order, competitive with the antiferromagnetic phase at intermediate coupling. However, its stability is strongly cluster and lattice size dependent, explaining controver- sies in earlier work. Our results support the paramagnetic insulator as being a metastable, rather than a true, intermediate phase, in the thermodynamic limit

Explicitly correlated plane waves: Accelerating convergence in periodic wavefunction expansions

We present an investigation into the use of an explicitly correlated plane wave basis for periodic wavefunction expansions at the level of second-order M{\o}ller-Plesset perturbation theory (MP2). The convergence of the electronic correlation energy with respect to the one-electron basis set is investigated and compared to conventional MP2 theory in a finite homogeneous electron gas model. In addition to the widely used Slater-type geminal correlation factor, we also derive and investigate a novel correlation factor that we term Yukawa-Coulomb. The Yukawa-Coulomb correlation factor is motivated by analytic results for two electrons in a box and allows for a further improved convergence of the correlation energies with respect to the employed basis set. We find the combination of the infinitely delocalized plane waves and local short-ranged geminals provides a complementary, and rapidly convergent basis for the description of periodic wavefunctions. We hope that this approach will expand the scope of discrete wavefunction expansions in periodic systems.Comment: 15 pages, 13 figure

Energy-weighted density matrix embedding of open correlated chemical fragments

We present a multi-scale approach to efficiently embed an ab initio correlated chemical fragment described by its energy-weighted density matrices, and entangled with a wider mean-field many-electron system. This approach, first presented in Phys. Rev. B, 98, 235132 (2018), is here extended to account for realistic long-range interactions and broken symmetry states. The scheme allows for a systematically improvable description in the range of correlated fluctuations out of the fragment into the system, via a self-consistent optimization of a coupled auxiliary mean-field system. It is discussed that the method has rigorous limits equivalent to existing quantum embedding approaches of both dynamical mean-field theory, as well as density matrix embedding theory, to which this method is compared, and the importance of these correlated fluctuations is demonstrated. We derive a self-consistent local energy functional within the scheme, and demonstrate the approach for Hydrogen rings, where quantitative accuracy is achieved despite only a single atom being explicitly treated.Comment: 14 pages, 8 figure

Spectroscopic accuracy directly from quantum chemistry: application to ground and excited states of beryllium dimer

We combine explicit correlation via the canonical transcorrelation approach with the density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods to compute a near-exact beryllium dimer curve, {\it without} the use of composite methods. In particular, our direct density matrix renormalization group calculations produce a well-depth of $D_e$=931.2 cm$^{-1}$ which agrees very well with recent experimentally derived estimates $D_e$=929.7$\pm 2$~cm$^{-1}$ [Science, 324, 1548 (2009)] and $D_e$=934.6~cm$^{-1}$ [Science, 326, 1382 (2009)]], as well the best composite theoretical estimates, $D_e$=938$\pm 15$~cm$^{-1}$ [J. Phys. Chem. A, 111, 12822 (2007)] and $D_e$=935.1$\pm 10$~cm$^{-1}$ [Phys. Chem. Chem. Phys., 13, 20311 (2011)]. Our results suggest possible inaccuracies in the functional form of the potential used at shorter bond lengths to fit the experimental data [Science, 324, 1548 (2009)]. With the density matrix renormalization group we also compute near-exact vertical excitation energies at the equilibrium geometry. These provide non-trivial benchmarks for quantum chemical methods for excited states, and illustrate the surprisingly large error that remains for 1$^1\Sigma^-_g$ state with approximate multi-reference configuration interaction and equation-of-motion coupled cluster methods. Overall, we demonstrate that explicitly correlated density matrix renormalization group and initiator full configuration interaction quantum Monte Carlo methods allow us to fully converge to the basis set and correlation limit of the non-relativistic Schr\"odinger equation in small molecules