Using the isotope effect to probe an aggregation induced emission mechanism: theoretical prediction and experimental validation.

Aggregation-induced emission (AIE) has become a hot topic for a variety of potential applications, but the understanding of its working mechanism is still under scrutiny. Herein, we proposed the use of the isotope effect (IE) to identify the AIE mechanism: under the restriction of an internal motion mechanism, the IE is pronouncedly different in excited-state decay rates when contrasting AIE luminogens (AIEgens) and non-AIEgens in theoretical calculations. For the complete deuteration of AIEgens, the IE of nonradiative decay rate in solution (<-10%) is much weaker than that (-65% to -95%) in aggregate, because the former stems from the overall results of competitive vibronic coupling and the severe mixing of low-frequency modes while the latter mainly comes from the vibronic coupling only. The experimental results confirm the isotopic "jump" behaviors in AIEgens well. However, non-AIEgens exhibit equivalent IEs (-40% to -90%) in both solution and solid phases. Further partial deuteration schemes for the 6-ring AIE analogues show positional dependence.

Introduction

Molecular aggregation in general tends to reduce luminescence efficiency: intermolecular charge transfer or energy transfer to quenching sites, or the formation of a dipole-forbidden lowest excited state.1 However, a series of aggregation induced emissive luminogens (AIEgens) are found to be poorly luminescent in dilute solution but highly emissive in an aggregate or crystalline state, demonstrating potential application in organic solid-state lighting and cell imaging.2 Despite the varieties of experimental efforts to clarify the AIE mechanism, there are still controversies regarding the microscopic origin of AIE: restriction of intramolecular motion (RIM),3 J-aggregation formation,3b,c hydrogen-bond induced excimer-emission,3 intramolecular planarization,3 vibration induced emission,3 and so on. To better understand AIE, here we propose to use the isotope effect (IE) to probe the influence of molecular aggregation on the excited-state non-radiative decay from the lowest excited singlet state (S1) to the ground state (S0) through comparative computational studies for AIEgens with RIM mechanisms and conventional luminogens.

Deuteration has been widely applied to probe excited-state decay processes.4 For conventional fluorophores, it is known that deuteration always causes a decrease in the nonradiative decay rate while hardly changing the radiative decay rate, leading to an increase in luminescence efficiency. This point is easily understood from the internal conversion rate under the displaced harmonic oscillator model based on the Fermi golden rule, where the kic is exponentially proportional to –S, the Huang–Rhys factor of the j-th normal mode, which is determined by the normal mode reorganization energy λ = ħSω.5 According to the four-point method in the potential energy surface, the total reorganization energy summated over all the normal modes is also defined as λg(e) = Eg(e)(S1(0)-geometry) – Eg(e)(S0(1)-geometry) in the ground (excited) state. Thus, it is independent of isotopic substitution since the equilibrium geometry and electronic-state energy are the same for different isotopomers.6 The lower frequency ω induced by the deuteration of the modes implies the increase of the Huang–Rhys factor S, which reduces the kic with the displaced approximation.

For the flexible AIEgens, mixing between the low-frequency normal modes is significant and the kic cannot be correctly described under the displaced harmonic oscillator model. Therefore, we have developed a multimode coupled kic formalism including the Duschinsky rotation effect (DRE)7 (see Part III of ESI†). The Duschinsky rotation matrix (DRM) M correlates the normal coordinates of the ground state S0 (Qg) and the excited state S1 (Qe) as Qe = MQg + De. The DRE occurs most notably for low-frequency modes and becomes more remarkable when more modes with lower frequency are activated, which significantly increases the kic according to our previous investigations.7 Hence, the frequency reduction from isotope substitution gives rise to two competitive effects on kic: a negative effect through increasing the Huang–Rhys factor and a positive effect via strengthening of the DRE to enhance inter-mode mixing (Chart 1). The IE on kic can be defined as,where kDic and kHic are the internal conversion rates for the deuterated and hydrogen system, respectively. As the low-frequency modes are suppressed in the solid/rigid matrix for AIEgens, the DRE becomes minor, and the positive effect diminishes. Then, a more remarkable IE can be anticipated in the solid state compared to in solution. While for the non-AIEgens, since the nonradiative decay is mainly controlled by the vibronic coupling of high-frequency vibration modes with little DRE, the IE is expected to be normal and behaves similarly in the solid state and in solution. In this work, we aim to make such a proposal to use IE for elucidating the microscopic mechanism of the AIE phenomena through theoretical prediction followed by experimental validation.

