Theoretical perspective on fine-tuning ISC efficiency in D-A BODIPY photosensitizers through heteroatom incorporation and nitrogen ring strain modulation

Jianfang Cao Xue Ma Xinyu Chen Tianci Zhang Wen Sun

Citation:  Jianfang Cao, Xue Ma, Xinyu Chen, Tianci Zhang, Wen Sun. Theoretical perspective on fine-tuning ISC efficiency in D-A BODIPY photosensitizers through heteroatom incorporation and nitrogen ring strain modulation[J]. Chinese Chemical Letters, 2025, 36(10): 110897. doi: 10.1016/j.cclet.2025.110897 shu

Theoretical perspective on fine-tuning ISC efficiency in D-A BODIPY photosensitizers through heteroatom incorporation and nitrogen ring strain modulation

English

  • The standard treatments for cancer often come with significant side effects and risks, particularly from surgical procedures, and they typically show limited effectiveness against large tumors [1,2]. As a result, photodynamic therapy (PDT) has recently emerged as a promising alternative for cancer treatment, attracting increased attention [36]. PDT operates on the principle of administering a photosensitizing agent, either intravenously or topically, followed by irradiating the cancerous tissue with specific wavelengths of light. This initiates a photoexcitation process where the photosensitizer transitions from a ground state (S0) to an excited singlet state (Sn). This is followed by intersystem crossing (ISC) to an excited triplet state (Tn), leading to reactions that produce reactive oxygen species (ROS) to destroy the tumor. Additionally, the transition back to the S0 state involves fluorescence emission, suggesting a trade-off between ROS generation and fluorescence intensity in photosensitizers [7].

    A key challenge in developing photosensitizers is optimizing the balance between fluorescence quantum yield (Φf) and singlet oxygen quantum yield (ΦΔ). Donor-acceptor (D-A) type photosensitizers [5,811] manage this by tuning the electron-donating and accepting properties of their components. In D-A systems without heavy atoms [12,13], the ISC often involves a charge transfer (CT) state, leading to the 1CT → 3LE (locally excited) process [14,15]. In earlier research [1619], the identification of charge transfer recombination in heavy atom-free D-A type distyryl-BODIPY photosensitizers was primarily established by altering solvent polarity and confirmed through femtosecond and nanosecond transient absorption spectroscopy. Unfortunately, the aforementioned studies do not provide a clear microscopic description of the charge transfer-recombination process from the 1CT to the 3LE.

    To bridge the existing knowledge gaps regarding charge transfer-recombination in photosensitizers, our research centered on a systematic study of experimentally synthesized M1-M3 [18,19] in Scheme 1. These molecules display maximum absorption and emission in the near-infrared range when dissolved in toluene. M3, modified with a phenothiazine electron-donating group attached to distyryl-BODIPY, exhibits increased Φf and decreased ΦΔ. In contrast, M2, which introduces a fluorine group at the para position of the benzene ring in M3, boasts enhanced electron-withdrawing capacity and significantly higher ΦΔ. M1 with a benzene donor group, serves to contrast the electron-donating strength less robust than phenothiazine, underpinning the pivotal role of donor characteristics in modulating singlet oxygen production. Donors with multiple heteroatoms like oxygen, sulfur, and nitrogen demonstrate enhanced electron-donating capabilities [5,10,1719]. We suggest that the lone-pair electron on these heteroatoms can integrate into the donor's π-conjugated system, boosting electron donation. Based on this, we designed molecules M5-M7 with donor groups containing both thiophene [20] and nitrogen rings [21]. Adjusting the sizes of the nitrogen rings, we modulated the donors' strain energy [22], fine-tuning their oxidation potentials and electron-donating properties, which is crucial for reducing the ΔEST and improving ISC efficiency. We also introduced a comparative molecule, M4, combining triphenylamine as the donor with distyryl-BODIPY as the acceptor [2325]. This comparison will enhance our understanding of sulfur's role in modulating excited states. To further probe the ISC process in D-A type photosensitizers transitioning from 1CT to 3LE, we applied computational chemistry (detailed in Supporting information) to analyze geometric configurations, oxidation potentials, strain energy, energy level diagrams, and electron-hole orbital interactions for molecules M1-M7.

    Scheme 1

    Scheme 1.  Chemical structure and photophysical properties of reported M1-M3 [18,19] and predicted M4-M7 in this work (acceptor in blue, donor in yellow).

    To investigate the impact of conformation on the photophysical behavior of the seven molecules, we analyzed their major bond lengths, bond angles, and dihedral angles (Table S1 in Supporting information). In the S0, S1, and T1 states, the dihedral angles between the donor and acceptor (C1-C2-C3-C4) for molecules M1 through M7, illustrated in Figs. S1-S7 (Supporting information), consistently measure around 90°. This consistency suggests minimal conformational variation in these states. However, a significant divergence is observed in the S2 state, where molecules M2-M4 exhibit a twisted configuration, contrasting with the orthorhombic configuration of M5-M7. Subsequent computational analysis identifies the S2 state as a charge-transfer state, highlighting the significant impact of conformation on photophysical properties. The orthorhombic configuration of M5-M7 in S2 state enhances charge transfer efficiency relative to the twisted configuration of M2-M4. This discrepancy affects the π-conjugation, influencing molecular orbital angular momentum and necessitating adjustments in spin angular momentum to preserve total angular momentum conservation [26]. These adjustments are crucial as they support ISC, a key process for efficient transitions between singlet and triplet states and for enhancing singlet oxygen generation.

    To assess the impact of nitrogen-containing aliphatic rings on the electron-donating abilities of molecules M5-M7, we first calculated their strain energies using a homodesmotic reaction method [22], where higher strain energy indicates lower stability, with M5 being the least stable (Table S2 in Supporting information). Next, we evaluated their oxidation potentials with Eq. S2 (Supporting information). Data from Table S3 (in Supporting information) shows that higher oxidation potentials (EOX) suggest weaker electron donation [27], following the order: M1 > M4 > M2 = M3 > M6 > M7 > M5. Combining the results from the strain energy and oxidation potential calculations, we deduce that M5 possesses the strongest electron-donating capability among molecules M5, M6, and M7.

    This suggests that a lower number of heteroatoms is associated with weaker electron donation. Molecules M5-M7, having lower oxidation potentials compared to M2 and M3, exhibit stronger electron-donating properties, attributed to their thiophene and nitrogen-containing aliphatic rings. Additionally, the near-vertical configuration of M5-M7 in the S2 state works synergistically with the enhanced electron donation to promote charge separation in the CT states, facilitating ISC through the 1CT → 3LE pathway. Thus, the observed consistency between strain energy and oxidation potentials confirms that precise control of strain energy is crucial for fine-tuning the ISC rate.

    To probe the absorption and emission behavior, we computed the absorption and emission wavelengths and frontier orbitals of molecules M1-M7 (Tables S4 and S5 in Supporting information). Discrepancies between calculated and experimental values are due to theoretical approximations, but trends are consistent with experimental data [28]. The molecular orbitals of M1-M7 (Fig. 1, Figs. S8-S10 in Supporting information) show that HOMO → LUMO and HOMO-1 → LUMO transitions correspond to S0 → S1 and S0 → S2 transitions, respectively, with oscillator strengths over 0.01, indicating allowed transitions. The significant HOMO-LUMO overlap suggests enhanced radiative decay from S1 to S0, potentially resulting in fluorescence emission from all M1-M7 [29]. M1 displays the highest oscillator strength (f = 1.1282) among the molecules. According to Eq. S3 (Supporting information), we calculated the radiation transition rate coefficients (fluorescence) kr for seven molecules (Table S6 in Supporting information). Among M1-M3, M1 exhibits the highest radiative transition rate coefficient, aligning with its highest experimentally observed fluorescence quantum yield. M2-M7 show kr values near 2.00 × 108 s-1, indicating their potential to emit fluorescence as visualizers.

    Figure 1

    Figure 1.  Frontier orbitals of the optimized S0 state and the S1 state geometries for M5.

    In the hole-electron diagram, the electronic excitation characteristics are clearly depicted [3032]. The Dhe value, which measures the distance of the center-of-mass between the hole and the electron, and the t value (t is stated in more detail in Eqs. S4-S6 in Supporting information), which measures the degree of separation between the hole and the electron, are both critical parameters. Electron density difference plots are essential for analyzing the electron migration trajectory between two electronic states. As shown in Fig. S11 (Supporting information), the S1 state of all molecules have small Dhe value and negative t value, the holes and electrons are distributed on the acceptor moiety which indicates that the S1 state of M1-M7 is a 1LE state with no significant change in the distribution region before and after the electronic excitation. As seen in Fig. S12 (Supporting information), the S2 state of M1 has both hole and electron delocalized domains throughout the molecule one can categorize the S2 state as a 1LE state which can also be judged by the small Dhe value and the negative t value. We speculate that the inability of the M1 molecule to generate 1O2 (singlet oxygen) is primarily due to its failure to produce a 1CT state and inefficient ISC.

    In our investigation of the electronic excitation characteristics of molecules M2, M3, and M5-M7, we observed that the distances between sulfur (S) atoms and acceptors significantly influence the electronic properties. Specifically, in the S2 state, the S-BODIPY distances in M2 and M3 measure approximately 5.1 Å, whereas in M5-M7, these distances are substantially shorter, around 2.8 Å (Figs. S1-S7 in Supporting information). This reduction in distance in M5-M7 enhances electronic coupling, which likely facilitates more efficient charge transfer and supports the formation of CT excited states. This hypothesis is corroborated by the larger net charge transfer observed between donors and acceptors in M5-M7 compared to M2 and M3, as detailed in Table S7 (Supporting information). Further analysis of the electronic states of molecules M2-M7 revealed that in the S2 state, a large Dhe value and a positive t value (Fig. 2A and Figs. S13A-S17A in Supporting information) indicate that the holes are predominantly located on the donors and the electrons on the acceptors, classifying the S2 state as a 1CT state. Conversely, in the T2 state, Dhe values decrease and t values become significantly negative, suggesting that the T2 state is a 3LE state with holes and electrons primarily located on the acceptor. Significant electron recombination from the S2 to the T2 state is evident in the electron density difference plots (Fig. 2B and Figs. S13B-S17B in Supporting information) of M2-M7, confirming a 1CT → 3LE transition. This transition, involving charge transfer and recombination, could potentially enhance ISC. Notably, M5 exhibits the highest electron-donating capability among these molecules, which is accompanied by the greatest hole-electron separation, as reflected by its largest t value in the S2 state and a maximum net charge transfer of 0.876 between donor and acceptor. M4 contains triphenylamine, which has a smaller oxidation potential than M1 which also reflected in the fact that the net charge transfer between its donor and acceptor is much larger than that of M1, indicating that it is a stronger electron donor than the benzene ring in M1 can induce larger Dhe and positive t values in S2 state. The donor oxidation potential of M4 is slightly higher than that of M2, M3, M5-M7, indicating a weaker electron-donating capacity than them, which is also reflected in the fact that the net charge transfer between donor and acceptor in S2 state of M4 is slightly smaller than that of M2, M3, M5-M7, resulting in marginally smaller Dhe and t than M2, M3, M5-M7. These minor variations are due to the absence of sulfur atoms in M4's donor, but they do not alter the fact that M4 exhibits similar electronic excitation characteristics to M2, M3, and M5-M7.

    Figure 2

    Figure 2.  (A) The hole (blue)-electron (yellow) diagram and (B) the electron density difference between the S2 state and the T2 state diagram of M5 (blue represents a decrease in electron density and yellow represents an increase in electron density).

    To examine how the introduction of S atoms affects the excited states, we have analyzed the excitation types of S2 and T2 states of M2, M3, and M5-M7 in Table S8 (Supporting information) by natural transition orbital (NTOs). The results show that the excitation modes for both the S2 and T2 states of these five molecules are ππ*. When enhancing the value of the occupied NTO isosurface for the S2 state during rendering, the molecular NTO isosurface clearly retains its π-orbital characteristic. In this setup, the n-electrons from sulfur atoms participate in π-conjugation with the donors in both the S2 and T2 states, without leading to n → π* excitation modes. To characterize that the addition of heteroatoms can expand the π-conjugation, we analyze the π-electron distribution of M1-M7 (Fig. S18 in Supporting information), where the green color represents the distribution of π-electrons, and it can be seen that there are π-electrons distributed on the heteroatoms. This facilitates electron transfer between donors and acceptors. In addition to this, S has a different hybridization pattern in thiophene and phenothiazine. The S atom in thiophene forms a conjugated system with the carbon atom through sp2 hybridization, and its lone pair of electrons is highly delocalized from the π system of the entire five-membered ring. The structure of phenothiazine is more complex, with a bridging structure formed between the S and N atoms. We have analyzed the π-electron's net charges of the S atoms in M2, M3, and M5-M7 (Table S9 in Supporting information), and although the S atoms are still involved in the conjugated system of phenothiazines, the π-electron's net charges of the sulfur atom in M2 and M3 are lower than those of thiophene. These results reveal that the introduction of S can increase the degree of π-conjugation of the donor, and when S is involved in the formation of π-bonds in sp2 hybridization, the π-electronic is distributed more on the S, which contributes to the further formation of CT states.

    To investigate the ISC process in compounds M1-M7, we constructed photophysical process diagrams. These diagrams integrate the vertical energies of the S0, S1, S2, T1, and T2 states, SOC values, and ISC rate constants (kisc in Table S10 in Supporting information). Based on the electronic excitation characteristics, we identified that the ISC channel for M2-M7 is primarily S2 → T2, driven by the transition from 1CT to 3LE after excitation to S2. Once in the T2 state, the molecules undergo internal conversion to the T1 state. Furthermore, for effective ISC, photosensitizer molecules should display smaller energy gaps between singlet and triplet states (ΔEST), larger SOC values to enhance the kisc and a large energy gap between T1 and S0 to sensitize 3O2 (triplet oxygen) to 1O2.

    For molecules M2, M3, M5, M6, and M7, the ΔES2-T2 are 0.577, 0.554, 0.193, 0.157, and 0.163 eV respectively, as shown in Figs. 3A and B. Their corresponding spin-orbit coupling (SOC) values are 0.39, 0.41, 0.74, 0.73, and 0.73 cm−1. The intersystem crossing rates (kisc) for M5-M7 are significantly higher than those for M2 and M3, indicating that enhanced donor capacity facilitates ISC. The attachment of an electron-withdrawing fluorine atom to M2′s receptor boosts its electron-accepting ability, leading to a more pronounced electron push-pull effect and a higher kisc value compared to M3. This is consistent with M2′s higher experimentally measured quantum yield of singlet oxygen relative to M3. Overall, M5-M7 are predicted to exhibit stronger ISC capacity and potentially higher singlet oxygen production than M2 and M3. M1, with the largest ΔES2-T2 at 1.156 eV and the SOC value of 0.00 cm-1 (Fig. 3A), cannot undergo ISC, making it suitable only as a fluorescent dye. M4 primarily undergoes ISC through the S2 → T2 channel, with an energy gap of 0.422 eV and an SOC value of 0.02 cm−1 (Fig. S19 in Supporting information). The absence of sulfur in M4, as explained in Eqs. S7 and S8 (Supporting information), leads to lower SOC values and consequently a lower kisc compared to other molecules with multiple heteroatoms, this variation in SOC values among these molecules can be attributed to the dependency of the SOC value on the atomic number, where a higher nuclear charge strengthens the effective magnetic field, thereby intensifying the spin-orbit coupling effect.

    Figure 3

    Figure 3.  (A) Energy level diagram of M1-M3. (B) Energy level diagram of M5-M7.

    In this study, we examined the potential of D-A-type photosensitizers for dual-functional applications, specifically in fluorescence imaging and 1O2 generation. We focused on molecules M2-M7, which incorporate heteroatoms into the donor structure. This modification was found to decrease the oxidation potential and enhance the electron-donating ability, thereby facilitating efficient ISC. Additionally, the introduction of sulfur atoms into these molecules enhances the π conjugation of the donor, which not only provides more electrons but also strengthens the spin-orbit coupling. Electron density difference and hole-electron maps further confirm that ISC in M2-M7 is governed by the 1CT → 3LE mechanism via the S2 → T2 transition. The effectiveness of these molecules in fluorescence imaging is supported by the high overlap between the HOMO and LUMO in the S1 state and significant kr, indicating robust fluorescence capabilities across all these molecules. As a result, these molecules are effective as both therapeutic agents and visualization tools. Whereas M1 is limited to imaging due to its inability to undergo ISC. Specifically, in molecules M5-M7, we fine-tuned the ISC process by adjusting the strain energy in nitrogen-containing aliphatic rings attached to the thiophene unit. These modifications resulted in lower oxidation potentials compared to M1-M4, suggesting a stronger electron-donating capacity. Consequently, M5-M7 exhibit smaller ΔES2-T2, larger SOC constants, and higher kisc, will lead to more efficient triplet state formation and enhanced 1O2 generation, with M5 expected to be the most effective singlet oxygen generator. These findings provide valuable insights for designing dual-functional D-A-type photosensitizers tailored for therapeutic and imaging applications.

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Jianfang Cao: Writing – review & editing, Visualization, Project administration, Methodology, Funding acquisition, Conceptualization. Xue Ma: Writing – original draft, Visualization, Validation, Methodology, Formal analysis, Data curation. Xinyu Chen: Visualization, Methodology. Tianci Zhang: Visualization, Methodology. Wen Sun: Writing – review & editing, Conceptualization.

    This work was financially supported by the Fundamental Research Funds for the Central Universities (No. DUT20RC(3)076) and Natural Science Foundation of Liaoning Province (No. 2020-MS-293)

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2025.110897.


    1. [1]

      N.R. Datta, S. Gómez Ordóñez, U.S. Gaipl, et al., Cancer Treat. Rev. 41 (2015) 742-753. doi: 10.1016/j.ctrv.2015.05.009

    2. [2]

      J. Wang, Q.B. Gong, L.J. Jiao, E.H. Hao, Coord. Chem. Rev. 496 (2023) 215367. doi: 10.1016/j.ccr.2023.215367

    3. [3]

      W.Q. Jiang, M.K. Liang, Q.F. Lei, G.Z. Li, S. Wu, Cancers 15 (2023) 585. doi: 10.3390/cancers15030585

    4. [4]

      K.X. Teng, W.K. Chen, L.Y. Niu, et al., Angew. Chem. Int. Ed. 60 (2021) 19912-19920. doi: 10.1002/anie.202106748

    5. [5]

      K.K. Wen, H. Tan, Q. Peng, et al., Adv. Mater. 34 (2022) 2108146. doi: 10.1002/adma.202108146

    6. [6]

      R. Wang, M.M. He, Z.W. Zhang, et al., Smart Mol. 2 (2024) e20240010. doi: 10.1002/smo.20240010

    7. [7]

      J.F. Cao, T.C. Zhang, W. Sun, Dyes Pigm. 208 (2023) 110797. doi: 10.1016/j.dyepig.2022.110797

    8. [8]

      H. Gu, W. Sun, J.J. Du, J.L. Fan, X.J. Peng, Smart Mol. 2 (2024) e20230014. doi: 10.1002/smo.20230014

    9. [9]

      B. Lu, J.C. Xia, H. Quan, et al., Small 20 (2024) 2307664. doi: 10.1002/smll.202307664

    10. [10]

      V.N. Nguyen, Y.B. Yim, S.G. Kim, et al., Angew. Chem. Int. Ed. 59 (2020) 8957-8962. doi: 10.1002/anie.202002843

    11. [11]

      L.W. Zhang, S.D. Liu, X. Zhang, J.Y. Sun, L.X. Chen, Smart Mol. 2 (2024) e20230022. doi: 10.1002/smo.20230022

    12. [12]

      Y.Q. Hu, Y.Q. Hou, Z.J. Wang, Y.Q. Li, J.Z. Zhao, J. Chem. Phys. 153 (2020) 224304. doi: 10.1063/5.0025224

    13. [13]

      Y.T. Gao, X.X. Wang, X.W. He, et al., Adv. Funct. Mater. 29 (2019) 1902673. doi: 10.1002/adfm.201902673

    14. [14]

      L.F. Wang, Y. Qian, Biomater. Sci. 11 (2023) 1459-1469. doi: 10.1039/d2bm01838a

    15. [15]

      Y.Q. Hou, X. Zhang, K.P. Chen, et al., J. Mater. Chem. C 7 (2019) 12048-12074. doi: 10.1039/c9tc04285g

    16. [16]

      J. Deckers, T. Cardeynaels, H. Penxten, et al., Chem. Eur. J. 26 (2020) 15212-15225. doi: 10.1002/chem.202002549

    17. [17]

      J. Deckers, T. Cardeynaels, S. Doria, et al., J. Mater. Chem. C 10 (2022) 9344-9355. doi: 10.1039/d2tc01526a

    18. [18]

      Y.Q. Hou, Q.Y. Liu, J.Z. Zhao, Chem. Commun. 56 (2020) 1721-1724. doi: 10.1039/c9cc09058d

    19. [19]

      Y. Dong, A. Elmali, J.Z. Zhao, B. Dick, A. Karatay, ChemPhysChem 21 (2020) 1388-1401. doi: 10.1002/cphc.202000300

    20. [20]

      J.F. Cao, X.Y. Chen, X. Ma, T.C. Zhang, W. Sun, Phys. Chem. Chem. Phys. 26 (2024) 21520-21529. doi: 10.1039/d4cp01346h

    21. [21]

      L.J. Ma, Y.J. Geng, G.Y. Zhang, et al., Adv. Health. Mater. 12 (2023) 2300733. doi: 10.1002/adhm.202300733

    22. [22]

      T. Lu, Z.Y. Liu, Q.X. Chen, Chin. Phys. B 31 (2022) 126101. doi: 10.1088/1674-1056/ac873a

    23. [23]

      S.R. Zhou, R.P. Li, Y. Li, Y.X. Wang, L.H. Feng, J. Mater. Chem. B 10 (2022) 8003-8012. doi: 10.1039/d2tb01578a

    24. [24]

      Y. Chen, S. Zhang, T.J. Cheng, et al., J. Mater. Sci. Technol. 141 (2023) 135-148. doi: 10.3847/1538-4357/acd753

    25. [25]

      P.H. Xiao, Z.P. Shen, D.L. Wang, et al., Adv. Sci. 9 (2022) 2104079. doi: 10.1002/advs.202104079

    26. [26]

      D. Sasikumar, A.T. John, J. Sunny, M. Hariharan, Chem. Soc. Rev. 49 (2020) 6122-6140. doi: 10.1039/d0cs00484g

    27. [27]

      X.F. Zhang, N. Feng, Chem. Asian J. 12 (2017) 2447-2456. doi: 10.1002/asia.201700794

    28. [28]

      J.F. Cao, T.C. Zhang, X.Y. Chen, X. Ma, J.L. Fan, Smart Mol. 3 (2025) e20240023. doi: 10.1002/smo.20240023

    29. [29]

      K. Shizu, H. Tanaka, M. Uejima, et al., J. Phys. Chem. C 119 (2015) 1291-1297. doi: 10.1021/jp511061t

    30. [30]

      T. Lu, J. Chem. Phys. 161 (2024) 082503. doi: 10.1063/5.0216272

    31. [31]

      Z.Y. Liu, T. Lu, Q.X. Chen, Carbon 165 (2020) 461-467. doi: 10.1016/j.carbon.2020.05.023

    32. [32]

      T. Lu, F.W. Chen, J. Comput. Chem. 33 (2012) 580-592. doi: 10.1002/jcc.22885

  • Scheme 1  Chemical structure and photophysical properties of reported M1-M3 [18,19] and predicted M4-M7 in this work (acceptor in blue, donor in yellow).

    Figure 1  Frontier orbitals of the optimized S0 state and the S1 state geometries for M5.

    Figure 2  (A) The hole (blue)-electron (yellow) diagram and (B) the electron density difference between the S2 state and the T2 state diagram of M5 (blue represents a decrease in electron density and yellow represents an increase in electron density).

    Figure 3  (A) Energy level diagram of M1-M3. (B) Energy level diagram of M5-M7.

  • 加载中
计量
  • PDF下载量:  0
  • 文章访问数:  89
  • HTML全文浏览量:  0
文章相关
  • 发布日期:  2025-10-15
  • 收稿日期:  2024-11-10
  • 接受日期:  2025-01-20
  • 修回日期:  2025-01-12
  • 网络出版日期:  2025-01-21
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

/

返回文章