1
|
Knysh I, Lipparini F, Blondel A, Duchemin I, Blase X, Loos PF, Jacquemin D. Reference CC3 Excitation Energies for Organic Chromophores: Benchmarking TD-DFT, BSE/ GW, and Wave Function Methods. J Chem Theory Comput 2024. [PMID: 39237472 DOI: 10.1021/acs.jctc.4c00906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
To expand the QUEST database of highly accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly accurate excitation energies. For several singlet transitions, we have been able to determine CCSDT transition energies with a compact basis set, finding minimal deviations from the CC3 values for most states. Subsequently, we employ these reference values to benchmark a series of lower-order wave function approaches, including the popular ADC(2) and CC2 schemes, as well as time-dependent density-functional theory (TD-DFT), both with and without applying the Tamm-Dancoff approximation (TDA). At the TD-DFT level, we evaluate a large panel of global, range-separated, local, and double hybrid functionals. Additionally, we assess the performance of the Bethe-Salpeter equation (BSE) formalism relying on both G0W0 and evGW quasiparticle energies evaluated from various starting points. It turns out that CC2 and ADC(2.5) are the most accurate models among those with respective O ( N 5 ) and O ( N 6 ) scalings with system size. In contrast, CCSD does not outperform CC2. The best performing exchange-correlation functionals include BMK, M06-2X, M06-SX, CAM-B3LYP, ωB97X-D, and LH20t, with average deviations of approximately 0.20 eV or slightly below. Errors on vertical excitation energies can be further reduced by considering double hybrids. Both SOS-ωB88PP86 and SOS-ωPBEPP86 exhibit particularly attractive performances with overall quality on par with CC2, whereas PBE0-DH and PBE-QIDH are only slightly less efficient. BSE/evGW calculations based on Kohn-Sham starting points have been found to be particularly effective for singlet transitions, but much less for their triplet counterparts.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Aymeric Blondel
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
2
|
Hancock AC, Giudici E, Goerigk L. How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models. J Comput Chem 2024; 45:1667-1681. [PMID: 38553847 DOI: 10.1002/jcc.27351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 03/07/2024] [Accepted: 03/12/2024] [Indexed: 06/04/2024]
Abstract
Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- ω PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.
Collapse
Affiliation(s)
- Amy C Hancock
- School of Chemistry, The University of Melbourne, Parkville, Victoria, Australia
| | - Erica Giudici
- School of Chemistry, The University of Melbourne, Parkville, Victoria, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria, Australia
| |
Collapse
|
3
|
Singha R, Maity P, Samanta D. Chiral Induction in a Self-Assembled Pd 4 Coordination Cage with Chiral Guests. Chemistry 2024; 30:e202401013. [PMID: 38700019 DOI: 10.1002/chem.202401013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Revised: 05/03/2024] [Accepted: 05/03/2024] [Indexed: 05/05/2024]
Abstract
The dynamic interplay of coordination bonds within metal-organic cages offers a unique avenue for structural evolution in response to external stimuli, presenting a promising strategy for the construction of chiral assemblies. This adaptability is crucial for the selective synthesis of homochiral assemblies and advancement of asymmetric catalysis. In this study, we report the self-assembly of an achiral square-planar Pd(II) acceptor with a C2-symmetric tetrapyridyl donor resulted in the formation of a racemic mixture of the chiral octahedral cage Pd4L2. The existence of this racemic mixture was confirmed using circular dichroism spectroscopy as well as single crystal X-ray diffraction analysis. We encoded chiral information into the asymmetric cavity of the cage by encapsulating chiral aromatic guests through efficient π-π stacking and hydrophobic interactions in aqueous media. The inclusion of a chiral guest induces a preference for one enantiomeric conformation of the cage over the other, effectively shifting the equilibrium towards a single, enantiopure host-guest complex. While the concept of chiral guest recognition by a chiral host is well-established, this work constitutes a remarkable example of guest-mediated chirality transfer leading to the formation of a single enantiopure coordination complex from achiral building blocks.
Collapse
Affiliation(s)
- Raghunath Singha
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), An OCC of Homi Bhabha National Institute, Bhubaneswar, Jatni, Khurda, Odisha, 752050, India
| | - Pankaj Maity
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), An OCC of Homi Bhabha National Institute, Bhubaneswar, Jatni, Khurda, Odisha, 752050, India
| | - Dipak Samanta
- School of Chemical Sciences, National Institute of Science Education and Research (NISER), An OCC of Homi Bhabha National Institute, Bhubaneswar, Jatni, Khurda, Odisha, 752050, India
| |
Collapse
|
4
|
Montserrat R, Oliveira RR, Rocha AB. Total absorption spectrum of benzene aggregates obtained from two different approaches. J Mol Model 2024; 30:66. [PMID: 38345753 DOI: 10.1007/s00894-024-05859-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Accepted: 01/29/2024] [Indexed: 03/16/2024]
Abstract
CONTEXT The study of molecular aggregation effects on the electronic spectrum is essential for the development of optoelectronic devices. However, investigating the entire valence absorption spectrum of aggregates using quantum mechanical methods is a challenging task. In this work, we perform systematic simulations of the absorption spectrum of benzene molecular clusters up to 35 eV applying two approaches based on time-dependent density functional theory. The results show that depending on the dimer packing, different energy shifts occur for the symmetry allowed [Formula: see text] transition, in comparison to the monomer. The transition intensity increases for the band around 6 eV for larger aggregates from the monomer to dimers and tetramer, indicating the occurrence of the symmetry forbidden (in [Formula: see text] point group) [Formula: see text] [Formula: see text] transition. The benzene crystal exhibits a large redshift following the experimental spectrum. Also, the continuum regions of all spectra show a good agreement with the experiments both in gas and solid phases. METHODS Geometry optimization of the monomer was carried out with Gaussian 09 software using the PBE0/def2-TZVP level of theory. We used dimers and tetramer molecular geometries extracted from the experimental crystal structure. The absorption spectra were directly obtained by the Liouville-Lanczos TDDFT approach with plane waves basis set or indirectly by TDDFT pseudo-spectra calculated in a [Formula: see text] basis followed by analytic continuation procedure to obtain complex polarizability. The former is available at Quantum ESPRESSO, and the latter was calculated using Gaussian 09 with the post-processing performed with a code previously developed in our group.
Collapse
Affiliation(s)
- Ricardo Montserrat
- Departamento de Físico-Química, Universidade Federal do Rio de Janeiro, Avenida Athos da Silveira Ramos, 149, Rio de Janeiro, 21941-909, Rio de Janeiro, Brazil
| | - Ricardo R Oliveira
- Departamento de Físico-Química, Universidade Federal do Rio de Janeiro, Avenida Athos da Silveira Ramos, 149, Rio de Janeiro, 21941-909, Rio de Janeiro, Brazil.
| | - Alexandre B Rocha
- Departamento de Físico-Química, Universidade Federal do Rio de Janeiro, Avenida Athos da Silveira Ramos, 149, Rio de Janeiro, 21941-909, Rio de Janeiro, Brazil.
| |
Collapse
|
5
|
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
|
6
|
Sitkiewicz SP, Matito E, Luis JM, Zaleśny R. Pitfall in simulations of vibronic TD-DFT spectra: diagnosis and assessment. Phys Chem Chem Phys 2023; 25:30193-30197. [PMID: 37905423 DOI: 10.1039/d3cp04276f] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
In this Communication, we study the effect of spurious oscillations in the profiles of energy derivatives with respect to nuclear coordinates calculated with density functional approximations (DFAs) for formaldehyde, pyridine, and furan in their ground and electronic excited states. These spurious oscillations, which can only be removed using extensive integration grids that increase enormously the CPU cost of DFA calculations, are significant in the case of third- and fourth-order energy derivatives of the ground and excited states computed by M06-2X and ωB97X functionals. The errors in question propagate to anharmonic vibronic spectra computed under the Franck-Condon approximation, i.e., positions and intensities of vibronic transitions are affected to a large extent (shifts as significant as hundreds of cm-1 were observed). On the other hand, the LC-BLYP and CAM-B3LYP functionals show a much less pronounced effect due to spurious oscillations. Based on the results presented herein, we recommend either LC-BLYP or CAM-B3LYP with integration grids (250, 974) (or larger) for numerically stable simulations of vibronic spectra including anharmonic effects.
Collapse
Affiliation(s)
- Sebastian P Sitkiewicz
- Wrocław Centre for Networking and Supercomputing, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, Wrocław PL-50370, Poland.
| | - Eduard Matito
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, Donostia 20018, Euskadi, Spain
- Ikerbasque Foundation for Science, Bilbao 48011, Euskadi, Spain
| | - Josep M Luis
- Institute of Computational Chemistry and Catalysis and Department of Chemistry, University of Girona, Campus de Montilivi, 17003, Girona, Catalonia, Spain.
| | - Robert Zaleśny
- Faculty of Chemistry, Wroclaw University of Science and Technology, Wybrzeże Wyspiańskiego 27, Wrocław 50-370, Poland.
| |
Collapse
|
7
|
Aschi M, Palombi L, Amadei A. Theoretical-Computational Modeling of CD Spectra of Aqueous Monosaccharides by Means of Molecular Dynamics Simulations and Perturbed Matrix Method. Molecules 2023; 28:molecules28083591. [PMID: 37110825 PMCID: PMC10144652 DOI: 10.3390/molecules28083591] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2023] [Revised: 04/12/2023] [Accepted: 04/13/2023] [Indexed: 04/29/2023] Open
Abstract
The electronic circular dichroism (ECD) spectra of aqueous d-glucose and d-galactose were modeled using a theoretical-computational approach combining molecular dynamics (MD) simulations and perturbed matrix method (PMM) calculations, hereafter termed MD-PMM. The experimental spectra were reproduced with a satisfactory accuracy, confirming the good performances of MD-PMM in modeling different spectral features in complex atomic-molecular systems, as already reported in previous studies. The underlying strategy of the method was to perform a preliminary long timescale MD simulation of the chromophore followed by the extraction of the relevant conformations through essential dynamics analysis. On this (limited) number of relevant conformations, the ECD spectrum was calculated via the PMM approach. This study showed that MD-PMM was able to reproduce the essential features of the ECD spectrum (i.e., the position, the intensity, and the shape of the bands) of d-glucose and d-galactose while avoiding the rather computationally expensive aspects, which were demonstrated to be important for the final outcome, such as (i) the use of a large number of chromophore conformations; (ii) the inclusion of quantum vibronic coupling; and (iii) the inclusion of explicit solvent molecules interacting with the chromophore atoms within the chromophore itself (e.g., via hydrogen bonds).
Collapse
Affiliation(s)
- Massimiliano Aschi
- Dipartimento di Scienze Fisiche e Chimiche, Università dell'Aquila via Vetoio (Coppito 2), 67010 l'Aquila, Italy
| | - Laura Palombi
- Dipartimento di Scienze Fisiche e Chimiche, Università dell'Aquila via Vetoio (Coppito 2), 67010 l'Aquila, Italy
| | - Andrea Amadei
- Dipartimento di Scienze e Tecnologie Chimiche, Università di Roma "Tor Vergata", Via della Ricerca Scientifica 1, 00133 Roma, Italy
| |
Collapse
|
8
|
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
|
9
|
Kempfer-Robertson EM, Haase MN, Bersson JS, Avdic I, Thompson LM. Role of Exact Exchange in Difference Projected Double-Hybrid Density Functional Theory for Treatment of Local, Charge Transfer, and Rydberg Excitations. J Phys Chem A 2022; 126:8058-8069. [PMID: 36269072 DOI: 10.1021/acs.jpca.2c04338] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Difference approaches to the study of excited states have undergone a renaissance in recent years, with the development of a plethora of algorithms for locating self-consistent field approximations to excited states. Density functional theory is likely to offer the best balance of cost and accuracy for difference approaches, and yet there has been little investigation of how the parametrization of density functional approximations affects performance. In this work, we aim to explore the role of the global Hartree-Fock exchange parameter in tuning accuracy of different excitation types within the framework of the recently introduced difference projected double-hybrid density functional theory approach and contrast the performance with conventional time-dependent double-hybrid density functional theory. Difference projected double-hybrid density functional theory was demonstrated to give vertical excitation energies with average error and standard deviation with respect to multireference perturbation theory comparable to more expensive linear-response coupled cluster approaches ( J. Chem. Phys.2020, 153, 074103). However, despite benchmarking of local excitations, there has been no investigation of the methods performance for charge transfer or Rydberg excitations. In this work we report a new benchmark of charge transfer, Rydberg, and local excited state vertical excitation energies and examine how the exact Hartree-Fock exchange affects the benchmark performance to provide a deeper understanding of how projection and nonlocal correlation balance differing sources of error in the ground and excited states.
Collapse
Affiliation(s)
| | - Meagan N Haase
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Jonathan S Bersson
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Irma Avdic
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Lee M Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| |
Collapse
|
10
|
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
|
11
|
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
|
12
|
Santra G, Martin JML. Do Double-Hybrid Functionals Benefit from Regularization in the PT2 Term? Observations from an Extensive Benchmark. J Phys Chem Lett 2022; 13:3499-3506. [PMID: 35417181 PMCID: PMC9036584 DOI: 10.1021/acs.jpclett.2c00718] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We put to the test a recent suggestion [Shee, J., et al. J. Phys. Chem. Lett. 2021, 12 (50), 12084-12097] that MP2 regularization might improve the performance of double-hybrid density functionals. Using the very large and chemically diverse GMTKN55 benchmark, we find that κ-regularization is indeed beneficial at lower percentages of Hartree-Fock exchange, especially if spin-component scaling is not applied [such as in B2GP-PLYP or ωB97M(2)]. This benefit dwindles for DSD and DOD functionals and vanishes entirely in the ∼70% HF exchange region optimal for them.
Collapse
|
13
|
Niemeyer N, Caricato M, Neugebauer J. Origin invariant electronic circular dichroism in the length dipole gauge without London atomic orbitals. J Chem Phys 2022; 156:154114. [PMID: 35459317 DOI: 10.1063/5.0088922] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a method for obtaining origin-independent electronic circular dichroism (ECD) in the length-gauge representation LG(OI) without the usage of London atomic orbitals. This approach builds upon the work by Caricato [J. Chem. Phys. 153, 151101 (2020)] and is applied to rotatory strengths and ECD spectra from damped response theory. Numerical results are presented for time-dependent Hartree-Fock and density-functional theory, the second-order algebraic diagrammatic construction method, and linear-response coupled-cluster theory with singles and approximate doubles. We can support the finding that the common choice of placing the gauge origin in the center of mass of a molecule in conventional length-gauge calculations involving chiroptical properties might not be optimal and show that LG(OI) is a valuable alternative for the origin-independent calculation of ECD spectra. We show that, for a limited test set, the convergence of the rotatory strengths calculated with the LG(OI) approach toward the basis-set limit tends to be faster than for the established velocity gauge representation. Relationships between the sum-over-states expression of the optical rotation in the LG(OI) framework and its representation in terms of response functions are analyzed.
Collapse
Affiliation(s)
- Niklas Niemeyer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Marco Caricato
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, USA
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
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] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [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
|
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] [Received: 11/02/2021] [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
|
16
|
Santra G, Semidalas E, Martin JML. Surprisingly Good Performance of XYG3 Family Functionals Using a Scaled KS-MP3 Correlation. J Phys Chem Lett 2021; 12:9368-9376. [PMID: 34550706 DOI: 10.1021/acs.jpclett.1c02838] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
By adding a GLPT3 (third-order Görling-Levy perturbation theory, or KS-MP3) term E3 to the XYG7 form for a double hybrid, we are able to bring down WTMAD2 (weighted total mean absolute deviation) for the very large and chemically diverse GMTKN55 benchmark to an unprecedented 1.17 kcal/mol, competitive with much costlier composite wave function ab initio approaches. Intriguingly, (a) the introduction of E3 makes an empirical dispersion correction redundant; (b) generalized gradient approximation (GGA) or meta-GGA semilocal correlation functionals offer no advantage over the local density approximation (LDA) in this framework; (c) if a dispersion correction is retained, then simple Slater exchange leads to no significant loss in accuracy. It is possible to create a six-parameter functional with WTMAD2 = 1.42 that has no post-LDA density functional theory components and no dispersion correction in the final energy.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Emmanouil Semidalas
- 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
|
17
|
Meskers SCJ. Circular Polarization of Luminescence as a Tool To Study Molecular Dynamical Processes. CHEMPHOTOCHEM 2021. [DOI: 10.1002/cptc.202100154] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Stefan C. J. Meskers
- Molecular Materials and Nanosystems and Institute for Complex Molecular Systems Eindhoven University of Technology P.O. box 513 (STW 4.37) NL 5600 MB Eindhoven Netherlands
| |
Collapse
|
18
|
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
|
19
|
Mester D, Kállay M. Spin-Scaled Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:4211-4224. [PMID: 34152771 PMCID: PMC8280718 DOI: 10.1021/acs.jctc.1c00422] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Indexed: 11/28/2022]
Abstract
Our recently presented range-separated (RS) double-hybrid (DH) time-dependent density functional approach [J. Chem. Theory Comput. 17, 927 (2021)] is combined with spin-scaling techniques. The proposed spin-component-scaled (SCS) and scaled-opposite-spin (SOS) variants are thoroughly tested for almost 500 excitations including the most challenging types. This comprehensive study provides useful information not only about the new approaches but also about the most prominent methods in the DH class. The benchmark calculations confirm the robustness of the RS-DH ansatz, while several tendencies and deficiencies are pointed out for the existing functionals. Our results show that the SCS variant consistently improves the results, while the SOS variant preserves the benefits of the original RS-DH method reducing its computational expenses. It is also demonstrated that, besides our approaches, only the nonempirical functionals provide balanced performance for general applications, while particular methods are only suggested for certain types of excitations.
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
|
20
|
Mester D, Kállay M. A Simple Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:927-942. [PMID: 33400872 PMCID: PMC7884002 DOI: 10.1021/acs.jctc.0c01135] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Indexed: 01/12/2023]
Abstract
A simple and robust range-separated (RS) double-hybrid (DH) time-dependent density functional approach is presented for the accurate calculation of excitation energies of molecules within the Tamm-Dancoff approximation. The scheme can be considered as an excited-state extension of the ansatz proposed by Toulouse and co-workers [J. Chem. Phys. 2018, 148, 164105], which is based on the two-parameter decomposition of the Coulomb potential, for which both the exchange and correlation contributions are range-separated. A flexible and efficient implementation of the new scheme is also presented, which facilitates its extension to any combination of exchange and correlation functionals. The performance of the new approximation is tested for singlet excitations on several benchmark compilations and thoroughly compared to that of representative DH, RS hybrid, and RS DH functionals. The one-electron basis set dependence and computation times are also assessed. Our results show that the new approach improves on standard DHs in most cases, and it can provide a more robust and accurate alternative. In addition, on average, it noticeably surpasses the existing RS hybrid and RS 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, Budapest, H-1521, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| |
Collapse
|
21
|
Casanova-Páez M, Goerigk L. Global double hybrids do not work for charge transfer: A comment on "Double hybrids and time-dependent density functional theory: An implementation and benchmark on charge transfer excited states". J Comput Chem 2021; 42:528-533. [PMID: 33415788 DOI: 10.1002/jcc.26478] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Revised: 11/17/2020] [Accepted: 12/18/2020] [Indexed: 01/13/2023]
Abstract
We comment on the results published by Ottochian et al. in J. Comput Chem. 2020, 41, 1242. Therein, the authors claim that the second-order, perturbative correlation correction applied to the time-dependent version of the PBE-QIDH global double-hybrid functional approximation (DHDFA) enables the description of charge-transfer (CT) excitations. Herein, we point out some inadvertent oversights related to what had already been previously known and achieved in the field of time-dependent DHDFAs. Exemplified for the same four systems that Ottochian et al. have used to analyze intermolecular CT excitations, we show how a systematic and unacceptably large redshift in global DHDFAs is rectified when using the latest long-range corrected DHDFAs published earlier in J. Chem. Theory Comput. 2019, 15, 4735.
Collapse
Affiliation(s)
- Marcos Casanova-Páez
- School of Chemistry, The University of Melbourne, Melbourne, Victoria, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Melbourne, Victoria, Australia
| |
Collapse
|
22
|
Goerigk L, Casanova-Paéz M. The Trip to the Density Functional Theory Zoo Continues: Making a Case for Time-Dependent Double Hybrids for Excited-State Problems. Aust J Chem 2021. [DOI: 10.1071/ch20093] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This account is written for general users of time-dependent density functional theory (TD-DFT) methods as well as chemists who are unfamiliar with the field. It includes a brief overview of conventional TD-DFT approaches and recommendations for applications to organic molecules based on our own experience. The main emphasis of this work, however, lies in providing the first in-depth review of time-dependent double-hybrid density functionals. They were first established in 2007 with very promising follow-up studies in the subsequent four years before developments or applications became scarce. The topic has regained more interest since 2017, and this account reviews those latest developments led by our group. These developments have shown unprecedented robustness for a variety of different types of electronic excitations when compared to more conventional TD-DFT methods. In particular, time-dependent double hybrids do not suffer from artificial ghost states and are able to reproduce exciton-coupled absorption spectra. Our latest methods include range separation and belong to the currently best TD-DFT methods for singlet-singlet excitations in organic molecules. While there is still room for improvement and further development in this space, we hope that this account encourages users to adjust their computational protocols to such new methods to provide more real-life testing and scenarios.
Collapse
|
23
|
Giovannini T, Egidi F, Cappelli C. Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems. Phys Chem Chem Phys 2020; 22:22864-22879. [PMID: 33043930 DOI: 10.1039/d0cp04027d] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | | | | |
Collapse
|
24
|
Casanova-Páez M, Goerigk L. Assessing the Tamm–Dancoff approximation, singlet–singlet, and singlet–triplet excitations with the latest long-range corrected double-hybrid density functionals. J Chem Phys 2020; 153:064106. [DOI: 10.1063/5.0018354] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
25
|
Xu Q, Liu Y, Zhao X, Chen S, Li Q, Wang M, Yang C. Vibrationally resolved electronic circular dichroism and circularly polarized luminescence spectra of a boron-fused double helicene: A theoretical study. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020; 231:118132. [PMID: 32058915 DOI: 10.1016/j.saa.2020.118132] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Revised: 02/03/2020] [Accepted: 02/03/2020] [Indexed: 06/10/2023]
Abstract
In this work, we present the theoretical study of the vibrationally resolved absorption (ABS), electronic circular dichroism (ECD), emission (EMI), and circularly polarized luminescence (CPL) spectra of a boron-fused double helicene, with a detailed and complete discussion of the alternative possible approximate methods. Given the fact that few examples of CPL calculations exist, the potential energy surfaces (PESs) have been constructed and compared with Adiabatic (AH) and Vertical Hessian (VH) models. All the vibronic calculations have accounted for Duschinsky mixings, Franck-Condon (FC) effect and Herzberg-Teller (HT) contribution. Moreover, different HT expansions have been checked and compared, by computing the derivatives of the electric and magnetic dipole transition moments around the equilibrium geometries of the initial and final states. Our results show that both AH and VH models have well reproduced the experimental vibronic structures and VH model shows a better performance in the simulation of spectral lineshapes. They also show that HT effects dominate the shapes of EMI and CPL, tuning the relative heights of the different vibronic peaks, improving the agreement with the experiment for EMI. Moreover, HT effects are the main reason for the mirror-symmetry breaking between ECD and CPL spectra. Furthermore, interesting interference effects between FC and HT contributions have also been addressed.
Collapse
Affiliation(s)
- Qiushuang Xu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China; School of Physics Engineering, Qufu Normal University, 2673100 Qufu, Shandong, PR China
| | - Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China.
| | - Xian Zhao
- State Key Laboratory of Crystal Materials, Shandong University, 250100 Jinan, Shandong, PR China
| | - Shenghui Chen
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Quanjing Li
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Meishan Wang
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China.
| | - Chuanlu Yang
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| |
Collapse
|
26
|
Lonsdale DR, Goerigk L. The one-electron self-interaction error in 74 density functional approximations: a case study on hydrogenic mono- and dinuclear systems. Phys Chem Chem Phys 2020; 22:15805-15830. [PMID: 32458849 DOI: 10.1039/d0cp01275k] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The self-interaction error (SIE), i.e. unphysical interactions of electrons with themselves, has plagued developers and users of Density Functional Approximations (DFAs) since the inception of Density Functional Theory (DFT). Formally, it can be separated into the one-electron and many-electron SIE; herein we present one of the most comprehensive studies of the first. While we focus mostly on the total SIE, we also make use of two different decompositions. The first is a separation into functional and density-driven errors as championed by Sim, Burke and co-workers [J. Phys. Chem. Lett., 2018, 9, 6385-6392]; the second separates the error into exchange, correlation, and one-electron components, with the latter being a density error that has not been discussed in this form before. After investigating the familiar hydrogen atom and dihydrogen cation, we establish a relationship between the SIE and the nuclear charge with the help of a series of heavier hydrogenic analogues. For the mononuclear systems and the diatomics at the dissociation limit, this relationship is linear in nature with prominent exceptions, mostly belonging to the Minnesota and range-separated (double-)hybrid DFAs. For the first time, we also show how the magnitude of the SIE depends on the underlying atomic-orbital basis set and how DFAs that rely on a popular van-der-Waals DFT type London-dispersion term exhibit "self-dispersion". We find that range separation is not a panacea for solving the one-electron SIE. DFAs that have been developed to be one-electron SIE free for one system, such as the hydrogen atom, show larger errors for heavier hydrogenic systems. Often, one-electron SIE-free DFAs rely on fortuitous error cancellation between their exchange and correlation components. An analysis of the most robust methods for general applications to date reveals that they suffer moderately from the one-electron SIE, while DFAs that are nearly SIE-free do not perform well in applications. Implicit in the continued existence of the one-electron SIE is that well-performing DFAs continue to suffer insufficiencies at their fundamental levels that are being compensated for by the SIE. Our analysis includes more than 250 000 datapoints, resulting in multiple insights that may drive future developments of new DFAs or SIE corrections.
Collapse
Affiliation(s)
- Dale R Lonsdale
- School of Chemistry, The University of Melbourne, Parkville, Australia.
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Australia.
| |
Collapse
|
27
|
Salim F, Yunus YM, Anouar EH, Awang K, Langat M, Cordell GA, Ahmad R. Absolute Configuration of Alkaloids from Uncaria longiflora through Experimental and Computational Approaches. JOURNAL OF NATURAL PRODUCTS 2019; 82:2933-2940. [PMID: 31686505 DOI: 10.1021/acs.jnatprod.8b00380] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The structure elucidation of three new alkaloids named isoformosaninol (1), formosaninol (2), and longiflorine (3), isolated from the leaves of Uncaria longiflora var. pteropoda (Miq.) Ridsdale, along with their biosynthetic pathways are discussed. Their absolute structures were determined through a combination of physical data interpretation and quantum chemical calculations using the time-dependent density functional theory (TDDFT) method.
Collapse
Affiliation(s)
- Fatimah Salim
- Atta-ur-Rahman Institute for Natural Product Discovery (AuRIns) , Universiti Teknologi MARA Selangor Branch , Puncak Alam Campus , 42300 Bandar Puncak Alam , Selangor , Malaysia
| | - Yusri Mohd Yunus
- Faculty of Applied Sciences , Universiti Teknologi MARA , 40450 Shah Alam , Selangor , Malaysia
| | - El Hassane Anouar
- Department of Chemistry, College of Sciences and Humanites , Prince Sattam bin Abdulaziz University , P.O. Box 83, 11942 , Al Kharj , Saudi Arabia
| | - Khalijah Awang
- Department of Chemistry, Faculty of Science , University Malaya , 50603 Kuala Lumpur , Malaysia
| | - Moses Langat
- Natural Product Research Group, Department of Chemistry, Faculty of Engineering and Physical Sciences , University of Surrey , Guildford GU2 7XH , United Kingdom
- Jodrell Laboratory, Natural Capital and Plant Health Department, Royal Botanic Gardens , Kew, Richmond, Surrey TW9 3DS , United Kingdom
| | - Geoffrey A Cordell
- Natural Products Inc. , Evanston , Illinois 60203 , United States
- College of Pharmacy , University of Florida , Gainesville , Florida 32610 , United States
| | - Rohaya Ahmad
- Faculty of Applied Sciences , Universiti Teknologi MARA , 40450 Shah Alam , Selangor , Malaysia
| |
Collapse
|
28
|
Casanova-Páez M, Dardis MB, Goerigk L. ωB2PLYP and ωB2GPPLYP: The First Two Double-Hybrid Density Functionals with Long-Range Correction Optimized for Excitation Energies. J Chem Theory Comput 2019; 15:4735-4744. [PMID: 31298850 DOI: 10.1021/acs.jctc.9b00013] [Citation(s) in RCA: 86] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Double-hybrid density functionals are currently the most accurate density functionals for ground-state properties and electronic excitations. Nevertheless, the lack of a long-range correction scheme makes them unreliable when it comes to long-range excitations. For this reason, we propose the first two time-dependent double-hybrid functionals with correct asymptotic long-range behavior named ωB2PLYP and ωB2GPPLYP. Herein, we demonstrate their excellent performance and show that they are the most accurate and robust time-dependent density functional theory methods for electronic excitation energies. They provide a balanced description of local-valence, Rydberg, and charge-transfer states. They are also able to tackle the difficult first two transitions in polycyclic aromatic hydrocarbons and show very promising results in a preliminary study on transition-metal compounds, exemplified for titanium dioxide clusters. This work shows that double hybrids can be systematically improved also for excitation energies, and further work in this field is warranted.
Collapse
Affiliation(s)
| | - Michael B Dardis
- School of Chemistry , The University of Melbourne , Victoria 3010 , Australia
| | - Lars Goerigk
- School of Chemistry , The University of Melbourne , Victoria 3010 , Australia
| |
Collapse
|
29
|
Wappett DA, Goerigk L. Toward a Quantum-Chemical Benchmark Set for Enzymatically Catalyzed Reactions: Important Steps and Insights. J Phys Chem A 2019; 123:7057-7074. [DOI: 10.1021/acs.jpca.9b05088] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
30
|
Mester D, Kállay M. Combined Density Functional and Algebraic-Diagrammatic Construction Approach for Accurate Excitation Energies and Transition Moments. J Chem Theory Comput 2019; 15:4440-4453. [PMID: 31265275 DOI: 10.1021/acs.jctc.9b00391] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A composite of time-dependent density functional theory (TDDFT) and the second-order algebraic-diagrammatic construction [ADC(2)] approach is presented for efficient calculation of spectral properties of molecules. Our method can be regarded as a new excited-state double-hybrid (DH) approach or a dressed TDDFT scheme, but it can also be interpreted as an empirically tuned ADC(2) model. Several combinations of exchange-correlation functionals and spin-scaling schemes are explored. Our best-performing method includes the Perdew, Burke, and Ernzerhof exchange and Perdew's 1986 correlation functional and employs the scaled-opposite-spin approximation for the higher-order terms. The computation time of the new method scales as the fourth power of the system size, and an efficient cost-reduction approach is also presented, which further speeds up the calculations. Our benchmark calculations show that the proposed model outperforms not only the existing DH approaches and ADC(2) variants but also the considerably more expensive coupled-cluster methods.
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
|
31
|
Michaelis M, Hildebrand N, Meißner RH, Wurzler N, Li Z, Hirst JD, Micsonai A, Kardos J, Delle Piane M, Colombi Ciacchi L. Impact of the Conformational Variability of Oligopeptides on the Computational Prediction of Their CD Spectra. J Phys Chem B 2019; 123:6694-6704. [DOI: 10.1021/acs.jpcb.9b03932] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- M. Michaelis
- Faculty of Production Engineering, Bremen Center for Computational Materials Science, Center for Environmental Research and Sustainable Technology (UFT), and MAPEX Center for Materials and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, Bremen 28359, Germany
- Biomolecular and Materials Interface Research Group, Interdisciplinary Biomedical Research Centre, School of Science and Technology, Nottingham Trent University, Clifton Lane, Nottingham NG11 8NS, United Kingdom
| | - N. Hildebrand
- Faculty of Production Engineering, Bremen Center for Computational Materials Science, Center for Environmental Research and Sustainable Technology (UFT), and MAPEX Center for Materials and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, Bremen 28359, Germany
| | - R. H. Meißner
- Faculty of Production Engineering, Bremen Center for Computational Materials Science, Center for Environmental Research and Sustainable Technology (UFT), and MAPEX Center for Materials and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, Bremen 28359, Germany
| | - N. Wurzler
- Faculty of Production Engineering, Bremen Center for Computational Materials Science, Center for Environmental Research and Sustainable Technology (UFT), and MAPEX Center for Materials and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, Bremen 28359, Germany
| | - Z. Li
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - J. D. Hirst
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - A. Micsonai
- Department of Biochemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/C, Budapest H-1117, Hungary
| | - J. Kardos
- Department of Biochemistry, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/C, Budapest H-1117, Hungary
| | - M. Delle Piane
- Faculty of Production Engineering, Bremen Center for Computational Materials Science, Center for Environmental Research and Sustainable Technology (UFT), and MAPEX Center for Materials and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, Bremen 28359, Germany
| | - L. Colombi Ciacchi
- Faculty of Production Engineering, Bremen Center for Computational Materials Science, Center for Environmental Research and Sustainable Technology (UFT), and MAPEX Center for Materials and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, Bremen 28359, Germany
| |
Collapse
|
32
|
Bedi A, Gidron O. Chiroptical Properties of Twisted Acenes: Experimental and Computational Study. Chemistry 2019; 25:3279-3285. [PMID: 30614087 DOI: 10.1002/chem.201805728] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Indexed: 01/29/2023]
Abstract
Acenes that are twisted out of planarity are expected to display chiroptical properties. However, the effect of twisting on the chiroptical properties of acenes has not been investigated computationally or experimentally. Herein, we present a computational investigation of the chiroptical properties of anthracenes to pentacenes, combined with an experimental study using a series of helically locked acenes, twisted to different torsional angles in their enantiopure form. The lowest energy transition, which is relatively weak in acenes, becomes dominant in their circular dichroism spectra upon twisting. We find that the rotational strength of acenes consistently increases with increasing twist. The experimental data obtained from enantiopure tethered twistacenes show the same trend as the calculated result, with a strong Cotton effect and anisotropy factor, rendering twisted acenes as excellent chiroptical materials.
Collapse
Affiliation(s)
- Anjan Bedi
- Institute of Chemistry, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Jerusalem, Israel
| | - Ori Gidron
- Institute of Chemistry, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Jerusalem, Israel
| |
Collapse
|
33
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 128] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
34
|
Miyahara T, Nakatsuji H. Accuracy of Td-DFT in the Ultraviolet and Circular Dichroism Spectra of Deoxyguanosine and Uridine. J Phys Chem A 2018; 122:100-118. [PMID: 29190101 DOI: 10.1021/acs.jpca.7b09733] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accuracy of the time-dependent density functional theory (Td-DFT) was examined for the ultraviolet (UV) and circular dichroism (CD) spectra of deoxyguanosine (dG) and uridine, using 11 different DFT functionals and two different basis sets. The Td-DFT results of the UV and CD spectra were strongly dependent on the functionals used. The basis-set dependence was observed only for the CD spectral calculations. For the UV spectra, the B3LYP and PBE0 functionals gave relatively good results. For the CD spectra, the B3LYP and PBE0 with 6-311G(d,p) basis gave relatively permissible result only for dG. The results of other functionals were difficult to be used for the studies of the UV and CD spectra, though the symmetry adapted cluster-configuration interaction (SAC-CI) method reproduced well the experimental spectra of these molecules. To obtain valuable information from the theoretical calculations of the UV and CD spectra, the theoretical tool must be able to reproduce correctly both of the intensities and peak positions of the UV and CD spectra. Then, we can analyze the reasons of the changes of the intensity and/or the peak position to clarify the chemistry involved. It is difficult to recommend Td-DFT as such tools of science, at least from the examinations using dG and uridine.
Collapse
Affiliation(s)
- Tomoo Miyahara
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16 , 14 Yoshida Kawara-machi, Sakyou-ku, Kyoto 606-8305, Japan
| | - Hiroshi Nakatsuji
- Quantum Chemistry Research Institute, Kyoto Technoscience Center 16 , 14 Yoshida Kawara-machi, Sakyou-ku, Kyoto 606-8305, Japan
| |
Collapse
|
35
|
Azarias C, Pascal S, Siri O, Jacquemin D. Central substitution of azacalixphyrins: a strategy towards acidochromic NIR dyes. Phys Chem Chem Phys 2018; 20:20056-20069. [DOI: 10.1039/c8cp03789b] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
A new substitution pattern of azacalixphyrins is proposed and assessed with theoretical tools.
Collapse
Affiliation(s)
- Cloé Azarias
- Chimie Et Interdisciplinarité, Synthèse, Analyse, Modélisation (CEISAM)
- UMR CNRS no. 6230
- BP 92208, Université de Nantes, 2, Rue de la Houssinière
- 44322 Nantes
- France
| | - Simon Pascal
- Aix Marseille Univ., CNRS-CINAM
- Campus de Luminy
- Case 913
- 13288 Marseille Cedex 09
- France
| | - Olivier Siri
- Aix Marseille Univ., CNRS-CINAM
- Campus de Luminy
- Case 913
- 13288 Marseille Cedex 09
- France
| | - Denis Jacquemin
- Chimie Et Interdisciplinarité, Synthèse, Analyse, Modélisation (CEISAM)
- UMR CNRS no. 6230
- BP 92208, Université de Nantes, 2, Rue de la Houssinière
- 44322 Nantes
- France
| |
Collapse
|
36
|
Brémond É, Savarese M, Pérez-Jiménez ÁJ, Sancho-García JC, Adamo C. Speed-Up of the Excited-State Benchmarking: Double-Hybrid Density Functionals as Test Cases. J Chem Theory Comput 2017; 13:5539-5551. [DOI: 10.1021/acs.jctc.7b00627] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Éric Brémond
- CompuNet, Istituto
Italiano di Tecnologia, via Morego
30, I-16163 Genoa, Italy
| | - Marika Savarese
- CompuNet, Istituto
Italiano di Tecnologia, via Morego
30, I-16163 Genoa, Italy
| | | | | | - Carlo Adamo
- Chimie
ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris IRCP, F-75005 Paris, France
| |
Collapse
|
37
|
Schwabe T, Goerigk L. Time-Dependent Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling. J Chem Theory Comput 2017; 13:4307-4323. [PMID: 28763220 DOI: 10.1021/acs.jctc.7b00386] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
For the first time, we combine time-dependent double-hybrid density functional approximations (TD-DHDFAs) for the calculation of electronic excitation energies with the concepts of spin-component and spin-opposite scaling (SCS/SOS) of electron-pair contributions to their nonlocal correlation components. Different flavors of this idea, ranging from standard SCS parameters to fully fitted parameter sets, are presented and tested on six different parent DHDFAs. For cross-validation, we assess those methods on three benchmark sets that cover small- to medium-sized chromophores (up to 78 atoms) and different excitation types. For this purpose, we also introduce new CC3 reference values for the popular Gordon benchmark set that we recommend using in future studies. Our results confirm that already the (unscaled) parent TD-DHDFAs are accurate and outperform some wave function methods. Further introduction of SCS/SOS eliminates extreme outliers, reduces deviation spans from reference values by up to 0.5 eV, aligns the performance of the Tamm-Dancoff approximation (TDA) to that of full TD calculations, and also enables a more balanced description of different excitation types. The best-performing TD-based methods in our cross validation have mean absolute deviations as low as 0.14 eV compared to the time- and resource-intensive CC3 approach. A very important finding is that we also obtained SOS variants with excellent performance, contrary to wave function based methods. This opens a future pathway to highly efficient methods for the optimization of excited-state geometries, particularly when paired with computing strategies such as the Laplace transform. We recommend our SCS- and SOS-based variants for further testing and subsequent applications.
Collapse
Affiliation(s)
- Tobias Schwabe
- Center for Bioinformatics and Physical Chemistry Institute, University of Hamburg , Bundesstraße 43, D-20146 Hamburg, Germany
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne , Parkville, Victoria 3010, Australia
| |
Collapse
|
38
|
Alipour M, Karimi N. Dissecting the accountability of parameterized and parameter-free single-hybrid and double-hybrid functionals for photophysical properties of TADF-based OLEDs. J Chem Phys 2017. [PMID: 28641443 DOI: 10.1063/1.4986777] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
Organic light emitting diodes (OLEDs) based on thermally activated delayed fluorescence (TADF) emitters are an attractive category of materials that have witnessed a booming development in recent years. In the present contribution, we scrutinize the accountability of parameterized and parameter-free single-hybrid (SH) and double-hybrid (DH) functionals through the two formalisms, full time-dependent density functional theory (TD-DFT) and Tamm-Dancoff approximation (TDA), for the estimation of photophysical properties like absorption energy, emission energy, zero-zero transition energy, and singlet-triplet energy splitting of TADF molecules. According to our detailed analyses on the performance of SHs based on TD-DFT and TDA, the TDA-based parameter-free SH functionals, PBE0 and TPSS0, with one-third of exact-like exchange turned out to be the best performers in comparison to other functionals from various rungs to reproduce the experimental data of the benchmarked set. Such affordable SH approximations can thus be employed to predict and design the TADF molecules with low singlet-triplet energy gaps for OLED applications. From another perspective, considering this point that both the nonlocal exchange and correlation are essential for a more reliable description of large charge-transfer excited states, applicability of the functionals incorporating these terms, namely, parameterized and parameter-free DHs, has also been evaluated. Perusing the role of exact-like exchange, perturbative-like correlation, solvent effects, and other related factors, we find that the parameterized functionals B2π-PLYP and B2GP-PLYP and the parameter-free models PBE-CIDH and PBE-QIDH have respectable performance with respect to others. Lastly, besides the recommendation of reliable computational protocols for the purpose, hopefully this study can pave the way toward further developments of other SHs and DHs for theoretical explorations in the field of OLEDs technology.
Collapse
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, College of Sciences, Shiraz University, Shiraz, Iran
| | - Niloofar Karimi
- Department of Chemistry, College of Sciences, Shiraz University, Shiraz, Iran
| |
Collapse
|
39
|
Kołodziejska R, Kwit M, Studzińska R, Jelecki M. Enantio- and diastereoselective acylation of prochiral hydroxyl group of pyrimidine acyclonucleosides. ACTA ACUST UNITED AC 2016. [DOI: 10.1016/j.molcatb.2016.07.013] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
40
|
Diedrich D, Moita AJR, Rüther A, Frieg B, Reiss GJ, Hoeppner A, Kurz T, Gohlke H, Lüdeke S, Kassack MU, Hansen FK. α-Aminoxy Oligopeptides: Synthesis, Secondary Structure, and Cytotoxicity of a New Class of Anticancer Foldamers. Chemistry 2016; 22:17600-17611. [PMID: 27573537 DOI: 10.1002/chem.201602521] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2016] [Indexed: 11/11/2022]
Abstract
α-Aminoxy peptides are peptidomimetic foldamers with high proteolytic and conformational stability. To gain an improved synthetic access to α-aminoxy oligopeptides we used a straightforward combination of solution- and solid-phase-supported methods and obtained oligomers that showed a remarkable anticancer activity against a panel of cancer cell lines. We solved the first X-ray crystal structure of an α-aminoxy peptide with multiple turns around the helical axis. The crystal structure revealed a right-handed 28 -helical conformation with precisely two residues per turn and a helical pitch of 5.8 Å. By 2D ROESY experiments, molecular dynamics simulations, and CD spectroscopy we were able to identify the 28 -helix as the predominant conformation in organic solvents. In aqueous solution, the α-aminoxy peptides exist in the 28 -helical conformation at acidic pH, but exhibit remarkable changes in the secondary structure with increasing pH. The most cytotoxic α-aminoxy peptides have an increased propensity to take up a 28 -helical conformation in the presence of a model membrane. This indicates a correlation between the 28 -helical conformation and the membranolytic activity observed in mode of action studies, thereby providing novel insights in the folding properties and the biological activity of α-aminoxy peptides.
Collapse
Affiliation(s)
- Daniela Diedrich
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Ana J Rodrigues Moita
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Anja Rüther
- Institut für Pharmazeutische Wissenschaften, Albert-Ludwigs-Universität Freiburg, 79104, Freiburg, Germany
| | - Benedikt Frieg
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Guido J Reiss
- Institut für Anorganische Chemie und Strukturchemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Astrid Hoeppner
- X-Ray Facility and Crystal Farm, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Thomas Kurz
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Holger Gohlke
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Steffen Lüdeke
- Institut für Pharmazeutische Wissenschaften, Albert-Ludwigs-Universität Freiburg, 79104, Freiburg, Germany
| | - Matthias U Kassack
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| | - Finn K Hansen
- Institut für Pharmazeutische und Medizinische Chemie, Heinrich-Heine-Universität Düsseldorf, 40225, Düsseldorf, Germany
| |
Collapse
|
41
|
Sanna D, Sciortino G, Ugone V, Micera G, Garribba E. Nonoxido V(IV) Complexes: Prediction of the EPR Spectrum and Electronic Structure of Simple Coordination Compounds and Amavadin. Inorg Chem 2016; 55:7373-87. [PMID: 27399275 DOI: 10.1021/acs.inorgchem.6b00409] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Density functional theory (DFT) calculations of the (51)V hyperfine coupling (HFC) tensor A have been completed for 20 "bare" V(IV) complexes with different donor sets, electric charges, and coordination geometries. Calculations were performed with ORCA and Gaussian software, using functionals BP86, TPSS0, B1LYP, PBE0, B3LYP, B3P, B3PW, O3LYP, BHandHLYP, BHandH, and B2PLYP. Among the basis sets, 6-311g(d,p), 6-311++g(d,p), VTZ, cc-pVTZ, def2-TZVPP, and the "core properties" CP(PPP) were tested. The experimental Aiso and Ai (where i = x or z, depending on the geometry and electronic structure of V(IV) complex) were compared with the values calculated by DFT methods. The results indicated that, based on the mean absolute percentage deviation (MAPD), the best functional to predict Aiso or Ai is the double hybrid B2PLYP. With this functional and the basis set VTZ, it is possible to predict the Aiso and Az of the EPR spectrum of amavadin with deviations of -1.1% and -2.0% from the experimental values. The results allowed us to divide the spectra of nonoxido V(IV) compounds in three types-called "type 1", "type 2", and "type 3", characterized by different composition of the singly occupied molecular orbital (SOMO) and relationship between the values of Ax, Ay, and Az. For "type 1" spectra, Az ≫ Ax ≈ Ay and Az is in the range of (135-155) × 10(-4) cm(-1); for "type 2" spectra, Ax ≈ Ay ≫ Az and Ax ≈ Ay are in the range of (90-120) × 10(-4) cm(-1); and for the intermediate spectra of "type 3", Az > Ay > Ax or Ax > Ay > Az, with Az or Ax values in the range of (120-135) × 10(-4) cm(-1). The electronic structure of the V(IV) species was also discussed, and the results showed that the values of Ax or Az are correlated with the percent contribution of V-dxy orbital in the SOMO. Similarly to V(IV)O species, for amavadin the SOMO is based mainly on the V-dxy orbital, and this accounts for the large experimental value of Az (153 × 10(-4) cm(-1)).
Collapse
Affiliation(s)
- Daniele Sanna
- Istituto di Chimica Biomolecolare, Consiglio Nazionale delle Ricerche, UOS di Sassari, Trav. La Crucca 3, I-07040 Sassari, Italy
| | - Giuseppe Sciortino
- Dipartimento di Chimica e Farmacia, Università di Sassari , Via Vienna 2, I-07100 Sassari, Italy
| | - Valeria Ugone
- Dipartimento di Chimica e Farmacia, Università di Sassari , Via Vienna 2, I-07100 Sassari, Italy
| | - Giovanni Micera
- Dipartimento di Chimica e Farmacia, Università di Sassari , Via Vienna 2, I-07100 Sassari, Italy
| | - Eugenio Garribba
- Dipartimento di Chimica e Farmacia, Università di Sassari , Via Vienna 2, I-07100 Sassari, Italy
| |
Collapse
|
42
|
Liu Y, Cerezo J, Mazzeo G, Lin N, Zhao X, Longhi G, Abbate S, Santoro F. Vibronic Coupling Explains the Different Shape of Electronic Circular Dichroism and of Circularly Polarized Luminescence Spectra of Hexahelicenes. J Chem Theory Comput 2016; 12:2799-819. [DOI: 10.1021/acs.jctc.6b00109] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Yanli Liu
- State
Key Laboratory of Crystal Materials, Shandong University, 250100 Jinan, Shandong, People’s Republic of China
- CNR − Consiglio
Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo
Metallici (ICCOM-CNR), UOS di Pisa,
Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Javier Cerezo
- CNR − Consiglio
Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo
Metallici (ICCOM-CNR), UOS di Pisa,
Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Giuseppe Mazzeo
- Dipartimento
di Medicina Molecolare e Traslazionale, Universitá di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Na Lin
- State
Key Laboratory of Crystal Materials, Shandong University, 250100 Jinan, Shandong, People’s Republic of China
| | - Xian Zhao
- State
Key Laboratory of Crystal Materials, Shandong University, 250100 Jinan, Shandong, People’s Republic of China
| | - Giovanna Longhi
- Dipartimento
di Medicina Molecolare e Traslazionale, Universitá di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Sergio Abbate
- Dipartimento
di Medicina Molecolare e Traslazionale, Universitá di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Fabrizio Santoro
- CNR − Consiglio
Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo
Metallici (ICCOM-CNR), UOS di Pisa,
Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
43
|
On the performance of time-dependent double-hybrid density functionals for description of absorption and emission spectra of heteroaromatic compounds. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1838-4] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
44
|
Komjáti B, Urai Á, Hosztafi S, Kökösi J, Kováts B, Nagy J, Horváth P. Systematic study on the TD-DFT calculated electronic circular dichroism spectra of chiral aromatic nitro compounds: A comparison of B3LYP and CAM-B3LYP. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2016; 155:95-102. [PMID: 26583523 DOI: 10.1016/j.saa.2015.11.002] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/22/2015] [Revised: 10/05/2015] [Accepted: 11/01/2015] [Indexed: 06/05/2023]
Abstract
B3LYP is one of the most widely used functional for the prediction of electronic circular dichroism spectra, however if the studied molecule contains aromatic nitro group computations may fail to produce reliable results. A test set of molecules of known stereochemistry were synthesized to study this phenomenon in detail. Spectra were computed by B3LYP and CAM-B3LYP functionals with 6-311++G(2d,2p) basis set. It was found that the range separated CAM-B3LYP gives better predictions than B3LYP for all test molecules. Fragment population analysis revealed that the nitro groups form highly localized molecule orbitals but the exact composition depends on the functional. CAM-B3LYP allows sufficient spatial overlap between the nitro group and distant parts of the molecule, which is necessary for the accurate description of excited states especially for charge transfer states. This phenomenon and the synthesized test molecules can be used to benchmark theoretical methods as well as to help the development of new functionals intended for spectroscopical studies.
Collapse
Affiliation(s)
- Balázs Komjáti
- Department of Organic Chemistry and Technology, Budapest University of Technology and Economics, Szent Gellért tér 4, Budapest 1111, Hungary.
| | - Ákos Urai
- Department of Pharmaceutical Chemistry, Semmelweis University, Hőgyes Endre utca 7, Budapest 1092, Hungary
| | - Sándor Hosztafi
- Department of Pharmaceutical Chemistry, Semmelweis University, Hőgyes Endre utca 7, Budapest 1092, Hungary
| | - József Kökösi
- Department of Pharmaceutical Chemistry, Semmelweis University, Hőgyes Endre utca 7, Budapest 1092, Hungary
| | - Benjámin Kováts
- Department of Organic Chemistry and Technology, Budapest University of Technology and Economics, Szent Gellért tér 4, Budapest 1111, Hungary
| | - József Nagy
- Department of Organic Chemistry and Technology, Budapest University of Technology and Economics, Szent Gellért tér 4, Budapest 1111, Hungary
| | - Péter Horváth
- Department of Pharmaceutical Chemistry, Semmelweis University, Hőgyes Endre utca 7, Budapest 1092, Hungary
| |
Collapse
|
45
|
Gajewy J, Szymkowiak J, Kwit M. Fine tuning of molecular and supramolecular properties of simple trianglimines – the role of the functional group. RSC Adv 2016. [DOI: 10.1039/c6ra06095a] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Readily available chiral, triangular poly-azamacrocycles can act as receptors, low molecular weight supergelators and ligands in asymmetric synthesis.
Collapse
Affiliation(s)
- J. Gajewy
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznan
- Poland
- Wielkopolska Center for Advanced Technologies (WCAT)
| | - J. Szymkowiak
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznan
- Poland
- Wielkopolska Center for Advanced Technologies (WCAT)
| | - M. Kwit
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznan
- Poland
- Wielkopolska Center for Advanced Technologies (WCAT)
| |
Collapse
|
46
|
Markworth PB, Adamson BD, Coughlan NJA, Goerigk L, Bieske EJ. Photoisomerization action spectroscopy: flicking the protonated merocyanine–spiropyran switch in the gas phase. Phys Chem Chem Phys 2015; 17:25676-88. [DOI: 10.1039/c5cp01567g] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular ions drifting through buffer gas are exposed to tunable laser light causing a detectable change in their collision cross section and arrival time at the ion detector. The effect can be exploited to obtain action electronic spectra for selected molecular isomers in the gas phase.
Collapse
Affiliation(s)
| | - Brian D. Adamson
- School of Chemistry
- The University of Melbourne
- Melbourne
- Australia
| | | | - Lars Goerigk
- School of Chemistry
- The University of Melbourne
- Melbourne
- Australia
| | - Evan J. Bieske
- School of Chemistry
- The University of Melbourne
- Melbourne
- Australia
| |
Collapse
|
47
|
Tang HY, Zhang Q, Gao YQ, Zhang AL, Gao JM. Miniolins A–C, novel isomeric furanones induced by epigenetic manipulation of Penicillium minioluteum. RSC Adv 2015. [DOI: 10.1039/c4ra11712c] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
Cultivation of Penicillium minioluteum with a DNA methyltransferase inhibitor led to the isolation of a novel type of aspertetronin dimer, miniolins A–C, along with their precursor. The dimeric absolute configurations were assigned by chiral HPLC and ECD calculations.
Collapse
Affiliation(s)
- Hao-Yu Tang
- College of Science
- Northwest A&F University
- Yangling 712100
- China
| | - Qiang Zhang
- College of Science
- Northwest A&F University
- Yangling 712100
- China
| | - Yu-Qi Gao
- College of Science
- Northwest A&F University
- Yangling 712100
- China
| | - An-Ling Zhang
- College of Science
- Northwest A&F University
- Yangling 712100
- China
| | - Jin-Ming Gao
- College of Science
- Northwest A&F University
- Yangling 712100
- China
| |
Collapse
|
48
|
Prusinowska N, Bendzińska-Berus W, Szymkowiak J, Warżajtis B, Gajewy J, Jelecki M, Rychlewska U, Kwit M. Double helicity induction in chiral bis(triphenylacetamides). RSC Adv 2015. [DOI: 10.1039/c5ra15554a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Helicity induction in a highly flexible trityl chromophore is not only due to the presence of the neighboring stereogenic center(s).
Collapse
Affiliation(s)
| | | | - Joanna Szymkowiak
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznań
- Poland
- Wielkopolska Center for Advanced Technologies (WCAT)
| | - Beata Warżajtis
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznań
- Poland
| | - Jadwiga Gajewy
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznań
- Poland
- Wielkopolska Center for Advanced Technologies (WCAT)
| | - Maciej Jelecki
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznań
- Poland
- Wielkopolska Center for Advanced Technologies (WCAT)
| | | | - Marcin Kwit
- Department of Chemistry
- Adam Mickiewicz University
- 61 614 Poznań
- Poland
- Institute of Organic Chemistry
| |
Collapse
|
49
|
Milde B, Leibeling M, Pawliczek M, Grunenberg J, Jones PG, Werz DB. π-Helicenes Truncated to a Minimum: Access Through a Domino Approach Involving Multiple Carbopalladations and a Stille Coupling. Angew Chem Int Ed Engl 2014; 54:1331-5. [DOI: 10.1002/anie.201408637] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2014] [Revised: 11/07/2014] [Indexed: 11/10/2022]
|
50
|
Milde B, Leibeling M, Pawliczek M, Grunenberg J, Jones PG, Werz DB. π-Helicene auf ein Minimum reduziert: Zugang mithilfe einer mehrfachen Domino-Carbopalladierungs-Stille-Sequenz. Angew Chem Int Ed Engl 2014. [DOI: 10.1002/ange.201408637] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|