1
|
Gałyńska M, Boguslawski K. Benchmarking Ionization Potentials from pCCD Tailored Coupled Cluster Models. J Chem Theory Comput 2024; 20:4182-4195. [PMID: 38752491 PMCID: PMC11137826 DOI: 10.1021/acs.jctc.4c00172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 05/04/2024] [Accepted: 05/07/2024] [Indexed: 05/29/2024]
Abstract
The ionization potential (IP) is an important parameter providing essential insights into the reactivity of chemical systems. IPs are also crucial for designing, optimizing, and understanding the functionality of modern technological devices. We recently showed that limiting the CC ansatz to the seniority-zero sector proves insufficient in predicting reliable and accurate ionization potentials within an IP equation-of-motion coupled-cluster formalism. Specifically, the absence of dynamical correlation in the seniority-zero pair coupled cluster doubles (pCCD) model led to unacceptably significant errors of approximately 1.5 eV. In this work, we aim to explore the impact of dynamical correlation and the choice of the molecular orbital basis (canonical vs localized) in CC-type methods targeting 230 ionized states in 70 molecules, comprising small organic molecules, medium-sized organic acceptors, and nucleobases. We focus on pCCD-based approaches as well as the conventional IP-EOM-CCD and IP-EOM-CCSD. Their performance is compared to the CCSD(T) or CCSDT equivalent and experimental reference data. Our statistical analysis reveals that all investigated frozen-pair coupled cluster methods exhibit similar performance, with differences in errors typically within chemical accuracy (1 kcal/mol or 0.05 eV). Notably, the effect of the molecular orbital basis, such as canonical Hartree-Fock or natural pCCD-optimized orbitals, on the IPs is marginal if dynamical correlation is accounted for. Our study suggests that triple excitations are crucial in achieving chemical accuracy in IPs when modeling electron detachment processes with pCCD-based methods.
Collapse
Affiliation(s)
- Marta Gałyńska
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| |
Collapse
|
2
|
Schraivogel T, Kats D. Two determinant distinguishable cluster. J Chem Phys 2024; 160:124109. [PMID: 38526108 DOI: 10.1063/5.0199274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 03/06/2024] [Indexed: 03/26/2024] Open
Abstract
A two reference determinant version of the distinguishable cluster with singles and doubles (DCSD) has been developed. We have implemented the two determinant distinguishable cluster (2D-DCSD) and the corresponding traditional 2D-CCSD method in a new open-source package written in Julia called ElemCo.jl. The methods were benchmarked on singlet and triplet excited states of valence and Rydberg character, as well as for singlet-triplet gaps of diradicals. It is demonstrated that the distinguishable cluster approximation improves the accuracy of 2D-CCSD.
Collapse
Affiliation(s)
- Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
3
|
Hait D, Head-Gordon M. When Is a Bond Broken? The Polarizability Perspective. Angew Chem Int Ed Engl 2023; 62:e202312078. [PMID: 37713599 DOI: 10.1002/anie.202312078] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 09/15/2023] [Accepted: 09/15/2023] [Indexed: 09/17/2023]
Abstract
The question of when a chemical bond can be said to be broken is of fundamental chemical interest but has not been widely studied. Herein we propose that the maxima of static polarizability along bond dissociation coordinates naturally define cutoff points for bond rupture, as they represent the onset of localization of shared electron density into constituent fragments. Examples of computed polarizability maxima over the course of bond cleavage in main-group and transition metal compounds are provided, across covalent, dative and charge-shift bonds. The behavior along reaction paths is also considered. Overall, the static polarizability is found to be a sensitive reporter of electronic structure reorganization associated with bond stretching, and thus can serve as a metric for describing bond cleavage (or diagnose the absence of a chemical bond).
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, CA 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, CA 94305, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, Berkeley, CA 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| |
Collapse
|
4
|
Tecmer P, Boguslawski K. Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry. Phys Chem Chem Phys 2022; 24:23026-23048. [PMID: 36149376 DOI: 10.1039/d2cp02528k] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.
Collapse
Affiliation(s)
- Paweł Tecmer
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland.
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
5
|
Abstract
The pair coupled cluster doubles (pCCD) ansatz represents an inexpensive but accurate single-reference method to describe multi-reference problems. By construction, pCCD remains, however, applicable to closed-shell systems. For the first time, we present extensions to pCCD that allow us to target open-shell molecules with up to 4 unpaired electrons. Although requiring only modest computational cost, our methods approach chemical accuracy for some challenging cases, while their performance is comparable to more expensive models like DMRG or CCSD(T).
Collapse
Affiliation(s)
- Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziądzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
6
|
Chattopadhyay S. Single-Root Multireference Brillouin-Wigner Perturbative Approach to Excitation Energies. ACS OMEGA 2021; 6:1668-1686. [PMID: 33490826 PMCID: PMC7818614 DOI: 10.1021/acsomega.0c05714] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 12/28/2020] [Indexed: 06/12/2023]
Abstract
The state-specific Brillouin-Wigner multireference perturbation theory [which employs Jeziorski-Monkhorst parametrization of the wave function] using improved virtual orbitals, denoted as IVO-BWMRPT, is applied to calculate excitation energies (EEs) for methylene, ethylene, trimethylenemethane, and benzyne systems exhibiting various degrees of diradical character. In IVO-BWMRPT, all of the parameters appearing in the wave function ansatz are optimized for a specific electronic state. For these systems, the IVO-BWMRPT method provides EEs that are in close agreement with the benchmark results and experiments, where available, indicating that the method does not introduce imbalance in the target-specific treatment of closed- and open-shell states involved. The good performance of the present methodology is primarily related to structural compactness of the formalism. Overall, present findings are encouraging for both further development of the approach and chemical applications on the energy differences of strongly correlated systems.
Collapse
Affiliation(s)
- Sudip Chattopadhyay
- Department of Chemistry, Indian Institute of Engineering Science and Technology,
Shibpur, Howrah 711103, India
| |
Collapse
|
7
|
Bartlett RJ, Park YC, Bauman NP, Melnichuk A, Ranasinghe D, Ravi M, Perera A. Index of multi-determinantal and multi-reference character in coupled-cluster theory. J Chem Phys 2020; 153:234103. [PMID: 33353328 DOI: 10.1063/5.0029339] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
A full configuration interaction calculation (FCI) ultimately defines the innate molecular orbital description of a molecule. Its density matrix and the natural orbitals obtained from it quantify the difference between having N-dominantly occupied orbitals in a reference determinant for a wavefunction to describe N-correlated electrons and how many of those N-electrons are left to the remaining virtual orbitals. The latter provides a measure of the multi-determinantal character (MDC) required to be in a wavefunction. MDC is further split into a weak correlation part and a part that indicates stronger correlation often called multi-reference character (MRC). If several virtual orbitals have high occupation numbers, then one might argue that these additional orbitals should be allowed to have a larger role in the calculation, as in MR methods, such as MCSCF, MR-CI, or MR-coupled-cluster (MR-CC), to provide adequate approximations toward the FCI. However, there are problems with any of these MR methods that complicate the calculations compared to the uniformity and ease of application of single-reference CC calculations (SR-CC) and their operationally single-reference equation-of-motion (EOM-CC) extensions. As SR-CC theory is used in most of today's "predictive" calculations, an assessment of the accuracy of SR-CC at some truncation of the cluster operator would help to quantify how large an issue MRC actually is in a calculation, and how it might be alleviated while retaining the convenient SR computational character of CC/EOM-CC. This paper defines indices that identify MRC situations and help assess how reliable a given calculation is.
Collapse
Affiliation(s)
- Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Nicholas P Bauman
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ann Melnichuk
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Duminda Ranasinghe
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Moneesha Ravi
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
8
|
Ranasinghe DS, Margraf JT, Perera A, Bartlett RJ. Vertical valence ionization potential benchmarks from equation-of-motion coupled cluster theory and QTP functionals. J Chem Phys 2019; 150:074108. [DOI: 10.1063/1.5084728] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Duminda S. Ranasinghe
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Johannes T. Margraf
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
9
|
Lutz JJ, Nooijen M, Perera A, Bartlett RJ. Reference dependence of the two-determinant coupled-cluster method for triplet and open-shell singlet states of biradical molecules. J Chem Phys 2018; 148:164102. [DOI: 10.1063/1.5025170] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Jesse J. Lutz
- Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio 45433, USA
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
10
|
Margraf JT, Perera A, Lutz JJ, Bartlett RJ. Single-reference coupled cluster theory for multi-reference problems. J Chem Phys 2017; 147:184101. [DOI: 10.1063/1.5003128] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Johannes T. Margraf
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Jesse J. Lutz
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
11
|
Bazante AP, Perera A, Bartlett RJ. Towards core-excitation spectra in attosecond spectroscopy: A coupled-cluster study of ClF. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.05.017] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
12
|
Zimmerman PM. Singlet–Triplet Gaps through Incremental Full Configuration Interaction. J Phys Chem A 2017; 121:4712-4720. [DOI: 10.1021/acs.jpca.7b03998] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Paul M. Zimmerman
- Department of Chemistry, University of Michigan 930 North University
Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
13
|
A comparative DFT study of interactions of Au and small gold clusters Aun (n = 2–4) with CH3S and CH2 radicals. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.01.019] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
14
|
Lefrancois D, Rehn DR, Dreuw A. Accurate adiabatic singlet-triplet gaps in atoms and molecules employing the third-order spin-flip algebraic diagrammatic construction scheme for the polarization propagator. J Chem Phys 2016; 145:084102. [DOI: 10.1063/1.4961298] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Dirk R. Rehn
- Departments of Physics, Chemistry and Biology, Linköping University, SE-581 83 Linköping, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
15
|
Sinha Ray S, Ghosh A, Chattopadhyay S, Chaudhuri RK. Taming the Electronic Structure of Diradicals through the Window of Computationally Cost Effective Multireference Perturbation Theory. J Phys Chem A 2016; 120:5897-916. [DOI: 10.1021/acs.jpca.6b03211] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Suvonil Sinha Ray
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Anirban Ghosh
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Sudip Chattopadhyay
- Department of Chemistry, Indian Institute of Engineering Science and Technology, Shibpur, Howrah 711103, India
| | - Rajat K. Chaudhuri
- Theoretical Physics, Indian Institute of Astrophysics, Bangalore 560034, India
| |
Collapse
|
16
|
Evaluation of modern DFT functionals and G3n-RAD composite methods in the modelization of organic singlet diradicals. J Mol Model 2016; 22:76. [DOI: 10.1007/s00894-016-2950-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2015] [Accepted: 02/29/2016] [Indexed: 11/26/2022]
|
17
|
Hansen JA, Bauman NP, Shen J, Borden WT, Piecuch P. Ab initiocoupled-cluster and multi-reference configuration interaction studies of the low-lying electronic states of 1,2,3,4-cyclobutanetetraone. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1112926] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
18
|
Liu J, Bao X, Hrovat DA, Borden WT. Theoretical Analysis of the Fragmentation of (CO)5: A Symmetry-Allowed Highly Exothermic Reaction that Follows a Stepwise Pathway. J Org Chem 2015; 80:11788-93. [PMID: 26375997 DOI: 10.1021/acs.joc.5b01546] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
B3LYP and CCSD(T) calculations, using an aug-cc-pVTZ basis set, have been carried out on the fragmentation of 1,2,3,4,5-cyclopentanepentone, (CO)(5), to five molecules of CO. Although this reaction is calculated to be highly exothermic and is allowed to be concerted by the Woodward-Hoffmann rules, our calculations find that the D(5h) energy maximum is a multidimensional hilltop on the potential energy surface. This D(5h) hilltop is 16-20 kcal/mol higher in energy than a C(2) transition structure for the endothermic cleavage of (CO)(5) to (CO)(4) + CO and 11-15 kcal/mol higher than a C(s) transition structure for the loss of two CO molecules. The reasons for the very high energy of the D(5h) hilltop are discussed, and the geometries of the two lower energy transition structures are rationalized on the basis of mixing of the e(2)' HOMO and the a(2)″ LUMO of the hilltop.
Collapse
Affiliation(s)
- Jiajun Liu
- College of Chemistry, Chemical Engineering and Materials Science, Soochow University , 199 Ren-Ai Road, Suzhou Industrial Park, Suzhou, Jiangsu 215123, China
| | - Xiaoguang Bao
- College of Chemistry, Chemical Engineering and Materials Science, Soochow University , 199 Ren-Ai Road, Suzhou Industrial Park, Suzhou, Jiangsu 215123, China
| | - David A Hrovat
- Department of Chemistry and the Center for Advanced Scientific Computing and Modeling, University of North Texas , 1155 Union Circle, #305070, Denton, Texas 76203-5017, United States
| | - Weston Thatcher Borden
- Department of Chemistry and the Center for Advanced Scientific Computing and Modeling, University of North Texas , 1155 Union Circle, #305070, Denton, Texas 76203-5017, United States
| |
Collapse
|
19
|
|
20
|
Borden WT. The Unpredictability of Research Directions and the Synergy between Theory and Experiment in Physical-Organic Chemistry. Isr J Chem 2015. [DOI: 10.1002/ijch.201400170] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
21
|
Chen B, Hrovat DA, West R, Deng SHM, Wang XB, Borden WT. The Negative Ion Photoelectron Spectrum of Cyclopropane-1,2,3-Trione Radical Anion, (CO)3•– — A Joint Experimental and Computational Study. J Am Chem Soc 2014; 136:12345-54. [DOI: 10.1021/ja505582k] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bo Chen
- Department
of Chemistry and the Center for Advanced, Scientific Computing and
Modeling, University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5070, United States
| | - David A. Hrovat
- Department
of Chemistry and the Center for Advanced, Scientific Computing and
Modeling, University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5070, United States
| | - Robert West
- Department
of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Shihu H. M. Deng
- Physical
Sciences Division, Pacific Northwest National Laboratory, P.O. Box
999, MS K8-88, Richland, Washington 99352, United States
| | - Xue-Bin Wang
- Physical
Sciences Division, Pacific Northwest National Laboratory, P.O. Box
999, MS K8-88, Richland, Washington 99352, United States
| | - Weston Thatcher Borden
- Department
of Chemistry and the Center for Advanced, Scientific Computing and
Modeling, University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5070, United States
| |
Collapse
|