Anamika Raya,
Haridas Palb and
Sumanta Bhattacharya*a
aDepartment of Chemistry, The University of Burdwan, Golapbag, Burdwan-713 104, West Bengal, India. E-mail: sum_9974@rediffmail.com
bMolecular Photochemistry Section, Radiation & Photochemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400 085, India
First published on 12th March 2015
This communication reports UV-vis, fluorescence and quantum chemical investigations on the supramolecular interactions of a porphyrazine derivative, namely, 2,7,12,17-tetra-tert-butyl-5,10,15,20-tetraaza-21H,23H-porphine (1) with C60 and C70 in toluene and dichlorobenzene. A synergistic combination of rigidity in the host molecule and a shape-selective binding motif give rise to the strongest binding of C60 with porphyrazine.
Fig. 1 Structure of various tetrapyrrole macrocycles along with 1, namely, (a) porphyrin, (b) phthalocyanine, (c) porphyrazine and (d) 1. |
The Pz, namely, 1 and fullerenes (C60 & C70) are obtained from Aldrich, USA and used without any further purifications. The solvents, namely, toluene and 1,2-dichlorobenzene (DCB) are of UV-vis spectroscopic grade and used after proper distillation. Evidence in favor of interactions between fullerenes and 1 first comes from UV-vis spectroscopic measurements. Addition of 1 to C70 in DCB causes appearance of additional absorption peak in the visible region. This observation is a good fingerprint in favor of complexation between fullerenes and 1, leading us to postulate electronic interactions between these two chromophores, in the ground state. Identification of the charge transfer (CT) spectra in the UV-vis spectroscopic experiment of fullerene–1 complexes is quite interesting as this may be the first example to track the CT behaviour in solution for such complexation process. Fig. 2 nicely demonstrates the colour changes for various CT complexes of 1 with different electron acceptors, i.e., o-chloranil, p-chloranil, tetracyanoethylene (TCNE), 2,3-dichloro-5,6-dicyano-p-benzoquinone (DDQ) and C70 in comparison to pristine 1 and C70. No significant colour change is noticed for C60–1 complex in DCB. Fig. 3 shows the electronic absorption spectra of 1 and mixtures containing 1 + C70, 1 + TCNE, 1 + o-chloranil and 1 + DDQ, in DCB. Fig. S1 (in ESI‡) shows the electronic absorption spectra of remaining mixture containing 1 with p-chloranil in DCB. It should be mentioned at this point that high dielectric constant value of DCB (ε = 9.98) induces CT character as we failed to detect such absorption bands in toluene (ε = 2.39). Moreover, the solubility of fullerene, i.e., C60 in toluene and DCB are estimated to be 2.8 and 27 mg mL−1, respectively.41 Therefore, solvent polarity associated with solvophobic effect are very much responsible behind CT interaction between fullerene and porphyrazine measured in DCB.42 The CT absorption spectra have been analyzed by fitting to the Gaussian function y = y0 + [A/w√(π/2)]exp[−2(x − xc)2/w2], where x and y denote wavelength and absorbance, respectively. The wavelengths at these new absorption maxima (λmax = xc) and the corresponding CT transition energies (hνCT) are summarized in Table 1. CT transition energies of these complexes are related to the vertical electron affinity (EvA) of the electron acceptors by the relation, 2C1 + hνCT = (1/IvD)C1(C1 + hνCT) + {(C2/IvD) + IvD}.43 The vertical electron affinities of C70, p-chloranil, TCNE, DDQ and o-chloranil are collected from literature.44–47 A plot (Fig. S2, in ESI‡) of 2C1 + hνCT versus C1(C1 + hνCT) for a given donor and various electron acceptors yields a slope of 1/IvD from which the value of 1/IvD is determined to be 7.87 eV. Utilizing the values of C1 and C2, degrees of CT (α) may be calculated as follows: α = (C2/2)/[(IvD − EvA + C1)2 + (C2/2)].43,45 The values of α (listed in Table 1) are small and indicate that very little amount of charge transfer occurs in the ground state. Variation of α with the electron affinity of the acceptors is found to be non-linear (Fig. 4). This is quite expected as the model we use here to determine various electronic parameters are based upon interaction of a common donor species with structurally dissimilar electron acceptors. From the CT absorption spectra, we can enumerate the magnitude of the oscillator strength (f). The oscillator strength is estimated using the formula: f = 4.32 × 10−9∫εCTdν, where the term, ∫εCTdν denotes the area under the curve of the extinction coefficient of the absorption band in question vs. wavenumber. The observed oscillator strengths of the CT complexes of 1 are summarized in Table 1. The extinction coefficient is related to the transition dipole moment (μEN) by the relationship: μEN = 0.0952(εmaxΔν1/2/Δν)1/2, where μEN is defined as −e∫ψex(∑ri)ψgdτ. Large value of μEN (Table 1) for C70–1 complex substantiates the CT character in ground state. Another important physicochemical parameter, namely, resonance energy (RN),48 is estimated for the CT complexes of 1 with different electron acceptors in present work and reported in Table 1. RN is determined according to the equation: εmax = 7.7 × 104/{(hνCT/RN) − 3.5}. Moderate value of RN for the C70–1 complex provides very good support in favour of strong electronic interaction between these two species in solution.
System | λCT, nm | hνCT, eV | EvA, eV | α | f | μEN,D | RN, eV |
---|---|---|---|---|---|---|---|
C70–1 | 655 | 1.88 | 2.59 | 0.026 | 0.038 | 7.37 | 0.217 |
o-Chloranil–1 | 643 | 1.92 | 2.87 | 0.020 | 0.020 | 5.33 | 0.205 |
p-Chloranil–1 | 473 | 2.62 | 1.37 | 0.043 | 0.033 | 4.28 | 0.115 |
DDQ–1 | 702 | 1.76 | 3.27 | 0.011 | 0.0028 | 1.85 | 0.022 |
TCNE–1 | 638 | 1.94 | 3.17 | 0.014 | 0.0023 | 1.71 | 0.024 |
It is observed that both fullerenes and 1 attract each other spontaneously in toluene and DCB as evidenced from steady state fluorescence investigations; fluorescence intensity of 1 diminishes gradually by the addition of varying concentration of C60 and C70 upon excitation at its corresponding Soret (S2 ← S0) absorption peak. The steady state fluorescence titration experiments have been performed at constant concentration of 1. The spectral changes finally reach a plateau, indicating that the fluorescence quenching is induced by the complexation. It should be mentioned at this point that a purely diffusion-driven process is ruled out, on the basis of the applied fullerene concentration. The decrease of the fluorescence intensity of 1 suggests a static quenching event inside the well defined fullerene–1 supramolecular complexes. On the basis of the aforementioned results, we reach the conclusions that, in the fullerene–1 complexes, the fluorescence state of 1 is quenched by the addition of electron-accepting C60 and C70. The binding constant (K) of various fullerene–1 complexes have been determined according to Benesi–Hildebrand (BH) fluorescence type equation49 and reported in Table 2. An excellent correlation factor of 0.99 is observed for all the systems studied in the present work. One typical fluorescence quenching curve for C60–1 system in toluene along with fluorescence induced curve and BH plot are shown in Fig. 5. Similar sort of fluorescence spectral representations for C60–1 system in DCB and C70–1 system in both the solvents are demonstrated in Fig. S3–S5, respectively (in ESI‡). Noteworthy is that, ZnPc, a structural analogue of 1, is reported to exhibit ∼3.1 times lower magnitude36 of K than that of 1 during complexation with C60 in toluene. However, both 1 and ZnPc exhibit comparable value of K with C70 as reported in Table 2 (i.e., KC70–1 = 37000 M−1) and literature (i.e., KC70–ZnPc = 25690 M−1).36 It is interesting to note that the increase in magnitude of K led to increase of the fluorescence quenching efficiency. This can be viewed in terms of the fact that the rigidity of 1 introduces tight fixation of C60 compared to C70 unit in fullerene–1 host–guest complex. Table 2 also suggests that C60 always gives larger magnitude of K with 1 in both toluene and DCB in comparison to C70. This is quite unnatural as structural analogue of Pz, i.e., Znporphyrin22,50,51 and ZnPc,36 exhibit stronger binding with C70. As we use the Soret absorption band as our source of excitation wavelength in fluorescence experiment, the second excited singlet state of 1 has been deactivated by the following mechanism:
System | K, M−1 | ΔH0f, kcal mol−1 | τs, ns | ksCS, s−1 | φsCS |
---|---|---|---|---|---|
1 | — | — | 2.86 (2.135) | — | — |
C60–1 | 73000 (44585) | −1.090 | 2.84 (2.05) | 2.40 × 106 (1.9 × 107) | 0.0068 (0.037) |
C70–1 | 37000 (27570) | −0.660 (end-on) | 2.85 (2.11) | 4.2 × 105 (5.6 × 106) | 0.0015 (0.0126) |
−0.100 (side-on) |
Fig. 5 Steady state fluorescence spectral variation of 1 (2.60 × 10−6 M) in presence of C60 (7.15 × 10−6 to 4.60 × 10−5 M) recorded in toluene at 298 K; the inset of Fig. 5 shows fluorescence induced curve and BH plot. λex = 337 nm; λem = 630 nm. |
A similar sort of experimental observations have already been reported by Yin et al.52 and Guldi et al.53 We are very much thankful to one of the learned reviewers for his/her comments, as it creates the opportunity to find out the value of one important photophysical parameter, namely, the Stern–Volmer constants (KSV) for the C60–1 and C70–1 systems in steady state fluorescence titration experiments. If we make a plot of F0/F (where F0 and F are the steady state fluorescence intensity of 1 in absence and presence of fullerenes, respectively) vs. concentration of the quencher (here fullerene) at dilute concentration range of the quencher, excellent linear plots are obtained for all the systems studied in present work. However, at higher concentrations of quencher, the plots generate an upward curvature. From the slope of such plots, KSV may be enumerated. The values of KSV for C60–1 and C70–1 systems in toluene (and DCB) are estimated to be 69620 M−1 (56600 M−1) and 65410 M−1 (51200 M−1), respectively. Typical fluorescence spectral variation for determination of SV plot for C60–1 system in DCB is shown Fig. 6. SV plots for C60–1 system in toluene and C70–1 system in both the solvents are provided as Fig. S6–S8, respectively (in ESI‡). The bimolecular quenching constant, kq, calculated from the linear segment of the SV plots revealed that the kq values are 3 order of magnitude higher than that expected for diffusion controlled processes in DCB and toluene, respectively. The linear nature of SV plots at dilute concentration range of quencher along with order of magnitude of kq suggest that static type quenching is the main process for quenching of fluorescence of 1 in presence of fullerenes.54,55 The KSV values measured in this manuscript fall into the range between pure diffusion controlled quenching and hydrogen bonding facilitated quenching in solutions where fullerene is a fluorescent quencher as suggested by one of the learned reviewer.56–58 The above results, therefore, provide enough evidences in favour of the complexation between fullerene and 1 in both toluene and DCB.
Fig. 6 Steady state fluorescence spectral variation of 1 (2.6 × 10−6 M) in presence of C60 (1.0 × 10−6 to 1.4 × 10−5 M) recorded in DCB at 298 K; the inset of Fig. 6 shows SV plot. λex = 340 nm; λem = 632 nm. |
Another interesting aspect of the present investigations is that 1 did not exhibit significant selectivity in binding toward C60 or C70. The selectivities of C60 over C70 are found to be ∼2.0 and 1.6 in toluene and DCB, respectively. Although, the present selectivity ratio of porphyrazine is comparable to that of azacalix[m]arene[n]pyridine (∼1.94),59 it is much lower than those of calixarene diporphyrin (∼4.3),34 cyclic dimers of Zn–porphyrins (∼25.5)60 and corresponding H2–diporphyrin (∼32)61 complexes.
In our present work, a detailed conformational analysis of the individual components as well as the fullerene–1 complexes have been performed in vacuo by HF/3-21G calculations. The geometric parameters of the complexes are obtained after complete energy minimization. To limit the number of basis functions, we have replaced tert-butyl group with methyl in the structure of 1. The strongest binding of C60 towards 1 is reflected in terms of larger negative value of heat of formation (ΔHf°) for this complex compared to C70–1 estimated in vacuo (Table 2). Comparing the relative stability of C60–1 and C70–1 complexes in terms of ΔHf° value (listed in Table 2), it is revealed that C70–1 complex is supposed to be less stable in side-on orientation of C70 (Fig. S9, in ESI‡) compared to its end-on orientation (Fig. 7(b)) by 0.560 kcal mol−1. We anticipate that the contact area of 1 with C60 (Fig. 7(a)) would not differ much with respect to that of C70. Estimation of the surface area of the C70 complexes of 1 (side-on) and 1 (end-on) are estimated to be 857.32 and 857.23, Å2, respectively, which is very much close to that of C60–1 complex, i.e., 836.12 Å2. Thus, unlike fullerene–porphyrin interaction,22 computations of the preferred conformations of fullerene–1 complexes confirm that van der Waals attractive forces play negligible role in present case. The frontier highest occupied molecular orbitals (HOMO) and lowest unoccupied molecular orbitals (LUMO) for the investigated supramolecules have been obtained by DFT/B3LYP/6-31G* method and are demonstrated in Fig. 8 and S10 (in ESI‡) for C60–1 and C70–1 systems, respectively. In the studied systems, the electron distribution of HOMO and LUMO is exclusively found to be located on the 1 entity and fullerene spheroid, respectively. These results suggest weak charge-transfer interaction between 1 and fullerene in the ground state and at the same time, prove the existence of charge-transfer transition from 1 entity to fullerene spheroid. It is also observed that, HOMO of both C60–1 (−5.5945 eV) and C70–1 (−5.5964 eV) complexes corroborate fairly well with that of uncomplexed 1 (−5.6013 eV) while LUMO of these complexes, i.e., −3.0820 eV and −3.0912 eV lie very close to that of C60 (−3.1469 eV) and C70 (−3.1427 eV), respectively. The small HOMO–LUMO gap gives formidable support in favour of weak CT from 1 to fullerene as evidenced from low values of α tabulated in Table 1. Molecular electrostatic potential (MEP) maps have also been generated for the investigated supramolecules in present work with the help of DFT calculations and visualized in Fig. 9. Fig. 9 gives strong indication of electronic redistribution in donor (here 1) and acceptor (here fullerene) during non-covalent interaction. It may be mentioned here that validity of molecular orbitals generated by density functional methods is already being recognized.62 The accuracy of these methods, especially B3LYP, is being well demonstrated by D'Souza54 and Schaefer et al.63 and, more recently, on covalently linked C60–Znporphyrin system64 and fullerene-based photoactive layers for the construction of heteroconjugation solar cells.65
Fig. 7 Ab initio (HF/3-21G) optimized geometric structures of (a) C60–1 complexes and (b) C70–1 (end-on orientation of C70) done in vacuo. The calculations are done using SPARTAN '06 software. |
Fig. 8 DFT (B3LYP/6-31G*) calculated frontier (a) HOMO and (b) LUMO for the C60–1 complex done in vacuo. The calculations are done using SPARTAN '14 software. |
Fig. 9 MEP maps of (a) C60–1 and (b) C70–1 (end-on orientation of C70) systems done by DFT (B3LYP/6-31G*) calculations in vacuo. The calculations are done using SPARTAN '14 software. |
Evidence in favor of end-on orientation of C70 toward the plane of the 1 gets strong experimental support from variable temperature 13C NMR experiment of C60–1 (Fig. 10) and C70–1 systems (Fig. 11) recorded in toluene-d8 medium. The 13C NMR spectra of C60 and C70 shown in Fig. 10 and 11, respectively, are not decoupled. Similar to the case of C60, the 13C NMR signals due to C70 in C70–1 complex are upfield-shifted from those of uncomplexed C70. Since such an upfield shift is more pronounced for the pole carbon signals than for those located at the equator, C70 appears to adopt an end-on conformation with respect to 1 moiety in order to maximize the host–guest interaction. Both theoretical and experimental observations reveal that such a binding motif of C70 toward 1 takes place to maximize the electrostatic interactions. As shown in Fig. 10, 13C-enriched C60 and 1 in Tol-d8 shows a single resonance for C60 at 298 K due to time-averaging of bound and unbound C60. Upon cooling, this peak gradually broadens and at 243 K, it appears as separate singlet signal corresponding to bound C60 (140.3 ppm) and unbound C60 (142.8 ppm). The upfield shifts of the fullerene in the bound states are the result of ring currents effects of the porphyrazine π-system on the fullerene. Similar sort of observations is noticed by Boyd et al. for calix[4]arene-linked bisporphyrin host for fullerene.66
Apart from steady state fluorescence measurements, we have performed detailed nanosecond time-resolved fluorescence experiment for the fullerene–1 systems in toluene and DCB using 374 nm laser pulse. The experiment has been carried out at a fixed concentration of 1. The time-resolved fluorescence measurement shows mono-exponential decay for uncomplexed 1; lifetime value of the singlet excited state (τs) of 1 is measured to be 2.86 (Fig. S11 in ESI‡) and 2.135 ns (Fig. 12(a)) in toluene and DCB, respectively. When the experiment is carried out in toluene, no significant change in the value of τs is observed for C60 (2.84 ns, Fig. S11 in ESI‡) and C70 (2.85 ns, Fig. S11 in ESI‡). However, slight reduction in τs value is observed for 1 in presence of C70 (2.11 ns, Fig. 12(b)) and C60 (2.05 ns, Fig. 12(c)) measured in DCB. Utilizing the value of τs in DCB, rate constant of charge-separation of the excited singlet (ksCS) and quantum yield of charge-separation for the excited singlet (φsCS)67–69 have been evaluated and listed in Table 2. The most fascinating thing is that in DCB, C60–1 complex exhibits little charge-separation as φsCS value is estimated to be 0.037. The lifetime experiment, therefore, suggests enhancement in electrostatic interaction between 1 and fullerene in solvent having high dielectric constant value. To validate the presence of the static quenching mechanism in our present investigations, we have performed detailed nanosecond time-resolved fluorescence measurements for fullerene–1 complexes in both toluene and DCB. The titration experiment is carried out at a fixed concentration of 1 and a variable concentration of C60 and C70. It is observed that, upon the gradual addition of fullerenes C60 and C70 in toluene (and DCB), the value of τs suffers little change, and mono-exponential decay is followed. Plot of τs/τ vs. concentration of quencher resulting from time-resolved titration experiment for C60–1 and C70–1 systems recorded in DCB are shown in Fig. S12(a) & S12(b) respectively (in ESI‡); here, τs and τ denote the lifetime value of 1 in absence and presence of quencher, respectively. Since the value of τs does not change much in the presence of the quencher, there is no question of observing the lifetime value in another time range like the picosecond or femtosecond region. This important photophysical observation establishes that diffusion controlled mechanism is not operative behind quenching of the fluorescence intensity of 1 in presence of C60 and C70 in both toluene and DCB. Similar sort of observations is noticed when the experiments are performed in toluene. Scanning the emission feature in longer wavelength regions (700–800 nm), it is revealed a weak emission band at 763 nm corresponding to the singlet emission of C60 in DCB (marked as red colour in Fig. 13). The fluorescence intensity of this band for a fixed concentration of C60 is found to remain almost same as that obtained in presence of 1. This result indicates that in DCB, there is a relaxation pathway from the excited singlet state of 1 to that of the C60. Similar sort of observations have been made for C60–1 system in toluene (Fig. S13‡) studied in present work.
Fig. 13 Steady state fluorescence spectral variation of C60 (2.0 × 10−5 M, red colour line) in presence of 1 (2.5 × 10−6 M, blue colour line) recorded in DCB at 298 K; λex = 340 nm; λem = 763 nm. |
In conclusion, 1 has been demonstrated to be good shape-selective host molecule for fullerene guests. The fullerene–1 supramolecular recognition element may find potential application to construct discrete host–guest complexes in near future. Further studies on photophysical characterization and single crystal X-ray analysis of the complexes of porphyrazine with various fullerene derivatives are in progress.
Footnotes |
† The authors wish to dedicate this work to Late Mr S. B. Banerjee for his constant encouragement in pursuing research work. |
‡ Electronic supplementary information (ESI) available: CT analysis of p-chloranil–1 complex in DCB (Fig. S1), plot for determination of IvD of 1 (Fig. S2), steady state fluorescence spectra, induced curves and BH plots for the determination of the binding constants of 1 with C60 and C70 in different solvents (Fig. S3–S5), steady state fluorescence spectral variation for the determinations of SV constants for C60–1 complex in toluene (Fig. S6) and C70–1 complex in toluene (Fig. S7) and DCB (Fig. S8), optimized structure of C70–1 complex at side-on orientation of C70 (Fig. S9), B3LYP/6-31G* optimized electronic, i.e., HOMO and LUMO, structure of C70–1 system in end-on orientation of C70 (Fig. S10), decay profile of 1 in absence and presence of C60 and C70 recorded in toluene (Fig. S11), plot of τs/τ vs. concentration of quencher for C60–1 and C70–1 systems measured in DCB (Fig. S12), and finally, fluorescence spectral variation of C60 in presence of 1 recorded in toluene (Fig. S13). See DOI: 10.1039/c5ra02003d |
This journal is © The Royal Society of Chemistry 2015 |