1
|
Song M, Alavi A, Li Manni G. Permutation symmetry in spin-adapted many-body wave functions. Faraday Discuss 2024. [PMID: 39158096 DOI: 10.1039/d4fd00061g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/20/2024]
Abstract
In the domain of exchange-coupled polynuclear transition-metal (PNTM) clusters, local emergent symmetries exist which can be exploited to greatly increase the sparsity of the configuration interaction (CI) eigensolutions of such systems. Sparsity of the CI secular problem is revealed by exploring the site permutation space within spin-adapted many-body bases, and highly compressed wave functions may arise by finding optimal site orderings. However, the factorial cost of searching through the permutation space remains a bottleneck for clusters with a large number of metal centers. In this work, we explore ways to reduce the factorial scaling, by combining permutation and point group symmetry arguments, and using commutation relations between cumulative partial spin and the Hamiltonian operators, . Certain site orderings lead to commuting operators, from which more sparse wave functions arise. Two graphical strategies will be discussed, one to rapidly evaluate the commutators of interest, and one in the form of a tree search algorithm to predict how many and which distinct site permutations are to be analyzed, eliminating redundancies in the permutation space. Particularly interesting is the case of the singlet spin states for which an additional reversal symmetry can be utilized to further reduce the number of distinct site permutations.
Collapse
Affiliation(s)
- Maru Song
- Electronic Structure Theory Department, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany.
| | - Ali Alavi
- Electronic Structure Theory Department, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany.
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| | - Giovanni Li Manni
- Electronic Structure Theory Department, Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany.
| |
Collapse
|
2
|
Greiner J, Gianni I, Nottoli T, Lipparini F, Eriksen JJ, Gauss J. MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces. J Chem Theory Comput 2024; 20:4663-4675. [PMID: 38809011 DOI: 10.1021/acs.jctc.4c00388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals.
Collapse
Affiliation(s)
- Jonas Greiner
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz 55128, Germany
| | - Ivan Gianni
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Tommaso Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Janus J Eriksen
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg. 206, Kgs. Lyngby 2800, Denmark
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, Mainz 55128, Germany
| |
Collapse
|
3
|
Cheng Y, Ma H. Renormalized-Residue-Based Multireference Configuration Interaction Method for Strongly Correlated Systems. J Chem Theory Comput 2024; 20:1988-2009. [PMID: 38380619 DOI: 10.1021/acs.jctc.3c01247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2024]
Abstract
The implementation of multireference configuration interaction (MRCI) methods in quantum systems with large active spaces is hindered by the expansion of configuration bases or the intricate handling of reduced density matrices (RDMs). In this work, we present a spin-adapted renormalized-residue-based MRCI (RR-MRCI) approach that leverages renormalized residues to effectively capture the entanglement between active and inactive orbitals. This approach is reinforced by a novel efficient algorithm, which also facilitates an efficient deployment of spin-adapted matrix product state MRCI (MPS-MRCI). The RR-MRCI framework possesses several advantages: (1) It considers the orbital entanglement and utilizes highly compressed MPS structure, improving computational accuracy and efficiency compared with internally contracted (ic) MRCI. (2) Utilizing small-sized buffer environments of a few external orbitals as probes based on quantum information theory, it enhances computational efficiency over MPS-MRCI and offers potential application to large molecular systems. (3) The RR framework can be implemented in conjunction with ic-MRCI, eliminating the need for high-rank RDMs, by using distinct renormalized residues. We evaluated this method across nine diverse molecular systems, including Cu2O22+ with an active space of (24e,24o) and two complexes of lanthanide and actinide with active space (38e,36o), demonstrating the method's versatility and efficacy.
Collapse
Affiliation(s)
- Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, School of Chemistry and Chemical Engineering, Shandong University, Qingdao, Shandong 266237, China
| |
Collapse
|
4
|
Crisci L, Coppola F, Petrone A, Rega N. Tuning ultrafast time-evolution of photo-induced charge-transfer states: A real-time electronic dynamics study in substituted indenotetracene derivatives. J Comput Chem 2024; 45:210-221. [PMID: 37706600 DOI: 10.1002/jcc.27231] [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: 07/20/2023] [Revised: 08/31/2023] [Accepted: 09/05/2023] [Indexed: 09/15/2023]
Abstract
Photo-induced charge transfer (CT) states are pivotal in many technological and biological processes. A deeper knowledge of such states is mandatory for modeling the charge migration dynamics. Real-time time-dependent density functional theory (RT-TD-DFT) electronic dynamics simulations are employed to explicitly observe the electronic density time-evolution upon photo-excitation. Asymmetrically substituted indenotetracene molecules, given their potential application as n-type semiconductors in organic photovoltaic materials, are here investigated. Effects of substituents with different electron-donating characters are analyzed in terms of the overall electronic energy spacing and resulting ultrafast CT dynamics through linear response (LR-)TD-DFT and RT-TD-DFT based approaches. The combination of the computational techniques here employed provided direct access to the electronic density reorganization in time and to its spatial and rational representation in terms of molecular orbital occupation time evolution. Such results can be exploited to design peculiar directional charge dynamics, crucial when photoactive materials are used for light-harvesting applications.
Collapse
Affiliation(s)
- Luigi Crisci
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Normale Superiore di Pisa, Pisa, Italy
| | | | - Alessio Petrone
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| | - Nadia Rega
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| |
Collapse
|
5
|
Safari AA, Anderson RJ, Manni GL. Toward a Stochastic Complete Active Space Second-Order Perturbation Theory. J Phys Chem A 2024; 128:281-291. [PMID: 38154124 PMCID: PMC10788896 DOI: 10.1021/acs.jpca.3c05109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 11/14/2023] [Accepted: 11/14/2023] [Indexed: 12/30/2023]
Abstract
In this work, an internally contracted stochastic complete active space second-order perturbation theory, stochastic-CASPT2, is reported. The method relies on stochastically sampled reduced density matrices (RDMs) up to rank four and contractions thereof with the generalized Fock matrix. A new protocol for calculating higher-order RDMs in full configuration interaction quantum Monte Carlo (FCIQMC) has been designed based on (1) restricting sampling of the corresponding excitations to a deterministic subspace, (2) averaging the RDMs from independent dynamics and (3) projecting them onto the closest positive semi-definite matrix. Our protocol avoids previously encountered numerical conditioning problems in the orthogonalization of the perturber overlap matrix stemming from numerical noise. The chromium dimer CASSCF(12,12)/CASPT2 binding curve is computed as a proof of concept.
Collapse
Affiliation(s)
- Arta A. Safari
- Max-Planck-Institute for Solid State
Research, 70569 Stuttgart, Germany
| | | | - Giovanni Li Manni
- Max-Planck-Institute for Solid State
Research, 70569 Stuttgart, Germany
| |
Collapse
|
6
|
Zhai H, Larsson HR, Lee S, Cui ZH, Zhu T, Sun C, Peng L, Peng R, Liao K, Tölle J, Yang J, Li S, Chan GKL. Block2: A comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond. J Chem Phys 2023; 159:234801. [PMID: 38108484 DOI: 10.1063/5.0180424] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 11/16/2023] [Indexed: 12/19/2023] Open
Abstract
block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency and to support integration with external numerical packages. Here, we explain the design principles and currently supported features and present numerical examples in a range of applications.
Collapse
Affiliation(s)
- Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Henrik R Larsson
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Tianyu Zhu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chong Sun
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Linqing Peng
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ruojing Peng
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ke Liao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Johannes Tölle
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Junjie Yang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Shuoxue Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
7
|
Phung QM, Nam HN, Ghosh A. Local Oxidation States in {FeNO} 6-8 Porphyrins: Insights from DMRG/CASSCF-CASPT2 Calculations. Inorg Chem 2023. [PMID: 38010736 DOI: 10.1021/acs.inorgchem.3c03689] [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/2023]
Abstract
A first DMRG/CASSCF-CASPT2 study of a series of paradigmatic {FeNO}6, {FeNO}7, and {FeNO}8 heme-nitrosyl complexes has led to substantial new insight as well as uncovered key shortcomings of the DFT approach. By virtue of its balanced treatment of static and dynamic correlation, the calculations have provided some of the most authoritative information available to date on the energetics of low- versus high-spin states of different classes of heme-nitrosyl complexes. Thus, the calculations indicate low doublet-quartet gaps of 1-4 kcal/mol for {FeNO}7 complexes and high singlet-triplet gaps of ≳20 kcal/mol for both {FeNO}6 and {FeNO}8 complexes. In contrast, DFT calculations yield widely divergent spin state gaps as a function of the exchange-correlation functional. DMRG-CASSCF calculations also help calibrate DFT spin densities for {FeNO}7 complexes, pointing to those obtained from classic pure functionals as the most accurate. The general picture appears to be that nearly all the spin density of Fe[P](NO) is localized on the Fe, while the axial ligand imidazole (ImH) in Fe[P](NO)(ImH) pushes a part of the spin density onto the NO moiety. An analysis of the DMRG-CASSCF wave function in terms of localized orbitals and of the resulting configuration state functions in terms of resonance forms with varying NO(π*) occupancies has allowed us to address the longstanding question of local oxidation states in heme-nitrosyl complexes. The analysis indicates NO(neutral) resonance forms [i.e., Fe(II)-NO0 and Fe(III)-NO0] as the major contributors to both {FeNO}6 and {FeNO}7 complexes. This finding is at variance with the common formulation of {FeNO}6 hemes as Fe(II)-NO+ species but is consonant with an Fe L-edge XAS analysis by Solomon and co-workers. For the {FeNO}8 complex {Fe[P](NO)}-, our analysis suggests a resonance hybrid description: Fe(I)-NO0 ↔ Fe(II)-NO-, in agreement with earlier DFT studies. Vibrational analyses of the compounds studied indicate an imperfect but fair correlation between the NO stretching frequency and NO(π*) occupancy, highlighting the usefulness of vibrational data as a preliminary indicator of the NO oxidation state.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Ho Ngoc Nam
- Institute of Materials Innovation, Institutes of Innovation for Future Society, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
- Department of Chemical Systems Engineering, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8603, Japan
| | - Abhik Ghosh
- Department of Chemistry, UiT the Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
8
|
Kim SY, Park JW. Approximate Excited-State Geometry Optimization with the State-Averaged Adaptive Sampling Configuration Interaction Algorithm with Large Active Spaces. J Chem Theory Comput 2023; 19:7260-7272. [PMID: 37800852 DOI: 10.1021/acs.jctc.3c00808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
The selected configuration interaction (SCI) wave function is a useful approximation to the full configuration interaction (FCI) one. The adaptive sampling CI (ASCI) method is a deterministic SCI method. By combining ASCI and orbital optimization, the ASCI self-consistent field (ASCI-SCF) method, which is an approximation of the complete active space self-consistent field (CASSCF) method, can be formulated as well. However, their applicability has been tested mainly on the systems in their electronically ground states. In this work, we implement the state-average (SA) ansatz in ASCI-SCF calculations to calculate excited states. We also derive expressions for the approximate analytical gradient and implement them as a computer program. We demonstrate the applicability of the current method for calculating vertical and adiabatic excitation energies and optimizing the molecular geometries of thermally activated delayed fluorescence (TADF) molecules.
Collapse
Affiliation(s)
- So Yeon Kim
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| | - Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
9
|
Li S, Misiewicz JP, Evangelista FA. Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory. J Chem Phys 2023; 159:114106. [PMID: 37712785 DOI: 10.1063/5.0159403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 08/14/2023] [Indexed: 09/16/2023] Open
Abstract
Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.
Collapse
Affiliation(s)
- Shuhang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Jonathon P Misiewicz
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
10
|
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
|
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
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
13
|
Magoulas I, Shen J, Piecuch P. Addressing strong correlation by approximate coupled-pair methods with active-space and full treatments of three-body clusters. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2057365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
- Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
14
|
Benediktsson B, Bjornsson R. Analysis of the Geometric and Electronic Structure of Spin-Coupled Iron-Sulfur Dimers with Broken-Symmetry DFT: Implications for FeMoco. J Chem Theory Comput 2022; 18:1437-1457. [PMID: 35167749 PMCID: PMC8908755 DOI: 10.1021/acs.jctc.1c00753] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
![]()
The open-shell electronic
structure of iron–sulfur clusters
presents considerable challenges to quantum chemistry, with the complex
iron–molybdenum cofactor (FeMoco) of nitrogenase representing
perhaps the ultimate challenge for either wavefunction or density
functional theory. While broken-symmetry density functional theory
has seen some success in describing the electronic structure of such
cofactors, there is a large exchange–correlation functional
dependence in calculations that is not fully understood. In this work,
we present a geometric benchmarking test set, FeMoD11, of synthetic
spin-coupled Fe–Fe and Mo–Fe dimers, with relevance
to the molecular and electronic structure of the Mo-nitrogenase FeMo
cofactor. The reference data consists of high-resolution crystal structures
of metal dimer compounds in different oxidation states. Multiple density
functionals are tested on their ability to reproduce the local geometry,
specifically the Fe–Fe/Mo–Fe distance, for both antiferromagnetically
coupled and ferromagnetically coupled dimers via the broken-symmetry
approach. The metal–metal distance is revealed not only to
be highly sensitive to the amount of exact exchange in the functional
but also to the specific exchange and correlation functionals. For
the antiferromagnetically coupled dimers, the calculated metal–metal
distance correlates well with the covalency of the bridging metal–ligand
bonds, as revealed via the corresponding orbital analysis, Hirshfeld
S/Fe charges, and Fe–S Mayer bond order. Superexchange via
bridging ligands is expected to be the dominant interaction in these
dimers, and our results suggest that functionals that predict accurate
Fe–Fe and Mo–Fe distances describe the overall metal–ligand
covalency more accurately and in turn the superexchange of these systems.
The best performing density functionals of the 16 tested for the FeMoD11
test set are revealed to be either the nonhybrid functionals r2SCAN and B97-D3 or hybrid functionals with 10–15% exact
exchange: TPSSh and B3LYP*. These same four functionals are furthermore
found to reproduce the high-resolution X-ray structure of FeMoco well
according to quantum mechanics/molecular mechanics (QM/MM) calculations.
Almost all nonhybrid functionals systematically underestimate Fe–Fe
and Mo–Fe distances (with r2SCAN and B97-D3 being
the sole exceptions), while hybrid functionals with >15% exact
exchange
(including range-separated hybrid functionals) overestimate them.
The results overall suggest r2SCAN, B97-D3, TPSSh, and
B3LYP* as accurate density functionals for describing the electronic
structure of iron–sulfur clusters in general, including the
complex FeMoco cluster of nitrogenase.
Collapse
Affiliation(s)
- Bardi Benediktsson
- Science Institute, University of Iceland, Dunhagi 3, 107 Reykjavik, Iceland
| | - Ragnar Bjornsson
- Science Institute, University of Iceland, Dunhagi 3, 107 Reykjavik, Iceland.,Max-Planck Institute for Chemical Energy Conversion, Stiftstrasse 34-36, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
15
|
Vitillo JG, Cramer CJ, Gagliardi L. Multireference Methods are Realistic and Useful Tools for Modeling Catalysis. Isr J Chem 2022. [DOI: 10.1002/ijch.202100136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jenny G. Vitillo
- Department of Science and High Technology and INSTM Università degli Studi dell'Insubria Via Valleggio 9 I-22100 Como Italy
| | - Christopher J. Cramer
- Underwriters Laboratories Inc. 333 Pfingsten Road Northbrook Illinois 60602 United States
| | - Laura Gagliardi
- Department of Chemistry Pritzker School of Molecular Engineering James Franck Institute University of Chicago Chicago Illinois 60637 United States
| |
Collapse
|
16
|
Cheng Y, Xie Z, Ma H. Post-Density Matrix Renormalization Group Methods for Describing Dynamic Electron Correlation with Large Active Spaces. J Phys Chem Lett 2022; 13:904-915. [PMID: 35049302 DOI: 10.1021/acs.jpclett.1c04078] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The ab initio density matrix renormalization group (DMRG) method has been well-established and has become one of the most accurate numerical methods for the precise electronic structure solution of large active spaces. In the past few years, to capture the missing dynamic correlation, various post-DMRG approaches have been proposed through the combination of DMRG and multireference quantum chemical methods or density functional theory. With this in mind, this work provides a brief overview of ab initio DMRG principles and the new developments within post-DMRG methods. For clarity, post-DMRG methods are classified into two main categories depending on whether high-order n-electron reduced density matrices are used, and their merits and disadvantages are properly discussed. Finally, we conclude by discussing unsolved bottlenecks and giving development perspectives of post-DMRG approaches, which are expected to yield quantitative descriptions of complex electronic structures in large strongly correlated molecules and materials.
Collapse
Affiliation(s)
- Yifan Cheng
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Jiangsu Key Laboratory of Vehicle Emissions Control, Nanjing University, Nanjing 210023, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing 210023, China
| |
Collapse
|
17
|
Phung QM, Muchammad Y, Yanai T, Ghosh A. A DMRG/CASPT2 Investigation of Metallocorroles: Quantifying Ligand Noninnocence in Archetypal 3d and 4d Element Derivatives. JACS AU 2021; 1:2303-2314. [PMID: 34984418 PMCID: PMC8717376 DOI: 10.1021/jacsau.1c00417] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/22/2021] [Indexed: 05/03/2023]
Abstract
Hybrid density functional theory (B3LYP) and density matrix renormalization group (DMRG) theory have been used to quantitatively compare the degree of ligand noninnocence (corrole radical character) in seven archetypal metallocorroles. The seven complexes, in decreasing order of corrole noninnocent character, are Mn[Cor]Cl > Fe[Cor]Cl > Fe[Cor](NO) > Mo[Cor]Cl2 > Ru[Cor](NO) ≈ Mn[Cor]Ph ≈ Fe[Cor]Ph ≈ 0, where [Cor] refers to the unsubstituted corrolato ligand. DMRG-based second-order perturbation theory calculations have also yielded detailed excited-state energetics data on the compounds, shedding light on periodic trends involving middle transition elements. Thus, whereas the ground state of Fe[Cor](NO) (S = 0) is best described as a locally S = 1/2 {FeNO}7 unit antiferromagnetically coupled to a corrole A' radical, the calculations confirm that Ru[Cor](NO) may be described as simply {RuNO}6-Cor3-, that is, having an innocent corrole macrocycle. Furthermore, whereas the ferromagnetically coupled S = 1{FeNO}7-Cor•2- state of Fe[Cor](NO) is only ∼17.5 kcal/mol higher than the S = 0 ground state, the analogous triplet state of Ru[Cor](NO) is higher by a far larger margin (37.4 kcal/mol) relative to the ground state. In the same vein, Mo[Cor]Cl2 exhibits an adiabatic doublet-quartet gap of 36.1 kcal/mol. The large energy gaps associated with metal-ligand spin coupling in Ru[Cor](NO) and Mo[Cor]Cl2 reflect the much greater covalent character of 4d-π interactions relative to analogous interactions involving 3d orbitals. As far as excited-state energetics is concerned, DMRG-CASPT2 calculations provide moderate validation for hybrid density functional theory (B3LYP) for qualitative purposes, but underscore the possibility of large errors (>10 kcal/mol) in interstate energy differences.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Yasin Muchammad
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Takeshi Yanai
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Abhik Ghosh
- Department
of Chemistry, UiT-The Arctic University
of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
18
|
Kollmar C, Sivalingam K, Guo Y, Neese F. An efficient implementation of the NEVPT2 and CASPT2 methods avoiding higher-order density matrices. J Chem Phys 2021; 155:234104. [PMID: 34937355 DOI: 10.1063/5.0072129] [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/14/2022] Open
Abstract
A factorization of the matrix elements of the Dyall Hamiltonian in N-electron valence state perturbation theory allowing their evaluation with a computational effort comparable to the one needed for the construction of the third-order reduced density matrix at the most is presented. Thus, the computational bottleneck arising from explicit evaluation of the fourth-order density matrix is avoided. It is also shown that the residual terms arising in the case of an approximate complete active space configuration interaction solution and containing even the fifth-order density matrix for two excitation classes can be evaluated with little additional effort by choosing again a favorable factorization of the corresponding matrix elements. An analogous argument is also provided for avoiding the fourth-order density matrix in complete active space second-order perturbation theory. Practical calculations indicate that such an approach leads to a considerable gain in computational efficiency without any compromise in numerical accuracy or stability.
Collapse
Affiliation(s)
- Christian Kollmar
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Kantharuban Sivalingam
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, China
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
19
|
Jara-Cortés J, Leal-Sánchez E, Francisco E, Pérez-Pimienta JA, Martín Pendás Á, Hernández-Trujillo J. Implementation of the interacting quantum atom energy decomposition using the CASPT2 method. Phys Chem Chem Phys 2021; 23:27508-27519. [PMID: 34874377 DOI: 10.1039/d1cp02837e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Edith Leal-Sánchez
- Departamento de Física y Química Teórica, Facultad de Química, UNAM, México City 04510, Mexico
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | - José A Pérez-Pimienta
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | | |
Collapse
|
20
|
Golub P, Antalik A, Veis L, Brabec J. Machine Learning-Assisted Selection of Active Spaces for Strongly Correlated Transition Metal Systems. J Chem Theory Comput 2021; 17:6053-6072. [PMID: 34570505 DOI: 10.1021/acs.jctc.1c00235] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Active space quantum chemical methods could provide very accurate description of strongly correlated electronic systems, which is of tremendous value for natural sciences. The proper choice of the active space is crucial but a nontrivial task. In this article, we present a neural network-based approach for automatic selection of active spaces, focused on transition metal systems. The training set has been formed from artificial systems composed of one transition metal and various ligands, on which we have performed the density matrix renormalization group and calculated the single-site entropy. On the selected set of systems, ranging from small benchmark molecules up to larger challenging systems involving two metallic centers, we demonstrate that our machine learning models could predict the active space orbitals with reasonable accuracy. We also tested the transferability on out-of-the-model systems, including bimetallic complexes and complexes with ligands, which were not involved in the training set. Also, we tested the correctness of the automatically selected active spaces on a Fe(II)-porphyrin model, where we studied the lowest states at the DMRG level and compared the energy difference between spin states or the energy difference between conformations of ferrocene with recent studies.
Collapse
Affiliation(s)
- Pavlo Golub
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| | - Andrej Antalik
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| | - Libor Veis
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| | - Jiri Brabec
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
21
|
Dobrautz W, Weser O, Bogdanov NA, Alavi A, Li Manni G. Spin-Pure Stochastic-CASSCF via GUGA-FCIQMC Applied to Iron-Sulfur Clusters. J Chem Theory Comput 2021; 17:5684-5703. [PMID: 34469685 PMCID: PMC8444347 DOI: 10.1021/acs.jctc.1c00589] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Indexed: 11/28/2022]
Abstract
In this work, we demonstrate how to efficiently compute the one- and two-body reduced density matrices within the spin-adapted full configuration interaction quantum Monte Carlo (FCIQMC) method, which is based on the graphical unitary group approach (GUGA). This allows us to use GUGA-FCIQMC as a spin-pure configuration interaction (CI) eigensolver within the complete active space self-consistent field (CASSCF) procedure and hence to stochastically treat active spaces far larger than conventional CI solvers while variationally relaxing orbitals for specific spin-pure states. We apply the method to investigate the spin ladder in iron-sulfur dimer and tetramer model systems. We demonstrate the importance of the orbital relaxation by comparing the Heisenberg model magnetic coupling parameters from the CASSCF procedure to those from a CI-only (CASCI) procedure based on restricted open-shell Hartree-Fock orbitals. We show that the orbital relaxation differentially stabilizes the lower-spin states, thus enlarging the coupling parameters with respect to the values predicted by ignoring orbital relaxation effects. Moreover, we find that, while CASCI results are well fit by a simple bilinear Heisenberg Hamiltonian, the CASSCF eigenvalues exhibit deviations that necessitate the inclusion of biquadratic terms in the model Hamiltonian.
Collapse
Affiliation(s)
- Werner Dobrautz
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Oskar Weser
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Nikolay A. Bogdanov
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Giovanni Li Manni
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany
| |
Collapse
|
22
|
Khedkar A, Roemelt M. Modern multireference methods and their application in transition metal chemistry. Phys Chem Chem Phys 2021; 23:17097-17112. [PMID: 34355719 DOI: 10.1039/d1cp02640b] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.
Collapse
Affiliation(s)
- Abhishek Khedkar
- Lehrstuhl für theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany.
| | | |
Collapse
|
23
|
Guo Y, Sivalingam K, Kollmar C, Neese F. Approximations of density matrices in N-electron valence state second-order perturbation theory (NEVPT2). II. The full rank NEVPT2 (FR-NEVPT2) formulation. J Chem Phys 2021; 154:214113. [PMID: 34240984 DOI: 10.1063/5.0051218] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In Paper I, the performances of pre-screening (PS), extended PS (EPS), and cumulant (CU) approximations to the fourth-order density matrix were examined in the context of second-order N-electron valence state perturbation theory (NEVPT2). It has been found that the CU, PS, and even EPS approximations with loose thresholds may introduce intruder states. In the present work, the origin of these "false intruder" states introduced by approximated density matrices is discussed. Canonical NEVPT2 implementations employ a rank reduction trick. By analyzing its residual error, we find that the omission of the rank reduction leads to a more stable multireference perturbation theory for incomplete active space reference wave functions. Such a full rank (FR)-NEVPT2 formulation is equivalent to the conventional NEVPT2 method for the complete active space self-consistent field/complete active space configuration interaction reference wave function. A major drawback of the FR-NEVPT2 formulation is the necessity of the fifth-order density matrix. To avoid the construction of the high-order density matrices, the combination of the FR-NEVPT2 with the CU approximation is studied. However, we find that the CU approximation remains problematic as it still introduces intruder states. The question of how to robustly and efficiently perform internally contracted multireference perturbation theories with approximate densities remains a challenging field of investigation.
Collapse
Affiliation(s)
- Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, China
| | - Kantharuban Sivalingam
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Christian Kollmar
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
24
|
Guo Y, Sivalingam K, Neese F. Approximations of density matrices in N-electron valence state second-order perturbation theory (NEVPT2). I. Revisiting the NEVPT2 construction. J Chem Phys 2021; 154:214111. [PMID: 34240991 DOI: 10.1063/5.0051211] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Over the last decade, the second-order N-electron valence state perturbation theory (NEVPT2) has developed into a widely used multireference perturbation method. To apply NEVPT2 to systems with large active spaces, the computational bottleneck is the construction of the fourth-order reduced density matrix. Both its generation and storage become quickly problematic beyond the usual maximum active space of about 15 active orbitals. To reduce the computational cost of handling fourth-order density matrices, the cumulant approximation (CU) has been proposed in several studies. A more conventional strategy to address the higher-order density matrices is the pre-screening approximation (PS), which is the default one in the ORCA program package since 2010. In the present work, the performance of the CU, PS, and extended PS (EPS) approximations for the fourth-order density matrices is compared. Following a pedagogical introduction to NEVPT2, contraction schemes, as well as the approximations to density matrices, and the intruder state problem are discussed. The CU approximation, while potentially leading to large computational savings, virtually always leads to intruder states. With the PS approximation, the computational savings are more modest. However, in conjunction with conservative cutoffs, it produces stable results. The EPS approximation to the fourth-order density matrices can reproduce very accurate NEVPT2 results without any intruder states. However, its computational cost is not much lower than that of the canonical algorithm. Moreover, we found that a good indicator of intrude states problems in any approximation to high order density matrices is the eigenspectra of the Koopmans matrices.
Collapse
Affiliation(s)
- Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, China
| | - Kantharuban Sivalingam
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
28
|
Anderson RJ, Shiozaki T, Booth GH. Efficient and stochastic multireference perturbation theory for large active spaces within a full configuration interaction quantum Monte Carlo framework. J Chem Phys 2020; 152:054101. [DOI: 10.1063/1.5140086] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Robert J. Anderson
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Toru Shiozaki
- Quantum Simulation Technologies, Inc., Cambridge, Massachusetts 02139, USA
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| |
Collapse
|
29
|
Backhouse OJ, Nusspickel M, Booth GH. Wave Function Perspective and Efficient Truncation of Renormalized Second-Order Perturbation Theory. J Chem Theory Comput 2020; 16:1090-1104. [PMID: 31951406 DOI: 10.1021/acs.jctc.9b01182] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present an approach to a renormalized second-order Green's function perturbation theory (GF2), which avoids all dependency on continuous variables, grids, or explicit Green's functions and is instead formulated entirely in terms of static quantities and wave functions. Correlation effects from MP2 diagrams are iteratively incorporated to modify the underlying spectrum of excitations by coupling the physical system to fictitious auxiliary degrees of freedom, allowing for single-particle orbitals to delocalize into this additional space. The overall approach is shown to be rigorously O[N5], after an appropriate compression of this auxiliary space. This is achieved via a novel scheme, which ensures that a desired number of moments of the underlying occupied and virtual spectra are conserved in the compression, allowing a rapid and systematically improvable convergence to the limit of the effective dynamical resolution. The approach is found to then allow for the qualitative description of stronger correlation effects, avoiding the divergences of MP2, as well as its orbital-optimized version. On application to the G1 test set, we find that modification up to only the third spectral moment of the underlying spectrum from which the double excitations are built are required for accurate energetics, even in strongly correlated regimes. This is beyond simple self-consistent changes to the density matrix of the system but far from requiring a description of the full dynamics of the frequency-dependent self-energy.
Collapse
Affiliation(s)
- Oliver J Backhouse
- Department of Physics , King's College London , Strand , London WC2R 2LS , U.K
| | - Max Nusspickel
- Department of Physics , King's College London , Strand , London WC2R 2LS , U.K
| | - George H Booth
- Department of Physics , King's College London , Strand , London WC2R 2LS , U.K
| |
Collapse
|
30
|
Mahajan A, Blunt NS, Sabzevari I, Sharma S. Multireference configuration interaction and perturbation theory without reduced density matrices. J Chem Phys 2019; 151:211102. [DOI: 10.1063/1.5128115] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ankit Mahajan
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Nick S. Blunt
- Department of Chemistry, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Iliya Sabzevari
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| |
Collapse
|
31
|
Freitag L, Ma Y, Baiardi A, Knecht S, Reiher M. Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method. J Chem Theory Comput 2019; 15:6724-6737. [PMID: 31670947 DOI: 10.1021/acs.jctc.9b00969] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.
Collapse
Affiliation(s)
- Leon Freitag
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Yingjin Ma
- Computer Network Information Center , Chinese Academy of Sciences , Beijing 100190 , China.,Center of Scientific Computing Applications & Research, Chinese Academy of Sciences , Beijing 100190 , China
| | - Alberto Baiardi
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Stefan Knecht
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| |
Collapse
|
32
|
Nakagaki M, Nakatani N, Sakaki S. How to understand very weak Cr-Cr double bonds and negative spin populations in trinuclear Cr complexes: theoretical insight. Phys Chem Chem Phys 2019; 21:22976-22989. [PMID: 31599298 DOI: 10.1039/c9cp03645h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Trinuclear Cr(ii) complex [Cr3(dpa)4Cl2] 1 (Hdpa = dipyridylamine) has two Cr-Cr double bonds linked with each other. DMRG-CASPT2 calculations reproduced its symmetrical structure. The Cr-Cr effective bond order (EBO) was evaluated to be only 0.59 based on the density matrix based on localized orbitals from DMRG-CASSCF orbitals. The CASCI calculations showed a significantly large α-spin population on the terminal Cr atoms as expected but a significantly large β-spin population on the central Cr atom against expectations. The very small EBO and the presence of a large β-spin population are not consistent with the simple understanding that 1 has two Cr-Cr double bonds and a quintet ground state, which requests correct understanding of 1 from the viewpoint of chemical bond theory. Comparison of 1 with the allene molecule and allyl radical disclosed that the linked Cr-Cr bonds of 1 resembled the C-C bond of the allyl radical but completely differed from the linked C-C double bonds of allene despite the similar molecular structure. Its N3 analogue [Cr3(dpa)4(N3)2] 2 has non-symmetrical structure with shorter Cr1-Cr2 and longer Cr2-Cr3 bonds unlike 1, indicating that 2 is a valence tautomer of 1. DMRG-CASPT2 could reproduce its non-symmetrical structure but DFT/B3PW91 could not. In 2, the EBO is 0.95 for the shorter Cr1-Cr2 bond and 0.47 for the longer Cr2-Cr3 one. The terminal Cr3 has a very large α spin population, and the other terminal Cr1 has a somewhat large α spin population, but the central Cr2 has a considerably large β spin population. These results indicate that the Cr1-Cr2 bond conjugates with the Cr2-Cr3 bond, which is inconsistent with the simple understanding that 2 has a quadruple bond between Cr1 and Cr2 and no bond between Cr2 and Cr3. The symmetrical structure has a stronger Cr-X coordinate bond (X = Cl or N3) but less stable Cr3 core than does the non-symmetrical one. The relative stabilities of the symmetrical and the non-symmetrical structures are determined by the balance between stabilization energies from the Cr3 core and the Cr-X coordinate bond. All these findings show that electronic structures and Cr-Cr bonds of 1 and 2 are interesting from the viewpoint of molecular science.
Collapse
Affiliation(s)
- Masayuki Nakagaki
- Fukui Institute for Fundamental Chemistry, Kyoto University, Takano-Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan.
| | | | | |
Collapse
|
33
|
Fdez. Galván I, Vacher M, Alavi A, Angeli C, Aquilante F, Autschbach J, Bao JJ, Bokarev SI, Bogdanov NA, Carlson RK, Chibotaru LF, Creutzberg J, Dattani N, Delcey MG, Dong SS, Dreuw A, Freitag L, Frutos LM, Gagliardi L, Gendron F, Giussani A, González L, Grell G, Guo M, Hoyer CE, Johansson M, Keller S, Knecht S, Kovačević G, Källman E, Li Manni G, Lundberg M, Ma Y, Mai S, Malhado JP, Malmqvist PÅ, Marquetand P, Mewes SA, Norell J, Olivucci M, Oppel M, Phung QM, Pierloot K, Plasser F, Reiher M, Sand AM, Schapiro I, Sharma P, Stein CJ, Sørensen LK, Truhlar DG, Ugandi M, Ungur L, Valentini A, Vancoillie S, Veryazov V, Weser O, Wesołowski TA, Widmark PO, Wouters S, Zech A, Zobel JP, Lindh R. OpenMolcas: From Source Code to Insight. J Chem Theory Comput 2019; 15:5925-5964. [DOI: 10.1021/acs.jctc.9b00532] [Citation(s) in RCA: 399] [Impact Index Per Article: 79.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Ignacio Fdez. Galván
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
- Department of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Morgane Vacher
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Ali Alavi
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Celestino Angeli
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Ferrara, Via Luigi Borsari 46, 44121 Ferrara, Italy
| | - Francesco Aquilante
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Jie J. Bao
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Sergey I. Bokarev
- Institut für Physik, Universität Rostock, Albert-Einstein-Straße 23-24, 18059 Rostock, Germany
| | - Nikolay A. Bogdanov
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Rebecca K. Carlson
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Liviu F. Chibotaru
- Theory of Nanomaterials Group, University of Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Joel Creutzberg
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Nike Dattani
- Harvard Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138, United States
| | - Mickaël G. Delcey
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Sijia S. Dong
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205 A, 69120 Heidelberg, Germany
| | - Leon Freitag
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Luis Manuel Frutos
- Departamento de Química Analítica, Química Física e Ingeniería Química, and Instituto de Investigación Química “Andrés M. del Río”, Universidad de Alcalá, E-28871 Alcalá de Henares, Madrid, Spain
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Frédéric Gendron
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Angelo Giussani
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
- Instituto de Ciencia Molecular, Universitat de València, Apartado 22085, ES-46071 Valencia, Spain
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Gilbert Grell
- Institut für Physik, Universität Rostock, Albert-Einstein-Straße 23-24, 18059 Rostock, Germany
| | - Meiyuan Guo
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Chad E. Hoyer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Marcus Johansson
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Sebastian Keller
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Stefan Knecht
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Goran Kovačević
- Division of Materials Physics, Ruđer Bošković Institute, P.O.B. 180, Bijenička 54, HR-10002 Zagreb, Croatia
| | - Erik Källman
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Giovanni Li Manni
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Yingjin Ma
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sebastian Mai
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - João Pedro Malhado
- Department of Chemistry, Imperial College London, London SW7 2AZ, United Kingdom
| | - Per Åke Malmqvist
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205 A, 69120 Heidelberg, Germany
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study (NZIAS), Massey University Albany, Private Bag
102904, Auckland 0632, New Zealand
| | - Jesper Norell
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, via A. Moro 2, 53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
- USIAS and Institut de Physique et Chimie des Matériaux de Strasbourg, Université de Strasbourg-CNRS, 67034 Strasbourg, France
| | - Markus Oppel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Quan Manh Phung
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven 3001, Belgium
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Andrew M. Sand
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Igor Schapiro
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Prachi Sharma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Stein
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Lasse Kragh Sørensen
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Mihkel Ugandi
- Department of Chemistry − Ångström Laboratory, Uppsala University, P.O. Box 538, SE-751 21 Uppsala, Sweden
| | - Liviu Ungur
- Department of Chemistry, National University of Singapore, 117543 Singapore
| | - Alessio Valentini
- Theoretical Physical Chemistry, Research Unit MolSys, Allée du 6 Août, 11, 4000 Liège, Belgium
| | - Steven Vancoillie
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Valera Veryazov
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Oskar Weser
- Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Tomasz A. Wesołowski
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Per-Olof Widmark
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Sebastian Wouters
- Brantsandpatents, Pauline van Pottelsberghelaan 24, 9051 Sint-Denijs-Westrem, Belgium
| | - Alexander Zech
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - J. Patrick Zobel
- Division of Theoretical Chemistry, Kemicentrum, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Roland Lindh
- Department of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
- Uppsala Center for Computational Chemistry (UC3), Uppsala University, P.O. Box 596, SE-751 24 Uppsala, Sweden
| |
Collapse
|
34
|
Phung QM, Pierloot K. Electronic Structure of N-Bridged High-Valent Diiron-Oxo. Chemistry 2019; 25:12491-12496. [PMID: 31313861 DOI: 10.1002/chem.201902766] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2019] [Revised: 07/14/2019] [Indexed: 11/09/2022]
Abstract
Density functional theory (DFT) and an advanced ab initio technique based on density matrix renormalization group (DMRG-CASPT2) were employed to investigate a reactive N-bridged high-valent diiron-oxo species involved in H-abstraction reactions. We studied in detail two important doublet states, the ground state with two iron(IV) centers and a mixed valence FeV -FeIV excited state. We found that the latter state is low-lying. Furthermore, its electronic structure and spin density imply that it has significantly higher H-abstraction reactivity than the ground state. This low-lying excited state might be the reason behind the high oxidation reactivity of this diiron-oxo species towards methane.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven, 3001, Belgium.,Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Chikusa, Nagoya, 464-8602, Japan
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven, 3001, Belgium
| |
Collapse
|
35
|
Phung QM, Martín-Fernández C, Harvey JN, Feldt M. Ab Initio Calculations for Spin-Gaps of Non-Heme Iron Complexes. J Chem Theory Comput 2019; 15:4297-4304. [DOI: 10.1021/acs.jctc.9b00370] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry, KU Leuven, Celestijnenlaan 200f, Box 2404, 3001 Leuven, Belgium
| | | | - Jeremy N. Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200f, Box 2404, 3001 Leuven, Belgium
| | - Milica Feldt
- Department of Chemistry, KU Leuven, Celestijnenlaan 200f, Box 2404, 3001 Leuven, Belgium
| |
Collapse
|
36
|
Reaction mechanism between Cu(II)-enolate complex and O2 as a test case for methodology used in DFT computational studies. J Mol Model 2019; 25:122. [DOI: 10.1007/s00894-019-3998-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Accepted: 03/15/2019] [Indexed: 11/26/2022]
|
37
|
Phung QM, Pierloot K. Low-Lying Electromeric States in Chloro-Ligated Iron(IV)-Oxo Porphyrin as a Model for Compound I, Studied with Second-Order Perturbation Theory Based on Density Matrix Renormalization Group. J Chem Theory Comput 2019; 15:3033-3043. [PMID: 30995039 DOI: 10.1021/acs.jctc.9b00166] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Employing second-order perturbation theory based on the density matrix renormalization group (DMRG-CASPT2), this work aims at providing a quantitative description of the spin state energetics of a chloro-ligated iron(IV)-oxo porphyrin as a model for the cytochromes P450 active species, also known as compound I (Cpd I). We explored DMRG-CASPT2 to its full extent with an extensive active space (up to 31 active orbitals) as well as a large number of renormalized states m (up to 10000). Different flavors of DMRG-CASPT2, using either the costly exact 4-particle reduced density matrix (4-RDM) or the cheaper cumulant approximated 4-RDM (cu(4)), were analyzed. All flavors essentially converge to similar relative energies between different spin states. Including a correction for the protein environment, we found a quartet FeIVO ground state and, more importantly, a thermally accessible doublet FeVO excited state that might directly contribute to the reactivity of this iron-oxo species. Our results also showed that cheaper approaches, such as CASPT2 based on a smaller active space or the cumulant approximation DMRG-cu(4)-CASPT2, are capable of accurately describing the spin state energetics of this species.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , Leuven 3001 , Belgium.,Institute of Transformative Bio-Molecules (WPI-ITbM) , Nagoya University , Chikusa, Nagoya 464-8602 , Japan
| | - Kristine Pierloot
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , Leuven 3001 , Belgium
| |
Collapse
|
38
|
Vogiatzis KD, Polynski MV, Kirkland JK, Townsend J, Hashemi A, Liu C, Pidko EA. Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities. Chem Rev 2019; 119:2453-2523. [PMID: 30376310 PMCID: PMC6396130 DOI: 10.1021/acs.chemrev.8b00361] [Citation(s) in RCA: 214] [Impact Index Per Article: 42.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Indexed: 12/28/2022]
Abstract
Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.
Collapse
Affiliation(s)
| | | | - Justin K. Kirkland
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Jacob Townsend
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ali Hashemi
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Chong Liu
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Evgeny A. Pidko
- TheoMAT
group, ITMO University, Lomonosova 9, St. Petersburg 191002, Russia
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| |
Collapse
|
39
|
Földvári D, Tóth Z, Surján PR, Szabados Á. Geminal perturbation theory based on the unrestricted Hartree-Fock wavefunction. J Chem Phys 2019; 150:034103. [PMID: 30660159 DOI: 10.1063/1.5060731] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A perturbative correction exploiting natural orbitals and the pair function structure of the unrestricted Hartree-Fock (UHF) wavefunction is devised. The method offers a simple framework for describing multireference systems where static correlation is captured by UHF. The UHF wavefunction is built of two-electron fragments (geminals), involving both singlet and triplet (ms = 0) parts. At order zero of the perturbative treatment, configuration interaction coefficients of UHF geminals are relaxed. The zero order Hamiltonian is of the Dyall-type, including explicit two-electron interaction within geminals and leading to a formal 6th power scaling. Adopting an effective one-electron zero order Hamiltonian term for the subset of virtual orbitals reduces scaling of the correction step to 4th power. Formal properties of the proposed schemes are discussed. Energetic data and natural occupation numbers of illustrative test systems are used to assess the new approach. The cases where the wavefunction becomes essentially spin pure at the level of reference show good performance. Spin contamination remaining at order zero is found to undermine the perturbative correction.
Collapse
Affiliation(s)
- D Földvári
- Laboratory of Theoretical Chemistry, Faculty of Science, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Zs Tóth
- Laboratory of Theoretical Chemistry, Faculty of Science, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - P R Surján
- Laboratory of Theoretical Chemistry, Faculty of Science, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| | - Á Szabados
- Laboratory of Theoretical Chemistry, Faculty of Science, Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
40
|
Nakatani N, Hada M. Can large active-space CASSCF calculation make sense to the reaction analysis of iron complex? A benchmark study of methane oxidation reaction by FeO . J Comput Chem 2019; 40:414-420. [PMID: 30351477 DOI: 10.1002/jcc.25640] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Revised: 09/14/2018] [Accepted: 09/14/2018] [Indexed: 01/04/2023]
Abstract
A methane oxidation reaction by FeO+ cation was theoretically investigated based on the density functional theory (DFT) and the complete active-space self-consistent field (CASSCF) method as well as the coupled-cluster singles, doubles, and perturbative triples (CCSD(T)) to explore the active-space dependency to computational analyses in such strongly correlated reaction systems. A small active-space CASSCF(5e in 5o) calculation, which only includes five 3d orbitals of the Fe atom in the active-space, showed remarkable difference both in energy and geometry compared to those computed by the DFT and CCSD(T) methods. Interestingly, a large active-space CASSCF(17e in 17o) calculation, which includes almost all the valence orbitals gives a qualitative agreement with either the DFT or the CCSD(T) results in the first half part of the reaction, although it varies from them in the latter half part. Therefore, it is indicated that the active-space dependency is serious in some part of the reaction and the small active-space CASSCF might lead a wrong discussion. We further investigated the optimized geometry of the intermediate complex with the small and the large active-space CASSCF methods as well as the CCSD(T) method, and found that the CASSCF(5e in 5o)-optimized geometry is considerably different from the others. In consequence, a small active-space CASSCF/CASPT2 calculation does not really work for such a strongly correlated reaction system even qualitatively, and a sophisticated assessment using the large active-space CASSCF/CASPT2 method will be indispensable. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Naoki Nakatani
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa Hachioji-shi, Tokyo 192-0397, Japan
| | - Masahiko Hada
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, 1-1 Minami-Osawa Hachioji-shi, Tokyo 192-0397, Japan
| |
Collapse
|
41
|
Dong G, Phung QM, Pierloot K, Ryde U. Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods. Inorg Chem 2018; 57:15289-15298. [PMID: 30500163 DOI: 10.1021/acs.inorgchem.8b02590] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
[NiFe] hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons. This seemingly simple reaction has attracted much attention because of the prospective use of H2 as a clean fuel. In this paper, we have studied the full reaction mechanism of this enzyme with various computational methods. Geometries were obtained with combined quantum mechanical and molecular mechanics (QM/MM) calculations. To get more accurate energies and obtain a detailed account of the surroundings, we performed big-QM calculations with 819 atoms in the QM region. Moreover, QM/MM thermodynamic cycle perturbation calculations were performed to obtain free energies. Finally, density matrix renormalisation group complete active space self-consistent field calculations were carried out to study the electronic structures of the various states in the reaction mechanism. Our calculations indicate that the Ni-L state is not involved in the reaction mechanism. Instead, the Ni-C state is reduced by one electron and then the bridging hydride ion is transferred to the sulfur atom of Cys546 as a proton and the two electrons transfer to the Ni ion. This step turned out to be rate-determining with an energy barrier of 58 kJ/mol, which is consistent with the experimental rate of 750 ± 90 s-1 (corresponding to ∼52 kJ/mol). The cleavage of the H-H bond is facile with an energy barrier of 33 kJ/mol, according to our calculations. We also find that the reaction energies are sensitive to the size of the QM system, the basis set, and the density functional theory method, in agreement with previous studies.
Collapse
Affiliation(s)
- Geng Dong
- Department of Theoretical Chemistry, Chemical Centre , Lund University , P.O. Box 124, SE-221 00 Lund , Sweden
- Department of Biochemistry and Molecular Biology , Shantou University Medical College , Shantou 514041 , Guangdong , PR China
| | - Quan Manh Phung
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Kristine Pierloot
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre , Lund University , P.O. Box 124, SE-221 00 Lund , Sweden
| |
Collapse
|
42
|
Sharma P, Bernales V, Knecht S, Truhlar DG, Gagliardi L. Density matrix renormalization group pair-density functional theory (DMRG-PDFT): singlet-triplet gaps in polyacenes and polyacetylenes. Chem Sci 2018; 10:1716-1723. [PMID: 30842836 PMCID: PMC6368241 DOI: 10.1039/c8sc03569e] [Citation(s) in RCA: 61] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2018] [Accepted: 11/23/2018] [Indexed: 11/21/2022] Open
Abstract
The density matrix renormalization group (DMRG) is a powerful method to treat static correlation.
The density matrix renormalization group (DMRG) is a powerful method to treat static correlation. Here we present an inexpensive way to calculate correlation energy starting from a DMRG wave function using pair-density functional theory (PDFT). We applied this new approach, called DMRG-PDFT, to study singlet–triplet gaps in polyacenes and polyacetylenes that require active spaces larger than the feasibility limit of the conventional complete active-space self-consistent field (CASSCF) method. The results match reasonably well with the most reliable literature values and have only a moderate dependence on the compression of the initial DMRG wave function. Furthermore, DMRG-PDFT is significantly less expensive than other commonly applied ways of adding additional correlation to DMRG, such as DMRG followed by multireference perturbation theory or multireference configuration interaction.
Collapse
Affiliation(s)
- Prachi Sharma
- Department of Chemistry , Chemical Theory Center , Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , USA . ;
| | - Varinia Bernales
- Department of Chemistry , Chemical Theory Center , Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , USA . ;
| | - Stefan Knecht
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , CH-8093 Zürich , Switzerland .
| | - Donald G Truhlar
- Department of Chemistry , Chemical Theory Center , Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , USA . ;
| | - Laura Gagliardi
- Department of Chemistry , Chemical Theory Center , Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , USA . ;
| |
Collapse
|
43
|
Schlimgen AW, Mazziotti DA. Analytical gradients of variational reduced-density-matrix and wavefunction-based methods from an overlap-reweighted semidefinite program. J Chem Phys 2018; 149:164111. [DOI: 10.1063/1.5043104] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anthony W. Schlimgen
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - David A. Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
44
|
Lischka H, Nachtigallová D, Aquino AJA, Szalay PG, Plasser F, Machado FBC, Barbatti M. Multireference Approaches for Excited States of Molecules. Chem Rev 2018; 118:7293-7361. [DOI: 10.1021/acs.chemrev.8b00244] [Citation(s) in RCA: 197] [Impact Index Per Article: 32.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hans Lischka
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 16610 Prague 6, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Palacký University, 78371 Olomouc, Czech Republic
| | - Adélia J. A. Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, P.R. China
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
- Institute for Soil Research, University of Natural Resources and Life Sciences Vienna, Peter-Jordan-Strasse 82, A-1190 Vienna, Austria
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Laboratory of Theoretical Chemistry, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary
| | - Felix Plasser
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
- Department of Chemistry, Loughborough University, Leicestershire LE11 3TU, United Kingdom
| | - Francisco B. C. Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | | |
Collapse
|
45
|
Guo S, Li Z, Chan GKL. A Perturbative Density Matrix Renormalization Group Algorithm for Large Active Spaces. J Chem Theory Comput 2018; 14:4063-4071. [DOI: 10.1021/acs.jctc.8b00273] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sheng Guo
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Zhendong Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
46
|
Guo S, Li Z, Chan GKL. Communication: An efficient stochastic algorithm for the perturbative density matrix renormalization group in large active spaces. J Chem Phys 2018; 148:221104. [DOI: 10.1063/1.5031140] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sheng Guo
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Zhendong Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
47
|
Schilling M, Luber S. Computational Modeling of Cobalt-Based Water Oxidation: Current Status and Future Challenges. Front Chem 2018; 6:100. [PMID: 29721491 PMCID: PMC5915471 DOI: 10.3389/fchem.2018.00100] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Accepted: 03/20/2018] [Indexed: 12/19/2022] Open
Abstract
A lot of effort is nowadays put into the development of novel water oxidation catalysts. In this context, mechanistic studies are crucial in order to elucidate the reaction mechanisms governing this complex process, new design paradigms and strategies how to improve the stability and efficiency of those catalysts. This review is focused on recent theoretical mechanistic studies in the field of homogeneous cobalt-based water oxidation catalysts. In the first part, computational methodologies and protocols are summarized and evaluated on the basis of their applicability toward real catalytic or smaller model systems, whereby special emphasis is laid on the choice of an appropriate model system. In the second part, an overview of mechanistic studies is presented, from which conceptual guidelines are drawn on how to approach novel studies of catalysts and how to further develop the field of computational modeling of water oxidation reactions.
Collapse
Affiliation(s)
- Mauro Schilling
- Department of Chemistry, University of Zürich, Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Zurich, Switzerland
| |
Collapse
|
48
|
Phung QM, Feldt M, Harvey JN, Pierloot K. Toward Highly Accurate Spin State Energetics in First-Row Transition Metal Complexes: A Combined CASPT2/CC Approach. J Chem Theory Comput 2018; 14:2446-2455. [PMID: 29614218 DOI: 10.1021/acs.jctc.8b00057] [Citation(s) in RCA: 82] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
In previous work on the performance of multiconfigurational second-order perturbation theory (CASPT2) in describing spin state energetics in first-row transition metal systems [ Pierloot et al. J. Chem. Theory Comput. 2017 , 13 , 537 - 553 ], we showed that standard CASPT2 works well for valence correlation but does not describe the metal semicore (3s3p) correlation effects accurately. This failure is partially responsible for the well-known bias toward high-spin states of CASPT2. In this paper, we expand our previous work and show that this bias could be partly removed with a combined CASPT2/CC approach: using high-quality CASPT2 with extensive correlation-consistent basis sets for valence correlation and low-cost CCSD(T) calculations with minimal basis sets for the metal semicore (3s3p) correlation effects. We demonstrate that this approach is efficient by studying the spin state energetics of a series of iron complexes modeling important intermediates in oxidative catalytic processes in chemistry and biochemistry. On the basis of a comparison with bare CCSD(T) results from this and previous work, the average error of the CASPT2/CC approach is estimated at around 2 kcal mol-1 in favor of high spin states.
Collapse
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Milica Feldt
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Jeremy N Harvey
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| | - Kristine Pierloot
- Department of Chemistry , KU Leuven , Celestijnenlaan 200F , B-3001 Leuven , Belgium
| |
Collapse
|
49
|
Roemelt M, Krewald V, Pantazis DA. Exchange Coupling Interactions from the Density Matrix Renormalization Group and N-Electron Valence Perturbation Theory: Application to a Biomimetic Mixed-Valence Manganese Complex. J Chem Theory Comput 2017; 14:166-179. [PMID: 29211960 DOI: 10.1021/acs.jctc.7b01035] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The accurate description of magnetic level energetics in oligonuclear exchange-coupled transition-metal complexes remains a formidable challenge for quantum chemistry. The density matrix renormalization group (DMRG) brings such systems for the first time easily within reach of multireference wave function methods by enabling the use of unprecedentedly large active spaces. But does this guarantee systematic improvement in predictive ability and, if so, under which conditions? We identify operational parameters in the use of DMRG using as a test system an experimentally characterized mixed-valence bis-μ-oxo/μ-acetato Mn(III,IV) dimer, a model for the oxygen-evolving complex of photosystem II. A complete active space of all metal 3d and bridge 2p orbitals proved to be the smallest meaningful starting point; this is readily accessible with DMRG and greatly improves on the unrealistic metal-only configuration interaction or complete active space self-consistent field (CASSCF) values. Orbital optimization is critical for stabilizing the antiferromagnetic state, while a state-averaged approach over all spin states involved is required to avoid artificial deviations from isotropic behavior that are associated with state-specific calculations. Selective inclusion of localized orbital subspaces enables probing the relative contributions of different ligands and distinct superexchange pathways. Overall, however, full-valence DMRG-CASSCF calculations fall short of providing a quantitative description of the exchange coupling owing to insufficient recovery of dynamic correlation. Quantitatively accurate results can be achieved through a DMRG implementation of second order N-electron valence perturbation theory (NEVPT2) in conjunction with a full-valence metal and ligand active space. Perspectives for future applications of DMRG-CASSCF/NEVPT2 to exchange coupling in oligonuclear clusters are discussed.
Collapse
Affiliation(s)
- Michael Roemelt
- Lehrstuhl für Theoretische Chemie, Ruhr-University Bochum , 44780 Bochum, Germany.,Max Planck Institute for Coal Research, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Vera Krewald
- Department of Chemistry, University of Bath , Bath BA2 7AY, United Kingdom
| | - Dimitrios A Pantazis
- Max Planck Institute for Chemical Energy Conversion, Stiftstraße 34-36, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
50
|
Phung QM, Domingo A, Pierloot K. Dinuclear Iron(II) Spin-Crossover Compounds: A Theoretical Study. Chemistry 2017; 24:5183-5190. [DOI: 10.1002/chem.201704441] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2017] [Indexed: 11/06/2022]
Affiliation(s)
- Quan Manh Phung
- Department of Chemistry; KU Leuven; Celestijnenlaan 200F Leuven 3001 Belgium
| | - Alex Domingo
- Department of Chemistry; KU Leuven; Celestijnenlaan 200F Leuven 3001 Belgium
| | - Kristine Pierloot
- Department of Chemistry; KU Leuven; Celestijnenlaan 200F Leuven 3001 Belgium
| |
Collapse
|