Natalia S.
Soldatova
a,
Pavel S.
Postnikov
*ab,
Daniil M.
Ivanov
ac,
Oleg V.
Semyonov
a,
Olga S.
Kukurina
a,
Olga
Guselnikova
ad,
Yusuke
Yamauchi
de,
Thomas
Wirth
f,
Viktor V.
Zhdankin
g,
Mekhman S.
Yusubov
a,
Rosa M.
Gomila
h,
Antonio
Frontera
i,
Giuseppe
Resnati
*aj and
Vadim Yu.
Kukushkin
*c
aResearch School of Chemistry and Applied Biomedical Sciences, Tomsk Polytechnic University, Tomsk 634034, Russian Federation. E-mail: postnikov@tpu.ru
bDepartment of Solid State Engineering, Institute of Chemical Technology, Prague 16628, Czech Republic
cInstitute of Chemistry, Saint Petersburg State University, Saint Petersburg 199034, Russian Federation. E-mail: v.kukushkin@spbu.ru
dJST-ERATO Yamauchi Materials Space-Tectonics Project, National Institute for Materials Science (NIMS), 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
eAustralian Institute for Bioengineering and Nanotechnology (AIBN), The University of Queensland, Brisbane, QLD 4072, Australia
fSchool of Chemistry, Cardiff University, Park Place, Cardiff, UK
gDepartment of Chemistry and Biochemistry, University of Minnesota, Duluth, MN 55812, USA
hServeis Científico-Tècnics, Universitat de les Illes Balears, Crta. de Valldemossa Km 7.5, 07122 Palma de Mallorca, Spain
iDepartament de Química, Universitat de les Illes Balears, Crta. de Valldemossa Km 7.5, 07122 Palma de Mallorca, Spain
jNFMLab, Department of Chemistry, Materials and Chemical Engineering “Giulio Natta”; Politecnico di Milano, via Mancinelli 7, I-20131 Milano, Italy. E-mail: giuseppe.resnati@polimi.it
First published on 12th April 2022
Porous architectures characterized by parallel channels arranged in honeycomb or rectangular patterns are identified in two polymorphic crystals of a zwitterionic 4-(aryliodonio)-benzenesulfonate. The channels are filled with disordered water molecules which can be reversibly removed on heating. Consistent with the remarkable strength and directionality of the halogen bonds (XBs) driving the crystal packing formation, the porous structure is stable and fully preserved on almost quantitative removal and readsorption of water. The porous systems described here are the first reported cases of one-component 3D organic frameworks whose assembly is driven by XB only (XOFs). These systems are a proof of concept for the ability of zwitterionic aryliodonium tectons in affording robust one-component 3D XOFs. The high directionality and strength of the XBs formed by these zwitterions and the geometrical constraints resulting from the tendency of their hypervalent iodine atoms to act as bidentate XB donors might be key factors in determining this ability.
Very recently, the modern palette of porous materials has been enriched by novel systems wherein the structural units are assembled via relatively weak interactions, with supramolecular organic frameworks18–20 (SOFs) and hydrogen-bonded organic frameworks21–26 (HOFs) being the most common representatives. However, the development of these materials is challenging as a consequence of the lower directionality of weak noncovalent bonds in comparison to coordination or covalent bonds affording MOFs and COFs.27,28 In this context, halogen bonding (XB) is becoming a novel tool for the controlled assembly of molecular tectons into highly ordered frameworks29,30 and is here employed for the obtainment of 3D and one-component purely halogen-bonded organic frameworks (XOFs). It was considered that the high directionality of the interaction might act as an advantage of the approach. Notably, a one-component 2D XOF was observed in crystals of tetrabromobenzene-1,4-dicarboxylic acid.30
Among the wealth of porous molecular crystals wherein XB contributes to sustaining the architecture, only two examples that can be considered as 3D XOFs were obtained via (i) cocrystallization of the tetrahedral pyridine-based XB acceptor and 1,3,4,5-tetraiododifluorobenzene to form a di-component XOF29 and (ii) assembly of a heterotopic tecton to give an unstable 3D XOF.31
The well-developed strategy of assembly of XB donors with nitrogen-based XB acceptors furnishes,32,33 in particular, halogen-bonded 3D structures (exhibiting diamondoid-like cage entities), which are promising porous materials; these structures have not been tested as adsorbents.34–38 In some other cases, assembly of iodine species with N-based XB acceptors provided 2D or 1D architectures, featuring cavities.29,39–41
Other reported examples of XB-assembled SOFs are built by the simultaneous action of XB and hydrogen bonding (HB).42,43
Diaryliodonium salts44,45 can function as useful building blocks in supramolecular chemistry46 due to their ability to form discrete adducts and 1D chains or ribbons.47–50 The iodine typically functions as a bidentate XB donor51–56 and the XB acceptor is, by far in most cases, the anion. Despite the wide applicability of iodonium salts as halogen bond donors in catalysis,54–57 the formation of 3D structures with well-defined pores is still unknown. Among other iodonium salts, the zwitterionic compounds can be promising for the directed design of 3D networks due to the presence of two binding sites. Nevertheless, the number of zwitterionic systems characterized via single-crystal X-ray diffractometry (XRD) is too small to enable statistically meaningful statements, but it is interesting to observe that although 1D chains and ribbons58–62 are the most common halogen-bonded systems formed by these iodonium salts, the tendency to form multidimensional nets, e.g., 2D59,63 networks, seems definitely higher than when starting from diaryliodonium salts wherein the positive and negative sites are in different molecular entities.
We decided to verify whether zwitterionic iodonium species can indeed function as promising tectons for the reliable design and synthesis of multidimensional supramolecular architectures in general and new XOFs in particular. In this study, we synthesized novel zwitterionic 4-(aryliodonio)benzene-sulfonate species (Scheme 1) and we report that crystallization of one of them leads to porous XOFs. We reasoned that the strong electrostatic attraction between the cationic and anionic sites in the zwitterion may act synergistically with the remarkable directionality of the XB and enable the self-assembly of robust XOFs. The design and preparation of XOFs based on zwitterionic iodonium salts may thus represent a new opportunity in the challenging and rewarding field of novel porous organic materials. Density functional theory (DFT) calculations provide an insight into the relative intensity of both σ-holes located at the extensions of both C–I bonds and the energetic features associated with the XBs.
Scheme 1 Synthesis of zwitterionic iodonium salts 1 (Ar = 2,4,6-Me3C6H2; 74%), 2 (2,5-Me2C6H3; 60%), and 3 (4-ClC6H4; 98%). |
The overall crystal packings of iodonium derivatives 1–3 are quite different from each other, but strict similarities exist in the attractive interactions driving zwitterion binding. Some C–H⋯π and C–H⋯O HBs are present in crystals 1–3 but the overall crystal packing of these systems is controlled by the short, directional, and charge assisted I⋯O XBs (Table 1). Iodine atoms function as bidentate XB donors and sulfonate moieties as bidentate acceptors, with two out of the three sulfonate O atoms acting as monodentate donors of electron density. As typical for XBs, oxygen atoms approach iodine atoms approximately on the extension of C–I covalent bonds. C–I⋯O angles span the range 167–175° but in two interactions involving the disordered sulfonated group in 1b where greater deviations from linearity are observed.
Compound | C1–I⋯O | d(I⋯O), Å | ∠C1–I⋯O, ° | C7–I⋯O | d(I⋯O), Å | ∠C7–I⋯O, ° |
---|---|---|---|---|---|---|
a Probably due to rotation of one sulfonate group of chain 1 around the C–S covalent bond, the oxygen atoms of the corresponding sulfonate group are disordered and have been modeled over two positions; geometric parameters of the XBs associated with the second position are C1–I1⋯O3, d(I⋯O) 2.988(3) Å, ∠C1–I1⋯O3 171.39(7)°; C7B–I1B⋯O1, d(I⋯O) 3.231(3) Å, C7B–I1B⋯O1 149.22(9)°. | ||||||
1a | C1–I1⋯O3 | 2.861(4) | 168.39(10) | C7–I1⋯O2 | 2.877(2) | 170.31(10) |
1b chain 1 | C1–I1⋯O3D | 3.267(3) | 166.70(8) | C7–I1⋯O1B | 2.8305(16) | 171.08(7) |
C1B–I1B⋯O3B | 3.0577(18) | 169.14(7) | C7B–I1B⋯O1D | 2.755(3) | 155.97(10) | |
1b chain 2 | C1A–I1A⋯O2A | 2.949(2) | 173.56(9) | C7A–I1A⋯O1A | 2.8173(19) | 169.60(6) |
2 | C1–I1⋯O2 | 3.007(5) | 171.34(17) | C7–I1⋯O1 | 2.864(4) | 169.41(12) |
3 chain 1 | C1–I1⋯O3 | 2.900(4) | 175.07(16) | C7–I1⋯O1 | 2.734(4) | 171.94(16) |
3 chain 2 | C1A–I1A⋯O3A | 2.910(4) | 171.90(16) | C7A–I1A⋯O1A | 2.763(4) | 173.24(15) |
The I⋯O distances, typically varying in the range 2.73–3.05 Å, are similar to those in crystalline 2-(aryliodonio)benzene-sulfonates60 and are substantially shorter than the sum of Bondi vdW radii69 (ΣvdW(O + I) is 3.5 Å). Deviations from these most common values are observed once again for interactions involving the disordered sulfonate group in 1b, the I⋯O1 and I⋯O3D distances increasing to 3.231 and 3.267 Å. Interestingly, in 1b, 2, and 3, the C–I⋯O XBs opposite the phenylsulfonyl moiety are 0.1–0.2 Å longer than the C–I⋯O XBs opposite the other aryl group at iodine, despite the σ-holes opposite the former moiety being more positive than those opposite the second one. A possible rationalization is suggested by the modelling of the interactional landscape present in the packings and is discussed later on. In 1a, the difference between the two C–I⋯O separations is negligible and within 3σ.
A halogen-bonded dimer formed on antiparallel pairing of two zwitterions is a common motif present in all crystals 1–3. A similar motif is present also in the zwitterion 4-(CF3SO2NSO2)C6H4IPh.59 These dimers are further connected to each other via XBs and linear ribbons (1D-network) are present in 1b, 2, and 3 (Fig. 1c and S1†). To the best of our knowledge, a similar 1D-net has never been reported for any other zwitterionic iodonium salts.59–61 One ribbon is found in the lattice of 2, while two distinct but very similar ribbons are present in 1b and 3. In 1a, any dimer formed on antiparallel pairing of two zwitterions is halogen-bonded to four other analogous dimers (Fig. 1a) and a 3D network is formed.
Consistent with the fact that the density of 2 and 3 (1.942 and 2.061 g cm−3, respectively) is much greater than that of 1a and 1b (1.492 and 1.399 g cm−3, respectively), the crystal packings of the former two systems present no voids while the packings of the latter two systems are porous and show parallel channels, along the c and a axes, respectively (Fig. 2). Water is present in these channels but could not be identified through crystallography due to its mobility.
Channels in both crystals have a cylindrical shape with a diameter of ∼10 Å in 1a and ∼12 Å in 1b; the percentage of empty volume is 22.5 in 1a and 25.0 in 1b. Before SQUEEZE in Platon70 was applied, we identified numerous residual peaks, ranging from 1.01 to 2.70 e, in the channel structure of 1a. We believe that the observed residual peaks are related to adsorbed disordered water that was also identified by FTIR and NMR spectroscopy; the NMR spectra of 1a do not exhibit any signals of MeOH. Additional evidence favoring the presence of water was obtained from the IR spectra, where the water peaks were found at 3420 cm−1 (Fig. 3e). Moreover, NMR experiments (residual DMSO in DMSO-d6 was used as an internal standard) clearly demonstrated a lower water content in a dried sample in comparison with the pristine sample and a sample after water adsorption (Fig. S3, ESI†).
Parallel 1D-channels filled with solvent molecules and similar to those found in 1a and 1b have been observed in systems formed on assembly of pyridine based XB acceptors with 1,4-diiodoperfluorobenzene and 1,3,5-triiodoperfluorobenzene.24,52 In 1a, the overall arrangement of the 3D-net adopts a honeycomb-like structure (Fig. 1b) and channels display a hexagonal pattern, while in 1b the 1D-nets assemble in a rectangular geometry (Fig. 1d) which translates into a channel array with analogous topology. In contrast to 1a, the XRD pattern of 1b does not agree with the simulated pattern because of the decomposition occurring on drying (Fig. S2†).
Key features for the further harnessing of XOF 1a are its thermal stability and the possibility to remove the solvent while preserving the original crystal structure. The TG/DSC analysis revealed that crystals were stable up to 180 °C, then exothermic decomposition occurred (Fig. 3d). According to the TG data, the prepared XOF 1a can be activated at least at 140 °C without either destruction or phase transition. The initial weight-loss is tentatively associated with the removal of solvent from the voids. To assess the crystal structure preservation on solvent removal, the crystalline powder of 1a was dried for 1.5 hours until the weight of the sample remained constant. The powder XRD study of this batch (Fig. 3b) showed a pattern very similar to that of the pristine material, indicating the original crystal structure.
To verify the functional potential of XOF 1a, we tested the reversibility of water adsorption/desorption from the vapor phase. The weighed sample of dried 1a was placed in a desiccator saturated with water vapor and the weight was monitored. The amount of chemisorbed water was 5.69 mass% (calculated with respect to the mass of the pristine sample of 1a) and this value agrees well with the mass-loss from the pristine powder (5.85 mass%) and mass-loss TG data (Fig. 3d).
The consequence of water adsorption and desorption by crystals of 1a is fully reversible (Fig. 3a). Thus, we conducted a gravimetric determination of the adsorption capacity, along with XRD and spectroscopic monitoring, and found that the structure of 1a is fully preserved after at least 5 cycles of adsorption followed by desorption performed on drying. These experiments indicate that XB stabilizes the XOF crystal structure that is preserved even under relatively harsh conditions of the repeated adsorption and desorption cycles.
We also performed an experiment with the adsorption of D2O (instead of H2O) by the dehydrated sample. The IR-spectra (Fig. 3e) revealed the appearance of νas,s(D2O) stretching at approximately 2530 cm−1 and these data additionally proved the adsorption of D2O from the gas phase by XOF 1a (the assignment of the IR bands is provided in ESI Table S2†).
To get additional data illustrating the stability of XOF 1a, we also carried out variable-temperature Raman spectroscopy experiments. We observed that these spectra remain unchanged even after heating XOF 1a at 140 °C for 1 h (Fig. 3c). The Raman spectrum of 1a, after 5 cycles of water adsorption–desorption, demonstrated the preservation of the crystal structure.
The porous structure of 1a was also studied using nitrogen adsorption–desorption isotherm measurements. Fig. 4a shows a gradual adsorption of nitrogen gas with a small hysteresis loop. The BET surface area of the XOF exhibits an expectedly low value of 17.7 m2 g−1, probably because of the presence of water in the pores. In view of the known difficulties of the Barrett–Joyner–Halenda (BJH) method,71 the pore size distribution curve was obtained by the more accurate Non-Local Density Functional Theory (NLDFT) approach (Fig. 4b). We found that the microporous structure is indeed formed and its total pore volume is 0.017 cm3 g−1. Fig. 4b shows largest pore size of about 7 Å, which corresponds to a cylindrical shaped pore with a diameter of approx. 10 Å; the latter value matches well with our X-ray crystallography data (Fig. 2). Other additional pores (approx. 1.7 nm) were also observed and they correspond to inter-crystallite spaces.
Fig. 4 (a) Nitrogen adsorption–desorption isotherm of the XOF; (b) pore size distribution curve obtained by the NLDFT approach (the enlarged image in the range from 3 to 6 Å is shown as an inset). |
In the context of the discovery of the stable 3D XOF, it is noteworthy that Uriel and coworkers31 have also employed a heterotopic tecton (featuring two σ-hole donating iodine centers, which are perpendicular to two XB accepting N sites in the same neutral molecule) for the design of XB-based XOFs. The assembly of this tecton led to porous XOFs with empty volumes of 20–24%. The obtained architectures were quite unstable and they irreversibly collapsed even on keeping in air at RT. This instability is believed to be associated31 with the weakness of the XBs involved in the construction of these XOFs. In our study, we generated XOFs of a higher stability and, moreover, we demonstrated the reversibility of water adsorption–desorption associated with these materials. Most likely, the stability is determined by a stronger XB formed between ionic sites of our iodonium zwitterions.
We have also analyzed the energetic features of some supramolecular motifs observed in the prepared systems. Crystals 1b, 2, and 3 present very similar tetrameric and dimeric motifs which are presented in Fig. S4, 6 and S5,† respectively. The tetrameric arrays are assembled via XBs, HBs, and π–π interactions and their formation energies are very large, ranging from ∼−156 to ∼−184 kcal mol−1 (Table 2). The largest value is for compound 3 in agreement with its more positive σ-holes evidenced by the MEP surface computation. Two distinct dimeric motifs can be extracted from the tetramers present in 1b, 2, and 3 and they are denoted here as dimer A and dimer B. In all three crystals, dimers A are assembled via two XBs and one π-stacking interaction (Fig. S4b, 6b, and S5b†). This π-stacking is expected to be highly stabilizing since the iodonium benzenesulfonate moiety is very polarized and probably its antiparallel pairing is highly exoenergetic. Indeed, the dimerization energies computed for dimers A range from ∼−61 to ∼−75 kcal mol−1 (Table 2). The combined QTAIM/NCIplot analysis reveals that the two equivalent XBs assembling dimers A (denoted as XB1) are characterized by a bond critical point (CP, red sphere of bond CP#1) and a bond path interconnecting the I and O atoms. Moreover, a blue NCIplot isosurface is located between the O and I atoms, thus showing a strong interaction. The π–π stacking is characterized by an extended and green isosurface that embraces both aromatic systems. The formation of dimers B is sustained by one single XB (denoted as XB2 and resulting in bond CP#2; Fig. S4c, 6c, and S5c†) and some ancillary contacts like C–H⋯O HBs (resulting in bond CP#3) and O⋯π interactions. It is thus no surprise that the corresponding dimerization energies are significantly smaller (in absolute value) than those of dimer A (e.g., −26.3 vs. −66.9 kcal mol−1 for compound 2, Table 2). Each interaction is characterized also by the corresponding bond path and NCIplot isosurface. The blue color of the NCIplot isosurface that characterizes XB2 is less intense than that of XB1, thus suggesting that XB1 is stronger, in line with the shorter XB experimental distance (Table 1).
Compound | ΔEtetramer | ΔEdimer A | ΔEdimer B |
---|---|---|---|
1b | −155.7 | −61.2 | −25.8 |
2 | −167.5 | −66.9 | −26.3 |
3 | −183.9 | −74.7 | −28.2 |
As described above, the solid-state architecture of compound 1a is different from that of 1b, 2, and 3 and the tetrameric motif observed in these three latter crystals is not present in 1a, while the two dimeric motifs A and B are present with minor differences in the interaction pattern with respect to 1b, 2, and 3 (Fig. 7). QTAIM/NCIplot analysis reveals that dimer A in 1a is assembled via two XBs (denoted as XB1) and one π-stacking, while dimer B is assembled via one single XB (denoted as XB2) and a weak ancillary C–H⋯O HB between one methyl group and one O atom of the sulfonate group. Consistent with these interaction patterns, the stabilization energy of dimer A (−62.3 kcal mol−1) is much greater than that of dimer B (−28.1 kcal mol−1).
Finally, taking into consideration that the σ-hole opposite the benzenesulfonate group is more positive than the other one, the fact that XB1 is shorter and stronger than XB2 can probably be explained considering the favorable interplay between the π-stacking and the XBs in dimers A. That is, the co-existing π-stacking and XB interactions reinforce each other.
Specifically, the single crystal X-ray analyses of 1a and 1b, the two polymorphic crystals formed by one of the zwitterionic 4-(aryliodonium)benzenesulfonates described here, reveal the presence of porous architectures wherein parallel channels are arranged in honeycomb and rectangular patterns, respectively. These two polymorphs are the first cases of one-component 3D XOFs reported up to now. The diameter of the channels spans the range 1.0–1.2 nanometers and the percentage of empty volume varies from 22.5 to 25.0. The powder X-ray analysis shows that no change of the crystalline architecture occurs when the solvent is removed from the cavities or when it refills the emptied cavities. Computation of MEP surfaces of the zwitterions confirms the relevance of the iodine σ-holes in determining the charge assisted XBs driving the formation of the observed architectures, and QTAIM/NCIplot analyses prove the robustness of these interactions. Importantly, while the porosity of 1a and 1b is not particularly high, the relevance of the obtainment of these systems consists in the fact that their crystalline architectures afford a proof of concept for the potential of zwitterionic aryliodonium tectons in affording robust one-component 3D XOFs. Greater porosity might be obtained by modifying the structure of the used zwitterionic aryliodonium tectons, e.g., by employing zwitterions with an increased separation between the anionic and cationic sites.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2096180, 2096183, 2096184 and 2096193. For ESI and crystallographic data in CIF or other electronic format see https://doi.org/10.1039/d2sc00892k |
This journal is © The Royal Society of Chemistry 2022 |