1
|
Mahajan A, Kurian JS, Lee J, Reichman DR, Sharma S. Response properties in phaseless auxiliary field quantum Monte Carlo. J Chem Phys 2023; 159:184101. [PMID: 37937933 DOI: 10.1063/5.0171996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 10/11/2023] [Indexed: 11/09/2023] Open
Abstract
We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo by applying automatic differentiation (AD). Biases and statistical efficiency of the resulting estimators are discussed. Our approach demonstrates that AD enables the calculation of reduced density matrices with the same computational cost scaling per sample as energy calculations, accompanied by a cost prefactor of less than four in our numerical calculations. We investigate the role of self-consistency and trial orbital choice in property calculations. We find that orbitals obtained using density functional theory perform well for the dipole moments of selected molecules compared to those optimized self-consistently.
Collapse
Affiliation(s)
- Ankit Mahajan
- Department of Chemistry, Columbia University, New York, New York 10027, USA
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Jo S Kurian
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| |
Collapse
|
2
|
Xiao ZY, Shi H, Zhang S. Interfacing Branching Random Walks with Metropolis Sampling: Constraint Release in Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2023; 19:6782-6795. [PMID: 37661928 DOI: 10.1021/acs.jctc.3c00521] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
We present an approach to interface branching random walks with Markov chain Monte Carlo sampling and to switch seamlessly between the two. The approach is discussed in the context of auxiliary-field quantum Monte Carlo (AFQMC) but can be applied to other Monte Carlo calculations or simulations. In AFQMC, the formulation of branching random walks along imaginary-time is needed to realize a constraint to control the sign or phase problem. The constraint is derived from an exact gauge condition and is in practice implemented approximately with a trial wave function or trial density matrix, which can break exactness in the algorithm. We use the generalized Metropolis algorithm to sample a selected portion of the imaginary-time path after it has been produced by the branching random walk. This interfacing allows a constraint release to follow seamlessly from constrained-path sampling, which can reduce the systematic error from the latter. It also provides a way to improve the computation of correlation functions and observables that do not commute with the Hamiltonian. We illustrate the method in atoms and molecules, where improvements in accuracy can be clearly quantified and near-exact results are obtained. We also discuss the computation of the variance of the Hamiltonian and propose a convenient way to evaluate it stochastically without changing the scaling of AFQMC.
Collapse
Affiliation(s)
- Zhi-Yu Xiao
- Department of Physics, College of William & Mary, Williamsburg, Virginia 23187, United States
| | - Hao Shi
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
3
|
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
|
4
|
Malone FD, Mahajan A, Spencer JS, Lee J. ipie: A Python-Based Auxiliary-Field Quantum Monte Carlo Program with Flexibility and Efficiency on CPUs and GPUs. J Chem Theory Comput 2023; 19:109-121. [PMID: 36503227 DOI: 10.1021/acs.jctc.2c00934] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We report the development of a python-based auxiliary-field quantum Monte Carlo (AFQMC) program, ipie, with preliminary timing benchmarks and new AFQMC results on the isomerization of [Cu2O2]2+. We demonstrate how implementations for both central and graphical processing units (CPUs and GPUs) are achieved in ipie. We show an interface of ipie with PySCF as well as a straightforward template for adding new estimators to ipie. Our timing benchmarks against other C++ codes, QMCPACK and Dice, suggest that ipie is faster or similarly performing for all chemical systems considered on both CPUs and GPUs. Our results on [Cu2O2]2+ using selected configuration interaction trials show that it is possible to converge the ph-AFQMC isomerization energy between bis(μ-oxo) and μ-η2:η2 peroxo configurations to the exact known results for small basis sets with 105-106 determinants. We also report the isomerization energy with a quadruple-zeta basis set with an estimated error less than a kcal/mol, which involved 52 electrons and 290 orbitals with 106 determinants in the trial wave function. These results highlight the utility of ph-AFQMC and ipie for systems with modest strong correlation and large-scale dynamic correlation.
Collapse
Affiliation(s)
- Fionn D Malone
- Google Research, Venice, California 90291, United States
| | - Ankit Mahajan
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, United States
| | | | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
5
|
Greene SM, Webber RJ, Smith JET, Weare J, Berkelbach TC. Full Configuration Interaction Excited-State Energies in Large Active Spaces from Subspace Iteration with Repeated Random Sparsification. J Chem Theory Comput 2022; 18:7218-7232. [PMID: 36345915 DOI: 10.1021/acs.jctc.2c00435] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is based on an asymmetric variant of subspace iteration, avoids the use of dot products of random vectors and instead relies upon trial vectors to maintain orthogonality and estimate eigenvalues. By leveraging recent advances, we apply our method to calculate ground- and excited-state energies of challenging molecular systems in large active spaces, including the carbon dimer with 8 electrons in 108 orbitals (8e,108o), an oxo-Mn(salen) transition metal complex (28e,28o), ozone (18e,87o), and butadiene (22e,82o). In the majority of these test cases, our approach yields total excited-state energies that agree with those from state-of-the-art methods─including heat-bath CI, the density matrix renormalization group approach, and FCIQMC─to within sub-milliHartree accuracy. In all cases, estimated excitation energies agree to within about 0.1 eV.
Collapse
Affiliation(s)
- Samuel M Greene
- Department of Chemistry, Columbia University, New York, New York10027, United States
| | - Robert J Webber
- Courant Institute of Mathematical Sciences, New York University, New York, New York10012, United States
| | - James E T Smith
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York10010, United States
| | - Jonathan Weare
- Courant Institute of Mathematical Sciences, New York University, New York, New York10012, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York10027, United States.,Center for Computational Quantum Physics, Flatiron Institute, New York, New York10010, United States
| |
Collapse
|
6
|
Lee J, Pham HQ, Reichman DR. Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking. J Chem Theory Comput 2022; 18:7024-7042. [PMID: 36255074 DOI: 10.1021/acs.jctc.2c00802] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) approach from a computational quantum chemistry perspective and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQMC for general main-group chemistry applications. For systems where single determinant wave functions are qualitatively accurate, we expect the accuracy of ph-AFQMC in conjunction with a single-determinant trial wave function to be between that of coupled-cluster with singles and doubles (CCSD) and CCSD with perturbative triples (CCSD(T)). For these applications, ph-AFQMC should be a method of choice when canonical CCSD(T) is too expensive to run. For systems where multireference (MR) wave functions are needed for qualitative accuracy, ph-AFQMC is far more accurate than MR perturbation theory methods and competitive with MR configuration interaction (MRCI) methods. Due to the computational efficiency of ph-AFQMC compared to MRCI, we recommended ph-AFQMC as a method of choice for handling dynamic correlation in MR problems. We conclude with a discussion of important directions for future development of the ph-AFQMC approach.
Collapse
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Q Pham
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
7
|
Weber JL, Vuong H, Devlaminck PA, Shee J, Lee J, Reichman DR, Friesner RA. A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2022; 18:3447-3459. [PMID: 35507769 DOI: 10.1021/acs.jctc.2c00111] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium- to large-sized molecules because of its accuracy and favorable polynomial scaling with system size. Unfortunately, the memory footprints of standard energy evaluation algorithms are nontrivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low-rank structure of the Coulombic integrals have been successful but exhibit high prefactors, making their utility limited to very large systems. Here we present a complementary cubic-scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to ph-AFQMC. We show that the error due to this approximation, which we term localized-orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable and that the method is computationally favorable even for small systems. We present results demonstrating robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized-orbital-based methods, including the singlet-triplet gaps of the polyacenes benzene through pentacene, the heats of formation for a set of Platonic hydrocarbon cages, and the total energy of ferrocene, Fe(Cp)2. Finally, we reproduce our previous result for the gas-phase ionization energy of Ni(Cp)2, agreeing with full ph-AFQMC to within statistical error while using less than 1/15th of the computer time.
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
8
|
Mahajan A, Lee J, Sharma S. Selected configuration interaction wave functions in phaseless auxiliary field quantum Monte Carlo. J Chem Phys 2022; 156:174111. [DOI: 10.1063/5.0087047] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present efficient algorithms for using selected configuration interaction (sCI) trial wave functions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC). These advances, geared towards optimizing computational performance for longer CI expansions, allow us to use up to a million configurations in the trial state for ph-AFQMC. In one example, we found the cost of ph-AFQMC per sample to increase only by a factor of about 3 for a calculation with 104 configurations compared to that with a single one, demonstrating the tiny computational overhead due to a longer expansion. This favorable scaling allows us to study the systematic convergence of the phaseless bias in AFQMC calculations with an increasing number of configurations and provides a means to gauge the accuracy of ph-AFQMC with other trial states. We also show how the scalability issues of sCI trial states for large system sizes could be mitigated by restricting them to a moderately sized orbital active space and leveraging the near-cancellation of out of active space phaseless errors.
Collapse
Affiliation(s)
- Ankit Mahajan
- University of Colorado at Boulder, United States of America
| | - Joonho Lee
- Chemistry, Columbia University, United States of America
| | - Sandeep Sharma
- University of Colorado at Boulder, United States of America
| |
Collapse
|
9
|
Lee J, Malone FD, Morales MA, Reichman DR. Spectral Functions from Auxiliary-Field Quantum Monte Carlo without Analytic Continuation: The Extended Koopmans' Theorem Approach. J Chem Theory Comput 2021; 17:3372-3387. [PMID: 33983735 DOI: 10.1021/acs.jctc.1c00100] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
We explore the extended Koopmans' theorem (EKT) within the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method. The EKT allows for the direct calculation of electron addition and removal spectral functions using reduced density matrices of the N-particle system and avoids the need for analytic continuation. The lowest level of EKT with AFQMC, called EKT1-AFQMC, is benchmarked using atoms, small molecules, 14-electron and 54-electron uniform electron gas supercells, and a minimal unit cell model of diamond at the Γ-point. Via comparison with numerically exact results (when possible) and coupled-cluster methods, we find that EKT1-AFQMC can reproduce the qualitative features of spectral functions for Koopmans-like charge excitations with errors in peak locations of less than 0.25 eV in a finite basis. We also note the numerical difficulties that arise in the EKT1-AFQMC eigenvalue problem, especially when back-propagated quantities are very noisy. We show how a systematic higher-order EKT approach can correct errors in EKT1-based theories with respect to the satellite region of the spectral function. Our work will be of use for the study of low-energy charge excitations and spectral functions in correlated molecules and solids where AFQMC can be reliably performed for both energy and back propagation.
Collapse
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
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
10
|
Shi H, Zhang S. Some recent developments in auxiliary-field quantum Monte Carlo for real materials. J Chem Phys 2021; 154:024107. [DOI: 10.1063/5.0031024] [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)
- Hao Shi
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, 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
|
11
|
Tucholska AM, Moszynski R. Molecular properties from the explicitly connected expressions of the response functions within the coupled-cluster theory. ADVANCES IN QUANTUM CHEMISTRY 2021. [DOI: 10.1016/bs.aiq.2021.05.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
12
|
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
|
13
|
Lee J, Reichman DR. Stochastic resolution-of-the-identity auxiliary-field quantum Monte Carlo: Scaling reduction without overhead. J Chem Phys 2020; 153:044131. [DOI: 10.1063/5.0015077] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [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, 3000 Broadway, New York, New York 10027, USA
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| |
Collapse
|
14
|
Malone FD, Zhang S, Morales MA. Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Units. J Chem Theory Comput 2020; 16:4286-4297. [DOI: 10.1021/acs.jctc.0c00262] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [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, Livermore, California 94550, United States
| | - Shuai Zhang
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
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
|
18
|
Motta M, Shee J, Zhang S, Chan GKL. Efficient Ab Initio Auxiliary-Field Quantum Monte Carlo Calculations in Gaussian Bases via Low-Rank Tensor Decomposition. J Chem Theory Comput 2019; 15:3510-3521. [DOI: 10.1021/acs.jctc.8b00996] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Mario Motta
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - James Shee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
19
|
Landinez Borda EJ, Gomez J, Morales MA. Non-orthogonal multi-Slater determinant expansions in auxiliary field quantum Monte Carlo. J Chem Phys 2019; 150:074105. [DOI: 10.1063/1.5049143] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - John Gomez
- Applied Physics Program and Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
20
|
Malone FD, Zhang S, Morales MA. Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting. J Chem Theory Comput 2018; 15:256-264. [DOI: 10.1021/acs.jctc.8b00944] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, United States
| | - Shuai Zhang
- 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
|
21
|
Rusakov AA, Iskakov S, Tran LN, Zgid D. Self-Energy Embedding Theory (SEET) for Periodic Systems. J Chem Theory Comput 2018; 15:229-240. [DOI: 10.1021/acs.jctc.8b00927] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alexander A. Rusakov
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Sergei Iskakov
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Lan Nguyen Tran
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
- Ho Chi Minh City Institute of Physics, Vietnam Academy of Science and Technology (VAST), Ho Chi Minh City 70000, Vietnam
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
- Center for Computational Quantum Physics, The Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
22
|
Zhang S, Malone FD, Morales MA. Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide. J Chem Phys 2018; 149:164102. [DOI: 10.1063/1.5040900] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shuai Zhang
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Fionn D. Malone
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
23
|
Dash M, Moroni S, Scemama A, Filippi C. Perturbatively Selected Configuration-Interaction Wave Functions for Efficient Geometry Optimization in Quantum Monte Carlo. J Chem Theory Comput 2018; 14:4176-4182. [PMID: 29953810 PMCID: PMC6096455 DOI: 10.1021/acs.jctc.8b00393] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We
investigate the performance of a class of compact and systematically
improvable Jastrow–Slater wave functions for the efficient
and accurate computation of structural properties, where the determinantal
component is expanded with a perturbatively selected configuration
interaction scheme (CIPSI). We concurrently optimize the molecular
ground-state geometry and full wave function—Jastrow factor,
orbitals, and configuration interaction coefficients—in variational
Monte Carlo (VMC) for the prototypical case of 1,3-trans-butadiene, a small yet theoretically challenging π-conjugated
system. We find that the CIPSI selection outperforms the conventional
scheme of correlating orbitals within active spaces chosen by chemical
intuition: it gives significantly better variational and diffusion
Monte Carlo energies for all but the smallest expansions, and much
smoother convergence of the geometry with the number of determinants.
In particular, the optimal bond lengths and bond-length alternation
of butadiene are converged to better than 1 mÅ with just a few
thousand determinants, to values very close to the corresponding CCSD(T)
results. The combination of CIPSI expansion and VMC optimization represents
an affordable tool for the determination of accurate ground-state
geometries in quantum Monte Carlo.
Collapse
Affiliation(s)
- Monika Dash
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| | - Saverio Moroni
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, and SISSA Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , I-34136 Trieste , Italy
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques , Université de Toulouse , CNRS, UPS, 118 Route de Narbonne , 31077 Toulouse , France
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology , University of Twente , P.O. Box 217, 7500 AE Enschede , The Netherlands
| |
Collapse
|
24
|
Samanta PK, Blunt NS, Booth GH. Response Formalism within Full Configuration Interaction Quantum Monte Carlo: Static Properties and Electrical Response. J Chem Theory Comput 2018; 14:3532-3546. [DOI: 10.1021/acs.jctc.8b00454] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pradipta Kumar Samanta
- Institut für Theoretische Chemie, Universität Stuttgart, D-70569 Stuttgart, Germany
- Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| | - Nick S. Blunt
- University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| |
Collapse
|
25
|
Shee J, Arthur EJ, Zhang S, Reichman DR, Friesner RA. Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units. J Chem Theory Comput 2018; 14:4109-4121. [DOI: 10.1021/acs.jctc.8b00342] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrödinger
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
26
|
Ab initio computations of molecular systems by the auxiliary‐field quantum Monte Carlo method. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1364] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
27
|
Motta M, Zhang S. Communication: Calculation of interatomic forces and optimization of molecular geometry with auxiliary-field quantum Monte Carlo. J Chem Phys 2018; 148:181101. [DOI: 10.1063/1.5029508] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mario Motta
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Shiwei Zhang
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, USA
| |
Collapse
|