1
|
Weymuth T, Unsleber JP, Türtscher PL, Steiner M, Sobez JG, Müller CH, Mörchen M, Klasovita V, Grimmel SA, Eckhoff M, Csizi KS, Bosia F, Bensberg M, Reiher M. SCINE-Software for chemical interaction networks. J Chem Phys 2024; 160:222501. [PMID: 38857173 DOI: 10.1063/5.0206974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/09/2024] [Indexed: 06/12/2024] Open
Abstract
The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.
Collapse
Affiliation(s)
- Thomas Weymuth
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P Unsleber
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L Türtscher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan-Grimo Sobez
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Charlotte H Müller
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Veronika Klasovita
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A Grimmel
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Eckhoff
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Katja-Sophia Csizi
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Francesco Bosia
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Moritz Bensberg
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
2
|
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
|
3
|
Xiang C, Jia W, Fang WH, Li Z. Distributed Multi-GPU Ab Initio Density Matrix Renormalization Group Algorithm with Applications to the P-Cluster of Nitrogenase. J Chem Theory Comput 2024; 20:775-786. [PMID: 38198503 DOI: 10.1021/acs.jctc.3c01228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
The presence of many degenerate d/f orbitals makes polynuclear transition-metal compounds, such as iron-sulfur clusters in nitrogenase, challenging for state-of-the-art quantum chemistry methods. To address this challenge, we present the first distributed multi-graphics processing unit (GPU) ab initio density matrix renormalization group (DMRG) algorithm suitable for modern high-performance computing (HPC) infrastructures. The central idea is to parallelize the most computationally intensive part─the multiplication of O(K2) operators with a trial wave function, where K is the number of spatial orbitals, by combining operator parallelism for distributing the workload with a batched algorithm for performing contractions on GPU. With this new implementation, we are able to reach an unprecedentedly large bond dimension D = 14,000 on 48 GPUs (NVIDIA A100 80 GB SXM) for an active space model (114 electrons in 73 active orbitals) of the P-cluster, which is nearly 3 times larger than the bond dimensions reported in previous DMRG calculations for the same system using only central processing units (CPUs).
Collapse
Affiliation(s)
- Chunyang Xiang
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Weile Jia
- State Key Lab of Processors, Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100190, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
4
|
Friesecke G, Barcza G, Legeza Ö. Predicting the FCI Energy of Large Systems to Chemical Accuracy from Restricted Active Space Density Matrix Renormalization Group Calculations. J Chem Theory Comput 2024; 20:87-102. [PMID: 38109339 DOI: 10.1021/acs.jctc.3c01001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [J. Phys. Chem. A 126, 9709] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the chromium dimer, dicarbon up to a large cc-pVQZ basis and even a large chemical complex such as the FeMoco with significantly lower computational demands than those of previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry and more generally for the computation of strong correlations in nuclear and condensed matter physics.
Collapse
Affiliation(s)
- Gero Friesecke
- Department of Mathematics, Technical University of Munich, München 85748, Germany
| | - Gergely Barcza
- Strongly Correlated Systems Lendület Research Group, Wigner Research Centre for Physics, Budapest H-1525, Hungary
| | - Örs Legeza
- Strongly Correlated Systems Lendület Research Group, Wigner Research Centre for Physics, Budapest H-1525, Hungary
- Institute for Advanced Study, Technical University of Munich, Germany, Lichtenbergstrasse 2a, Garching 85748, Germany
| |
Collapse
|
5
|
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
|
6
|
King DS, Truhlar DG, Gagliardi L. Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2023; 19:8118-8128. [PMID: 37905518 DOI: 10.1021/acs.jctc.3c00792] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.
Collapse
Affiliation(s)
- Daniel S King
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Group, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| |
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
|
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
|
9
|
Liao K, Zhai H, Christlmaier EM, Schraivogel T, Ríos PL, Kats D, Alavi A. Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules. J Chem Theory Comput 2023; 19:1734-1743. [PMID: 36912635 DOI: 10.1021/acs.jctc.2c01207] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2023]
Abstract
We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.
Collapse
Affiliation(s)
- Ke Liao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Huanchen Zhai
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Pablo López Ríos
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany.,Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
10
|
Xie Z, Song Y, Peng F, Li J, Cheng Y, Zhang L, Ma Y, Tian Y, Luo Z, Ma H. Kylin 1.0: An ab-initio density matrix renormalization group quantum chemistry program. J Comput Chem 2023; 44:1316-1328. [PMID: 36809661 DOI: 10.1002/jcc.27085] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 01/15/2023] [Accepted: 01/19/2023] [Indexed: 02/23/2023]
Abstract
The accurate evaluation of electron correlations is highly necessary for the proper descriptions of the electronic structures in strongly correlated molecules, ranging from bond-dissociating molecules, polyradicals, to large conjugated molecules and transition metal complexes. For this purpose, in this paper, a new ab-initio quantum chemistry program Kylin 1.0 for electron correlation calculations at various quantum many-body levels, including configuration interaction (CI), perturbation theory (PT), and density matrix renormalization group (DMRG), is presented. Furthermore, fundamental quantum chemistry methods such as Hartree-Fock self-consistent field (HF-SCF) and the complete active space SCF (CASSCF) are also implemented. The Kylin 1.0 program possesses the following features: (1) a matrix product operator (MPO) formulation-based efficient DMRG implementation for describing static electron correlation within a large active space composed of more than 100 orbitals, supporting both U 1 n × U 1 S z $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{U}{(1)}_{{\mathrm{S}}_{\mathrm{z}}} $$ and U 1 n × SU 2 S $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{SU}{(2)}_{\mathrm{S}} $$ symmetries; (2) an efficient second-order DMRG-self-consistent field (SCF) implementation; (3) an externally contracted multi-reference CI (MRCI) and Epstein-Nesbet PT with DMRG reference wave functions for including the remaining dynamic electron correlation outside the large active spaces. In this paper, we introduce the capabilities and numerical benchmark examples of the Kylin 1.0 program.
Collapse
Affiliation(s)
- Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yinxuan Song
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Fangwen Peng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Jianhao Li
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Lingzhi Zhang
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Yingqi Tian
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Zhen Luo
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China.,Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, 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
|
Li Z. Time-reversal symmetry adaptation in relativistic density matrix renormalization group algorithm. J Chem Phys 2023; 158:044119. [PMID: 36725514 DOI: 10.1063/5.0127621] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
In the nonrelativistic Schrödinger equation, the total spin S and spin projection M are good quantum numbers. In contrast, spin symmetry is lost in the presence of spin-dependent interactions, such as spin-orbit couplings in relativistic Hamiltonians. Therefore, the relativistic density matrix renormalization group algorithm (R-DMRG) only employing particle number symmetry is much more expensive than nonrelativistic DMRG. In addition, artificial breaking of Kramers degeneracy can happen in the treatment of systems with an odd number of electrons. To overcome these issues, we propose time-reversal symmetry adaptation for R-DMRG. Since the time-reversal operator is antiunitary, this cannot be simply achieved in the usual way. We introduce a time-reversal symmetry-adapted renormalized basis and present strategies to maintain the structure of basis functions during the sweep optimization. With time-reversal symmetry adaptation, only half of the renormalized operators are needed, and the computational costs of Hamiltonian-wavefunction multiplication and renormalization are reduced by half. The present construction of the time-reversal symmetry-adapted basis also directly applies to other tensor network states without loops.
Collapse
Affiliation(s)
- Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
13
|
Beran P, Pernal K, Pavošević F, Veis L. Projection-Based Density Matrix Renormalization Group in Density Functional Theory Embedding. J Phys Chem Lett 2023; 14:716-722. [PMID: 36648273 PMCID: PMC10017021 DOI: 10.1021/acs.jpclett.2c03298] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 01/11/2023] [Indexed: 06/17/2023]
Abstract
The density matrix renormalization group (DMRG) method has already proved itself as a very efficient and accurate computational method, which can treat large active spaces and capture the major part of strong correlation. Its application on larger molecules is, however, limited by its own computational scaling as well as demands of methods for treatment of the missing dynamical electron correlation. In this work, we present the first step in the direction of combining DMRG with density functional theory (DFT), one of the most employed quantum chemical methods with favorable scaling, by means of the projection-based wave function (WF)-in-DFT embedding. On two proof-of-concept but important molecular examples, we demonstrate that the developed DMRG-in-DFT approach provides a very accurate description of molecules with a strongly correlated fragment.
Collapse
Affiliation(s)
- Pavel Beran
- J.
Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223Prague 8, Czech Republic
- Faculty
of Mathematics and Physics, Charles University, 121 16Prague, Czech Republic
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005Lodz, Poland
| | - Fabijan Pavošević
- Center
for Computational Quantum Physics, Flatiron
Institute, 162 Fifth Avenue, New York, 10010New York, United
States
| | - Libor Veis
- J.
Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223Prague 8, Czech Republic
| |
Collapse
|
14
|
Tian Y, Xie Z, Luo Z, Ma H. Mixed-Precision Implementation of the Density Matrix Renormalization Group. J Chem Theory Comput 2022; 18:7260-7271. [PMID: 36322159 DOI: 10.1021/acs.jctc.2c00632] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The mixed-precision optimization is an effective emerging technique for quantum chemistry methods to obtain better computational performance and maintain the chemical accuracy. Here, we developed a two-level mixed-precision implementation for the density matrix renormalization group (DMRG) method. This implementation is based on the idea that the DMRG process is an iterative process. Therefore, the first several iteration steps can be executed in single precision. A feasible single-precision DMRG may generate moderate accuracy, and when a few double-precision cleanup sweeps are added, the double-precision accuracy will be recovered. In the double-precision sweeps, we developed a mixed-precision diagonalization method that can run the most time-consuming step in single precision and maintain the double-precision accuracy. By combining these two mixed-precision schemes, we implemented our mixed-precision DMRG method. The benchmark result shows that our mixed-precision implementation achieved a good performance. A speed-up of up to 2.31 is achieved, and the accuracy is preserved within 0.01 kcal/mol.
Collapse
Affiliation(s)
- Yingqi Tian
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Zhen Luo
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
15
|
Drwal D, Beran P, Hapka M, Modrzejewski M, Sokół A, Veis L, Pernal K. Efficient Adiabatic Connection Approach for Strongly Correlated Systems: Application to Singlet-Triplet Gaps of Biradicals. J Phys Chem Lett 2022; 13:4570-4578. [PMID: 35580342 PMCID: PMC9150121 DOI: 10.1021/acs.jpclett.2c00993] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 05/09/2022] [Indexed: 06/15/2023]
Abstract
Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The ACn method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, ACn is more efficient than existing ab initio multireference dynamic correlation methods. ACn affords excellent results for singlet-triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.
Collapse
Affiliation(s)
- Daria Drwal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Pavel Beran
- J.
Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
- Faculty
of Mathematics and Physics, Charles University, Ke Karlovu 2027/3, 12116 Prague 2, Czech Republic
| | - Michał Hapka
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Adam Sokół
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Libor Veis
- J.
Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Katarzyna Pernal
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|
16
|
Li W, Ren J, Yang H, Shuai Z. On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:254003. [PMID: 35378514 DOI: 10.1088/1361-648x/ac640e] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2022] [Accepted: 04/04/2022] [Indexed: 06/14/2023]
Abstract
Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, includingab initioelectronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term 'on the fly swapping' (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model,ab initioelectronic structure of N2molecule, and the S1/S2internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.
Collapse
Affiliation(s)
- Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Hengrui Yang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing, People's Republic of China
| |
Collapse
|