Chun-Lin
Sun
a,
Qing
Liao
b,
Ting
Li
c,
Jun
Li
a,
Jian-Qiao
Jiang
a,
Zhen-Zhen
Xu
b,
Xue-Dong
Wang
b,
Rong
Shen
d,
De-Cheng
Bai
d,
Qiang
Wang
a,
Sheng-Xiang
Zhang
c,
Hong-Bing
Fu
be and
Hao-Li
Zhang
*a
aState Key Laboratory of Applied Organic Chemistry (SKLAOC), College of Chemistry and Chemical Engineering, Lanzhou University, Lanzhou 73000, P. R. China. E-mail: haoli.zhang@lzu.edu.cn
bBeijing National Laboratory for Molecular Sciences (BNLMS) Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
cSchool of Life Sciences, Lanzhou University, Lanzhou 73000, P. R. China
dSchool of Basic Medical Sciences, Lanzhou University, Lanzhou 730000, P. R. China
eDepartment of Chemistry, Capital Normal University, Beijing 100048, P. R. China
First published on 7th October 2014
Small organic dyes with large two-photon absorption (TPA) cross sections (δ) are more desirable in many applications compared with large molecules. Herein, we proposed a facile theoretical method for the fast screening of small organic molecules as potential TPA dyes. This method is based on a theoretical analysis to the natural transition orbitals (NTOs) directly associated with the TPA transition. Experimental results on the small indolic squaraine dyes (ISD) confirmed that their TPA cross sections is strongly correlated to the delocalization degree of the NTOs of the S2 excited states. Aided by this simple and intuitive method, we have successfully designed and synthesized a small indolic squaraine dye (ISD) with a remarkable δ value above 8000 GM at 780 nm. The ISD dye also exhibits a high singlet oxygen generation quantum yield about 0.90. The rationally designed TPA dye was successfully applied in both two-photon excited fluorescence cell imaging and in vivo cerebrovascular blood fluid tracing.
Previous design of TPA dyes has generally been based on using extended π-conjugated frameworks with enforced coplanarity and addition of donor and acceptor groups.7 By using these strategies, most of the reported TPA dyes with high δ consist of large π-systems and various donor–acceptor units.8 However, large conjugated molecules are not always appreciated in practice because they may meet difficulties in synthesis, purification and solubility.9 Importantly, large conjugated chromophores tend to have small S1–S0 energy gaps, and rapid S1–S0 nonradiative internal conversion, resulting in low-energy short-lived excited states, and low quantum yield, which are undesirable for fluorescence or photochemical activations and photodynamic therapy (PDT).7a,10 For applications such as TPA imaging, photoinitiator and optical limiting, ideal TPA dyes should exhibit peak absorption matching the commercial laser source.11 Only in a few special applications, such as TPA at telecommunication wavelengths, a narrow optical gap is preferred.12 Therefore, maximizing the δ is not the only consideration when designing a TPA dye. Besides, it is often useful to pack the maximum effect into the smallest possible chromosphere. Anderson has proposed that δ/Ne (where Ne is the number of π electrons) can be used as a useful parameter to evaluate how efficiently the π-conjugated system contributes to the δ.7a To date, there has been very few systems achieving δ/Ne > 200 GM, and these reported molecules tend to be either too complicated for synthesis13 or have too low energy gap.8a,14 In TPA microscopy, small TPA dyes that have peak TPA absorption in the 700–900 nm window with high δ and are most favorable,15 while an ideal TPA dye should absorb at the shortest possible wavelengths to enable the best possible resolution.16
Developing theoretical methods to guide the design of smaller but more efficient TPA dyes is an attractive idea.8b,16 Because the TPA process is a third-order nonlinear phenomenon and δ is related to the imaginary part of the third nonlinear polarization, theoretical calculation of δ is very complicated. Brédas et al. have studied the relationship of δ with single photon absorption transition dipole moment using semi-empirical AM1 and INDO calculations.17 Full ab initio calculation of the third nonlinear polarization has been studied by response theory, including the sum-over-states technique18 and the finite field method.19 Unfortunately, even with the state-of-the-art calculation tools, the theoretically predicted δ often showed large discrepancy with the experimental results, in some cases, by orders of magnitude.7b,20 Another problem with these complicated theoretical methods is that the structure–property relationship is hidden by the complex wave function calculation, and can hardly provide any intuitive information for molecular design. To date, new TPA materials are typically discovered experimentally, with theoretical characterization providing mostly post facto justification.21 Therefore, an easy-to-use theoretical method that can guide the practicable design of TPA materials is critically needed.
In this work, we report a simple theory-assisted method for the fast screening of organic molecules as potential TPA dyes. This method is based on morphology-analysis of the final excited state of the TPA process, which provides qualitative information that is useful for both pre-synthesis predication and post facto analysis. By applying this method on a series of indolic squaraine dye (ISD), we have successfully obtained a small ISD molecule with an exceptionally large δ value. The peak TPA absorption of the obtained dye is around 780 nm, which matches well with the laser source of the TPA microscope. The newly designed TPA dye shows potential in both in vitro and in vivo bio-imaging applications.
To unveil the structure–property relation, we have designed seven small ISD molecules, as shown in Fig. 1. Two positions were readily used in structural modification, i.e., the 5-position of benzene ring and the N-position of the pyrrole ring. The seven molecules can be classified into two groups. From ISD-1 to ISD-4, the maximum number of π electron in conjugation is maintained at 22. The Br, n-hexyl and benzyl substitutions to the molecules are designed to study the effects of heavy atoms, non-conjugated alkyl and aromatic substitutions, respectively. ISD-1 is very similar to that reported in previous literature,25 and is used as a reference. From ISD-5 to ISD-7, the π system is expanded by directly attaching aromatic structures to the indole part (Table 1). Some other designed molecules are presented in the ESI† (part 6).
ISD-1 | ISD-2 | ISD-3 | ISD-4 | ISD-5 | ISD-6 | ISD-7 | |
---|---|---|---|---|---|---|---|
a N e are the maximum number of conjugated π electrons in the planar part with the total number of π electrons in the molecules given in parentheses. b λ abs (nm) is the maximum absorption wavelength. c λ em (nm) is the maximum emission wavelength. d τ (ns) is the fluorescence lifetime. e Φ is the fluorescence quantum yield. f k r (ns−1) is the rate constant for radiative deactivation from S1 to S0. g k nr (ns−1) is the rate constant for non-radiative deactivation. h δ (GM) is the TPA cross section. i Φ Δ is the singlet oxygen quantum yield. j δ/Ne is δmax normalized by the π electron number. k δΦ Δ presents the TPA singlet oxygen generation capability. | |||||||
N e | 22 (22) | 22 (22) | 22 (34) | 22 (34) | 38 (42) | 38 (54) | 22 (34) |
λ abs | 560 | 561 | 562 | 561 | 561 | 578 | 575 |
λ em | 580 | 581 | 582 | 582 | 582 | 598 | 600 |
τ | 1.58 | 1.65 | 1.53 | 1.57 | 1.64 | 2.67 | 1.47 |
Φ | 0.701 | 0.499 | 0.505 | 0.301 | 0.676 | 0.580 | 0.288 |
k r | 0.444 | 0.302 | 0.330 | 0.192 | 0.412 | 0.218 | 0.196 |
k nr | 0.189 | 0.304 | 0.324 | 0.445 | 0.198 | 0.157 | 0.484 |
δ max | 330 | 1260 | 730 | 701 | 712 | 645 | 8019 |
Φ Δ | 0.30 | 0.62 | 0.50 | 0.77 | 0.70 | 0.64 | 0.90 |
δ/Nej | 15 (15) | 57 (57) | 33 (21) | 32 (21) | 19 (17) | 17 (12) | 365 (236) |
δΦ Δ | 99.7 | 785.8 | 365.1 | 543.0 | 501.9 | 411.4 | 7217.1 |
Previous studies have shown that using an extended π-conjugated framework and donor–acceptor groups are likely to enhance the δ.7 Meanwhile, it is important to be aware that the extended conjugation structure or added groups need to participate the TPA process in order to contribute to the δ. In a TPA process, the ground state molecule simultaneously absorbs two photons to be excited to a excited state through a virtual state (Fig. 2). Both the virtual state and the final state should have strong influence on the δ. Accurate description of TPA excited state need to include electron-correlation effects known as configuration interaction (CI)26 and thus requires complex and time-consuming calculations. Moreover, a precise calculation of the virtual state is difficult to achieve at the moment. To simplify the calculation, we focused on the final excited state of the TPA process. The electronic properties of the excited state can be derived by the calculation of natural transition orbitals (NTOs),27 which can be used to evaluate the participation or influence of a certain structural unit in the molecule. It is known that a high degree of delocalization of the excited states may enlarge the oscillator strength of a particular absorption and increase the corresponding transition moment.7a In the TPA process, a highly delocalized final excited state is more likely to give a large δ than a localized one. Herein, we demonstrate that the visualization of the NTOs not only can help to understand the excited state associated to the TPA process but also help to screen TPA materials with high δ.
It is known that the TPA process depends on the symmetry of the chromophore, and the final state is different for noncentrosymmetric or centrosymmetric chromophores.7a For dipolar noncentrosymmetric chromophores, their one- and two-photon transitions are the same, which implies that the maximum TPA transition energy is half of the energy needed in the one-photon counterpart. However, the one and two photon transitions are different for the symmetric quadrupolar chromophores, in which the TPA transition is from the ground state (S0, 1Ag) to the lowest excited state with Ag symmetry (S2, 2Ag), so that the S2 state is the final state in TPA process.7 Because the ISD compounds designed in this work are all centrosymmetric, we studied their NTOs of the S2 states.
Fig. 2 illustrates the electron transition of centrosymmetric chromophores in a TPA process. As illustrated in Fig. 2, the final TPA excited state (S2 state) of the ISD molecules can be described by the particle–hole pair of the NTO-2. To study the excited states, TD-DFT calculations were conducted using the Gaussian 09 package.28 The detailed simulation results are collected in the ESI (Fig. S1†). Structural optimization suggests that the indolic squaraine structure, including the four-membered ring and two indolic units, has a coplanar conformation, which is ideal for an extended excited state. The planarity of ISD-5 and ISD-6 extends to the phenylethynyl groups. The NTOs of S1 and S2 were studied to visualize the excited state in the form of “electron” and “hole” distributions. We paid particular attention to the delocalization of the NTOs of S2, as S2 is the final state after absorbing two photons. The calculated “holes” of S2 states of the representative ISDs are shown in Fig. 3, while the other NTOs are provided in Fig. S2.†
![]() | ||
Fig. 3 The natural transition orbitals for the “hole” in the second excited state (S2) of seven ISD molecules. |
Fig. 3 reveals that the excited states are distinctly different for the ISD molecules, despite their structural similarity. Importantly, not all parts of the conjugated structures take part in the TPA final states. For ISD-1, which is the alkyl substituted ISD, the “hole” of the NTO-2 is mainly restricted to the four-membered squaraine ring, and the benzene rings do not contribute to the excited state. When the hexyl units are replaced by the benzyl groups (such as ISD-3), the situation is very different. The excited sate of ISD-3 is mainly localized on the two indole rings. However, the four-membered ring does not participate in the excited state, so that the NTO is divided into two separated parts. The Br atoms on the indole rings are involved in the NTO-2 but they do not change the NTO morphology very much, and the NTOs of ISD-2 and ISD-4 are very similar to that of ISD-3. Expanding the π system was achieved by attaching phenylethynyl moieties to the 5-position of the indole benzene (ISD-5 and ISD-6). Seen from the “hole” of S2, the orbital of the excited state indeed extends to contain both the indole and the phenylethynyl units. However, the NTOs are still separated by the four-member ring, which gives limited delocalization of the TPA excited state. The most different NTO is observed on ISD-7, which has a benzene group attached on the N-position of each indole ring. It is very interesting that the “hole” of ISD-7 extended to cover both the squaraine core and the two indole units, though the two benzene rings attached on the N positions are not involved in the S2 state due to non-planar conformation.
As mentioned above, the S2 state is the final excited state of the TPA process. A reasonable presumption is that a high degree of delocalization of the final excited state is beneficial for the corresponding absorption process. The variations of NTO of different molecules are associated to their energy levels and orbital symmetry, and cannot be simply predicted from their molecular structures. Our NTO analysis helps to visualize the excited state and hence provides useful information to qualitatively predict their δ values. For compound ISD-1, its NTO of the S2 state is very small and locates only in the four-membered ring of the squaraine core. Therefore, it is not surprising that this type of molecules has relatively small δ, as reported in previous studies.25 In the case of ISD-2 to ISD-6, their NTO of the excited states are larger than that of the ISD-1, so that larger δ values are expected. However, their NTOs are separated into two parts, which is not the most ideal situation for TPA process. In contrast to all the other molecules, the calculation shows that the ISD-7 have the most delocalized excited state, and hence is expected to have the largest TPA cross section.
To test the above hypothesis, the compounds ISD-1 to ISD-7 were synthesized and their absorption and emission properties have been characterized. The experimental one-photon absorption and emission spectra of the seven ISD molecules in diluted THF solution are shown in Fig. 4a and b, respectively; while the key parameters are summarized in Table 1. The one-photon absorption spectra have also been simulated by time-dependent density functional theory (TD-DFT) methods (Fig. 4c), which helps to give precise assignments to the absorption peaks.
![]() | ||
Fig. 4 (a) Steady-state one-photon absorption and (b) fluorescence of ISD molecules in THF, and (c) the simulated absorption spectra. |
The one-photon absorption spectrum of the ISD-1 shows a typical strong and sharp absorption band around 560 nm, corresponding to the S0–S1 transition. The absorption bands of ISD-2 to ISD-4 are nearly unaffected by the substitution. The main absorption bands of ISD-6 and ISD-7 show bathochromic-shifts for 18 and 15 nm, respectively. Experimental results are generally consistent with the calculated results, although certain deviations exist, which is common in TD-DFT calculation.29 For instance, the transition energy of S0–S1 is underestimated for ISD-5 and ISD-6 because of the presence of the C–C triple bond.30
The one-photon absorption spectra also indicate that the S0–S2 transition intensity (shown in the Fig. 4a) is much weaker than that of the S0–S1 transition. From the calculated results, the oscillator strengths of the S0–S2 transitions are all zero for the seven ISDs, which can be attributed to the symmetry forbidden one-photon absorption process. In practice, the S0–S2 transition intensity comes from the vibronic coupling (also called Herzberg–Teller effect) by which the symmetry-forbidden rule is broken.31 The selection rule is reversed between two-photon and one-photon absorption processes in general to centrosymmetric chromophores,7 and the S0–S2 transitions are in favor of the TPA process for our ISDs. The emission spectra of ISD-1 to ISD-5 are nearly the same, while the emissions of ISD-6 and ISD-7 are slightly red-shifted, consistent with their absorption spectra (Fig. 4b). Fig. S2† shows that the colors of ISD-1 to 7 are very similar with each other under the one-photon absorption condition, and are almost the same in their fluorescence emission. The similar color indicates very small differences exist in the one-photon absorption processes among these seven ISDs.
The fluorescent quantum yields (Φ) of these seven ISDs were measured and compared in Table 1. The alkyl substituted ISD-1 gives the highest Φ of 0.701, while introduction of all the other substituent groups tends to reduce the Φ. The Φ of the Br-containing molecules (ISD-2 and ISD-4) were suppressed compared with their non-Br counterparts (ISD-1 and ISD-3) due to the heavy-atom effect,32 which results in reduced radiative rate constant of kr and concomitant enhanced non-radiative deactivation rate of knr, as derived from the time-resolved fluorescence spectra (Fig. S3†). The Φ values of benzyl-containing ISDs (ISD-3, ISD-4 and ISD-6) were lower compared with the corresponding non-benzyl-containing molecules (ISD-1, ISD-2 and ISD-5). The Φ reduction also occurs upon phenylethynyl substitutions. The benzene modification at the N-position brings about the smallest Φ for ISD-7, which is only 0.288. The Φ suppression after attaching different aromatic groups probably arises from the intra-molecular rotation of the substituent groups, which provides additional channels for non-radiative de-excitation.33
The TPA spectra between 730 and 840 nm are shown in Fig. 5, and the maximum TPA δ is shown in Table 1. The maximum TPA absorption occurs around 780 nm for all the studied molecules, corresponding well with the two-photon absorption for the S0–S2 transition. The TPA δ of ISD-1 was measured to be 330 GM under our condition, which is in good agreement with previously reported ISD molecules with similar structure (450 GM).25 For the substituted molecules from ISD-2 to ISD-6, the δ values are in the range of 600–1200 GM, much larger than that of the model molecule ISD-1. This result is consistent with our theoretical calculation results, where ISD-2 to ISD-6 have larger NTOs of the S2 states compared with that of ISD-1. It is very interesting that the TPA δ of ISD-2 is about 1200, three times larger than that of ISD-1. As discussed before, the halogenation effect to the photoluminescence process is generally very complicated.29b,34 From the NTOs of ISD-2, such TPA δ enhancement in ISD-2 can be attributed to the participation of the Br-atom in the TPA process.
It was discovered that the δ of ISD-7 is exceptionally high, reaching 8019 GM, which is more than 20 times higher than that of ISD-1. Indeed, two-photon excitation fluorescence shows that the ISD-7 is significantly brighter than ISD-1 under the same condition (inset of Fig. 5). Such a high δ value is in excellent agreement with our theoretical calculation that ISD-7 has the most extended NTO of the S2 states. It is noted that the δ value of ISD-7 measured under our condition is about 17 times higher than the largest δ value ever reported for indolic squaraine derivatives.16,25 If we only consider the π electrons in the planar part, the δ/Ne is as high as 365 GM, which is almost the highest for all the known organic TPA dyes with maximum absorption below 800 nm.7a,13 The successful design of ISD-7 strongly suggests that the proposed NTO morphology analysis method is an efficient method for the rational design of organic dyes with large TPA cross sections.
The above experimental results confirmed that extension of the π-system is not necessarily a useful way to increase δ.7a For instance, the molecules with the largest number of π electrons in this work, i.e. ISD-5 and ISD-6, did not show very high δ values, indicating that the expansion of the π system does not necessarily enlarge δ. The π-expansion strategy is only effective when the additional π structures participate in the TPA process. We also found that the non-planar substitution, such as the phenyl groups in ISD-7, may facilitate the δ enlargement by extending the NTOs associated with TPA absorption. In addition, we note that one should be cautious to directly use the linear optical properties, such as one-photon absorption, fluorescence spectra, fluorescence lifetime and quantum yield, to design TPA dyes.
The large δ of molecule ISD-7 is remarkable considering the simple molecular structure. Importantly, its TPA absorption peaks at 780 nm, which matches well the efficient output of the Ti:sapphire laser and the optical window in biological tissue.16 Moreover, ISD-7 exhibits low cytotoxicity (Fig. S4†). To examine its luminescence imaging capability in vitro, fibroblast cells were incubated with ISD-7, and then imaged using two-photon laser confocal scanning microscopy (TPLCSM). Intense intracellular luminescence was observed with femtosecond laser pluses at 780 nm (Fig. 6a), under which excitation no background fluorescence from the fibroblast cells would interfere. Benefited from the high resolution of the TPA microscopic image, much detailed intracellular structural information can be obtained. It is found that the luminescence was mainly from the phospholipid bilayer membrane of cytoplasm whereas the nucleus gave no luminescence emission. Though there are very few ISD-7 molecules in the membrane, the membrane structure can be observed in high resolution (inset of Fig. 6a) because of the high δ value.
![]() | ||
Fig. 6 (a) Two-photon laser confocal scanning microscopy (TPLCSM) image of fibroblast cells incubated with ISD-7 upon excitation at 780 nm; the inset shows the high resolution image. (b) TPLCSM image of the blood plasma labelled with the tracer of ISD-7 emulsion in blood vessels; the image shows the flow trajectory by merging the dynamic photo micrographs (the video is included in the ESI, Video-S2†). The blood vessels are indicated by the dashed lines, and the arrows present the flow direction. |
Though NIR excitation is favorable for biological imaging due to the weak absorption by tissue in this wavelength region, the two-photon excited fluorescence (TPEF) intensity diminishes exponentially as the focus deepens into the tissue. Hence fluorescent dyes with large TPA cross section are in great need for in vivo imaging. The fluorescent emulsion of ISD-7 was injected into the mice blood vessel by tail vein for in vivo tracing of the blood flow. First, the fluorescent emulsion was monitored by the one-photon excited fluorescence microscopy, see video in ESI (Video-S1†). Compared with traditional dyes, the bright red emulsion flows smoothly in the brain vessel in the microscopic video, providing unprecedented contrast to observe both the vessel distribution and the inside fluid situation simultaneously. To track the movements of individual emulsion particles, TPLCSM imaging was performed. The tracking is shown in the ESI (Video-S2†). Twelve dynamic TPLCSM photos are merged in Fig. 6b to indicate the blood flowing behavior. In Fig. 6b, the movements of three separated emulsion particles are marked, which are located in one stem vessel and two branch vessels, respectively. From the tracer of the emulsion particles, it was found that the blood fluid in right branch vessel is faster than the left one and the stem vessel. Such data may provide information for understanding the brain blood flow in different vessels. This result demonstrates that the fluorescent emulsion of ISD-7 can be successfully applied for in vivo TPLCSM tracing.
As mentioned above, the designed ISD molecules have a relatively large S1–S0 energy gap, which is beneficial for photochemical activation. One possible application to take advantage of this property is to activate oxygen and promote generation of singlet oxygen species. To verify this proposal, the singlet oxygen quantum yields (ΦΔ) of the ISDs were measured by the singlet oxygen luminescent method, using Rose Bengal (RB) as a standard. The singlet oxygen 1Δg emission spectra are shown in Fig. 7 and the determined ΦΔ values are collected in Table 1. The ΦΔ of ISD-1 is about 0.30, which is normal for this type of dyes. Table 1 shows that all the other ISDs show enhanced ΦΔ compared to ISD-1. The ΦΔ of Br-containing ISD dyes are larger than their non-Br counterparts, presumably because the heavy Br atom promotes excited-state intersystem crossing.35 The ΦΔ of ISD-4 is similar to that of other good photosensitizers such as RB. The ΦΔ of ISD-7 is above 0.90, which is significantly larger than the other molecules, even larger than RB. Such high ΦΔ is very rare among metal-free sensitizers.
Chromophores that have both high δ and ΦΔ are being considered for two-photon photodynamic therapy (PDT) applications. PDT is a promising noninvasive treatment for cancers and other diseases through singlet oxygen produced by the photosensitizer. Two-photon PDT is advantageous over the traditional one-photon counterpart, because of its deeper penetration capability into living tissues and 3D selectivity. For two-photon PDT applications, the product of δΦΔ describes the efficiency of molecules for two photon singlet oxygen generation. As shown Table 1, the δΦΔ of ISD-7 is as high as 7217, which is the largest for known squaraine dyes. This result indicates that ISD-7 has potential application for two-photon PDT. The relevant biological research is currently under way.
For the measurements of the fluorescence time profiles, the time-correlated single-photon counting (TCSPC) method using a nanosecond pulsed LED sources (464 nm, FWHM ca. 800 ps) with 40 MHz repetition rate was employed. A photomultiplier tube and a counting board (PicoQuanta, PicoHarp 300, Germany) were used for signal detection.
It is noted that one should not be confused by the different values of Φ and ΦΔ as they corresponding to different processes and were measured under different conditions. The solution concentration for the ΦΔ measurement was 0.2 mM, which was much higher than that used for the Φ measurement. Besides, all the ΦΔ were measured under air-saturation condition. According to the literature,38 uncertainties for the quantum yields based on relative luminescence method were 5–10%.
For in vivo time-lapse imaging of emulsion flowing, snap (1024 × 1024 pixel arrays) were taken with the emulsion flowing (Video-S2†).
We also realize that there are still some limitations of the NTO analysis methodology. For example, the information from the simulation is qualitative and the relationship between the NTOs and the molecule structure needs more definitive description. Our further endeavors will focus on addressing these issues. After all, it is believed that this new theory-assisted molecular design strategy demonstrated in this work will significantly accelerate the discovery of new organic TPA dyes with even further improved performance.
Footnote |
† Electronic supplementary information (ESI) available: Synthetic procedures and structural characterization of compounds, all the NTOs, time-resolved fluorescence spectra and the video of microscopy (Video-S1 and S2). See DOI: 10.1039/c4sc02165g |
This journal is © The Royal Society of Chemistry 2015 |