1
|
Li C, Mao S, Huang R, Evangelista FA. Frozen Natural Orbitals for the State-Averaged Driven Similarity Renormalization Group. J Chem Theory Comput 2024; 20:4170-4181. [PMID: 38747709 DOI: 10.1021/acs.jctc.4c00152] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
We present a reduced-cost implementation of the state-averaged driven similarity renormalization group (SA-DSRG) based on the frozen natural orbital (FNO) approach. The natural orbitals (NOs) are obtained by diagonalizing the one-body reduced density matrix from SA-DSRG second-order perturbation theory (SA-DSRG-PT2). We consider three criteria to truncate the virtual NOs for the subsequent electron correlation treatment beyond SA-DSRG-PT2. An additive second-order correction is applied to the SA-DSRG Hamiltonian to reintroduce correlation effects from the discarded orbitals. The FNO SA-DSRG method is benchmarked on 35 small organic molecules in the QUEST database. When keeping 98-99% of the cumulative occupation numbers, the mean absolute error in the vertical transition energies due to FNO is less than 0.01 eV. Using the same FNO threshold, we observe a speedup of 9 times compared to the conventional SA-DSRG implementation for nickel carbonyl with a quadruple-ζ basis set. The FNO approach enables nonperturbative SA-DSRG computations on chloroiron corrole [FeCl(C19H11N4)] with more than 1000 basis functions, surpassing the current limit of a conventional implementation.
Collapse
Affiliation(s)
- Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Shuxian Mao
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Renke Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
2
|
Li L, Prindle CR, Shi W, Nuckolls C, Venkataraman L. Radical Single-Molecule Junctions. J Am Chem Soc 2023; 145:18182-18204. [PMID: 37555594 DOI: 10.1021/jacs.3c04487] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
Radicals are unique molecular systems for applications in electronic devices due to their open-shell electronic structures. Radicals can function as good electrical conductors and switches in molecular circuits while also holding great promise in the field of molecular spintronics. However, it is both challenging to create stable, persistent radicals and to understand their properties in molecular junctions. The goal of this Perspective is to address this dual challenge by providing design principles for the synthesis of stable radicals relevant to molecular junctions, as well as offering current insight into the electronic properties of radicals in single-molecule devices. By exploring both the chemical and physical properties of established radical systems, we will facilitate increased exploration and development of radical-based molecular systems.
Collapse
Affiliation(s)
- Liang Li
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Claudia R Prindle
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Wanzhuo Shi
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Colin Nuckolls
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Latha Venkataraman
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Department of Applied Physics and Applied Mathematics, Columbia University, New York, New York 10027, United States
| |
Collapse
|
3
|
Krstić M, Fink K, Sharapa DI. The Adsorption of Small Molecules on the Copper Paddle-Wheel: Influence of the Multi-Reference Ground State. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27030912. [PMID: 35164179 PMCID: PMC8840508 DOI: 10.3390/molecules27030912] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 01/26/2022] [Accepted: 01/26/2022] [Indexed: 11/16/2022]
Abstract
We report a theoretical study of the adsorption of a set of small molecules (C2H2, CO, CO2, O2, H2O, CH3OH, C2H5OH) on the metal centers of the “copper paddle-wheel”—a key structural motif of many MOFs. A systematic comparison between DFT of different rungs, single-reference post-HF methods (MP2, SOS–MP2, MP3, DLPNO–CCSD(T)), and multi-reference approaches (CASSCF, DCD–CAS(2), NEVPT2) is performed in order to find a methodology that correctly describes the complicated electronic structure of paddle-wheel structure together with a reasonable description of non-covalent interactions. Apart from comparison with literature data (experimental values wherever possible), benchmark calculations with DLPNO–MR–CCSD were also performed. Despite tested methods show qualitative agreement in the majority of cases, we showed and discussed reasons for quantitative differences as well as more fundamental problems of specific cases.
Collapse
Affiliation(s)
- Marjan Krstić
- Institute for Theoretical Solid State Physics (TFP), Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany;
| | - Karin Fink
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany;
| | - Dmitry I. Sharapa
- Institute of Catalysis Research and Technology (IKFT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Correspondence:
| |
Collapse
|
4
|
Arikawa S, Shimizu A, Shiomi D, Sato K, Shintani R. Synthesis and Isolation of a Kinetically Stabilized Crystalline Triangulene. J Am Chem Soc 2021; 143:19599-19605. [PMID: 34767718 DOI: 10.1021/jacs.1c10151] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The synthesis and isolation of hydrocarbons with a triplet ground state in crystalline forms have been sought in materials science. Triangulene is one of the most famous triplet-ground-state benzenoid hydrocarbons. Its unique electronic structure and highly symmetric structure have prompted many scientists to synthesize and isolate triangulene and its derivatives, but all attempts so far to isolate them as crystals have been unsuccessful. Herein we report the synthesis and isolation of a kinetically stabilized crystalline triangulene for the first time. The key to success is the introduction of bulky substituents onto the reactive zigzag edges. Its highly symmetric structure was confirmed by X-ray crystallography, and its fundamental properties, including the triplet ground state, were revealed. The achievement here will open the door for the synthesis and isolation of other hydrocarbons with higher spin multiplicity.
Collapse
Affiliation(s)
- Shinobu Arikawa
- Division of Chemistry, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Akihiro Shimizu
- Division of Chemistry, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Daisuke Shiomi
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Kazunobu Sato
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Ryo Shintani
- Division of Chemistry, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
5
|
Lang J, Antalík A, Veis L, Brandejs J, Brabec J, Legeza Ö, Pittner J. Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T). J Chem Theory Comput 2020; 16:3028-3040. [PMID: 32275424 DOI: 10.1021/acs.jctc.0c00065] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Sciences, Charles University, Albertov 6, 128 00 Prague 2, Czech Republic
| | - Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jan Brandejs
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
6
|
Saitow M, Yanai T. A multireference coupled-electron pair approximation combined with complete-active space perturbation theory in local pair-natural orbital framework. J Chem Phys 2020; 152:114111. [PMID: 32199413 DOI: 10.1063/1.5142622] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Complete-Active Space Second-order Perturbation Theory (CASPT2) has been one of the most widely-used methods for reliably calculating electronic structures of multireference systems. Because of its lowest level treatment of dynamic correlation, it has a high computational feasibility; however, its accuracy in some cases falls short of needs. Here, as a simple yet higher-order alternative, we introduce a hybrid theory of the CASPT2 and a multireference variant of the Coupled-Electron Pair Approximation (CEPA), which is a class of high level correlation theory. A central feature of our theory (CEPT2) is to use the two underlying theories for describing different divisions of correlation components based on the full internal contraction framework. The external components, which usually give a major contribution to the dynamic correlation, are intensively described using the CEPA Ansatz, while the rests are treated at the CASPT2 level. Furthermore, to drastically reduce the computational demands, we have incorporated the pair-natural orbital (PNO) method into our multireference implementations. This development, thus, requires highly complex derivations and coding, while it has been largely facilitated with an automatic expression and code generation technique. To highlight the accuracy of the CEPT2 approach and to assess the errors caused by the PNO truncation, benchmark calculations are shown on small- to medium-size molecules, illustrating the high accuracy of the present CEPT2 model. By tightening the truncation thresholds, the PNO-CEPT2 energy converges toward the canonical counterpart and is more accurate than that of PNO-CASPT2 as long as the same truncation thresholds are used.
Collapse
Affiliation(s)
- Masaaki Saitow
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
7
|
Zhou R, Liang Q, Yang J. Complete OSV-MP2 Analytical Gradient Theory for Molecular Structure and Dynamics Simulations. J Chem Theory Comput 2019; 16:196-210. [DOI: 10.1021/acs.jctc.9b00806] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Ruiyi Zhou
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| | - Qiujiang Liang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| | - Jun Yang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| |
Collapse
|
8
|
Miyagawa K, Kawakami T, Suzuki Y, Isobe H, Shoji M, Yamanaka S, Okumura M, Nakajima T, Yamaguchi K. Domain-based local pair natural orbital CCSD(T) calculations of strongly correlated electron systems: Examination of dynamic equilibrium models based on multiple intermediates in S1 state of photosystem II. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1666171] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- K. Miyagawa
- Institute for Scientific and Industrial Research, Osaka University, Osaka 567-0047, Japan
| | - T. Kawakami
- Graduate School of Science, Osaka University, Toyonaka 560-0043, Japan
- RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Y. Suzuki
- Graduate School of Science, Osaka University, Toyonaka 560-0043, Japan
| | - H. Isobe
- Graduate School of Natural Science and Technology, Okayama University, Okayama 700-8530, Japan
| | - M. Shoji
- Center of Computational Sciences, Tsukuba University, Tsukuba, Ibaraki 305-8577, Japan
| | - S. Yamanaka
- Graduate School of Science, Osaka University, Toyonaka 560-0043, Japan
| | - M. Okumura
- Graduate School of Science, Osaka University, Toyonaka 560-0043, Japan
| | - T. Nakajima
- RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - K. Yamaguchi
- Institute for Scientific and Industrial Research, Osaka University, Osaka 567-0047, Japan
- RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
- Insitute for Nanoscience Design, Osaka University, Toyonaka 560-0043, Japan
| |
Collapse
|
9
|
Antalík A, Veis L, Brabec J, Demel O, Legeza Ö, Pittner J. Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen). J Chem Phys 2019; 151:084112. [PMID: 31470730 DOI: 10.1063/1.5110477] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.
Collapse
Affiliation(s)
- Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research Group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
10
|
Pastorczak E, Hapka M, Veis L, Pernal K. Capturing the Dynamic Correlation for Arbitrary Spin-Symmetry CASSCF Reference with Adiabatic Connection Approaches: Insights into the Electronic Structure of the Tetramethyleneethane Diradical. J Phys Chem Lett 2019; 10:4668-4674. [PMID: 31356083 DOI: 10.1021/acs.jpclett.9b01582] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The recently proposed approach to multireference dynamic correlation energy based on the adiabatic connection (AC) is extended to an arbitrary spin symmetry of the reference state. We show that both the spin-free AC approach and its computationally inexpensive approximation, AC0, when combined with a complete active space wave function, constitute viable alternatives to the perturbation-based and density-functional-based multiconfiguration methods. In particular, the AC0 approach, thanks to its favorable scaling with the system size and the size of the active space, allows for treating larger systems than its perturbation-based counterparts while maintaining comparable accuracy. We show the method's robustness on illustrative chemical systems, including the elusive tetramethyleneethane (TME) diradical, potential energy surfaces of which present a challenge to most computational approaches. For the latter system, AC0 outperforms other methods, staying in close agreement with the full configuration interaction quantum Monte Carlo benchmark. A careful analysis of the contributions to the correlation energy of TME's lowest singlet and triplet states reveals the subtle interplay of the dynamic and static correlation as the key to understanding the shape of the diradical's potential energy surfaces.
Collapse
Affiliation(s)
- Ewa Pastorczak
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| | - Michał Hapka
- Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , 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
|