Open Access Article
Shuhai
Qiu†
a,
Li
Zhang†
a,
Da-Hui
Qu
*a and
Zhaohui
Wang
*ab
aLaboratory for Advanced Materials and Joint International Research Laboratory of Precision Chemistry and Molecular Engineering, Feringa Nobel Prize Scientist Joint Research Center, Frontiers Science Center for Materiobiology and Dynamic Chemistry, Institute of Fine Chemicals, School of Chemistry and Molecular Engineering, East China University of Science and Technology, Shanghai, 200237, P. R. China. E-mail: dahui_qu@ecust.edu.cn; wangzhaohui@mail.tsinghua.edu.cn
bKey Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, P. R. China
First published on 26th September 2025
Cycloalkanoid molecules are of fundamental interest due to their unique geometric structures and physical properties. However, studies on the molecular dynamics of these molecules tend to be difficult owing to labile conformations. Herein, we report a series of “supercycloalkanes”, namely dihydropyrazine-embedded macrocycles, that show puckered structures imitating those of cyclobutane, cyclopentane and cyclohexane, as confirmed by crystal structures. 1H NMR spectrum analyses and theoretical calculations reveal that the macrocycle containing four dihydropyrazine corners undergoes a ring-flapping process in solution, and the one with five dihydropyrazine corners involves a pseudorotation process. Photophysical results indicate that these macrocycles show donor (D)–acceptor (A)-type structures, resulting in solvent-dependent emissions. Moreover, each dihydropyrazine corner could effectively interact with fullerene molecules owing to the V-shaped π-surface. This study provides a series of cycloalkanoids with a “superstructure” and discloses the molecular dynamics of cycloalkanoid macrocycles.
Cycloalkanes, which contain a cyclic pure carbon skeleton, have been well studied by chemists in the early stages.11 One of the hottest topics around these molecules focuses on molecular conformations and dynamics.12,13 Due to the ideal angle of 109.5° among sp3-hybridized carbon atoms, cycloalkanes tend to adopt puckered geometries to alleviate strain. Specifically, cyclobutane and cyclohexane have the lowest energy conformer with a folded and chair geometry, respectively, while cyclopentane has two types of conformers with nearly equal energy, the Cs symmetric envelope conformation with four carbon atoms in a plane and the C2 symmetric half-chair form with three carbon atoms in a plane.14,15 Theoretical calculations and infrared spectroscopic measurements revealed that cyclobutane underwent an inversion process with an energy barrier of less than 1.5 kcal mol−1 (Fig. 1b),16–19 while cyclopentane underwent a pseudorotation process between the envelope and half-chair conformers (Fig. 1c).20–22 However, the labile conformations of cyclobutane and cyclopentane preclude in-depth experimental studies, and only cyclohexane and higher cycloalkanes could be directly investigated using nuclear magnetic resonance (NMR) spectroscopy.23,24 Therefore, cyclic molecules with sufficient conformational stability as well as similar molecular dynamics to cycloalkanes should be of interest. In the past decades, pioneering studies on cyclic oligomers constituted by ortho-linked aromatics, including benzene,25,26 thiophene,27 pyrene,28 and porphyrins,29 have been reported. Among them, most tetramers adopt a folded conformation in the solid state similar to that of cyclobutene, while higher aromatic oligomers, such as pentamers and hexamers, showing dynamic behaviors that resemble those of cyclopentane and cyclohexane, remain largely unexplored. Inspired by the less strained cycloalkanes, we envision that the installation of adaptive π-corners into the rigid skeleton of cyclo-para-phenylenes (CPPs)30–33 should circumvent the high strain by regulating the corner angles and rotating around the diphenylene linkers; therefore, these polygonal macrocycles can act as “supercycloalkanes” that are structurally akin to cycloalkenes and display interesting molecular dynamics (Fig. 1d).
1,4-Dihydropyrazine is an 8π-electron antiaromatic ring,34 and its π-extended derivatives, such as dihydrophenazines, have shown intriguing properties and functions, such as photoredox catalysis, tunable emissions, redox activity, and conformational dynamics.35–38 On the one hand, by virtue of the antiaromatic character of dihydropyrazine, the hexagonal ring tends to adopt a V-shaped geometry in the ground state, especially when sterically hindered polycyclic aromatic units are fused.39,40 On the other hand, the introduction of nitrogen atoms allows for the equatorial and axial exchanges on the substituents, providing more flexibility for conformational relaxation. Overall, these factors make dihydropyrazine derivatives ideal candidates as corners to construct π-conjugated polygonal macrocycles. Herein, we report the one-pot synthesis of dihydropyrazine-embedded macrocycles uniformly containing four to six dipyrene-fused dihydropyrazine corners, respectively. Instead of forming ideal polygonal tube-like structures, these macrocyclic molecules adopt puckered geometries as the lowest energy conformations similar to those of cycloalkanes, as unambiguously manifested by X-ray structures. The conformational analyses of these cyclic molecules are also studied using 1H NMR spectra and theoretical calculations to obtain insight into the molecular dynamics of these flexible macrocycles. Moreover, these macrocycles exhibit intriguing properties, including D–A type electronic structures, and supramolecular interactions with fullerenes.
Single crystals suitable for X-ray diffraction analyses were grown through slow evaporation of methanol into a chloroform solution for DPP-Q and DPP-H and a mesitylene solution for DPP-P. As shown in Fig. 2a, DPP-Q showed a puckered geometry with C2 symmetry, slightly deviating from the ideal D2d symmetric geometry for cyclobutane. The edge length of the folded structure, that is the diphenyl diamine linker, is around 10.0 Å. The distances between two centroids of the diagonal 1,4-dihydropyrazine rings are 12.7 and 14.9 Å, respectively. Correspondingly, the bending angles are measured to be 75.0° and 86.4°. Such a folded geometry could efficiently alleviate the strain, which is in line with the calculated strain energy of 6.1 kcal mol−1 according to hypothetical homodesmotic reactions (Scheme S1). This value is significantly lower than that of [12]CPP (48.1 kcal mol−1),42,43 indicating the high conformational flexibility of 1,2-dihydropyrazine corners. Interestingly, similar to sp3-hybridized carbon atoms in cyclobutene, each dipyreneo[c,e]dihydrophenazine corner of DPP-Q has one pyrene unit at the pseudoequatorial position and the other one at the pseudoaxial position (Fig. 2d). Similarly, DPP-P and DPP-H also adopted puckered configurations to minimize the strain energy in the crystalline state. Specifically, only the Cs symmetric envelope configuration was observed for DPP-P, even though the C2 symmetric half-chair conformation was reported to have a low energy close to that of the envelope form in the case of cyclopentane.15DPP-H adopted a chair-shaped geometry, in line with the lowest-energy conformation of cyclohexane.44,45 The bending angle was measured to be 40.5°for DPP-P and 58.7° for DPP-H, respectively. Correspondingly, the strain energy was calculated to be 8.6 kcal mol−1 for DPP-P and 2.3 kcal mol−1 for DPP-H. Notably, the strain energy is closely correlated with the angle at each dihydropyrazine corner (θ), which can be defined as the sum of the angles between the N–N axis and phenyl groups at the N atoms (θa and θb, Fig. 2). The corner angles θ of both DPP-Q and DPP-H are in the range of 227.2–235.8°, which are slightly larger than that of the strainless precursor DPP-2Br (226.5°).41 This is in agreement with the small strain energy of DPP-Q and DPP-H. However, in DPP-P, both the corner angles on the short edge (12.4 Å) of the trapezoid plane are 247.1°, much larger than that on the long edge side (θ = 223.2°). Such a large angle deviation from that of the strainless DPP unit leads to a higher strain than that of DPP-Q and DPP-H. As a result, the diphenyl linker in the short edge adopts a bent geometry, which is reminiscent of that of highly strained [n]CPPs.30
The molecular structures of DPP-Q, DPP-P and DPP-H in solution were investigated by variable-temperature (VT) 1H NMR measurements. Considering the folded structure of DPP-Q in the crystalline state, two sets of 1H NMR signals deriving from the pyrene units in both the pseudoequatorial and pseudoaxial positions were expected to appear. Interestingly, only one set of proton signals (δ = 8.06, 8.20 and 8.67 ppm) belonging to the pyrene unit of DPP-Q was observed in the aromatic region at 283 K (Fig. 3a), indicating that the protons of pyrene subunits at the pseudoequatorial positions are equal to those at the pseudoaxial positions. These results imply that the folded conformers probably undergo a fast inversion process to afford a time-average D4h symmetric structure. As the temperature decreased from 283 K to 188 K, the signals gradually broadened due to a slowdown of conformational motions. When the temperature was lowered to the instrumental temperature limit of 173 K, each singlet signal belonging to the protons of pyrene subunits fully divided into two sets of peaks, suggesting that the inversion process is completely restricted, and the molecular structure adopts a folded D2d symmetric geometry, in which the protons belonging to the pyrenes at the pseudoequatorial positions differ from those at the pseudoaxial positions. Based on line-shape analysis of proton signals, the exchange rate constant k at different temperatures was estimated according to the reported method developed by Gasparro and Kolodny46 (see details in the SI). By fitting the parameter ln(k/T) versus the reciprocal number of temperature (1/T) to the Eyring equation (Fig. 3a, inset), the thermodynamic parameters ΔH* = 8.33 kcal mol−1 and ΔS* = −2.54 kcal mol−1 K−1 were obtained. Accordingly, the inversion energy barrier (ΔG*) at the coalescence temperature (188 K) was determined to be 8.81 kcal mol−1. This value matches well with theoretical calculations, which indicated that the inversion of the folded conformers overcame an energy barrier of 9.2 kcal mol−1via a coplanar D4h symmetric structure in the transition state (Fig. 3c). It is worth noting that this energy barrier is significantly higher than that of cyclobutane (1.45 kcal mol−1),19 which involves a rapid inversion process beyond the measurable range by NMR. Similar to DPP-Q, DPP-P also exhibited one set of proton signals belonging to the pyrene subunits (Fig. 3b), suggesting fast interconversion amongst all the conformers at room temperature. Further lowering the measurement temperature led to broadening of the proton signals, indicating the increased restrictions of conformational changes. Unfortunately, the conformational interconversion still occurred even at the temperature limit (173 K) as no signal split was observed. These results suggest that the conformational interconversion process of DPP-P has a relatively low energy barrier, lower than 5 kcal mol−1, which is also reminiscent of the pseudorotation process of cyclopentane. In contrast to the inversion process with a high energy barrier of 5 kcal mol−1,47 pseudorotation of cyclopentane generally proceeds without passing through a significant energy barrier, and the energy difference between the half-chair and envelope conformations is quite small with the envelope form being more stable by 0.5 kcal mol−1.48 However, owing to the labile conformations, experimental investigations on the pseudorotation process are still difficult.
In order to further understand the conformational dynamics of DPP-P, the conformational interconversion process was investigated by theoretical calculations. First, a ring-flipping process similar to that of DPP-Q was proposed (Fig. 3d). The D5h symmetric transition state structure was optimized by relaxed scan to confine the N–N axis of 1,4-dihydropyrazine rings in coplanar positions. However, the high energy barrier of 24.4 kcal mol−1, which is inconsistent with the results obtained from VT 1H NMR measurements, precluded this hypothesis. Therefore, we assumed that the conformational interconversion of DPP-P might undergo a pseudorotation process, similar to that of cyclopentane.49 Optimization of the molecular structures of DPP-P indicates that the energy of the envelope configuration is slightly higher than that of the half-chair one by 0.4 kcal mol−1, which is in line with the slight energy difference between the two conformers in the case of cyclopentane. These results indicate that, instead of the inversion among the envelope conformers, a pseudorotation process involving both envelope and half-chair conformers is more feasible for DPP-P. In the case of DPP-H, multiple signals in 1H NMR spectra were observed at room temperature (see Fig. S18), implying a high energy barrier in the conformational interconversion process compared to DPP-Q and DPP-P. However, further investigations on the molecular dynamics were impeded due to poor solubility even at a high temperature up to 413 K (Fig. S19).
To capture the conformational motions of these molecules, X-ray snapshots of different conformations were obtained by growing single crystals under various conditions. Single crystals of DPP-Qvia slow evaporation of methanol into a 1,2-dichloroethane solution revealed a tightened conformation compared to that from the chloroform solution (Fig. 3e). The distances between two centroids of the diagonal 1,4-dihydropyrazine rings were shortened by 1.5–3.6 Å, and the bending angle was decreased by ∼33°, which further supports a ring-flipping process as proposed in Fig. 3c.
The optical properties of these macrocyclic molecules were explored in solution (Fig. 4). Three molecules exhibited nearly identical absorptions in dichloromethane with the absorbance maximum at ∼387 nm, which is also commonly observed in [n]CPPs.50 All compounds exhibited green emission at different wavelengths (DPP-Q: 472 nm, DPP-P: 487 nm, and DPP-H: 476 nm), and the fluorescence quantum yields were measured to be 5.8% for DPP-Q, 6.8% for DPP-P and 5.3% for DPP-H. The fluorescence decay time determined using the nanosecond time-correlated single photon counting technique revealed a short lifetime of 2.0 ns for these molecules in dichloromethane (Fig. S7). The weak fluorescence and short emission lifetimes could be attributed to the fast conformational changes in solution. To further support this explanation, variable temperature fluorescence of DPP-Q and DPP-P in toluene was measured by gradually lowering the temperature from room temperature to 78 K (Fig. S10). As the temperature was lowered to 113 K, only a small blue-shift in the emission (∼5 nm for DPP-Q and ∼10 nm for DPP-P) as well as limited enhancement of intensity was observed, which could be explained by the restrictions on the conformational inversions of the puckered structures, as supported by VT 1H NMR results. When the temperature was further decreased to 78 K, the emission intensity dramatically increased, and a distinct blue-shift (∼25 nm for DPP-Q and ∼17 nm for DPP-P) was observed. These phenomena indicated a complete suppression of the conformation vibration of the DPP unit. These results also are consistent with reported π-extended dihydrophenazines, which involved bent-to-planar excited-state dynamics upon photoexcitation.37,51
The optical properties of these macrocycles in different solvents were also investigated. By modulating the polarity of the solvent, the absorptions showed negligible changes (Fig. S8). In contrast, the emission was distinctly red-shifted by 40 nm for DPP-Q, 30 nm for DPP-P and 52 nm for DPP-H (Fig. S9), which suggests that an intramolecular charge transfer process occurs upon excitation. These results are further supported by the disjoint distributions of the calculated frontier molecular orbitals (Fig. 4b–c, e, f, h and i). The highest occupied molecular orbitals (HOMOs) of DPP-Q localize on the whole dihydropyrazine-embedded CPP subunit, while the lowest unoccupied molecular orbitals (LUMOs) distribute on the pyrene units. In comparison to DPP-Q, the HOMOs of DPP-P and DPP-H distribute partially on the N-doped CPP skeleton, and the LUMOs localize on the partial pyrene units.
Considering the large cavity as well as the electron-rich structure, host–guest interactions of these macrocycles with fullerenes (C60 and C70) were studied in solution. Fluorescence titration was performed by gradually adding fullerene molecules into a solution of DPP-Q or DPP-P in toluene (Fig. S11–S12). Fittings by monitoring the fluorescence changes (Fig. 5a–d) suggested a 1
:
1 binding mode for DPP-Q and DPP-P.52 The association constant (Ka) of DPP-Q with C60 was determined to be 2.50 × 104 M−1, while the Ka value with C70 was significantly higher at 1.86 × 105 M−1 in comparison with that with C60. These results indicated a selective complexation of C70 over C60 for DPP-Q. In contrast, DPP-P exhibited comparable Ka values of 3.13 × 104 M−1 with C60 and 2.02 × 104 M−1 with C70, respectively, suggesting a similar binding affinity toward C60 and C70.
In order to clarify the origins of supramolecular interactions with fullerenes, single crystals of the complexes were grown. Fortunately, single crystals of the complex DPP-Q@C60 suitable for X-ray analysis were obtained by slow evaporation of methanol into the DPP-Q/C60 mixture in toluene. As shown in Fig. 5e, one C60 molecule was located at the outside of each corner of DPP-Q with the distances of 3.22–3.35 Å, indicating strong interactions between the V-shaped π-surface and the convex of fullerenes. Notably, no C60 molecules were encapsulated in the cavity of DPP-Q, which might be attributed to the hindrance from the pyrene units at the openings. Therefore, the recognition selectivity of C70 over C60 can be explained by the larger π-surface of C70 molecules. Moreover, four C60 molecules aggregated together among DPP-Q molecules, forming a highly ordered supramolecular assembly (Fig. 5f and S13). In the case of DPP-P, the overcrowded configuration as well as the rapid conformational changes might restrict the interactions with fullerenes, thus affording moderate Ka values. It is worth noting that electron donor–acceptor systems have been proved to be essential for applications including photovoltaics and molecular electronics;53,54 thus, these complexes of the electron-rich macrocycles with fullerenes might hold potential for applications in materials science in the future.
The data supporting this article have been included as part of the SI. Supplementary information: synthetic details, spectroscopic measurements and theoretical calculations. See DOI: https://doi.org/10.1039/d5sc06636k.
Footnote |
| † These authors contributed equally to this work. |
| This journal is © The Royal Society of Chemistry 2025 |