1
|
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
|
2
|
Papastathopoulos-Katsaros A, Henderson TM, Scuseria GE. Linear Combinations of Cluster Mean-Field States Applied to Spin Systems. J Chem Theory Comput 2024; 20:3697-3705. [PMID: 38695526 DOI: 10.1021/acs.jctc.4c00171] [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
We present an innovative cluster-based method employing linear combinations of diverse cluster mean-field states and apply it to describe the ground state of strongly correlated spin systems. In cluster mean-field theory, the ground state wave function is expressed as a factorized tensor product of optimized cluster states. While our prior work concentrated on a single cMF tiling, this study removes that constraint by combining different tilings of cMF states. Selection criteria, including translational symmetry and spatial proximity, guide this process. We present benchmark calculations for the one- and two-dimensional J1 - J2 and XXZ Heisenberg models. Our findings highlight two key aspects. First, the method offers a semiquantitative description of the 0.4 ≲ J2/J1 ≲ 0.6 regime of the J1 - J2 model─a particularly challenging regime for existing methods. Second, our results demonstrate the capability of our method to provide qualitative descriptions for all the models and regimes considered, establishing it as a valuable reference. However, the inclusion of additional (weak) correlations is necessary for quantitative agreement, and we explore methods to incorporate these extra correlations.
Collapse
Affiliation(s)
| | - Thomas M Henderson
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
| | - Gustavo E Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
3
|
Vu N, Mejia-Rodriguez D, Bauman NP, Panyala A, Mutlu E, Govind N, Foley JJ. Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory. J Chem Theory Comput 2024; 20:1214-1227. [PMID: 38291561 PMCID: PMC10876286 DOI: 10.1021/acs.jctc.3c01207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 12/21/2023] [Accepted: 12/26/2023] [Indexed: 02/01/2024]
Abstract
Polariton chemistry has attracted great attention as a potential route to modify chemical structure, properties, and reactivity through strong interactions among molecular electronic, vibrational, or rovibrational degrees of freedom. A rigorous theoretical treatment of molecular polaritons requires the treatment of matter and photon degrees of freedom on equal quantum mechanical footing. In the limit of molecular electronic strong or ultrastrong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of ab initio quantum chemistry, yielding an approach we refer to as ab initio cavity quantum electrodynamics, where the photon degrees of freedom are treated at the level of cavity quantum electrodynamics. Here, we present an approach called Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction theory to provide ground- and excited-state polaritonic surfaces with a balanced description of strong correlation effects among electronic and photonic degrees of freedom. This method provides a platform for ab initio cavity quantum electrodynamics when both strong electron correlation and strong light-matter coupling are important and is an important step toward computational approaches that yield multiple polaritonic potential energy surfaces and couplings that can be leveraged for ab initio molecular dynamics simulations of polariton chemistry.
Collapse
Affiliation(s)
- Nam Vu
- Department
of Chemistry, University of North Carolina
Charlotte, 9201 University City Blvd., Charlotte, North Carolina 28223, United States
| | - Daniel Mejia-Rodriguez
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Nicholas P. Bauman
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Ajay Panyala
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Erdal Mutlu
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jonathan J. Foley
- Department
of Chemistry, University of North Carolina
Charlotte, 9201 University City Blvd., Charlotte, North Carolina 28223, United States
| |
Collapse
|
4
|
Gao H, Imamura S, Kasagi A, Yoshida E. Distributed Implementation of Full Configuration Interaction for One Trillion Determinants. J Chem Theory Comput 2024; 20:1185-1192. [PMID: 38314701 PMCID: PMC10867839 DOI: 10.1021/acs.jctc.3c01190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 12/24/2023] [Accepted: 01/12/2024] [Indexed: 02/07/2024]
Abstract
Full configuration interaction (FCI) can provide an exact molecular ground-state energy within a given basis set and serve as a benchmark for approximate methods in quantum chemical calculations, including the emerging variational quantum eigensolver. However, its exponential computational and memory requirements easily exceed the capability of a single server and limit its applicability to large molecules. In this paper, we present a distributed FCI implementation employing a hybrid parallelization scheme with multithreading and multiprocessing to expand FCI's applicability. We optimize this scheme to minimize the bottlenecks arising from interprocess communications and interthread data management. Our implementation achieves higher scalability than the naive combination of prior works and successfully calculates the exact energy of C3H8/STO-3G with 1.31 trillion determinants, which is the largest FCI calculation to the best of our knowledge. Furthermore, we provide a comprehensive list of FCI results with 136 combinations of molecules and basis sets for future evaluation and development of approximate methods.
Collapse
Affiliation(s)
- Hong Gao
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| | - Satoshi Imamura
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| | - Akihiko Kasagi
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| | - Eiji Yoshida
- Computing Laboratory, Fujitsu Laboratories, Fujitsu Limited, Kawasaki City, Kanagawa, 211-0053, Japan
| |
Collapse
|
5
|
Coe JP. Analytic Non-adiabatic Couplings for Selected Configuration Interaction via Approximate Degenerate Coupled Perturbed Hartree-Fock. J Chem Theory Comput 2023; 19:8053-8065. [PMID: 37939698 PMCID: PMC10687870 DOI: 10.1021/acs.jctc.3c00601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Revised: 10/24/2023] [Accepted: 10/25/2023] [Indexed: 11/10/2023]
Abstract
We use degenerate perturbation theory and assume that for degenerate pairs of orbitals, the coupled perturbed Hartree-Fock coefficients are symmetric in the degenerate basis to show [Formula: see text] is the only modification needed in the original molecular orbital basis. This enables us to develop efficient and accurate analytic nonadiabatic couplings between electronic states for selected configuration interactions (CIs). Even when the states belong to different irreducible representations, degenerate orbital pairs cannot be excluded by symmetry. For various excited states of carbon monoxide and trigonal planar ammonia, we benchmark the method against the full CI and find it to be accurate. We create a semi-numerical approach and use it to show that the analytic approach is correct even when a high-symmetry structure is distorted to break symmetry so that near degeneracies in orbitals occur. For a range of geometries of trigonal planar ammonia, we find that the analytic non-adiabatic couplings for selected CI can achieve sufficient accuracy using a small fraction of the full CI space.
Collapse
Affiliation(s)
- Jeremy P. Coe
- Institute of Chemical Sciences, School
of Engineering and Physical Sciences, Heriot-Watt
University, Edinburgh EH14 4AS, U.K.
| |
Collapse
|
6
|
Dey M, Ghosh D. Machine Learning the Quantum Mechanical Wave Function. J Phys Chem A 2023; 127:9159-9166. [PMID: 37906959 DOI: 10.1021/acs.jpca.3c05322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
Strongly correlated systems have been challenging to computational chemists for a long time. To solve these systems, multireference methods have been developed over the years. Recently, with the fast development of machine learning and artificial intelligence methods, these methods have also influenced the quest for optimal wave function ansatz. Machine learning approaches have been used in many different flavors. From this perspective, we will discuss the different milestones achieved in the use of machine learning for solving the quantum many body problem.
Collapse
Affiliation(s)
- Mandira Dey
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
7
|
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
|
8
|
Papastathopoulos-Katsaros A, Henderson TM, Scuseria GE. Symmetry-projected cluster mean-field theory applied to spin systems. J Chem Phys 2023; 159:084107. [PMID: 37610020 DOI: 10.1063/5.0155765] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Accepted: 07/11/2023] [Indexed: 08/24/2023] Open
Abstract
We introduce Sz spin-projection based on cluster mean-field theory and apply it to the ground state of strongly correlated spin systems. In cluster mean-fields, the ground state wavefunction is written as a factorized tensor product of optimized cluster states. In previous work, we have focused on unrestricted cluster mean-field, where each cluster is Sz symmetry adapted. We here remove this restriction by introducing a generalized cluster mean-field (GcMF) theory, where each cluster is allowed to access all Sz sectors, breaking Sz symmetry. In addition, a projection scheme is used to restore global Sz, which gives rise to the Sz spin-projected generalized cluster mean-field (SzGcMF). Both of these extensions contribute to accounting for inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and 2D J1 - J2 and XXZ Heisenberg models. Our results indicate that the new methods (GcMF and SzGcMF) provide a qualitative and semi-quantitative description of the Heisenberg lattices in the regimes considered, suggesting them as useful references for further inter-cluster correlations, which are discussed in this work.
Collapse
Affiliation(s)
| | - Thomas M Henderson
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| | - Gustavo E Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| |
Collapse
|
9
|
Zhang H, Zou J, Ren X, Li S. Equation-of-Motion Block-Correlated Coupled Cluster Method for Excited Electronic States of Strongly Correlated Systems. J Phys Chem Lett 2023:6792-6799. [PMID: 37478417 DOI: 10.1021/acs.jpclett.3c01474] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2023]
Abstract
An equation-of-motion block-correlated coupled cluster method based on the generalized valence bond wave function (EOM-GVB-BCCC) is proposed to describe low-lying excited states for strongly correlated systems. The EOM-GVB-BCCC2b method with up to two-pair correlation has been implemented and tested for a few strongly correlated systems. For a water hexamer with stretched O-H bonds, which is beyond the capability of the CASSCF method, EOM-GVB-BCCC2b provides very close results as the density matrix renormalization group (DMRG). For four conjugated diradical species with triplet ground states, we found that their vertical S-T gaps from EOM-GVB-BCCC2b are also quite consistent with the DMRG results. This new method is expected to be a promising theoretical tool for describing the low-lying excited states of strongly correlated systems with large active spaces.
Collapse
Affiliation(s)
- Haodong Zhang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Jingxiang Zou
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Xiaochuan Ren
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, People's Republic of China
| |
Collapse
|
10
|
Ma Y, Li Z, Chen X, Ding B, Li N, Lu T, Zhang B, Suo B, Jin Z. Machine-learning assisted scheduling optimization and its application in quantum chemical calculations. J Comput Chem 2023; 44:1174-1188. [PMID: 36648254 DOI: 10.1002/jcc.27075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/16/2022] [Accepted: 12/23/2022] [Indexed: 01/18/2023]
Abstract
Easy and effective usage of computational resources is crucial for scientific calculations, both from the perspectives of timeliness and economic efficiency. This work proposes a bi-level optimization framework to optimize the computational sequences. Machine-learning (ML) assisted static load-balancing, and different dynamic load-balancing algorithms can be integrated. Consequently, the computational and scheduling engine of the ParaEngine is developed to invoke optimized quantum chemical (QC) calculations. Illustrated benchmark calculations include high-throughput drug suit, solvent model, P38 protein, and SARS-CoV-2 systems. The results show that the usage rate of given computational resources for high throughput and large-scale fragmentation QC calculations can primarily profit, and faster accomplishing computational tasks can be expected when employing high-performance computing (HPC) clusters.
Collapse
Affiliation(s)
- Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - ZhiYing Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Xin Chen
- ShenZhen Bay Laboratory, Shenzhen, China
| | - Bowen Ding
- Institute of Chemistry, Chinese Academy of Sciences, Beijing, China
| | - Ning Li
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
- College of Chemistry and Materials Engineering, Wenzhou University, Wen Zhou, China
| | - Teng Lu
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Baohua Zhang
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - BingBing Suo
- Department of Physics, Northwest University, Xi'an, China
| | - Zhong Jin
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
11
|
Iino T, Shiozaki T, Yanai T. Algorithm for analytic nuclear energy gradients of state averaged DMRG-CASSCF theory with newly derived coupled-perturbed equations. J Chem Phys 2023; 158:054107. [PMID: 36754810 DOI: 10.1063/5.0130636] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
We present an algorithm for evaluating analytic nuclear energy gradients of the state-averaged density matrix renormalization group complete-active-space self-consistent field (SA-DMRG-CASSCF) theory based on the newly derived coupled-perturbed (CP) DMRG-CASSCF equations. The Lagrangian for the conventional SA-CASSCF analytic gradient theory is extended to the SA-DMRG-CASSCF variant that can fully consider a whole set of constraints on the parameters of multi-root canonical matrix product states formed at all the DMRG block configurations. An efficient algorithm to solve the CP-DMRG-CASSCF equations for determining the multipliers was developed. The complexity of the resultant analytic gradient algorithm is overall the same as that of the unperturbed SA-DMRG-CASSCF algorithm. In addition, a reduced-scaling approach was developed to directly compute the SA reduced density matrices (SA-RDMs) and their perturbed ones without calculating separate state-specific RDMs. As part of our implementation scheme, we neglect the term associated with the constraint on the active orbitals in terms of the active-active rotation in the Lagrangian. Thus, errors from the true analytic gradients may be caused in this scheme. The proposed gradient algorithm was tested with the spin-adapted implementation by checking how accurately the computed analytic energy gradients reproduce numerical gradients of the SA-DMRG-CASSCF energies using a common number of renormalized bases. The illustrative applications show that the errors are sufficiently small when using a typical number of the renormalized bases, which is required to attain adequate accuracy in DMRG's total energies.
Collapse
Affiliation(s)
- Tsubasa Iino
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Toru Shiozaki
- Quantum Simulation Technologies, Inc., Boston, Massachusetts 02135, USA
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
12
|
Abstract
We develop analytic gradients for selected configuration interaction wave functions. Despite all pairs of molecular orbitals now potentially having to be considered for the coupled perturbed Hartree-Fock equations, we show that degenerate orbital pairs belonging to different irreducible representations in the largest abelian subgroup do not need to be included and instabilities due to degeneracies are avoided. We introduce seminumerical gradients and use them to validate the analytic approach even when near degeneracies are present due to high-symmetry geometries being slightly distorted to break symmetry. The method is applied to carbon monoxide, ammonia, square planar H4, hexagonal planar H6, and methane for a range of bond lengths where we demonstrate that analytic gradients for selected configuration interaction can approach the quality of full configuration interaction yet only use a very small fraction of its determinants.
Collapse
|
13
|
Dang DK, Kammeraad JA, Zimmerman PM. Advances in Parallel Heat Bath Configuration Interaction. J Phys Chem A 2023; 127:400-411. [PMID: 36580361 DOI: 10.1021/acs.jpca.2c07949] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Heat-bath configuration interaction (HCI) is a deterministic method that approaches the full CI limit at greatly reduced computational cost. In this work, computational improvements to the HCI algorithm are introduced targeting speed, parallel efficiency, and memory requirements. The new implementation introduces a hash function to distribute determinants and takes advantage of MPI and OpenMP for parallelism allowing for a (22e,168o) active space to be studied, which explicitly includes 2.39 × 107 variational determinants and 8.95 × 1010 perturbative determinants. Benchmarks show up to 86% parallel efficiency of the perturbative step on 32 nodes (4096 cores) and a total efficiency of 74%. The new HCI implementation is benchmarked for accuracy against prior results and applied to study the triplet-quintet gap in the challenging [FeO(NH3)5]2+ complex.
Collapse
Affiliation(s)
- Duy-Khoi Dang
- Department of Chemistry, University of Michigan 930 North University Avenue Ann Arbor, Michigan 48109, United States
| | - Joshua A Kammeraad
- Department of Chemistry, University of Michigan 930 North University Avenue Ann Arbor, Michigan 48109, United States
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan 930 North University Avenue Ann Arbor, Michigan 48109, United States
| |
Collapse
|
14
|
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
|
15
|
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
|
16
|
Coe JP, Moreno Carrascosa A, Simmermacher M, Kirrander A, Paterson MJ. Efficient Computation of Two-Electron Reduced Density Matrices via Selected Configuration Interaction. J Chem Theory Comput 2022; 18:6690-6699. [PMID: 36198067 PMCID: PMC9648180 DOI: 10.1021/acs.jctc.2c00738] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We create an approach to efficiently calculate two-electron
reduced
density matrices (2-RDMs) using selected configuration interaction
wavefunctions. This is demonstrated using the specific example of
Monte Carlo configuration interaction (MCCI). The computation of the
2-RDMs is accelerated by using ideas from fast implementations of
full configuration interaction (FCI) and recent advances in implementing
the Slater–Condon rules using hardware bitwise operations.
This method enables a comparison of MCCI and truncated CI 2-RDMs with
FCI values for a range of molecules, which includes stretched bonds
and excited states. The accuracy in energies, wavefunctions, and 2-RDMs
is seen to exhibit a similar behavior. We find that MCCI can reach
sufficient accuracy of the 2-RDM using significantly fewer configurations
than truncated CI, particularly for systems with strong multireference
character.
Collapse
Affiliation(s)
- Jeremy P Coe
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, EdinburghEH14 4AS, U.K
| | | | - Mats Simmermacher
- EaStCHEM, School of Chemistry, University of Edinburgh, EdinburghEH9 3FJ, U.K
| | - Adam Kirrander
- EaStCHEM, School of Chemistry, University of Edinburgh, EdinburghEH9 3FJ, U.K
| | - Martin J Paterson
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, EdinburghEH14 4AS, U.K
| |
Collapse
|
17
|
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
|
18
|
Papastathopoulos-Katsaros A, Jiménez-Hoyos CA, Henderson TM, Scuseria GE. Coupled Cluster and Perturbation Theories Based on a Cluster Mean-Field Reference Applied to Strongly Correlated Spin Systems. J Chem Theory Comput 2022; 18:4293-4303. [PMID: 35729717 DOI: 10.1021/acs.jctc.2c00338] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We introduce perturbation and coupled-cluster theories based on a cluster mean-field reference for describing the ground state of strongly correlated spin systems. In cluster mean-field, the ground state wave function is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allow for a straightforward improvement which uses perturbation theory and coupled-cluster to account for intercluster correlations. We present benchmark calculations on the 1D chain and 2D square J1-J2 Heisenberg model, using cluster mean-field, perturbation theory, and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies accurately. Our results indicate that, with sufficiently large clusters, the correlated methods (cPT2, cPT4, and cCCSD) can provide a relatively accurate description of the Heisenberg model in the regimes considered, which suggests that the methods presented can be used for other strongly correlated systems. Some ways to improve upon the methods presented in this work are discussed.
Collapse
Affiliation(s)
| | - Carlos A Jiménez-Hoyos
- Department of Chemistry, Wesleyan University, Middletown, Connecticut 06459, United States
| | - Thomas M Henderson
- Department of Chemistry, Rice University, Houston, Texas 77005, United States.,Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
| | - Gustavo E Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005, United States.,Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
19
|
Larsson HR, Zhai H, Gunst K, Chan GKL. Matrix Product States with Large Sites. J Chem Theory Comput 2022; 18:749-762. [DOI: 10.1021/acs.jctc.1c00957] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Henrik R. Larsson
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Klaas Gunst
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
- Center for Molecular Modeling, Ghent University, Technologiepark 46, B-9052 Zwijnaarde, Belgium
- Department of Physics and Astronomy, Ghent University, Krijgslaan 281, S9, B-9000 Ghent, Belgium
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
20
|
Leszczyk A, Máté M, Legeza Ö, Boguslawski K. Assessing the Accuracy of Tailored Coupled Cluster Methods Corrected by Electronic Wave Functions of Polynomial Cost. J Chem Theory Comput 2021; 18:96-117. [PMID: 34965121 DOI: 10.1021/acs.jctc.1c00284] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various coupled cluster methods tailored by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled cluster (fpCC) methods, which are tailored by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N2, F2, C2, CN+, CO, BN, BO+, and Cr2) and molecules (ammonia, ethylene, cyclobutadiene, benzene, hydrogen chains, rings, and cuboids) for which the conventional single-reference coupled cluster singles and doubles (CCSD) method is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to accurate theoretical and/or experimental reference data. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.
Collapse
Affiliation(s)
- Aleksandra Leszczyk
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Mihály Máté
- Strongly Correlated Systems "Lendület" Research Group, Wigner Research Center for Physics, H-1525 Budapest, Hungary.,Department of Physics of Complex Systems, Eötvös Loránd University, Pf. 32, H-1518 Budapest, Hungary
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research Group, Wigner Research Center for Physics, H-1525 Budapest, Hungary.,Institute for Advanced Study, Technical University of Munich, 80333 Munich, Germany
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| |
Collapse
|
21
|
Jeong W, Gaggioli CA, Gagliardi L. Active Learning Configuration Interaction for Excited-State Calculations of Polycyclic Aromatic Hydrocarbons. J Chem Theory Comput 2021; 17:7518-7530. [PMID: 34787422 PMCID: PMC8675132 DOI: 10.1021/acs.jctc.1c00769] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Indexed: 11/30/2022]
Abstract
We present the active learning configuration interaction (ALCI) method for multiconfigurational calculations based on large active spaces. ALCI leverages the use of an active learning procedure to find important electronic configurations among the full configurational space generated within an active space. We tested it for the calculation of singlet-singlet excited states of acenes and pyrene using different machine learning algorithms. The ALCI method yields excitation energies within 0.2-0.3 eV from those obtained by traditional complete active-space configuration interaction (CASCI) calculations (affordable for active spaces up to 16 electrons in 16 orbitals) by including only a small fraction of the CASCI configuration space in the calculations. For larger active spaces (we tested up to 26 electrons in 26 orbitals), not affordable with traditional CI methods, ALCI captures the trends of experimental excitation energies. Overall, ALCI provides satisfactory approximations to large active-space wave functions with up to 10 orders of magnitude fewer determinants for the systems presented here. These ALCI wave functions are promising and affordable starting points for the subsequent second-order perturbation theory or pair-density functional theory calculations.
Collapse
Affiliation(s)
- WooSeok Jeong
- Department
of Chemistry, Nanoporous Materials Genome Center, Chemical Theory
Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Carlo Alberto Gaggioli
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, 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, Lemont, Illinois 60439, United States
| |
Collapse
|
22
|
Abraham V, Mayhall NJ. Revealing the Contest between Triplet-Triplet Exchange and Triplet-Triplet Energy Transfer Coupling in Correlated Triplet Pair States in Singlet Fission. J Phys Chem Lett 2021; 12:10505-10514. [PMID: 34677988 DOI: 10.1021/acs.jpclett.1c03217] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Understanding the separation of the correlated triplet pair state 1(TT) intermediate is critical for leveraging singlet fission to improve solar cell efficiency. This separation mechanism is dominated by two key interactions: (i) the exchange interaction (K) between the triplets which leads to the spin splitting of the biexciton state into 1(TT),3(TT) and 5(TT) states, and (ii) the triplet-triplet energy transfer integral (t) which enables the formation of the spatially separated (but still spin entangled) state 1(T···T). We develop a simple ab initio technique to compute both the biexciton exchange (K) and biexciton transfer coupling. Our key findings reveal new conditions for successful correlated triplet pair state dissociation. The biexciton exchange interaction needs to be ferromagnetic or negligible to the triplet energy transfer for favorable dissociation. We also explore the effect of chromophore packing to reveal geometries where these conditions are achieved for tetracene.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24060, United States
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24060, United States
| |
Collapse
|
23
|
Goings JJ, Hu H, Yang C, Li X. Reinforcement Learning Configuration Interaction. J Chem Theory Comput 2021; 17:5482-5491. [PMID: 34423637 DOI: 10.1021/acs.jctc.1c00010] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Selected configuration interaction (sCI) methods exploit the sparsity of the full configuration interaction (FCI) wave function, yielding significant computational savings and wave function compression without sacrificing the accuracy. Despite recent advances in sCI methods, the selection of important determinants remains an open problem. We explore the possibility of utilizing reinforcement learning approaches to solve the sCI problem. By mapping the configuration interaction problem onto a sequential decision-making process, the agent learns on-the-fly which determinants to include and which to ignore, yielding a compressed wave function at near-FCI accuracy. This method, which we call reinforcement-learned configuration interaction, adds another weapon to the sCI arsenal and highlights how reinforcement learning approaches can potentially help solve challenging problems in electronic structure theory.
Collapse
Affiliation(s)
- Joshua J Goings
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hang Hu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Chao Yang
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
24
|
Abraham V, Mayhall NJ. Cluster many-body expansion: A many-body expansion of the electron correlation energy about a cluster mean field reference. J Chem Phys 2021; 155:054101. [PMID: 34364343 DOI: 10.1063/5.0057752] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. In the past, applications of MBE to correlation energy have been unfeasible for large systems, but recent improvements to computing resources have sparked renewed interest in capturing the correlation energy using the generalized nth order Bethe-Goldstone equation. In this work, we extend this approach, originally proposed for a Slater determinant, to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean field reference TPS configuration and includes the correlation contribution of the excited TPSs using the MBE. This method, named cluster MBE (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from a RHF reference. We present numerical results for strongly correlated systems, such as the one- and two-dimensional Hubbard models and the chromium dimer. The performance of the cMBE method is also tested by partitioning the extended π space of several large π-conjugated systems, including a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 1066 determinants for the exact FCI solution.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24060, USA
| | | |
Collapse
|
25
|
Li Manni G, Dobrautz W, Bogdanov NA, Guther K, Alavi A. Resolution of Low-Energy States in Spin-Exchange Transition-Metal Clusters: Case Study of Singlet States in [Fe(III) 4S 4] Cubanes. J Phys Chem A 2021; 125:4727-4740. [PMID: 34048648 PMCID: PMC8201447 DOI: 10.1021/acs.jpca.1c00397] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
![]()
Polynuclear transition-metal
(PNTM) clusters owe their catalytic
activity to numerous energetically low-lying spin states and stable
oxidation states. The characterization of their electronic structure
represents one of the greatest challenges of modern chemistry. We
propose a theoretical framework that enables the resolution of targeted
electronic states with ease and apply it to two [Fe(III)4S4] cubanes. Through direct access to their many-body
wave functions, we identify important correlation mechanisms and their
interplay with the geometrical distortions observed in these clusters,
which are core properties in understanding their catalytic activity.
The simulated magnetic coupling constants predicted by our strategy
allow us to make qualitative connections between spin interactions
and geometrical distortions, demonstrating its predictive power. Moreover,
despite its simplicity, the strategy provides magnetic coupling constants
in good agreement with the available experimental ones. The complexes
are intrinsically frustrated anti-ferromagnets, and the obtained spin
structures together with the geometrical distortions represent two
possible ways to release spin frustration (spin-driven Jahn–Teller
distortion). Our paradigm provides a simple, yet rigorous, route to
uncover the electronic structure of PNTM clusters and may be applied
to a wide variety of such clusters.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Department of Electronic Structure Theory, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Werner Dobrautz
- Department of Electronic Structure Theory, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Nikolay A Bogdanov
- Department of Electronic Structure Theory, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Kai Guther
- Department of Electronic Structure Theory, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Department of Electronic Structure Theory, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany.,Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
26
|
Meitei OR, Mayhall NJ. Spin-Flip Pair-Density Functional Theory: A Practical Approach To Treat Static and Dynamical Correlations in Large Molecules. J Chem Theory Comput 2021; 17:2906-2916. [PMID: 33861603 DOI: 10.1021/acs.jctc.1c00121] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a practical approach to treat static and dynamical correlation accurately in large multiconfigurational systems. The static correlation is taken into account by using the spin-flip approach, which is well-known for capturing static correlation accurately at low-computational expense. Unlike previous approaches to add dynamical correlation to spin-flip models which use perturbation theory or coupled-cluster theory, we explore the ability to use the on-top pair-density functional theory approaches recently developed by Gagliardi and co-workers (J. Comput. Theor. Chem., 2014, 10, 3669). External relaxations are performed in the spin-flip calculations through a restricted active space framework for which a truncation scheme for the orbitals used in the external excitation is presented. The performance of the approach is demonstrated by computing energy gaps between ground and excited states for diradicals, triradicals, and linear polyacene chains ranging from naphthalene to dodecacene. Accurate results are obtained using the new approach for these challenging open-shell molecular systems.
Collapse
Affiliation(s)
- Oinam Romesh Meitei
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
27
|
Dutta R, Chen GP, Henderson TM, Scuseria GE. Construction of linearly independent non-orthogonal AGP states. J Chem Phys 2021; 154:114112. [DOI: 10.1063/5.0045006] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rishab Dutta
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
| | - Guo P. Chen
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
| | - Thomas M. Henderson
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| |
Collapse
|
28
|
Ghosh SK, Rano M, Ghosh D. Configuration interaction trained by neural networks: Application to model polyaromatic hydrocarbons. J Chem Phys 2021; 154:094117. [PMID: 33685176 DOI: 10.1063/5.0040785] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The main bottleneck of a stochastic or deterministic configuration interaction method is determining the relative weights or importance of each determinant or configuration, which requires large scale matrix diagonalization. Therefore, these methods can be improved significantly from a computational standpoint if the relative importance of each configuration in the ground and excited states of molecular/model systems can be learned using machine learning techniques such as artificial neural networks (ANNs). We have used neural networks to train the configuration interaction coefficients obtained from full configuration interaction and Monte Carlo configuration interaction methods and have tested different input descriptors and outputs to find the more efficient training techniques. These ANNs have been used to calculate the ground states of one- and two-dimensional Heisenberg spin chains along with Heisenberg ladder systems, which are good approximations of polyaromatic hydrocarbons. We find excellent efficiency of training and the model this trained was used to calculate the variational ground state energies of the systems.
Collapse
Affiliation(s)
- Sumanta K Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Madhumita Rano
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
29
|
Mniszewski SM, Dub PA, Tretiak S, Anisimov PM, Zhang Y, Negre CFA. Reduction of the molecular hamiltonian matrix using quantum community detection. Sci Rep 2021; 11:4099. [PMID: 33602988 PMCID: PMC7892829 DOI: 10.1038/s41598-021-83561-x] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Accepted: 02/04/2021] [Indexed: 11/09/2022] Open
Abstract
Quantum chemistry is interested in calculating ground and excited states of molecular systems by solving the electronic Schrödinger equation. The exact numerical solution of this equation, frequently represented as an eigenvalue problem, remains unfeasible for most molecules and requires approximate methods. In this paper we introduce the use of Quantum Community Detection performed using the D-Wave quantum annealer to reduce the molecular Hamiltonian matrix in Slater determinant basis without chemical knowledge. Given a molecule represented by a matrix of Slater determinants, the connectivity between Slater determinants (as off-diagonal elements) is viewed as a graph adjacency matrix for determining multiple communities based on modularity maximization. A gauge metric based on perturbation theory is used to determine the lowest energy cluster. This cluster or sub-matrix of Slater determinants is used to calculate approximate ground state and excited state energies within chemical accuracy. The details of this method are described along with demonstrating its performance across multiple molecules of interest and bond dissociation cases. These examples provide proof-of-principle results for approximate solution of the electronic structure problem using quantum computing. This approach is general and shows potential to reduce the computational complexity of post-Hartree-Fock methods as future advances in quantum hardware become available.
Collapse
Affiliation(s)
- Susan M Mniszewski
- Computer, Computational and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Pavel A Dub
- Chemistry Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Petr M Anisimov
- Accelerator Operations and Technology Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Christian F A Negre
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| |
Collapse
|
30
|
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
|