# Controlling the fluorescence quantum yields of benzothiazole-difluoroborates by optimal substitution

Patryk Rybczyński, Manon H E Bousquet, Anna Kaczmarek-Kędziera, Beata

Jędrzejewska, Denis Jacquemin, Borys Ośmialowski

► To cite this version:

Patryk Rybczyński, Manon H E Bousquet, Anna Kaczmarek-Kędziera, Beata Jędrzejewska, Denis Jacquemin, et al. Controlling the fluorescence quantum yields of benzothiazole-difluoroborates by optimal substitution. Chemical Science, 2022, 13 (45), pp.13347 - 13360. 10.1039/d2sc05044g. hal-04945274

HAL Id: hal-04945274

https://hal.science/hal-04945274v1

Submitted on 13 Feb 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
EDGE ARTICLE

View Article Online

View Journal | View Issue

<figure><img src="image_2.png" /></figure>

Check for updates

Cite this: *Chem. Sci.*, 2022, 13, 13347

All publication charges for this article have been paid for by the Royal Society of Chemistry

# Controlling the fluorescence quantum yields of benzothiazole-difluoroborates by optimal substitution†

Patryk Rybczyński,‡ª Manon H. E. Bousquet,‡ᵇ Anna Kaczmarek-Kędzieria,ª Beata Jędrzejewska, ID c Denis Jacquemin ID *bd and Borys Ośmiałowski ID *a

Precise tuning of the fluorescence quantum yield, vital for countless applications of fluorophores, remains exceptionally challenging due to numerous factors affecting energy dissipation phenomena often leading to its counterintuitive behavior. In contrast to the absorption and emission wavelength which can be precisely shifted to the desired range by simple structural changes, no general strategy exists for controllable modification of the fluorescence quantum yield. The rigidification of the molecular skeleton is known to usually enhance the emission and can be practically realized via the limiting molecular vibrations by aggregation. However, the subtle balance between the abundant possible radiative and non-radiative decay pathways makes the final picture exceptionally sophisticated. In the present study, a series of nine fluorophores obtained by peripheral substitution with two relatively mild electron donating and electron withdrawing groups are reported. The obtained fluorescence quantum yields range from dark to ultra-bright and the extreme values are obtained for the isomeric molecules. These severe changes in emission efficiency have been shown to arise from the complex relationship between the Franck-Condon excited state and conical intersection position. The experimental findings are rationalized by the advanced quantum chemical calculations delivering good correlation between the measured emission parameters and theoretical radiative and internal conversion rate constants. Therefore, the described substituent exchange provides a method to rigorously adjust the properties of molecular probes structurally similar to thioflavin T.

Received 9th September 2022

Accepted 21st October 2022

DOI: 10.1039/d2sc05044g

rsc.li/chemical-science

## 1 Introduction

The fine-tuning of light emission of organic dyes is a subject of intensive research due to the rapid development of myriads of fluorescence- and phosphorescence-based technologies. Efforts are focused on optimizing the synthetic methods, improving the purification processes, and, of course, performing structural modifications allowing tailoring the dyes to meet the market

demands. Starting from a known emissive core, such an optimization can be achieved through various structural changes including the introduction of side substituent(s), the elongation of the π-conjugation path, atom exchange in the emitter, doubling the acceptor core, etc.¹⁻<fnref n="4" /> Among the available fluorescent cores, the BODIPY family, that contains a BF₂ group providing both improved rigidity and electron-withdrawing character, is clearly one of the most popular platforms.⁵,⁶ This stems from the high fluorescence quantum yield of the parent BODIPY molecule (close to 100%), and the many synthetic approaches that have been developed for modifying BODIPYs. It is therefore not surprising that BODIPY derivatives are now used in countless areas such as chemosensing,⁷,⁸ bioimaging,⁹,¹⁰ photodynamic therapy,¹¹,¹² and redox batteries.¹³ For BODIPYs, like for other dyes, it is generally rather straightforward to predict the impact of chemical modifications on the absorption and emission wavelengths. For instance, longer π-conjugation paths typically redshift the fluorescence,¹⁴,¹⁵ while the same may be realized for absorption and emission peaks by introducing an intramolecular charge-transfer (ICT) state caused by an electron acceptor (A) and electron donor (D) attached at two extremities of a dye. In contrast, intuitively predicting the variations of the fluorescence quantum yield (φ<sub>f</sub>) remains very

<figure><img src="image_3.png" /></figure>

<footnote>ªFaculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina Street 7, 87-100 Toruń, Poland. E-mail: borys.osmiałowski@umk.pl</footnote>

<footnote>ᵇNantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France. E-mail: Denis.Jacquemin@univ-nantes.fr</footnote>

<footnote>ᶜBydgoszcz University of Science and Technology, Faculty of Chemical Technology and Engineering, Seminaryjna 3, 85-326 Bydgoszcz, Poland</footnote>

<footnote n="4">⁴Institut Universitaire de France (IUF), Paris, FR-75005, France</footnote>

<footnote>†Electronic supplementary information (ESI) available: Fluorescence quantum yield in various media, fluorescence lifetime, experimental radiative and non-radiative rate constants, fluorescence spectra in glycerol-methanol mixtures, NMR and mass spectra, barycenters for charge transfer, Wiberg bond indexes, dipole moments in the ground and excited states, comparison of various vibronic scheme approaches, Cartesian coordinates of the S₀ and S₁ minima. See DOI: https://doi.org/10.1039/d2sc05044g</footnote>

<footnote>‡These authors contributed equally to this work.</footnote>
<figure><img src="image_4.png" /></figure>

challenging. For instance, it is known that more rigid molecules exhibit higher $\phi_f$ and that redshifting the emission, for instance through a D-A strategy, typically yields a smaller $\phi_f$ (the so-called “energy-gap law”), yet these rules-of-thumb remain qualitative, often fail, and are typically not very helpful for practical design. This is because $\phi_f$ depends on all non-radiative pathways and their balance,¹⁶ so that subtle chemical changes might open (or close) some specific routes or, more often, tilt the radiative/non-radiative balance in an unexpected direction. In this framework, it is interesting to stress that first-principles calculations face exactly the same challenges as chemical intuition: whilst emission energies and band shapes can typically be predicted with sufficient reliability for dye design,¹⁷ accurate *ab initio* predictions of $\phi_f$ are in their infancy.¹⁸⁻²⁶ Likewise, one might wish to turn towards machine learning approaches, but yet again, they suffer from the same bias: reliable predictions can be obtained for wavelengths, much less so for $\phi_f$.²⁷,²⁸

Among the ICT dyes, thioflavin T (TT) attracted attention as a dye used for investigating amyloid β (Aβ) aggregation²⁹⁻³¹ but its usage is not limited to amyloids.³² TT is fluorescent but with a low quantum yield in water,³³ while its brightness significantly improves after aggregation with amyloids, with a $\phi_f$ reaching ca. 45%.³⁴ The aggregation induced emission (AIE) mechanism³⁵,³⁶ responsible for this behavior is related to the restriction of molecular motions that inhibit the emission in solution.³⁷ Indeed, the non-radiative processes in TT are driven not only by solute-solvent interactions but also by the rotation around a single bond present in its structure,³⁸,³⁹ and this rotation is hampered in the aggregated state (see below). From the structural point of view the common feature of many molecules used in studies of Aβ and its aggregation is the presence of a benzothiazole core.⁴⁰⁻⁴² This applies to TT that contains a positively charged benzothiazole moiety coupled to a strong electron donating group, namely the N,N-dimethylamino moiety. In several other series of benzothiazole-based fluorophores a specific twisting mode controls (and often limits) the $\phi_f$.⁴³ In TT, and in many other compounds, the introduction of the strong electron-donating NMe₂ moiety advantageously maximizes the ICT and redshifts the emission bands, but unfortunately also reduces $\phi_f$ in most solvents,³³ in part due to the enhanced molecular motions. In contrast, weaker donors might not deliver the desired absorption and emission in the red part of the spectrum. In the same vein, it is also known that most compounds containing the strong electron-accepting NO₂ group fluoresce weakly,⁴⁴ besides few exceptions.⁴⁵ Therefore, the elimination of the strongest D/A groups from the fluorophore together with the use of multiple substitutions could *a priori* be a smart strategy to bathochromically shift the emission while maintaining a relatively large $\phi_f$. Yet again, such “deduction” remains of qualitative nature only.

The current investigation tackles the impact of double substitution (Fig. 1) on a structural core greatly inspired by TT. The BF₂ group was introduced to rigidify the skeleton, delivering a core structure combining the advantages of both the BODIPY core and TT. Still, one “twistable” single bond is

<figure><img src="image_5.png" /><figcaption>Fig. 1 Structure of the nine compounds investigated herein.</figcaption></figure>

present in these molecules, which might allow significant relaxation in the excited-state.

