1
|
Greiner J, Gauss J, Eriksen JJ. Error Control and Automatic Detection of Reference Active Spaces in Many-Body Expanded Full Configuration Interaction. J Phys Chem A 2024. [PMID: 39099303 DOI: 10.1021/acs.jpca.4c04056] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
We present a wide-reaching revamp of the generalized many-body expanded full configuration interaction (MBE-FCI) method. First, we outline how to automatize the selection of reference active spaces, whereby the inherent bias introduced through a manual identification is reduced, also within the context of traditional complete active space methods. Second, we allow for the use of compact orbital clusters as expansion objects, which works to circumvent the unfavorable scaling with the number of orbitals included in the space complementary to the reference orbitals. Finally, we present a new algorithm for ensuring that many-body expansions can be efficiently terminated while conservatively accounting for resulting errors. These developments are all tested on a variety of molecular systems and different orbital representations to illustrate the abilities of our algorithm to produce correlation energies within predetermined error bounds, significantly broadening the overall applicability of the MBE-FCI method.
Collapse
Affiliation(s)
- Jonas Greiner
- Department Chemie, Johannes Gutenberg-Universität Mainz Duesbergweg 10-14, 55128 Mainz, Germany
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz Duesbergweg 10-14, 55128 Mainz, Germany
| | - Janus J Eriksen
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg. 206, 2800 Kgs. Lyngby, Denmark
| |
Collapse
|
2
|
Glaser N, Baiardi A, Lieberherr AZ, Reiher M. Vibrational Entanglement through the Lens of Quantum Information Measures. J Phys Chem Lett 2024; 15:6958-6965. [PMID: 38940632 PMCID: PMC11247497 DOI: 10.1021/acs.jpclett.4c01298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024]
Abstract
We introduce a quantum information analysis of vibrational wave functions to understand complex vibrational spectra of molecules with strong anharmonic couplings and vibrational resonances. For this purpose, we define one- and two-modal entropies to guide the identification of strongly coupled vibrational modes and to characterize correlations within modal basis sets. We evaluate these descriptors for multiconfigurational vibrational wave functions which we calculate with the n-mode vibrational density matrix renormalization group algorithm. Based on the quantum information measures, we present a vibrational entanglement analysis of the vibrational ground and excited states of CO2, which display strong anharmonic effects due to the symmetry-induced and accidental (near-) degeneracies. We investigate the entanglement signature of the Fermi resonance and discuss the maximally entangled state arising from the two degenerate bending modes.
Collapse
Affiliation(s)
- Nina Glaser
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Alberto Baiardi
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Annina Z Lieberherr
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
3
|
Glaser N, Baiardi A, Reiher M. Flexible DMRG-Based Framework for Anharmonic Vibrational Calculations. J Chem Theory Comput 2023; 19:9329-9343. [PMID: 38060309 PMCID: PMC10753801 DOI: 10.1021/acs.jctc.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 10/16/2023] [Accepted: 10/23/2023] [Indexed: 12/08/2023]
Abstract
We present a novel formulation of the vibrational density matrix renormalization group (vDMRG) algorithm tailored to strongly anharmonic molecules described by general, high-dimensional model representations of potential energy surfaces. For this purpose, we extend the vDMRG framework to support vibrational Hamiltonians expressed in the so-called n-mode second-quantization formalism. The resulting n-mode vDMRG method offers full flexibility with respect to both the functional form of the PES and the choice of the single-particle basis set. We leverage this framework to apply, for the first time, vDMRG based on an anharmonic modal basis set optimized with the vibrational self-consistent field algorithm on an on-the-fly constructed PES. We also extend the n-mode vDMRG framework to include excited-state-targeting algorithms in order to efficiently calculate anharmonic transition frequencies. We demonstrate the capabilities of our novel n-mode vDMRG framework for methyloxirane, a challenging molecule with 24 coupled vibrational modes.
Collapse
Affiliation(s)
- Nina Glaser
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Alberto Baiardi
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Department of Chemistry
and Applied Biosciences, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
4
|
Xu Y, Cheng Y, Song Y, Ma H. New Density Matrix Renormalization Group Approaches for Strongly Correlated Systems Coupled with Large Environments. J Chem Theory Comput 2023. [PMID: 37471519 DOI: 10.1021/acs.jctc.2c01316] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
Thanks to the high compression of the matrix product state (MPS) form of the wave function and the efficient site-by-site iterative sweeping optimization algorithm, the density matrix normalization group (DMRG) and its time-dependent variant (TD-DMRG) have been established as powerful computational tools in accurately simulating the electronic structure and quantum dynamics of strongly correlated molecules with a large number (101-2) of quantum degrees of freedom (active orbitals or vibrational modes). However, the quantitative characterization of the quantum many-body behaviors of realistic strongly correlated systems requires a further consideration of the interaction between the embedded active subsystem and the remaining correlated environment, e.g., a larger number (102-3) of external orbitals in electronic structure or infinite condensed-phase phononic modes in nucleus dynamics. To this end, we introduced three new post-DMRG and TD-DMRG approaches, namely (1) DMRG2sCI-MRCI and DMRG2sCI-ENPT by the reconstruction of selected configuration interaction (sCI) type of compact reference function from DMRG coefficients and the use of externally contracted MRCI (multireference configuration interaction) and Epstein-Nesbet perturbation theory (ENPT), without recourse to the expensive high order n-electron reduced density matrices (n-RDMs). (2) DMRG combined with RR-MRCI (renormalized residue-based MRCI), which improves the computational accuracy and efficiency of internally contracted (ic) MRCI by renormalizing the contracted bases with small-sized buffer environment(s) of a few external orbitals as probes based on quantum information theory. (3) HM (hierarchical mapping)-TD-DMRG in which a large environment is reduced to a small number of renormalized environmental modes (which accounts for the most vital system-environment interactions) through stepwise mapping transformation. These advances extend the efficacy of highly accurate DMRG/TD-DMRG computations to the quantitative characterization of the electronic structure and quantum dynamics in realistic strongly correlated systems coupled with large environments and are reviewed in this paper.
Collapse
Affiliation(s)
- Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yinxuan Song
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| |
Collapse
|
5
|
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
|
6
|
Abstract
We introduce DMRG[FEAST], a new method for optimizing excited-state many-body wave functions with the density matrix renormalization group (DMRG) algorithm. Our approach applies the FEAST algorithm, originally designed for large-scale diagonalization problems, to matrix product state wave functions. We show that DMRG[FEAST] enables the stable optimization of both low- and high-energy eigenstates, therefore overcoming the limitations of state-of-the-art excited-state DMRG algorithms. We demonstrate the reliability of DMRG[FEAST] by calculating anharmonic vibrational excitation energies of molecules with up to 30 fully coupled degrees of freedom.
Collapse
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Anna Klára Kelemen
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
7
|
Ren JJ, Wang YH, Li WT, Jiang T, Shuai ZG. Time-dependent density matrix renormalization group coupled with n-mode representation potentials for the excited state radiationless decay rate: Formalism and application to azulene. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2108138] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jia-jun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Yuan-heng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Wei-tang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Zhi-gang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| |
Collapse
|
8
|
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
|
9
|
Wang T, Ma Y, Zhao L, Jiang J. Portably parallel construction of a configuration-interaction wave function from a matrix-product state using the Charm++ framework. J Comput Chem 2020; 41:2707-2721. [PMID: 32986283 DOI: 10.1002/jcc.26424] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2020] [Revised: 08/16/2020] [Accepted: 09/02/2020] [Indexed: 11/10/2022]
Abstract
The construction of configuration-interaction (CI) expansions from a matrix product state (MPS) involves numerous matrix operations and the skillful sampling of important configurations in a large Hilbert space. In this work, we present an efficient procedure for constructing CI expansions from MPS employing the parallel object-oriented Charm++ programming framework, upon which automatic load-balancing and object migrating facilities can be employed. This procedure was employed in the MPS-to-CI utility (Moritz et al., J. Chem. Phys. 2007, 126, 224109), the sampling-reconstructed complete active-space algorithm (SR-CAS, Boguslawski et al., J. Chem. Phys. 2011, 134, 224101), and the entanglement-driven genetic algorithm (EDGA, Luo et al., J. Chem. Theory Comput. 2017, 13, 4699). It enhances productivity and allows the sampling programs to evolve to their population-expansion versions, for example, EDGA with population expansion (PE-EDGA). Further, examples of 1,2-dioxetanone and firefly dioxetanone anion (FDO- ) molecules demonstrated the following: (a) parallel efficiencies can be persistently improved by simply by increasing the proportions of the asynchronous executions and (b) a sampled CAS-type CI wave function of a bi-radical-state FDO- molecule utilizing the full valence (30e,26o) active space can be constructed within a few hours with using thousands of cores.
Collapse
Affiliation(s)
- Ting Wang
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China.,Center of Scientific Computing Applications and Research, Chinese Academy of Sciences, Beijing, China.,School of Computer Science, Shaanxi Normal University, Xi'an, China
| | - Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China.,Center of Scientific Computing Applications and Research, Chinese Academy of Sciences, Beijing, China
| | - Lian Zhao
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China.,Center of Scientific Computing Applications and Research, Chinese Academy of Sciences, Beijing, China
| | - Jinrong Jiang
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China.,Center of Scientific Computing Applications and Research, Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
10
|
Pierloot K, Phung QM, Ghosh A. Electronic Structure of Neutral and Anionic Iron–Nitrosyl Corrole. A Multiconfigurational and Density Matrix Renormalization Group Investigation. Inorg Chem 2020; 59:11493-11502. [DOI: 10.1021/acs.inorgchem.0c01312] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Quan Manh Phung
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Chikusa, Nagoya 464-8602, Japan
| | - Abhik Ghosh
- Department of Chemistry, UiT-The Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
11
|
Lang J, Antalík A, Veis L, Brandejs J, Brabec J, Legeza Ö, Pittner J. Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T). J Chem Theory Comput 2020; 16:3028-3040. [PMID: 32275424 DOI: 10.1021/acs.jctc.0c00065] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Sciences, Charles University, Albertov 6, 128 00 Prague 2, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
12
|
Brandejs J, Višňák J, Veis L, Maté M, Legeza Ö, Pittner J. Toward DMRG-tailored coupled cluster method in the 4c-relativistic domain. J Chem Phys 2020; 152:174107. [DOI: 10.1063/1.5144974] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Jakub Višňák
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
- Czech Academic City in Erbil, Yassin Najar Street, Kurani Ankawa, Erbil, Kurdistan, Region of Iraq
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Mihály Maté
- Strongly Correlated Systems “Lendület” Research Group, Institute for Solid State Physics and Optics, MTA Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 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, Institute for Solid State Physics and Optics, MTA Wigner Research Centre for Physics, Konkoly-Thege Miklós út 29-33, H-1121 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
13
|
Stein CJ, Reiher M. Semiclassical Dispersion Corrections Efficiently Improve Multiconfigurational Theory with Short-Range Density-Functional Dynamic Correlation. J Phys Chem A 2020; 124:2834-2841. [PMID: 32186877 DOI: 10.1021/acs.jpca.0c02130] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Multiconfigurational wave functions are known to describe the electronic structure across a Born-Oppenheimer surface qualitatively correct. However, for quantitative reaction energies, dynamic correlation originating from the many configurations involving excitations out of the restricted orbital space, the active space, must be considered. Standard procedures involve approximations that eventually limit the ultimate accuracy achievable (most prominently, multireference perturbation theory). At the same time, the computational cost increases dramatically due to the necessity to obtain higher-order reduced density matrices. It is this disproportion that leads us here to propose an MC-srDFT-D hybrid approach of semiclassical dispersion (D) corrections to cover long-range dynamic correlation in a multiconfigurational (MC) wave function theory, which includes short-range (sr) dynamic correlation by density functional theory (DFT) without double counting. We demonstrate that the reliability of this approach is very good (at negligible cost), especially when considering that standard second-order multireference perturbation theory usually overestimates dispersion interactions.
Collapse
Affiliation(s)
- Christopher J Stein
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
14
|
Antalík A, Veis L, Brabec J, Demel O, Legeza Ö, Pittner J. Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen). J Chem Phys 2019; 151:084112. [PMID: 31470730 DOI: 10.1063/1.5110477] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.
Collapse
Affiliation(s)
- Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research Group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
15
|
Baiardi A, Stein CJ, Barone V, Reiher M. Optimization of highly excited matrix product states with an application to vibrational spectroscopy. J Chem Phys 2019; 150:094113. [DOI: 10.1063/1.5068747] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alberto Baiardi
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Christopher J. Stein
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
16
|
Luo Z, Ma Y, Wang X, Ma H. Externally-Contracted Multireference Configuration Interaction Method Using a DMRG Reference Wave Function. J Chem Theory Comput 2018; 14:4747-4755. [PMID: 30052433 DOI: 10.1021/acs.jctc.8b00613] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The recent development of the density matrix renormalization group (DMRG) method in multireference quantum chemistry makes it practical to evaluate static correlation in a large active space, while dynamic correlation provides a critical correction to the DMRG reference for strong-correlated systems and is usually obtained using multireference perturbation (MRPT) or configuration interaction (MRCI) methods with internal contraction (ic) approximation. These methods can use an active space scalable to relatively larger size references than has previously been possible. However, they are still hardly applicable to systems with an active space larger than 30 orbitals and/or a large basis set because of high computation and storage costs of high-order reduced density matrices (RDMs) and the crucial dependence of the MRCI Hamiltonian dimension on the number of virtual orbitals. In this work, we propose a new effective implementation of DMRG-MRCI, in which we use reconstructed CASCI-type configurations from DMRG wave function via the entropy-driving genetic algorithm (EDGA) [ Luo et al. J. Chem. Theory Comput. 2017 , 13 , 4699 - 4710 . ] and integrate it with MRCI by an external contraction (ec) scheme. This bypasses the bottleneck of computing high-order RDMs in traditional DMRG dynamic correlation methods with ic approximation, and the number of MRCI configurations is not dependent on the number of virtual orbitals. Therefore, the DMRG-ec-MRCI method is promising for dealing with a larger active space than 30 orbitals and large basis sets. We demonstrate the capability of our DMRG-ec-MRCI method in several benchmark applications, including the evaluation of the potential energy curve of Cr2, single-triplet gaps of higher n-acene molecules, and the energy of the Eu-BTBP(NO3)3 complex.
Collapse
Affiliation(s)
- Zhen Luo
- School of Chemistry and Chemical Engineering , Nanjing University , Nanjing 210023 , China
| | - Yingjin Ma
- Department of High Performance Computing Technology and Application Development , Computer Network Information Center, Chinese Academy of Sciences , Beijing 100190 , China.,Center of Scientific Computing Applications & Research , Chinese Academy of Sciences , Beijing 100190 , China
| | - Xicun Wang
- Nanjing Foreign Language School, Nanjing 210008 , China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering , Nanjing University , Nanjing 210023 , China
| |
Collapse
|
17
|
Luo Z, Ma Y, Liu C, Ma H. Efficient Reconstruction of CAS-CI-Type Wave Functions for a DMRG State Using Quantum Information Theory and a Genetic Algorithm. J Chem Theory Comput 2017; 13:4699-4710. [DOI: 10.1021/acs.jctc.7b00439] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Zhen Luo
- Key
Laboratory of Mesoscopic Chemistry of MOE, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Yingjin Ma
- Department
of High Performance Computing Technology and Application Development,
Computer Network Information Center, Chinese Academy of Sciences, Beijing 100190, China
- Center of Scientific Computing Applications & Research, Chinese Academy of Sciences, Beijing 100190, China
| | - Chungen Liu
- Key
Laboratory of Mesoscopic Chemistry of MOE, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Key
Laboratory of Mesoscopic Chemistry of MOE, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
18
|
Baiardi A, Stein CJ, Barone V, Reiher M. Vibrational Density Matrix Renormalization Group. J Chem Theory Comput 2017; 13:3764-3777. [DOI: 10.1021/acs.jctc.7b00329] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alberto Baiardi
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Christopher J. Stein
- Laboratorium
für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Vincenzo Barone
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Markus Reiher
- Laboratorium
für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
19
|
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Zürich, Switzerland
| |
Collapse
|
20
|
Veis L, Antalík A, Brabec J, Neese F, Legeza Ö, Pittner J. Coupled Cluster Method with Single and Double Excitations Tailored by Matrix Product State Wave Functions. J Phys Chem Lett 2016; 7:4072-4078. [PMID: 27682626 DOI: 10.1021/acs.jpclett.6b01908] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
In the past decade, the quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as the method of choice for calculations of strongly correlated molecular systems. Despite its favorable scaling, it is in practice not suitable for computations of dynamic correlation. We present a novel method for accurate "post-DMRG" treatment of dynamic correlation based on the tailored coupled cluster (CC) theory in which the DMRG method is responsible for the proper description of nondynamic correlation, whereas dynamic correlation is incorporated through the framework of the CC theory. We illustrate the potential of this method on prominent multireference systems, in particular, N2 and Cr2 molecules and also oxo-Mn(Salen), for which we have performed the first post-DMRG computations in order to shed light on the energy ordering of the lowest spin states.
Collapse
Affiliation(s)
- Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic , v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic , v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic , v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Frank Neese
- Max Planck Institut für Chemische Energiekonversion , Stiftstrasse 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research Group, Wigner Research Centre for Physics , H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic , v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
21
|
Sivalingam K, Krupicka M, Auer AA, Neese F. Comparison of fully internally and strongly contracted multireference configuration interaction procedures. J Chem Phys 2016; 145:054104. [DOI: 10.1063/1.4959029] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kantharuban Sivalingam
- Max-Planck Institute of Chemical Energy Conversion, Stiftstrasse 34, 45470 Mülheim an der Ruhr, Germany
| | - Martin Krupicka
- Max-Planck Institute of Chemical Energy Conversion, Stiftstrasse 34, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A. Auer
- Max-Planck Institute of Chemical Energy Conversion, Stiftstrasse 34, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck Institute of Chemical Energy Conversion, Stiftstrasse 34, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
22
|
Keller S, Dolfi M, Troyer M, Reiher M. An efficient matrix product operator representation of the quantum chemical Hamiltonian. J Chem Phys 2015; 143:244118. [DOI: 10.1063/1.4939000] [Citation(s) in RCA: 106] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
|
23
|
Keller S, Boguslawski K, Janowski T, Reiher M, Pulay P. Selection of active spaces for multiconfigurational wavefunctions. J Chem Phys 2015; 142:244104. [DOI: 10.1063/1.4922352] [Citation(s) in RCA: 74] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sebastian Keller
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Katharina Boguslawski
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Tomasz Janowski
- Department of Chemistry and Biochemistry, Fulbright College of Arts and Sciences, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Peter Pulay
- Department of Chemistry and Biochemistry, Fulbright College of Arts and Sciences, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
24
|
Chattopadhyay S, Chaudhuri RK, Mahapatra US. State-specific multireference perturbation theory with improved virtual orbitals: Taming the ground state of F2, Be2,and N2. J Comput Chem 2015; 36:907-25. [DOI: 10.1002/jcc.23873] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2014] [Revised: 01/09/2015] [Accepted: 01/11/2015] [Indexed: 11/07/2022]
Affiliation(s)
- Sudip Chattopadhyay
- Department of Chemistry; Indian Institute of Engineering Science and Technology; Shibpur Howrah 711103 India
| | | | - Uttam Sinha Mahapatra
- Department of Physics; Maulana Azad College; 8 Rafi Ahmed Kidwai Road Kolkata 700013 India
| |
Collapse
|
25
|
Olivares-Amaya R, Hu W, Nakatani N, Sharma S, Yang J, Chan GKL. Theab-initiodensity matrix renormalization group in practice. J Chem Phys 2015; 142:034102. [DOI: 10.1063/1.4905329] [Citation(s) in RCA: 231] [Impact Index Per Article: 25.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
| | - Weifeng Hu
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Naoki Nakatani
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
- Catalysis Research Center, Hokkaido University, Kita 21 Nishi 10, Sapporo, Hokkaido 001-0021, Japan
| | - Sandeep Sharma
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Jun Yang
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Garnet Kin-Lic Chan
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
26
|
Kurashige Y, Chalupský J, Lan TN, Yanai T. Complete active space second-order perturbation theory with cumulant approximation for extended active-space wavefunction from density matrix renormalization group. J Chem Phys 2014; 141:174111. [DOI: 10.1063/1.4900878] [Citation(s) in RCA: 100] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yuki Kurashige
- Department of Theoretical and Computational Molecular Science,
Institute for Molecular Science, Okazaki, Aichi 444-8585,
Japan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
- Japan Science and Technology Agency,
PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Jakub Chalupský
- Department of Theoretical and Computational Molecular Science,
Institute for Molecular Science, Okazaki, Aichi 444-8585,
Japan
| | - Tran Nguyen Lan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Takeshi Yanai
- Department of Theoretical and Computational Molecular Science,
Institute for Molecular Science, Okazaki, Aichi 444-8585,
Japan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
27
|
Tecmer P, Boguslawski K, Johnson PA, Limacher PA, Chan M, Verstraelen T, Ayers PW. Assessing the Accuracy of New Geminal-Based Approaches. J Phys Chem A 2014; 118:9058-68. [PMID: 24745368 DOI: 10.1021/jp502127v] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Paweł Tecmer
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, L8S 4M1, Hamilton, Ontario, Canada
| | - Katharina Boguslawski
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, L8S 4M1, Hamilton, Ontario, Canada
| | - Paul A. Johnson
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, L8S 4M1, Hamilton, Ontario, Canada
| | - Peter A. Limacher
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, L8S 4M1, Hamilton, Ontario, Canada
| | - Matthew Chan
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, L8S 4M1, Hamilton, Ontario, Canada
| | - Toon Verstraelen
- Center
for Molecular Modeling, QCMM Alliance Ghent-Brussels, Ghent University, Technologiepark
903, Zwijnaarde 9052, Belgium
| | - Paul W. Ayers
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, L8S 4M1, Hamilton, Ontario, Canada
| |
Collapse
|
28
|
Kelly TP, Perera A, Bartlett RJ, Greer JC. Monte Carlo configuration interaction with perturbation corrections for dissociation energies of first row diatomic molecules: C2, N2, O2, CO, and NO. J Chem Phys 2014; 140:084114. [DOI: 10.1063/1.4866609] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
Ma Y, Ma H. Assessment of various natural orbitals as the basis of large active space density-matrix renormalization group calculations. J Chem Phys 2013; 138:224105. [DOI: 10.1063/1.4809682] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Sousa C, Tosoni S, Illas F. Theoretical Approaches to Excited-State-Related Phenomena in Oxide Surfaces. Chem Rev 2012. [DOI: 10.1021/cr300228z] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Carmen Sousa
- Departament de Química
Física and Institut de Química Teòrica i Computacional
(IQTCUB), Universitat de Barcelona, C/Martí
i Franquès 1, 08028 Barcelona, Spain
| | - Sergio Tosoni
- Departament de Química
Física and Institut de Química Teòrica i Computacional
(IQTCUB), Universitat de Barcelona, C/Martí
i Franquès 1, 08028 Barcelona, Spain
- Departamento de Química, Universidad de Las Palmas de Gran Canaria, Campus Universitario
de Tafira, 35017 Las Palmas de Gran Canaria, Spain
| | - Francesc Illas
- Departament de Química
Física and Institut de Química Teòrica i Computacional
(IQTCUB), Universitat de Barcelona, C/Martí
i Franquès 1, 08028 Barcelona, Spain
| |
Collapse
|
31
|
Coe JP, Taylor DJ, Paterson MJ. Calculations of potential energy surfaces using Monte Carlo configuration interaction. J Chem Phys 2012. [DOI: 10.1063/1.4767052] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
32
|
Boguslawski K, Tecmer P, Legeza Ö, Reiher M. Entanglement Measures for Single- and Multireference Correlation Effects. J Phys Chem Lett 2012; 3:3129-3135. [PMID: 26296018 DOI: 10.1021/jz301319v] [Citation(s) in RCA: 103] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Electron correlation effects are essential for an accurate ab initio description of molecules. A quantitative a priori knowledge of the single- or multireference nature of electronic structures as well as of the dominant contributions to the correlation energy can facilitate the decision regarding the optimum quantum chemical method of choice. We propose concepts from quantum information theory as orbital entanglement measures that allow us to evaluate the single- and multireference character of any molecular structure in a given orbital basis set. By studying these measures we can detect possible artifacts of small active spaces.
Collapse
Affiliation(s)
- Katharina Boguslawski
- †ETH Zurich, Laboratory of Physical Chemistry, Wolfgang-Pauli-Str. 10, CH-8093 Zurich, Switzerland
| | - Pawel̷ Tecmer
- †ETH Zurich, Laboratory of Physical Chemistry, Wolfgang-Pauli-Str. 10, CH-8093 Zurich, Switzerland
| | - Örs Legeza
- ‡MTA-WRCP Strongly Correlated Systems "Lendület" Research Group, H-1525, Budapest, Hungary
| | - Markus Reiher
- †ETH Zurich, Laboratory of Physical Chemistry, Wolfgang-Pauli-Str. 10, CH-8093 Zurich, Switzerland
| |
Collapse
|
33
|
Zimmerman PM, Bell F, Goldey M, Bell AT, Head-Gordon M. Restricted active space spin-flip configuration interaction: Theory and examples for multiple spin flips with odd numbers of electrons. J Chem Phys 2012; 137:164110. [DOI: 10.1063/1.4759076] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
|
34
|
Boguslawski K, Marti KH, Legeza O, Reiher M. Accurate ab Initio Spin Densities. J Chem Theory Comput 2012; 8:1970-1982. [PMID: 22707921 PMCID: PMC3373140 DOI: 10.1021/ct300211j] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2012] [Indexed: 01/12/2023]
Abstract
We present an approach for the calculation of spin density
distributions
for molecules that require very large active spaces for a qualitatively
correct description of their electronic structure. Our approach is
based on the density-matrix renormalization group (DMRG) algorithm
to calculate the spin density matrix elements as a basic quantity
for the spatially resolved spin density distribution. The spin density
matrix elements are directly determined from the second-quantized
elementary operators optimized by the DMRG algorithm. As an analytic
convergence criterion for the spin density distribution, we employ
our recently developed sampling-reconstruction scheme [J.
Chem. Phys.2011, 134, 224101]
to build an accurate complete-active-space configuration-interaction
(CASCI) wave function from the optimized matrix product states. The
spin density matrix elements can then also be determined as an expectation
value employing the reconstructed wave function expansion. Furthermore,
the explicit reconstruction of a CASCI-type wave function provides
insight into chemically interesting features of the molecule under
study such as the distribution of α and β electrons in
terms of Slater determinants, CI coefficients, and natural orbitals.
The methodology is applied to an iron nitrosyl complex which we have
identified as a challenging system for standard approaches [J. Chem. Theory Comput.2011, 7, 2740].
Collapse
|