Combining the Bethe–Salpeter Formalism with Time-Dependent DFT Excited-State Forces to Describe Optical Signatures: NBO Fluoroborates as Working Examples

Double Hybrid Density Functionals for the Electronic Excitation Energies of Linear Cyanines

The lowest bright electronic excitations of seven model linear cyanines (CN3-CN15) using 28 double-hybrid (DH) density functionals are benchmarked against accurate and recent CC3 results. Some of these DH functionals are recently designed specifically for excited electronic state calculations. In addition, CIS, CIS(D), SCS-CIS(D), and SOS-CIS(D) were also tested. Four different basis sets were used for the vertical electronic excitation calculations: cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis. Augmented basis sets (e.g. aug-cc-pVDZ and aug-cc-pVTZ) are found to be required for accurate and consistent results using DH functionals. The DH functionals tested in this work are classified into four main groups: global double-hybrids (GDH), range-separated double-hybrids (RSDH), spin-component and spin-opposite scaling global double-hybrids (SCS/SOS-GDH), and spin-component and spin-opposite scaling range-separated double-hybrids (SCS/SOS-RSDH). Within these groups, the SCS/SOS-RSDH group of functionals is found to provide the lowest mean absolute error (MAE) values (in the range 0.020-0.148 eV) in comparison with the GDH group (0.195-0.441 eV), the RSDH group (0.186-0.511 eV), and the SCS/SOS-GDH group (0.079-0.461 eV). Of all the DH functionals and ab initio methods investigated in the present contribution, the following functionals are found to be the most accurate and consistent: SCS-ωB2GPPLYP (MAE = 0.036 eV), SOS-ωB2GPPLYP (MAE = 0.020 eV), SOS-ωB88PP86 (MAE = 0.035 eV), and SOS-ωPBEPP86 (MAE = 0.037 eV). In general, the ab initio methods tested here show mediocre performance as compared to many DH functionals.

Accurate Vertical Excitation Energies of BODIPY/Aza-BODIPY Derivatives from Excited-State Mean-Field Calculations

We report a benchmark study of vertical excitation energies and oscillator strengths for the HOMO → LUMO transitions of 17 boron–dipyrromethene (BODIPY) structures, showing a large variety of ring sizes and substituents. Results obtained at the time-dependent density functional theory (TDDFT) and at the delta-self-consistent-field (ΔSCF) by using 13 different exchange correlation kernels (within LDA, GGA, hybrid, and range-separated approximations) are benchmarked against the experimental excitation energies when available. It is found that the time-independent ΔSCF DFT method, when used in combination with hybrid PBE0 and B3LYP functionals, largely outperforms TDDFT and can be quite competitive, in terms of accuracy, with computationally more costly wave function based methods such as CC2 and CASPT2.

Accurate predictions of the electronic excited states of BODIPY based dye sensitizers using spin-component-scaled double-hybrid functionals: a TD-DFT benchmark study

The vertical excitation energies of 13 BODIPY based dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals from different DFT rungs. Most TD-DFT results were found to overestimate the excitation energies, and show mean absolute error (MAE) values in the range 0.2–0.5 eV. The dispersion-corrected, spin-component-scaled, double-hybrid (DSD) functionals DSD-BLYP and DSD-PBEP86 were found to have the smallest MAE values of 0.083 eV and 0.106 eV, respectively, which is close to the range of average errors found in the more expensive coupled-cluster methods. Moreover, DSD-BLYP and DSD-PBEP86 functionals show excellent consistency and quality of results (standard deviation = 0.048 eV and 0.069 eV respectively). However, the range separated hybrid (RSH) and the range separated double hybrid (RSDH) functionals were found to provide the best predictability (linear determination coefficient R² > 0.97 eV).

The excitation energies of 13 BODIPY dye sensitizers are benchmarked by means of TD-DFT, using 36 functionals. Spin-component-scaled double-hybrid (DSD) functionals are found to show the best performance.

