1
|
Wang Q, Agarawal V, Hermes MR, Motta M, Rice JE, Jones GO, Gagliardi L. Distinguishing homolytic vs heterolytic bond dissociation of phenylsulfonium cations with localized active space methods. J Chem Phys 2024; 161:014106. [PMID: 38949581 DOI: 10.1063/5.0215697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.
Collapse
Affiliation(s)
- Qiaohong Wang
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - Valay Agarawal
- Department of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, USA
| | - Matthew R Hermes
- Department of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, USA
| | - Mario Motta
- IBM Quantum, IBM T. J. Watson Research Center, Yorktown Heights, New York 1059, USA
| | - Julia E Rice
- IBM Quantum, IBM Research-Almaden, San Jose, California 95120, USA
| | - Gavin O Jones
- IBM Quantum, IBM Research-Almaden, San Jose, California 95120, USA
| | - Laura Gagliardi
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
- Department of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, USA
- James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
- Argonne National Laboratory, 9700 S, Cass Avenue, Lemont, Illinois 60439, USA
| |
Collapse
|
2
|
Greiner J, Gianni I, Nottoli T, Lipparini F, Eriksen JJ, Gauss J. MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces. J Chem Theory Comput 2024; 20:4663-4675. [PMID: 38809011 DOI: 10.1021/acs.jctc.4c00388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals.
Collapse
Affiliation(s)
- Jonas Greiner
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz 55128, Germany
| | - Ivan Gianni
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Tommaso Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Janus J Eriksen
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg. 206, Kgs. Lyngby 2800, Denmark
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz 55128, Germany
| |
Collapse
|
3
|
Agarawal V, King DS, Hermes MR, Gagliardi L. Automatic State Interaction with Large Localized Active Spaces for Multimetallic Systems. J Chem Theory Comput 2024; 20:4654-4662. [PMID: 38787596 DOI: 10.1021/acs.jctc.4c00376] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The localized active space self-consistent field method factorizes a complete active space wave function into an antisymmetrized product of localized active space wave function fragments. Correlation between fragments is then reintroduced through localized active space state interaction (LASSI), in which the Hamiltonian is diagonalized in a model space of LAS states. However, the optimal procedure for defining the LAS fragments and LASSI model space is unknown. We here present an automated framework to explore systematically convergent sets of model spaces, which we call LASSI[r, q]. This method requires the user to select only r, the number of electron hops from one fragment to another, and q, the number of fragment basis functions per Hilbert space, which converges to CASCI in the limit of r, q → ∞. Numerical tests of this method on the trimetal oxo-centered complexes [Fe(III)Al(III)Fe(II)(μ3-O)(HCOO)6] and [Fe(III)2Fe(II)(μ3-O)(HCOO)6] show efficient convergence to the CASCI limit with 4-10 orders of magnitude fewer states than CASCI.
Collapse
Affiliation(s)
- Valay Agarawal
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Daniel S King
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R Hermes
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
4
|
Pham HQ, Ouyang R, Lv D. Scalable Quantum Monte Carlo with Direct-Product Trial Wave Functions. J Chem Theory Comput 2024; 20:3524-3534. [PMID: 38700513 DOI: 10.1021/acs.jctc.3c00769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2024]
Abstract
The computational demand posed by applying multi-Slater determinant trials in phaseless auxiliary-field quantum Monte Carlo methods (MSD-AFQMC) is particularly significant for molecules exhibiting strong correlations. Here, we propose using direct-product wave functions as trials for MSD-AFQMC, aiming to reduce computational overhead by leveraging the compactness of multi-Slater determinant trials in direct-product form (DP-MSD). This efficiency arises when the active space can be divided into noncoupling subspaces, a condition we term "decomposable active space". By employing localized-active space self-consistent field wave functions as an example of such trials, we demonstrate our proposed approach across a range of molecular systems, each exhibiting varying degrees of complexity in their electronic structures. Our findings indicate that the compact DP-MSD trials can reduce computational costs substantially, by up to 36 times for the C2H6N4 molecule where the two double bonds between nitrogen N=N are clearly separated by a C-C single bond, while maintaining accuracy when active spaces are decomposable. In the case of larger systems such as the benzene dimer, characterized by weak coupling between the two monomers, we observed a decrease in computational cost compared to using a complete active space trial, yet we retained the same level of accuracy. However, for systems where these active subspaces strongly couple, a scenario we refer to as "strong subspace coupling", the method's accuracy decreases compared to that achieved with a complete active space approach. We anticipate that our method will be beneficial for systems with noncoupling to weakly coupling subspaces that require local multireference treatments.
Collapse
Affiliation(s)
- Hung Q Pham
- ByteDance Research, San Jose, California 95110, United States
| | | | - Dingshun Lv
- ByteDance Research, Zhonghang Plaza, No. 43, North third Ring West Road, 100098 Beijing, China
| |
Collapse
|
5
|
D'Cunha R, Otten M, Hermes MR, Gagliardi L, Gray SK. State Preparation in Quantum Algorithms for Fragment-Based Quantum Chemistry. J Chem Theory Comput 2024; 20:3121-3130. [PMID: 38607377 DOI: 10.1021/acs.jctc.3c01283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2024]
Abstract
State preparation for quantum algorithms is crucial for achieving high accuracy in quantum chemistry and competing with classical algorithms. The localized active space-unitary coupled cluster (LAS-UCC) algorithm iteratively loads a fragment-based multireference wave function onto a quantum computer. In this study, we compare two state preparation methods, quantum phase estimation (QPE) and direct initialization (DI), for each fragment. We test the two state preparation methods on three systems, ranging from a model system, a set of interacting hydrogen molecules, to more realistic chemical problems, like the C-C double bond breaking in transbutadiene and the spin ladder in a bimetallic system. We analyze the impact of QPE parameters, such as the number of ancilla qubits and Trotter steps, on the prepared state. We find a trade-off between the methods, where DI requires fewer resources for smaller fragments, while QPE is more efficient for larger fragments. Our resource estimates highlight the benefits of system fragmentation in state preparation for subsequent quantum chemical calculations. These findings have broad applications for preparing multireference quantum chemical wave functions on quantum circuits that can be used for realistic chemical applications.
Collapse
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew Otten
- HRL Laboratories LLC, Malibu, California 90265, United States
| | - Matthew R Hermes
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Argonne National Laboratory, Lemont, Illinois 60439, United States
- James Franck Institute, Chicago, Illinois 60637, United States
- Pritzker School of Molecular Engineering, Chicago, Illinois 60637, United States
| | - Stephen K Gray
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
6
|
Braunscheidel NM, Abraham V, Mayhall NJ. Generalization of the Tensor Product Selected CI Method for Molecular Excited States. J Phys Chem A 2023; 127:8179-8193. [PMID: 37733948 DOI: 10.1021/acs.jpca.3c03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/23/2023]
Abstract
In a recent paper [JCTC, 2020, 16, 6098], we introduced a new approach for accurately approximating full CI ground states in large electronic active-spaces called Tensor Product Selected CI (TPSCI). In TPSCI, a large orbital active space is first partitioned into disjoint sets (clusters) for which the exact, local many-body eigenstates are obtained. Tensor products of these locally correlated many-body states are taken as the basis for the full, global Hilbert space. By folding correlation into the basis states themselves, the low-energy eigenstates become increasingly sparse, creating a more compact selected CI expansion. While we demonstrated that this approach can improve accuracy for a variety of systems, there is even greater potential for applications to excited states, particularly those which have some excited-state character. In this paper, we report on the accuracy of TPSCI for excited states, including a far more efficient implementation in the Julia programming language. In traditional SCI methods that use a Slater determinant basis, accurate excitation energies are obtained only after a linear extrapolation and at a large computational cost. We find that TPSCI with perturbative corrections provides accurate excitation energies for several excited states of various polycyclic aromatic hydrocarbons with respect to the extrapolated result (i.e., near exact result). Further, we use TPSCI to report highly accurate estimates of the lowest 31 eigenstates for a tetracene tetramer system with an active space of 40 electrons in 40 orbitals, giving direct access to the initial bright states and the resulting 18 doubly excited (biexcitonic) states.
Collapse
Affiliation(s)
| | - Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, United States
| |
Collapse
|
7
|
Sekaran S, Bindech O, Fromager E. A unified density matrix functional construction of quantum baths in density matrix embedding theory beyond the mean-field approximation. J Chem Phys 2023; 159:034107. [PMID: 37466226 DOI: 10.1063/5.0157746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Accepted: 07/03/2023] [Indexed: 07/20/2023] Open
Abstract
The equivalence in one-electron quantum baths between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [S. Yalouz et al., J. Chem. Phys. 157, 214112 (2022)]. We prove mathematically that the equivalence remains valid when the density matrix is not idempotent anymore, thus allowing for the construction of correlated (one-electron) quantum baths. A density-matrix functional exactification of DMET is derived within the present unified quantum embedding formalism. Numerical examples reveal that the embedding cluster can be quite sensitive to the level of density-matrix functional approximation used for computing the reference density matrix.
Collapse
Affiliation(s)
- Sajanthan Sekaran
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| | - Oussama Bindech
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France
| |
Collapse
|
8
|
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
|
9
|
Otten M, Hermes MR, Pandharkar R, Alexeev Y, Gray SK, Gagliardi L. Localized Quantum Chemistry on Quantum Computers. J Chem Theory Comput 2022; 18:7205-7217. [PMID: 36346785 DOI: 10.1021/acs.jctc.2c00388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Quantum chemistry calculations of large, strongly correlated systems are typically limited by the computation cost that scales exponentially with the size of the system. Quantum algorithms, designed specifically for quantum computers, can alleviate this, but the resources required are still too large for today's quantum devices. Here, we present a quantum algorithm that combines a localization of multireference wave functions of chemical systems with quantum phase estimation (QPE) and variational unitary coupled cluster singles and doubles (UCCSD) to compute their ground-state energy. Our algorithm, termed "local active space unitary coupled cluster" (LAS-UCC), scales linearly with the system size for certain geometries, providing a polynomial reduction in the total number of gates compared with QPE, while providing accuracy above that of the variational quantum eigensolver using the UCCSD ansatz and also above that of the classical local active space self-consistent field. The accuracy of LAS-UCC is demonstrated by dissociating (H2)2 into two H2 molecules and by breaking the two double bonds in trans-butadiene, and resource estimates are provided for linear chains of up to 20 H2 molecules.
Collapse
Affiliation(s)
- Matthew Otten
- HRL Laboratories, LLC, 3011 Malibu Canyon Road, Malibu, California90265, United States
| | - Matthew R Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois60637, United States
| | - Riddhish Pandharkar
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois60637, United States
| | - Yuri Alexeev
- Computational Science Division, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Stephen K Gray
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, Illinois60439, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois60637, United States.,Argonne National Laboratory, Lemont, Illinois60439, United States
| |
Collapse
|
10
|
Ai Y, Sun Q, Jiang H. Efficient Multiconfigurational Quantum Chemistry Approach to Single-Ion Magnets Based on Density Matrix Embedding Theory. J Phys Chem Lett 2022; 13:10627-10634. [PMID: 36350106 DOI: 10.1021/acs.jpclett.2c02890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Density matrix embedding theory (DMET) provides a systematic framework to combine low-level (e.g., Hartree-Fock approximation) and high-level correlated quantum chemistry methods to treat strongly correlated systems. In this work, we propose an efficient quantum embedding approach that combines DMET with the complete active space self-consistent field and subsequent state interaction treatment of spin-orbit coupling (CASSI-SO) and apply it to a theoretical description of single-ion magnets (SIMs). We have developed a novel regularized direct inversion of iterative subspace (R-DIIS) technique that ensures restricted open-shell Hartree-Fock converging to a physically correct ground state, which is found to be crucial for the efficacy of subsequent CASSI-SO calculation. The DMET+CASSI-SO approach can produce reliable zero-field splitting parameters in typical 3d-SIMs with dramatically reduced computational cost compared to its all-electron counterpart. This work therefore demonstrates the great potential of DMET-based multiconfigurational approaches for efficient ab initio study of magneto-structural correlations in complex molecular magnetic systems.
Collapse
Affiliation(s)
- Yuhang Ai
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Qiming Sun
- Axiomquant Investment Management LLC, Beijing 100871, China
| | - Hong Jiang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
11
|
Pandharkar R, Hermes MR, Cramer CJ, Gagliardi L. Localized Active Space-State Interaction: a Multireference Method for Chemical Insight. J Chem Theory Comput 2022; 18:6557-6566. [PMID: 36257065 DOI: 10.1021/acs.jctc.2c00536] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Multireference electronic structure methods, like the complete active space (CAS) self-consistent field model, have long been used to characterize chemically interesting processes. Important work has been done in recent years to develop modifications having a lower computational cost than CAS, but typically these methods offer no more chemical insight than that from the CAS solution being approximated. In this paper, we present the localized active space-state interaction (LASSI) method that can be used not only to lower the intrinsic cost of the multireference calculation but also to improve interpretability. The localized active space (LAS) approach utilizes the local nature of the electron-electron correlation to express a composite wave function as an antisymmetrized product of unentangled wave functions in local active subspaces. LASSI then uses these LAS states as a basis from which to express complete molecular wave functions. This not only makes the molecular wave function more compact but also permits flexibility in choosing those states to be included in the basis. Such selective inclusion of states translates to the selective inclusion of specific types of interactions, thereby allowing a quantitative analysis of these interactions. We demonstrate the use of LASSI to study charge migration and spin-flip excitations in multireference organic molecules. We also compute the J coupling parameter for a bimetallic compound using various LAS bases to construct the Hamiltonian to provide insights into the coupling mechanism.
Collapse
Affiliation(s)
- Riddhish Pandharkar
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois60637, United States.,Argonne National Laboratory, Lemont, Illinois60439, USA
| | - Matthew R Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois60637, United States
| | - Christopher J Cramer
- Underwriters Laboratories Inc., 333 Pfingsten Road., Northbrook, Illinois60062, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois60637, United States.,Argonne National Laboratory, Lemont, Illinois60439, USA
| |
Collapse
|
12
|
Mitra A, Hermes MR, Cho M, Agarawal V, Gagliardi L. Periodic Density Matrix Embedding for CO Adsorption on the MgO(001) Surface. J Phys Chem Lett 2022; 13:7483-7489. [PMID: 35939641 PMCID: PMC9393885 DOI: 10.1021/acs.jpclett.2c01915] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Accepted: 08/01/2022] [Indexed: 05/19/2023]
Abstract
The adsorption of simple gas molecules to metal oxide surfaces is a primary step in many heterogeneous catalysis applications. Quantum chemical modeling of these reactions is a challenge in terms of both cost and accuracy, and quantum-embedding methods are promising, especially for localized chemical phenomena. In this work, we employ density matrix embedding theory (DMET) for periodic systems to calculate the adsorption energy of CO to the MgO(001) surface. Using coupled-cluster theory with single and double excitations and second-order Møller-Plesset perturbation theory as quantum chemical solvers, we perform calculations with embedding clusters up to 266 electrons in 306 orbitals, with the largest embedding models agreeing to within 1.2 kcal/mol of the non-embedding references. Moreover, we present a memory-efficient procedure of storing and manipulating electron repulsion integrals in the embedding space within the framework of periodic DMET.
Collapse
Affiliation(s)
- Abhishek Mitra
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Minsik Cho
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Valay Agarawal
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Laura Gagliardi
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
- Argonne
National Laboratory 9700
South Cass Avenue Lemont, Illinois 60439, United
States
- E-mail:
| |
Collapse
|
13
|
Abraham V, Mayhall NJ. Coupled Electron Pair-Type Approximations for Tensor Product State Wave Functions. J Chem Theory Comput 2022; 18:4856-4864. [PMID: 35878319 DOI: 10.1021/acs.jctc.2c00589] [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
Size extensivity, defined as the correct scaling of energy with system size, is a desirable property for any many-body method. Traditional configuration interaction (CI) methods are not size extensive, hence the error increases as the system gets larger. Coupled electron pair approximation (CEPA) methods can be constructed as simple extensions of a truncated CI that ensures size extensivity. One of the major issues with the CEPA and its variants is that singularities arise in the amplitude equations when the system starts to be strongly correlated. In this work, we extend the traditional Slater determinant based coupled electron pair approaches like CEPA-0, averaged coupled-pair functional, and average quadratic coupled-cluster to a new formulation based on tensor product states (TPS). We show that a TPS basis can often be chosen such that it removes the singularities that commonly destroy the accuracy of CEPA based methods. A suitable TPS representation can be formed by partitioning the system into separate disjoint clusters and forming the final wave function as the tensor product of the many body states of these clusters. We demonstrate the application of these methods on simple bond breaking systems such as CH4 and F2 where determinant based CEPA methods fail. We further apply the TPS-CEPA approach to stillbene isomerization and few planar π-conjugated systems. Overall, the results show that the TPS-CEPA method can remove the singularities and provide improved numerical results compared to common electronic structure methods.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, United States
| |
Collapse
|
14
|
Zhou C, Hermes MR, Wu D, Bao JJ, Pandharkar R, King DS, Zhang D, Scott TR, Lykhin AO, Gagliardi L, Truhlar DG. Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory. Chem Sci 2022; 13:7685-7706. [PMID: 35865899 PMCID: PMC9261488 DOI: 10.1039/d2sc01022d] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Accepted: 06/01/2022] [Indexed: 12/25/2022] Open
Abstract
Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals. This feature article overviews recent work on active spaces, matrix product reference states, treatment of quasidegeneracy, hybrid theory, density-coherence functionals, machine-learned functionals, spin–orbit coupling, gradients, and dipole moments.![]()
Collapse
Affiliation(s)
- Chen Zhou
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Matthew R. Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Dihua Wu
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Jie J. Bao
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Riddhish Pandharkar
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Daniel S. King
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Dayou Zhang
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| | - Thais R. Scott
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Aleksandr O. Lykhin
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, IL 60637, USA
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, MN 55455-0431, USA
| |
Collapse
|
15
|
Pandharkar R, Hermes MR, Cramer CJ, Truhlar DG, Gagliardi L. Localized Active Space Pair-Density Functional Theory. J Chem Theory Comput 2021; 17:2843-2851. [PMID: 33900078 DOI: 10.1021/acs.jctc.1c00067] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.
Collapse
Affiliation(s)
- Riddhish Pandharkar
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Matthew R Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, 5735 S Ellis Avenue, Chicago, Illinois 60637, United States
| |
Collapse
|
16
|
Senjean B, Sen S, Repisky M, Knizia G, Visscher L. Generalization of Intrinsic Orbitals to Kramers-Paired Quaternion Spinors, Molecular Fragments, and Valence Virtual Spinors. J Chem Theory Comput 2021; 17:1337-1354. [PMID: 33555866 DOI: 10.1021/acs.jctc.0c00964] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Localization of molecular orbitals finds its importance in the representation of chemical bonding (and antibonding) and in the local correlation treatments beyond mean-field approximation. In this paper, we generalize the intrinsic atomic and bonding orbitals [G. Knizia, J. Chem. Theory Comput. 2013, 9, 11, 4834-4843] to relativistic applications using complex and quaternion spinors, as well as to molecular fragments instead of atomic fragments only. By performing a singular value decomposition, we show how localized valence virtual orbitals can be expressed on this intrinsic minimal basis. We demonstrate our method on systems of increasing complexity, starting from simple cases such as benzene, acrylic acid, and ferrocene molecules, and then demonstrate the use of molecular fragments and inclusion of relativistic effects for complexes containing heavy elements such as tellurium, iridium, and astatine. The aforementioned scheme is implemented into a standalone program interfaced with several different quantum chemistry packages.
Collapse
Affiliation(s)
- Bruno Senjean
- Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands
| | - Souloke Sen
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Gerald Knizia
- Department of Chemistry, The Pennsylvania State University, University Park, State College, Pennsylvania 16802, United States
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
17
|
Head-Marsden K, Flick J, Ciccarino CJ, Narang P. Quantum Information and Algorithms for Correlated Quantum Matter. Chem Rev 2020; 121:3061-3120. [PMID: 33326218 DOI: 10.1021/acs.chemrev.0c00620] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Johannes Flick
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Christopher J Ciccarino
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Prineha Narang
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
18
|
Abraham V, Mayhall NJ. Selected Configuration Interaction in a Basis of Cluster State Tensor Products. J Chem Theory Comput 2020; 16:6098-6113. [DOI: 10.1021/acs.jctc.0c00141] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, United States
| | - Nicholas J. Mayhall
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, United States
| |
Collapse
|