M. Saeed Mirzaei†
a,
Saber Mirzaei†ab,
Hormoz Khosravia and
Raúl Hernández Sánchez
*abc
aDepartment of Chemistry, Rice University, 6100 Main St., Houston, TX 77005, USA. E-mail: raulhs@rice.edu
bDepartment of Chemistry, University of Pittsburgh, 219 Parkman Ave., Pittsburgh, PA 15260, USA
cRice Advanced Materials Institute, Rice Sustainability Institute, and Rice WaTER Institute, Rice University, Houston, TX 77005, USA
First published on 14th August 2025
Access to highly strained molecules remains a challenge. We report the synthesis of a bench-stable and highly strained (∼135 kcal mol−1) deep cavitand (1). Synthesis of 1 follows a two-step protocol from known starting materials. The rigid structure of 1 contains four arch-shaped biphenyelenes with a large cavity (∼500 Å3) to host fullerenes.
![]() | ||
Fig. 1 Highly strained conjugated aromatic structures. (a)–(f) Previous examples reported in the literature. (g) This work, deep cavitand 1. Values shown correspond to the DFT-calculated strain energies. SE in (f) was calculated using StrainViz.29 |
Our group and others have developed a synthetic approach towards contorted macrocycles based on macrocyclic arene species, e.g., resorcin[4]arenes.22–25 In our work, the resulting tubular species develop a sizable strain in the last synthetic step reaching values of ca. 90 kcal mol−1 at poor isolated yields of <1–3%. We hypothesized that low yields could be overcome by moving away from intermolecular reactions in exchange for intramolecular couplings. The success of this approach was demonstrated recently with the synthesis of 1m,26 a non-strained analogue of 1, and has also been exemplified by others in recent reports.27 Following an intramolecular Ni-mediated homocoupling reaction, and comparing it to our previous reports employing intermolecular cross couplings, led us to improve our yields by 20-fold in the strain-forming step.
The synthesis of 1 – a highly strained deep-cavity molecular container – consists on derivatizing 328 through a Suzuki – Miyaura cross coupling reaction with para-chlorophenyl boronic acid resulting in compound 2 in 75% yield (Fig. 2a). The yield obtained for 2 is remarkable as eight C–C bonds are formed. Subsequently, a Ni-mediated Yamamoto cross coupling leads to 1 as a white solid in 23% yield. In comparison, the intermolecular Suzuki–Miyaura cross coupling between 3 and 4,4′-biphenyldiboronic acid bis(pinacol) ester does not produce 1 in detectable quantities (Fig. 2a). Using meta-chlorophenyl boronic acid in this reaction sequence leads to 1m, which is substantially less rigid than 1 (SE = ∼9 kcal mol−1, SI).26
Analysis of 2 via 1H NMR displays six aromatic resonances indicating an ideal C4-symmetric structure in solution (Fig. S1). Formation of 1 results in eight aromatic resonances suggesting that a species of similar symmetry to 2 is formed yet with additional intricacies (Fig. 2b). This observation led us to propose that phenylene ring “P” cannot freely rotate making resonances “a” and “c” no longer equivalent as in 2; similarly, “b” and “d” are unique resonances in 1. COSY and NOESY NMR supports the assignment of 1 (Fig. S5–S7). Note that MALDI-MS of 1 match perfectly the simulated pattern of its cationic molecular ion peak, [M]+ (Fig. 2c). Variable temperature 1H NMR of 1 in 1,1,2,2-tetrachloroethane-d2 (TeCA-d2) from 25 to 115 °C showcases non-coalescing resonances confirming the rigid nature of “P” (Fig. S8). DFT calculations indicate a rotational barrier of ∼25 kcal mol−1 for “P”, similar to other rigid systems.30
The top-rim of 1 contains twelve phenylenes connected in a unique and alternating 2:
1 para
:
meta fashion. We did not observe decomposition, color change, or the formation of new products during the purification of 1 using silica gel or after prolonged time on the benchtop, as shown through MALDI-MS and 1H NMR. In contrast, other strained systems are known to decompose over time.9,16 We grew colorless crystals by slow diffusion of MeCN into a chlorobenzene solution of 1. The structure of 1 is shown in Fig. 2d providing a clear view of the biphenyl contortion with an average displacement angle of 12.4(3) degrees, similar to [6]CPP (12.6 deg.),12 and a torsional angle (α) of 10(3) deg. which is smaller than [5]CPP and [6]CPP at 12(2)9 and 26.4(7) deg.,31 respectively. Last, biphenyl contortion produces a large and unexpected shielding of the “g” resonance as it shifts from 6.74 ppm in 2 to 4.82 ppm in 1, suggesting an anisotropic effect is induced by the biphenyl moieties despite the distance of ∼5.2 Å to the centroid of “P”.
The SE of 1′ (R = Me) was calculated considering two different homodesmotic reactions (Fig. S9) and six different DFT functionals (Table S2). From all twelve calculations, the SE of 1′ ranges from 128 to 147 kcal mol−1. Data obtained using B3LYP/6-31G(d) provides a SE for 1′ of ∼135 kcal mol−1, which allow us to compare with literature values. For instance, the SE of [6]CPP is 97 kcal mol−1 making an average strain of ∼16 kcal mol−1 per aryl ring.12 In comparison, the same metric for 1′ is ∼17 kcal mol−1 considering only the aromatic rings within the four distorted biphenyls. In stark contrast, the meta-connected species 1m′ has a minor contortion of ∼9 kcal mol−1 of SE or 1.1 kcal mol−1 per aryl ring in the connecting biphenyls (Fig. S10).
We examined the electronic properties of 1 via UV-vis absorption, fluorescence emission, and supported our findings through TD-DFT calculations. Its lowest energy absorption band is observed at λmax of 331 nm in CH2Cl2 at room temperature, which is red shifted by ∼80 nm relative to 1m (Fig. 3a). A fluorescence red shift of 81 nm is seen between 1 (λem = 425 nm) and 1m (λem = 344 nm). A visual comparison of the fluorescence shift is shown in Fig. S11. We attribute these differences to greater π conjugation in 1. The fluorescence quantum yield (ϕf) of 1 is 0.06 as determined by the relative method described by Williams et al.32 using anthracene as standard. For reference, [5]CPP and [6]CPP are not emissive.9,10,31 The smallest member in the [n]CPP family with measurable fluorescence is [7]CPP with ϕf of 0.0007.33 Based on TD-DFT calculations, we attribute the major absorption band in 1 at λmax of 331 nm to H → L + 1 and H → L + 2 transitions, where H and L stands for HOMO and LUMO, respectively. The H → L transition is forbidden (Table S3), similar to [n]CPPs.11 Despite heavy contortion of 1, its HOMO displays an overall benzenoid structure distributed over all twelve phenylenes atop the resorcin[4]arene (Fig. 3b). DFT results on 1′ indicate a HOMO–LUMO gap of 3.99 eV, while the unstrained analogue 1m′ has a larger gap of 4.73 eV, which is similar to the trend in [n]CPPs where the HOMO–LUMO gap increases as strain decreases.34
![]() | ||
Fig. 3 (a) UV-vis absorption and emission spectra of 1 (blue traces) and 1m (black traces) collected in CH2Cl2 at room temperature. (b) HOMO and LUMO density plots of 1 (±0.02 au). |
The rigidity of 1 produces a large shape-persistent cavity. The crystal structure of 1 showcases cofacial packing producing a shared cavity space for twelve MeCN molecules. Using the molecular volume of MeCN (87.8 Å3)35 and the void space calculator of Olex2,36 we determine the size of this cavity to be between 465 and 527 Å3 (Fig. S13). Molecular containers display cavities from ∼90 to ∼400 Å3,37–44 except for a recent report describing a species with a cavity of ∼800 Å3.45 Despite these reports, hosts with large pseudospherical cavities are rare.46
We hypothesized 1 can host fullerenes, similar to 1m.26 1H NMR titration experiments showed that addition of C60 or C70 to 1 generates a new set of resonances that co-exist with those of free host 1 (Fig. 4) suggesting a large association constant (Ka) for adduct formation. No further spectral changes are observed after addition of 1 equivalent of C60. The benzal proton in 1, resonance “i”, experiences an upfield shift from 5.96 to 5.87 ppm and points directly at the π surface of C60. Computational data demonstrates that shielding effects occur above the hexagonal rings in C60 and C70, and deshielding over the pentagons.47 We conclude that on average the benzal proton in C60 ⊂ 1 points at hexagons in C60. Titration of C70 to 1 also results in an upfield shift of the benzal proton to 5.87 ppm (Fig. S15); thus, nesting of C60 and C70 within 1 occurs in a similar fashion. Last, the diffusion properties of 1 and its fullerene adducts were probed via DOSY NMR. We expect 1 to have a larger diffusion coefficient (D) than its fullerene adducts (D = 1.34(1) × 10−10 m2 s−1, Fig. S16); however, it was surprising to find that adducts C60 ⊂ 1 and C70 ⊂ 1 also differ in its D from 1.26(1) × 10−10 to 1.10(1) × 10−10 m2 s−1 (Fig. S17 and S18), respectively.
Binding data for fullerene ⊂ 1 was obtained from fluorescence quenching experiments. Titration of C60/C70 into 1 were performed in TeCA (Fig. S19 and S20). Fitting of the data using Bindfit48 reveals Kas of 5.7(13) × 105 and 5.5(19) × 105 M−1 for C60 ⊂ 1 and C70 ⊂ 1, respectively, which are lower than the non-strained analogues C60 ⊂ 1m and C70 ⊂ 1m.26 These Kas are medium-to-large,49 with the exception of a few examples displaying higher binding affinities.50,51 Last, we calculated the non-covalent interaction surface through the Hirshfeld partition of molecular density (IGMH) to visualize the nesting between C60 or C70 and 1.52 The IGMH isosurface shows the interaction between the contorted biphenyls, π contacts, and the benzal protons of 1 with fullerene's π surface (Fig. S21).
In summary, we developed the synthesis of 1 – a highly strained deep cavitand built atop an alkyl resorcin[4]arene scaffold. The final product contains four highly contorted arch-shaped biphenylenes. The overall strain energy of 1 is ∼135 kcal mol−1. The rigid structure of 1 produces a large internal void space enough to accommodate fullerenes. Our approach will allow the synthesis of other highly strained molecules pushing the known boundaries in contorted aromatic chemistry.
This research was partially supported by NSF CAREER CHE-2302628 and funds by Rice University. R. H. S. acknowledges the support of the Robert A. Welch Foundation Young Investigator Award and Welch Foundation grant C-2142-20230405. We thank Dr Christopher L. Pennington for assistance with HRMS. S. M. acknowledges the support from the Dietrich School of Arts & Sciences Graduate Fellowship and the Andrew Mellon Pre-doctoral Fellowship. We thank the support from the CRC and instrumentation made available through SEA at Rice University.
CCDC 2393731 contains the supplementary crystallographic data for this paper.53
Footnote |
† These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2025 |