Calculated linear and nonlinear optical absorption spectra of phosphine-ligated gold clusters

Although prediction of optical excitations of ligated gold clusters by time-dependent density functional theory (TDDFT) is relatively well-established, limitations still exist, for example in the choice of the exchange-correlation functional....

Optical and Electronic Properties of Organic NIR-II Fluorophores by Time-Dependent Density Functional Theory and Many-Body Perturbation Theory: GW-BSE Approaches

Organic-molecule fluorophores with emission wavelengths in the second near-infrared window (NIR-II, 1000-1700 nm) have attracted substantial attention in the life sciences and in biomedical applications because of their excellent resolution and sensitivity. However, adequate theoretical levels to provide efficient and accurate estimations of the optical and electronic properties of organic NIR-II fluorophores are lacking. The standard approach for these calculations has been time-dependent density functional theory (TDDFT). However, the size and large excitonic energies of these compounds pose challenges with respect to computational cost and time. In this study, we used the GW approximation combined with the Bethe-Salpeter equation (GW-BSE) implemented in many-body perturbation theory approaches based on density functional theory. This method was used to perform calculations of the excited states of two NIR molecular fluorophores (BTC980 and BTC1070), going beyond TDDFT. In this study, the optical absorption spectra and frontier molecular orbitals of these compounds were compared using TDDFT and GW-BSE calculations. The GW-BSE estimates showed excellent agreement with previously reported experimental results.

Assessment of local coupled cluster methods for excited states of BODIPY/Aza-BODIPY families

It was previously reported that Laplace transformed local CC2 (LCC2*) provided the best agreement (MAE = 0.145 eV) when comparing vertical excitation energies to experimental λ_max for a benchmark set of 17 BODIPY/Aza-BODIPY molecules. However, these energies did not agree with values obtained from canonical CC2. Here we report LCC2* computations of vertical excitation energies on the same benchmark set of molecules using a newly implemented treatment of the ground state. Comparison with resolution-of-identity approximate coupled cluster to second-order (RI-CC2) results demonstrate that the new LCC2* results agree quantitatively. Furthermore, these values can easily be corrected empirically to also provide excellent agreement with the experiment. We show that the local algebraic diagrammatic construction to second-order (LADC(2)) method exhibits the same differences between implementations as seen for LCC2. The source of the difference is traced to an improved treatment of the ground state in the local methods, which decreases agreement with the experiment (as attributed to a fortuitous cancellation of errors) but significantly improves agreement with RI-CC2. While the absolute vertical excitation energies now show larger deviations, there remains a strong linear correlation between the LCC2* results and the experiment. For the 17 BODIPY/Aza-BODIPY molecules vertical excitation energies are determined using DLPNO-STEOM-CCSD and shown to have excellent agreement with experimental λ_max (MAE = 0.145 eV), which is the best of all the single-reference methods. The vertical excitation energies are determined using LCC2*, empirically corrected LCC2*, and RI-CC2 for a series of eight large BODIPYs and Aza-BODIPYs.

Density Functional Studies on Photophysical Properties of Boron-Pyridyl-Imino-Isoindoline Dyes: Effect of the Fusion

In this work, to make out the aryl-fusion effect on the photophysical properties of boron-pyridyl-imino-isoindoline dyes, compounds 1-5 were theoretically studied through analyses of their geometric and electronic structures, optical properties, transport abilities, and radiative (k r) and non-radiative decay rate (k nr) constants. The highest occupied molecular orbitals of aryl-fused compounds 2-5 are higher owing to the extended conjugation. Interestingly, aryl fusion in pyridyl increases the lowest unoccupied molecular orbital (LUMO) level, while isoindoline decreases the LUMO level; thus, 4 and 5 with aryl fusion both in pyridyl and isoindoline exhibit a similar LUMO to 1. Compounds 4 and 5 show relatively low ionization potentials and high electron affinities, suggesting a better ability to inject holes and electrons. Importantly, the aryl fusion is conducive to the decrease of k IC. The designed compound 5 exhibits a red-shifted emission maximum, low λh, and low k IC, which endow it with great potential for applications in organic electronics. Our investigation provides an in-depth understanding of the aryl-fusion effect on boron-pyridyl-imino-isoindoline dyes at molecular levels and demonstrates that it is achievable.