Our global target is to understand how the photophysical properties can be controlled by substitution with a rather mild donor (OMe, which is much weaker than NMe₂) and an induc-tive acceptor (CF₃, which is also much less potent than NO₂; CF₃ has typically a trifling impact on absorption and emission wavelengths). In more detail, our main goals are: (i) to determine to what extent extremely strong donors or acceptors are mandatory for tuning the photophysical properties of dyes; (ii) to check if the substitution pattern in isomeric molecules significantly impacts their emissive properties; and (iii) to test the molecules topologically similar to TT for their use in various environments. The OMe and CF₃ groups are introduced systematically in the two substitution points (R and R' in Fig. 1), namely the benzothiazole fragment (R) and the *para* position of the phenyl ring (R'). The simple exchange of the selected groups between R and R' in Fig. 1 provides a set of nine compounds exhibiting, as we detail below, non-intuitive patterns in terms of photophysical properties.

## 2 Results and discussion

### 2.1 Experimental

The synthesis of all derivatives is described in the Experimental details section. The key photophysical properties of all dyes measured in chloroform are collected in Table 1. The corresponding absorption (black, solid), emission (red) and the mirrored emission spectra (black, dashed) are shown in Fig. 2.

All compounds exhibit a main absorption band between 340 and 370 nm. It is noteworthy that the H-to-CF₃ substitution has always a very mild impact on both the absorption position and the band shape, whereas the use of the donating OMe group induces noticeable changes, as it tends to bathochromically move the absorption by *ca.* 20–30 nm (Fig. 2). Interestingly, dye 1, that has a DAD quadrupolar topology, presents the most redshifted absorption band, the traditional push-pull strategy represented by compounds 3 and 7 being slightly less effective. This confirms that the BF₂-substituted benzothiazole is a strongly electron-withdrawing moiety. Notably the intensity of the absorption, as given by $ε$, can be maximized by plugging a methoxy group on the phenyl side (at R') irrespective of the group substituting the benzothiazole. Nevertheless the $ε$ always remains in the 30–80 × 10³ M⁻¹ cm⁻¹ regime, values that are
<table><caption>**Table 1** Photophysical properties of dyes **1–9** in CHCl<sub>3</sub>: longest wavelength of maximal absorption (nm), molar absorption coefficient (M<sup>-1</sup> cm<sup>-1</sup>), emission wavelength (nm), Stokes shift (cm<sup>-1</sup>), fluorescence quantum yield (%), full width at half maximum of the absorption and emission bands (cm<sup>-1</sup>), and dye topology. Radiative and non-radiative rates are provided in the ESI</caption><thead><tr><th>Compound (R/R')</th><th>λ<sup>abs</sup><sub>max</sub></th><th>ε</th><th>λ<sup>em</sup><sub>max</sub></th><th>Δ<sub>SS</sub></th><th>φ<sub>f</sub></th><th>FWHM<sub>abs</sub></th><th>FWHM<sub>em</sub></th><th>Topology</th></tr></thead><tbody><tr><td>1 (OMe/OMe)</td><td>368.0</td><td>76 900</td><td>454.5</td><td>5172</td><td>13.1</td><td>4246</td><td>4519</td><td>DAD</td></tr><tr><td>2 (OMe/H)</td><td>361.0</td><td>28 700</td><td>458.5</td><td>5891</td><td>0.9</td><td>4691</td><td>4586</td><td>DA-</td></tr><tr><td>3 (OMe/CF<sub>3</sub>)</td><td>367.0</td><td>26 000</td><td>491.5</td><td>6902</td><td>0.4</td><td>5046</td><td>4229</td><td>DAA</td></tr><tr><td>4 (H/OMe)</td><td>364.0</td><td>46 700</td><td>419.5</td><td>3635</td><td>75.6</td><td>3722</td><td>4253</td><td>-AD</td></tr><tr><td>5 (H/H)</td><td>343.0</td><td>32 100</td><td>439.5</td><td>6401</td><td>6.9</td><td>4269</td><td>5245</td><td>-A-</td></tr><tr><td>6 (H/CF<sub>3</sub>)</td><td>345.0</td><td>27 100</td><td>447.0</td><td>6614</td><td>0.4</td><td>4602</td><td>5130</td><td>-AA</td></tr><tr><td>7 (CF<sub>3</sub>/OMe)</td><td>362.5</td><td>63 500</td><td>412.0</td><td>3314</td><td>98.0</td><td>3465</td><td>3622</td><td>AAD</td></tr><tr><td>8 (CF<sub>3</sub>/H)</td><td>339.5</td><td>43 700</td><td>406.0</td><td>4825</td><td>28.0</td><td>4159</td><td>4972</td><td>AA-</td></tr><tr><td>9 (CF<sub>3</sub>/CF<sub>3</sub>)</td><td>341.0</td><td>30 300</td><td>434.5</td><td>6311</td><td>4.7</td><td>4361</td><td>5347</td><td>AAA</td></tr></tbody></table>

typical of conjugated organic dyes. The substitution of OMe by H and then by CF<sub>3</sub> in the sub-series (1–3, 4–6, 7–9, see Table 1) always delivers broader absorption bands.

Let us now turn to emission. While, in general, $\lambda_{\max}^{em}$ values are more variable across the series than their absorption counterparts, the changes in fluorescence spectra remain rather mild for the nine compounds studied herein (Fig. 2). As for the absorption, the influence of the CF<sub>3</sub> group is weak, whereas OMe has a more clearly detectable (bathochromic) effect, especially when added to the benzothiazole side (R). The presence of the methoxy also partially washes out the vibronic progression. It is 3 that delivers the largest $\lambda_{\max}^{em}$. Nevertheless, one should be cautious in interpreting directly the $\lambda_{\max}^{em}$ listed in Table 1. Indeed, as can be seen in Fig. 2, the relative heights of the various vibronic peaks may differ from one compound to another, indicating that the $\lambda_{\max}^{em}$, related to the most intense vibronic transition, could significantly change, while the actual energy shift between the electronic states (the bands) is mild. In any case, Fig. 2 indicates that the emission and absorption bands are close to mirror-shape hinting that the geometry of the emissive state should not (too) strongly differ from its absorption counterpart. We note that this mirror-shape trend is especially close to perfection for 7 that also displays a notably smaller Stokes shift than the other derivatives. On the other hand, the largest value of the Stokes shift is noted for 3, the only molecule that shows a broader absorption spectrum than the emission. For all dyes, and for both absorption and emission, the full width at half maximum (FWHM) values are in the 3500–5500 cm<sup>-1</sup> range, which is quite typical of such systems.

Until now, we have not unravelled any specific trends departing from “chemical expectations”, as one could clinically summarize the above findings by stating that all nine dyes share similar absorption and emission features but for slight/significant bathochromic/hyperchromic shifts induced by the methoxy donor. The picture is however much more blurred when investigating the fluorescence quantum yields listed in Table 1. Indeed $\phi_f$ is almost 100% for one of the push-pull dyes, 7, but close to zero for 3, presenting inverted R/R' groups. From the data of Table 1, one can deduce that adding a methoxy substituent at R' is beneficial in terms of brightness, irrespective of the group at R, whereas, in contrast, CF<sub>3</sub> group at R position increases the $\phi_f$. It is however not easy to interpret such a trend.

A first explanation could be that the methoxy group favours a more quinoidal form on the side phenyl group in which the single bond separating the phenyl and benzothiazole moieties gains a double character, what, in turn decreases deleterious excited-state motions. Such changes of bonding pattern would however need to be quite drastic to explain the extreme changes of $\phi_f$ noted in the series. In such a scenario, one would also expect the FWHMs, which are roughly related to structural differences between the ground and excited states to significantly change, an unseen outcome. Indeed, the linear determination coefficient ($R^2$) between $\phi_f$ and FWHM<sub>em</sub> is 0.50 only. Likewise, the “energy-gap law” is not very effective here, since the $R^2$ obtained by a $\phi_f$ to $\lambda_{\max}^{em}$ relationship attains 0.44. The only experimental parameter nicely correlating with $\phi_f$ seems to be the Stokes shift ($R^2 = 0.88$), but, as stated above, one should consider the $\lambda_{\max}^{em}$ with caution.

To obtain more insights into the relationship between the $\phi_f$ and molecular motions, changing the viscosity of the solution is a rather straightforward experimental approach. We therefore report in Fig. 3 and Table 2 a comparison between the $\phi_f$ measured in methanol and glycerol. We have used mixtures of methanol and glycerol to obtain media of controllable viscosity. This pair of solvents was selected as both are able to form hydrogen bonds due to their OH groups, allowing a smooth evolution of the solvents' physicochemical signatures. It is easy to notice that an increased viscosity positively influences the emission properties. Indeed as can be seen in Fig. S2† $\phi_f$ increases with the growing viscosity, giving good linear correlation coefficients for the doubly logarithmic plot of $\phi_f = f(\eta)$ except for 7, which exhibits practically a constant (close to 100%) fluorescence quantum yield irrespective of the medium. The observed relationships are in agreement with the ones of known viscosity-sensitive probes, in which the slope of the linear function is defined as a response factor of the probe.⁴⁶ We therefore underline that, for the current series, the largest slopes are logically obtained for the weakest emitters in non-viscous solvents, i.e., 2 and 6. In contrast, as stated above, 7 enjoys a fluorescence dependent neither upon the modification of the hydrogen bonding nature of the solvent (CHCl<sub>3</sub> to MeOH) nor upon the medium viscosity (MeOH vs. glycerol). This clearly suggests that several non-radiative mechanisms are at play in the series.
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 2** Normalized, corrected absorption (solid black line), fluorescence (solid red line) and mirror-image of fluorescence spectra (dashed black line) shifted to match the 0–0 point in CHCl<sub>3</sub>.</figcaption></figure>

From the points of view of the photophysical properties of
the dyes, we note that φ<sub>f</sub> is very similar in chloroform (non-
viscous, weakly hydrogen-bonding solvent) and glycerol (very
viscous and strong hydrogen bonding solvent) for 4, 7, and 8. In

contrast, a very low φ<sub>f</sub> is recorded for 2, 3, and 6 in both chlo-
roform and methanol but a reasonable increase is noticed in
glycerol for these systems. The structural features for those
compounds clearly differ from those of 4, 7, and 8, with a donor
<figure><img src="image_8.png" /><figcaption>Fig. 3 Heatmaps of the fluorescence quantum yields (%) measured in CHCl₃ (top left), methanol (top right), glycerol (bottom left). The bottom right figure provides the ratio between the glycerol and methanol $\phi_f$.</figcaption></figure>

<table><caption>Table 2 Key emissive properties of 1–9 in MeOH (m), glycerol (g) and response factor (last column, slope), see the caption of Table 1 and text for more details</caption><thead><tr><th>Compound (R/R')</th><th>φ<sub>f</sub> (m)</th><th>φ<sub>f</sub>(g)/φ<sub>f</sub>(m)</th><th>FWHM<sub>em</sub> (m)</th><th>FWHM<sub>em</sub> (g)</th><th>Log(φ<sub>f</sub>) = f(log(η)) (slope)</th></tr></thead><tbody><tr><td>1 (OMe/OMe)</td><td>2.4</td><td>11.8</td><td>5320</td><td>5022</td><td>0.350</td></tr><tr><td>2 (OMe/H)</td><td>0.4</td><td>27.1</td><td>5732</td><td>5202</td><td>0.441</td></tr><tr><td>3 (OMe/CF<sub>3</sub>)</td><td>0.8</td><td>13.0</td><td>5394</td><td>5149</td><td>0.376</td></tr><tr><td>4 (H/OMe)</td><td>26.5</td><td>2.7</td><td>5093</td><td>4570</td><td>0.148</td></tr><tr><td>5 (H/H)</td><td>2.0</td><td>9.2</td><td>6735</td><td>5502</td><td>0.350</td></tr><tr><td>6 (H/CF<sub>3</sub>)</td><td>0.6</td><td>23.3</td><td>6391</td><td>5559</td><td>0.463</td></tr><tr><td>7 (CF<sub>3</sub>/OMe)</td><td>99.8</td><td>1.0</td><td>3585</td><td>3445</td><td>0.000</td></tr><tr><td>8 (CF<sub>3</sub>/H)</td><td>4.0</td><td>9.0</td><td>6303</td><td>5364</td><td>0.306</td></tr><tr><td>9 (CF<sub>3</sub>/CF<sub>3</sub>)</td><td>1.6</td><td>15.3</td><td>6631</td><td>5533</td><td>0.392</td></tr></tbody></table>

on the benzothiazole and an electron acceptor on the phenyl or both. This means that reversing the substitution of the acceptor and donor allows going from normal compounds to structures behaving like molecular rotors sensitive to viscosity. Finally, the remaining 1, 5, and 9 are fluorescent in chloroform, very weakly emissive in methanol and their emission is much higher in glycerol. They are simultaneously sensitive to the hydrogen bonding/polar environment and viscosity (see Fig. S1†).

## 2.2 Theory

We have used ab initio calculations to probe the excited state (ES) features of all nine compounds. We have decided to

consider CHCl₃ as the medium for these simulations as the theoretical protocol relies on a continuum model for the environmental effects (see Experimental details section for technical specification). The careful inspection of the geometry of the investigated systems in the two electronic states confirms that the excitation promotes the *p*-quinoid form of the terminal phenyl ring. This can be deduced from the respective shortcomings and elongations of the alternate C–C bonds to form the =C(–CH=CH–)₂C= moiety. These changes are similar to those found previously in NMe₂-substituted benzothiazole derivatives.⁴⁷ Since one of the envisaged deactivation pathways involves the rotation of the phenyl ring, (TD-)DFT calculations
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>Fig. 4 (a) Barriers for rotation around the C15–C17 bond computed at the PCM(CHCl₃)-(TD-)MN15/6-311++G(d,p)) level for the ground and excited state geometries, (b) corresponding heatmaps.</figcaption></figure>

have been performed to quantify the related rotation barrier. Full DFT geometry optimizations of the transition state for the rotation around the C–C bond (C15–C17 in Fig. S52†) provide a ground-state picture consistent with chemical intuition. Indeed, as can be seen in Fig. 4, increasing the electron-donating character of the substituent at R' causes increase of the rotation barrier by *ca.* 2 kcal mol⁻¹, irrespective of the group at R, which fits with the quinoidal picture leading to a strengthening of the inter-ring bond upon addition of a mesomeric donor at R'. We also note that increasing the electron-accepting character of the substituent at the benzo-thiazole moiety induces a slight increase of the barrier, but the effect (*ca.* 0.5 kcal mol⁻¹) is insignificant for a DFT calculation. The largest ground-state barrier, 10.4 kcal mol⁻¹, is observed for 7, the lowest one for 3, 7.68 kcal mol⁻¹. If such a pattern seems to fit the observed *φ*<sub>*f*</sub> ranking (larger barriers, larger emission yields), one has to remember that these values are valid for the ground electronic state, *i.e.*, they are not representative of the structures of the emissive species.

In contrast, the ES picture does not follow a clear systematic order nor correlates with the experimental φ<sub>F</sub> (see Fig. 4). Interestingly, the ES rotation barriers are significantly higher than in the ground state, 8-out-of-9 exceeding 17 kcal mol<sup>-1</sup>, 7 being the only compound presenting a rather low ES barrier

(13.4 kcal mol⁻¹). If this smaller change of bond strength between the ground and excited states for 7 is consistent with the trends in the measured Stokes shift and almost perfect mirror symmetry of absorption and emission spectra, it is clearly not intuitive that the dye showing, by far, the lowest excited-state rotation barrier also displays the largest φ<sub>F</sub>. We note that the increase of barriers after photon absorption is consistent with chemical intuition (π-π* excitation decreasing the bond length alternation) and also with the Wiberg bond indexes being larger by 0.1 to 0.2 in the excited than in the ground state, a statement holding for all molecules (see Table S4†).

Let us now turn towards the charge transfer properties, *i.e.*, the transferred charge, *q*<sub>CT</sub>, the distance between the barycenters of density gain and depletion, *d*<sub>CT</sub>, and the variation of the dipole moment between the states, Δμ. Table 3 reports these parameters together with density difference plots. The gain of density is always localized on the fluoroborate ring, whereas the loss of density varies from one compound to another, the role of the OMe group being clear. This is confirmed by the positions of the barycenters (see the ESI†): R<sub>+</sub> remains almost unaffected by the variation of the substituents whereas R<sub>-</sub> depends on the considered dye. Among the analyzed fluorophores, the strongest photoinduced charge transfer appears for 3 which
<table><caption>**Table 3** Electron density difference (EDD, 0.002 au contour) plots and charge-transfer parameters⁴⁸ determined at the PCM(CHCl₃)-TD-MN15/6-311++G(d,p) level on the ground state optimized geometry. Blue (red) surface marks the electron density depletion (gain) upon excitation. *d*<sub>CT</sub> is given in Å, *q*<sub>CT</sub> in a.u., Δμ = μ<sub>ES</sub> – μ<sub>GS</sub> and dipole moments in D. The dipole moments of the optimized geometry of the excited state can be found in the ESI</caption><thead><tr><th>Compound (R/R<sup>1</sup>)</th><th>EDD</th><th>d<sub>CT</sub></th><th>q<sub>CT</sub></th><th>μ<sub>GS</sub></th><th>μ<sub>ES</sub></th><th>Δμ</th></tr></thead><tbody><tr><td>1 (OMe/OMe)</td><td></td><td>0.880</td><td>0.575</td><td>4.881</td><td>5.411</td><td>0.530</td></tr><tr><td>2 (OMe/H)</td><td></td><td>2.547</td><td>0.606</td><td>6.605</td><td>11.036</td><td>4.432</td></tr><tr><td>3 (OMe/CF<sub>3</sub>)</td><td></td><td>3.271</td><td>0.651</td><td>8.376</td><td>17.181</td><td>8.806</td></tr><tr><td>4 (H/OMe)</td><td></td><td>0.904</td><td>0.545</td><td>6.729</td><td>7.816</td><td>1.088</td></tr><tr><td>5 (H/H)</td><td></td><td>1.147</td><td>0.514</td><td>4.887</td><td>6.147</td><td>1.259</td></tr><tr><td>6 (H/CF<sub>3</sub>)</td><td></td><td>2.119</td><td>0.540</td><td>6.644</td><td>11.427</td><td>4.783</td></tr><tr><td>7 (CF<sub>3</sub>/OMe)</td><td></td><td>1.658</td><td>0.551</td><td>6.427</td><td>10.609</td><td>4.181</td></tr><tr><td>8 (CF<sub>3</sub>/H)</td><td></td><td>0.412</td><td>0.495</td><td>5.736</td><td>5.501</td><td>-0.235</td></tr><tr><td>9 (CF<sub>3</sub>/CF<sub>3</sub>)</td><td></td><td>1.479</td><td>0.503</td><td>4.056</td><td>6.091</td><td>2.036</td></tr></tbody></table>

undergoes an 8.7 D dipole increase in the ES, corresponding to a displacement of 0.65e over 3.2 Å, quite large values in Le Bahers' model.⁴⁹ There is a weak positive correlation between *d*<sub>CT</sub> and λ<sub>max</sub><sup>em</sup> with an R<sup>2</sup> of 0.61. Nevertheless, as for the barriers above, the CT parameters do not provide reasonable insights

into the emission yields. Indeed, the R<sup>2</sup> between the measured φ<sub>F</sub> and *d*<sub>CT</sub>, *q*<sub>CT</sub>, and Δμ are all tiny: 0.12, 0.03, and 0.04, respectively. In other words, the magnitude of the charge-transfer does not control the fluorescence efficiency.

<figure><img src="image_11.png" /></figure>
<table><caption>**Table 4** Singlet-triplet gap computed on the relaxed S<sub>1</sub> geometry (eV), with the corresponding SOC (cm<sup>-1</sup>). Energy difference between the most accessible MECP and the FC point calculated at the PCM(CHCl<sub>3</sub>)-SF-TD-MN15/6-311+G(d,p) level (eV). Radiative and non-radiative rates as obtained by theory (10<sup>8</sup> s<sup>-1</sup>). Theoretical fluorescence quantum yields (%)</caption><thead><tr><th>Compound (R/R<sup>1</sup>)</th><th>ΔE<sup>S<sub>1</sub>-T<sub>1</sub></sup></th><th>SOC</th><th>ΔE<sup>MECP-S<sub>1</sub>FC</sup></th><th>k<sub>r</sub></th><th>k<sub>ic</sub></th><th>k<sub>MECP</sub></th><th>φ<sup>no-MECPa</sup><sub>f</sub></th><th>φ<sup>MECPb</sup><sub>f</sub></th></tr></thead><tbody><tr><td>1 (OMe/OMe)</td><td>0.46</td><td>0.04</td><td>0.06</td><td>3.09</td><td>0.83</td><td>5060</td><td>79</td><td>0.0</td></tr><tr><td>2 (OMe/H)</td><td>0.66</td><td>0.03</td><td>-0.14</td><td>2.18</td><td>7.87</td><td>∞</td><td>22</td><td>0.0</td></tr><tr><td>3 (OMe/CF<sub>3</sub>)</td><td>0.73</td><td>0.04</td><td>-0.27</td><td>2.02</td><td>2.16</td><td>∞</td><td>48</td><td>0.0</td></tr><tr><td>4 (H/OMe)</td><td>0.28</td><td>0.03</td><td>0.47</td><td>4.76</td><td>0.36</td><td>0.00</td><td>93</td><td>93</td></tr><tr><td>5 (H/H)</td><td>0.44</td><td>0.03</td><td>0.12</td><td>3.21</td><td>0.57</td><td>682</td><td>85</td><td>0.1</td></tr><tr><td>6 (H/CF<sub>3</sub>)</td><td>0.52</td><td>0.03</td><td>-0.07</td><td>2.70</td><td>1.50</td><td>∞</td><td>64</td><td>0.0</td></tr><tr><td>7 (CF<sub>3</sub>/OMe)</td><td>0.24</td><td>0.04</td><td>0.86</td><td>5.61</td><td>0.26</td><td>0.00</td><td>96</td><td>96</td></tr><tr><td>8 (CF<sub>3</sub>/H)</td><td>0.36</td><td>0.04</td><td>0.56</td><td>3.92</td><td>0.28</td><td>0.00</td><td>93</td><td>93</td></tr><tr><td>9 (CF<sub>3</sub>/CF<sub>3</sub>)</td><td>0.44</td><td>0.04</td><td>0.10</td><td>2.78</td><td>2.14</td><td>1250</td><td>56</td><td>0.0</td></tr></tbody></table>

<sup>a</sup> Computed considering only the radiative and internal conversion rates. <sup>b</sup> Considering the radiative, internal and MECP pathways. <sup>c</sup> SOC calculation did not converge.

Given the above rather unfruitful attempts to rationalize the measured $\phi_f$, we decided to search for the most accessible pathways allowing to directly go back from the excited to the ground electronic state through non-radiative decay. First, given the presence of a sulfur atom, we evaluated the possibility of intersystem crossing (ISC), by evaluating the singlet-triplet gap and the spin-orbit coupling (SOC) at the relaxed excited state geometry. The results are listed in Table 4. Only one triplet is below S<sub>1</sub>, and as can be seen in that table the gaps are typically quite large (at least 0.24 eV) while the SOC matrix elements are truly trifling (<0.05 cm<sup>-1</sup>). This strongly hints that ISC is not a major pathway for the considered systems.

Next, we used spin-flip (SF) TD-DFT to locate the most accessible minimum energy crossing point (MECP) between the ground and excited states for all dyes. Indeed, several previous studies have demonstrated a direct relationship between the measured emission yield and the relative energy of the lowest MECP.<sup>18,21,50,51</sup> Using the procedure described in the Experimental section, we found that this most accessible MECP corresponds, in all cases, to a strongly deformed geometry with an almost 90° kink between the two parts of the molecule (see Fig. 5). The relative SF-TD-DFT energies of this MECP and the S<sub>1</sub> FC point are given in Table 4. Despite several attempts, we could not locate a proper transition state corresponding to the path to this MECP, likely indicating that the path is smooth and that these relative energies are indicators

of the accessibility of this MECP. Interestingly, we note that the three dyes having an MECP more stable than the FC point are 2, 3, and 6, that is, the three experimentally “dark” compounds, i.e., compounds showing $\phi_f$ smaller than 1% (see Table 1). This is consistent with a downhill path to the MECP. Likewise, the three compounds for which reaching the MECP implies an energetic penalty larger than 0.3 eV, namely, 4, 8, and 7 are the only three showing rather large (>25%) $\phi_f$ in the experiment. We note that this is the first clear correlation between theoretical estimates and experimental $\phi_f$ found in this work. It is worth noting that the $R^2 = 0.81$ for the correlation between these energies and the measured $\phi_f$ (in glycerol); a very similar value ($R^2 = 0.78$) is obtained in chloroform which was used in measurements and calculations. The slightly higher value for glycerol may be due to the higher rigidity of the molecules in that solvent, a feature of well-defined geometry as in calculations limiting additional dynamic effects associated with energy dissipation.

To go further, we have used Peng and Shuai's TVCF formalism to model both the radiative and internal conversion rates of all dyes on an equal footing (again see the Experimental section for details).<sup>52,53</sup> Such calculations involve the calculations of vibronic couplings on the basis of harmonic vibrational frequencies, as well as the selection of models for the coupling between the potential energy surfaces, the coordinates, etc. After some tests, we have chosen the most refined model (so-called

<figure><img src="image_12.png" /><figcaption>**Fig. 5** Two views of the MECP between the ground and the lowest singlet excited states as obtained at the PCM(CHCl<sub>3</sub>)-SF-TD-MN15/6-311+G(d,p) level of theory for 9. Similar MECPs are found for all dyes.</figcaption></figure>

<figure><img src="image_13.png" /></figure>
adiabatic Hessian), in a time-dependent framework and selected Cartesian coordinates – the computed rates are listed in Table 4. In Fig. 6, one can find comparisons between the theoretical and experimental emission lineshapes for the four

brightest derivatives and the agreement gives confidence that theory correctly describes the vibronic couplings for the present class of derivatives. The *ab initio* radiative rates are in the $10^8 \text{ s}^{-1}$ range, which fits the experimental trends (see the ESI†), the

<figure><img src="image_14.png" /><figcaption>Fig. 6 Left: computed (colored lines) and measured (black lines and grey background) fluorescence spectra for the four brightest dyes. Right: plots showing the dependence between the internal conversion rates and the energies for the same four systems, the vertical line gives the computed adiabatic energy.</figcaption></figure>
largest value being computed for 7, 5.61 × 10⁸ s⁻¹, a result perfectly consistent with the measurement, 5.93 × 10⁸ s⁻¹. The R² between experiment and theory for the k_f attains a very satisfying 0.87 value. We note nevertheless that the absolute magnitudes of the theoretical radiative rates significantly overestimate their experimental counterparts for the poorly emissive dyes, but for these almost dark compounds the experimental error bar becomes (much) larger, so it is not so straightforward to pinpoint if theory or experiment is faulty in that case. The internal conversion rates, k_ic, measure the non-radiative deactivation through non-adiabatic couplings between the excited and ground states. The computed values are in the 2 × 10⁷ to 8 × 10⁸ s⁻¹ range (see Table 4), i.e., cover a quite broad panel of values, yet these k_ic have in general no direct experimental counterparts. In addition, we wish to stress the inherent limits of theory in evaluation of this rate, as it is rather common that k_ic strongly depends on the selected broadening, in contrast to k_f.²² In the ESI Tables S6 and S7,† we provide comparisons between IC rates computed with different vibronic models, and one notices that the IC rate varies significantly depending on the selected model, but those changes do not impact the relatively poor correlation with experimental data.

By computing the φ_f on the basis of radiative and internal conversion rates only, one can estimate the emission yield if no other non-radiative processes take place, meaning, in the present case, that the MECPs would always be inaccessible. In this way, too large φ_f values are obtained but one notices that a value of 96%, consistent with the measurements, is determined for 7 which also presents the least accessible MECP (see above). Finally, it is possible to compute the non-radiative rate constant related to the MECP simply using

$$k_{\text{MECP}} = \frac{k_b T}{h} \exp\left(-\frac{\Delta E^{\text{MECP-S}_1^{\text{FC}}}}{RT}\right) \quad (1)$$

in an approach equivalent to the one used previously for BODIPYs.²¹ As can be seen in Table 4 such a more complete approach would split the nine dyes into two groups: (highly) emissive and non-fluorescent, and these groups are close to the experimental ones. In fact, taking the experimental φ_f as the reference, the errors are smaller than 0.1, but for three compounds, namely: **1** which is incorrectly predicted to be totally dark whereas the experimental yield is 13%, **4** for which theory overshoots the φ_f (93 versus 76%), and **8** which is the worst failure to our views, as theory grossly overestimates the yield (93 versus 28%). These deviations are likely significant revealing the inherent limitations of the theoretical model (both electronic structure and vibronic calculations). Alternatively, one could indeed compute k_MECP using as reference the energy of the relaxed S₁ geometry instead of the FC point. This has been investigated and the results are displayed in Fig. S60 in the ESI† showing that the correlation with the experimental quantum yields is not improved in this case. For compound **1**, computing the MECP rate from the FC point leads to k_MECP = 5.06 × 10¹¹ s⁻¹ and hence a negligible fluorescence yield, whereas from the relaxed S₁ one obtains k_MECP = 3.24 × 10⁵ s⁻¹ and a very large fluorescence yield. Obviously, the experimental truth stands in

between, and this illustrates the limits of static approaches: one can estimate high and low limit values for the MECP rate, whereas a dynamical approach would likely be needed to obtain a more balanced description.

## 3 Conclusions and outlook

The series of benzothiazoles considered here shows that the fluorescence quantum yield may be tuned systematically and strongly by exchanging two (relatively mild) groups located on the periphery of the structures only. From the intuitive point of view this seems to be realised by modifying the height of the rotation barrier around the only single bond in the structure, as a viscous environment is beneficial for the brightness of the synthesized dyes, indicating the presence of significant detrimental motions. However, theoretical calculations indicate that these motions correspond to the kinking of the structure around the single bond, rather than a rotation. Indeed, the difference of excited state energy between the FC point and the minimum energy crossing point between ground and excited states strongly correlates with experimental fluorescence quantum yield. This correlation is even higher in viscous environments. We conclude that simple substituent exchange by no means induces an intuitive effect though such an approach delivers significant changes in the emission. The current series undoubtedly shows that substituents being not particularly potent ones in light of their electronic effects may be applied to cover the full range of emission quantum yields. This feature may be further adapted in the design of rotor-based/conical intersection-tuned emitters capable of sensing the viscosity and aggregation, and thus in molecular probes for a variety of needs including biochemistry-oriented ones.

## 4 Experimental details

### 4.1 Synthesis and structure confirmation

All substrates for the synthesis (amines and benzoyl chlorides) and solvents were obtained from commercial sources. The studied series of compounds was obtained by a two-step synthesis. The amide: to a solution of 2-aminobenzothiazole (1 eq.) and triethylamine (2 eq.) in dry THF (20 ml, -78 °C, under an inert gas) the respective benzoyl chloride (1 eq.) was added dropwise (in dry THF). The mixture was allowed to warm to room temperature and stirred overnight. The solvent was evaporated and the residue extracted using DCM/NaHCO₃ water solution. After evaporation of the organic layer the compound was purified by column chromatography (eluent hexane/ethyl acetate). In the second step, the appropriate amide was complexed with BF₃·OEt₂. The etherate was added dropwise to a mixture of amide (1 eq.) and diisopropylethylamine (3 eq.) in DCM (20 ml, room temperature). The reaction mixture was stirred overnight and then the DCM was evaporated. The compound was purified by column chromatography (SiO₂, eluent DCM).

The structure of the synthesized compounds was confirmed by NMR and MS techniques. All of the NMR spectra were recorded at 400 MHz or 700 MHz on a Bruker spectrometer at
room temperature in CDCl₃. HRMS data were obtained with a Waters spectrometer in the EI mode. The melting points were measured using a Stuart SMP50 Digital Melting Point Instrument.

**4.1.1 Compound 1.** ((Z)-[(Difluoroboryloxy)(p-methoxyphenyl)methylene](6-methoxy-1,3-benzothiazol-2-yl)amine, OMe/OMe) yield 56.7%. Mp 215–217 °C, green powder ¹H NMR (700 MHz, from TMS, CDCl₃): δ (ppm) 8.36 (d, 2H, *J* = 9.03 Hz), 7.95 (d, 1H, *J* = 9.03 Hz), 7.26 (d, 1H, *J* = 2.52 Hz), 7.18 (dd, 1H, *J* = 9.10, 2.52 Hz), 7.00 (d, 2H, *J* = 8.96 Hz), 3.93 (s, 3H), 3.91 (s, 3H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 0.576 (t). ¹³C (100 MHz, from TMS, DMSO-d₆): δ (ppm): 172.3, 167.4, 164.8, 134.1, 132.7, 128.4, 123.3, 119.2, 116.8, 114.1, 105.3, 56.0, 55.6. ¹⁹F (376 MHz, from CFCl₃, DMSO-d₆): δ (ppm) –137.6. EI HRMS (*m/z*) [*M*]⁺ 362.0703 cal. 362.0708.

**4.1.2 Compound 2.** ((Z)-[(Difluoroboryloxy)phenyl-methylene](6-methoxy-1,3-benzothiazol-2-yl)amine, OMe/H) yield 63.4%. Mp 201–202 °C, green powder, ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.40 (d, 2H, *J* = 8.50 Hz), 7.66 (t, 1H, *J* = 7.56 Hz), 7.53 (dd, 2H, *J* = 7.82 Hz), 7.28 (d, 1H, *J* = 2.65 Hz), 7.20 (dd, 1H, *J* = 9.04, 2.16 Hz), 3.92 (s, 3H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 1.04 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 172.2, 167.7, 158.6, 134.2, 134.0, 131.1, 130.9, 128.6, 119.4, 117.1, 105.2, 56.0. ¹⁹F (376 MHz, from CFCl₃, DMSO-d₆): δ (ppm) –137.0. EI HRMS (*m/z*) [*M*]⁺ 332.0605 cal. 332.0602.

**4.1.3 Compound 3.** ((Z)-[(Difluoroboryloxy)[p-(trifluoromethyl)phenyl]methylene(6-methoxy-1,3-benzothiazol-2-yl)amine, OMe/CF₃) yield 80.2%. Mp 251–252 °C, green powder ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.51 (d, 2H, *J* = 8.08 Hz), 8.01 (d, 1H, *J* = 9.08 Hz), 7.80 (d, 2H, *J* = 8.20 Hz), 7.31 (d, 1H, *J* = 2.36 Hz), 7.24 (dd, 1H, *J* = 9.12, 2.52 Hz), 3.94 (s, 3H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 0.64 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 171.7, 166.1, 158.8, 135.3, 135.1, 134.5, 133.9, 130.5, 129.1, 125.6, 124.4, 122.8, 119.7, 117.5, 105.1, 56.0. ¹⁹F (376 MHz, from CFCl₃, DMSO-d₆): δ (ppm) –64.1, –136.6. EI HRMS (*m/z*) [*M*]⁺ 400.0478 cal. 400.0476.

**4.1.4 Compound 4.** ((Z)-[(Difluoroboryloxy)(p-methoxyphenyl)methylene]-1,3-benzothiazol-2-ylamine, H/OMe) yield 82.7%. Mp 215–216 °C, green powder ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.39 (d, 2H, *J* = 9.02 Hz), 8.06 (d, 1H, *J* = 8.08 Hz), 7.81 (d, 1H, *J* = 8.36 Hz), 7.61 (m, 1H), 7.49 (m, 1H), 7.02 (d, 2H, *J* = 9.02 Hz), 3.94 (s, 3H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 1.03 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 174.0, 168.1, 165.0, 140.1, 132.9, 131.7, 128.2, 126.2, 123.1, 122.1, 118.4, 114.1, 55.6. ¹⁹F (376 MHz, from CFCl₃, CDCl₃): δ (ppm) –137.4. EI HRMS (*m/z*) [*M*]⁺ 332.0601 cal. 332.0602.

**4.1.5 Compound 5.** ((Z)-[(Difluoroboryloxy)phenyl-methylene]-1,3-benzothiazol-2-ylamine, H/H) yield 66.8%. Mp 197–198 °C, white powder ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.42 (dd, 2H, *J* = 8.55 Hz), 8.09 (d, 1H, *J* = 9.08 Hz), 7.83 (d, 1H, *J* = 8.68 Hz), 7.67 (m, 1H), 7.63 (m, 1H), 7.54 (m, 2H) 7.51 (m, 1H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 1.10 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 174.0, 168.5, 140.1, 134.5, 131.0, 128.7, 128.4, 127.2, 126.6, 122.2, 118.7. ¹⁹F

(376 MHz, from CFCl₃, CDCl₃): δ (ppm) –136.7. EI HRMS (*m/z*) [*M*]⁺ 302.0501 cal. 302.0497.

**4.1.6 Compound 6.** ((Z)-[(Difluoroboryloxy)[p-(trifluoromethyl)phenyl]methylene-1,3-benzothiazol-2-ylamine, H/CF₃) yield 90.1%. Mp 176–177 °C, white powder ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.52 (d, 2H, *J* = 8.15 Hz), 8.12 (d, 1H, *J* = 8.17 Hz), 7.86 (d, 2H, *J* = 8.20 Hz), 7.80 (d, 2H, *J* = 2.36 Hz), 7.67 (dd, 1H, *J* = 7.58 Hz), 7.67 (dd, 1H, *J* = 7.32 Hz). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 1.10 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 166.8, 139.9, 135.6, 135.1, 134.3, 130.6, 128.6, 125.6, 125.5, 125.3, 122.2, 121.7, 118.8. ¹⁹F (376 MHz, from CFCl₃, DMSO-d₆): δ (ppm) –64.1, –136.6. EI HRMS (*m/z*) [*M*]⁺ 370.0375 cal. 370.0371.

**4.1.7 Compound 7.** ((Z)-[(Difluoroboryloxy)(p-methoxyphenyl)methylene](6-(trifluoromethyl)-1,3-benzothiazol-2-yl)amine, CF₃/OMe) yield 53.3%. Mp 288–289 °C, green powder ¹H NMR (700 MHz, from TMS, CDCl₃): δ (ppm) 8.38 (d, 2H, *J* = 8.89 Hz), 8.12 (d, 1H, *J* = 8.61 Hz), 8.08 (s, 1H), 7.83 (dd, 1H, *J* = 8.70, 1.42), 7.00 (d, 1H, *J* = 9.10 Hz), 3.93 (s, 3H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 0.66 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 175.6, 169.2, 165.2, 142.5, 133.3, 128.7, 127.2, 125.3, 122.8, 119.7, 118.7, 114.3, 55.65. ¹⁹F (376 MHz, from CFCl₃, DMSO-d₆): δ (ppm) –62.85. EI HRMS (*m/z*) [*M*]⁺ 400.0478 cal. 400.0476.

**4.1.8 Compound 8.** ((Z)-[(Difluoroboryloxy)phenyl-methylene](6-(trifluoromethyl)-1,3-benzothiazol-2-yl)amine, CF₃/H) yield 77.9%. Mp 249–250 °C, white powder ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.44 (d, 2H, *J* = 8.08 Hz), 8.19 (d, 1H, *J* = 9.08 Hz), 8.14 (m, 1H), 7.88 (dd, 1H, *J* = 8.78, 1.31 Hz), 7.71 (m, 1H), 7.57 (m, 2H), 3.94 (s, 3H). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 0.73 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 175.8, 169.5, 142.3, 135.1, 130.8, 130.6, 129.1, 128.8, 128.7, 127.4, 125.6, 125.3, 121.7, 119.9, 119.0, 118.1. ¹⁹F (376 MHz, from CFCl₃, DMSO-d₆): δ (ppm) –62.9, –136.6. EI HRMS (*m/z*) [*M*]⁺ 370.0375 cal. 370.0371.

**4.1.9 Compound 9.** ((Z)-[(Difluoroboryloxy)[p-(trifluoromethyl)phenyl]methylene](6-(trifluoromethyl)-1,3-benzothiazol-2-yl)amine, CF₃/CF₃) yield 85.0%. Mp 237–238 °C, white powder ¹H NMR (400 MHz, from TMS, CDCl₃): δ (ppm) 8.53 (d, 2H, *J* = 8.23 Hz), 8.20 (d, 1H, *J* = 8.83 Hz), 7.80 (d, 1H, *J* = 8.20 Hz), 7.31 (m, 1H), 7.24 (dd, 1H, *J* = 8.77, 1.19 Hz), 7.80 (d, 2H, *J* = 8.49 Hz). ¹¹B (128 MHz, from BF₃·Et₂O, CDCl₃): δ (ppm) 0.65 (t). ¹³C (100 MHz, from TMS, CDCl₃): δ (ppm): 175.6, 168.0, 142.2, 136.1, 134.0, 131.0, 129.1, 127.7, 125.8, 124.8, 122.1, 120.0, 119.3. ¹⁹F (376 MHz, from CFCl₃): δ (ppm) –64.2, –62.9, –135.8. EI HRMS (*m/z*) [*M*]⁺ 438.0239 cal. 438.0244.

## 4.2 Optical measurements

UV-vis absorption spectra of *ca.* 10⁻⁵ M solutions of dyes were recorded in quartz cells (1 cm) using a Shimadzu UV-1900 spectrometer. Fluorescence spectra of *ca.* 10⁻⁶ M solutions of complexes were recorded in quartz cells (1 cm) using an FS5 fluorimeter (Edinburgh Instruments). Fluorescence quantum yield measurements were carried out on the same equipment with the use of an integrating sphere. Time-correlated single-photon counting measurements were performed with an Edinburgh
<figure><img src="image_15.png" /></figure>

Analytical Instruments F920P spectrometer. Samples were excited at 373 nm using a laser diode. The emission intensity was recorded at fluorescence maximum wavelength. A solution of colloidal silica was used to obtain the instrument response function (IRF). Fluorescence lifetimes were calculated using the FAST software.

## 4.3 Calculations

Geometry optimizations of the ground and excited state have been performed within the (TD-)MN15 (ref. 54) approach known to deliver especially accurate energies for compounds of the family investigated here⁵⁵ using either the 6-311++G(d,p) or 6-311+G(d,p) basis sets, the former being applied for CT parameters and rotational barrier, the (more compact) latter for (the more expensive) vibronic calculations. The solvent was modelled through the well-known Polarizable Continuum Mode (PCM),⁵⁶ selecting chloroform. The character of the stationary points on the potential energy surface has been confirmed with the harmonic vibrational analysis. The rotation barriers have been estimated by the optimization of the geometry of the transition state for the phenyl ring rotation around the C15–C17 bond (for atom numbering see Fig. S52†). All these calculations have been carried out with Gaussian16.⁵⁷ The SOC matrix elements were determined at the M06-2X/def2-TZVP level⁵⁸ using ORCA5.0.1 (ref. 59) with CHCl₃ as the solvent as modelled with the SMD solvation model. The RIJCOSX method was used to accelerate the calculation, and DefGrid3 and TightSCF settings were applied, whereas TDA was turned off. The S–T gaps were evaluated using SCS-CC2/aug-cc-pVTZ (in the gas phase)⁶⁰,⁶¹ because this method performs especially well for these gaps.⁶² These SCS-CC2 calculations were performed with Turbomole 7.3.⁶³ Given that ISC typically occurs after geometry relaxation in organic dyes (free of heavy elements), these calculations have been performed on the ES minima. For the MECP search, we used the Q-Chem 5.4 code,⁶⁴ applying the same level of theory as above, that is PCM-MN15/6-311+G(d,p), using the spin-flip approach. We considered four roots, applied the branching plane approach to locate the MECP, and searched for the MECP between the two lowest singlet states, defined with a spin-contamination threshold of 1.2. Default Q-Chem parameters were otherwise applied. In practice the MECP search was first performed with the non-substituted 5, exploring a large number of possible deformations as the starting point (kinked structures twisting various moieties, hydrogen atoms moved out of plane, elongated bonds, etc.). It is indeed our experience that starting from the FC geometry does not necessarily yield the lowest MECP, especially for quasiplanar systems. This allowed the determination of the most accessible MECP. Starting from there, the MECPs of all other compounds were searched for. The energy difference between the MECP and the excited-state minimum was determined through PCM-SF-TD-MN15/6-311+G(d,p) single point calculations performed on the optimal excited-state structure found by conventional TD-DFT (see above) and on the SF-TD-DFT optimized MECP. All the vibronic calculations have been carried out using FCClasses 3 (v3-0.1-216-g29ba701).⁶⁵ Both the adiabatic and vertical Hessian PES approximations (AH and VH, respectively) have been

originally used along with the Time-Dependent (TD)⁶⁶ and Time-Independent (TI)⁶⁷,⁶⁸ formalisms (see also below). Considering the large oscillator strength values ($f > 0.5$, typically ca. 0.8–1.2), the Herzberg–Teller effects were neglected. We used geometries, gradients, and Hessians computed at the PCM-MN15/6-311+G(d,p) level, and further corrected the energies and transition dipole moment using gas-phase CC2/aug-cc-pVTZ results. These CC2 calculations were performed in an effort to improve the accuracy and were made with Turbomole. Regarding the pre-screening parameters of FCClasses we used $C_1^{max} = 30$, $C_2^{max} = 25$ for the maximum number of integrals in both $C_1$ and $C_2$ classes and we set the threshold to select the thermally excited initial states to 0.8. The number of maximum integrals was set to $10^6$ for the routine TI calculation. All the emission spectra were convoluted with a Gaussian broadening of $450 \text{ cm}^{-1}$ to ensure matching with experimental data. The Internal Conversion (IC) spectra use a $10 \text{ cm}^{-1}$ Lorentzian broadening, a typical choice.⁶⁹ In order to select the appropriate vibronic model we carried out some tests on 4, H/OMe (the corresponding spectra shown in the ESI†). Since a better convergence of the TI spectra was obtained within the AH approximation, as well as the expected shape of the $k_{IC}$ versus energy plot, we pursued the work using TD AH|FC formalism, which we report in the body of the text.

## Data availability

Additional computational data are available upon request to DJ while experimental ones from BO.

## Author contributions

PR: investigation, formal analysis, validation, visualization, writing – original draft preparation; MHEB: investigation, formal analysis, validation, visualization, writing – original draft preparation; AKK: investigation, visualization, writing – original draft preparation; BJ: investigation, formal analysis, writing – review & editing; DJ: conceptualization, methodology, resources, writing – original draft preparation, supervision; BO: conceptualization, project administration, methodology, resources, writing – original draft preparation, supervision.

## Conflicts of interest

There are no conflicts to declare.

## Acknowledgements

MHEB and DJ are indebted to Dr R. Zaleśny (Wrocław) for fruitful discussions on the convergence of vibronic calculations and to Dr F. Santoro and J. Cerezo for making their great code available to them. We are very indebted for the financial support of research from the National Science Centre (project no. 2021/43/B/ST5/00753). Wrocław Centre for Networking and Supercomputing and Nantes' Glicid/CCIPL computational centers are gratefully acknowledged for the generous allotment of computational resources.
<figure><img src="image_16.png" /></figure>

## Notes and references

1. B. Jędrzejewska, A. Skotnicka, A. D. Laurent, M. Pietrzak, D. Jacquemin and B. Ośmiałowski, *J. Org. Chem.*, **2018**, *83*, 7779-7788.

2. R. Hu, E. Lager, A. Aguilar-Aguilar, J. Liu, J. W. Y. Lam, H. H. Y. Sung, I. D. Williams, Y. Zhong, K. S. Wong, E. Peña-Cabrera and B. Z. Tang, *J. Phys. Chem. C*, **2009**, *113*, 15845-15853.

3. Y. Kage, S. Kang, S. Mori, M. Mamada, C. Adachi, D. Kim, H. Furuta and S. Shimizu, *Chem.-Eur. J.*, **2021**, *27*, 5259-5267.

4. Z. Shi, X. Han, W. Hu, H. Bai, B. Peng, L. Ji, Q. Fan, L. Li and W. Huang, *Chem. Soc. Rev.*, **2020**, *49*, 7533-7567.

5. R. Ziessel, G. Ulrich and A. Harriman, *New J. Chem.*, **2007**, *31*, 496-501.

6. A. Loudet and K. Burgess, *Chem. Rev.*, **2007**, *107*, 4891-4932.

7. P. Kaur and K. Singh, *J. Mater. Chem. C*, **2019**, *7*, 11361-11405.

8. Q. Li, Y. Guo and S. Shao, *Sens. Actuators, B*, **2012**, *171-172*, 872-877.

9. Z. Shi, X. Han, W. Hu, H. Bai, B. Peng, L. Ji, Q. Fan, L. Li and W. Huang, *Chem. Soc. Rev.*, **2020**, *49*, 7533-7567.

10. S. Kolemen and E. U. Akkaya, *Coord. Chem. Rev.*, **2018**, *354*, 121-134.

11. J. Zou, Z. Yin, K. Ding, Q. Tang, J. Li, W. Si, J. Shao, Q. Zhang, W. Huang and X. Dong, *ACS Appl. Mater. Interfaces*, **2017**, *9*, 32475-32481.

12. A. Kamkaew, S. H. Lim, H. B. Lee, L. V. Kiew, L. Y. Chung and K. Burgess, *Chem. Soc. Rev.*, **2013**, *42*, 77-88.

13. A. M. Kosswattaarachchi, A. E. Friedman and T. R. Cook, *ChemSusChem*, **2016**, *9*, 3317-3323.

14. E. Deniz, G. C. Isbasar, Ö. A. Bozdemir, L. T. Yildirim, A. Siemiarczuk and E. U. Akkaya, *Org. Lett.*, **2008**, *10*, 3401-3403.

15. P. Verwilst, H.-R. Kim, J. Seo, N.-W. Sohn, S.-Y. Cha, Y. Kim, S. Maeng, J.-W. Shin, J. H. Kwak, C. Kang and J. S. Kim, *J. Am. Chem. Soc.*, **2017**, *139*, 13393-13403.

16. B. Valeur, *Molecular Fluorescence: Principles and Applications*, Wiley-VCH, Weinheim, 2002.

17. F. Santoro and D. Jacquemin, *Wiley Interdiscip. Rev.: Comput. Mol. Sci.*, **2016**, *6*, 460-486.

18. A. Prlj, A. Fabrizio and C. Corminboeuf, *Phys. Chem. Chem. Phys.*, **2016**, *18*, 32668-32672.

19. A. Prlj, L. Vannay and C. Corminboeuf, *Helv. Chim. Acta*, **2017**, *100*, e1700093.

20. Z. Lin, A. W. Kohn and T. Van Voorhis, *J. Phys. Chem. C*, **2020**, *124*, 3925-3938.

21. Q. Ou, Q. Peng and Z. Shuai, *J. Phys. Chem. Lett.*, **2020**, *11*, 7790-7797.

22. A. Humeniuk, M. Bužančić, J. Hoche, J. Cerezo, R. Mitrić, F. Santoro and V. Bonačić-Koutecký, *J. Chem. Phys.*, **2020**, *152*, 054107.

23. K. Veys and D. Escudero, *J. Phys. Chem. A*, **2020**, *124*, 7228-7237.

24. J. S. A. Badaro, B. Koszarna, M. H. E. Bousquet, E. T. Ouellette, D. Jacquemin and D. T. Gryko, *Org. Chem. Front.*, **2022**, *9*, 1861-1874.

25. K. V. Vygranenko, Y. M. Poronik, M. H. E. Bousquet, O. Vakuliuk, D. Jacquemin and D. T. Gryko, *Chem. Commun.*, **2022**, *58*, 1542-1545.

26. M. Krzeszewski, S. Modrzycka, M. H. E. Bousquet, D. Jacquemin, M. Drąg and D. T. Gryko, *Org. Lett.*, **2022**, *24*, 5602-5607.

27. C.-W. Ju, H. Bai, B. Li and R. Liu, *J. Chem. Inf. Model.*, **2021**, *61*, 1053-1065.

28. J. F. Joung, M. Han, J. Hwang, M. Jeong, D. H. Choi and S. Park, *JACS Au*, **2021**, *1*, 427-438.

29. H. Levine III, *Protein Sci.*, **1993**, *2*, 404-410.

30. M. Biancalana and S. Koide, *Biochim. Biophys. Acta*, **2010**, *1804*, 1405-1412.

31. H. Levine III, *Arch. Biochem. Biophys.*, **1997**, *342*, 306-316.

32. S. Verma, V. Ravichandiran and N. Ranjan, *Biochimie*, **2021**, *190*, 111-123.

33. N. Amdursky, Y. Erez and D. Huppert, *Acc. Chem. Res.*, **2012**, *45*, 1548-1557.

34. A. I. Sulatskaya, A. A. Maskevich, I. M. Kuznetsova, V. N. Uversky and K. K. Turoverov, *PLoS One*, **2010**, *5*, 1-7.

35. G. R. Suman, M. Pandey and A. S. J. Chakravarthy, *Mater. Chem. Front.*, **2021**, *5*, 1541-1584.

36. J. Mei, N. L. C. Leung, R. T. K. Kwok, J. W. Y. Lam and B. Z. Tang, *Chem. Rev.*, **2015**, *115*, 11718-11940.

37. L. S. Wolfe, M. F. Calabrese, A. Nath, D. V. Blaho, A. D. Miranker and Y. Xiong, *Proc. Natl. Acad. Sci. U. S. A.*, **2010**, *107*, 16863-16868.

38. V. I. Stsiapura, A. A. Maskevich, S. A. Tikhomirov and O. V. Buganov, *J. Phys. Chem. A*, **2010**, *114*, 8345-8350.

39. V. I. Stsiapura, A. A. Maskevich, V. A. Kuzmitsky, V. N. Uversky, I. M. Kuznetsova and K. K. Turoverov, *J. Phys. Chem. B*, **2008**, *112*, 15893-15902.

40. H. F. Kung, *ACS Med. Chem. Lett.*, **2012**, *3*, 265-267.

41. N. S. Mason, C. A. Mathis and W. E. Klunk, *J. Labelled Compd. Radiopharm.*, **2013**, *56*, 89-95.

42. S. Noël, S. Cadet, E. Gras and C. Hureau, *Chem. Soc. Rev.*, **2013**, *42*, 7747-7762.

43. S. Gao, Y. Ma and W. Lin, *Anal. Methods*, **2019**, *11*, 2626-2629.

44. M.-C. Chen, D.-G. Chen and P.-T. Chou, *ChemPlusChem*, **2021**, *86*, 11-27.

45. K. Skonieczny, I. Papadopoulos, D. Thiel, K. Gutkowski, P. Haines, P. M. McCosker, A. D. Laurent, P. A. Keller, T. Clark, D. Jacquemin, D. M. Guldi and D. T. Gryko, *Angew. Chem., Int. Ed.*, **2020**, *59*, 16104-16113.

46. S. Ye, H. Zhang, J. Fei, C. H. Wolstenholme and X. Zhang, *Angew. Chem., Int. Ed.*, **2021**, *60*, 1339-1346.

47. A. M. Grabarz, B. Jędrzejewska, A. Skotnicka, N. A. Murugan, F. Patalas, W. Bartkowiak, D. Jacquemin and B. Ośmiałowski, *Dyes Pigment*, **2019**, *170*, 107481.

48. T. Le Bahers, C. Adamo and I. Ciofini, *J. Chem. Theory Comput.*, **2011**, *7*, 2498-2506.

49. I. Ciofini, T. Le Bahers, C. Adamo, F. Odobel and D. Jacquemin, *J. Phys. Chem. C*, **2012**, *116*, 11946-11955.
<figure><img src="image_17.png" /></figure>

50 Y. Harabuchi, T. Taketsugu and S. Maeda, *Phys. Chem. Chem. Phys.*, 2015, **17**, 22561-22565.

51 J. Shi, M. A. Izquierdo, S. Oh, S. Y. Park, B. Milián-Medina, D. Roca-Sanjuán and J. Gierschner, *Org. Chem. Front.*, 2019, **6**, 1948-1954.

52 Q. Peng, Y. Niu, C. Deng and Z. Shuai, *Chem. Phys.*, 2010, **370**, 215-222.

53 Z. Shuai, *Chin. J. Chem.*, 2020, **38**, 1223-1232.

54 H. S. Yu, X. He, S. L. Li and D. G. Truhlar, *Chem. Sci.*, 2016, **7**, 5032-5051.

55 A. M. Grabarz and B. Ośmiałowski, *Molecules*, 2021, **26**, 7434.

56 J. Tomasi, B. Mennucci and R. Cammi, *Chem. Rev.*, 2005, **105**, 2999-3094.

57 M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery Jr, J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman and D. J. Fox, *Gaussian 16 Revision C.01*, Gaussian Inc., Wallingford CT, 2016.

58 Y. Zhao and D. G. Truhlar, *Theor. Chem. Acc.*, 2008, **120**, 215-241.

59 F. Neese, F. Wennmohs, U. Becker and C. Riplinger, *J. Chem. Phys.*, 2020, **152**, 224108.

60 O. Christiansen, H. Koch and P. Jørgensen, *Chem. Phys. Lett.*, 1995, **243**, 409-418.

61 A. Hellweg, S. A. Grün and C. Hättig, *Phys. Chem. Chem. Phys.*, 2008, **10**, 4119-4127.

62 A. Pershin, D. Hall, V. Lemaur, J.-C. Sancho-Garcia, L. Muccioli, E. Zysman-Colman, D. Beljonne and Y. Olivier, *Nat. Commun.*, 2019, **10**, 597.

63 S. G. Balasubramani, G. P. Chen, S. Coriani, M. Diedenhofen, M. S. Frank, Y. J. Franzke, F. Furche, R. Grotjahn, M. E. Harding, C. Hättig, A. Hellweg, B. Helmich-Paris, C. Holzer, U. Huniar, M. Kaupp, A. Marefat Khah, S. Karbalaei Khani, T. Müller, F. Mack, B. D. Nguyen, S. M. Parker, E. Perlt, D. Rappoport, K. Reiter, S. Roy, M. Rückert, G. Schmitz, M. Sierka, E. Tapavicza, D. P. Tew, C. van Wüllen, V. K. Voora, F. Weigend, A. Wodyński and J. M. Yu, *J. Chem. Phys.*, 2020, **152**, 184107.

64 E. Epifanovsky, A. T. B. Gilbert, X. Feng, J. Lee, Y. Mao, N. Mardirossian, P. Pokhilko, A. F. White, M. P. Coons, A. L. Dempwolff, Z. Gan, D. Hait, P. R. Horn, L. D. Jacobson, I. Kaliman, J. Kussmann, A. W. Lange, K. U. Lao, D. S. Levine, J. Liu, S. C. McKenzie, A. F. Morrison, K. D. Nanda, F. Plasser, D. R. Rehn, M. L. Vidal, Z.-Q. You, Y. Zhu, B. Alam, B. J. Albrecht, A. Aldossary, E. Alguire, J. H. Andersen, V. Athavale, D. Barton, K. Begam, A. Behn, N. Bellonzi, Y. A. Bernard, E. J. Berquist, H. G. A. Burton, A. Carreras, K. Carter-Fenk, R. Chakraborty, A. D. Chien, K. D. Closser, V. Cofer-Shabica, S. Dasgupta, M. de Wergifosse, J. Deng, M. Diedenhofen, H. Do, S. Ehlert, P.-T. Fang, S. Fatehi, Q. Feng, T. Friedhoff, J. Gayvert, Q. Ge, G. Gidofalvi, M. Goldey, J. Gomes, C. E. González-Espinoza, S. Gulania, A. O. Gunina, M. W. D. Hanson-Heine, P. H. P. Harbach, A. Hauser, M. F. Herbst, M. H. Vera, M. Hodecker, Z. C. Holden, S. Houck, X. Huang, K. Hui, B. C. Huynh, M. Ivanov, Á. Jász, H. Ji, H. Jiang, B. Kaduk, S. Kähler, K. Khistyaev, J. Kim, G. Kis, P. Klunzinger, Z. Koczor-Benda, J. H. Koh, D. Kosenkov, L. Koulias, T. Kowalczyk, C. M. Krauter, K. Kue, A. Kunitsa, T. Kus, I. Ladjánszki, A. Landau, K. V. Lawler, D. Lefrancois, S. Lehtola, R. R. Li, Y.-P. Li, J. Liang, M. Liebenthal, H.-H. Lin, Y.-S. Lin, F. Liu, K.-Y. Liu, M. Loipersberger, A. Luenser, A. Manjanath, P. Manohar, E. Mansoor, S. F. Manzer, S.-P. Mao, A. V. Marenich, T. Markovich, S. Mason, S. A. Maurer, P. F. McLaughlin, M. F. S. J. Menger, J.-M. Mewes, S. A. Mewes, P. Morgante, J. W. Mullinax, K. J. Oosterbaan, G. Paran, A. C. Paul, S. K. Paul, F. Pavošević, Z. Pei, S. Prager, E. I. Proynov, Á. Rák, E. Ramos-Cordoba, B. Rana, A. E. Rask, A. Rettig, R. M. Richard, F. Rob, E. Rossomme, T. Scheele, M. Scheurer, M. Schneider, N. Sergueev, S. M. Sharada, W. Skomorowski, D. W. Small, C. J. Stein, Y.-C. Su, E. J. Sundstrom, Z. Tao, J. Thirman, G. J. Tornai, T. Tsuchimochi, N. M. Tubman, S. P. Veccham, O. Vydrov, J. Wenzel, J. Witte, A. Yamada, K. Yao, S. Yeganeh, S. R. Yost, A. Zech, I. Y. Zhang, X. Zhang, Y. Zhang, D. Zuev, A. Aspuru-Guzik, A. T. Bell, N. A. Besley, K. B. Bravaya, B. R. Brooks, D. Casanova, J.-D. Chai, S. Coriani, C. J. Cramer, G. Cserey, A. E. DePrince, R. A. DiStasio, A. Dreuw, B. D. Dunietz, T. R. Furlani, W. A. Goddard, S. Hammes-Schiffer, T. Head-Gordon, W. J. Hehre, C.-P. Hsu, T.-C. Jagau, Y. Jung, A. Klamt, J. Kong, D. S. Lambrecht, W. Liang, N. J. Mayhall, C. W. McCurdy, J. B. Neaton, C. Ochsenfeld, J. A. Parkhill, R. Peverati, V. A. Rassolov, Y. Shao, L. V. Slipchenko, T. Stauch, R. P. Steele, J. E. Subotnik, A. J. W. Thom, A. Tkatchenko, D. G. Truhlar, T. Van Voorhis, T. A. Wesolowski, K. B. Whaley, H. L. Woodcock, P. M. Zimmerman, S. Faraji, P. M. W. Gill, M. Head-Gordon, J. M. Herbert and A. I. Krylov, *J. Chem. Phys.*, 2021, **155**, 084801.

65 F. Santoro and J. Cerezo, *FCclasses3*, a code for vibronic calculations. Version 1.0 (v3-0.1-177-g4b1514a), August 2021, visit, http://www.pi.iccom.cnr.it/fcclasse.

66 F. J. Avila Ferrer, J. Cerezo, J. Soto, R. Improta and F. Santoro, *Comput. Theor. Chem.*, 2014, **1040-1041**, 328-337.

67 F. Santoro, R. Improta, A. Lami, J. Bloino and V. Barone, *J. Chem. Phys.*, 2007, **126**, 084509.

68 F. Santoro, R. Improta, A. Lami, J. Bloino and V. Barone, *J. Chem. Phys.*, 2007, **126**, 169903.

69 A. Humeniuk, M. Bužančić, J. Hoche, J. Cerezo, R. Mitrić, F. Santoro and V. Bonačić-Koutecký, *J. Chem. Phys.*, 2020, **152**, 054107.