1
|
Zielinski P, Black JA, Köhn A. Performance Tests of the Second-Order Approximate Internally Contracted Multireference Coupled-Cluster Singles and Doubles Method icMRCC2. J Chem Theory Comput 2023; 19:8671-8688. [PMID: 37991987 PMCID: PMC10720349 DOI: 10.1021/acs.jctc.3c00969] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2023] [Revised: 11/03/2023] [Accepted: 11/06/2023] [Indexed: 11/24/2023]
Abstract
Benchmark results are presented for the second-order approximation of the internally contracted multireference coupled-cluster method with single and double excitations, icMRCC2 [Köhn, Bargholz, J. Chem. Phys. 2019, 151, 041106], which was designed as a multireference analogue of the single-reference second-order approximate coupled-cluster method CC2 [Christiansen, Koch, Jørgensen, Chem. Phys. Lett. 1995, 243, 409-418]. Vertical excitation energies of various small to medium-sized organic molecules are investigated based on established test sets from the literature. Additionally, the spectroscopic constants of ground and excited states of diatomics and the geometric parameters of excited triatomic molecules were determined and compared to the experimental data. The results show that the method clearly extends the applicability of single-reference CC2, including doubly excited states, and also artifacts of CC2 like too low Rydberg excitations and too weak multiple bonds are eliminated. The method is computationally more demanding than standard multireference second-order perturbation theories but improves significantly in accuracy, as shown by the benchmark results. In addition, it is demonstrated that small active spaces are often sufficient to obtain accurate energies with icMRCC2. Example applications like the automerization of cyclobutadiene, the deactivation pathway of ethylene, and the excited states of an iron complex with a noninnocent nitrosyl ligand demonstrate the potential of icMRCC2 in cases with strong multireference character.
Collapse
Affiliation(s)
- Patrik Zielinski
- Institute for Theoretical
Chemistry, University of Stuttgart, Paffenwaldring 55, Stuttgart D-70569, Germany
| | | | - Andreas Köhn
- Institute for Theoretical
Chemistry, University of Stuttgart, Paffenwaldring 55, Stuttgart D-70569, Germany
| |
Collapse
|
2
|
Hancock AC, Goerigk L. Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models. RSC Adv 2023; 13:35964-35984. [PMID: 38090083 PMCID: PMC10712016 DOI: 10.1039/d3ra07381e] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 11/21/2023] [Indexed: 05/12/2024] Open
Abstract
Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.
Collapse
Affiliation(s)
- Amy C Hancock
- School of Chemistry, The University of Melbourne Parkville Australia +61-(0)3-8344 6784
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne Parkville Australia +61-(0)3-8344 6784
| |
Collapse
|
3
|
Jacquemin D, Kossoski F, Gam F, Boggio-Pasqua M, Loos PF. Reference Vertical Excitation Energies for Transition Metal Compounds. J Chem Theory Comput 2023. [PMID: 37965941 DOI: 10.1021/acs.jctc.3c01080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.
Collapse
Affiliation(s)
- Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Franck Gam
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| |
Collapse
|
4
|
Mester D, Kállay M. Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level. J Chem Theory Comput 2023. [PMID: 37326360 DOI: 10.1021/acs.jctc.3c00363] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
5
|
Schrader T, Perlt E, Fritz T, Sierka M. Performance of Common Density Functionals for Excited States of Tetraphenyldibenzoperiflanthene. J Phys Chem A 2023; 127:3265-3273. [PMID: 37037005 DOI: 10.1021/acs.jpca.2c06715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
Abstract
Time-dependent density functional theory is the method of choice to efficiently calculate excitation spectra with the functional and basis set choice allowing one to compromise between accuracy and computational cost. In this work, the performance of different functionals as well as the second-order approximate coupled cluster singles and doubles model CC2 is evaluated by comparing the results to experimental results of the example molecule tetraphenyldibenzoperiflanthene (DBP). The choice of the functional has a significant impact on the calculated spectrum of DBP. The performance of a number of different functionals was evaluated, quantified, and, where possible, discussed. The best functional, tuned-CAM-B3LYP, is used to investigate DBP on a surface of hexagonal boron nitride (h-BN). The resulting spectrum shows excellent agreement with experimental results for a monolayer of DBP on h-BN.
Collapse
Affiliation(s)
- Tim Schrader
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Eva Perlt
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Torsten Fritz
- Institute of Solid State Physics, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Marek Sierka
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| |
Collapse
|
6
|
Mester D, Kállay M. Double-Hybrid Density Functional Theory for Core Excitations: Theory and Benchmark Calculations. J Chem Theory Comput 2023; 19:1310-1321. [PMID: 36721871 PMCID: PMC9979613 DOI: 10.1021/acs.jctc.2c01222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.
Collapse
Affiliation(s)
- Dávid Mester
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,
| |
Collapse
|
7
|
Sandoval-Salinas ME, Brémond E, Pérez-Jiménez AJ, Adamo C, Sancho-García JC. Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions. J Chem Phys 2023; 158:044105. [PMID: 36725511 DOI: 10.1063/5.0134946] [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/25/2023] Open
Abstract
A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet 1La and 1Lb states of linear oligoacenes as well as the triplet 3La and 3Lb states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S0-T1 singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.
Collapse
Affiliation(s)
- M E Sandoval-Salinas
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - E Brémond
- ITODYS, CNRS, Université Paris Cité, F-75006 Paris, France
| | - A J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - C Adamo
- Chimie ParisTech, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), UMR8060, PSL Research University, F-75005 Paris, France
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| |
Collapse
|
8
|
Helal W. Double Hybrid Density Functionals for the Electronic Excitation Energies of Linear Cyanines. J Phys Chem A 2023; 127:131-141. [PMID: 36537875 DOI: 10.1021/acs.jpca.2c07192] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The lowest bright electronic excitations of seven model linear cyanines (CN3-CN15) using 28 double-hybrid (DH) density functionals are benchmarked against accurate and recent CC3 results. Some of these DH functionals are recently designed specifically for excited electronic state calculations. In addition, CIS, CIS(D), SCS-CIS(D), and SOS-CIS(D) were also tested. Four different basis sets were used for the vertical electronic excitation calculations: cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis. Augmented basis sets (e.g. aug-cc-pVDZ and aug-cc-pVTZ) are found to be required for accurate and consistent results using DH functionals. The DH functionals tested in this work are classified into four main groups: global double-hybrids (GDH), range-separated double-hybrids (RSDH), spin-component and spin-opposite scaling global double-hybrids (SCS/SOS-GDH), and spin-component and spin-opposite scaling range-separated double-hybrids (SCS/SOS-RSDH). Within these groups, the SCS/SOS-RSDH group of functionals is found to provide the lowest mean absolute error (MAE) values (in the range 0.020-0.148 eV) in comparison with the GDH group (0.195-0.441 eV), the RSDH group (0.186-0.511 eV), and the SCS/SOS-GDH group (0.079-0.461 eV). Of all the DH functionals and ab initio methods investigated in the present contribution, the following functionals are found to be the most accurate and consistent: SCS-ωB2GPPLYP (MAE = 0.036 eV), SOS-ωB2GPPLYP (MAE = 0.020 eV), SOS-ωB88PP86 (MAE = 0.035 eV), and SOS-ωPBEPP86 (MAE = 0.037 eV). In general, the ab initio methods tested here show mediocre performance as compared to many DH functionals.
Collapse
Affiliation(s)
- Wissam Helal
- Department of Chemistry, The University of Jordan, Amman11942, Jordan
| |
Collapse
|
9
|
Santra G, Calinsky R, Martin JML. Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights. J Phys Chem A 2022; 126:5492-5505. [PMID: 35930677 PMCID: PMC9393870 DOI: 10.1021/acs.jpca.2c03922] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
To better understand the thermochemical kinetics and
mechanism
of a specific chemical reaction, an accurate estimation of barrier
heights (forward and reverse) and reaction energies is vital. Because
of the large size of reactants and transition state structures involved
in real-life mechanistic studies (e.g., enzymatically catalyzed reactions),
density functional theory remains the workhorse for such calculations.
In this paper, we have assessed the performance of 91 density functionals
for modeling the reaction energies and barrier heights on a large
and chemically diverse data set (BH9) composed of 449 organic chemistry
reactions. We have shown that range-separated hybrid functionals perform
better than the global hybrids for BH9 barrier heights and reaction
energies. Except for the PBE-based range-separated nonempirical double
hybrids, range separation of the exchange term helps improve the performance
for barrier heights and reaction energies. The 16-parameter Berkeley
double hybrid, ωB97M(2), performs remarkably well for both properties.
However, our minimally empirical range-separated double hybrid functionals
offer marginally better accuracy than ωB97M(2) for BH9 barrier
heights and reaction energies.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Rivka Calinsky
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
10
|
Loos PF, Lipparini F, Matthews DA, Blondel A, Jacquemin D. A Mountaineering Strategy to Excited States: Revising Reference Values with EOM-CC4. J Chem Theory Comput 2022; 18:4418-4427. [PMID: 35737466 DOI: 10.1021/acs.jctc.2c00416] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the framework of the computational determination of highly accurate vertical excitation energies in small organic compounds, we explore the possibilities offered by the equation-of-motion formalism relying on the approximate fourth-order coupled-cluster (CC) method, CC4. We demonstrate, using an extended set of more than 200 reference values based on CC including up to quadruples excitations (CCSDTQ), that CC4 is an excellent approximation to CCSDTQ for excited states with a dominant contribution from single excitations with an average deviation as small as 0.003 eV. We next assess the accuracy of several additive basis set correction schemes, in which vertical excitation energies obtained with a compact basis set and a high-order CC method are corrected with lower-order CC calculations performed in a larger basis set. Such strategies are found to be overall very beneficial, though their accuracy depends significantly on the actual scheme. Finally, CC4 is employed to improve several theoretical best estimates of the QUEST database for molecules containing between four and six (nonhydrogen) atoms, for which previous estimates were computed at the CCSDT level.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Aymeric Blondel
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
11
|
Van Dijk J, Casanova-Páez M, Goerigk L. Assessing Recent Time-Dependent Double-Hybrid Density Functionals on Doublet-Doublet Excitations. ACS PHYSICAL CHEMISTRY AU 2022; 2:407-416. [PMID: 36855692 PMCID: PMC9955292 DOI: 10.1021/acsphyschemau.2c00014] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This work is the first thorough investigation of time-dependent double-hybrid density functionals (DHDFs) for the calculation of doublet-doublet excitation energies. It sheds light on the current state-of-the-art techniques in the field and clarifies if there is still room for future improvements. Overall, 29 hybrid functionals and DHDFs are investigated. We separately analyze the individual impacts of the Tamm-Dancoff approximation (TDA), range separation, and spin-component/opposite scaling (SCS/SOS) on 45 doublet-doublet excitations in 23 radicals before concluding with an overarching analysis that includes and excludes challenging excitations with double-excitation or multireference character. Our results show again that so-called "nonempirical" DHDFs are outperformed by semiempirical ones. While the best assessed functionals are DHDFs, some of the worst are also DHDFs and outperformed by all assessed hybrids. SCS/SOS is particularly beneficial for range-separated DHDFs. Spin-scaled, range-separated DHDFs paired with the TDA belong to the best tested methods here, and we particularly highlight SCS-ωB2GP-PLYP, SOS-ωB2PLYP, SOS-ωB2GP-PLYP, SOS-ωB88PP86, SOS-RSX-QIDH, and SOS-ωPBEPP86. When comparing our functional rankings with previous studies on singlet-singlet and singlet-triplet excitations, we recommend TDA-SOS-ωB88PP86 and TDA-SOS-ωPBEPP86 as robust methods for excitation energies in general until further improvements have been achieved that surpass the chemical accuracy threshold for challenging open-shell excitations without increasing the computational effort.
Collapse
Affiliation(s)
- Joshua Van Dijk
- School
of Chemistry, The University of Melbourne, Victoria 3010, Australia
| | - Marcos Casanova-Páez
- School
of Chemistry, The University of Melbourne, Victoria 3010, Australia,Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Lars Goerigk
- School
of Chemistry, The University of Melbourne, Victoria 3010, Australia,. Phone: +61 3 834 46784
| |
Collapse
|
12
|
Hancock AC, Goerigk L. Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models. RSC Adv 2022; 12:13014-13034. [PMID: 35520129 PMCID: PMC9062889 DOI: 10.1039/d2ra01703b] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Accepted: 04/12/2022] [Indexed: 01/21/2023] Open
Abstract
Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.
Collapse
Affiliation(s)
- Amy C Hancock
- School of Chemistry, The University of Melbourne Parkville Australia +61-3-8344-6784
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne Parkville Australia +61-3-8344-6784
| |
Collapse
|
13
|
Prokopiou G, Hartstein M, Govind N, Kronik L. Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals. J Chem Theory Comput 2022; 18:2331-2340. [PMID: 35369687 PMCID: PMC9009176 DOI: 10.1021/acs.jctc.2c00082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Indexed: 11/29/2022]
Abstract
We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.
Collapse
Affiliation(s)
- Georgia Prokopiou
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Michal Hartstein
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Leeor Kronik
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
14
|
Mester D, Kállay M. Charge-Transfer Excitations within Density Functional Theory: How Accurate Are the Most Recommended Approaches? J Chem Theory Comput 2022; 18:1646-1662. [PMID: 35200021 PMCID: PMC8908740 DOI: 10.1021/acs.jctc.1c01307] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The performance of
the most recent density functionals is assessed
for charge-transfer (CT) excitations using comprehensive intra- and
intermolecular CT benchmark sets with high-quality reference values.
For this comparison, the state-of-the-art range-separated (RS) and
long-range-corrected (LC) double hybrid (DH) approaches are selected,
and global DH and LC hybrid functionals are also inspected. The correct
long-range behavior of the exchange–correlation (XC) energy
is extensively studied, and various CT descriptors are compared as
well. Our results show that the most robust performance is attained
by RS-PBE-P86/SOS-ADC(2), as it is suitable to describe both types
of CT excitations with outstanding accuracy. Furthermore, concerning
the intramolecular transitions, unexpectedly excellent results are
obtained for most of the global DHs, but their limitations are also
demonstrated for bimolecular complexes. Despite the outstanding performance
of the LC-DH methods for common intramolecular excitations, serious
deficiencies are pointed out for intermolecular CT transitions, and
the wrong long-range behavior of the XC energy is revealed. The application
of LC hybrids to such transitions is not recommended in any respect.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
15
|
Alipour M, Izadkhast T. Do any types of double-hybrid models render the correct order of excited state energies in inverted singlet–triplet emitters? J Chem Phys 2022; 156:064302. [DOI: 10.1063/5.0077722] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| | - Tahereh Izadkhast
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| |
Collapse
|
16
|
Mester D, Kállay M. Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach. J Chem Theory Comput 2022; 18:865-882. [PMID: 35023739 PMCID: PMC8830052 DOI: 10.1021/acs.jctc.1c01100] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Our second-order algebraic-diagrammatic construction [ADC(2)]-based double-hybrid (DH) ansatz (J. Chem. Theory Comput. 2019, 15, 4440. DOI: 10.1021/acs.jctc.9b00391) is combined with range-separation techniques. In the present scheme, both the exchange and the correlation contributions are range-separated, while spin-scaling approaches are also applied. The new methods are thoroughly tested for the most popular benchmark sets including 250 singlet and 156 triplet excitations, as well as 80 oscillator strengths. It is demonstrated that the range separation for the correlation contributions is highly recommended for both the genuine and the ADC(2)-based DH approaches. Our results show that the latter scheme slightly but consistently outperforms the former one for single excitation dominated transitions. Furthermore, states with larger fractions of double excitations are assessed as well, and challenging charge-transfer excitations are also discussed, where the recently proposed spin-scaled long-range corrected DHs fail. The suggested iterative fourth-power scaling RS-PBE-P86/SOS-ADC(2) method, using only three adjustable parameters, provides the most robust and accurate excitation energies within the DH theory. In addition, the relative error of the oscillator strengths is reduced by 65% compared to the best genuine DH functionals.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
17
|
Alkhatib Q, Helal W, Marashdeh A. Accurate predictions of the electronic excited states of BODIPY based dye sensitizers using spin-component-scaled double-hybrid functionals: a TD-DFT benchmark study. RSC Adv 2022; 12:1704-1717. [PMID: 35425182 PMCID: PMC8978916 DOI: 10.1039/d1ra08795a] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Accepted: 01/01/2022] [Indexed: 12/31/2022] Open
Abstract
The vertical excitation energies of 13 BODIPY based dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals from different DFT rungs. Most TD-DFT results were found to overestimate the excitation energies, and show mean absolute error (MAE) values in the range 0.2–0.5 eV. The dispersion-corrected, spin-component-scaled, double-hybrid (DSD) functionals DSD-BLYP and DSD-PBEP86 were found to have the smallest MAE values of 0.083 eV and 0.106 eV, respectively, which is close to the range of average errors found in the more expensive coupled-cluster methods. Moreover, DSD-BLYP and DSD-PBEP86 functionals show excellent consistency and quality of results (standard deviation = 0.048 eV and 0.069 eV respectively). However, the range separated hybrid (RSH) and the range separated double hybrid (RSDH) functionals were found to provide the best predictability (linear determination coefficient R2 > 0.97 eV). The excitation energies of 13 BODIPY dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals. Spin-component-scaled double-hybrid (DSD) functionals are found to show the best performance.![]()
Collapse
Affiliation(s)
- Qabas Alkhatib
- Department of Chemistry, The University of Jordan Amman 11 942 Jordan
| | - Wissam Helal
- Department of Chemistry, The University of Jordan Amman 11 942 Jordan
| | - Ali Marashdeh
- Department of Chemistry, Al-Balqa Applied University 19 117 Al-Salt Jordan.,Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University P. O. Box 9502 2300 RA Leiden The Netherlands
| |
Collapse
|
18
|
Loos PF, Jacquemin D. A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Bicyclic Systems. J Phys Chem A 2021; 125:10174-10188. [PMID: 34792354 DOI: 10.1021/acs.jpca.1c08524] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Pursuing our efforts to define highly accurate estimates of the relative energies of excited states in organic molecules, we investigate, with coupled-cluster methods including iterative triples (CC3 and CCSDT), the vertical excitation energies of 10 bicyclic molecules (azulene, benzoxadiazole, benzothiadiazole, diketopyrrolopyrrole, furofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide aug-cc-pVTZ reference vertical excitation energies for 91 excited states of these relatively large systems. We use these reference values to benchmark various wave function methods, i.e., CIS(D), EOM-MP2, CC2, CCSD, STEOM-CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, ADC(2), ADC(2.5), and ADC(3), as well as some spin-scaled variants of both CC2 and ADC(2). These results are compared to those obtained previously on smaller molecules. It turns out that while the accuracy of some methods is almost unaffected by system size, e.g., CIS(D) and CC3, the performance of others can significantly deteriorate as the systems grow, e.g., EOM-MP2 and CCSD, whereas others, e.g., ADC(2) and CC2, become more accurate for larger derivatives.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, F-31062, Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
19
|
Wang J, Durbeej B. Thermal Fluctuations in Conjugation and their Effect on Calculated Excitation Energies: A Case Study on the Astaxanthin Carotenoid. CHEMPHOTOCHEM 2021. [DOI: 10.1002/cptc.202100178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Jun Wang
- Jiangsu Key Laboratory for Chemistry of Low-Dimensional Materials Jiangsu Engineering Laboratory for Environment Functional Materials School of Chemistry and Chemical Engineering Huaiyin Normal University No. 111 West Changjiang Road 223300 Huaian Jiangsu Province China
| | - Bo Durbeej
- Division of Theoretical Chemistry IFM Linköping University 581 83 Linköping Sweden
| |
Collapse
|
20
|
Casanova-Páez M, Goerigk L. Time-Dependent Long-Range-Corrected Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling: A Comprehensive Analysis of Singlet-Singlet and Singlet-Triplet Excitation Energies. J Chem Theory Comput 2021; 17:5165-5186. [PMID: 34291643 DOI: 10.1021/acs.jctc.1c00535] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Following the work on spin-component and spin-opposite scaled (SCS/SOS) global double hybrids for singlet-singlet excitations by Schwabe and Goerigk [ J. Chem. Theory Comput. 2017, 13, 4307-4323] and our own works on new long-range corrected (LC) double hybrids for singlet-singlet and singlet-triplet excitations [ J. Chem. Theory Comput. 2019, 15, 4735-4744 and J. Chem. Phys. 2020, 153, 064106], we present new LC double hybrids with SCS/SOS that demonstrate further improvement over previously published results and methods. We introduce new unscaled and scaled versions of different global and LC double hybrids based on Becke88 or PBE exchange combined with LYP, PBE, or P86 correlation. For singlet-singlet excitations, we cross-validate them on six benchmark sets that cover small to medium-sized chromophores with different excitation types (local-valence, Rydberg, and charge transfer). For singlet-triplet excitations, we perform the cross-validation on three different benchmark sets following the same analysis as in our previous work in 2020. In total, 203 excitations are analyzed. Our results confirm and extend those of Schwabe and Goerigk regarding the superior performance of SCS and SOS variants compared to their unscaled parents by decreasing mean absolute deviations, root-mean-square deviations, or error spans by more than half and bringing absolute mean deviations closer to zero. Our SCS/SOS variants are shown to be highly efficient and robust for the computation of vertical excitation energies, which even outperform specialized double hybrids that also contain an LC in their perturbative part. In particular, our new SCS/SOS-ωPBEPP86 and SCS/SOS-ωB88PP86 functionals are four of the most accurate and robust methods tested in this work, and we fully recommend them for future applications. However, if the relevant SCS and SOS algorithms are not available to the user, we suggest ωPBEPP86 as the best unscaled method in this work.
Collapse
Affiliation(s)
- Marcos Casanova-Páez
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|