The Bethe–Salpeter equation in chemistry: relations with TD-DFT, applications and challenges

We review the Bethe–Salpeter formalism and analyze its performances for the calculation of the excited state properties of molecular systems.

Benchmark of Bethe-Salpeter for Triplet Excited-States

We have evaluated the accuracy of the Bethe-Salpeter singlet-triplet transition energies as well as singlet-triplet and triplet-triplet splittings for 20 organic molecules, using as reference the CC3 values determined by Thiel and co-workers with both the TZVP and aug-cc-pVTZ atomic basis sets. Our excitation energies are obtained on the basis of GW quasiparticle energy levels that are self-consistently converged with respect to the starting DFT eigenvalues. In its current form, BSE/GW is often unable to provide a balanced description of both singlet and triplet excited-states. While the singlet-singlet and triplet-triplet energy separations are obtained accurately, triplets are located too close in energy from the ground-state, by typically -0.55 eV when using standard functionals to generate the starting eigenstates. Applying the Tamm-Dancoff approximation upshifts the BSE triplet energies and allows reducing this error to ca. -0.40 eV, while using M06-HF eigenstates allows a further increase and hence a reduction of the error for triplet states, but at the cost of larger errors for the singlet excited-states. At this stage, the most accurate TD-DFT estimates therefore remain competitive for computing singlet-triplet transition energies. Indeed, with M06-2X, irrespective of the application or not of the Tamm-Dancoff approximation and of the selected atomic basis set, the deviations obtained with TD-DFT are rather small.

Bethe-Salpeter study of cationic dyes: Comparisons with ADC(2) and TD-DFT

We present a theoretical investigation of the excited-state properties of a large series of structurally diverse arylcarbonium derivatives that are known to be challenging for theoretical models. More specifically, we compare the pros and cons of TD-DFT (TD-M06-2X), ADC(2), and BSE/GW approaches for a large panel of compounds, using two different solvent models. Both 0-0 and vertical transition energies are considered and compared to the experimental values. All approaches reasonably reproduce the auxochromic and acidochromic shifts, although in most cases both TD-DFT and BSE/GW return larger correlation with experimental values than ADC(2) for these shifts. In contrast, the absolute transition energies obtained with ADC(2) tend to be closer to the measurements, TD-DFT using the M06-2X functional largely overestimating the experimental references (by ca. 0.5 eV), and BSE/GW providing intermediate values. In addition, we show that the selected solvent model has a significant impact on the results, the corrected linear-response approach providing larger transition energies than its linear-response counterpart.

Investigating the optical properties of BOIMPY dyes using ab initio tools

Using a computational approach combining Time-Dependent Density Functional Theory (TD-DFT) and second-order Coupled Cluster (CC2) approaches, we investigate the spectral properties of a large panel of BOIMPY dyes.

TD-DFT Benchmark on Inorganic Pt(II) and Ir(III) Complexes

We report in the present paper a comprehensive investigation of representative Pt(II) and Ir(III) complexes with special reference to their one-photon absorption spectra employing methods rooted in density functional theory and its time dependent extension. We have compared nine different functionals ranging from generalized gradient approximation (GGA) to global or range-separated hybrids, and two different basis sets, including pseudopotentials for 4 iridium and 7 platinum complexes. It turns out that hybrid functionals with the same exchange part give comparable results irrespective of the specific correlation functional (i.e., B3LYP is very close to B3PW91 and PBE0 is very close to MPW1PW91). More recent functionals, such as CAM-B3LYP and M06-2X, overestimate excitation energies, whereas local functionals (BP86 -GGA-, M06-L -Meta GGA-) strongly underestimate transition energies with respect to experimental results. As expected, basis set effects are weak, and the use of a triple-ζ polarized (def2-TZVP) basis set does not significantly improve the computed excitation energies with respect to a classical double-ζ basis set (LANL2DZ) augmented by polarization functions, but it significantly raises the computational effort.

