Open Access Article
Biying Liu†
a,
Ying Wang†
b,
Miao She†a,
Tianfu Liu
cd,
Manxin Luoa,
Yongxin Xiea,
Defei Liu
a,
Ziyi Dua,
Liyun Huang*a and
Wenbing Yuan
*a
aSchool of Environmental and Chemical Engineering, Foshan University, Foshan 528000, P. R. China. E-mail: huangliyun@fosu.edu.cn; hnyuanwb@126.com
bSchool of Environmental and Chemical Engineering, Wuyi University, Jiangmen 529020, P. R. China
cState Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, P. R. China
dUniversity of the Chinese Academy of Sciences, Beijing 100049, China
First published on 9th March 2026
The separation of benzene and cyclohexane is a long-standing challenge in the petrochemical industry, yet conventional distillation remains highly energy-intensive due to the nearly identical boiling points of the two molecules. Adsorptive separation offers a far more energy-efficient alternative; however, the development of adsorbents that simultaneously combine stable lightweight structures, high adsorption capacity, and strong selectivity remains difficult. To overcome these limitations, we present a robust hydrogen-bonded organic framework, H8ETTOB-HOF, constructed from a polyphenyl octacarboxylic acid linker. Its precisely engineered pore architecture and tailored electrostatic environment work synergistically to enhance benzene affinity, affording exceptional benzene/cyclohexane selectivity (16.9) and a high adsorption capacity (3.9 mmol g−1 at 298 K, 13 kPa), surpassing all previously reported HOFs. Density functional theory shows that the synergistic C–H⋯O and C–H⋯π interactions facilitate benzene capture. Experimental investigations employing dynamic columns demonstrate efficient one-step benzene purification from equimolar mixtures, showcasing excellent recyclability. This work highlights a strategic approach to constructing stable, multifunctional HOFs for challenging hydrocarbon separations.
Against this backdrop, both metal–organic frameworks (MOFs) and hydrogen-bonded organic frameworks (HOFs) have emerged as highly promising classes of materials. MOFs often achieve ultra-selective recognition of targeted molecules through precise pore engineering and functionalization,6–8 such as JNU-81-C,9 UiO-66/MFM-300 (ref. 10) and Zn-Ade-TCPE.11 In contrast, HOFs possess distinct advantages due to their mild synthesis conditions, readily tunable pore structures, and inherent molecular recognition capabilities. HOFs are porous crystalline materials composed of molecular building blocks connected through intermolecular hydrogen bonds.12–15 HOFs possess significant advantages such as light weight, high crystallinity, processability in solution, exceptional recyclability, and easy purification, endowing them with promising potential for gas adsorption and separation, proton conduction, drug delivery, and so on.16,17 Recent advancements in HOFs have demonstrated their potential in the selective adsorption of aromatic compounds like benzene. For instance, HOF-8 displays a benzene uptake of ∼2.5 mmol g−1 at 298 K, with its ordered π-stacked pores favoring aromatic adsorption via van der Waals and π⋯π interactions.18 The B ← N-based framework of BNH-1 shows ∼2.2 mmol g−1 toluene uptake and effectively separates it from methylcyclohexane, enabled by directional C–H⋯π interactions and a self-healing hydrogen-bonded network.19 Nevertheless, achieving efficient discrimination between structurally similar aromatic and aliphatic molecules—such as benzene and cyclohexane—remains challenging for HOFs, due to limitations in pore rigidity, binding site definition, and framework persistence. These challenges underscore the need for rational design of HOFs featuring enhanced connectivity, shape-selective cavities, and stable hydrogen-bonded architectures to achieve high selectivity and structural robustness for aromatic-aliphatic separations.20,21
For HOFs, the well-known inherent weakness of hydrogen bonds often renders them fragile, making it challenging to sustain their permanent porosity, which significantly impedes the use of these materials. To date, considerable endeavors have been undertaken to construct stable HOFs with accessible permanent porosity.22 Among them, increasing the number of connection sites in the building blocks improves connectivity and promotes interpenetration. This structural reinforcement contributes to greater framework robustness. Furthermore, strengthening the interlocking interactions between interpenetrating layers can further improve the chemical and thermal stability of HOFs. Moreover, the combination of rigid and flexible organic polycarboxylic acids can increase the connectivity of the constructed HOF, and it is also likely to increase the degree of interpenetration, which is a very effective way to improve the stability of HOFs. Currently, the vast majority of ligands in HOFs are at most six-linked. Bu et al. reported rare 8-connected carboxyl hydrogen-bonded organic frameworks, NKM–HOF–1 and NKM–HOF–2, for white light emission.23 Generally, HOFs are “soft” and always in a balanced state of robustness and flexibility.24 HOFs constructed from rigid and flexible organic compounds with more than six potential multi-linker sites happen to have this softness or flexibility, which increases the interpenetration degree, and their inherent pores have a more sensitive corresponding adaptability to the storage of guest molecules.
Thus, designing a new HOF should focus on organic polyphenyl linkers, multi-connectivity, and a rigid framework to maintain narrow windows, allowing for effective host-guest van der Waals interactions, especially with benzene over cyclohexane. Herein, we designed and synthesized a novel tetraphenylethylene (TPE)-based octacarboxylic HOF, designated as H8ETTOB-HOF, constructed from 4,4′,4″,4‴,4′′′′,4′′′′′′,4′′′′′′,4′′′′′′–((ethene-1,2,2 tetrayltetrakis (benzene-4,1-diyl))tetrakis- (azanetriyl))octabenzoic acid (H8ETTOB, Fig. 1a) as the molecular building block. This stable framework with π-conjugated aromatic linkers or aromatic pore surfaces can preferentially interact with aromatic guests compared to aliphatic ones. The judicious arrangement of these components facilitates stable assembly, exhibiting potential for selective adsorption and separation of Bz from Cy.
776.3 (5) Å3, and a density of 1011 g cm−3 at 100 K (Table S2). The H8ETTOB framework is stabilized by eight carboxylic groups via O–H⋯O hydrogen bonds, forming a robust three-dimensional (3D) network.26 The O–H⋯O bond lengths range from 2.586 to 2.653 Å with angles near 170°, indicating strong and directional hydrogen bonding. These values are consistent with reported data for similar HOFs (Fig. 1 and S3, Table S3).27,28
In this structure, adjacent layers are interconnected by additional O–H⋯O hydrogen bonds (2.601 Å – 2.699 Å), reinforcing the stability of the 3D framework (Fig. 1b, c and S3). The C37–H37⋯O12 and C56–H56⋯O4 interactions (3.455 Å and 3.400 Å, respectively, Table S3) further enhance the structural integrity.29,30 These supramolecular interactions stabilize the framework while imparting flexibility due to their weak nature. Fig. 1d, e and S4 illustrate the hydrogen bonding interactions within and between layers, forming a robust 3D interlocked network. This structural arrangement, featuring intra- and inter-layer hydrogen bonds, generates stable yet flexible porosity, making H8ETTOB-HOF promising for gas adsorption and separation.31
Notably, H8ETTOB acts as a building block for H8ETTOB-HOF, which exhibits an unexpected coordination number of six, forming unconventional “double link” connections (Fig. 1d, e and S4). For accurate topological analysis, we deconstructed each 8-c building block into one four-connected (4-c) vertex and four trigonal (3-c) vertices, revealing a 6-c topology (Fig. 1d).32 The layers feature alternating large and small hexagonal structures. The slight misalignment between adjacent large hexagons creates zigzag channels, with a similar pattern observed in the smaller hexagons. This results in two distinct one-dimensional (1D) channels with window sizes of approximately 12 × 8 Å2 and 9 × 6 Å2 (Fig. S5), sufficiently large to accommodate Bz and Cy molecules (Table S1).
The network topology of H8ETTOB-HOF, as analyzed by the ToposPro program, features fully interdigitated symmetry elements, classifying it as a twofold interpenetrated framework.27,33 The framework exhibits a 6-connected net with the Schläfli symbol 412·63. As shown in Fig. S5, H8ETTOB-HOF displays a pcu (primitive cubic) topology with alpha-Po-like charateristics, integrating the geometric features of both nets, which likely endows the material with distinctive physicochemical properties. The H8ETTOB building block exhibits flexible and adaptable structural characteristics facilitated by bond rotation, the presence of a helical axis (Fig. 1a), and variations in bond angles, to form a stable and functional framework. The single-layer structure of H8ETTOB-HOF is connected by hydrogen bonds, displaying a “double-spiral upward figure-eight” pattern along the b-axis (Fig. S6). Correspondingly, two layers of these “spiral eight-shaped chains” of H8ETTOB-HOF become intertwined and risen. These robust interlayer hydrogen bonding interactions and the structural feature of two helices, combined with the well-defined porous channels, contribute to the stability and functionality of the 3D framework of H8ETTOB-HOF.
The powder X-ray diffraction (PXRD) patterns of H8ETTOB-HOF were carefully examined under various conditions, including as-synthesized and thermally activated (473 K) (Fig. 2a). The experimental patterns match well with the simulated one, confirming framework integrity upon thermal treatment. To evaluate chemical stability, H8ETTOB-HOF was immersed in concentrated hydrochloric acid (12 M), NaOH aqueous solution (pH = 11), H2O, methanol, ethyl acetate, and acetone for 3 days, and PXRD showed that H8ETTOB-HOF maintained its crystallinity (Fig. 2a and S7c). Optical microscopy and scanning electron microscopy (SEM) results show that the crystal form does not change under various soaking conditions (Fig. S1 and S2). These observations demonstrate that H8ETTOB-HOF maintains its morphological integrity under chemically aggressive conditions. Furthermore, thermal stability was investigated using variable-temperature X-ray diffraction (VT-XRD): the diffraction patterns remain consistent from 298 K to 473 K, confirming that the crystal structure of H8ETTOB-HOF is stable up to high temperatures (Fig. S7d). Complementary FTIR spectra (Fig. S7a) reinforce the stability of the framework, with strong absorption bands in the –COOH region (1700–1500 cm−1) affirming the integrity of carboxylic groups within the HOF. Notably, the TGA curves (Fig. S7b) for both untreated and post-activation H8ETTOB-HOF reveal comparable weight loss patterns, with significant decomposition commencing around 573 K, attributed to the breakdown of organic ligands. Compared to similar hexacarboxylic acid-based HOFs reported in the literature,34 H8ETTOB-HOF exhibits superior thermal stability, likely due to its high connectivity and robust framework.25 Additionally, nitrogen adsorption–desorption isotherms (Fig. 2b) of H8ETTOB-HOF activated at 473 K reveal a Brunauer–Emmett–Teller (BET) surface area of 555.11 m2 g−1 and a pore volume of 0.26 cm3 g−1. Nonlocal density functional theory (NLDFT) pore size distribution indicated a main pore at 6.4 Å, with additional pores distributed in the range of 12.0 to 18.5 Å (Fig. 2b). Notably, the H8ETTOB-HOF exhibits remarkable chemical stability, as evidenced by its BET surface areas (Fig. 2c) of 509.18 and 453.73 m2 g−1 after immersion in concentrated hydrochloric acid (12 M) and water for 3 days, respectively. These results are consistent with the PXRD data, further confirming the structural integrity of H8ETTOB-HOF under harsh conditions.
The demonstration of permanent porosity in H8ETTOB-HOF promoted an investigation into the single-component adsorption isotherms of Bz and Cy at 298, 318, and 373 K, at pressures up to 13 kPa. As shown in Fig. 3a, S6a and S7, H8ETTOB-HOF exhibits a pronounced preference for Bz over Cy at all tested temperatures. At 298 K and 13 kPa, the Bz uptake reaches 3.9 mmol g−1, significantly higher than that of Cy (1.2 mmol g−1), resulting in a high Bz/Cy uptake ratio. Ideal adsorbed solution theory (IAST) was applied to calculate the adsorption selectivity of H8ETTOB-HOF towards binary Bz/Cy mixtures (50
:
50, v/v) at 298, 318, and 373 K. As shown in Fig. 3b and S9, the IAST selectivity of H8ETTOB-HOF for benzene over cyclohexane exhibits a notable pressure-dependent increase at various temperatures. At 298 K, the selectivity reached 16.9 at 14 kPa, underscoring its pronounced affinity for benzene under ambient conditions. As the temperature increased, the selectivity decreased to approximately 10 at 318 K, indicating that elevated thermal energy reduced the differential adsorption affinity between benzene and cyclohexane. Subsequently, the selectivity rose significantly to 20 at 373 K. This recovery is attributed to a pronounced thermal contraction of the HOF pore apertures, which imposes a stronger steric restriction on the bulkier cyclohexane molecule.35 In contrast, the planar benzene molecule maintains efficient adsorption through resilient π–π interactions with the aromatic pore walls, thereby preserving high separation efficacy even at elevated temperature interactions.36 The nearly linear progression of selectivity with pressure at each temperature suggests a stable adsorption behavior, maintaining significant Bz/Cy discrimination even at elevated pressures. The PXRD patterns confirm the structural stability of H8ETTOB-HOF after Bz/Cy uptake (Fig. S8a).
![]() | ||
| Fig. 3 (a) Adsorption isotherms of Bz and Cy on H8ETTOB-HOF at 298 K, 318 K and 373 K; (b) IAST selectivity of H8ETTOB-HOF for 50/50 Bz/Cy at different temperatures. (c) Isosteric heat of adsorption (Qst) plots for the adsorption of Bz and Cy by H8ETTOB-HOF. (d) Adsorption kinetic curves of Bz and Cy on H8ETTOB-HOF at 298 K under 200 Pa. (e) Trends in performance for the HOFs.45,46 (f) Performance comparison of IAST selectivity between H8ETTOB-HOF and classical adsorbents.47–52 | ||
The adsorption heat (Qst) for Bz was calculated to be approximately 42 kJ mol−1 at zero coverage, higher than that for Cy (∼30 kJ mol−1), indicating stronger interactions with Bz (Fig. 3c and S11).37,38 Additionally, the gentle upward trend in the Qst curve over a broad range of adsorbed amounts suggests additional intermolecular interactions as gas molecules progressively stack within the pore channels until saturation.39 The exceptional performance of H8ETTOB-HOF in Bz capture is attributed to its fine-tuned pore structure, the polyphenyl environment of its pore walls, and its strong binding affinity for Bz molecules. The adsorption isotherms of benzene in H8ETTOB-HOF display a distinct profile, indicating strong interactions and efficient accommodation within the framework's channels (Fig. S10). Notably, the adsorption data were effectively modelled using the dual-site Langmuir–Freundlich (DSL-F) fitting, which revealed the presence of two types of adsorption sites within the H8ETTOB-HOF framework. The high-affinity sites, as indicated by the Langmuir model, suggest strong and uniform adsorption, likely driven by π–π interactions and hydrogen bonding with benzene.40 In contrast, the low-affinity sites, modelled by the Freundlich approach, are indicative of a more heterogeneous surface, supporting multilayer adsorption and weaker van der Waals interactions, which are more characteristic of cyclohexane.41 This dual-site nature elucidates the observed higher selectivity and uptake of benzene, particularly at lower temperatures and higher pressures.
To further elucidate the adsorption kinetics of benzene Bz and Cy within H8ETTOB-HOF, we measured the kinetic adsorption curves at 298 K and 200 Pa. As depicted in Fig. 3d, Bz reached adsorption equilibrium in approximately 160 min, whereas Cy required around 190 min to do so. The diffusion time constant (D/r2), obtained from the short-time approximation of the diffusion equation,42 reveals a kinetic selectivity of 21.4 for Bz/Cy in the H8ETTOB-HOF (Fig. S12 and Table S5–S6). The significantly higher value implies a greater relative diffusivity of Bz compared to Cy, likely attributable to the size-exclusion effect imposed by the pore channels on the bulkier Cy molecules.43 The faster adsorption rate of Bz compared to Cy suggests a more rapid intra-crystalline diffusion for Bz, indicating a higher affinity and interaction with the framework of H8ETTOB-HOF.44 Fig. 3e presents a radar chart comparison of key performance indicators—including IAST selectivity, Bz uptake, BET surface area, uptake ratio, density, and thermal stability—among H8ETTOB-HOF and representative HOF analogs such as X–HOF–4a45 and HOF-Tri-PMDI-Br.46 H8ETTOB-HOF demonstrates a balanced profile, with notably higher Bz uptake and IAST selectivity, indicating its superior adsorption capacity and affinity toward aromatic molecules under ambient conditions. Fig. 3f presents a comparative analysis of IAST selectivity for Bz over Cy across various porous materials. H8ETTOB-HOF exhibits a notably higher selectivity (16.9) than most classical COFs and MOFs. This highlights its superior molecular recognition capability among metal-free frameworks.
To verify the selective adsorption capabilities of H8ETTOB-HOF in mixed Bz/Cy, we conducted column breakthrough experiments at 298 K using a 50
:
50 Bz/Cy mixture. At a flow rate of 0.6 mL min−1, the breakthrough curve shows that Cy was eluted from the column first at 59.4 min·g−1, while Bz began to penetrate the column at 112.8 min·g−1, with a time interval of 53.4 min·g−1. This indicates that H8ETTOB-HOF has a higher selectivity for Bz than for Cy (Fig. 4a). At a flow rate of 1.2 mL min−1 for the Bz/Cy mixture, the separation factor of 3.35 (eqn (S6)) outperforms HOF-FJU-112a (2.28),36 HOF-202 (2.79)53 and HOF-303 (2.95).54 Moreover, as the flow rate increased from 0.6 to 0.9 and 1.2 mL min−1, the retention time of Bz on the H8ETTOB-HOF column gradually decreased from 53.4 to 43.6 and 28.9 min·g−1. The adsorption capacity of the material remained stable at approximately 3.1–3.6 mmol g−1 during three consecutive regeneration cycles under a flow rate of 1.2 mL min−1 (Fig. 4b and c), which further confirmed the stability of H8ETTOB-HOF. Besides, the samples after three penetration cycles still maintained consistency with the PXRD, FT-IR and 77 K N2 adsorption curves of the pristine samples, indicating that the samples were stable (Fig. S7a and S8).
To gain a deeper understanding of the selective adsorption behavior in H8ETTOB-HOF, density functional theory (DFT) calculations were employed to analyze the energy profiles and molecular interactions of Bz and Cy. The electrostatic potential (ESP) analysis of H8ETTOB-HOF, Bz, and Cy molecules reveals that the positive ESP on Bz's H atoms complements the negative ESP along H8ETTOB-HOF's pore channels (Fig. 5a and b).55 This electrostatic synergy drives the selective adsorption of Bz. Fig. 5c, d and S13 further show that Bz forms C–H⋯π interactions (2.99 Å) and C–H⋯O interactions (2.21 to 3.12 Å) with H8ETTOB-HOF, whereas Cy exhibits only C–H⋯O interactions (2.61 to 3.08 Å). The binding energy of Bz (−79 kJ mol−1) is higher than that of Cy (−62 kJ mol−1).56
The energy variations as Bz and Cy traverse different adsorption sites—transitioning from the initial state (IS) through the transition state (TS) to the final state (FS)—are depicted in Fig. 5e and S14. The lower TS energy of Bz (12.7 vs. 21.2 kJ mol−1) indicates its thermodynamic preference in H8ETTOB-HOF. Further analysis of the energy landscapes for Bz and Cy adsorption across different sites-outside, edge, aperture, and cavity—is presented in Fig. 5f. Bz exhibits relative energies of −37.6 kJ mol−1 (edge), −19.4 kJ mol−1 (aperture), and −8.7 kJ mol−1 (cavity), wheares Cy shows weaker relative energies: −15.8 kJ mol−1 (edge), 45.5 kJ mol−1 (aperture), and −3.6 kJ mol−1 (cavity). This stark contrast underscores Bz's superior affinity.57,58 The energy differences are attributed to the structural properties of Bz and Cy: Bz's planar and aromatic structure facilitates C–H⋯π and C–H⋯O interactions with the framework, enhancing its adsorption stability.59,60
Footnote |
| † These authors equally contributed to this work. |
| This journal is © The Royal Society of Chemistry 2026 |