Chart 1

Representation of isotope effect on kic.

We investigated the isotopic substitution effects on the luminescent properties in solution and aggregation of the AIE-active 1,1,2,3,4,5-hexaphenylsilole (HPS),8 1,1-bis(2′-thienyl)-2,3,4,5-tetraphenylsilole (BtTPS),9 1,2-diphenyl-3,4-bis (diphenylmethylene)-1-cyclobutene (HPDMCb),10cis,cis-1,2,3,4-tetraphenyl-1,3-butadiene (TPBD),11 and 2,2′-(6,12-diphenyltetracene-5,11-diyl)dithiophene (DPTDTP),12 in comparison with AIE-inactive diacenaphtho-[1,2-b;1′,2′-d]silole (BPS),13 perylene,14 distyrylbenzene (DSB),15 and anthracene14 (see Chart 2).

Chart 2

Molecular structures of the AIE and non-AIE compounds.

Results and discussion

The equilibrium geometries and frequencies at the S0(S1) state of the compounds were calculated at the (TD)-DFT/PBE0/6-31G(d) level, which has been proven to be an effective method to study the electronic structures of the organic molecules.16 The environmental effects are taken into account by employing a polarizable continuum model (PCM)17 for the solution phase and hybrid quantum mechanics/molecular mechanics (QM/MM)18 approach for the aggregates (computational models are given in Fig. 1 and Part I of ESI†). To better match the experimental absorption and emission maxima, we adopted the long-range separated CAM-B3LYP functional to recalculate the excitation energy of the optimized structure. The different PCM models, such as the linear-response (LR)19 and state-specific (SS)19 solvation methods were tested (see Part II and Tables S1–S9 of ESI†). As expected, the SS-PCM model performs better in describing the excitation energy, probably owing to the self-consistency between the solvent reaction field and the solute electrostatic potential, and has been widely used to investigate the optical properties of solution-phase organic molecules.20,21 The resultant vertical excitation energies and the comparisons between experimental and theoretical values are given in Table 1, and all deviations are within the acceptable range of 0.01–0.43 eV.22

Fig. 1

Setup of PCM (a) and QM/MM (b) models (taking BtTPS as an example).

Table 1

Calculated optical spectral properties (eV) for the AIE-active and AIE-inactive compounds in both solution and solid phases, as well as the experimental values in parenthesis for comparison

	HPS		BtTPS		HPDMCb
	Solution	Solid	Solution	Solid	Solution	Solid
Abs.	3.44 (3.39) a	3.39 (3.36) b	3.43 (3.33) d	3.44 (3.36) d	3.66 (—)	3.61 (3.50) e
Emi.	2.38 (2.48) a	2.51 (2.68) c	2.33 (—)	2.48 (2.59) d	2.38 (—)	2.49 (2.62) e

In acetone.8

In thin film.23

In crystal.23

In acetone and crystal.9

In crystal.10

In acetone and aggregate.11

In THF and aggregate.12

In CH2Cl2.13

In cyclohexane.23

In crystal.14

In CH2Cl2.15

In crystal.15

In cyclohexane.23

In crystal.14

We assume the deuterium substitutions hardly change the crystal structure and neglect the induced effects caused by environmental (MM part) deuteration. Normal mode analyses before and after deuteration were performed with the help of the DUSHIN24 program. Finally, we evaluated kic using the multimode coupled thermal vibration correlation function formalism realized in the home-built MOMAP program,25 which has successfully described the optical properties of many polyatomic molecules, including AIEgens and non-AIEgens.26 Both the distortion and DRE of the potential energy surfaces are taken into account in the kic formula (see Part III of ESI†).

The calculated results for kic as well as the IE results are presented in Table 2. It is seen that the kic decreases sharply by several orders of magnitude from the solution to solid phases for AIEgens but undergoes a slight change for non-AIEgens. We then plotted the IE results as shown in Fig. 2. It is clear that the IE is always negative. For the AIEgens, the IE is minor (less than ca. –10%) in solution but becomes remarkable in a solid phase (ca. –65% to –95%). Interestingly, the IE is strikingly different for the AIEgens compared to the non-AIEgens. For non-AIEgens, the IE results both in solution and solid phase are close to each other and fall within the range of ca. –40% to –90%.

