Vidisha Rai-Constapel and
Christel M. Marian*
Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitätsstr. 1, D-40225 Düsseldorf, Germany. E-mail: Christel.Marian@hhu.de
First published on 8th February 2016
High-level electronic structure methods and quantum chemistry programs have been employed for a thorough investigation of the photophysics of acridone in isolated and solvated states. A kaleidoscope of photophysical behavior results by varying the medium in which acridone finds itself. With the computed intersystem crossing rate constants of the order of 1011 s−1 in vacuum, the photoexcited species acts as an effective triplet sensitizer. In polar, aprotic media, the radiationless decay processes for the photoexcited singlet state are still predominant, however, a delayed fluorescence may also be observed. Such delayed fluorescence comes into play due to the near degeneracy of the 1(πHπ*L) and 3(nOπ*L) electronic states in such a solvent. As the polarity, as well as the proticity of the solvent increases, the energy shifts experienced by the involved electronic states cause the photoexcited molecule to relax primarily via fluorescence emission to the ground state with a rate constant of the order of about 107 s−1. In such solvents, acridone is a good fluorescence marker. Hence, tuning the polarity and proticity of the solvent should turn acridone into a dark (triplet formation) or a bright (fluorescence) sensitizer.
![]() | ||
Fig. 1 Molecular structure of the xanthone family of molecules: thioxanthone (X![]() ![]() ![]() |
A review article by Zhao et al. mentions the possible use of AC as triplet sensitizer in a triplet–triplet annihilation upconversion scheme based on the fact that intersystem crossing (ISC) from the singlet to the triplet state in photoexcited AC is feasible without the involvement of any heavy atoms, the ISC being promoted via n → π* type transitions.2 The triplet sensitizing properties of AC have also been exploited for exciting Eu3+ in a complex used as a luminescent marker in biological applications.3 AC itself is also applied as a fluorescent probe for DNA labeling.4 The solvent dependent luminescence of AC has been made use of to design a pair of donor–acceptor–donor triads with a potential for application as organic semi-conductors.5
Whether AC acts as a triplet sensitizer or as a fluorescent marker, depends upon the excited-state energy landscape. This landscape, in turn, varies with the type of medium (or the absence of it!) surrounding the molecule. Though experimental studies regarding the various applications of AC are numerous, an extensive theoretical analysis of its photophysics was not found in the literature. Previous theoretical works on the related molecules xanthone (X) and thioxanthone (TX), have established the importance of theory in order to gain a better insight and understanding of the excited-state processes.6–9 This paper extends such an analysis to AC and presents a detailed quantum chemical study of this molecule in vacuum, as well as four solvents, namely, chloroform, acetonitrile (AcN), methanol (MeOH) and water.
The conductor-like screening model (COSMO) was used to take account of the electrostatic interaction of AC with the solvent.16,17 The effect of specific bonding with solvent molecules was taken account of via micro-solvation. To this purpose, explicit solvent molecules were placed at appropriate hydrogen bonding sites of AC and the whole system was enveloped by the implicit solvation model COSMO (see Fig. 2).
Single-point calculations at the optimized structures were carried out using the DFT/MRCI method in order to obtain the vertical excitation energies, dipole moments and oscillator strengths.18 For the isolated molecule, the calculations were carried out in the C2v symmetry. Due to technical reasons, DFT/MRCI calculations involving COSMO had to be carried out in the C1 symmetry. In the DFT/MRCI step, the solvent effect is incorporated solely by adding the solute–solvent interaction terms to the Fock matrix of the reference state, i.e., electronic relaxation effects and the costs required for the polarization of the solvent are not accounted for.19 The ISC rate constants were calculated in the Condon approximation which requires the electronic spin–orbit coupling matrix elements and the vibrational overlaps. The program SPOCK developed in our group was employed to determine the spin–orbit matrix elements (SOMEs) between the correlated DFT/MRCI wavefunctions.20,21 Recently, a time-dependent method for computing ISC rate constants has been developed in our laboratory. This method can also take account of temperature effects on the rate constants if desired.22,23 Computation of the ISC rate constants in this work is based on this program.
A qualitative analysis of the photophysical processes occurring after the photoexcitation has been carried out by following the linearly interpolated path (LIP) between the minima of the electronic states involved. At each of the chosen points along the LIP, a single-point DFT/MRCI calculation was carried out to determine the energy profile of the singlet and triplet excited electronic states.
State | Electronic structurea | ΔE | |||||||
---|---|---|---|---|---|---|---|---|---|
Vacuum | CHCl3 | AcN | MeOH | Water (b) | Water (c) | Water (d) | |||
a The characterization of the electronic states is based on the gas-phase results.b The oscillator strengths are italicized.c AC monomer, band origin.24d AC in AcN, band maximum.25e AC in MeOH, band maximum.25f AC in water, band maximum.25 | |||||||||
S0 | 11A1 | (94.7) ground state | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 |
S1 | 11A2 | (75.4) nO → π*L | 3.46 | 3.67 | 3.75 | 3.97 | 4.11 | 3.85 | 3.87 |
0.00b | |||||||||
S2 | 21A1 | (80.1) πH → π*L | 3.50 (3.38c) | 3.43 | 3.40 (3.31d) | 3.31 (3.27e) | 3.32 (3.23f) | 3.33 | 3.33 |
T1 | 13A1 | (83.7) πH → π*L | 2.94 (2.90c) | 2.85 | 2.80 | 2.65 | 2.68 | 2.71 | 2.69 |
0.00 | |||||||||
T2 | 13A2 | (76.6) nO → π*L | 3.30 | 3.53 | 3.61 | 3.85 | 4.00 | 3.72 | 3.73 |
0.00 | |||||||||
T3 | 13B2 | (31.7) πH → π*L+1, (19.1) πH−1 → π*L | 3.58 | 3.58 | 3.58 | 3.54 | 3.55 | 3.56 | 3.56 |
(12.4) πH−3 → π*L, (11.1) πH−2 → π*L+2 | 0.00 | ||||||||
T4 | 23A1 | (57.8) πH−2 → π*L, (12.0) πH−3 → π*L+2 | 3.59 | 3.61 | 3.61 | 3.64 | 3.63 | 3.62 | 3.62 |
0.01 |
For the explicit solvent model (b) (see Fig. 2), the blue shift of the nOπ*L is enhanced as two H-bonds are formed at the carbonyl oxygen. We observe that the red shift experienced by the πHπ*L states remains pretty much the same irrelevant of the explicit solvation model. The major absorption band corresponding to S0 → S(πHπ*L) transition at 3.23 eV observed by Siegmund et al.25 for AC in aqueous solution agrees well with the calculated vertical excitation energy of 3.32 eV for the bright singlet S(πHπ*L) state at the S0 minimum. As in AcN and MeOH, the multiconfigurational triplet ππ* states are not much influenced by the solvent.
Geometry | Vacuum | CHCl3 | AcN | MeOH | Water | |||
---|---|---|---|---|---|---|---|---|
ΔEadia | ΔEest | ΔEadia | ΔEest | ΔEadia | ΔEest | ΔEadia | ΔEest | |
a Multiconfigurational triplet states, see Table 1. | ||||||||
1(nOπ*L) | 3.11 | 3.32 | 3.61 | 3.40 | — | 3.62 | — | 3.75 |
1(πHπ*L) | 3.41 | 3.34 | 3.30 | 3.31 | 3.17 | 3.22 | 3.20 | 3.23 |
3(πHπ*L) | 2.82 | 2.73 | 2.71 | 2.68 | — | 2.53 | 2.58 | 2.56 |
3(nOπ*L) | 3.01 | 3.30 | 3.50 | 3.31 | — | 3.56 | — | 3.71 |
3(ππ*)a | 3.55 | 3.55 | — | 3.55 | — | 3.51 | — | 3.52 |
The computed rate constant for radiative emission to the ground state is kf = 6.4 × 107 s−1. However, since the two singlets (S(πHπ*L) and S(nOπ*L)) are isoenergetic at the S(πHπ*L) minimum (see Fig. 4), the rate constant for the radiationless IC process is expected to be larger by some orders of magnitudes as compared to kf. Besides IC, the S(πHπ*L) T(nOπ*L) and S(nOπ*L)
T(πHπ*L) ISC channels are also feasible. The SOMEs for the two channels are i19.40 cm−1 and i42.38 cm−1, respectively. Both the channels are El-Sayed allowed channels, so why do the SOMEs for the two channels differ that much?
In order to understand this we took a closer look at the configuration of the involved electronic states at the minima of the two singlet states. We find that, at the adiabatic minimum of the S(πHπ*L) state, this state comprises 81% of πH → π*L electronic transition. The nO → π*L transition contributes 77% to the T(nOπ*L) configuration. Hence, both the states are highly pure. Now we go to the adiabatic minimum of the S(nOπ*L) and investigate the composition of this and the T(πHπ*L) state. The S(nOπ*L) is nearly pure with 77% contribution from the nO → π*L electronic transition. However, the T(πHπ*L) state shows a mixed configuration character. We find a 70% contribution from the πH → π*L and a 15% contribution from the πH−1 → π*L electronic transition for the T(πHπ*L) state at the S(nOπ*L) minimum. Looking at the involved orbitals (see Fig. 3), we see that the charge transfer character of the S(nOπ*L) T(πHπ*L) transition is enhanced due to the involvement of the πH−1 orbital. In the πH−1 → π*L transition, a marked increase in the electron density at the carbonyl group takes place, explaining the larger SOME for this process.
The scheme in Fig. 5 summarises the various ISC channels available in vacuum. Due to the energetically close vicinity of the (πHπ*L) and (nOπ*L) states, we expect the El-Sayed forbidden channels to be also feasible in addition to the allowed direct ISC path. A rate constant of kISC ≈ 1011 s−1 is computed for the S(πHπ*L) T(nOπ*L) El-Sayed allowed relaxation path. The rate constant of the vibronically allowed S(πHπ*L)
T(πHπ*L) channel is about 3 orders of magnitude smaller and may, hence, safely be ignored in the photophysical discussion. The dark S(nOπ*L) state is populated via IC from the S(πHπ*L) state. The rate constants for the El-Sayed allowed and forbidden ISC processes arising from this dark state are computed to be of the same order of magnitude, namely ≈1010 s−1, with the direct channel being slightly faster than the vibronically induced one.
Mitsui et al. conclude in their paper that in isolated AC an efficient formation of the triplet state occurs after photoexcitation.24 Our results agree and undermine the statement. They further deduce that the dominant nonradiative deactivating channel is S(ππ*) T(nπ*) followed by T(nπ*)
T(ππ*). Here we propose an alternative deactivation channel, namely, S(πHπ*L)
S(nOπ*L) followed by S(nOπ*L)
T(πHπ*L). As mentioned above, the isoenergetic position of the involved singlet states at the minimum of the bright state and the energetically lower-lying adiabatic minimum of the dark state lend credibility to the deactivation channel proposed by us.
Also, surprisingly, the SOME of ≈40 cm−1 reported by Mitsui et al. for the S(ππ*) T(nπ*) channel compares well with that computed by us for the S(nOπ*L)
T(πHπ*L) deactivation process.24
The ISC channel S(πHπ*L) T(nOπ*L) with a computed rate constant of ≈1010 s−1 is the competing channel to the IC process mentioned above. The adiabatic minima of these two states lie only 0.04 eV apart, however, the reverse ISC process T(nOπ*L)
S(πHπ*L) with a rate constant of kReISC = 4.1 × 108 s−1 is about an order of magnitude slower. Summarizing, the radiationless processes dominate the photophysics in the moderately polar solvent, quenching the radiative channels.
The rate constant for the fluorescence decay process is about the same as in vacuum, namely, 6.1 × 107 s−1. A comparison of the LIPs depicted in Fig. 4 and 7 shows clearly the effect of polarity on the PE curves of the various electronic states. Due to the polarity of the solvent, the PES of the (nOπ*L) states are blue shifted as compared to the vacuum and as a result, the 1(πHπ*L) now lies below the 1(nOπ*L) state at its adiabatic minimum. However, the PES of the two singlet states intersect each other at about 0.12 eV above the adiabatic minimum of the bright state. Furthermore, the adiabatic minimum of the 1(nOπ*L) state (at 3.39 eV) is nearly isoenergetic with that of the 1(πHπ*L) (at 3.30 eV), leading us to expect the IC between the two states to play a significant role in depleting the population of the latter state.
Moving along the LIP on the 1(πHπ*L) PES, another crossing is encountered, this time with the 3(nOπ*L) state. The PES of the bright and this triplet state cross each other around 3.36 eV, only 0.06 eV above the adiabatic minimum of the 1(πHπ*L) state. The rate constant for the population transfer via El-Sayed allowed ISC process, 1(πHπ*L) 3(nOπ*L), is computed as 5 × 109 s−1. The adiabatic minima of the involved states being isoenergetic (see Fig. 7), with the minimum of the triplet state lying marginally above that of the singlet state, and considering the fact that under such circumstances the vibrational density of states (VDOS) is not very large, one may not be surprised that this rate constant is about two orders of magnitude smaller than that for the same process in vacuum. Vibronic coupling effects may not be left out when the states with (nOπ*L) and (πHπ*L) character are in such close vicinity. In fact the largest ISC rate constant, kISC ≈ 1011 s−1, has been determined for the El-Sayed forbidden 1(nOπ*L)
3(nOπ*L) ISC process. The direct 1(nOπ*L)
3(πHπ*L) and the vibronically activated 1(πHπ*L)
3(πHπ*L) ISC channels are negligible, being at least an order of magnitude smaller than the other radiationless deactivation channels (see Fig. 8). Thus, we may summarize that the bright state depopulates via IC to 1(nOπ*L) and ISC to 3(nOπ*L). The 3(nOπ*L) is also populated through vibronically promoted ISC out of the 1(nOπ*L) state. This will then be followed by IC from the 3(nOπ*L) to the lowest-lying 3(πHπ*L) state.
It may appear that radiative decay plays no role in the photophysics, kf being a few orders of magnitude smaller than the rate constants for the radiationless transitions. However, the near isoenergeticity of the 1(πHπ*L) and 3(nOπ*L) states leads us to draw parallels with the case of TX in methanol.8,9 The back population of 1(πHπ*L) via reverse ISC (ReISC) from the 3(nOπ*L) is expected. A rate constant of kReISC = 5 × 1010 s−1 is determined computationally for the 3(nOπ*L) 1(πHπ*L) process. This de- and back population of the bright state allows us to predict that delayed fluorescence could be observed for AC in AcN. In the long run, however, the 3(πHπ*L) state would be populated via IC from the 3(nOπ*L) state. The 3(πHπ*L) state would be a long-lived state, as it does not couple strongly to the ground state. Such long-lived triplet states make AC a good triplet sensitizer.
With regard to the role of AC in triplet–triplet annihilation upconversion, we can safely consider AC as a good candidate in vacuum/apolar or slightly polar media. AcN, on the other hand, may already be too polar, though a stable 3(πHπ*L) formation would facilitate the triplet sensitizer function of AC.
As implied by the course of the PES (Fig. 10) and the relative positions of the adiabatic minima of the excited states (Table 2), radiationless deactivation of the 1(πHπ*L) state is no longer feasible. This is corroborated by the absence of a crossing between the 1(πHπ*L) state with other electronic states along the LIP (see Fig. 10). Hence, we conclude, that after photoexcitation, the molecule should relax predominantly via fluorescence. The computed rate constant for this process kf = 8.4 × 107 s−1 compares very well with that reported for AC in ethanol, kf = 7.4 × 107 s−1.28
Mory et al. measured the S1–Sn absorption of AC in ethanol and found a maximum at about 698 nm.29 Our calculations show this maximum to be at 721 nm, in excellent agreement with the experimental results. The estimated adiabatic energy of the 3(πHπ*L) at 2.53 eV agrees well with that reported by Shcherbo et al. for AC in ethanol, namely, 2.61 eV.30 In the literature, we found a work by Miyashita et al. who measured the triplet–triplet absorption of AC in ethanol at room temperature.31 They report two absorption bands at 571 nm and 357 nm, respectively. We calculated the T–Tn absorption spectrum in MeOH (comparable to ethanol) at the T(πHπ*L) minimum (see Fig. 9). The stick spectrum has been convoluted using Gaussian functions with a fwhm (full width at half maximum) of 50 nm. The calculations reproduce the experimental results nearly perfectly (compare Fig. 9 in this work and Fig. 6 in ref. 23). The strongest absorption at 546 nm (571 nm in experiment) in the present work is allocated to the π*L → π*L+2 transition. In the 420–350 nm (experimental value 357 nm) region, we find three transitions contributing to the small hump. Since the experiments were performed at room temperature, perhaps the individual lines could not be identified.
Boużyk et al. presented a theoretical study of the electronic absorption and emission of AC.26 Among other solvents, they have also carried out calculations for AC in AcN using the PM3/CI method. However, they compared these results for the polar, non-protic solvent with that of the results obtained by Miyashita et al. in ethanol which is a polar, protic solvent. We have already emphasised how strongly the proticity of the solvent affects the energy landscape. Hence, the good agreement with experiment achieved by Boużyk et al. is surprising. Since their theoretical results in AcN have been scaled in order to achieve a better fit with the experimental spectra, we have refrained from comparing these findings with our results for AcN.
As implied by the course of the PES (Fig. 10) and the relative positions of the adiabatic minima of the excited states (Table 2), radiationless deactivation of the 1(πHπ*L) state is no longer feasible. This is corroborated by the absence of crossings between the 1(πHπ*L) state with other electronic states along the LIP (see Fig. 10). Hence, we conclude, that after photoexcitation, the molecule should relax predominantly via fluorescence.
Mitsui et al. have discussed the effect of microscopic solvation with water on the photophysics of AC.24 From their analysis they conclude that the ISC in water between the bright singlet state and the triplet (nπ*) states is no longer feasible as the two states experience opposite energy shifts due to the formation of H-bonds.
The LIPs between the minima of the relevant electronic states in water (see Fig. 11) no longer show any crossings that would indicate facilitated radiationless relaxation channels for the photoexcited bright singlet state. Hence, fluorescence should be the dominant decay channel after photoexcitation with a calculated rate constant of kf = 6.2 × 107 s−1. González-Blanco et al. have reported a single exponential decay of AC in water with a time constant of 14.9 ns, which is in excellent agreement with our kf value.32
A detailed quantum chemical analysis of the photophysics of AC in the isolated and solvated states has been presented in this work. The minima and properties of the ground and relevant excited states have been determined using high level electronic structure methods. The effect of solvent polarity and specific bond formation on the photoexcitation processes has been studied in detail. Analog to what we found for related molecules TX and X, the photophysics is controlled by the rather large blue shift experienced by the nOπ*L states. The red shift undergone by the electronic states with 1(πHπ*L) is at least twice as small. The higher-lying multiconfigurational ππ* states are hardly affected by the surrounding medium.
These energy shifts cause the PES of the various electronic states to cross each other. The LIPs between the minima of the states of interest have been calculated in order to get a qualitative picture of the feasible decay processes that may occur after the initial photoexcitation. Besides, computations were carried out to determine the rate constants of the radiative and radiationless processes, wherever possible and plausible.
The rate constant for the decay via fluorescence from the 1(πHπ*L) state is nearly the same in all the cases. In vacuum and CHCl3 the radiative process is quenched by the much faster radiationless IC and ISC channels. Going on to the polar medium of AcN, we still find the IC and ISC processes to be predominant. However, there is an important difference as compared to the previous cases. In AcN we predict that the photoexcited molecule may also show delayed fluorescence. The delayed fluorescence in AcN is possible due to the isoenergeticity of the adiabatic minima of the 1(πHπ*L) and 3(nOπ*L). Such a constellation of energetic layout facilitates the process of reverse ISC from the triplet to the singlet state. Hence, at a larger timescale, fluorescence should be observed in AcN. The hydrogen bond formation occurring in the polar, protic solvents MeOH and water further enhances the energy gap between the πHπ*L and nOπ*L states. This leads to the non availability of the radiationless decay channels, since the states that could be suitable candidates are not reachable due to the large adiabatic energy differences. Hence, in polar, protic solvents the main decay channel is expected to be fluorescence to the ground state after photoexcitation.
AC is an unusual species, in that it may act as a triplet sensitizer (photorelaxation leads to triplet formation) or as a fluorescence marker (photorelaxation does not lead to triplet formation but to fluorescence) depending upon the medium surrounding it. Hence, the analysis of the photophysics as presented in this work may be employed as a guide to tune the surrounding medium in such a way, so as to turn on either the triplet sensitizer or the fluorescent marker properties of AC.
This journal is © The Royal Society of Chemistry 2016 |