Why do TD-DFT excitation energies of BODIPY/Aza-BODIPY families largely deviate from experiment? Answers from electron correlated and multireference methods

The vertical excitation energies of 17 boron-dipyrromethene (BODIPY) core structures with a variety of substituents and ring sizes are benchmarked using time-dependent density functional theory (TD-DFT) with nine different functionals combined with the cc-pVTZ basis set. When compared to experimental measurements, all functionals provide mean absolute errors (mean AEs) greater than 0.3 eV, larger than the 0.1-0.3 eV differences typically expected from TD-DFT. Due to the high linear correlation of TD-DFT results with experiment, most functionals can be used to predict excitation energies if corrected empirically. Using the CAM-B3LYP functional, 0-0 transition energies are determined, and while the absolute difference is improved (mean AE = 0.478 eV compared to 0.579 eV), the correlation diminishes substantially (R(2) = 0.961 to 0.862). Two very recently introduced charge transfer (CT) indices, q(CT) and d(CT), and electron density difference (EDD) plots demonstrate that CT does not play a significant role for most of the BODIPYs examined and, thus, cannot be the source of error in TD-DFT. To assess TD-DFT methods, vertical excitation energies are determined utilizing TD-HF, configuration interaction CIS and CIS(D), equation of motion EOM-CCSD, SAC-CI, and Laplace-transform based local coupled-cluster singles and approximate doubles LCC2* methods. Moreover, multireference CASSCF and CASPT2 vertical excitation energies were also obtained for all species (except CASPT2 was not feasible for the four largest systems). The SAC-CI/cc-pVDZ, LCC2*/cc-pVDZ, and CASPT2/cc-pVDZ approaches are shown to have the smallest mean AEs of 0.154, 0.109, and 0.100 eV, respectively; the utility of the LCC2* approach is demonstrated for eight extended BODIPYs and aza-BODIPYs. We found that the problems with TD-DFT arise from difficulties in dealing with the differential electron correlation (as assessed by comparing CCS, CC2, LR-CCSD, CCSDR(T), and CCSDR(3) vertical excitation energies for five compounds) and from contributions of multireference character and double excitations (from analysis of the CASSCF wave functions).

0-0 Energies Using Hybrid Schemes: Benchmarks of TD-DFT, CIS(D), ADC(2), CC2, and BSE/GW formalisms for 80 Real-Life Compounds

The 0-0 energies of 80 medium and large molecules have been computed with a large panel of theoretical formalisms. We have used an approach computationally tractable for large molecules, that is, the structural and vibrational parameters are obtained with TD-DFT, the solvent effects are accounted for with the PCM model, whereas the total and transition energies have been determined with TD-DFT and with five wave function approaches accounting for contributions from double excitations, namely, CIS(D), ADC(2), CC2, SCS-CC2, and SOS-CC2, as well as Green's function based BSE/GW approach. Atomic basis sets including diffuse functions have been systematically applied, and several variations of the PCM have been evaluated. Using solvent corrections obtained with corrected linear-response approach, we found that three schemes, namely, ADC(2), CC2, and BSE/GW allow one to reach a mean absolute deviation smaller than 0.15 eV compared to the measurements, the two former yielding slightly better correlation with experiments than the latter. CIS(D), SCS-CC2, and SOS-CC2 provide significantly larger deviations, though the latter approach delivers highly consistent transition energies. In addition, we show that (i) ADC(2) and CC2 values are extremely close to each other but for systems absorbing at low energies; (ii) the linear-response PCM scheme tends to overestimate solvation effects; and that (iii) the average impact of nonequilibrium correction on 0-0 energies is negligible.