Table 2

Calculated room-temperature kic (s–1) for non-deuterated (H-all) and fully-deuterated (D-all) isotopomers of the AIEgens and non-AIEgens in both solution and solid phases

	HPS		BtTPS		HPDMCb
	Solution	Solid	Solution	Solid	Solution	Solid
H-all	2.44 × 1011	8.60 × 106	2.20 × 1011	2.73 × 107	1.31 × 1011	2.26 × 107
D-all	2.22 × 1011	2.61 × 106	1.97 × 1011	6.89 × 106	1.27 × 1011	7.11 × 106
IE	–9.0%	–69.6%	–10.5%	–74.8%	–3.1%	–68.5%

Fig. 2

IE results for AIEgens and non-AIEgens.

These results can be well justified under the RIM mechanism. Since the internal conversion rate kic depends exponentially on –S, lowering ω implies an increased vibronic coupling strength since the total reorganization energy is unchanged upon deuteration. Thus, the negative IE on kic is generally considered as normal behavior. But for AIEgens in solution, the DRE caused by the low-frequency mode mixing plays a very important role in increasing the internal conversion rate.7 Lessening the frequency by isotope substitutions leads to a more pronounced DRE, which tends to sharply increase kic in solution. This counteracts the general IE-induced decrease in kic to a large extent. Thus, full deuteration leads to a much lesser IE for AIEgens in solution, distinct from AIEgens in solid phase and non-AIEgens in solution or solid phases.

We further introduce an effective frequency (ωeff) for a more quantitative interpretation,27which represents an effective mode contributing to the excited-state relaxation process. The ωeff results for the pristine systems are shown in Fig. 3. It is noted that the ωeff for AIEgens falls into the range of 100–160 cm–1 in solution since the low-frequency modes (<100 cm–1) contribute largely to the relaxation energy. But upon aggregation, the ωeff rises to 300–1000 cm–1 owing to the dominating high frequency vibrations. The ωeff (600–1200 cm–1) is close for all the non-AIEgens in solution and solid states. In addition, the deuteration effect for ωeff is given in Table S10† and the ωeff always decreases as expected. Detailed analyses of vibrational and structural relaxation modification induced by aggregation are presented in Part IV, Fig. S1–S7 and Tables S11–S18 of the ESI.†

Fig. 3

ω eff for AIEgens and non-AIEgens in both solution and solid phases.

The difference of the IE for solution-phase AIEgens compared to other cases mainly stems from the outstanding contributions of low-frequency normal modes during the excited-state decay process. The low-frequency modes arise from the flat potential energy surfaces, and the involved quantum number of these modes increases sharply with temperature. At room temperature, there are so many vibrational states which are close in energy. They tend to strongly mix among the different modes, namely, the DRE becomes very pronounced, leading to the increase of kic. The DRE can be directly characterized from the calculated DRM. The more the diagonal elements of the DRM deviate from 1.0, the more off-diagonal elements are non-zero, and the more significant the DRE is. The contour maps of the DRM for 20 normal modes with the lowest frequencies are presented in Fig. 4. For both solution-phase and solid-phase non-AIEgens, e.g.BPS, most diagonal elements are close to 1.0 and the DRE is tiny. A similar situation can be found for the AIEgens in the solid-phase e.g. for HPS, most of the lowest 20 modes are in the range of 100–200 cm–1 and the DRM elements gather in the vicinity of the diagonal line, indicating a very weak DRE. However, for AIEgens in the solution phase, the off-diagonal elements of the DRM become more outstanding corresponding to more considerable mixing between modes. When the frequencies of the modes decrease upon deuteration, the effect would become much more severe, this could sharply increase the kic and compensate for the normal negative IE on kic. To validate this assumption, we also calculate kic without the DRE for a typical AIEgen, HPS, and non-AIEgen, BPS, in both solution and solid phases (see Table S19†). IEs without the DRE are –88.4% (–71.4%) for HPS in solution (solid phase) and –47.8% (–48.0%) for BPS in solution (solid phase), respectively. When the DRE is considered, it is only –9.0% for solution-phase HPS. Therefore, isotopic substitution could be utilized as an effective tool to explore the aggregation effect on the nonradiative process. The DRE suppressed by molecular aggregation in AIEgens could result in a sudden “jump” or “drop” of the IE.

