1
|
Sukurma Z, Schlipf M, Humer M, Taheridehkordi A, Kresse G. Toward Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2024; 20:4205-4217. [PMID: 38750634 PMCID: PMC11137827 DOI: 10.1021/acs.jctc.4c00304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Revised: 04/29/2024] [Accepted: 05/06/2024] [Indexed: 05/29/2024]
Abstract
We report modifications of the ph-AFQMC algorithm that allow the use of large time steps and reliable time step extrapolation. Our modified algorithm eliminates size-consistency errors present in the standard algorithm when large time steps are employed. We investigate various methods to approximate the exponential of the one-body operator within the AFQMC framework, distinctly demonstrating the superiority of Krylov methods over the conventional Taylor expansion. We assess various propagators within AFQMC and demonstrate that the Split-2 propagator is the optimal method, exhibiting the smallest time-step errors. For the HEAT set molecules, the time-step extrapolated energies deviate on average by only 0.19 kcal/mol from the accurate small time-step energies. For small water clusters, we obtain accurate complete basis-set binding energies using time-step extrapolation with a mean absolute error of 0.07 kcal/mol compared to CCSD(T). Using large time-step ph-AFQMC for the N2 dimer, we show that accurate bond lengths can be obtained while reducing CPU time by an order of magnitude.
Collapse
Affiliation(s)
- Zoran Sukurma
- University
of Vienna, Faculty of Physics and Center for Computational Materials
Science, Kolingasse 14-16, A-1090 Vienna, Austria
- University
of Vienna, Faculty of Physics
& Vienna Doctoral School in Physics, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - Martin Schlipf
- VASP
Software GmbH, Berggasse
21/14, 1090 Vienna, Austria
| | - Moritz Humer
- University
of Vienna, Faculty of Physics and Center for Computational Materials
Science, Kolingasse 14-16, A-1090 Vienna, Austria
- University
of Vienna, Faculty of Physics
& Vienna Doctoral School in Physics, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - Amir Taheridehkordi
- University
of Vienna, Faculty of Physics and Center for Computational Materials
Science, Kolingasse 14-16, A-1090 Vienna, Austria
| | - Georg Kresse
- University
of Vienna, Faculty of Physics and Center for Computational Materials
Science, Kolingasse 14-16, A-1090 Vienna, Austria
- VASP
Software GmbH, Sensengasse
8, 1090 Vienna, Austria
| |
Collapse
|
2
|
Peng L, Zhang X, Chan GKL. Fermionic Reduced Density Low-Rank Matrix Completion, Noise Filtering, and Measurement Reduction in Quantum Simulations. J Chem Theory Comput 2023; 19:9151-9160. [PMID: 38095484 PMCID: PMC10753808 DOI: 10.1021/acs.jctc.3c00851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 11/20/2023] [Accepted: 11/27/2023] [Indexed: 12/27/2023]
Abstract
Fermionic reduced density matrices summarize the key observables in Fermionic systems. In electronic systems, the two-particle reduced density matrix (2-RDM) is sufficient to determine the energy and most physical observables of interest. Here, we consider the possibility of using matrix completion to reconstruct the two-particle reduced density matrix to chemical accuracy from partial information. We consider the case of noiseless matrix completion, where the partial information corresponds to a subset of the 2-RDM elements, as well as noisy completion, where the partial information corresponds to both a subset of elements and statistical noise in their values. Through experiments on a set of 24 molecular systems, we find that 2-RDM can be efficiently reconstructed from a reduced amount of information. In the case of noisy completion, this results in a multiple orders of magnitude reduction in the number of measurements needed to determine the 2-RDM with chemical accuracy. These techniques can be readily applied to both classical and quantum algorithms for quantum simulations.
Collapse
Affiliation(s)
- Linqing Peng
- Division of Chemistry and Chemical
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Xing Zhang
- Division of Chemistry and Chemical
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical
Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
3
|
Van Benschoten W, Petras HR, Shepherd JJ. Electronic Free Energy Surface of the Nitrogen Dimer Using First-Principles Finite Temperature Electronic Structure Methods. J Phys Chem A 2023; 127:6842-6856. [PMID: 37535315 PMCID: PMC10440793 DOI: 10.1021/acs.jpca.3c01741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Revised: 07/20/2023] [Indexed: 08/04/2023]
Abstract
We use full configuration interaction and density matrix quantum Monte Carlo methods to calculate the electronic free energy surface of the nitrogen dimer within the free-energy Born-Oppenheimer approximation. As the temperature is raised from T = 0, we find a temperature regime in which the internal energy causes bond strengthening. At these temperatures, adding in the entropy contributions is required to cause the bond to gradually weaken with increasing temperature. We predict a thermally driven dissociation for the nitrogen dimer between 22,000 to 63,200 K depending on symmetries and basis set. Inclusion of more spatial and spin symmetries reduces the temperature required. The origin of these observations is explored using the structure of the density matrix at various temperatures and bond lengths.
Collapse
Affiliation(s)
| | - Hayley R. Petras
- Department
of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
| | - James J. Shepherd
- Department
of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
4
|
Shen T, Barghathi H, Yu J, Del Maestro A, Rubenstein BM. Stable recursive auxiliary field quantum Monte Carlo algorithm in the canonical ensemble: Applications to thermometry and the Hubbard model. Phys Rev E 2023; 107:055302. [PMID: 37329093 DOI: 10.1103/physreve.107.055302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 03/07/2023] [Indexed: 06/18/2023]
Abstract
Many experimentally accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath or use projective algorithms which may suffer from nonoptimal scaling with system size or large algorithmic prefactors. In this paper, we introduce a highly stable, recursive auxiliary field quantum Monte Carlo approach that can directly simulate systems in the canonical ensemble. We apply the method to the fermion Hubbard model in one and two spatial dimensions in a regime known to exhibit a significant "sign" problem and find improved performance over existing approaches including rapid convergence to ground-state expectation values. The effects of excitations above the ground state are quantified using an estimator-agnostic approach including studying the temperature dependence of the purity and overlap fidelity of the canonical and grand canonical density matrices. As an important application, we show that thermometry approaches often exploited in ultracold atoms that employ an analysis of the velocity distribution in the grand canonical ensemble may be subject to errors leading to an underestimation of extracted temperatures with respect to the Fermi temperature.
Collapse
Affiliation(s)
- Tong Shen
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Hatem Barghathi
- Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37916, USA
| | - Jiangyong Yu
- Department of Physics, Brown University, Providence, Rhode Island 02912, USA
| | - Adrian Del Maestro
- Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37916, USA
- Min H. Kao Department of Electrical Engineering and Computer Science, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
- Department of Physics, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
5
|
Shee J, Weber JL, Reichman DR, Friesner RA, Zhang S. On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond. J Chem Phys 2023; 158:140901. [PMID: 37061483 PMCID: PMC10089686 DOI: 10.1063/5.0134009] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 - N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.
Collapse
Affiliation(s)
- James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, USA
| |
Collapse
|
6
|
Liu Y, Meitei OR, Chin ZE, Dutt A, Tao M, Chuang IL, Van Voorhis T. Bootstrap Embedding on a Quantum Computer. J Chem Theory Comput 2023; 19:2230-2247. [PMID: 37001026 DOI: 10.1021/acs.jctc.3c00012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
Abstract
We extend molecular bootstrap embedding to make it appropriate for implementation on a quantum computer. This enables solution of the electronic structure problem of a large molecule as an optimization problem for a composite Lagrangian governing fragments of the total system, in such a way that fragment solutions can harness the capabilities of quantum computers. By employing state-of-art quantum subroutines including the quantum SWAP test and quantum amplitude amplification, we show how a quadratic speedup can be obtained over the classical algorithm, in principle. Utilization of quantum computation also allows the algorithm to match─at little additional computational cost─full density matrices at fragment boundaries, instead of being limited to 1-RDMs. Current quantum computers are small, but quantum bootstrap embedding provides a potentially generalizable strategy for harnessing such small machines through quantum fragment matching.
Collapse
Affiliation(s)
- Yuan Liu
- Department of Physics, Co-Design Center for Quantum Advantage, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Oinam R. Meitei
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Zachary E. Chin
- Department of Physics, Co-Design Center for Quantum Advantage, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Arkopal Dutt
- Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Max Tao
- Department of Physics, Co-Design Center for Quantum Advantage, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Isaac L. Chuang
- Department of Physics, Co-Design Center for Quantum Advantage, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
7
|
Dornheim T, Tolias P, Moldabekov ZA, Cangi A, Vorberger J. Effective electronic forces and potentials from ab initio path integral Monte Carlo simulations. J Chem Phys 2022; 156:244113. [PMID: 35778089 DOI: 10.1063/5.0097768] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The rigorous description of correlated quantum many-body systems constitutes one of the most challenging tasks in contemporary physics and related disciplines. In this context, a particularly useful tool is the concept of effective pair potentials that take into account the effects of the complex many-body medium consistently. In this work, we present extensive, highly accurate ab initio path integral Monte Carlo (PIMC) results for the effective interaction and the effective force between two electrons in the presence of the uniform electron gas. This gives us a direct insight into finite-size effects, thereby, opening up the possibility for novel domain decompositions and methodological advances. In addition, we present unassailable numerical proof for an effective attraction between two electrons under moderate coupling conditions, without the mediation of an underlying ionic structure. Finally, we compare our exact PIMC results to effective potentials from linear-response theory, and we demonstrate their usefulness for the description of the dynamic structure factor. All PIMC results are made freely available online and can be used as a thorough benchmark for new developments and approximations.
Collapse
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany
| | - Panagiotis Tolias
- Space and Plasma Physics, Royal Institute of Technology (KTH), Stockholm SE-100 44, Sweden
| | | | - Attila Cangi
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| |
Collapse
|
8
|
Van Benschoten WZ, Shepherd JJ. Piecewise Interaction Picture Density Matrix Quantum Monte Carlo. J Chem Phys 2022; 156:184107. [DOI: 10.1063/5.0094290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC method (IP-DMQMC) which overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.
Collapse
|
9
|
Dornheim T, Vorberger J, Militzer B, Moldabekov ZA. Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization. Phys Rev E 2021; 104:055206. [PMID: 34942706 DOI: 10.1103/physreve.104.055206] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2021] [Accepted: 10/25/2021] [Indexed: 11/07/2022]
Abstract
We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.
Collapse
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany.,Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| | - Burkhard Militzer
- Department of Earth and Planetary Science, University of California, Berkeley, California 94720, USA.,Department of Astronomy, University of California, Berkeley, California 94720, USA
| | - Zhandos A Moldabekov
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany.,Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| |
Collapse
|
10
|
Petras HR, Van Benschoten WZ, Ramadugu SK, Shepherd JJ. The Sign Problem in Density Matrix Quantum Monte Carlo. J Chem Theory Comput 2021; 17:6036-6052. [PMID: 34546738 PMCID: PMC8515812 DOI: 10.1021/acs.jctc.1c00078] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Density matrix quantum Monte Carlo (DMQMC) is a recently developed method for stochastically sampling the N-particle thermal density matrix to obtain exact-on-average energies for model and ab initio systems. We report a systematic numerical study of the sign problem in DMQMC based on simulations of atomic and molecular systems. In DMQMC, the density matrix is written in an outer product basis of Slater determinants. In principle, this means that DMQMC needs to sample a space that scales in the system size, N, as O[(exp(N))2]. In practice, removing the sign problem requires a total walker population that exceeds a system-dependent critical walker population (Nc), imposing limitations on both storage and compute time. We establish that Nc for DMQMC is the square of Nc for FCIQMC. By contrast, the minimum Nc in the interaction picture modification of DMQMC (IP-DMQMC) is only linearly related to the Nc for FCIQMC. We find that this difference originates from the difference in propagation of IP-DMQMC versus canonical DMQMC: the former is asymmetric, whereas the latter is symmetric. When an asymmetric mode of propagation is used in DMQMC, there is a much greater stochastic error and is thus prohibitively expensive for DMQMC without the interaction picture adaptation. Finally, we find that the equivalence between IP-DMQMC and FCIQMC seems to extend to the initiator approximation, which is often required to study larger systems with large basis sets. This suggests that IP-DMQMC offers a way to ameliorate the cost of moving between a Slater determinant space and an outer product basis.
Collapse
Affiliation(s)
- Hayley R Petras
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242-1294, United States
| | | | - Sai Kumar Ramadugu
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242-1294, United States
| | - James J Shepherd
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242-1294, United States
| |
Collapse
|
11
|
Church MS, Rubenstein BM. Real-time dynamics of strongly correlated fermions using auxiliary field quantum Monte Carlo. J Chem Phys 2021; 154:184103. [PMID: 34241020 DOI: 10.1063/5.0049116] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| |
Collapse
|
12
|
Dornheim T, Vorberger J. Overcoming finite-size effects in electronic structure simulations at extreme conditions. J Chem Phys 2021; 154:144103. [DOI: 10.1063/5.0045634] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| |
Collapse
|
13
|
Lee J, Morales MA, Malone FD. A phaseless auxiliary-field quantum Monte Carlo perspective on the uniform electron gas at finite temperatures: Issues, observations, and benchmark study. J Chem Phys 2021; 154:064109. [DOI: 10.1063/5.0041378] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| |
Collapse
|
14
|
Abstract
We present a Perspective on what the future holds for full configuration interaction (FCI) theory, with an emphasis on conceptual rather than technical details. Upon revisiting the early history of FCI, a number of its key contemporary approximations are compared on as equal a footing as possible, using a recent blind challenge on the benzene molecule as a testbed [Eriksen et al., J. Phys. Chem. Lett., 2020 11, 8922]. In the process, we review the scope of applications for which FCI continues to prove indispensable, and the required traits in terms of robustness, efficacy, and reliability its modern approximations must satisfy are discussed. We close by conveying a number of general observations on the merits offered by the state-of-the-art alongside some of the challenges still faced to this day. While the field has altogether seen immense progress over the years-the past decade, in particular-it remains clear that our community as a whole has a substantial way to go in enhancing the overall applicability of near-exact electronic structure theory for systems of general composition and increasing size.
Collapse
Affiliation(s)
- Janus J Eriksen
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
15
|
Dornheim T, Invernizzi M, Vorberger J, Hirshberg B. Attenuating the fermion sign problem in path integral Monte Carlo simulations using the Bogoliubov inequality and thermodynamic integration. J Chem Phys 2020; 153:234104. [DOI: 10.1063/5.0030760] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany
| | - Michele Invernizzi
- Institute of Computational Sciences, Università Della Svizzera Italiana, 6900 Lugano, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, Università Della Svizzera Italiana, 6900 Lugano, Switzerland
- Department of Physics, ETH Zurich, 8092 Zurich, Switzerland
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| | - Barak Hirshberg
- Institute of Computational Sciences, Università Della Svizzera Italiana, 6900 Lugano, Switzerland
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland
| |
Collapse
|
16
|
Dornheim T, Vorberger J. Finite-size effects in the reconstruction of dynamic properties from ab initio path integral Monte Carlo simulations. Phys Rev E 2020; 102:063301. [PMID: 33466040 DOI: 10.1103/physreve.102.063301] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Accepted: 11/15/2020] [Indexed: 06/12/2023]
Abstract
We systematically investigate finite-size effects in the dynamic structure factor S(q,ω) of the uniform electron gas obtained via the analytic continuation of ab initio path integral Monte Carlo data for the imaginary-time density-density correlation function F(q,τ). Using the recent scheme by Dornheim et al. [Phys. Rev. Lett. 121, 255001 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.255001], we find that the reconstructed spectra are not afflicted with any finite-size effects for as few as N=14 electrons both at warm dense matter (WDM) conditions and at the margins of the strongly correlated electron liquid regime. Our results further corroborate the high quality of our current description of the dynamic density response of correlated electrons, which is of high importance for many applications in WDM theory and beyond.
Collapse
Affiliation(s)
- Tobias Dornheim
- Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany
| | - Jan Vorberger
- Helmholtz-Zentrum Dresden-Rossendorf (HZDR), D-01328 Dresden, Germany
| |
Collapse
|
17
|
Shen T, Liu Y, Yu Y, Rubenstein BM. Finite temperature auxiliary field quantum Monte Carlo in the canonical ensemble. J Chem Phys 2020; 153:204108. [DOI: 10.1063/5.0026606] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tong Shen
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Yuan Liu
- Center for Ultracold Atoms, Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Yang Yu
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | |
Collapse
|
18
|
Lawrence JE, Manolopoulos DE. Confirming the role of nuclear tunneling in aqueous ferrous–ferric electron transfer. J Chem Phys 2020; 153:154114. [DOI: 10.1063/5.0022678] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Joseph E. Lawrence
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David E. Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
19
|
Harsha G, Henderson TM, Scuseria GE. Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems. J Chem Phys 2020; 153:124115. [DOI: 10.1063/5.0022702] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gaurav Harsha
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| | - Thomas M. Henderson
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
| | - Gustavo E. Scuseria
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
| |
Collapse
|
20
|
Lee J, Malone FD, Reichman DR. The performance of phaseless auxiliary-field quantum Monte Carlo on the ground state electronic energy of benzene. J Chem Phys 2020; 153:126101. [DOI: 10.1063/5.0024835] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
21
|
White AF, Kin-Lic Chan G. Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems. J Chem Phys 2020; 152:224104. [DOI: 10.1063/5.0009845] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Affiliation(s)
- Alec F. White
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
22
|
Liu Y, Shen T, Zhang H, Rubenstein B. Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo. J Chem Theory Comput 2020; 16:4298-4314. [DOI: 10.1021/acs.jctc.0c00288] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yuan Liu
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Tong Shen
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Hang Zhang
- Department of Chemical Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Brenda Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
23
|
Kent PRC, Annaberdiyev A, Benali A, Bennett MC, Landinez Borda EJ, Doak P, Hao H, Jordan KD, Krogel JT, Kylänpää I, Lee J, Luo Y, Malone FD, Melton CA, Mitas L, Morales MA, Neuscamman E, Reboredo FA, Rubenstein B, Saritas K, Upadhyay S, Wang G, Zhang S, Zhao L. QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo. J Chem Phys 2020; 152:174105. [DOI: 10.1063/5.0004860] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- P. R. C. Kent
- Center for Nanophase Materials Sciences Division and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Abdulgani Annaberdiyev
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Anouar Benali
- Computational Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, USA
| | - M. Chandler Bennett
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Edgar Josué Landinez Borda
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Peter Doak
- Center for Nanophase Materials Sciences Division and Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Hongxia Hao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Ilkka Kylänpää
- Computational Physics Laboratory, Tampere University, P.O. Box 692, 33014 Tampere, Finland
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Ye Luo
- Computational Science Division, Argonne National Laboratory, 9700 S. Cass Avenue, Lemont, Illinois 60439, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Cody A. Melton
- Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Lubos Mitas
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Fernando A. Reboredo
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Brenda Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Kayahan Saritas
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA
| | - Shiv Upadhyay
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Guangming Wang
- Department of Physics, North Carolina State University, Raleigh, North Carolina 27695-8202, USA
| | - Shuai Zhang
- Laboratory for Laser Energetics, University of Rochester, 250 E River Rd., Rochester, New York 14623, USA
| | - Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
24
|
Lee J, Malone FD, Morales MA. Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C36 Fullerene and an Iron Porphyrin Model Complex. J Chem Theory Comput 2020; 16:3019-3027. [DOI: 10.1021/acs.jctc.0c00055] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, United States
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, United States
| |
Collapse
|
25
|
Petras HR, Ramadugu SK, Malone FD, Shepherd JJ. Using Density Matrix Quantum Monte Carlo for Calculating Exact-on-Average Energies for ab Initio Hamiltonians in a Finite Basis Set. J Chem Theory Comput 2020; 16:1029-1038. [DOI: 10.1021/acs.jctc.9b01080] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | | | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, United States
| | | |
Collapse
|
26
|
Dornheim T. Fermion sign problem in path integral Monte Carlo simulations: Quantum dots, ultracold atoms, and warm dense matter. Phys Rev E 2019; 100:023307. [PMID: 31574603 DOI: 10.1103/physreve.100.023307] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2019] [Indexed: 11/07/2022]
Abstract
The ab initio thermodynamic simulation of correlated Fermi systems is of central importance for many applications, such as warm dense matter, electrons in quantum dots, and ultracold atoms. Unfortunately, path integral Monte Carlo (PIMC) simulations of fermions are severely restricted by the notorious fermion sign problem (FSP). In this paper, we present a hands-on discussion of the FSP and investigate in detail its manifestation with respect to temperature, system size, interaction-strength and -type, and the dimensionality of the system. Moreover, we analyze the probability distribution of fermionic expectation values, which can be non-Gaussian and fat-tailed when the FSP is severe. As a practical application, we consider electrons and dipolar atoms in a harmonic confinement, and the uniform electron gas in the warm dense matter regime. In addition, we provide extensive PIMC data, which can be used as a reference for the development of new methods and as a benchmark for approximations.
Collapse
Affiliation(s)
- T Dornheim
- Center for Advanced Systems Understanding (CASUS), Görlitz, Germany
| |
Collapse
|
27
|
He YY, Shi H, Zhang S. Reaching the Continuum Limit in Finite-Temperature Ab Initio Field-Theory Computations in Many-Fermion Systems. PHYSICAL REVIEW LETTERS 2019; 123:136402. [PMID: 31697528 DOI: 10.1103/physrevlett.123.136402] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Indexed: 06/10/2023]
Abstract
Finite-temperature grand-canonical computations based on field theory are widely applied in areas including condensed matter physics, ultracold atomic gas systems, and the lattice gauge theory. However, these calculations have computational costs scaling as N_{s}^{3} with the size of the lattice or basis set, N_{s}. We report a new approach based on systematically controllable low-rank factorization that reduces the scaling of such computations to N_{s}N_{e}^{2}, where N_{e} is the average number of fermions in the system. In any realistic calculations aiming to describe the continuum limit, N_{s}/N_{e} is large and needs to be extrapolated effectively to infinity for convergence. The method thus fundamentally changes the prospect for finite-temperature many-body computations in correlated fermion systems. Its application, in combination with frameworks to control the sign or phase problem as needed, will provide a powerful tool in ab initio quantum chemistry and correlated electron materials. We demonstrate the method by computing exact properties of the two-dimensional Fermi gas with zero-range attractive interaction as a function of temperature in both the normal and superfluid states.
Collapse
Affiliation(s)
- Yuan-Yao He
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, USA
| | - Hao Shi
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, USA
| |
Collapse
|
28
|
Lee J, Malone FD, Morales MA. An auxiliary-Field quantum Monte Carlo perspective on the ground state of the dense uniform electron gas: An investigation with Hartree-Fock trial wavefunctions. J Chem Phys 2019. [DOI: 10.1063/1.5109572] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- College of Chemistry, University of California, Berkeley, California 94720, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| |
Collapse
|
29
|
Dornheim T, Groth S, Filinov AV, Bonitz M. Path integral Monte Carlo simulation of degenerate electrons: Permutation-cycle properties. J Chem Phys 2019; 151:014108. [DOI: 10.1063/1.5093171] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- T. Dornheim
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstr. 15, Kiel, Germany
- Center for Advanced Systems Understanding (CASUS), Görlitz, Germany
| | - S. Groth
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstr. 15, Kiel, Germany
| | - A. V. Filinov
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstr. 15, Kiel, Germany
- Joint Institute for High Temperatures RAS, Izhorskaya Str. 13, Moscow, Russia
| | - M. Bonitz
- Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, Leibnizstr. 15, Kiel, Germany
| |
Collapse
|
30
|
White AF, Chan GKL. A Time-Dependent Formulation of Coupled-Cluster Theory for Many-Fermion Systems at Finite Temperature. J Chem Theory Comput 2018; 14:5690-5700. [DOI: 10.1021/acs.jctc.8b00773] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alec F. White
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|