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
|
Tecmer P, Gałyńska M, Szczuczko L, Boguslawski K. Geminal-Based Strategies for Modeling Large Building Blocks of Organic Electronic Materials. J Phys Chem Lett 2023; 14:9909-9917. [PMID: 37903084 PMCID: PMC10641881 DOI: 10.1021/acs.jpclett.3c02434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 10/09/2023] [Accepted: 10/24/2023] [Indexed: 11/01/2023]
Abstract
We elaborate on unconventional electronic structure methods based on geminals and their potential to advance the rapidly developing field of organic photovoltaics (OPVs). Specifically, we focus on the computational advantages of geminal-based methods over standard approaches and identify the critical aspects of OPV development. Examples are reliable and efficient computations of orbital energies, electronic spectra, and van der Waals interactions. Geminal-based models can also be combined with quantum embedding techniques and a quantum information analysis of orbital interactions to gain a fundamental understanding of the electronic structures and properties of realistic OPV building blocks. Furthermore, other organic components present in, for instance, dye-sensitized solar cells (DSSCs) represent another promising scope of application. Finally, we provide numerical examples predicting the properties of a small building block of OPV components and two carbazole-based dyes proposed as possible DSSC sensitizers.
Collapse
Affiliation(s)
- Paweł Tecmer
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Marta Gałyńska
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lena Szczuczko
- 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
|
3
|
Mamache S, Gałyńska M, Boguslawski K. Benchmarking ionization potentials using the simple pCCD model. Phys Chem Chem Phys 2023. [PMID: 37378457 DOI: 10.1039/d3cp01963b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
The electron-detachment energy is measured by the ionization potential (IP). As a result, it is a fundamental, observable and important molecular electronic signature in photoelectron spectroscopy. A precise theoretical prediction of electron-detachment energies or ionization potentials is essential for organic optoelectronic systems like transistors, solar cells, or light-emitting diodes. In this work, we benchmark the performance of the recently presented IP variant of the equation-of-motion pair coupled cluster doubles (IP-EOM-pCCD) model to determine IPs. Specifically, the predicted ionization energies are compared to experimental results and higher-order coupled cluster theories based on statistically assessing 201 electron-detached states of 41 organic molecules for three different molecular orbital basis sets and two sets of particle-hole operators. While IP-EOM-pCCD features a reasonable spread and skewness of ionization energies, its mean error and standard deviation differ by up to 1.5 eV from reference data. Our study, thus, highlights the importance of dynamical correlation to reliably predict IPs from a pCCD reference function in small organic molecules.
Collapse
Affiliation(s)
- Saddem Mamache
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - 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
|
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
|
Brzeski J, Jordan KD. Non-Valence Anions of Pyridine and the Diazines. J Phys Chem A 2022; 126:5310-5313. [PMID: 35920853 DOI: 10.1021/acs.jpca.2c04040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The dipole-bound anions of pyridine, pyridazine, and pyrimidine are characterized using equation of motion coupled cluster singles and doubles calculations. These calculations predict that the anions of pyridine, pyrimidine and pyridazine are bound in the Born-Oppenheimer approximation by 0.05, 0.8, and 19.0 meV, respectively. The binding energies of pyrimidine and pyridazine are large enough that the anions will remain bound even when allowing for corrections to the Born-Oppenheimer approximation, while that of pyridine is a borderline case. We were unable to find a stable non-valence correlation-bound anion for pyrazine, which has a zero dipole moment.
Collapse
Affiliation(s)
- Jakub Brzeski
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15218, United States.,Department of Bioinorganic Chemistry, Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Kenneth D Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15218, United States
| |
Collapse
|
6
|
Liebenthal M, Vu NH, DePrince E. Equation-of-motion cavity quantum electrodynamics coupled-cluster theory for electron attachment. J Chem Phys 2022; 156:054105. [DOI: 10.1063/5.0078795] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Nam Hoang Vu
- Chemistry & Biochemistry, Florida State University, United States of America
| | - Eugene DePrince
- Chemistry and Biochemistry, Florida State University, United States of America
| |
Collapse
|
7
|
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
|
8
|
Lu Y, Guo M, Wang Z, Wang F. Low-lying states of Tl2 and Nh2 with EOM-CC and FSCC methods. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
9
|
Guo M, Wang Z, Wang F. Equation-of-motion coupled-cluster theory for double electron attachment with spin-orbit coupling. J Chem Phys 2020; 153:214118. [PMID: 33291924 DOI: 10.1063/5.0032716] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report implementation of the equation-of-motion coupled-cluster (EOM-CC) method for double electron-attachment (DEA) with spin-orbit coupling (SOC) at the CC singles and doubles (CCSD) level using a closed-shell reference in this work. The DEA operator employed in this work contains two-particle and three-particle one-hole excitations, and SOC is included in post-Hartree-Fock treatment. Time-reversal symmetry and spatial symmetry are exploited to reduce computational cost. The EOM-DEA-CCSD method with SOC allows us to investigate SOC effects of systems with two-unpaired electrons. According to our results on atoms, double ionization potentials (DIPs), excitation energies (EEs), and SO splittings of low-lying states are calculated reliably using the EOM-DEA-CCSD method with SOC. Its accuracy is usually higher than that of EOM-CCSD for EEs or DIPs if the same target can be reached from single excitations by choosing a proper closed-shell reference. However, performance of the EOM-DEA-CCSD method with SOC on molecules is not as good as that for atoms. Bond lengths for the ground and the several lowest excited states of GaH, InH, and TlH are underestimated pronouncedly, although reasonable EEs are obtained, and splittings of the 3Σ- state from the π2 configuration are calculated to be too small with EOM-DEA-CCSD.
Collapse
Affiliation(s)
- Minggang Guo
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| | - Zhifan Wang
- College of Chemistry and Life Science, Chengdu Normal University, Chengdu 611130, People's Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| |
Collapse
|
10
|
Guo M, Wang Z, Wang F. Treating spin-orbit coupling at different levels in equation-of-motion coupled-cluster calculations. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1785029] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Minggang Guo
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People’s Republic of China
| | - Zhifan Wang
- College of Chemistry and Life Science, Chengdu Normal University, Chengdu, People’s Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People’s Republic of China
| |
Collapse
|
11
|
Zhao H, Wang Z, Guo M, Wang F. Splittings of d 8 configurations of late-transition metals with EOM-DIP-CCSD and FSCCSD methods. J Chem Phys 2020; 152:134105. [PMID: 32268764 DOI: 10.1063/1.5145077] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Multireference methods are usually required for transition metal systems due to the partially filled d electrons. In this work, the single-reference equation-of-motion coupled-cluster method at the singles and doubles level for double ionization potentials (EOM-DIP-CCSD) is employed to calculate energies of states from the d8 configuration of late-transition metal atoms starting from a closed-shell reference. Its results are compared with those from the multireference Fock-space coupled-cluster method at the CCSD level (FSCCSD) for DIP from the same closed-shell reference. Both scalar-relativistic effects and spin-orbit coupling are considered in these calculations. Compared with all-electron FSCCSD results with four-component Dirac-Coulomb Hamiltonian, FSCCSD with relativistic effective core potentials can provide reasonable results, except for atoms with unstable reference. Excitation energies for states in the (n - 1)d8ns2 configuration are overestimated pronouncedly with these two methods, and this overestimation is more severe than those in the (n - 1)d9ns1 configuration. Error of EOM-CCSD on these excitation energies is generally larger than that of FSCCSD. On the other hand, relative energies of most of the states in the d8 configuration with respect to the lowest state in the same configuration are predicted reliably with EOM-DIP-CCSD, except for the 3P0 state of Hg2+ and states in Ir+. FSCCSD can provide reasonable relative energies for the several lowest states, while its error tends to be larger for higher states.
Collapse
Affiliation(s)
- Hewang Zhao
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| | - Zhifan Wang
- School of Electronic Engineering, Chengdu Technological University, Chengdu 611730, China
| | - Minggang Guo
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| |
Collapse
|
12
|
Pathak H, Sasmal S, Talukdar K, Nayak MK, Vaval N, Pal S. Relativistic double-ionization equation-of-motion coupled-cluster method: Application to low-lying doubly ionized states. J Chem Phys 2020; 152:104302. [PMID: 32171231 DOI: 10.1063/1.5140988] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This article deals with the extension of the relativistic double-ionization equation-of-motion coupled-cluster (DI-EOMCC) method [H. Pathak et al. Phys. Rev. A 90, 010501(R) (2014)] for the molecular systems. The Dirac-Coulomb Hamiltonian with four-component spinors is considered to take care of the relativistic effects. The implemented method is employed to compute a few low-lying doubly ionized states of noble gas atoms (Ar, Kr, Xe, and Rn) and Cl2, Br2, HBr, and HI. Additionally, we presented results with two intermediate schemes in the four-component relativistic DI-EOMCC framework to understand the role of electron correlation. The computed double ionization spectra for the atomic systems are compared with the values from the non-relativistic DI-EOMCC method with spin-orbit coupling [Z. Wang et al. J. Chem. Phys. 142, 144109 (2015)] and the values from the National Institute of Science and Technology (NIST) database. Our atomic results are found to be in good agreement with the NIST values. Furthermore, the obtained results for the molecular systems agree well with the available experimental values.
Collapse
Affiliation(s)
- Himadri Pathak
- Electronic Structure Theory Group, Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune 411008, India
| | - Sudip Sasmal
- Electronic Structure Theory Group, Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune 411008, India
| | - Kaushik Talukdar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Nayana Vaval
- Electronic Structure Theory Group, Physical Chemistry Division, CSIR-National Chemical Laboratory, Pune 411008, India
| | - Sourav Pal
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
13
|
Cao Z, Wang F, Yang M. Spin-orbit coupling with approximate equation-of-motion coupled-cluster method for ionization potential and electron attachment. J Chem Phys 2016; 145:154110. [DOI: 10.1063/1.4964859] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Zhanli Cao
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People’s Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People’s Republic of China
| | - Mingli Yang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People’s Republic of China
| |
Collapse
|
14
|
Musiał M, Lupa Ł, Kucharski SA. Equation-of-motion coupled cluster method for the description of the high spin excited states. J Chem Phys 2016; 144:154105. [PMID: 27389207 DOI: 10.1063/1.4946031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The equation-of-motion (EOM) coupled cluster (CC) approach in the version applicable for the excitation energy(EE) calculations has been formulated for high spin components. The EE-EOM-CC scheme based on the restricted Hartree-Fock reference and standard amplitude equations as used in the Davidson diagonalization procedure yields the singlet states. The triplet and higher spin components require separate amplitude equations. In the case of quintets, the relevant equations are much simpler and easier to solve. Out of 26 diagrammatic terms contributing to the R1 and R2 singlet equations in the case of quintets, only R2 operator survives with 5 diagrammatic terms present. In addition all terms engaging three body elements of the similarity transformed Hamiltonian disappear. This indicates a substantial simplification of the theory. The implemented method has been applied to the pilot study of the excited states of the C2 molecule and quintet states of C and Si atoms.
Collapse
Affiliation(s)
- Monika Musiał
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Łukasz Lupa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | | |
Collapse
|
15
|
Wang Z, Hu S, Wang F, Guo J. Equation-of-motion coupled-cluster method for doubly ionized states with spin-orbit coupling. J Chem Phys 2015; 142:144109. [DOI: 10.1063/1.4917041] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Zhifan Wang
- College of Chemistry, Sichuan University, Chengdu 610064, People’s Republic of China
| | - Shu Hu
- Key Laboratory of Chemical Laser, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610064, People’s Republic of China
| | - Jingwei Guo
- Key Laboratory of Chemical Laser, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| |
Collapse
|