Fig. 4

Contour maps of the DRM for the lowest 20 modes in HPS solution (a), HPS solid phase (b), BPS solution (c), and BPS solid phase (d). The excited-state normal mode index is rearranged to make the largest element the diagonal.

To experimentally verify the above calculated isotopic characteristic of AIEgens, we synthesized 2,3,4,5-ring deuterated (D-) HPS (see Part V of the ESI†) and measured the radiative and nonradiative decay rates for a comparison with the pristine HPS as presented in Table 3. The experimental values are in good agreement with the calculated ones (computational details presented in Part VI and Table S20 of the ESI†). In solution, both the H- HPS and D- HPS are non-emissive as determined from both the calculated and experimental results because the nonradiative decay rates are far larger than the radiative ones, while they emit strong light in the solid state. Strikingly, the experimental IEs well reproduce the calculated results, in that they are very little in solution but remarkable in solid phase. In particular, the abnormal positive effect of the IE is observed in the experiment, which fully confirms the essential role of the DRE in the nonradiative decay process. In solid phase, a large normal negative IE is observed because the nonradiative decay rate is mainly determined by the vibronic coupling.

Table 3

The calculated and the measured excited-state decay rates of the non-deuterated (H-) and deuterated (D-)2,3,4,5-ring HPS in solution and solid phases at room temperature

HPS	Solution			Solid
	k r	k ic	IE	k r	k ic	IE
H-	6.54 × 107  a	2.44 × 1011  a		0.83 × 107  c	0.86 × 107  c
	(1.27 × 107) b	(1.05 × 109) b		(1.52 × 107) d	(1.29 × 107) d
D-	6.54 × 107  a	2.30 × 1011  a	–5.7%	0.83 × 107  c	0.29 × 107  c	–66.3%
	(1.64 × 107) b	(1.07 × 109) b	(1.9%) b	(1.37 × 107) d	(0.47 × 107) d	(–63.6%) d

Calculated in acetone solution.

Measured in THF solution.

Calculated in crystal.

Measured in microcrystal.

We further systematically investigated the isotope substituted position dependence as to provide more experimentally verifiable quantitative predictions, taking the 6-ring analogues HPS, BtTPS, and HPDMCb as examples. Partially-deuterated schemes are considered in both solution and solid phases, namely, 1,1-rings or 1,6-rings only deuterated (D-1,1), 2,5-rings only deuterated (D-2,5), and 3,4-rings only deuterated (D-3,4). Detailed data are presented in Part VII and Table S21 of the ESI.† It is obvious that the IEs in solid states are more remarkable than those in solution. We herein analyse the cases in solid phase as illustrated in Fig. 5. For HPS, the IE follows the order of D-3,4 > D-2,5 > D-1,1; for BtTPS, D-2,5 induces a much larger IE than D-3,4 and D-1,1; HPDMCb also shows position-dependent deuteration effects, with a major IE induced by D-1,1 and D-3,4 but a minor IE induced by D-2,5. Interestingly, the IE decreases in the order of the reduction of ωeff (see Table S22†) for the partially deuterated isotopomers, which further confirms the ωeff to be a good parameter to judge the IE. E.g. for BtTPS, D-2,5 undergoes a remarkable frequency decrease (–3.5%), while D-3,4 experiences a slight frequency change (–0.6%), and D-1,1 exhibits an almost unaffected frequency (–0.3%).

Fig. 5

IEs for all deuterated isotopomers of HPS, BtTPS and HPDMCb in solid phase.

Deuteration influences mostly the C–H bond stretching and C–H in-plane or out-of-plane bending vibrations. For complex polyatomic molecules, these vibrations are naturally mixed with other vibrations such as CCC in-plane bending, CC stretching or ring out-of-plane deformation or twisting vibration (see Tables S13 and S17†). The deuteration effect becomes non-trivial only when the corresponding relaxation energy is important. The vibration types of significant modes with major relaxation energies (see Fig. S2 and S8†) are detailed in Fig. S4, S9 and S10 and Tables S13 and S23.† E.g. for BtTPS, these vibrations are mainly from 2,5-rings, with minor contributions from 1,1,3,4-rings. To elaborate this more clearly, we project the total relaxation energy onto the geometry relaxation in internal coordinates. The contributions from the internal coordinates of the rings at the 1,1-positions or 1,6-positions (1,1), 2,5-positions (2,5) and 3,4-positions (3,4) to the total relaxation energy are depicted in Fig. 6 and listed in Tables S24–26.† It is clearly seen that the substitution positions contributing largely to the relaxation energy could induce remarkable deuteration effects, i.e., for HPS, D-3,4 and D-2,5 have obvious effects; for BtTPS, D-2,5 stands out; for HPDMCb, D-1,1 and D-3,4 have remarkable effects.

