Yang Liua,
Wangchang Lia,
Peiying Lia,
Yanmin Guoa,
Peng Cui*b and
Zhuxia Zhang*c
aCollege of Materials Science and Engineering, Key Laboratory of Interface Science and Engineering in Advanced Materials, Taiyuan University of Technology, Taiyuan 030024, China
bThe Novel Computer Architecture Laboratory, School of Information, Guizhou University of Finance and Economics, Guiyang 550025, China. E-mail: pcui@mail.gufe.edu.cn
cCollege of Chemistry, Taiyuan University of Technology, Taiyuan 030024, China. E-mail: zhangzhuxia@tyut.edu.cn
First published on 3rd February 2023
The encapsulation of fullerenes by carbon nanorings has gained increasing attention because of the unique molecular structure and special properties of the formed complexes. The host–guest interactions between the fullerenes and the carbon nanorings can influence the metal ion orientation and the molecular electronic structure. In this study, we hooped a series of carbide cluster metallofullerenes, namely Sc2C2@C2v(5)-C80, Sc2C2@C3v(8)-C82, and Sc2C2@D2d(23)-C84, with molecular carbon nanorings of [12]cycloparaphenylene ([12]CPP) and perfluoro[12]cycloparaphenylene (PF[12]CPP). The formed complexes were computationally studied via dispersion-corrected density functional theory calculations. The results showed that the deformation rate of PF[12]CPP after the formation of the fullerene-containing complexes was significantly smaller than that of [12]CPP. The binding energy and thermodynamic information showed that PF[12]CPP was more suitable for fullerene encapsulation. Moreover, charge population analysis showed that PF[12]CPP transferred more electrons to Sc2C2@C2n (n = 40, 41, and 42) compared with [12]CPP. Energy decomposition and real-space function analyses of host–guest interactions revealed the characteristics and nature of the noncovalent interactions in the supramolecules. These results provide theoretical support for the study of host–guest systems based on metallofullerenes.
One of the interesting approaches to this problem is the use of receptors based on the so-called concave–convex π–π interactions. Through the use of concave receptors capable of adapting to the convex surface of fullerenes, the number of atom contacts can be maximized at short distances, enabling the maximization of stabilization through dispersion.9 Considering the spherical form of metallofullerene, the circular carbon nanorings are a suitable host for accommodating metallofullerenes and altering their properties. Molecular carbon nanorings are distinct π systems resembling carbon nanotube (CNT) segments, and both can form stable supramolecular systems through concave–convex π–π interactions.10 The most representative examples of carbon nanorings are cyclophenylpropylenes (CPPs), which consist of conjugated multibenzene rings that represent the shortest armchair carbon nanotube segment. Since the discovery of C60 was trapped into [10]CPP by Yamago et al.,11 the following class of supramolecular systems has been developed: C59N⊂[10]CPP,12 C70⊂[10–11]CPP,13 La@C82⊂[11]CPP,14 Gd@C82⊂[11]CPP,15 Y2@C79N⊂[12]CPP,16 and Sc2C2@C82⊂[12]CPP17 etc.
Fluorine incorporation into organic compounds considerably modifies the properties of the compounds, including their polarity, solubility, and lipophilicity. Therefore, perfluorocycloparaphenylenes (PFCPPs), another class of carbon nanorings with a comparable structure to CPPs, have attracted extensive attention. PFCPPs are a class of highly strained ring-shaped perfluoroarenes in which all hydrogen atoms of the corresponding CPPs are replaced with fluorine atoms. The structural and electronic properties of PFCPPs are expected to change significantly following the replacement of hydrogen by fluorine.18 Thus, elucidating the implications of these changes in PFCPPs for supramolecular systems such as fullerene⊂CPP is vital. Moreover, clarifying the fullerene molecule accommodation property of PFCPPs as compared with CPPs is important. Recently, Itami et al. synthesized and isolated PFCPPs.19 However, studies on the combination of PFCPPs and fullerenes are few. In addition, experimentally elucidating the nature of these noncovalent interactions is difficult. Therefore, the computational groups have a good opportunity to explore in depth the fundamental properties of the π–π noncovalent interactions of these host–guest systems.
In this study, we hooped a series of carbide cluster metallofullerenes, namely Sc2C2@C2v(5)-C80, Sc2C2@C3v(8)-C82, and Sc2C2@D2d(23)-C84, with nanorings of [12]CPP and PF[12]CPP. Although Zhang and Lu et al.17,20 have used [12]CPP to encapsulate Sc2C2@C3v(8)-C82, its two isomers Sc2C2@C2v(5)-C80, and Sc2C2@D2d(23)-C84 have not been investigated. Thus, in the present study, a series of complexes formed by host molecules ([12]CPP and PF[12]CPP) and guest molecules Sc2C2@C2n (n = 40, 41, and 42) were computationally studied using the density functional theory (DFT). The molecular-level π–π interactions of two systems were theoretically assessed to explore the discrepancies between them. The results of this study can elucidate the structures and properties of these supramolecular complexes and contribute to the design of novel supramolecular nanoarchitectures for material chemistry applications.
In order to obtain accurate binding energy for supramolecular complexes, the ωB97M-V functional in combination with a large def2-TZVPP basis set were employed in single-point calculations.25,26 The Counterpoise correction was used to solve the problem of basis set superposition error (BSSE).27 The ωB97M-V is fairly accurate for evaluating intermolecular interactions, and the single-point calculations were conducted using the ORCA 5.0.3 program28 via the RIJCOSX technique29 to accelerate the calculations.
To determine the nature of weak interactions in the supramolecular complexes, the physical components of intermolecular interactions were quantitatively analyzed via symmetry-adapted perturbation theory (SAPT) analysis using the PSI4 code.30 Considering the size of the studied systems, we had to use the simplest SAPT level, SAPT0,31 in combination with a small def2-SVP basis set. The adopted SAPT0 included the following terms:
ESAPT0 = Eelest + Eexch + Eind + Edisp, | (1) |
(2) |
Electrostatic potential (ESP) and independent gradient model based on Hirshfeld partition (IGMH) analyses were all finished by Multiwfn 3.8(dev) code.32–37 Thermochemistry properties were calculated using the Shermo package.38 All isosurface and molecular structure maps were constructed using the VMD software.39
The [12]CPP featured negative regions associated with phenyl rings and positive regions on the rim of the ring, owing to hydrogen atoms (Fig. 1a). Therefore, the cavity was mostly negative and will favor the inclusion of electron-deficient species. The ESP of the van der Waals (vdW) surface of [12]CPP was distributed in the range of −16.79 to 16.22 kcal mol−1. The overall distribution characteristics of ESP can be found from the histogram of the area distribution of different ESP intervals. The relatively uniform distribution of the surface area in different ESP ranges for [12]CPP is shown in Fig. 1(a-right), and the positive and negative areas were calculated as 451.1 and 494.1 Å2, respectively, which correspond to 47.8% and 52.2% of the entire vdW surface, respectively. As the polarity of F atoms was much greater than that of H atoms, PF[12]CPP and [12]CPP considerably differed in ESP distribution. PF[12]CPP displayed positive regions in the internal region of the ring and negative regions on the ring rim (Fig. 1b). The ESP of the vdW surface of PF[12]CPP was distributed in the range −9.04 to 18.29 kcal mol−1, and the local minima and maxima of the ESP values were greater (more positive) than those of [12]CPP (−16.79 to 16.22 kcal mol−1). Fig. 1(b-right) displays the inhomogeneous surface area for different ESP ranges of PF[12]CPP compared with [12]CPP. The areas of the positive and negative parts of PF[12]CPP were respectively calculated as 587.7 and 518.1 Å2, corresponding to 46.9% and 53.1% of the entire vdW surface, respectively. The ESP of the vdW surface of PF[12]CPP was not distributed in the interval −18 to 12 kcal mol−1, whereas 18.3% of that of [12]CPP was distributed in this interval.
In previous studies, C60 and C70 showed almost zero ESP, with marginally positive regions, and C76 and C78 exhibited the same behavior.9 However, the occurrence of the endohedral unit in endofullerenes introduces larger asymmetries in the charge distribution, resulting in noticeable positive and negative ESP regions. Sc2C2@C80, Sc2C2@C82, and Sc2C2@C84 (Fig. 1c–e) showed a markedly positive region, where the Sc atom contacted the carbon cage, while the carbon cages near the C–C bonds of the Sc2C2 endohedral cluster exhibited a negative region. Moreover, the positive region of Sc2C2@C80 was more concentrated than those of Sc2C2@C82 and Sc2C2@C84, attributable to the small size of C80, which resulted in a dihedral angle of only 112.7° for Sc–C2–Sc in the Sc2C2 cluster, which is significantly smaller than those of C82 and C84 (126.0° and 180.0°, respectively; see Table S1† for specific geometric parameters). The ESPs of the vdW surface of Sc2C2@C80, Sc2C2@C82, and Sc2C2@C84 were distributed in the range −4.81 to 13.74 kcal mol−1, −6.07 to 13.87 kcal mol−1, and −2.35 to 11.96 kcal mol−1, respectively. The ESP of Sc2C2@C84 was not distributed in the interval of −7.0 to −3.5 kcal mol−1, whereas 2.9% and 7.9% of the ESPs of Sc2C2@C82 and Sc2C2@C84 were distributed in this interval. In addition, Sc2C2@C84 showed almost the largest marginally positive ESP regions (50.9% in the 0 to 3.5 kcal mol−1 interval); the corresponding percentages for Sc2C2@C80 and Sc2C2@C82 were 42.1% and 32.5%, respectively, consistent with the ESP diagram (Fig. 1e).
Furthermore, we explored the different orientations of guest molecules in the host molecule through the ESP study of host–guest molecules and according to the electrostatically complementary principle. We simulated two typical rotation processes by calculating the energy curves of Sc2C2@C80, Sc2C2@C82, and Sc2C2@C84 along their rotation path inside [12]CPP (Fig. 2). All isomers were optimized at the M06-2X-D3/6-31G(d) level (a total of 69 complexes). Path I suggests that two Sc atoms move perpendicularly along [12]CPP plane, while Path II suggests that two Sc atoms remained adhered to the wall of the nanoring during the rotation process, indicating that the positive region where Sc atoms contacted the carbon cage could overlap maximally with the negative region inside the [12]CPP cavity. Path II was the best path to rotate because its energy curve was below that of Path I (Fig. 2c and d). When Sc2C2@C82 rotated in [12]CPP, the barrier difference between Paths I and II was 4.75 kcal mol−1, while for Sc2C2@C84, the barrier difference was 4.27 kcal mol−1. In contrast to both complexes, for Sc2C2@C80 (Fig. 2b), the energy curve of Path II was not significantly lower than that of Path I because the positive region of Sc2C2@C80 was more concentrated (Fig. 1) than those of Sc2C2@C82 and Sc2C2@C84, as mentioned earlier. Path II was still the best rotation path because its barrier was 2.13 kcal mol−1 lower than that of Path I. When both Sc atoms were on parallel sites of the [12]CPP plane (θ = 180°), the conformer possessed minimum energy for Path I, because the positive region could overlap maximally with the negative region. The results of our study agree with the findings by Zhao and Lu et al. and explain why metallofullerenes such as Sc3N@C80 and Y2@C79N have the best rotation paths when rotated horizontally within [12]CPP.16,20 Likewise, for Sc2C2@C2n⊂PF[12]CPP (n = 40, 41, and 42), we positioned Sc atoms in parallel and perpendicular to F atoms, where the positive and negative regions overlapped to the maximum and minimum extents. The energy of the isomers was reduced when Sc atoms were parallel to F atoms. Section 3.2 focuses on the structures of the supramolecular complexes. The above results show that the molecular ESP is useful for identifying the stable structures of supramolecular complexes such as metallofullerene⊂CPP.
Species | Dc–ca (Å) | Long axis (Å) | Short axis (Å) | Aspect ratio | Rdb (%) | θc (deg) | Dipole (Debye) |
---|---|---|---|---|---|---|---|
a Dc–c is defined as the distance between the two centers of host and guest.b Rd (deformation ratio of the host ring) = (aspect ratio of the host in the complex – aspect ratio of the free host)/aspect ratio of the free host.c θ is defined as the dihedral angle range between ring units of the host [12]CPP and PF[12]CPP. | |||||||
Free [12]CPP | — | 16.407 | 16.283 | 1.007 | — | 28.1–39.1 | 0 |
Sc2C2@C80⊂[12]CPP | 1.519 | 18.000 | 14.384 | 1.251 | 24.2 | 27.5–47.3 | 1.647 |
Sc2C2@C82⊂[12]CPP | 1.110 | 17.537 | 15.009 | 1.168 | 16.0 | 18.2–40.5 | 1.970 |
Sc2C2@C84⊂[12]CPP | 1.102 | 17.278 | 15.160 | 1.140 | 13.2 | 25.2–36.2 | 1.160 |
Free PF[12]CPP | — | 15.955 | 14.960 | 1.069 | — | 46.1–60.1 | 0.014 |
Sc2C2@C80⊂PF[12]CPP | 0.608 | 16.866 | 14.929 | 1.130 | 5.7 | 44.1–55.7 | 1.297 |
Sc2C2@C82⊂PF[12]CPP | 0.530 | 16.776 | 14.905 | 1.126 | 5.3 | 42.7–56.2 | 0.898 |
Sc2C2@C84⊂PF[12]CPP | 0.674 | 16.400 | 15.667 | 1.047 | −2.1 | 40.5–54.0 | 0.536 |
As shown in Table 1, the free [12]CPP and PF[12]CPP were almost circular, with calculated aspect ratios of 1.007 and 1.069, respectively. Upon the encapsulation of the fullerenes, the hosts were deformed into an elliptical shape. The deformation ratio of the host ring (Rd), defined in the footnote of Table 1, is a parameter for describing the host ring deformation after the formation of the complexes. Rd decreased with increasing carbon cage size (Table 1).
The Rd values of Sc2C2@C80⊂[12]CPP, Sc2C2@C82⊂[12]CPP, and Sc2C2@C84⊂[12]CPP were 24.2%, 16.0%, and 13.2%, respectively. The Rd of Sc2C2@C2n⊂PF[12]CPP was significantly lower than that of Sc2C2@C2n⊂[12]CPP. The Rd of Sc2C2@C84⊂PF[12]CPP (−2.1%) was lower than that of Sc2C2@C80⊂PF[12]CPP (5.7%). The −2.1% Rd of Sc2C2@C84⊂PF[12]CPP indicates that the interaction between the PF[12]CPP host and Sc2C2@C84 guest makes PF[12]CPP closer to a circular geometry compared with the free PF[12]CPP. In addition, Table 2 calculates the deformation energy of the monomers. The deformation energy (ΔEdef) of PF[12]CPP is lower than that of [12]CPP. The dihedral angle (θ) between ring units of the host [12]CPP and PF[12]CPP is another necessary parameter for describing the structure and deformation of the host ring. The θ values of the free host PF[12]CPP were widely distributed (46.1–60.1°), whereas those of the free hosts [12]CPP were narrowly distributed (28.1–39.1°), because the steric repulsion between the F atoms was greater than that between the H atoms. Additionally, the dihedral angles in the complexes change to different extents, probably suggesting a larger deformation energy of the host and a very strong host–guest interaction based on repulsion deriving from steric energy. Furthermore, because of their quasi-symmetry, both the free [12]CPP and PF[12]CPP possessed small dipole moments nearly equal to zero. However, all of the host–guest complexes were polar supramolecules with a certain number of dipole moments, attributable to the electron-rich property of the host and the electron-deficient nature of the guest fullerene.
Complexes | ΔEa | BSSEa | ΔEcpa | ΔEdef | ΔG | ΔH | ΔS | |
---|---|---|---|---|---|---|---|---|
Full | Ring | |||||||
a Obtained at the ωB97M-V/TZVPP level//M06-2X-D3/6-31G(d) level; BSSE refers to basis set superposition error. | ||||||||
Sc2C2@C80⊂[12]CPP | −45.66 | 1.47 | −44.19 | 0.55 | 5.43 | −102.09 | −178.97 | −257.88 |
Sc2C2@C82⊂[12]CPP | −48.40 | 1.48 | −46.92 | 0.49 | 3.94 | −112.83 | −189.95 | −258.69 |
Sc2C2@C84⊂[12]CPP | −49.45 | 1.37 | −48.08 | 0.58 | 1.46 | −123.22 | −195.02 | −240.81 |
Sc2C2@C80⊂PF[12]CPP | −57.23 | 5.86 | −51.37 | 0.80 | 1.23 | −135.05 | −212.22 | −258.84 |
Sc2C2@C82⊂PF[12]CPP | −58.00 | 5.90 | −52.10 | 0.19 | 1.15 | −136.34 | −215.65 | −266.00 |
Sc2C2@C84⊂PF[12]CPP | −60.82 | 5.76 | −55.06 | 0.26 | 0.93 | −153.38 | −227.77 | −249.52 |
To accurately describe the position of the endofullerene guests in the [12]CPP and PF[12]CPP hosts, the distances (Dc–c) between the centers of host and guest were determined (Table 1). The Dc–c value is not zero, which indicates that the centers of the host and the guest did not completely overlap in all of the configurations of the complexes. A close relationship existed between the centers of the host and the guest (Dc–c) and the π–π interaction regions, and the shorter Dc–c was more advantageous to maximize the attractive interactions between the host and guests. The guest Sc2C2@C80 of Sc2C2@C80⊂[12]CPP, Sc2C2@C82 of Sc2C2@C82⊂[12]CPP, and Sc2C2@C84 of Sc2C2@C84⊂[12]CPP shifted to one side of the host ring by 1.519, 1.110, and 1.105 Å away from the [12]CPP center, respectively. The Dc–c values of Sc2C2@C80⊂PF[12]CPP, Sc2C2@C82⊂PF[12]CPP, and Sc2C2@C80⊂PF[12]CPP were 0.608, 0.530, and 0.674 Å, respectively. The Dc–c value of Sc2C2@C2n⊂PF[12]CPP was markedly smaller than that of Sc2C2@C2n⊂[12]CPP, suggesting that Sc2C2@C2n⊂PF[12]CPP exhibited stronger host–guest interaction.
All of the complexes of the Sc2C2@C2n⊂[12]CPP system exhibited similar general aspects (Fig. 3), with the fullerenes fitting within the [12]CPP cavity and establishing favorable interactions with all of the phenylene units. The fullerenes were oriented inside the cavity in such a way that the endohedral Sc atoms were in the plane of the nanoring which means that carbon atoms interacting with the endohedral unit were always oriented towards phenylene rings. All of the complexes of the Sc2C2@C2n⊂PF[12]CPP system also exhibited similar general features, with the fullerene suitably fitting within the PF[12]CPP cavity and establishing helpful interactions with F atoms of ring units. The Sc atoms were located close to F atoms in PF[12]CPP, which indicates that carbon atoms interacting with the endohedral unit were oriented toward the F atoms of ring units. Structures with other orientations relative to hosts were less stable. These behaviors were probably related to the asymmetry of the electron charge shown in Fig. 1. The ESP of endofullerenes was more positive in the regions where the endohedral unit contacted the carbon cage. Considering that the ESP of the [12]CPP cavity and the ring rim of PF[12]CPP were markedly negative, metallofullerenes were oriented such that their positive regions were directed toward those negative regions of hosts to facilitate the stabilization of electrostatic interactions. Fig. 3 also shows the interfacial distances (di) between hosts and guests. Here, di denotes the interfacial distance between the centroid of a ring unit of the host and the nearest centroid of a hexagon or pentagon of an endofullerene. The di values of Sc2C2@C80⊂[12]CPP, Sc2C2@C82⊂[12]CPP, and Sc2C2@C84⊂[12]CPP were 3.547–6.964 Å (Fig. 3a), 3.536–6.113 Å (Fig. 3b), and 3.622–5.812 Å (Fig. 3c), respectively, showing the wide distribution of di. This suggests that the attractive force of endofullerenes along the [12]CPP was nonuniform. A narrower interface will have a stronger force between Sc2C2@C2n and [12]CPP than a wider interface, because the interfacial distances were not significantly greater than 3.4 Å, the vdW distance between graphite sheets. The di values of Sc2C2@C80⊂PF[12]CPP, Sc2C2@C82⊂PF[12]CPP, and Sc2C2@C84⊂PF[12]CPP were 3.796–4.924 Å (Fig. 3d), 3.776–5.166 Å (Fig. 3e), and 3.619–5.369 Å (Fig. 3f), which shows a relatively narrow di distribution compared with the Sc2C2@C2n⊂[12]CPP system, respectively. Thus, the attractive force of endofullerenes along PF[12]CPP was relatively uniform compared with that along [12]CPP.
The effect of the host molecule on the endohedral unit is also of concern. The geometric parameters of Sc2C2 are presented in Table S1.† The cluster of Sc2C2@C80 was bent like a butterfly, with two tightly bonded C atoms in the cage center. The two Sc atoms existed between the C2 unit and the cage. The Sc–Sc distance in the Sc2C2@C80 cage was 3.692 Å, and the Sc–C2–Sc dihedral angle was 112.7°. Compared with the Sc2C2@C80 cage, the Sc2C2@C80⊂[12]CPP system exhibited a 0.006 Å higher Sc–Sc distance and a 0.4° lower Sc–C2–Sc dihedral angle, while the Sc2C2@C80⊂PF[12]CPP system exhibited a 0.015 Å lower Sc–Sc distance and 0.9° lower Sc–C2–Sc dihedral angle. Moreover, the Sc–Sc distance and dihedral angle of the Sc2C2 cluster of Sc2C2@C82 were 3.939 Å and 126.0°, greater than those of Sc2C2 cluster of Sc2C2@C80, which indicates that the values of both parameters increased with increasing carbon cage size. In addition, compared with Sc2C2@C82, the Sc–Sc distances of the Sc2C2@C82⊂[12]CPP and Sc2C2@C82⊂PF[12]CPP systems were reduced by 0.002 Å and 0.003 Å, respectively, while the Sc–C2–Sc dihedral angles were reduced by 0.4° and 0.2°, respectively. The Sc2C2 cluster of Sc2C2@C84 exhibited the largest Sc–Sc distance and dihedral angle: 4.473 Å and 180.0°, respectively. Compared with Sc2C2@C84, the Sc–Sc distances of the Sc2C2@C84⊂[12]CPP and Sc2C2@C84⊂PF[12]CPP systems were 0.004 Å and 0.003 Å higher, respectively, while the Sc–C2–Sc dihedral angles did not change. The C–C distance was almost constant for all complexes (around 1.265 Å). Overall, the weak interaction between the host and the guest did not significantly affect the endohedral clusters.
Sc2C2@C84 exhibited the greatest binding energy to the host molecule, attributable to the larger carbon cage size of C84 compared with C80 and C82, which allows for a large region of π–π contact between the host molecules ([12]CPP and PF[12]CPP) and C84. The binding energy of the Sc2C2@C2n⊂PF[12]CPP system was about 7 kcal mol−1 greater than that of the Sc2C2@C2n⊂[12]CPP system, which indicates that PF[12]CPP is more suitable for accommodating fullerenes than [12]CPP. In addition, to investigate the effect of Sc2C2 clusters on the noncovalent interactions between the host and guest, the binding of C80, C82, and C84 without Sc2C2 clusters to [12]CPP and PF[12]CPP was calculated (see Table S2†). The results show that Sc2C2@C2n binds more stable to [12]CPP and PF[12]CPP compared with C2n. This suggests that electrostatic interactions due to the introduction of endohedral clusters can promote the stabilization of supramolecular complexes. The thermodynamic information of all of the complexes at 298.15 K and 1 atm were obtained via DFT calculations at the M06-2X-D3/6-31G(d) level of theory (Table 2). The relative trends of the ΔG and ΔH of the Sc2C2@C2n⊂[12]CPP and Sc2C2@C2n⊂PF[12]CPP complexes were well consistent with those of ΔECP. The processes of Sc2C2@C80, Sc2C2@C82, and Sc2C2@C84 binding to [12]CPP and PF[12]CPP at 298.15 K were spontaneous, with ΔG values of −102.09 to −123.22 kJ mol−1 and −135.05 to −153.38 kJ mol−1 for Sc2C2@C2n⊂[12]CPP and Sc2C2@C2n⊂PF[12]CPP complexes, respectively. The ΔG values of the Sc2C2@C2n⊂[12]CPP complexes were considerably smaller (more positive) than those of the Sc2C2@C2n⊂PF[12]CPP complexes, indicating that the [12]CPP molecule exhibited a weaker thermodynamic spontaneity of endofullerenes binding than PF[12]CPP. According to the ΔH values, all host–guest binding reactions were highly exothermic. The free molecules decreased by half after the formation of the host–guest complexes, and correspondingly, the entropies of the six complexes decreased by −240.81 to −266.00 J mol−1 K−1. All of this thermodynamic information indicates that the binding of Sc2C2@C2n to [12]CPP and PF[12]CPP in the vacuum was enthalpy-driven and entropy-opposed.
It is well known that the possibility of intermolecular binding depends entirely on the Gibbs free energy (ΔGb), and it is closely related to environmental temperature. At a given temperature, a negative and positive ΔGb correspond to thermodynamically spontaneous and infeasible complex formation, respectively. Therefore, further research on ΔGb is necessary. In this work, we calculated the ΔGb of all of the complexes and decomposed it as ΔGb = ΔHb − TΔSb, where ΔHb and ΔSb represent the changes in enthalpy and entropy upon binding, respectively, to illustrate the thermodynamic nature of the formation of the complex (Fig. 4).
The ΔHb curve shows that the binding of endofullerenes to [12]CPP and PF[12]CPP is primarily driven by a significant reduction in enthalpy, and the extent of which is less dependent on temperature. The formation of the host–guest complexes reduced the motion freedom of the system, so that the entropy of the system was significantly decreased, resulting in a clear increase in the −TΔSb terms. Overall, the ΔGb value of all of the complexes increased with increasing temperature. Supramolecular complexes Sc2C2@C80⊂[12]CPP, Sc2C2@C82⊂[12]CPP, and Sc2C2@C84⊂[12]CPP were spontaneously formed below 700, 740, and 810 K, respectively, because the calculated ΔGb of the system was negative in this temperature range. Sc2C2@C80⊂PF[12]CPP, Sc2C2@C82⊂PF[12]CPP, and Sc2C2@C80⊂PF[12]CPP were spontaneously formed below 880, 890, and 980 K, respectively. The ΔGb values of the Sc2C2@C2n⊂PF[12]CPP systems were higher than those of the Sc2C2@C2n⊂[12]CPP systems, indicating that the thermodynamic spontaneity of the Sc2C2@C2n⊂PF[12]CPP system was stronger than that of Sc2C2@C2n⊂[12]CPP, as previously described.
Fig. 5 The frontier molecular orbital and the corresponding energy level diagram of the Sc2C2@C2n (n = 40, 41, and 42) and six complexes. |
To further investigate the electronic properties of the complexes, the ionization potential (IP) and the electron affinity (EA) of individual molecules and their complexes were determined from the energies of the corresponding charged and neutral species. Table 3 presents the computed IP and EA values of Sc2C2@C2n, [12]CPP, PF[12]CPP, and their complexes. The relative trends of the fundamental gap (i.e., IP minus EA, GapIP−EA) of the systems were well consistent with those of the frontier orbital gap (GapHOMO–LUMO, Table 3). Either the IP or EA of free PF[12]CPP was higher than that of the free [12]CPP, which indicates that [12]CPP exhibited a stronger electron-donating capacity but weaker electron-accepting capacity than PF[12]CPP in their isolated states. Additionally, the free host molecules [12]CPP and PF[12]CPP exhibited higher IPs but lower EAs than the free guest molecules Sc2C2@C80, Sc2C2@C82, and Sc2C2@C84. This indicates that the free endofullerenes featured stronger electron-donating and electron-accepting capacities than the free [12]CPP and PF[12]CPP. Therefore, the fundamental gaps and frontier orbital gaps of [12]CPP and PF[12]CPP were distinctly larger than those of the three metallofullerenes. A careful examination of the results indicates that the electrons in the Sc2C2@C82 were more bound than those in Sc2C2@C80 and Sc2C2@C84; however, Sc2C2@C82 had a weaker electron-capturing capacity, as its EA was 2.38 eV, which was 0.14 and 0.18 eV lower than those of Sc2C2@C80 and Sc2C2@C84, respectively.
Systems | IP | EA | GapIP−EA | GapHOMO–LUMO |
---|---|---|---|---|
[12]CPP | 6.86 | 0.50 | 6.36 | 5.58 |
PF[12]CPP | 8.68 | 1.58 | 7.10 | 6.30 |
Sc2C2@C80 | 6.61 | 2.52 | 4.09 | 3.05 |
Sc2C2@C82 | 6.83 | 2.38 | 4.45 | 3.42 |
Sc2C2@C84 | 6.66 | 2.56 | 4.10 | 3.11 |
Sc2C2@C80⊂[12]CPP | 6.29 | 2.39 | 3.90 | 3.04 |
Sc2C2@C82⊂[12]CPP | 6.51 | 2.24 | 4.27 | 3.42 |
Sc2C2@C84⊂[12]CPP | 6.44 | 2.45 | 3.99 | 3.09 |
Sc2C2@C80⊂PF[12]CPP | 6.52 | 2.68 | 3.84 | 3.05 |
Sc2C2@C82⊂PF[12]CPP | 6.79 | 2.51 | 4.28 | 3.40 |
Sc2C2@C84⊂PF[12]CPP | 6.69 | 2.66 | 4.03 | 3.09 |
As mentioned earlier, the frontier molecular orbitals of the six complexes were derived from those of endofullerenes and independent of [12]CPP and PF[12]CPP; that is, the IP or EA of the complexes was determined mainly by the concentration of endofullerenes Sc2C2@C2n. As expected, the IP and EA trends of the complexes were similar to those of the endofullerenes (Table 3). Among the Sc2C2@C2n⊂[12]CPP and Sc2C2@C2n⊂PF[12]CPP systems, Sc2C2@C82⊂[12]CPP and Sc2C2@C82⊂PF[12]CPP respectively exhibited both the highest IP and the lowest EA. The IP and EA of the Sc2C2@C2n⊂PF[12]CPP systems were higher than those of the Sc2C2@C2n⊂[12]CPP systems, which indicates that the Sc2C2@C2n⊂[12]CPP system possessed a stronger electron-donating capacity but weaker electron-accepting capacity than the Sc2C2@C2n⊂PF[12]CPP system.
Charge transfer commonly occurs in supramolecular systems. Charge transfer results in the additional stability of a complex in terms of interaction energy. In addition to the IP of the donor and the EA of the acceptor, the intermolecular distance and mutual orientation also highly influence the amplitude of charge transfer.41 Table 4 presents the amount of charge transfer as determined via Mulliken population analysis and natural population analysis (NPA). The Mulliken and NPA charge transfers of the Sc2C2@C2n⊂[12]CPP systems were around 0.07 and 0.04 e, respectively. The charge transfer amount of the Sc2C2@C2n⊂[12]CPP system as derived via Mulliken analysis was consistent with the charge transfer of the [11]CPP–La@C82 complex reported by Yamago et al.14 Sc2C2@C80⊂PF[12]CPP, Sc2C2@C82⊂PF[12]CPP, and Sc2C2@C84⊂PF[12]CPP exhibited Mulliken and NPA charge transfers of 0.19 and 0.10 e, 0.17 and 0.09 e, and 0.35 and 0.12 e, respectively, which were significantly greater than those of the Sc2C2@C2n⊂[12]CPP system. These results suggest that the amount of charge transfer between the host and guest molecules was more dependent on the species of the host molecule rather than the species and size of the fullerene.
Host–guest | Mulliken | NPA |
---|---|---|
Sc2C2@C80⊂[12]CPP | 0.08 | 0.04 |
Sc2C2@C82⊂[12]CPP | 0.07 | 0.04 |
Sc2C2@C84⊂[12]CPP | 0.07 | 0.04 |
Sc2C2@C80⊂PF[12]CPP | 0.19 | 0.10 |
Sc2C2@C82⊂PF[12]CPP | 0.17 | 0.09 |
Sc2C2@C84⊂PF[12]CPP | 0.35 | 0.12 |
Fragments | Eelest | Eind | Edisp | Eexch | ESAPT0 |
---|---|---|---|---|---|
[12]CPP-Sc2C2@C80 | −29.11 (22.19%) | −6.89 (5.25%) | −95.21 (72.56%) | 65.47 | −65.74 |
[12]CPP-Sc2C2@C82 | −29.42 (22.40%) | −6.82 (5.19%) | −95.09 (72.41%) | 64.96 | −66.37 |
[12]CPP-Sc2C2@C84 | −26.89 (21.11%) | −6.35 (4.99%) | −94.14 (73.90%) | 62.87 | −64.51 |
PF[12]CPP-Sc2C2@C80 | −30.55 (25.26%) | −6.33 (5.23%) | −84.06 (69.51%) | 73.59 | −47.35 |
PF[12]CPP-Sc2C2@C82 | −37.75 (27.68%) | −7.52 (5.51%) | −91.10 (66.80%) | 84.78 | −51.59 |
PF[12]CPP-Sc2C2@C84 | −66.75 (32.60%) | −13.38 (6.46%) | −126.25 (60.94%) | 167.61 | −39.57 |
Furthermore, for all of the six complexes, Edisp dominated the attractive component of the interaction energy, and for most of these cases, it even contributed more than 60% to this part. Even though electrostatic interaction is not as important as dispersion interaction in the formation of dimers, its effect on binding should not be ignored, with Eelest accounting for more than 20% of all complexes. As depicted in Table 4, a small amount of charge transfer existed between the host and guest molecules, and thus, Eind contributed the least to the interaction, only about 5%.
IGMH is a very intuitive method proposed to graphically depict the main interaction regions between specific fragments of compounds by using the isosurface of δginter function. Fig. 6 delineates the weak interaction regions of the host–guest complexes and the role played by each atom in the fragment. The atoms in the different fragments are colored according to their contributions to the dispersion interaction energy. The larger the contribution to the dispersion effect, the closer the bluer the hue. For all the complexes, the extended δginter isosurface in green suggests that the interaction between the host and guest molecules covered a wide spatial region and that the electron density on δginter isosurface was small (about 0.001 a.u). Thus, the interaction between host and guest molecules can be considered typical π–π stacking.7 The weak interaction region of the Sc2C2@C2n⊂[12]CPP system exhibited an open weak interaction region. The closer the interfacial distance (di) between the host [12]CPP and guest Sc2C2@C2n was to 3.4 Å (Fig. 3), the bluer the color of the atoms, indicating stronger π–π vdW interactions. In contrast, at an inter-atom distance considerably greater than 3.4 Å (corresponding to white atom color, the longest distance was >6 Å), vdW π–π interactions will no longer exist. Compared with the Sc2C2@C2n⊂[12]CPP system, the Sc2C2@C2n⊂PF[12]CPP system exhibited an almost closed weak interaction region, consistent with the relative trend of the binding energies. This trend was due to the relatively uniform distribution of the interfacial distance between PF[12]CPP and the fullerenes. IGMH not only visualizes the weak interaction regions but also identifies the weak interaction type as π–π stacking.
(1) ESP analysis showed that [12]CPP exhibited negative regions associated with phenyl rings and positive regions on the rim of the ring. However, PF[12]CPP exhibited positive regions in the internal region of the ring and negative regions on the ring rim because the polarity of F atoms was much greater than that of H atoms. In addition, the ESP can be useful for identifying the stabilized structure of the complexes.
(2) Three endofullerenes, namely Sc2C2@C80, Sc2C2@C82, and Sc2C2@C84, established favorable interactions with all of the phenylene units of the Sc2C2@C2n⊂[12]CPP system. For the Sc2C2@C2n⊂PF[12]CPP system, the three endofullerenes established favorable interactions with the F atoms of the ring units. This could be related to the greater polarity of the endofullerenes, which showed positive ESP regions around the carbon atoms contacting the endohedral units.
(3) With increasing carbon cage size, the binding energy of the complexes also increased, possibly because the larger carbon cage allowed for a large region of π–π contact between the host molecules ([12]CPP and PF[12]CPP) and guest molecules (Sc2C2@C2n). However, PF[12]CPP could more suitably accommodate Sc2C2@C2n than [12]CPP. Thermodynamic information indicated that the complexes of the two systems were enthalpy-driven and entropy-opposed. Moreover, the Sc2C2@C2n⊂PF[12]CPP system exhibited a stronger thermodynamic spontaneity than the Sc2C2@C2n⊂[12]CPP system.
(4) The frontier molecular orbitals of the six complexes were completely derived from those of the endofullerenes and independent of [12]CPP and PF[12]CPP. The IP and the EA showed that the electron-accepting capacities of the complexes were higher than those of the free [12]CPP and PF[12]CPP. Moreover, compared with [12]CPP, PF[12]CPP exhibited a stronger charge transfer with fullerenes.
(5) Energy decomposition analysis showed that the interaction between the host and guest molecules of the complexes was dominated by dispersion attraction. However, the intermolecular electrostatic interaction was considerable. The visualized weak interaction regions indicated that π–π vdW interaction would be nonexistent in the region with di much greater than 3.4 Å.
The elucidated structures and properties of those host–guest supramolecules, which are difficult to experimentally obtain, and the reported conclusions can provide useful theoretical guidance for materials chemistry applications.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2ra08153a |
This journal is © The Royal Society of Chemistry 2023 |