Gloria
Tabacchi
a,
Gion
Calzaferri
b and
Ettore
Fois
*a
aDepartment of Science and High Technology, University of Insubria, and INSTM, Via Valleggio 9, I-22100 Como, Italy. E-mail: ettore.fois@uninsubria.it
bDepartment of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
First published on 26th July 2016
Confined supramolecular architectures of chromophores are key components in artificial antenna composites for solar energy harvesting and storage. A typical fabrication process, based on the insertion of dye molecules into zeolite channels, is still unknown at the molecular level. We show that slipping of perylene diimide dyes into the one-dimensional channels of zeolite L and travelling inside is only possible because of steric-interaction-induced cooperative vibrational modes of the host and the guest. The funnel-like structure of the channel opening, larger at the entrance, along with a directionally asymmetric entrance–exit probability, ensures a favorable self-assembly process of the perylene units.
Scheme 1 Left: Sequential loading, leading after steps 1, 2, and 3 to the organized 3 color composite 4. Right: PDI with vdW diameter.1 |
Progress in theoretical modeling of zeolites has been enormous, fortunately.7 We unraveled the interaction and orientation of dyes inside ZL and we explained how stopper molecules irreversibly modify ZL at the channel entrances.8 Now, to investigate the slipping of a PDI dye across those entrances we build a first-principles model featuring the channel openings of ZL (see the ESI†), as previously done for the ZL entrance functionalization.8e Then, we select a symmetric PDI (Scheme 1, R = CH3) and position it right outside the opening, with its longest molecular axis parallel to the channel axis, as depicted in Fig. 1 – for PDI, this is clearly the only possible way to cross the channel entrance. The high temperatures needed to obtain PDI–ZL adducts indicate that the PDI loading is an activated process, with timescales normally not accessible to molecular simulations. To explore how PDIs penetrate through the channel, we need to augment the finite temperature sampling of first principles molecular dynamics9a using the metadynamics approach.9b To achieve this aim, we use as the collective variable (CV) the average distance of the PDI core atoms from the ZL entrance ring,9c,dFig. 1a. After 5.0 ps equilibration at a temperature of 200 °C, the metadynamics simulation starts with the PDI molecule in front of the entrance, as in Fig. 1a, and evolves as featured in the Movie (ESI†), to end up with the free energy profile shown in Fig. 1b.
Such a profile is characterized by a nearly flat initial region, followed by a deep well; then, a free energy increase occurs, followed by an even deeper well. This means that the transit of PDIs through the channel opening is a directionally asymmetric process:10 the entrance is favoured over the exit, and the molecule gets finally encapsulated into the porous matrix.
The first and nearly flat region of the curve (CV: from 0 Å to −5 Å, E ≈ −9 kcal mol−1) corresponds to the initial stages of the PDI encapsulation, Fig. 2a and b: only the relatively smaller and more polar part of the molecule, the diimide fragment, has passed through the channel entrance. In the first minimum region (CV ≈ −8 Å, E ≈ −48 kcal mol−1), roughly half of the molecule has slipped inside the channel: two of the PDI carbonyl groups interact with K+ cations of the first (and outermost) ZL cell, while the other two are still outside of the opening, Fig. 2d. In the second and deeper minimum (CV ≈ −16 Å, E ≈ −78 kcal mol−1), the PDI is fully encapsulated, and its carbonyl groups interact with K+ cations both in the 1st and the 2nd cell of the zeolite, Fig. 2f. Conversely, in the free energy maximum (CV ≈ −12 Å, E ≈ −4 kcal mol−1) the molecule moves from the 1st cell to the 2nd one and none of its carbonyl groups are involved in interactions with the zeolite K+ cations, Fig. 2e. The carbonyl–K+ interactions,8a therefore, trigger the slipping and may also rule the diffusion of PDIs inside the ZL channel system because of their substantial stabilizing effect: the interaction of a K+ cation with a PDI carbonyl oxygen amounts to 26.1 kcal mol−1. Another possible factor favouring the PDI inclusion could be a charge–π interaction, with the zeolite K+ cation positioned on top of the aromatic PDI plane. This interaction, although much weaker (9.9 kcal mol−1 (ESI,† Fig. S3)), may play a key role when the K+ cations are not bound to the PDI carbonyls, thus facilitating the transit of the molecule from the first to the second ZL cage (Fig. 2e).
Fig. 2 (a–f) Clockwise from the top left. They represent 6 relevant instantaneous configurations taken from the metadynamics simulation; the corresponding values of the CV parameter are in parentheses (see also the Movie, ESI†). Color code: The Si and Al of the framework are gray sticks. C, cyan; H, white; N, blue; carbonyl O, red; K+, purple. Coordination of K+ to the carbonyl oxygen atoms is highlighted by yellow sticks. |
The question then arises as to whether the porous inorganic framework is an active player in the slipping process.11 To find an answer, we take a closer look at the ZL structural changes as a function of the PDI penetration. The channel opening, defined by 2 sets of crystallographically different oxygens O1 and O2, is a 12-membered (12M) ring,6 which is slightly narrower along the O1–O1 directions (Fig. 3, left). Hence, the entrance process of a spherical object is controlled by the O1–O1 diameter. The average value of this quantity calculated for the entrance ring is larger compared to the ring positioned just below the entrance (Table 1, Fig. 3 and 4, and ESI,† Table S2). This result suggests a funnel-like structure of the ZL channel opening, which might well facilitate the PDI entrance, especially under the high temperature conditions adopted in the inclusion experiments.5d
Fig. 4 Evolution of the O1–O1 (a) and the O2–O2 (b) distances (in Å) during the PDI slipping in the metadynamics trajectory at 200 °C. From top to bottom, the distances refer to the 1st (the channel opening), the 2nd and the 3rd 12M ring (Fig. 3). In (b), the blue lines refer to the O2–O2 diameters along which the PDI slipping occurs (Fig. 5). |
Ring | O1–O1 | O2–O2 |
---|---|---|
1st | 10.338 (0.764) | 10.649 (0.795) |
2nd | 10.048 (0.735) | 10.547 (0.785) |
3rd | 10.067 (0.737) | 10.501 (0.780) |
∞ | 10.015 (0.732) | 10.475 (0.778) |
The temperature effects on the 12M ring opening are actually remarkable, and can be appreciated by looking at the variation of the O1–O1 and O2–O2 diameters along the 200 °C trajectory, depicted in Fig. 4. All distances show large temperature induced oscillations, of the order of 0.1 nm. As 200 °C corresponds to 329 cm−1, these oscillations should be ascribable mainly to the bending motion of the Si–O–Si (Al–O–Si) and O–Si–O (O–Al–O) angles (typical frequencies below 400 cm−1) and to collective vibrational modes of lower energy.12 What is to notice is that the oscillations do not depend solely on thermal effects, but are correlated with and driven by the slipping PDI. During the passage across the 1st and 2nd 12M rings, the O2–O2 distances which the PDI is aligned to show both the greatest average values and the largest amplitude oscillations, indicating that the dye, due to its size, has jostled its way to slip inside. This fact clearly emerges by observing that the structural changes of the 1st and 2nd rings along the process are indeed quite different from those of the 3rd and innermost ring, Fig. 4. Because the PDI never crossed that ring in the sampled configurations, the thermal oscillations are equally distributed among its O1–O1 and O2–O2 distances (Fig. 4a and b, bottom panels). On the other hand, in the first two rings, which are crossed by the transiting molecule, a large elongation of one O2–O2 diameter implies a shortening of the other diameters, hence an elliptic deformation of the ring. We therefore conclude that the channel gets locally deformed through a PDI-induced redistribution of the thermally activated – and otherwise random – oscillations of the zeolite framework structure (Fig. 5).
We have shown that the size and structure of PDI greatly influence the ZL framework oscillations, which are essential for the dye encapsulation. A closer inspection of the trajectory reveals that the vibrations of the PDI molecule also actively contribute to this process. Initially, the PDI enters the channel with its short axis oriented along one of the O1–O1 diameters, as shown in Fig. 2a and b, because the K+ sites facing the channel wall are located just below the O1 sites, Fig. 3. Due to the strong interaction with the PDI's carbonyl groups, these K+ cations draw the molecule to the entrance; however, then the dye reorients, and eventually slips inside the channel through the widest passage, i.e., along one of the O2–O2 diameters, Fig. 2c and d.
Such a molecular reorientation process, instrumental in the PDI insertion, is very far from being a rigid rotation of the aromatic backbone: it occurs, indeed, via an out-of-plane distortion of the molecule, Fig. 2c. This distortion corresponds just to the lowest energy normal mode of the PDI (ESI,† Fig. S4). Altogether, our results therefore indicate that both ZL and PDIs exploit their low energy vibrational modes, cooperatively, for reaching the innermost and deepest free energy minimum, thus accomplishing the encapsulation process.
In conclusion, the answer to the question of how PDIs can slip into ZL channels and travel inside despite their apparently too large width can be summarized by the following points: (i) the funnel-like shape of the channel opening does facilitate the initial phase of the PDI entrance; (ii) the stabilizing interactions among the PDI carbonyl oxygens and the ZL K+ cations provide the main energetic driving force; (iii) the cooperative vibrational modes of the host and the guest play so well to enable these processes, via a local elliptic deformation of the channel; and (iv) the asymmetry in the free energy profile clearly indicates that the PDI entrance process is favoured over the PDI exit.
Overall, the self-assembly of PDI–ZL adducts is possible, at high temperatures, because of the thermodynamically induced unidirectional motion of PDIs. This motion is combined with cooperative vibrational modes, which allow overcoming the barrier imposed by the structure of the channels. It is probable that similar mechanisms could govern the build-up and behaviour of other classes of host–guest composites.13
The present work was supported by the Italian MIUR through the projects ImPACT (FIRBRBFR12CLQD), INFOCHEM (PRIN2010CX2TLM_006), and Insubria University FAR2014. The CINECA supercomputing center (Italy) is acknowledged for computing time (ISCRA project DENARI, HP10B6JTDP).
Footnote |
† Electronic supplementary information (ESI) available: Computational set-up, tables and figures related to the computations, further details and a movie of the penetration process. See DOI: 10.1039/c6cc05303c |
This journal is © The Royal Society of Chemistry 2016 |