Benchmarking the Bethe-Salpeter Formalism on a Standard Organic Molecular Set

We perform benchmark calculations of the Bethe-Salpeter vertical excitation energies for the set of 28 molecules constituting the well-known Thiel’s set, complemented by a series of small molecules representative of the dye chemistry field. We show that Bethe-Salpeter calculations based on a molecular orbital energy spectrum obtained with non-self-consistent G0W0 calculations starting from semilocal DFT functionals dramatically underestimate the transition energies. Starting from the popular PBE0 hybrid functional significantly improves the results even though this leads to an average -0.59 eV redshift compared to reference calculations for Thiel’s set. It is shown, however, that a simple self-consistent scheme at the GW level, with an update of the quasiparticle energies, not only leads to a much better agreement with reference values, but also significantly reduces the impact of the starting DFT functional. On average, the Bethe-Salpeter scheme based on self-consistent GW calculations comes close to the best time-dependent DFT calculations with the PBE0 functional with a 0.98 correlation coefficient and a 0.18 (0.25) eV mean absolute deviation compared to TD-PBE0 (theoretical best estimates) with a tendency to be red-shifted. We also observe that TD-DFT and the standard adiabatic Bethe-Salpeter implementation may differ significantly for states implying a large multiple excitation character.

Taking up the cyanine challenge with quantum tools

Cyanine derivatives, named from the Greek word kyanos meaning dark-blue, were discovered more than 150 years ago and remain one of the most widely used classes of organic dyes with contemporary applications in photography (panchromatic emulsions), information storage (CD-R and DVD-R media) and biochemistry (DNA and protein labeling) fields. Cyanine chromogens consist of a charged π-conjugated segment containing an odd number of sp(2) carbon atoms with the chain capped at the extremities by two electronegative centers, typically nitrogen or oxygen atoms. Cyanines are characterized by a vanishing bond length alternation indicating nearly equal carbon-carbon bond lengths, as well as a very intense and sharp absorption band presenting a shoulder. This hallmark band undergoes a strong red shift when the chain is extended. This so-called vinyl shift is extremely large (ca. 100 nm for each pair of carbon atoms added in the π-conjugated path), making cyanines ideal building blocks for the design of devices with near-infrared applications. Numerous cyanines also exhibit emission bands with large quantum yields. These exceptional optical properties explain why both canonical cyanines and the corresponding fluoroborates (e.g., boron-dipyrromethene, BODIPY) remain the focus of an ever-growing body of experimental work. In turn, this popularity has stimulated quantum mechanical investigations aiming, on the one hand, at probing the specific electronic nature of cyanine dyes and, on the other hand, at helping to design new dyes. However, the adiabatic approximation to time-dependent density functional theory, the most widespread ab initio model for electronically excited states, fails to accurately reproduce the absorption spectra of cyanine derivatives: it yields a systematic and large underestimation of the experimental wavelengths irrespective of the details of the computational protocol. In contrast, highly correlated wave function approaches provide accurate transition energies for model systems but are hardly applicable to real-life cyanines and BODIPY. This indicates that setting up a computationally tractable theoretical protocol that provides both robust and accurate optical spectra for cyanine-based dyes is a major challenge that has only been taken up lately. In this Account, we compile the most recent advances in the field by considering both compact streptocyanines and large fluoroborates. For the former, we summarize the key results obtained with a large panel of theoretical approaches, allowing us not only to understand the origin of the cyanine challenge but also to pinpoint the schemes presenting the most promising accuracy/effort ratio. For the latter, we show via selected examples how theoretical models can be used to reproduce simultaneously experimental band shapes and transition energies, thus paving the way to an efficient in silico design of new compounds.