Fig. 6

Position dependent relaxation energy for solid-phase HPS, BtTPS and HPDMCb.

Conclusions

In summary, we investigated the isotope effect on the excited-state nonradiative decay rate for a series of AIEgens in both solution and solid phases in comparison with non-AIEgens through first-principles calculations. We found that deuteration always leads to a negative IE on the nonradiative decay rate in both solution and solid phases as expected. Interestingly, the IE is inconspicuous (ca. –10%) for AIEgens in solution in sharp contrast to that of the solid phase of ca. –65% to –95%. The deuteration effect has been observed in the experiments of a typical AIEgen, HPS, in solution and solid states. For non-AIEgens, the IE is comparable (ca. –40% to –90%) in solution and in solid states. The minor IE for AIEgens in solution stems from the competitive results of the negative effect induced by an increased vibronic coupling strength characterized by the Huang–Rhys factor and a positive effect caused by an enhanced DRE. The DRE mainly occurs for low-frequency rotating or twisting vibrations and it becomes weak for AIEgens in solid state compared to in solution. Therefore the IEs of AIEgens in solid phases are similar to those of non-AIEgens in solution and solid states because the nonradiative decay rates in these cases are mainly controlled by the vibronic coupling. Thus, isotopic substitution is shown to be a valuable tool to elucidate the nonradiative decay mechanism induced by aggregation. Partial deuteration effects are estimated for the 6-ring AIE analogues to evaluate the positional dependence. The positions involved in the vibrations with large relaxation energy are found to be responsible for the dependence. Agreement between the theoretical calculations and experimental measurements validates the RIM mechanism for AIE phenomena.

Finally, we should note that computational studies of excited-state dynamics are full of challenges.28 The present methodology focuses only on the intramolecular process by assuming a simple displaced and distorted harmonic oscillator model. As has been shown previously,26 the vibrational quanta decrease steadily with the molecular size. At the same time, the electronic adiabatic transition energy also decreases with size, implying a reduced anharmonic effect for a large system. For further developments by considering the anharmonicity, and excitonic effect as well as charge transfer delocalization, non-perturbation, etc. on the spectroscopy and nonradiative decay rates of multichromophoric aggregates, there is still a long way to go toward the quantitative prediction of light-emitting properties based on first-principles.

Experimental

Ab initio calculations

The solution-phase geometry optimizations and frequency calculations were performed through the equilibrium solvation method. Vertical transition energies at the optimized geometries were obtained via the non-equilibrium solvation method. Solid-phase geometry optimizations were performed on the QM region while the MM region was kept frozen. The centroid molecule of the cluster built from the X-ray crystal structure was chosen as the QM region, whereas the remaining molecules were treated as the MM region. The excitonic effect as well as the intermolecular charge transfer was ignored, since AIEgens in general present twisted structures with large intermolecular distances (more than 5 Å) without π–π interactions.16c,26c,26d For non-AIEgens, flat molecular planes tend to favor π–π stacking interactions, where the radiative process could be strongly affected by forming H- or J-aggregates or excimers. Moreover, the Coulomb coupling components of the excitonic couplings for the representative HPS and BPS were calculated (see Part VIII of the ESI†). The intermolecular excitonic couplings for HPS and BPS are at least one order of magnitude smaller than the intramolecular relaxation energy (Fig. S11 and S12 and Tables S27 and S28†), implying that the photophysical process is dominated by intramolecular motions. No symmetry constraint was imposed during the optimizations in both solution and solid phases. The absence of imaginary frequencies was carefully checked.

Materials and instruments

The absolute fluorescence quantum yields were determined by a calibrated integrating sphere of a Hamamatsu absolute PL quantum yield spectrometer C11347 Quantaurus-QY, excitation wavelength: 365 nm. Spectrophotometric-grade solvents were used in the measurements without further purification.

Click here for additional data file.