1
|
Grotjahn R, Purnomo J, Jin D, Lutfi N, Furche F. Chemically Accurate Singlet-Triplet Gaps of Arylcarbenes from Local Hybrid Density Functionals. J Phys Chem A 2024; 128:6046-6060. [PMID: 39012067 DOI: 10.1021/acs.jpca.4c02852] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/17/2024]
Abstract
Singlet-triplet (ST) gaps are key descriptors of carbenes, because their properties and reactivity are strongly spin-dependent. However, the theoretical prediction of ST gaps is challenging and generally thought to require elaborate correlated wave function methods or double-hybrid density functionals. By evaluating two recent test sets of arylcarbenes (AC12 and AC18), we show that local hybrid functionals based on the "common t" local mixing function (LMF) model achieve mean absolute errors below 1 kcal/mol at a computational cost only slightly higher than that of global hybrid functionals. An analysis of correlation contributions to the ST gaps suggests that the accuracy of the common t-LMF model is mainly due to an improved description of nondynamical correlation which, unlike exchange, is not additive in each spin-channel. Although spin-nonadditivity can be achieved using the local spin polarization alone, using the "common", i.e., spin-unresolved, iso-orbital indicator t for constructing the LMF is found to be critical for consistent accuracy in ST gaps of arylcarbenes. The results support the view of LHs as vehicles to improve the description of nondynamical correlation rather than sophisticated exchange mixing approaches.
Collapse
Affiliation(s)
- Robin Grotjahn
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Justin Purnomo
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Dayun Jin
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Nicolas Lutfi
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
2
|
Sülzner N, Hättig C. Role of Singles Amplitudes in ADC(2) and CC2 for Low-Lying Electronically Excited States. J Chem Theory Comput 2024; 20:2462-2474. [PMID: 38449383 DOI: 10.1021/acs.jctc.3c01355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/08/2024]
Abstract
The closely related second-order methods CC2 and ADC(2) usually perform very similarly for single excitations of organic molecules. However, as rationalized in this work, significant deviations between these two methods can arise if the ground state and a low-lying singly excited state arise from a strong coupling between their leading configurations. Such a configuration mixing is partially accounted for in CC2 through the ground-state singles amplitudes but is omitted in ADC(2). This can cause unusual deviations between the results obtained with these methods. In this work, we study how severe this effect can become at the example of two solvatochromic dyes: the negatively solvatochromic betaine dye N1-tBu and the positively solvatochromic bithiophene P1. These two dyes allow one to study the limits of both small and somewhat larger excitation energies and configuration mixing by tuning the S0 → S1 transition energy through the polarity of the environment. Higher-level calculations at the CC3 level provide information on the accuracy of ADC(2) and CC2 in these cases. The most extreme deviation between ADC(2) and CC2 is found for N1-tBu in vacuum, where the ADC(2) result is 0.45 eV below that of CC2. In this case, the methodical error of CC2 with respect to CC3 is only 0.05 eV. With increasing excitation energy in polar solvents, the CC2-ADC(2) deviation decreases and reaches a value of only 0.15 eV. For P1, which has larger excitation energies, these effects are reversed due to the opposite solvatochromism but also smaller in magnitude: the deviation increases from 0.08 eV in vacuum to 0.16 eV in the so-called conductor limit of the continuum solvation model. Although for these two dyes larger deviations are observed for smaller excitation energies, the extent of configuration mixing does not generally correlate with only the size of excitation energy. For example, s-triazine (0.15 eV), formamide (0.19 eV), and formaldehyde (0.23 eV) also show large deviations between CC2 and ADC(2) despite their much higher excitation energies compared to those of N1-tBu and P1.
Collapse
Affiliation(s)
- Niklas Sülzner
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
3
|
Ju CW, Shen Y, French EJ, Yi J, Bi H, Tian A, Lin Z. Accurate Electronic and Optical Properties of Organic Doublet Radicals Using Machine Learned Range-Separated Functionals. J Phys Chem A 2024. [PMID: 38382058 DOI: 10.1021/acs.jpca.3c07437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
Luminescent organic semiconducting doublet-spin radicals are unique and emergent optical materials because their fluorescent quantum yields (Φfl) are not compromised by the spin-flipping intersystem crossing (ISC) into a dark high-spin state. The multiconfigurational nature of these radicals challenges their electronic structure calculations in the framework of single-reference density functional theory (DFT) and introduces room for method improvement. In the present study, we extended our earlier development of ML-ωPBE [J. Phys. Chem. Lett., 2021, 12, 9516-9524], a range-separated hybrid (RSH) exchange-correlation (XC) functional constructed using the stacked ensemble machine learning (SEML) algorithm, from closed-shell organic semiconducting molecules to doublet-spin organic semiconducting radicals. We assessed its performance for a new test set of 64 doublet-spin radicals from five categories while placing all previously compiled 3926 closed-shell molecules in the new training set. Interestingly, ML-ωPBE agrees with the nonempirical OT-ωPBE functional regarding the prediction of the molecule-dependent range-separation parameter (ω), with a small mean absolute error (MAE) of 0.0197 a0-1, but saves the computational cost by 2.46 orders of magnitude. This result demonstrates an outstanding domain adaptation capacity of ML-ωPBE for diverse organic semiconducting species. To further assess the predictive power of ML-ωPBE in experimental observables, we also applied it to evaluate absorption and fluorescence energies (Eabs and Efl) using linear-response time-dependent DFT (TDDFT), and we compared its behavior with nine popular XC functionals. For most radicals, ML-ωPBE reproduces experimental measurements of Eabs and Efl with small MAEs of 0.299 and 0.254 eV, only marginally different from those of OT-ωPBE. Our work illustrates a successful extension of the SEML framework from closed-shell molecules to doublet-spin radicals and will open the venue for calculating optical properties for organic semiconductors using single-reference TDDFT.
Collapse
Affiliation(s)
- Cheng-Wei Ju
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, United States
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| | - Yili Shen
- Manning College of Information and Computer Sciences, University of Massachusetts, Amherst, Massachusetts 01003, United States
- Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Ethan J French
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, United States
- Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003, United States
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts 02129, United States
| | - Jun Yi
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, United States
- Department of Chemistry, Wake Forest University, Winston-Salem, North Carolina 27109, United States
| | - Hongshan Bi
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, United States
| | - Aaron Tian
- Manning College of Information and Computer Sciences, University of Massachusetts, Amherst, Massachusetts 01003, United States
- Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts 01003, United States
| | - Zhou Lin
- Department of Chemistry, University of Massachusetts, Amherst, Massachusetts 01003, United States
| |
Collapse
|
4
|
Jacquemin D, Kossoski F, Gam F, Boggio-Pasqua M, Loos PF. Reference Vertical Excitation Energies for Transition Metal Compounds. J Chem Theory Comput 2023. [PMID: 37965941 DOI: 10.1021/acs.jctc.3c01080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.
Collapse
Affiliation(s)
- Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Franck Gam
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| |
Collapse
|
5
|
Grotjahn R. Learning from the 4-(dimethylamino)benzonitrile twist: Two-parameter range-separated local hybrid functional with high accuracy for triplet and charge-transfer excitations. J Chem Phys 2023; 159:174102. [PMID: 37909451 DOI: 10.1063/5.0173701] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 10/06/2023] [Indexed: 11/03/2023] Open
Abstract
The recent ωLH22t range-separated local hybrid (RSLH) is shown to provide outstanding accuracy for the notorious benchmark problem of the two lowest excited-state potential energy curves for the amino group twist in 4-(dimethylamino)benzonitrile (DMABN). However, the design of ωLH22t as a general-purpose functional resulted in less convincing performance for triplet excitations, which is an important advantage of previous LHs. Furthermore, ωLH22t uses 8 empirical parameters to achieve broad accuracy. In this work, the RSLH ωLH23ct-sir is constructed with minimal empiricism by optimizing its local mixing function prefactor and range-separation parameter for only 8 excitation energies. ωLH23ct-sir maintains the excellent performance of ωLH22t for the DMABN twist and charge-transfer benchmarks but significantly improves the errors for triplet excitation energies (0.17 vs 0.24 eV). Additional test calculations for the AE6BH6 thermochemistry test set and large dipole moment and static polarizability test sets confirm that the focus on excitation energies in the optimization of ωLH23ct-sir has not caused any dramatic errors for ground-state properties. Although ωLH23ct-sir cannot replace ωLH22t as a general-purpose functional, it is preferable for problems requiring a universally good description of localized and charge-transfer excitations of both singlet and triplet multiplicity. Current limitations on the application of ωLH23ct-sir and other RSLHs to the study of singlet-triplet gaps of emitters for thermally activated delayed fluorescence are discussed. This work also includes the first systematic analysis of the influence of the local mixing function prefactor and the range-separation parameter in an RSLH on different types of excitations.
Collapse
Affiliation(s)
- Robin Grotjahn
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
6
|
Talbot JJ, Arias-Martinez JE, Cotton SJ, Head-Gordon M. Fantastical excited state optimized structures and where to find them. J Chem Phys 2023; 159:171102. [PMID: 37916588 DOI: 10.1063/5.0172015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Accepted: 10/16/2023] [Indexed: 11/03/2023] Open
Abstract
The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.
Collapse
Affiliation(s)
- Justin J Talbot
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Juan E Arias-Martinez
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Stephen J Cotton
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
7
|
Tapavicza E, Reutershan T, Thompson T. Ab Initio Simulation of the Ultrafast Circular Dichroism Spectrum of Provitamin D Ring-Opening. J Phys Chem Lett 2023; 14:5061-5068. [PMID: 37227143 PMCID: PMC10240533 DOI: 10.1021/acs.jpclett.3c00862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Accepted: 05/03/2023] [Indexed: 05/26/2023]
Abstract
We present a method to simulate ultrafast pump-probe time-resolved circular dichroism (TRCD) spectra based on time-dependent density functional theory trajectory surface hopping. The method is applied to simulate the TRCD spectrum along the photoinduced ring-opening of provitamin D. Simulations reveal that the initial decay of the signal is due to excited state relaxation, forming the rotationally flexible previtamin D. We further show that oscillations in the experimental TRCD spectrum arise from isomerizations between previtamin D rotamers with different chirality, which are associated with the helical conformation of the triene unit. We give a detailed description of the formation dynamics of different rotamers, playing a key role in the natural regulation of vitamin D photosynthesis. Going beyond the sole extraction of decay rates, simulations greatly increase the amount of information that can be retrieved from ultrafast TRCD, making it a sensitive tool to unravel details in the subpicosecond dynamics of photoinduced chirality changes.
Collapse
Affiliation(s)
- Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, United States
| | - Trevor Reutershan
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, United States
| | - Travis Thompson
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, United States
| |
Collapse
|
8
|
Bralick AK, Mitchell EC, Doner AC, Webb AR, Christianson MG, Turney JM, Rotavera B, Schaefer HF. Simulation of the VUV Absorption Spectra of Oxygenates and Hydrocarbons: A Joint Theoretical-Experimental Study. J Phys Chem A 2023; 127:3743-3756. [PMID: 37097841 PMCID: PMC10165657 DOI: 10.1021/acs.jpca.2c07743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Vacuum UV absorption spectroscopy is regularly used to provide unambiguous identification of a target species, insight into the electronic structure of molecules, and quantitative species concentrations. As molecules of interest have become more complex, theoretical spectra have been used in tandem with laboratory spectroscopic analysis or as a replacement when experimental data is unavailable. However, it is difficult to determine which theoretical methodologies can best simulate experiment. This study examined the performance of EOM-CCSD and 10 TD-DFT functionals (B3LYP, BH&HLYP, BMK, CAM-B3LYP, HSE, M06-2X, M11, PBE0, ωB97X-D, and X3LYP) to produce reliable vacuum UV absorption spectra for 19 small oxygenates and hydrocarbons using vertical excitation energies. The simulated spectra were analyzed against experiment using both a qualitative analysis and quantitative metrics, including cosine similarity, relative integral change, mean signed error, and mean absolute error. Based on our ranking system, it was determined that M06-2X was consistently the top performing TD-DFT method with BMK, CAM-B3LYP, and ωB97X-D also producing reliable spectra for these small combustion species.
Collapse
Affiliation(s)
- Addison K Bralick
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
- Center for Computational Quantum Chemistry, University of Georgia, 1004 Cedar Street, Athens, Georgia 30602, United States
| | - Erica C Mitchell
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
- Center for Computational Quantum Chemistry, University of Georgia, 1004 Cedar Street, Athens, Georgia 30602, United States
| | - Anna C Doner
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
| | - Annabelle R Webb
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
| | - Matthew G Christianson
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
| | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia, 1004 Cedar Street, Athens, Georgia 30602, United States
| | - Brandon Rotavera
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
- College of Engineering, University of Georgia, 597 D.W. Brooks Drive, Athens, Georgia 30602, United States
| | - Henry F Schaefer
- Department of Chemistry, University of Georgia, 302 East Campus Road, Athens, Georgia 30602, United States
- Center for Computational Quantum Chemistry, University of Georgia, 1004 Cedar Street, Athens, Georgia 30602, United States
| |
Collapse
|
9
|
Bertoni AI, Sánchez CG. Data-driven approach for benchmarking DFTB-approximate excited state methods. Phys Chem Chem Phys 2023; 25:3789-3798. [PMID: 36645084 DOI: 10.1039/d2cp04979a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this work we propose a chemically-informed data-driven approach to benchmark the approximate density-functional tight-binding (DFTB) excited state (ES) methods that are currently available within the DFTB+ suite. By taking advantage of the large volume of low-detail ES-data in the machine learning (ML) dataset, QM8, we were able to extract valuable insights regarding the limitations of the benchmarked methods in terms of the approximations made to the parent formalism, density-functional theory (DFT), while providing recommendations on how to overcome them. For this benchmark, we compared the first singlet-singlet vertical excitation energies (E1) predicted by the DFTB-approximate methods with predictions of less approximate methods from the reference ML-dataset. For the nearly 21800 organic molecules in the GDB-8 chemical space, we were able to identify clear trends in the E1 prediction error distributions, with respect to second-order approximate coupled cluster (CC2), showing a strong dependence on chemical identity.
Collapse
Affiliation(s)
- Andrés I Bertoni
- Instituto Interdisciplinario de Ciencias Básicas (ICB-CONICET), Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza 5502, Argentina.
| | - Cristián G Sánchez
- Instituto Interdisciplinario de Ciencias Básicas (ICB-CONICET), Universidad Nacional de Cuyo, Padre Jorge Contreras 1300, Mendoza 5502, Argentina.
| |
Collapse
|
10
|
Abou Taka A, Corzo HH, Pribram Jones A, Hratchian HP. Good Vibrations: Calculating Excited-State Frequencies Using Ground-State Self-Consistent Field Models. J Chem Theory Comput 2022; 18:7286-7297. [PMID: 36445860 DOI: 10.1021/acs.jctc.2c00672] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
The use of Δ-self-consistent field (SCF) approaches for studying excited electronic states has received a renewed interest in recent years. In this work, the use of this scheme for calculating excited-state vibrational frequencies is examined. Results from Δ-SCF calculations for a set of representative molecules are compared with those obtained using configuration interaction with single substitutions (CIS) and time-dependent density functional theory (TD-DFT) methods. The use of an approximate spin purification model is also considered for cases where the excited-state SCF solution is spin-contaminated. The results of this work demonstrate that an SCF-based description of an excited-state potential energy surface can be an accurate and cost-effective alternative to CIS and TD-DFT methods.
Collapse
Affiliation(s)
- Ali Abou Taka
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California95343, United States.,Combustion Research Facility, Sandia National Laboratories, Livermore, California94550, United States
| | - Hector H Corzo
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California95343, United States.,National Center for Computational Sciences, Oak Ridge Leadership Computing Facility, Oak Ridge National laboratory, Oak Ridge, Tennessee37831-6012, United States
| | - Aurora Pribram Jones
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California95343, United States
| | - Hrant P Hratchian
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, California95343, United States
| |
Collapse
|
11
|
Gronowski M, Kołos R. A DFT Study on the Excited Electronic States of Cyanopolyynes: Benchmarks and Applications. Molecules 2022; 27:molecules27185829. [PMID: 36144567 PMCID: PMC9500640 DOI: 10.3390/molecules27185829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/01/2022] [Accepted: 09/02/2022] [Indexed: 11/16/2022] Open
Abstract
Highly unsaturated chain molecules are interesting due to their potential application as nanowires and occurrence in interstellar space. Here, we focus on predicting the electronic spectra of polyynic nitriles HC2m+1N (m = 0–13) and dinitriles NC2n+2N (n = 0–14). The results of time-dependent density functional theory (TD-DFT) calculations are compared with the available gas-phase and noble gas matrix experimental data. We assessed the performance of fifteen functionals and five basis sets for reproducing (i) vibrationless electronic excitation energies and (ii) vibrational frequencies in the singlet excited states. We found that the basis sets of at least triple-ζ quality were necessary to describe the long molecules with alternate single and triple bonds. Vibrational frequency scaling factors are similar for the ground and excited states. The benchmarked spectroscopic parameters were shown to be acceptably reproduced with adequately chosen functionals, in particular ωB97X, CAM-B3LYP, B3LYP, B971, and B972. Select functionals were applied to study the electronic excitation of molecules up to HC27N and C30N2. It is demonstrated that optical excitation leads to a shift from the polyyne- to a cumulene-like electronic structure.
Collapse
|
12
|
Hehn AS, Sertcan B, Belleflamme F, Chulkov SK, Watkins MB, Hutter J. Excited-State Properties for Extended Systems: Efficient Hybrid Density Functional Methods. J Chem Theory Comput 2022; 18:4186-4202. [PMID: 35759470 PMCID: PMC9281608 DOI: 10.1021/acs.jctc.2c00144] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Time-dependent density functional theory has become state-of-the-art for describing photophysical and photochemical processes in extended materials because of its affordable cost. The inclusion of exact exchange was shown to be essential for the correct description of the long-range asymptotics of electronic interactions and thus a well-balanced description of valence, Rydberg, and charge-transfer excitations. Several approaches for an efficient treatment of exact exchange have been established for the ground state, while implementations for excited-state properties are rare. Furthermore, the high computational costs required for excited-state properties in comparison to ground-state computations often hinder large-scale applications on periodic systems with hybrid functional accuracy. We therefore propose two approximate schemes for improving computational efficiency for the treatment of exact exchange. Within the auxiliary density matrix method (ADMM), exact exchange is estimated using a relatively small auxiliary basis and the introduced basis set incompleteness error is compensated by an exchange density functional correction term. Benchmark results for a test set of 35 molecules demonstrate that the mean absolute error introduced by ADMM is smaller than 0.3 pm for excited-state bond lengths and in the range of 0.02-0.04 eV for vertical excitation, adiabatic excitation, and fluorescence energies. Computational timings for a series of covalent-organic frameworks demonstrate that a speed-up of at least 1 order of magnitude can be achieved for excited-state geometry optimizations in comparison to conventional hybrid functionals. The second method is to use a semiempirical tight binding approximation for both Coulomb and exchange contributions to the excited-state kernel. This simplified Tamm-Dancoff approximation (sTDA) achieves an accuracy comparable to approximated hybrid density functional theory when referring to highly accurate coupled-cluster reference data. We find that excited-state bond lengths deviate by 1.1 pm on average and mean absolute errors in vertical excitation, adiabatic excitation, and fluorescence energies are in the range of 0.2-0.5 eV. In comparison to ADMM-approximated hybrid functional theory, sTDA accelerates the computation of broad-band excitation spectra by 1 order of magnitude, suggesting its potential use for large-scale screening purposes.
Collapse
Affiliation(s)
- Anna-Sophia Hehn
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Beliz Sertcan
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Fabian Belleflamme
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Sergey K. Chulkov
- School
of Mathematics and Physics, University of
Lincoln, Brayford Pool, Lincoln LN67TS, United Kingdom
| | - Matthew B. Watkins
- School
of Mathematics and Physics, University of
Lincoln, Brayford Pool, Lincoln LN67TS, United Kingdom
| | - Jürg Hutter
- Department
of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| |
Collapse
|
13
|
Liang J, Feng X, Hait D, Head-Gordon M. Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four. J Chem Theory Comput 2022; 18:3460-3473. [PMID: 35533317 DOI: 10.1021/acs.jctc.2c00160] [Citation(s) in RCA: 58] [Impact Index Per Article: 29.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.
Collapse
Affiliation(s)
- Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States
| | - Xintian Feng
- Q-Chem Inc., Pleasanton, California 94588, United States
| | - Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
14
|
Alipour M, Izadkhast T. Do any types of double-hybrid models render the correct order of excited state energies in inverted singlet–triplet emitters? J Chem Phys 2022; 156:064302. [DOI: 10.1063/5.0077722] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| | - Tahereh Izadkhast
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| |
Collapse
|
15
|
Gupta A, Chakraborty S, Ghosh D, Ramakrishnan R. Data-driven modeling of S 0 → S 1 excitation energy in the BODIPY chemical space: High-throughput computation, quantum machine learning, and inverse design. J Chem Phys 2021; 155:244102. [PMID: 34972385 DOI: 10.1063/5.0076787] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Derivatives of BODIPY are popular fluorophores due to their synthetic feasibility, structural rigidity, high quantum yield, and tunable spectroscopic properties. While the characteristic absorption maximum of BODIPY is at 2.5 eV, combinations of functional groups and substitution sites can shift the peak position by ±1 eV. Time-dependent long-range corrected hybrid density functional methods can model the lowest excitation energies offering a semi-quantitative precision of ±0.3 eV. Alas, the chemical space of BODIPYs stemming from combinatorial introduction of-even a few dozen-substituents is too large for brute-force high-throughput modeling. To navigate this vast space, we select 77 412 molecules and train a kernel-based quantum machine learning model providing <2% hold-out error. Further reuse of the results presented here to navigate the entire BODIPY universe comprising over 253 giga (253 × 109) molecules is demonstrated by inverse-designing candidates with desired target excitation energies.
Collapse
Affiliation(s)
- Amit Gupta
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad 500107, India
| | - Sabyasachi Chakraborty
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad 500107, India
| | - Debashree Ghosh
- Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Raghunathan Ramakrishnan
- Centre for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad 500107, India
| |
Collapse
|
16
|
Wang J, Durbeej B. Thermal Fluctuations in Conjugation and their Effect on Calculated Excitation Energies: A Case Study on the Astaxanthin Carotenoid. CHEMPHOTOCHEM 2021. [DOI: 10.1002/cptc.202100178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Jun Wang
- Jiangsu Key Laboratory for Chemistry of Low-Dimensional Materials Jiangsu Engineering Laboratory for Environment Functional Materials School of Chemistry and Chemical Engineering Huaiyin Normal University No. 111 West Changjiang Road 223300 Huaian Jiangsu Province China
| | - Bo Durbeej
- Division of Theoretical Chemistry IFM Linköping University 581 83 Linköping Sweden
| |
Collapse
|
17
|
Grotjahn R, Kaupp M. Assessment of hybrid functionals for singlet and triplet excitations: Why do some local hybrid functionals perform so well for triplet excitation energies? J Chem Phys 2021; 155:124108. [PMID: 34598568 DOI: 10.1063/5.0063751] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
The performance of various hybrid density functionals is assessed for 105 singlet and 105 corresponding triplet vertical excitation energies from the QUEST database. The overall lowest mean absolute error is obtained with the local hybrid (LH) functional LH12ct-SsirPW92 with individual errors of 0.11 eV (0.11 eV) for singlet (triplet) n → π* excitations and 0.29 eV (0.17 eV) for π → π* excitations. This is slightly better than with the overall best performing global hybrid M06-2X [n → π*: 0.13 eV (0.17 eV), π → π*: 0.30 eV (0.20 eV)], while most other global and range-separated hybrids and some LHs suffer from the "triplet problem" of time-dependent density functional theory. This is exemplified by correlating the errors for singlet and triplet excitations on a state-by-state basis. The excellent performance of LHs based on a common local mixing function, i.e., an LMF constructed from the spin-summed rather than the spin-resolved semilocal quantities, is systematically investigated by the introduction of a spin-channel interpolation scheme that allows us to continuously modulate the fraction of opposite-spin terms used in the LMF. The correlation of triplet and singlet errors is systematically improved for the n → π* excitations when larger fractions of the opposite-spin-channel are used in the LMF, whereas this effect is limited for the π → π* excitations. This strongly supports a previously made hypothesis that attributes the excellent performance of LHs based on a common LMF to cross-spin-channel nondynamical correlation terms.
Collapse
Affiliation(s)
- Robin Grotjahn
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| |
Collapse
|
18
|
Sofferman DL, Konar A, Mastron JN, Spears KG, Cisneros C, Smith AC, Tapavicza E, Sension RJ. Probing the Formation and Conformational Relaxation of Previtamin D 3 and Analogues in Solution and in Lipid Bilayers. J Phys Chem B 2021; 125:10085-10096. [PMID: 34473504 DOI: 10.1021/acs.jpcb.1c04376] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The photosynthesis of vitamin D3 in mammalian skin results from UV-B irradiation of provitamin D3 (7-dehydrocholesterol, DHC) at ca. 290 nm. Upon return to the ground state, the hexatriene product, previtamin D3, undergoes a conformational equilibration between helical gZg and more planar tZg and tZt forms. The helical gZg forms provide a pathway for the formation of vitamin D3 via a [1,7]-sigmatropic hydrogen shift. Steady state photolysis and UV transient absorption spectroscopy are combined to explore the conformational relaxation of previtamin D3 formed from DHC in isotropic solution and confined to lipid bilayers chosen to model the biological cell membrane. The results are compared with measurements for two analogues: previtamin D2 formed from ergosterol (provitamin D2) and previtamin D3 acetate formed from DHC acetate. The resulting spectral dynamics are interpreted in the context of simulations of optical excitation energy and oscillator strength as a function of conformation. In solution, the relaxation dynamics and steady state product distributions of the three compounds are nearly identical, favoring tZg forms. When confined to lipid bilayers, the heterogeneity and packing forces alter the conformational distributions and enhance the population of a gZg conformer capable of vitamin D formation.
Collapse
Affiliation(s)
- Danielle L Sofferman
- Program in Applied Physics, University of Michigan, 450 Church Street, Ann Arbor, Michigan 48109-1040, United States
| | - Arkaprabha Konar
- Department of Physics, University of Michigan, 450 Church Street, Ann Arbor, Michigan 48109-1040, United States
| | - Joseph N Mastron
- Department of Physics, University of Michigan, 450 Church Street, Ann Arbor, Michigan 48109-1040, United States.,Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109-1055, United States
| | - Kenneth G Spears
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109-1055, United States
| | - Cecilia Cisneros
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840-9507, United States
| | - Adam C Smith
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840-9507, United States
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840-9507, United States
| | - Roseanne J Sension
- Department of Physics, University of Michigan, 450 Church Street, Ann Arbor, Michigan 48109-1040, United States.,Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109-1055, United States
| |
Collapse
|
19
|
Tao Z, Roy S, Schneider PE, Pavošević F, Hammes-Schiffer S. Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies. J Chem Theory Comput 2021; 17:5110-5122. [PMID: 34260237 DOI: 10.1021/acs.jctc.1c00454] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Patrick E Schneider
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
20
|
Tapavicza E, von Rudorff GF, De Haan DO, Contin M, George C, Riva M, von Lilienfeld OA. Elucidating an Atmospheric Brown Carbon Species-Toward Supplanting Chemical Intuition with Exhaustive Enumeration and Machine Learning. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2021; 55:8447-8457. [PMID: 34080853 DOI: 10.1021/acs.est.1c00885] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Brown carbon (BrC) is involved in atmospheric light absorption and climate forcing and can cause adverse health effects. Understanding the formation mechanisms and molecular structure of BrC is of key importance in developing strategies to control its environment and health impact. Structure determination of BrC is challenging, due to the lack of experiments providing molecular fingerprints and the sheer number of molecular candidates with identical mass. Suggestions based on chemical intuition are prone to errors due to the inherent bias. We present an unbiased algorithm, using graph-based molecule generation and machine learning, which can identify all molecular structures of compounds involved in biomass burning and the composition of BrC. We apply this algorithm to C12H12O7, a light-absorbing "test case" molecule identified in chamber experiments on the aqueous photo-oxidation of syringol, a prevalent marker in wood smoke. Of the 260 million molecular graphs, the algorithm leaves only 36,518 (0.01%) as viable candidates matching the spectrum. Although no unique molecular structure is obtained from only a chemical formula and a UV/vis absorption spectrum, we discuss further reduction strategies and their efficacy. With additional data, the method can potentially more rapidly identify isomers extracted from lab and field aerosol particles without introducing human bias.
Collapse
Affiliation(s)
- Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, United States
| | - Guido Falk von Rudorff
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - David O De Haan
- Department of Chemistry and Biochemistry, University of San Diego, 5998 Alcala Park, San Diego, California 92110, United States
| | - Mario Contin
- Facultad de Farmacia y Bioquímica, Departamento de Química Analitica y Fisicoquímica, Universidad de Buenos Aires, Junín 956, Buenos Aires C1113AAD, Argentina
| | - Christian George
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, IRCELYON, 69626 Villeurbanne, France
| | - Matthieu Riva
- Université Lyon, Université Claude Bernard Lyon 1, CNRS, IRCELYON, 69626 Villeurbanne, France
| | - O Anatole von Lilienfeld
- Faculty of Physics, University of Vienna, Kolingasse 14-16, AT-1090 Wien, Austria
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
21
|
Schieschke N, Bold BM, Dohmen PM, Wehl D, Hoffmann M, Dreuw A, Elstner M, Höfener S. Geometry dependence of excitonic couplings and the consequences for configuration-space sampling. J Comput Chem 2021; 42:1402-1418. [PMID: 33993548 DOI: 10.1002/jcc.26552] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Revised: 03/02/2021] [Accepted: 04/06/2021] [Indexed: 12/24/2022]
Abstract
Excitonic coupling plays a key role for the understanding of excitonic energy transport (EET) in, for example, organic photovoltaics. However, the calculation of realistic systems is often beyond the applicability range of accurate wavefunction methods so that lower-scaling semi-empirical methods are used to model EET events. In the present work, the distance and angle dependence of excitonic couplings of dimers of selected organic molecules are evaluated for the semi-empirical long-range corrected density functional based tight binding (LC-DFTB) method and spin opposite scaled second order approximate coupled cluster singles and doubles (SOS-CC2). While semi-empirically scaled methods can lead to slightly increased deviations for excitation energies, the excitonic couplings and their dependence on the dimer geometry are reproduced. LC-DFTB yields a similar accuracy range as density-functional theory (DFT) employing the ωB97X functional while the computation time is reduced by several orders of magnitude. The dependence of the exchange contributions to the excitonic couplings on the dimer geometry is analyzed assessing the calculation of Coulombic excitonic couplings from monomer local excited states only, which reduces the computational effort significantly. The present work is a necessary first step toward the simulation of excitonic energy transport using semi-empirical methods.
Collapse
Affiliation(s)
- Nils Schieschke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Philipp M Dohmen
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Daniel Wehl
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Marvin Hoffmann
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Heidelberg, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany.,Institute of Biological Interfaces (IGB2), Karlsruhe Institute for Technology (KIT), Karlsruhe, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
22
|
Sokolov M, Bold BM, Kranz JJ, Höfener S, Niehaus TA, Elstner M. Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition Energies. J Chem Theory Comput 2021; 17:2266-2282. [PMID: 33689344 DOI: 10.1021/acs.jctc.1c00095] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The absorption and emission of light is a ubiquitous process in chemical and biological processes, making a theoretical description inevitable for understanding and predicting such properties. Although ab initio and DFT methods are capable of describing excited states with good accuracy in many cases, the investigation of dynamical processes and the need to sample the phase space in complex systems often requires methods with reduced computational costs but still sufficient accuracy. In the present work, we report the derivation and implementation of analytical nuclear gradients for time-dependent long-range corrected density functional tight binding (TD-LC-DFTB) in the DFTB+ program. The accuracy of the TD-LC-DFTB potential-energy surfaces is benchmarked for excited-state geometries and adiabatic as well as vertical transition energies. The benchmark set consists of more than 100 organic molecules taken as subsets from available benchmark sets. The reported method yields a mean deviation of 0.31 eV for adiabatic excitation energies with respect to CC2. In order to study more subtle effects, seminumerical second derivatives based on the analytical gradients are employed to simulate vibrationally resolved UV/vis spectra. This extensive test exhibits few problematic cases, which can be traced back to the parametrization of the repulsive potential.
Collapse
Affiliation(s)
- Monja Sokolov
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Julian J Kranz
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Thomas A Niehaus
- Université Claude Bernard Lyon 1, Université Lyon, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany.,Institute of Biological Interfaces (IBG2), Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| |
Collapse
|
23
|
Véril M, Scemama A, Caffarel M, Lipparini F, Boggio‐Pasqua M, Jacquemin D, Loos P. QUESTDB
: A database of highly accurate excitation energies for the electronic structure community. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1517] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Mickaël Véril
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Michel Caffarel
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale University of Pisa Pisa Italy
| | - Martial Boggio‐Pasqua
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| | | | - Pierre‐François Loos
- Laboratoire de Chimie et Physique Quantiques Université de Toulouse, CNRS, UPS Toulouse France
| |
Collapse
|
24
|
Ghahghaey Z, Hekmati M, Darvish Ganji M. Theoretical investigation of phenol adsorption on functionalized graphene using DFT calculations for effective removal of organic contaminants from wastewater. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114777] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
25
|
Ren Y, Li MY, Song YX, Sui MY, Sun GY, Qu XC, Xie P, Lu JL. Refined standards for simulating UV–vis absorption spectra of acceptors in organic solar cells by TD-DFT. J Photochem Photobiol A Chem 2021. [DOI: 10.1016/j.jphotochem.2020.113087] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
26
|
Sarkar R, Boggio-Pasqua M, Loos PF, Jacquemin D. Benchmarking TD-DFT and Wave Function Methods for Oscillator Strengths and Excited-State Dipole Moments. J Chem Theory Comput 2021; 17:1117-1132. [DOI: 10.1021/acs.jctc.0c01228] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
27
|
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
|
28
|
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
|
29
|
Chrayteh A, Blondel A, Loos PF, Jacquemin D. Mountaineering Strategy to Excited States: Highly Accurate Oscillator Strengths and Dipole Moments of Small Molecules. J Chem Theory Comput 2020; 17:416-438. [DOI: 10.1021/acs.jctc.0c01111] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Amara Chrayteh
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Aymeric Blondel
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
30
|
Chrayteh A, Ewels CP, Jacquemin D. TD-DFT and CC2 insights into the dual-emissive behaviour of 2-(2'-hydroxyphenyl)oxazoles core and their derivatives. Phys Chem Chem Phys 2020; 22:25066-25074. [PMID: 33119009 DOI: 10.1039/d0cp04520a] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Two efficient excited state intramolecular proton transfer (ESIPT) dyes based on the hydroxyphenyl-oxazole core and containing one or two triphenylamine donor groups are explored with theoretical tools. These compounds are known to show clear experimental dual emission behaviour, leading to nearly pure white-light emission for one derivative. To probe the excited state properties, we use both Time Dependent Density Functional Theory (TD-DFT) and post Hartree-Fock methods [ADC(2) and CC2] coupled to different solvent models to describe polarisation effects. After validating our theoretical protocol on the two known systems, we design 14 new derivatives with different substitution patterns to quantify the impact of electron accepting and donating groups on the fluorescence spectrum and the ESIPT mechanism. We show that the selected protocol delivers accurate spectroscopic values for the two experimentally-characterised structures, and more importantly, that the relative stabilisation of the keto tautomer depends on the substitution side. Adding donor or acceptor groups to the ESIPT donor moiety favours the formation of the keto form, whereas when placed on the ESIPT accepting side, they tend to preclude ESIPT. Moreover, combining two donor or acceptor substituents generally results in similar ESIPT behaviour as single substitution on one of the two sides: simple additive rules do not apply.
Collapse
Affiliation(s)
- Amara Chrayteh
- Laboratoire CEISAM - UMR 6230 - CNRS - Université de Nantes, Nantes, France.
| | | | | |
Collapse
|
31
|
Barsuk I, Lainé PP, Maurel F, Brémond É. Triangulenium dyes: the comprehensive photo-absorption and emission story of a versatile family of chromophores. Phys Chem Chem Phys 2020; 22:20673-20684. [PMID: 32895673 DOI: 10.1039/d0cp02990d] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
The triangulenium dyes constitute a family of versatile chromophores whose impressive photo-absorption and emission properties are currently highlighted in numerous novel experimental applications. In this investigation, we provide a comprehensive TDDFT characterization of their spectroscopic properties elucidating the origin of their large and complex absorption and emission vibronic spectra spread over the (whole) visible region. More precisely, by benchmarking the performance of 10 commonly-used exchange-correlation density functionals belonging to different classes of approximation, we develop and validate a computational protocol allowing the accurate modeling of both the position and optical line-shape of their vibrationally-resolved absorption and emission band structures. We find that semilocal approximations provide the best estimate of the structure of the vibronic spectra, however they spuriously and strongly underestimate their position. We finally show that global-hybrid density functionals mixing between 20 and 30% of exact-like exchange are an excellent compromise to get a satisfactory estimate of both of these properties.
Collapse
Affiliation(s)
- Irina Barsuk
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France.
| | | | | | - Éric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France.
| |
Collapse
|
32
|
Ou Q, Peng Q, Shuai Z. Toward Quantitative Prediction of Fluorescence Quantum Efficiency by Combining Direct Vibrational Conversion and Surface Crossing: BODIPYs as an Example. J Phys Chem Lett 2020; 11:7790-7797. [PMID: 32787317 DOI: 10.1021/acs.jpclett.0c02054] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Accurate theoretical description of the electronic structure of boron dipyrromethene (BODIPY) molecules has been a challenge, let alone the prediction of fluorescence quantum efficiency. In this Letter, we show that the electronic structures of BODIPYs can be accurately evaluated via the spin-flip time-dependent density functional theory with the B3LYP functional. With the resulting electronic structures, the experimental spectral line shapes of representative BODIPYs are successfully reproduced by our previously developed thermal vibration correlation function method. Most importantly, a two-channel scheme is proposed to describe the internal conversion of S1 to S0 in BODIPYs: channel I via direct vibrational relaxation within the harmonic region and channel II via a distorted S0/S1 minimum energy crossing point well away from the harmonic region. The fluorescence quantum yields are accurately predicted within this two-channel scheme, which can therefore serve as a generalized method for predicting the photophysical parameters of organic fluorescent compounds.
Collapse
Affiliation(s)
- Qi Ou
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| | - Qian Peng
- CAS Key Laboratory of Organic Solids, Institute of Chemistry of the Chinese Academy of Sciences, Zhonguancun Beiyijie 2, Beijing 100190, China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, China
| |
Collapse
|
33
|
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
|
34
|
Grotjahn R, Kaupp M. Validation of Local Hybrid Functionals for Excited States: Structures, Fluorescence, Phosphorescence, and Vibronic Spectra. J Chem Theory Comput 2020; 16:5821-5834. [DOI: 10.1021/acs.jctc.0c00520] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Robin Grotjahn
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
35
|
Liu J, Hättig C, Höfener S. Analytical nuclear gradients for electron-attached and electron-detached states for the second-order algebraic diagrammatic construction scheme combined with frozen-density embedding. J Chem Phys 2020; 152:174109. [DOI: 10.1063/5.0002851] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jing Liu
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, D-76049 Karlsruhe, Germany
| |
Collapse
|
36
|
Wang J, Durbeej B. How accurate are TD‐DFT excited‐state geometries compared to DFT ground‐state geometries? J Comput Chem 2020; 41:1718-1729. [DOI: 10.1002/jcc.26213] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 04/07/2020] [Accepted: 04/09/2020] [Indexed: 01/18/2023]
Affiliation(s)
- Jun Wang
- Division of Theoretical Chemistry, IFMLinköping University Linköping Sweden
- Institut de Química Computacional i Catàlisi, Facultat de CiènciesUniversitat de Girona Girona Spain
| | - Bo Durbeej
- Division of Theoretical Chemistry, IFMLinköping University Linköping Sweden
| |
Collapse
|
37
|
Loos PF, Scemama A, Jacquemin D. The Quest for Highly Accurate Excitation Energies: A Computational Perspective. J Phys Chem Lett 2020; 11:2374-2383. [PMID: 32125872 DOI: 10.1021/acs.jpclett.0c00014] [Citation(s) in RCA: 87] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
We provide an overview of the successive steps that made it possible to obtain increasingly accurate excitation energies with computational chemistry tools, eventually leading to chemically accurate vertical transition energies for small- and medium-size molecules. First, we describe the evolution of ab initio methods employed to define benchmark values, with the original Roos CASPT2 method, then the CC3 method as in the renowned Thiel set, and more recently the resurgence of selected configuration interaction methods. The latter method has been able to deliver consistently, for both single and double excitations, highly accurate excitation energies for small molecules, as well as medium-size molecules with compact basis sets. Second, we describe how these high-level methods and the creation of representative benchmark sets of excitation energies have allowed the fair and accurate assessment of the performance of computationally lighter methods. We conclude by discussing possible future theoretical and technological developments in the field.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
38
|
Narsaria AK, Ruijter JD, Hamlin TA, Ehlers AW, Guerra CF, Lammertsma K, Bickelhaupt FM. Performance of TDDFT Vertical Excitation Energies of Core-Substituted Naphthalene Diimides. J Comput Chem 2020; 41:1448-1455. [PMID: 32142173 PMCID: PMC7317478 DOI: 10.1002/jcc.26188] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Revised: 02/10/2020] [Accepted: 02/22/2020] [Indexed: 01/15/2023]
Abstract
We have evaluated the performance of various density functionals, covering generalized gradient approximation (GGA), global hybrid (GH) and range‐separated hybrid (RSH), using time dependent density functional theory (TDDFT) for computing vertical excitation energies against experimental absorption maximum (λmax) for a set of 10 different core‐substituted naphthalene diimides (cNDI) recorded in dichloromethane. The computed excitation in case of GH PBE0 is most accurate while the trend is most systematic with RSH LCY‐BLYP compared to λmax. We highlight the importance of including solvent effects for optimal agreement with the λmax. Increasing the basis set size from TZ2P to QZ4P has a negligible influence on the computed excitation energies. Notably, RSH CAMY‐B3LYP gave the least error for charge‐transfer excitation. The poorest agreement with λmax is obtained with semi‐local GGA functionals. Use of the optimally‐tuned RSH LCY‐BLYP* is not recommended because of the high computational cost and marginal improvement in results.
Collapse
Affiliation(s)
- Ayush K Narsaria
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Julian D Ruijter
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Trevor A Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Andreas W Ehlers
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands.,Department of Chemistry, University of Johannesburg, Johannesburg, South Africa
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, Leiden, The Netherlands
| | - Koop Lammertsma
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Department of Chemistry, University of Johannesburg, Johannesburg, South Africa
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Institute of Molecules and Materials, Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
39
|
Loos PF, Lipparini F, Boggio-Pasqua M, Scemama A, Jacquemin D. A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules. J Chem Theory Comput 2020; 16:1711-1741. [PMID: 31986042 DOI: 10.1021/acs.jctc.9b01216] [Citation(s) in RCA: 107] [Impact Index Per Article: 26.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Moruzzi 3, 56124 Pisa, Italy
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques, CNRS et Université Toulouse III - Paul Sabatier, 118 route de Narbonne, 31062 Toulouse, France
| | - Denis Jacquemin
- CEISAM Lab, UMR 6230, Université de Nantes, CNRS, F-44000 Nantes, France
| |
Collapse
|
40
|
Loos PF, Jacquemin D. Is ADC(3) as Accurate as CC3 for Valence and Rydberg Transition Energies? J Phys Chem Lett 2020; 11:974-980. [PMID: 31913639 DOI: 10.1021/acs.jpclett.9b03652] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The search for new models rapidly delivering accurate excited-state energies and properties is one of the most active research lines of theoretical chemistry. Along with these developments, the performance of known methods is constantly reassessed on the basis of new benchmark values. In this Letter, we show that the third-order algebraic diagrammatic construction, ADC(3), does not yield transition energies of the same quality as the third-order coupled cluster method, CC3. This is demonstrated by extensive comparisons with several hundred high-quality vertical transition energies obtained with FCI, CCSDTQ, and CCSDT. Direct comparisons with experimental 0-0 energies of small- and medium-size molecules support the same conclusion, which holds for both valence and Rydberg transitions. Considering these results, we introduce a composite approach, ADC(2.5), which consists of averaging the ADC(2) and ADC(3) excitation energies. Although ADC(2.5) does not match the CC3 accuracy, it significantly improves the ADC(3) results, especially for vertical energies.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), CNRS , Université de Toulouse , 31077 Toulouse , France
| | - Denis Jacquemin
- Laboratoire CEISAM UMR UN-CNRS 6230 , Université de Nantes , F-44000 Nantes , France
| |
Collapse
|
41
|
Wang Y, Verma P, Zhang L, Li Y, Liu Z, Truhlar DG, He X. M06-SX screened-exchange density functional for chemistry and solid-state physics. Proc Natl Acad Sci U S A 2020; 117:2294-2301. [PMID: 31953258 PMCID: PMC7007546 DOI: 10.1073/pnas.1913699117] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Screened-exchange hybrid density functionals are especially recommended for solid-state systems because they combine the advantages of hybrid functionals with the correct physics and lower computational cost associated with the attenuation of Hartree-Fock exchange at long range. We present a screened-exchange hybrid functional, M06-SX, that combines the functional form of the local revM06-L functional with a percentage of short-range nonlocal Hartree-Fock exchange. The M06-SX functional gives good results not only for a large set of training data but also for several databases quite different from the training data. The mean unsigned error (MUE) of the M06-SX functional is 2.85 kcal/mol for 418 atomic and molecular energies (AME418) in Minnesota Database 2019, which is better than all five other screened-exchange hybrid functionals tested in this work. The M06-SX functional also gives especially good results for semiconductor band gaps, molecular dissociation energies, noncovalent interactions, barrier heights, and electronic excitation energies excluding long-range charge transfer excitations. For the LC18 lattice constants database, the M06-SX functional gives an MUE of only 0.034 Å. Therefore, the M06-SX functional is well suited for studying molecular chemistry as well as solid-state physics.
Collapse
Affiliation(s)
- Ying Wang
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Pragya Verma
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455-0431
- Chemical Theory Center, University of Minnesota, Minneapolis, MN 55455-0431
- Nanoporous Materials Genome Center, University of Minnesota, Minneapolis, MN 55455-0431
- Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Lujia Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| | - Yaqi Li
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
| | - Zhonghua Liu
- The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
| | - Donald G Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, MN 55455-0431;
- Chemical Theory Center, University of Minnesota, Minneapolis, MN 55455-0431
- Nanoporous Materials Genome Center, University of Minnesota, Minneapolis, MN 55455-0431
- Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
42
|
Shao Y, Mei Y, Sundholm D, Kaila VRI. Benchmarking the Performance of Time-Dependent Density Functional Theory Methods on Biochromophores. J Chem Theory Comput 2020; 16:587-600. [PMID: 31815476 PMCID: PMC7391796 DOI: 10.1021/acs.jctc.9b00823] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Quantum chemical calculations are important for elucidating light-capturing mechanisms in photobiological systems. The time-dependent density functional theory (TDDFT) has become a popular methodology because of its balance between accuracy and computational scaling, despite its problems in describing, for example, charge transfer states. As a step toward systematically understanding the performance of TDDFT calculations on biomolecular systems, we study here 17 commonly used density functionals, including seven long-range separated functionals, and compare the obtained results with excitation energies calculated at the approximate second order coupled-cluster theory level (CC2). The benchmarking set includes the first five singlet excited states of 11 chemical analogues of biochromophores from the green fluorescent protein, rhodopsin/bacteriorhodopsin (Rh/bR), and the photoactive yellow protein. We find that commonly used pure density functionals such as BP86, PBE, M11-L, and hybrid functionals with 20-25% of Hartree-Fock (HF) exchange (B3LYP, PBE0) have a tendency to consistently underestimate vertical excitation energies (VEEs) relative to the CC2 values, whereas hybrid density functionals with around 50% HF exchange such as BHLYP, PBE50, and M06-2X and long-range corrected functionals such as CAM-B3LYP, ωPBE, ωPBEh, ωB97X, ωB97XD, BNL, and M11 overestimate the VEEs. We observe that calculations using the CAM-B3LYP and ωPBEh functionals with 65% and 100% long-range HF exchange, respectively, lead to an overestimation of the VEEs by 0.2-0.3 eV for the benchmarking set. To reduce the systematic error, we introduce here two new empirical functionals, CAMh-B3LYP and ωhPBE0, for which we adjusted the long-range HF exchange to 50%. The introduced parameterization reduces the mean signed average (MSA) deviation to 0.07 eV and the root mean square (rms) deviation to 0.17 eV as compared to the CC2 values. In the present study, TDDFT calculations using the aug-def2-TZVP basis sets, the best performing functionals relative to CC2 are ωhPBE0 (rms = 0.17, MSA = 0.06 eV); CAMh-B3LYP (rms = 0.16, MSA = 0.07 eV); and PBE0 (rms = 0.23, MSA = -0.14 eV). For the popular range-separated CAM-B3LYP functional, we obtain an rms value of 0.31 eV and an MSA value of 0.25 eV, which can be compared with the rms and MSA values of 0.37 and -0.31 eV, respectively, as obtained at the B3LYP level.
Collapse
Affiliation(s)
- Yihan Shao
- Department of Chemistry and Biochemistry , University of Oklahoma , Norman , Oklahoma 73019 , United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science , East China Normal University , Shanghai 200062 , China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China
| | - Dage Sundholm
- Department of Chemistry, Faculty of Science , University of Helsinki , P.O. Box 55, Helsinki FIN-00014 , Finland
| | - Ville R I Kaila
- Department Chemie , Technische Universität München (TUM) , Lichtenbergstrasse 4 , Garching D-85747 , Germany
- Department of Biochemistry and Biophysics , Stockholm University , Stockholm SE-10691 , Sweden
| |
Collapse
|
43
|
Patil V, Sekar N, Padalkar VS, Rajput J, Patil SR, Patil SV. Molecular properties of 5-(1H-Benzo[D]Oxa, thia, imid azole-2-Yl)-2-methyl quinazolin-4-ol fluorescent brighteners: Theoretical and experimental approach. J Mol Struct 2020. [DOI: 10.1016/j.molstruc.2019.126984] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
44
|
Batra K, Zahn S, Heine T. Benchmark of Simplified Time‐Dependent Density Functional Theory for UV–Vis Spectral Properties of Porphyrinoids. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201900192] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Kamal Batra
- Technische Universität DresdenTheoretische Chemie Bergstraße 66c 01062 Dresden Germany
| | - Stefan Zahn
- Leibniz Institut für Oberflächenmodifizierung (IOM) Permoserstraße 15 04318 Leipzig Germany
| | - Thomas Heine
- Technische Universität DresdenTheoretische Chemie Bergstraße 66c 01062 Dresden Germany
- Helmholtz‐Zentrum Dresden‐Rossendorf, Forschungsstelle Leipzig Permoserstraße 15 04318 Leipzig Germany
| |
Collapse
|
45
|
Grathwol CW, Wössner N, Swyter S, Smith AC, Tapavicza E, Hofstetter RK, Bodtke A, Jung M, Link A. Azologization and repurposing of a hetero-stilbene-based kinase inhibitor: towards the design of photoswitchable sirtuin inhibitors. Beilstein J Org Chem 2019; 15:2170-2183. [PMID: 31598174 PMCID: PMC6774072 DOI: 10.3762/bjoc.15.214] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 08/29/2019] [Indexed: 12/31/2022] Open
Abstract
The use of light as an external trigger to change ligand shape and as a result its bioactivity, allows the probing of pharmacologically relevant systems with spatiotemporal resolution. A hetero-stilbene lead resulting from the screening of a compound that was originally designed as kinase inhibitor served as a starting point for the design of photoswitchable sirtuin inhibitors. Because the original stilbenoid structure exerted unfavourable photochemical characteristics it was remodelled to its heteroarylic diazeno analogue. By this intramolecular azologization, the shape of the molecule was left unaltered, whereas the photoswitching ability was improved. As anticipated, the highly analogous compound showed similar activity in its thermodynamically stable stretched-out (E)-form. Irradiation of this isomer triggers isomerisation to the long-lived (Z)-configuration with a bent geometry causing a considerably shorter end-to-end distance. The resulting affinity shifts are intended to enable real-time photomodulation of sirtuins in vitro.
Collapse
Affiliation(s)
- Christoph W Grathwol
- Institute of Pharmacy, University of Greifswald, Friedrich-Ludwig-Jahn-Str. 17, 17489 Greifswald, Germany
| | - Nathalie Wössner
- Institute of Pharmaceutical Sciences, University of Freiburg, Albertstr. 25, 79104 Freiburg, Germany
| | - Sören Swyter
- Institute of Pharmaceutical Sciences, University of Freiburg, Albertstr. 25, 79104 Freiburg, Germany
| | - Adam C Smith
- Department of Chemistry and Biochemistry, California State University Long Beach, 1250 Bellflower Boulevard, Long Beach, CA, 90840 USA
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University Long Beach, 1250 Bellflower Boulevard, Long Beach, CA, 90840 USA
| | - Robert K Hofstetter
- Institute of Pharmacy, University of Greifswald, Friedrich-Ludwig-Jahn-Str. 17, 17489 Greifswald, Germany
| | - Anja Bodtke
- Institute of Pharmacy, University of Greifswald, Friedrich-Ludwig-Jahn-Str. 17, 17489 Greifswald, Germany
| | - Manfred Jung
- Institute of Pharmaceutical Sciences, University of Freiburg, Albertstr. 25, 79104 Freiburg, Germany
| | - Andreas Link
- Institute of Pharmacy, University of Greifswald, Friedrich-Ludwig-Jahn-Str. 17, 17489 Greifswald, Germany
| |
Collapse
|
46
|
Fang C, Durbeej B. Calculation of Free-Energy Barriers with TD-DFT: A Case Study on Excited-State Proton Transfer in Indigo. J Phys Chem A 2019; 123:8485-8495. [DOI: 10.1021/acs.jpca.9b05163] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Changfeng Fang
- Center for Optics Research and Engineering (CORE), Shandong University, Qingdao 266237, China
| | - Bo Durbeej
- Division of Theoretical Chemistry, IFM, Linköping University, SE-581 83 Linköping, Sweden
| |
Collapse
|
47
|
Grotjahn R, Furche F, Kaupp M. Development and Implementation of Excited-State Gradients for Local Hybrid Functionals. J Chem Theory Comput 2019; 15:5508-5522. [DOI: 10.1021/acs.jctc.9b00659] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
48
|
Kaur J, Ospadov E, Staroverov VN. What Is the Accuracy Limit of Adiabatic Linear-Response TDDFT Using Exact Exchange-Correlation Potentials and Approximate Kernels? J Chem Theory Comput 2019; 15:4956-4964. [PMID: 31386366 DOI: 10.1021/acs.jctc.9b00618] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Calculation of vertical excitation energies by the adiabatic linear-response time-dependent density-functional theory (TDDFT) requires static Kohn-Sham potentials and exchange-correlation kernels. When these quantities are derived from standard density-functional approximations (DFA), mean absolute errors (MAE) of the method are known to range from 0.2 eV to over 1 eV, depending on the functional and type of excitation. We investigate how the performance of TDDFT varies when increasingly accurate exchange-correlation potentials derived from Hartree-Fock (HF) and post-HF wavefunctions are combined with different approximate kernels. The lowest MAEs obtained in this manner for valence excitations are about 0.15-0.2 eV, which appears to be the practical limit of the accuracy of TDDFT that can be achieved by improving the Kohn-Sham potentials alone. These findings are consistent with previous reports on the benefits of accurate exchange-correlation potentials in TDDFT, but provide new insights and afford more definitive conclusions.
Collapse
Affiliation(s)
- Jaspreet Kaur
- Department of Chemistry , The University of Western Ontario , London , Ontario N6A 5B7 , Canada
| | - Egor Ospadov
- Department of Chemistry , The University of Western Ontario , London , Ontario N6A 5B7 , Canada
| | - Viktor N Staroverov
- Department of Chemistry , The University of Western Ontario , London , Ontario N6A 5B7 , Canada
| |
Collapse
|
49
|
Suellen C, Freitas RG, Loos PF, Jacquemin D. Cross-Comparisons between Experiment, TD-DFT, CC, and ADC for Transition Energies. J Chem Theory Comput 2019; 15:4581-4590. [DOI: 10.1021/acs.jctc.9b00446] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Cinthia Suellen
- Departamento de Quimica, Laboratorio Computacional de Materiais, Universidade Federal de Mato Grosso, Cuiaba, Mato Grosso 78060, Brazil
| | - Renato Garcia Freitas
- Departamento de Quimica, Laboratorio Computacional de Materiais, Universidade Federal de Mato Grosso, Cuiaba, Mato Grosso 78060, Brazil
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Cedex 9 Toulouse, France
| | - Denis Jacquemin
- Laboratoire CEISAM—UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Cedex 3 Nantes, France
| |
Collapse
|
50
|
Affiliation(s)
- Pierre‐François Loos
- Laboratoire de Chimie et Physique QuantiquesUniversité de Toulouse, CNRS, UPS France
| | - Denis Jacquemin
- Laboratoire CEISAM – UMR CNRS 6230Université de Nantes 2 Rue de la Houssinière BP 92208, 44322 Nantes Cedex 3 France
| |
Collapse
|