1
|
Kats D, Christlmaier EMC, Schraivogel T, Alavi A. Orbital optimisation in xTC transcorrelated methods. Faraday Discuss 2024. [PMID: 39072553 DOI: 10.1039/d4fd00036f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Additionally, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The accuracy of these methods in comparison to standard xTC methods is demonstrated, and the advantages and disadvantages of the orbital optimisation are discussed.
Collapse
Affiliation(s)
- Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| | | | - Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany.
| |
Collapse
|
2
|
Liao K, Ding L, Schilling C. Quantum Information Orbitals (QIO): Unveiling Intrinsic Many-Body Complexity by Compressing Single-Body Triviality. J Phys Chem Lett 2024; 15:6782-6790. [PMID: 38913404 DOI: 10.1021/acs.jpclett.4c01105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
The simultaneous treatment of static and dynamic correlations in strongly correlated electron systems is a critical challenge. In particular, finding a universal scheme for identifying a single-particle orbital basis that minimizes the representational complexity of the many-body wave function is a formidable and longstanding problem. As a contribution toward its solution, we show that the total orbital correlation actually reveals and quantifies the intrinsic complexity of the wave function, once it is minimized via orbital rotations. To demonstrate the power of this concept in practice, an iterative scheme is proposed to optimize the orbitals by minimizing the total orbital correlation calculated by the tailored coupled cluster singles and doubles (TCCSD) ansatz. The optimized orbitals enable the limited TCCSD ansatz to capture more nontrivial information on the many-body wave function, indicated by the improved wave function and energy. An initial application of this scheme shows great improvement of TCCSD in predicting the singlet ground state potential energy curves of the strongly correlated C2 and Cr2 molecule.
Collapse
Affiliation(s)
- Ke Liao
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Lexin Ding
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Christian Schilling
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| |
Collapse
|
3
|
Dobrautz W, Sokolov IO, Liao K, Ríos PL, Rahm M, Alavi A, Tavernelli I. Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method. J Chem Theory Comput 2024; 20:4146-4160. [PMID: 38723159 PMCID: PMC11137825 DOI: 10.1021/acs.jctc.4c00070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 04/15/2024] [Accepted: 04/16/2024] [Indexed: 05/29/2024]
Abstract
Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.
Collapse
Affiliation(s)
- Werner Dobrautz
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, 41296 Gothenburg, Sweden
| | - Igor O. Sokolov
- IBM
Quantum, IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| | - Ke Liao
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Pablo López Ríos
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Martin Rahm
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, 41296 Gothenburg, Sweden
| | - Ali Alavi
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, U.K.
| | - Ivano Tavernelli
- IBM
Quantum, IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| |
Collapse
|
4
|
Urban L, Laqua H, Thompson TH, Ochsenfeld C. Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term. J Chem Theory Comput 2024; 20:3706-3718. [PMID: 38626443 PMCID: PMC11099969 DOI: 10.1021/acs.jctc.4c00193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Revised: 03/23/2024] [Accepted: 03/25/2024] [Indexed: 04/18/2024]
Abstract
We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 μEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.
Collapse
Affiliation(s)
- Lars Urban
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
- Max
Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| | - Henryk Laqua
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Travis H. Thompson
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
- Max
Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| |
Collapse
|
5
|
Safari AA, Anderson RJ, Manni GL. Toward a Stochastic Complete Active Space Second-Order Perturbation Theory. J Phys Chem A 2024; 128:281-291. [PMID: 38154124 PMCID: PMC10788896 DOI: 10.1021/acs.jpca.3c05109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 11/14/2023] [Accepted: 11/14/2023] [Indexed: 12/30/2023]
Abstract
In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.
Collapse
Affiliation(s)
- Arta A. Safari
- Max-Planck-Institute for Solid State
Research, 70569 Stuttgart, Germany
| | | | - Giovanni Li Manni
- Max-Planck-Institute for Solid State
Research, 70569 Stuttgart, Germany
| |
Collapse
|
6
|
Hou L, Irons TJP, Wang Y, Furness JW, Wibowo-Teale AM, Sun J. Capturing the electron-electron cusp with the coupling-constant averaged exchange-correlation hole: A case study for Hooke's atoms. J Chem Phys 2024; 160:014103. [PMID: 38180252 DOI: 10.1063/5.0173370] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 11/15/2023] [Indexed: 01/06/2024] Open
Abstract
In density-functional theory, the exchange-correlation (XC) energy can be defined exactly through the coupling-constant (λ) averaged XC hole n̄xc(r,r'), representing the probability depletion of finding an electron at r' due to an electron at r. Accurate knowledge of n̄xc(r,r') has been crucial for developing XC energy density-functional approximations and understanding their performance for molecules and materials. However, there are very few systems for which accurate XC holes have been calculated since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of λ while the electron density remains fixed at the physical (λ = 1) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron-electron cusp using finite basis sets. Focusing on Hooke's atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the λ-integrated XC hole. Our results indicate that, for Hooke's atoms, the error introduced by the description of the electron-electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle-, and coupling-constant-averaged XC holes are also calculated and provide a benchmark against which the Perdew-Burke-Ernzerhof and local density approximation XC hole models are assessed.
Collapse
Affiliation(s)
- Lin Hou
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Tom J P Irons
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Yanyong Wang
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - James W Furness
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| | - Andrew M Wibowo-Teale
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, N-0315 Oslo, Norway
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, USA
| |
Collapse
|
7
|
Neufeld VA, Berkelbach TC. Highly Accurate Electronic Structure of Metallic Solids from Coupled-Cluster Theory with Nonperturbative Triple Excitations. PHYSICAL REVIEW LETTERS 2023; 131:186402. [PMID: 37977636 DOI: 10.1103/physrevlett.131.186402] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Accepted: 09/28/2023] [Indexed: 11/19/2023]
Abstract
Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T))-often considered the "gold standard" of main-group quantum chemistry-is inapplicable to three-dimensional metals due to an infrared divergence, preventing its application to many important problems in materials science. We study the full, nonperturbative inclusion of triple excitations (CCSDT) and propose a new, iterative method, which we call ring-CCSDT, that resums the essential triple excitations with the same N^{7} run-time scaling as CCSD(T). CCSDT and ring-CCSDT are used to calculate the correlation energy of the uniform electron gas at metallic densities and the structural properties of solid lithium. Inclusion of connected triple excitations is shown to be essential to achieving high accuracy. We also investigate semiempirical CC methods based on spin-component scaling and the distinguishable cluster approximation and find that they enhance the accuracy of their parent ab initio methods.
Collapse
Affiliation(s)
- Verena A Neufeld
- 1Department of Chemistry, Columbia University, New York, New York 10027, USA
- 2Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | | |
Collapse
|
8
|
Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J Chem Theory Comput 2023; 19:6933-6991. [PMID: 37216210 PMCID: PMC10601490 DOI: 10.1021/acs.jctc.3c00182] [Citation(s) in RCA: 56] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Indexed: 05/24/2023]
Abstract
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ignacio Fdez. Galván
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Ali Alavi
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Yusuf Hamied
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Flavia Aleotti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Francesco Aquilante
- Theory and
Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
| | - Davide Avagliano
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Alberto Baiardi
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jie J. Bao
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Stefano Battaglia
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Letitia Birnoschi
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Alejandro Blanco-González
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Sergey I. Bokarev
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
- Chemistry
Department, School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Ria Broer
- Theoretical
Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Roberto Cacciari
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Paul B. Calio
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Rebecca K. Carlson
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Rafael Carvalho Couto
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luis Cerdán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
- Instituto
de Óptica (IO−CSIC), Consejo
Superior de Investigaciones Científicas, 28006, Madrid, Spain
| | - Liviu F. Chibotaru
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Nicholas F. Chilton
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | | | - Irene Conti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Juliana Cuéllar-Zuquin
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Razan E. Daoud
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Nike Dattani
- HPQC Labs, Waterloo, N2T 2K9 Ontario Canada
- HPQC College, Waterloo, N2T 2K9 Ontario Canada
| | - Piero Decleva
- Istituto
Officina dei Materiali IOM-CNR and Dipartimento di Scienze Chimiche
e Farmaceutiche, Università degli
Studi di Trieste, I-34121 Trieste, Italy
| | - Coen de Graaf
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
- ICREA, Pg. Lluís
Companys 23, 08010 Barcelona, Spain
| | - Mickaël
G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luca De Vico
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Werner Dobrautz
- Chalmers
University of Technology, Department of Chemistry
and Chemical Engineering, 41296 Gothenburg, Sweden
| | - Sijia S. Dong
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Chemical Biology, Department of Physics, and Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Rulin Feng
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, Fudan University, Shanghai 200433, China
| | - Nicolas Ferré
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | | | - Laura Gagliardi
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Garavelli
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Yafu Guan
- State Key
Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Meiyuan Guo
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Matthew R. Hennefarth
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chad E. Hoyer
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Miquel Huix-Rotllant
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | - Vishal Kumar Jaiswal
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Andy Kaiser
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Danil S. Kaliakin
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Marjan Khamesian
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Daniel S. King
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Vladislav Kochetov
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Marek Krośnicki
- Institute
of Theoretical Physics and Astrophysics, Faculty of Mathematics, Physics
and Informatics, University of Gdańsk, ul Wita Stwosza 57, 80-952, Gdańsk, Poland
| | | | - Ernst D. Larsson
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Susi Lehtola
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 University of Helsinki, Finland
| | - Marie-Bernadette Lepetit
- Condensed
Matter Theory Group, Institut Néel, CNRS UPR 2940, 38042 Grenoble, France
- Theory
Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409-1061, United States
| | - Pablo López Ríos
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department
of Chemistry − Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| | - Dongxia Ma
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | | | - Francesco Montorsi
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Maximilian Mörchen
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Artur Nenov
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Vu Ha Anh Nguyen
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Yoshio Nishimoto
- Graduate
School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Meagan S. Oakley
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Markus Oppel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Riddhish Pandharkar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Felix Plasser
- Department
of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K.
| | - Gerardo Raggi
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Quantum
Materials and Software LTD, 128 City Road, London, EC1V 2NX, United Kingdom
| | - Elisa Rebolini
- Scientific
Computing Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Markus Reiher
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ivan Rivalta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Thies Romig
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Arta Anushirwan Safari
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aitor Sánchez-Mansilla
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
| | - Andrew M. Sand
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Igor Schapiro
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Thais R. Scott
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Javier Segarra-Martí
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Francesco Segatta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Dumitru-Claudiu Sergentu
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Laboratory
RA-03, RECENT AIR, A. I. Cuza University of Iaşi, RA-03 Laboratory (RECENT AIR), Iaşi 700506, Romania
| | - Prachi Sharma
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, USA
| | - Yinan Shu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Jakob K. Staab
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | | | - Bruno Nunes Cabral Tenorio
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Liviu Ungur
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Nantes
Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Valera Veryazov
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Torben Arne Voß
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oskar Weser
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Dihua Wu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Xuchun Yang
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - David Yarkony
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chen Zhou
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Roland Lindh
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala
Center for Computational Chemistry (UC3), Uppsala University, PO Box 576, SE-751 23 Uppsala. Sweden
| |
Collapse
|
9
|
Ammar A, Scemama A, Giner E. Transcorrelated selected configuration interaction in a bi-orthonormal basis and with a cheap three-body correlation factor. J Chem Phys 2023; 159:114121. [PMID: 37732558 DOI: 10.1063/5.0163831] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Accepted: 08/21/2023] [Indexed: 09/22/2023] Open
Abstract
In this work, we develop a mathematical framework for a selected configuration interaction (SCI) algorithm within a bi-orthogonal basis for transcorrelated (TC) calculations. The bi-orthogonal basis used here serves as the equivalent of the standard Hartree-Fock (HF) orbitals. However, within the context of TC, it leads to distinct orbitals for the left and right vectors. Our findings indicate that the use of such a bi-orthogonal basis allows for a proper definition of the frozen core approximation. In contrast, the use of HF orbitals results in bad error cancellations for ionization potentials and atomization energies (AE). Compared to HF orbitals, the optimized bi-orthogonal basis significantly reduces the positive part of the second-order energy (PT2), thereby facilitating the use of standard extrapolation techniques of hermitian SCI. While we did not observe a significant improvement in the convergence of the SCI algorithm, this is largely due to the use in this work of a simple three-body correlation factor introduced in a recent study. This correlation factor, which depends only on atomic parameters, eliminates the need for re-optimization of the correlation factor for molecular systems, making its use straightforward and user-friendly. Despite the simplicity of this correlation factor, we were able to achieve accurate results on the AE of a series of 14 molecules on a triple-zeta basis. We also successfully broke a double bond until the full dissociation limit while maintaining the size consistency property. This work thus demonstrates the potential of the BiO-TC-SCI approach in handling complex molecular systems.
Collapse
Affiliation(s)
- Abdallah Ammar
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
10
|
Lee N, Thom AJW. Studies on the Transcorrelated Method. J Chem Theory Comput 2023; 19:5743-5759. [PMID: 37640393 PMCID: PMC10500994 DOI: 10.1021/acs.jctc.3c00046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Indexed: 08/31/2023]
Abstract
We investigate the possibility of using a transcorrelated (TC) Hamiltonian to describe electron correlation. A method to obtain TC wavefunctions was developed based on the mathematical framework of the bi-variational principle. This involves the construction of an effective TC Hamiltonian matrix, which can be solved in a self-consistent manner. This was optimized using a method we call second-order-moment minimization and demonstrate that it is possible to obtain highly accurate energies for some closed-shell atoms and helium-like ions. The effects of certain correlator terms on the description of electron-electron and electron-nuclear cusps were also examined graphically, and some TC wavefunctions were compared against near-exact Hylleraas wavefunctions.
Collapse
Affiliation(s)
- Nicholas Lee
- Department
of Chemistry, Physical and Theoretical Chemistry
Laboratory, South Parks
Road, Oxford OX1 3QZ, U.K.
| | - Alex J. W. Thom
- Yusuf
Hamied Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, U.K.
| |
Collapse
|
11
|
Ammar A, Scemama A, Giner E. Biorthonormal Orbital Optimization with a Cheap Core-Electron-Free Three-Body Correlation Factor for Quantum Monte Carlo and Transcorrelation. J Chem Theory Comput 2023. [PMID: 37390472 DOI: 10.1021/acs.jctc.3c00257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/02/2023]
Abstract
We introduce a novel three-body correlation factor that is designed to vanish in the core region around each nucleus and approach a universal two-body correlation factor for valence electrons. The transcorrelated Hamiltonian is used to optimize the orbitals of a single Slater determinant within a biorthonormal framework. The Slater-Jastrow wave function is optimized on a set of atomic and molecular systems containing both second-row elements and 3d transition metal elements. The optimization of the correlation factor and the orbitals, along with an increase in the basis set, results in a systematic lowering of the variational Monte Carlo energy for all systems tested. Importantly, the optimal parameters of the correlation factor obtained for atomic systems can be transferred to molecules. Additionally, the present correlation factor is computationally efficient and uses a mixed analytical-numerical integration scheme that reduces the costly numerical integration from R6 to R3.
Collapse
Affiliation(s)
- Abdallah Ammar
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique, Sorbonne Université and CNRS, F-75005 Paris, France
| |
Collapse
|
12
|
Haupt JP, Hosseini SM, López Ríos P, Dobrautz W, Cohen A, Alavi A. Optimizing Jastrow factors for the transcorrelated method. J Chem Phys 2023; 158:2895246. [PMID: 37290083 DOI: 10.1063/5.0147877] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Accepted: 05/25/2023] [Indexed: 06/10/2023] Open
Abstract
We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods, such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than those obtained by minimizing the variational energy. We compute all-electron atomization energies for the challenging first-row molecules C2, CN, N2, and O2 and find that the TC method yields chemically accurate results using only the cc-pVTZ basis set, roughly matching the accuracy of non-TC calculations with the much larger cc-pV5Z basis set. We also investigate an approximation in which pure three-body excitations are neglected from the TC-FCIQMC dynamics, saving storage and computational costs, and show that it affects relative energies negligibly. Our results demonstrate that the combination of tailored real-space Jastrow factors with the multi-configurational TC-FCIQMC method provides a route to obtaining chemical accuracy using modest basis sets, obviating the need for basis-set extrapolation and composite techniques.
Collapse
Affiliation(s)
- J Philip Haupt
- Max-Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | | | - Pablo López Ríos
- Max-Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Werner Dobrautz
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 41296 Gothenburg, Sweden
| | - Aron Cohen
- DeepMind, 6 Pancras Square, London N1C 4AG, United Kingdom
| | - Ali Alavi
- Max-Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
13
|
Ludovicy J, Dahl R, Lüchow A. Toward Compact Selected Configuration Interaction Wave Functions with Quantum Monte Carlo─A Case Study of C 2. J Chem Theory Comput 2023; 19:2792-2803. [PMID: 37130194 DOI: 10.1021/acs.jctc.2c01229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
The 1Σg+ ground state of C2 is investigated using truncated CIPSI-Jastrow CSF wave functions with Hartree-Fock orbitals within the framework of variational and diffusion quantum Monte Carlo. The truncation is performed based on the absolute value of the CI coefficients, and the Jastrow, molecular orbitals, and CI parameters are either partially or fully reoptimized with respect to the variational energy. Excellent absolute as well as bond dissociation energies are obtained at DMC level with very compact, fully optimized wave functions. By studying the expansions in more detail, we observe a change in the CI picture when reoptimizing the antisymmetric part of the CIPSI-Jastrow wave functions. Furthermore, we demonstrate that a decrease in the VMC energy as well as an improvement of the nodal surface quality can be achieved─with the same expansion size─if the CSFs are selected in the presence of a Jastrow correlation function, laying the foundation for a Jastrow selected CI scheme with quantum Monte Carlo.
Collapse
Affiliation(s)
- Jil Ludovicy
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52074 Aachen, Germany
| | - Robin Dahl
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52074 Aachen, Germany
| | - Arne Lüchow
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52074 Aachen, Germany
| |
Collapse
|
14
|
Li X, Li Z, Chen J. Ab initio calculation of real solids via neural network ansatz. Nat Commun 2022; 13:7895. [PMID: 36550157 PMCID: PMC9780243 DOI: 10.1038/s41467-022-35627-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Accepted: 12/13/2022] [Indexed: 12/24/2022] Open
Abstract
Neural networks have been applied to tackle many-body electron correlations for small molecules and physical models in recent years. Here we propose an architecture that extends molecular neural networks with the inclusion of periodic boundary conditions to enable ab initio calculation of real solids. The accuracy of our approach is demonstrated in four different types of systems, namely the one-dimensional periodic hydrogen chain, the two-dimensional graphene, the three-dimensional lithium hydride crystal, and the homogeneous electron gas, where the obtained results, e.g. total energies, dissociation curves, and cohesive energies, reach a competitive level with many traditional ab initio methods. Moreover, electron densities of typical systems are also calculated to provide physical intuition of various solids. Our method of extending a molecular neural network to periodic systems can be easily integrated into other neural network structures, highlighting a promising future of ab initio solution of more complex solid systems using neural network ansatz, and more generally endorsing the application of machine learning in materials simulation and condensed matter physics.
Collapse
Affiliation(s)
- Xiang Li
- ByteDance Inc, Zhonghang Plaza, No. 43, North 3rd Ring West Road, Haidian District, Beijing, China
| | - Zhe Li
- ByteDance Inc, Zhonghang Plaza, No. 43, North 3rd Ring West Road, Haidian District, Beijing, China
| | - Ji Chen
- grid.11135.370000 0001 2256 9319School of Physics, Interdisciplinary Institute of Light-Element Quantum Materials, Frontiers Science Center for Nano-Optoelectronics, Peking University, Beijing, 100871 P. R. China
| |
Collapse
|
15
|
Jiang T, Fang W, Alavi A, Chen J. General Analytical Nuclear Forces and Molecular Potential Energy Surface from Full Configuration Interaction Quantum Monte Carlo. J Chem Theory Comput 2022; 18:7233-7242. [PMID: 36326847 DOI: 10.1021/acs.jctc.2c00440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe dynamics at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy surfaces (PES) from FCIQMC. In this work, we derive the general formula for nuclear forces from FCIQMC, and clarify different contributions of the total force. This method to obtain FCIQMC forces eliminates previous restrictions and can be used with frozen core approximation and free selection of orbitals, making it promising for more efficient nuclear forces calculations. After some numerical checks of this procedure on the binding curve of N2 molecule, we use the FCIQMC energy and force to obtain the full-dimensional ground state PES of the water molecule via Gaussian processes regression. The new water FCIQMC PES can be used as the basis for H2O ground state nuclear dynamics, structure optimization, and rotation-vibrational spectrum calculation.
Collapse
Affiliation(s)
- Tonghuan Jiang
- School of Physics, Peking University, Beijing100871, P. R. China
| | - Wei Fang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian116023, P. R. China.,Department of Chemistry, Fudan University, Shanghai200438, P. R. China
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569Stuttgart, Germany.,University of Cambridge, Lensfield Road, CambridgeCB2 1EW, United Kingdom
| | - Ji Chen
- School of Physics, Peking University, Beijing100871, P. R. China.,Collaborative Innovation Center of Quantum Matter, Beijing100871, P. R. China.,Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing100871, P. R. China.,Frontiers Science Center for Nano-Optoelectronics, Peking University, Beijing100871, P. R. China
| |
Collapse
|
16
|
Ammar A, Giner E, Scemama A. Optimization of Large Determinant Expansions in Quantum Monte Carlo. J Chem Theory Comput 2022; 18:5325-5336. [PMID: 35997484 DOI: 10.1021/acs.jctc.2c00556] [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/30/2022]
Abstract
We present a new method for the optimization of large configuration interaction (CI) expansions in the quantum Monte Carlo (QMC) framework. The central idea here is to replace the nonorthogonal variational optimization of CI coefficients performed in usual QMC calculations by an orthogonal non-Hermitian optimization thanks to the so-called transcorrelated (TC) framework, the two methods yielding the same results in the limit of a complete basis set. By rewriting the TC equations as an effective self-consistent Hermitian problem, our approach requires the sampling of a single quantity per Slater determinant, leading to minimal memory requirements in the QMC code. Using analytical quantities obtained from both the TC framework and the usual CI-type calculations, we also propose improved estimators which reduce the statistical fluctuations of the sampled quantities by more than an order of magnitude. We demonstrate the efficiency of this method on wave functions containing 105-106 Slater determinants, using effective core potentials or all-electron calculations. In all the cases, a sub-milli-Hartree convergence is reached within only two or three iterations of optimization.
Collapse
Affiliation(s)
- Abdallah Ammar
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse - CNRS, 118 route de Narbonne, 31062 Toulouse cedex 09, France
| | - Emmanuel Giner
- Laboratoire de Chimie Théorique (UMR 7616), Université Paris Sorbonne - CNRS, 4 place Jussieu, 75052 Paris cedex, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse - CNRS, 118 route de Narbonne, 31062 Toulouse cedex 09, France
| |
Collapse
|
17
|
Luo H, Alavi A. Perturbation Calculation of the Uniform Electron Gas with a Transcorrelated Hamiltonian. J Chem Phys 2022; 157:074105. [DOI: 10.1063/5.0101776] [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
With a transcorrelated Hamiltonian, we perform a many body perturbation (MBPT) calculation on the uniform electron gas in the high density regime. By using a correlation factor optimised for a single determinant Jastrow ansatz, the second order correlation energy is calculated as $\frac{1-\ln2}{\pi^{2}}\ln(r_{s})-0.05075$. This already reproduces the exact logarithmic term of the random phase approximation (RPA) result, while the constant term is roughly $7\%$ larger than the RPA one. The close agreement with the RPA method demonstrates that the transcorrelated method offers a viable and potentially efficient method for treating metallic systems.
Collapse
Affiliation(s)
- Hongjun Luo
- Max-Planck-Institute for Solid State Research, Germany
| | - Ali Alavi
- Max-Planck-Institute for Solid State Research, Germany
| |
Collapse
|
18
|
Baiardi A, Lesiuk M, Reiher M. Explicitly Correlated Electronic Structure Calculations with Transcorrelated Matrix Product Operators. J Chem Theory Comput 2022; 18:4203-4217. [PMID: 35666238 DOI: 10.1021/acs.jctc.2c00167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we present the first implementation of the transcorrelated electronic Hamiltonian in an optimization procedure for matrix product states by the density matrix renormalization group (DMRG) algorithm. In the transcorrelation ansatz, the electronic Hamiltonian is similarity-transformed with a Jastrow factor to describe the cusp in the wave function at electron-electron coalescence. As a result, the wave function is easier to approximate accurately with the conventional expansion in terms of one-particle basis functions and Slater determinants. The transcorrelated Hamiltonian in first quantization comprises up to three-body interactions, which we deal with in the standard way by applying robust density fitting to two- and three-body integrals entering the second-quantized representation of this Hamiltonian. The lack of hermiticity of the transcorrelated Hamiltonian is taken care of along the lines of the first work on transcorrelated DMRG [ J. Chem. Phys. 2020, 153, 164115] by encoding it as a matrix product operator and optimizing the corresponding ground state wave function with imaginary-time time-dependent DMRG. We demonstrate our quantum chemical transcorrelated DMRG approach at the example of several atoms and first-row diatomic molecules. We show that transcorrelation improves the convergence rate to the complete basis set limit in comparison to conventional DMRG. Moreover, we study extensions of our approach that aim at reducing the cost of handling the matrix product operator representation of the transcorrelated Hamiltonian.
Collapse
Affiliation(s)
- Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Michał Lesiuk
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland.,Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
19
|
Weser O, Liebermann N, Kats D, Alavi A, Li Manni G. Spin Purification in Full-CI Quantum Monte Carlo via a First-Order Penalty Approach. J Phys Chem A 2022; 126:2050-2060. [PMID: 35298155 PMCID: PMC8978180 DOI: 10.1021/acs.jpca.2c01338] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
![]()
In this article,
we demonstrate that a first-order spin penalty
scheme can be efficiently applied to the Slater determinant based
Full-CI Quantum Monte Carlo (FCIQMC) algorithm, as a practical route
toward spin purification. Two crucial applications are presented to
demonstrate the validity and robustness of this scheme: the 1Δg ← 3Σg vertical excitation in O2 and
key spin gaps in a [Mn3(IV)O4] cluster.
In the absence of a robust spin adaptation/purification technique,
both applications would be unattainable by Slater determinant based
ground state methods, with any starting wave function collapsing into
the higher-spin ground state during the optimization. This strategy
can be coupled to other algorithms that use the Slater determinant
based FCIQMC algorithm as configuration interaction eigensolver, including
the Stochastic Generalized Active Space, the similarity-transformed
FCIQMC, the tailored-CC, and second-order perturbation theory approaches.
Moreover, in contrast to the GUGA-FCIQMC technique, this strategy
features both spin projection and total spin adaptation, making it
appealing when solving anisotropic Hamiltonians. It also provides
spin-resolved reduced density matrices, important for the investigation
of spin-dependent properties in polynuclear transition metal clusters,
such as the hyperfine-coupling constants.
Collapse
Affiliation(s)
- Oskar Weser
- Max-Planck-Institute for Solid State Research, 70569 Stuttgart, Germany
| | - Niklas Liebermann
- Max-Planck-Institute for Solid State Research, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max-Planck-Institute for Solid State Research, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max-Planck-Institute for Solid State Research, 70569 Stuttgart, Germany.,Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Giovanni Li Manni
- Max-Planck-Institute for Solid State Research, 70569 Stuttgart, Germany
| |
Collapse
|
20
|
Affiliation(s)
- Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie Vladimir-Prelog-Weg 2 8093 Zürich Switzerland
| |
Collapse
|
21
|
Schraivogel T, Cohen AJ, Alavi A, Kats D. Transcorrelated coupled cluster methods. J Chem Phys 2021; 155:191101. [PMID: 34800963 DOI: 10.1063/5.0072495] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Transcorrelated coupled cluster and distinguishable cluster methods are presented. The Hamiltonian is similarity transformed with a Jastrow factor in the first quantization, which results in up to three-body integrals. The coupled cluster with singles and doubles equations on this transformed Hamiltonian are formulated and implemented. It is demonstrated that the resulting methods have a superior basis set convergence and accuracy to the corresponding conventional and explicitly correlated methods. Additionally, approximations for three-body integrals are suggested and tested.
Collapse
Affiliation(s)
- Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aron J Cohen
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
22
|
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
|
23
|
Guther K, Cohen AJ, Luo H, Alavi A. Binding curve of the beryllium dimer using similarity-transformed FCIQMC: Spectroscopic accuracy with triple-zeta basis sets. J Chem Phys 2021; 155:011102. [PMID: 34241381 DOI: 10.1063/5.0055575] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We demonstrate how similarity-transformed full configuration interaction quantum Monte Carlo (FCIQMC) based on the transcorrelated Hamiltonian can be applied to make highly accurate predictions for the binding curve of the beryllium dimer, marking the first case study of a molecular system with this method. In this context, the non-Hermitian transcorrelated Hamiltonian, resulting from a similarity transformation with a Jastrow factor, serves the purpose to effectively address dynamic correlation beyond the used basis set and thus allows for obtaining energies close to the complete basis set limit from FCIQMC already with moderate basis sets and computational effort. Building on results from other explicitly correlated methods, we discuss the role of the Jastrow factor and its functional form, as well as potential sources for size consistency errors, and arrive at Jastrow forms that allow for high accuracy calculations of the vibrational spectrum of the beryllium dimer.
Collapse
Affiliation(s)
- Kai Guther
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| | - Aron J Cohen
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| | - Hongjun Luo
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Federal Republic of Germany
| |
Collapse
|
24
|
Callahan JM, Lange MF, Berkelbach TC. Dynamical correlation energy of metals in large basis sets from downfolding and composite approaches. J Chem Phys 2021; 154:211105. [PMID: 34240964 DOI: 10.1063/5.0049890] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Coupled-cluster theory with single and double excitations (CCSD) is a promising ab initio method for the electronic structure of three-dimensional metals, for which second-order perturbation theory (MP2) diverges in the thermodynamic limit. However, due to the high cost and poor convergence of CCSD with respect to basis size, applying CCSD to periodic systems often leads to large basis set errors. In a common "composite" method, MP2 is used to recover the missing dynamical correlation energy through a focal-point correction, but the inadequacy of finite-order perturbation theory for metals raises questions about this approach. Here, we describe how high-energy excitations treated by MP2 can be "downfolded" into a low-energy active space to be treated by CCSD. Comparing how the composite and downfolding approaches perform for the uniform electron gas, we find that the latter converges more quickly with respect to the basis set size. Nonetheless, the composite approach is surprisingly accurate because it removes the problematic MP2 treatment of double excitations near the Fermi surface. Using this method to estimate the CCSD correlation energy in the combined complete basis set and thermodynamic limits, we find that CCSD recovers 85%-90% of the exact correlation energy at rs = 4. We also test the composite approach with the direct random-phase approximation used in place of MP2, yielding a method that is typically (but not always) more cost effective due to the smaller number of orbitals that need to be included in the more expensive CCSD calculation.
Collapse
Affiliation(s)
- James M Callahan
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Malte F Lange
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
25
|
Jiang T, Chen Y, Bogdanov NA, Wang E, Alavi A, Chen J. A full configuration interaction quantum Monte Carlo study of ScO, TiO, and VO molecules. J Chem Phys 2021; 154:164302. [PMID: 33940817 DOI: 10.1063/5.0046464] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate ab initio calculations of 3d transition metal monoxide molecules have attracted extensive attention because of their relevance in physical and chemical science as well as theoretical challenges in treating strong electron correlation. Meanwhile, recent years have witnessed the rapid development of the full configuration interaction quantum Monte Carlo (FCIQMC) method to tackle electron correlation. In this study, we carry out FCIQMC simulations to ScO, TiO, and VO molecules and obtain accurate descriptions of 13 low-lying electronic states (ScO 2Σ+, 2Δ, 2Π; TiO 3Δ, 1Δ, 1Σ+, 3Π, 3Φ; VO 4Σ-, 4Φ, 4Π, 2Γ, 2Δ), including states that have significant multi-configurational character. The FCIQMC results are used to assess the performance of several other wave function theory and density functional theory methods. Our study highlights the challenging nature of the electronic structure of transition metal oxides and demonstrates FCIQMC as a promising technique going forward to treat more complex transition metal oxide molecules and materials.
Collapse
Affiliation(s)
- Tonghuan Jiang
- School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Yilin Chen
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Nikolay A Bogdanov
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Enge Wang
- International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, People's Republic of China
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Ji Chen
- School of Physics, Peking University, Beijing 100871, People's Republic of China
| |
Collapse
|
26
|
Eriksen JJ, Gauss J. Incremental treatments of the full configuration interaction problem. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1525] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Jürgen Gauss
- Department Chemie Johannes Gutenberg‐Universität Mainz Mainz Germany
| |
Collapse
|
27
|
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
|
28
|
Ghanem K, Guther K, Alavi A. The adaptive shift method in full configuration interaction quantum Monte Carlo: Development and applications. J Chem Phys 2020; 153:224115. [DOI: 10.1063/5.0032617] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Khaldoon Ghanem
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Kai Guther
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
29
|
Abstract
We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG with the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
30
|
Bircher MP, Villard J, Rothlisberger U. Efficient Treatment of Correlation Energies at the Basis-Set Limit by Monte Carlo Summation of Continuum States. J Chem Theory Comput 2020; 16:6550-6559. [PMID: 32915565 PMCID: PMC7584365 DOI: 10.1021/acs.jctc.0c00724] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Martin P. Bircher
- Computational and Soft Matter Physics, Universität Wien, Sensengasse 8/9, A-1090 Wien, Austria
| | - Justin Villard
- Laboratory of Computational Chemistry and Biochemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Av. F.A. Forel 2, CH-1015 Lausanne, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Av. F.A. Forel 2, CH-1015 Lausanne, Switzerland
| |
Collapse
|
31
|
Guther K, Anderson RJ, Blunt NS, Bogdanov NA, Cleland D, Dattani N, Dobrautz W, Ghanem K, Jeszenszki P, Liebermann N, Manni GL, Lozovoi AY, Luo H, Ma D, Merz F, Overy C, Rampp M, Samanta PK, Schwarz LR, Shepherd JJ, Smart SD, Vitale E, Weser O, Booth GH, Alavi A. NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods. J Chem Phys 2020; 153:034107. [DOI: 10.1063/5.0005754] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
- Kai Guther
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Robert J. Anderson
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Nick S. Blunt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Nikolay A. Bogdanov
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | | | - Nike Dattani
- Department of Electrical and Computer Engineering, University of Waterloo, 200 University Avenue, Waterloo, Ontario N2L 3G1, Canada
| | - Werner Dobrautz
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Khaldoon Ghanem
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Peter Jeszenszki
- Centre for Theoretical Chemistry and Physics, NZ Institute for Advanced Study, Massey University, Auckland, New Zealand
- Dodd-Walls Centre for Photonic and Quantum Technologies, P.O. Box 56, Dunedin 9056, New Zealand
| | - Niklas Liebermann
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Giovanni Li Manni
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Alexander Y. Lozovoi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Hongjun Luo
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Dongxia Ma
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Florian Merz
- Lenovo HPC and AI Innovation Center, Meitnerstr. 9, 70563 Stuttgart, Germany
| | - Catherine Overy
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Markus Rampp
- Max Planck Computing and Data Facility (MPCDF), Gießenbachstr. 2, 85748 Garching, Germany
| | - Pradipta Kumar Samanta
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Lauretta R. Schwarz
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - James J. Shepherd
- Department of Chemistry and Informatics Institute, University of Iowa, Iowa City, Iowa 52242, USA
| | - Simon D. Smart
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Eugenio Vitale
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Oskar Weser
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
32
|
Needs RJ, Towler MD, Drummond ND, López Ríos P, Trail JR. Variational and diffusion quantum Monte Carlo calculations with the CASINO code. J Chem Phys 2020; 152:154106. [DOI: 10.1063/1.5144288] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- R. J. Needs
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - M. D. Towler
- University College London, London WC1E 6BT, United Kingdom
| | - N. D. Drummond
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
| | - P. López Ríos
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - J. R. Trail
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
33
|
Baiardi A, Reiher M. The density matrix renormalization group in chemistry and molecular physics: Recent developments and new challenges. J Chem Phys 2020; 152:040903. [DOI: 10.1063/1.5129672] [Citation(s) in RCA: 85] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
34
|
Neufeld VA, Thom AJW. Accelerating Convergence in Fock Space Quantum Monte Carlo Methods. J Chem Theory Comput 2020; 16:1503-1510. [DOI: 10.1021/acs.jctc.9b01023] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Verena A. Neufeld
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
35
|
Ghanem K, Lozovoi AY, Alavi A. Unbiasing the initiator approximation in full configuration interaction quantum Monte Carlo. J Chem Phys 2019; 151:224108. [DOI: 10.1063/1.5134006] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Khaldoon Ghanem
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Alexander Y. Lozovoi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
36
|
Petras HR, Graham DS, Ramadugu SK, Goodpaster JD, Shepherd JJ. Fully Quantum Embedding with Density Functional Theory for Full Configuration Interaction Quantum Monte Carlo. J Chem Theory Comput 2019; 15:5332-5342. [DOI: 10.1021/acs.jctc.9b00571] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Hayley R. Petras
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
- University of Iowa Informatics Initiative, University of Iowa, Iowa City, Iowa 52242, United States
| | - Daniel S. Graham
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Sai Kumar Ramadugu
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
- University of Iowa Informatics Initiative, University of Iowa, Iowa City, Iowa 52242, United States
| | - Jason D. Goodpaster
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - James J. Shepherd
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
- University of Iowa Informatics Initiative, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
37
|
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
|
38
|
Cohen AJ, Luo H, Guther K, Dobrautz W, Tew DP, Alavi A. Similarity transformation of the electronic Schrödinger equation via Jastrow factorization. J Chem Phys 2019. [DOI: 10.1063/1.5116024] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Aron J. Cohen
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Hongjun Luo
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Kai Guther
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Werner Dobrautz
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - David P. Tew
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
39
|
Spencer JS, Blunt NS, Choi S, Etrych J, Filip MA, Foulkes WMC, Franklin RST, Handley WJ, Malone FD, Neufeld VA, Di Remigio R, Rogers TW, Scott CJC, Shepherd JJ, Vigor WA, Weston J, Xu R, Thom AJW. The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up. J Chem Theory Comput 2019; 15:1728-1742. [DOI: 10.1021/acs.jctc.8b01217] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- James S. Spencer
- Department of Physics, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
- Department of Materials, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Nick S. Blunt
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- St. John’s College, St. John’s Street, Cambridge, CB2 1TP, United Kingdom
| | - Seonghoon Choi
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Jiří Etrych
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Maria-Andreea Filip
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - W. M. C. Foulkes
- Department of Physics, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
- Department of Physics, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Ruth S. T. Franklin
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Will J. Handley
- Astrophysics Group, Cavendish Laboratory, Cambridge, CB3 OHE, United Kingdom
- Kavli Institute for Cosmology, Madingley Road, Cambridge, CB3 0HA, United Kingdom
- Gonville & Caius College, Trinity Street, Cambridge, CB2 1TA, United Kingdom
| | - Fionn D. Malone
- Department of Physics, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Verena A. Neufeld
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø—The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Thomas W. Rogers
- Department of Physics, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Charles J. C. Scott
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | | | - William A. Vigor
- Department of Chemistry, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - Joseph Weston
- Department of Physics, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - RuQing Xu
- Department of Modern Physics, University of Science and Technology, Hefei, Anhui 230026, China
| | - Alex J. W. Thom
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom
- Department of Chemistry, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
40
|
Neufeld VA, Thom AJW. Exciting Determinants in Quantum Monte Carlo: Loading the Dice with Fast, Low-Memory Weights. J Chem Theory Comput 2018; 15:127-140. [DOI: 10.1021/acs.jctc.8b00844] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Verena A. Neufeld
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|