1. Introduction
Crystallization in the form of hydrates is inherent to many types of organic materials. These compositions are particularly significant to pharmaceutical science since the hydrated forms of active species may be more suited for applications with respect to long-term stability, solubility and drug performance (Morris, 1999). The general interest in organic hydrates is reflected in the increasing research interest in the crystal chemistry of such materials (Sanii et al., 2021). A recently reported data-driven and topological mapping approach allows the prediction of hydrate crystal structures from an anhydrous framework (Hong et al., 2022). It is believed that up to one third of organic compounds can afford hydrates (Stahly, 2007), while a recent analysis identified as many as 23698 unique hydrate structures within 286752 entries matching the search criteria (Werner & Swift, 2020). The existence of hydrated and anhydrous structures for a given compound is also remarkable. Only 1476 hydrates had a corresponding anhydrate in the database, although this number was possibly underestimated due to a bias of the Cambridge Structural Database (CSD; Groom et al., 2016) toward structures that have only been crystallized once (Cruz-Cabeza et al., 2015). Such an underestimation could be even more significant when considering the existence of hydrated and anhydrous structures as a special issue of a more general case of poly-stoichiometry of hydrates. Therefore, organic compositions with different substrate-to-aqua ratios are apparently uncommon, as revealed by the analysis of Basford (2021). Only seven compounds within a set of 6575 entries adopted at least three different hydrate stoichiometries. One may suppose that the number of such systems is still limited primarily at the expense of metastable and elusive higher hydrate forms, the existence of which can be foreseen with a degree of certainty for many cases.
Poly-stoichiometry of hydrates and the isolation of higher hydrates are particularly interesting in view of nucleation and crystal growth, since the crystallization of anhydrous or less hydrated forms inevitably succeeds significant desolvation of the dissolved species. Following the well-known Ostwald's rule of stages, one can intuitively assume initial crystallization of metastable highly hydrated phases. They may also be a primary outcome of crash crystallization. Such higher hydrates may inherit some features of the substrate–water interactions in solution, as a sufficient amount of the coordinated solvent remains an integral part of the crystal (Nangia & Desiraju, 1999). This situation is illustrated by the crystal structure of hexamethylenetetramine hexahydrate, bearing a close resemblance to the solution environment of the substrate (Burton et al., 2009). Although crystallizations of hydrates often do not obey the rule of stages (Tian et al., 2010), such a scenario may be more realistic for compounds that are most prone to forming hydrates. For example, 2,4-dihydroxybenzoic acid usually forms a hemihydrate, but fast cooling of hot solutions (crash crystallization) led to a very unstable monohydrate (Braun et al., 2011).
The propensity of organic compounds for hydrate formation was associated with a (dis)balance of the available hydrogen-bond donor and acceptor sites (Desiraju, 1991), while a subsequent study suggested correlation rather with a total number of such functionalities within the molecule (Infantes et al., 2007). With the aim of isolating high hydrates, we have explored the behaviour of 3,3′,5,5′-tetranitro-4,4′-bipyrazole (H2Tnbpz), which is an excellent candidate for providing a set of valuable pre-requisites. First, this substrate perfectly fits the above criteria, in view of the exceptionally rich hydrogen-bonding functionality accompanied by a striking mismatch between the numbers of available donors and acceptors, i.e. two N—H donors and ten O- and N-atom acceptors. Second, conformational flexibility of the bipyrazole molecule is largely beneficial for polymorphism and pseudopolymorphism, which are well known for the 3,3′,5,5′-tetramethyl-4,4′-bipyrazole prototype (H2Me4bpz; Boldog et al., 2003). An even more important functional feature is that the highly acidic N—H groups are very strong donors, while the NO2 groups, as well as the weakly basic dinitropyrazole N atoms, are only poor acceptors of hydrogen bonding. In accordance with Etter's hydrogen-bonding rule (Etter, 1990), the need for the bonding of strong proton donors and acceptors here are to be implemented preferably in the case of the substrate–aqua pair. Therefore, unlike the N—H⋯N hydrogen-bonded patterns seen for the H2Me4bpz systems (Boldog et al., 2001), rather N—H⋯OH2 interactions may occur in the case of H2Tnbpz. One can identify this substrate as being particularly prone either to hydrate formation or to hydrate poly-stoichiometry. Following these inputs, we succeeded in the crystallization of two metastable high hydrates and report their structures here. 3,3′,5,5′-Tetranitro-4,4′-bipyrazole and its ionic derivatives have attracted attention as perspective energetic materials, which combine sufficient performance and low sensitivity for safe applications (Gospodinov et al., 2024).
3. Results and discussion
The anhydrate H2Tnbpz (CSD refcode PITGEH) crystallized with difficulty, using the Schlenk technique, from hot high-boiling aromatic solvents after azeotropic removal of any dissolved water (Domasevitch et al., 2019). A strong trend for the hydration is best reflected by the formation of stable H2Tnbpz·H2O (refcode PITGIL), which dominated crystallization of the substrate from every examined common solvent under contact with ambient air. In spite of such invariant isolation of the monohydrate, further insights into this presumably rich hydrate system are possible when exploiting the special protolytic properties of H2Tnbpz. Due to the appreciable acidity of 3,5-dinitropyrazole groups [pKa = 3.14 for 3,5-dinitropyrazole; Janssen et al., 1973], the substrate readily forms nitropyrazolate salts, such as Li2(Tnbpz)(H2O)4 (Domasevitch et al., 2023). Singly charged hydrogen bipyrazolate anions, [H(Tnbpz)]−, are also known (Gospodinov et al., 2020). Aqueous solutions of MI[H(Tnbpz)] may be prepared by neutralization with appropriate amounts of alkali metal carbonates, but the stability of such salts is rather different (Fig. 1). Hydrogen bipyrazolates were isolated with the largest Rb and Cs cations (Domasevitch & Ponomarova, 2021), but in the case of MI = Li–K, they readily dismutated to afford soluble normal salts and the excess amount of the acid was deposited in the form of H2Tnbpz·H2O. However, we have found that the outcome of the crystallization depends on the cooling rate, and fast cooling of hot solutions allows the preparation of the metastable highly hydrated species (1) and (2).
| Figure 1 Crystallization of H2Tnbpz·nH2O as a result of prototropic dismutation of hydrogen bipyrazolates (M = Na). The postulated aggregation of the H(Tnpbz)− anions by extensive mutual NO2/NO2 interactions and hydration of the pyrazole sites are well reflected by the structures of the metastable high hydrates (1) and (2). |
The molecular structures of the title compounds are shown in Figs. 2 and 3. The main geometries of the organic frames agree well with the parameters for H2Tnbpz and its monohydrate (Domasevitch et al., 2019). In particular, clear differentiation of the angles at the ring N atoms suggests neutral pyrazole structures with localized and immobile H atoms [N—N(H)—C = 110.35 (15)–111.34 (14)° and N(H)—N—C = 103.74 (14)–105.11 (15)°]. Certain conjugation between the nitro groups and heterocycles is indicated by the nearly flat structure of the dinitropyrazole fragments, with corresponding NO2/ring dihedral angles in a range 1.36 (9)–13.52 (11)° [mean 5.09 (14) and 6.97 (12)° for (1) and (2), respectively]. The dihedral angles between two pyrazole rings, however, are significant, being 66.95 (7)° for (1) and 60.38 (7)° for (2). For H2Tnbpz·H2O, this angle was even more appreciable [78.99 (6)°]. In spite of such a twisted conformation of the organic molecules, in each of the hydrates, the components afford relatively dense packing, as is indicated by packing indices of 73.6 for the monohydrate, 74.2 for (1) and 73.0 for (2), which are slightly higher than the value of 72.2 for the anhydrate structure. These values approach the upper limit of the 65–75% range expected for organic solids (Dunitz, 1995).
| Figure 2 The molecular structure of H2Tnbpz tetrahydrate, (1), showing the atom and ring labelling, and with displacement ellipsoids drawn at the 50% probability level. Dotted lines indicate hydrogen bonding. |
| Figure 3 The molecular structure of H2Tnbpz pentahydrate, (2), showing the atom and ring labelling, and with displacement ellipsoids drawn at the 50% probability level. One water molecule is unequally disordered (0.86:0.14) and the minor contribution (O2WA) is isotropic. Dotted lines indicate hydrogen bonding. |
Progressive hydration results in a very illustrative evolution of crystal patterns (Fig. 4). H2Tnbpz·H2O was a genuine pocket hydrate incorporating isolated water molecules, whereas the two present metastable materials are channel hydrates showing extensive aqua–aqua hydrogen bonding. An increase in the molar fraction of water molecules does not lead to their uniform distribution in the lattice, but primarily causes clustering and the formation of extended hydrate networks. An immediate result of the clustering is the assembly of one-dimensional (1D) hydrate tapes in tetrahydrate (1), but with a higher number of water molecules in pentahydrate (2), the system develops an assembly of two-dimensional (2D) hydrate layers. In both structures, the mutual bonding of water molecules is directional and follows a standard geometry with typical O⋯O separations, which are 2.685 (2)–2.949 (3) Å for (1) (Table 2) and 2.766 (2)–2.939 (2) Å for (2) (Table 3). The longer O2W—H4⋯O5Wi hydrogen bond in (1) [O⋯O = 2.949 (3) Å and H⋯O = 2.33 (3) Å; symmetry code: (i) −x, −y + 1, −z + 1], however, exists as a shorter branch of a bifurcated interaction with aqua and nitro acceptors.
D—H⋯A | D—H | H⋯A | D⋯A | D—H⋯A | N1—H1N⋯O1W | 0.96 (2) | 1.68 (2) | 2.634 (2) | 177 (3) | N5—H2N⋯O4Wi | 0.97 (2) | 1.76 (2) | 2.723 (2) | 172 (3) | O1W—H1⋯O3W | 0.87 (2) | 2.00 (2) | 2.866 (2) | 178 (3) | O1W—H2⋯O2W | 0.87 (2) | 1.97 (2) | 2.771 (2) | 153 (3) | O2W—H3⋯N6ii | 0.87 (2) | 2.10 (2) | 2.942 (2) | 164 (3) | O2W—H4⋯O4Wiii | 0.88 (2) | 2.07 (2) | 2.939 (2) | 168 (3) | O3W—H5⋯O3iv | 0.86 (2) | 2.31 (3) | 3.116 (2) | 156 (4) | O3W—H6⋯O2Wiii | 0.88 (2) | 1.90 (2) | 2.766 (2) | 170 (3) | O4W—H7⋯O3W | 0.88 (2) | 1.92 (2) | 2.779 (2) | 166 (3) | O4W—H8⋯O4v | 0.88 (2) | 2.62 (4) | 3.136 (2) | 118 (3) | O4W—H8⋯O6iv | 0.88 (2) | 2.39 (4) | 3.051 (2) | 132 (4) | Symmetry codes: (i) ; (ii) ; (iii) ; (iv) ; (v) . | |
D—H⋯A | D—H | H⋯A | D⋯A | D—H⋯A | N1—H1N⋯O1W | 0.92 (3) | 1.70 (3) | 2.610 (2) | 173 (3) | N5—H2N⋯O4Wi | 0.87 (2) | 1.85 (2) | 2.723 (2) | 179 (2) | O1W—H1⋯O2W | 0.87 (2) | 1.99 (2) | 2.833 (2) | 163 (2) | O1W—H2⋯O5Wii | 0.86 (2) | 1.84 (2) | 2.685 (2) | 170 (3) | O2W—H3⋯N6iii | 0.88 (2) | 2.21 (3) | 3.000 (2) | 151 (3) | O2W—H4⋯O5Wiv | 0.83 (2) | 2.33 (3) | 2.949 (3) | 131 (3) | O2W—H4⋯O1v | 0.83 (2) | 2.57 (3) | 3.241 (3) | 139 (3) | O3W—H5⋯N2 | 0.85 (2) | 2.07 (2) | 2.885 (2) | 160 (2) | O3W—H6⋯O2W | 0.83 (2) | 1.99 (2) | 2.802 (2) | 166 (2) | O4W—H7⋯O3Wiv | 0.88 (2) | 1.87 (2) | 2.735 (2) | 173 (2) | O4W—H8⋯O3W | 0.85 (2) | 1.91 (2) | 2.749 (2) | 176 (3) | O5W—H9⋯O4W | 0.85 (2) | 1.94 (2) | 2.784 (2) | 175 (3) | O5W—H10⋯O3 | 0.86 (2) | 2.19 (3) | 2.910 (2) | 141 (3) | O5W—H10⋯O8iv | 0.86 (2) | 2.67 (3) | 3.065 (2) | 109 (2) | Symmetry codes: (i) ; (ii) ; (iii) ; (iv) ; (v) . | |
| Figure 4 Evolution of the crystal patterns adopted by H2Tnbpz·nH2O as a response to the increased number of water molecules. (a) The pocket hydrate PITGIL (n = 1), (b) the channel hydrate (1) (n = 4) embedding 1D aqua tapes and (c) the channel hydrate (2) (n = 5) with 2D aqua layers. For (b) and (c), the infinite aqua connectivities are orthogonal to the drawing plane. |
The infinite hydrate tapes in (1), running along the b direction in the crystal, consist of distorted seven-membered rings sharing three water molecules between adjacent units, which is a rare T7(3) hydrate topology according to Infantes & Motherwell (2002). The H2Tnbpz molecules are integrated into the tape by double hydrogen bonding, but this scheme is relevant only for rings B, which are N—H⋯O hydrogen-bond donors and O—H⋯N hydrogen-bond acceptors (Fig. 5). The former bond is particularly strong [N5⋯O4Wi = 2.723 (2) Å], whereas relative weakness of the latter interaction [O2W⋯N6ii = 2.942 (2) Å; symmetry codes: (i) −x, −y + 1, −z + 1; (ii) x + , −y + , −z + 1] reflects the low basicity of the dinitropyrazole N atom. The interactions of the A rings are more complicated and only strong N—H⋯OH2 hydrogen bonds are retained [N1⋯O1W = 2.634 (2) Å]. The resulting dinitropyrazole/aqua ensemble is bound to the adjacent dinitropyrazole ring as a triple donor of lone pair–π-hole interactions (Fig. 5). In particular, the N2 atom is situated nearly above the ring centroid, with an N2⋯Cg(ring A)iii distance of 3.2492 (16) Å, while nitro–nitro O3⋯N4iii and aqua–nitro O1W⋯N3iii separations are as short as 3.031 (2) and 3.029 (2) Å, respectively [symmetry code: (iii) −x + , y + , z.]
| Figure 5 (a) Projection of the structure of (1) on the ab plane, showing the topology of the aqua tape (which is highlighted with a blue strip) and the accommodation of the H2Tnbpz molecules through bonding of the B rings. (b) The bonding mode of the A rings with a set of lone pair–π-hole interactions involving also the O1W water molecule. [Symmetry codes: (i) −x, −y + 1, −z + 1; (iii) −x + , y + , z; (vi) x − , −y + , −z + 1.] |
The hydrogen-bonded pattern of the H2Tnbpz molecules in (2) is simpler due to a uniform function of both pyrazole rings. Similar to the B rings in the tetrahydrate, strong N—H⋯O hydrogen bonds to water molecules [N⋯O = 2.610 (2) and 2.723 (2) Å] are accompanied by weaker O—H⋯N interactions with aqua donors [N⋯O = 2.885 (2) and 3.000 (2) Å]. With such different donor and acceptor bonding, each pyrazole ring is installed at the hydrate layer (Fig. 6) and with two pairs of such bonds, the H2Tnbpz molecules are embedded between two hydrate layers separated by 11.2164 (5) Å, which is the a parameter of the unit cell.
| Figure 6 (a) Projection of the structure of (2) on the bc plane, showing the 2D aqua connectivity and how it integrates the organic molecules by multiple hydrogen bonding. (b) Two H2Tnbpz molecules embedded between two successive hydrate layers, with the aqua shell wrapping the hydrophilic dinitropyrazole sites. [Symmetry codes: (ii) −x + 1, y − , −z + ; (iv) x, −y + , z − .] |
Corrugation of the layer is conditioned by the need for the most effective interactions between the water molecules and the hydrophilic pyrazole sites. In fact, the water molecules tend to envelop dinitropyrazole groups, forming a certain sphere segment around them [Fig. 6(b)]. In this view, the array is reminiscent of the structures of clathrate hydrates incorporating the developed water frameworks. Such a structure could be assumed as a further possible step of the pattern evolution, with a higher water-to-substrate ratio. For example, the layered structure of 3-pyrroline trihydrate is similar to (2), but in the case of the pentahydrate, a genuine clathrate hydrate was observed (Rzepiński et al., 2019).
However, unlike many relatively small molecules, such as 3-pyrroline, total encapsulation within the hydrate shell is hardly possible due to the very extensive and energetically favourable interactions of the H2Tnbpz molecules themselves. Mutual stacking of the NO2 groups, as well as the stacking of NO2 and pyrazole groups, which are kinds of peculiar lone pair–π-hole interactions (Bauzá et al., 2015), are particularly prevalent and such bonding may be traced either for the anhydrate or any of the three hydrates. With a cut-off of 3.30 Å, as many as five short N⋯O contacts involving the nitro groups and also one O2N⋯OH2 contact of this type, assemble the molecules into tapes along the b axis in (1) (Fig. 7). We note the local motif with four doubly stacked NO2 groups and O7 atoms involved in three N⋯O contacts in the range 3.008 (2)–3.277 (2) Å (Table 4). In the case of (2) (Fig. 8), infinite double stacks along the c axis are generated by (N4/O3/O4)/(N4/O3/O4)vii interactions, with N⋯O = 3.019 (2) Å [symmetry code: (vii) x, −y + , z + ]. These separations are comparable to the short NO2/NO2 stacking seen in Li2(Tnbpz)(H2O)4 [N⋯O = 3.0349 (15)–3.0887 (15) Å; Domasevitch et al., 2023], although they are slightly longer than similar parameters for H2Tnbpz itself [2.911 (2) Å; Domasevitch et al., 2019]. However, all these interactions are highly directional, as is indicated by the nearly orthogonal orientation of the N⋯O axes with respect to the planes of the acceptor NO2 groups (Table 4). The NO2/pyrazole stacking in every case is selective to the N(H) sites, with the N⋯O separations down to 3.063 (2) Å for (1).
Compound | O-atom donor | Group | O⋯N | O⋯plane | φ | O⋯N (nitro) | | | | | | (1) | O3 | (C3/N4/O3/O4)iii | 3.031 (2) | 2.8316 (16) | 69.10 (3) | | O6 | (C6/N8/O7/O8)vii | 3.002 (2) | 2.8835 (18) | 73.85 (3) | | O7 | (C4/N7/O5/O6)ix | 3.008 (2) | 2.8201 (19) | 69.64 (4) | | O7 | (C1/N3/O1/O2)ix | 3.277 (2) | 3.098 (2) | 70.98 (5) | | O7 | (C6/N8/O7/O8)viii | 3.085 (2) | 2.951 (2) | 73.05 (5) | | O1W | (C1/N3/O1/O2)iii | 3.029 (2) | 2.8845 (19) | 72.23 (4) | (2) | O3 | (C3/N4/O3/O4)iv | 3.034 (2) | 2.9979 (16) | 81.16 (2) | | O4 | (C3/N4/O3/O4)vii | 3.019 (2) | 2.9866 (16) | 81.65 (2) | | O6 | (C6/N8/O7/O8)viii | 3.217 (2) | 2.911 (2) | 64.81 (4) | | O7 | (C1/N3/O1/O2)ix | 3.175 (2) | 3.0653 (18) | 74.89 (3) | | O1W | (C1/N3/O1/O2)vi | 3.132 (2) | 2.9153 (19) | 68.56 (4) | | | | | | | O⋯N(H) (ring) | | | | | | (1) | O2 | (N5/N6/C4/C5/C6)i | 3.144 (2) | 3.0726 (16) | 77.77 (3) | | O5 | (N5/N6/C4/C5/C6)xiii | 3.063 (2) | 2.770 (3) | 64.73 (5) | (2) | O1 | (N5/N6/C4/C5/C6)iv | 3.087 (2) | 2.738 (3) | 62.49 (5) | | O5 | (N5/N6/C4/C5/C6)x | 3.0857 (19) | 2.781 (2) | 64.32 (5) | | O8 | (N1/N2/C1/C2/C3)ix | 3.190 (2) | 3.1830 (19) | 86.20 (4) | Notes: O⋯plane is a distance of the O-atom donor to the mean plane of the nitro (pyrazole) group and φ is the angle of the O⋯N axis to the plane of the nitro (pyrazole) group. [Symmetry codes for (1): (iii) −x + , y + , z; (vii) x, y + 1, z; (viii) −x, −y, −z + 1; (ix) x, y − 1, z, (xiii) −x, y + , −z + . Symmetry codes for (2): (iv) x, −y + , z?1/2 [PLEASE CLARIFY]; (vi) −x + 1, −y, −z + 1; (vii) x, −y + , z + ; (viii) x, y, z − 1; (ix) x, y, z + 1; (x) −x, −y, −z + 1.] | |
| Figure 7 Lone pair–π-hole interactions of the NO2 groups in the structure of (1), with the O7 atom acting as a triple donor of such bonds. The grey circle indicates a special kind of bonding involving A rings, which is further detailed in Fig. 5(b). [Symmetry codes: (iii) −x + , y + , z; (vii) x, y + 1, z; (viii) −x, −y, −z + 1; (ix) x, y − 1, z.] |
| Figure 8 Lone pair–π-hole interactions of the NO2 groups in the structure of (2), arranging the molecules into 1D stacks along the c direction. [Symmetry codes: (iv) x, −y + , z − ; (vi) −x + 1, −y, −z + 1; (vii) x, −y + , z + ; (viii) x, y, z − 1; (ix) x, y, z + 1.] |
Dense NO2/NO2 networks give rise also to less favourable close O⋯O contacts. For (1), they are only slightly below the sum of the van der Waals radii, e.g. O7⋯O7viii = 2.970 (3) Å [symmetry code: (viii) −x, −y, −z + 1], whereas (2) reveals an exceptionally short nonbonded contact O4⋯O6vii of 2.5963 (18) Å [symmetry code: (vii) x, −y + , z + ]. This may be compared with inter-polyhedral O⋯O contacts of 2.687 (2) Å in the much more robust covalent framework of NiSO4 (Wildner, 1990). Such an arrangement is likely essential for the instability of the present hydrate.
One can note that the NO2/NO2 stacking tolerates well the progressively increased intensity of the substrate–aqua interactions in line with the increased number of the water molecules (H2Tnbpz·nH2O; n = 0, 1, 4 or 5). This lone pair–π bonding is highly competitive to weak hydrogen bonds with nitro groups (Bauzá et al., 2017) and therefore different kinds of such stacks remain intact even in the present high hydrates. In spite of the gradually increased number of identified O—H⋯O(NO2) hydrogen bonds, which is two for the monohydrate, three for (1) and four for (2) (Table 5), these interactions represent only the weakest hydrogen bonds in the structures [O⋯O = 2.910 (2)–3.241 (3) Å] and remain in the shadow of the apparently stronger mutual interactions of the nitro groups. This is contrary to the mutual hydrogen bonding of H2Tnbpz molecules. Thus, four N—H⋯O/O⋯H—N hydrogen bonds per molecule in H2Tnbpz itself were completely eliminated even in the case of the monohydrate, whereas two N—H⋯N/N⋯H—N hydrogen bonds still retained in H2Tnbpz·H2O, but disappear irrevocably upon further hydration, in favour of stronger N—H⋯OH2 and O—H⋯N hydrogen bonds.
Bond type | PITGEH | PITGIL | (1) | (2) | NH⋯N/N⋯HN | 2 | 2 | – | – | NH⋯O/O⋯HN | 4 | – | – | – | NH⋯OW | – | 1 | 2 | 2 | N⋯HOW | – | 1 | 2 | 2 | O⋯HOW | – | 2 | 3 | 4 | N⋯O/O⋯N (NO2) | 8 | 4 | 7 | 6 | N⋯O/O⋯N (ring) | 10 | 10 | 4 | 6 | O2N⋯OW | – | 1 | 1 | 1 | π⋯OW | – | – | – | 1 | O(N)⋯π/π⋯O(N) | 2 | 2 | 4 | – | N⋯NO2/O2N⋯N | 2 | – | – | – | O⋯O (NO2)a | 4 | 4 | 6 | 3 | Note: (a) the number of O⋯O contacts below 3.05 Å adopted by nitro groups. | |
The supramolecular interactions in the H2Tnbpz systems were further assessed by Hirshfeld surface analysis (Spackman & Byrom, 1997; McKinnon et al., 2004; Hirshfeld, 1977) performed with CrystalExplorer17 (Turner et al., 2017). The 2D fingerprint plots for the anhydrate PITGEH [Fig. 9(a)] and individual H2Tnbpz molecules in the monohydrate PITGIL [Fig. 9(b)] and the present tetra- [Fig. 9(c)] and pentahydrates [Fig. 9(d)] indicate the prevalence of the hydrogen-bond and O⋯N/N⋯O interactions, while general features for the partial contributions of the different contacts are very informative for the entire series (Table 6). The exceedingly large fraction of O⋯O contacts for H2Tnbpz (32.9%) indicates essential contraction upon hydration, down to 29.7% for the monohydrate, 17.9% for (1) and 21.1% for (2). However, in the latter two cases, the packing generated shorter O⋯O contacts. This is particularly relevant for (2) and is reflected by a short spike pointing to the lower left, at di + de = 2.6 Å [Fig. 6(d)]. The possible reasons for the ease of hydration of H2Tnbpz, as well as for the instability of the present high hydrates, may be at least in part associated with the unfavourable O⋯O contacts between the nitro groups. A similar trend for decreased contributions of O⋯N/N⋯O contacts occurs primarily at the expense of less favourable NO2/pyrazole interactions, since the fractions of entirely O⋯N/N⋯O (nitro) contacts themselves are almost insensitive to hydration (7.8–8.9%). The nature of the latter contacts is identical across the series and they are identified in the plots by symmetrical (about the diagonal where di = de) pairs of features. Therefore, both donor and acceptor sites of such bonds, namely, the O and N atoms of the nitro groups, appear within the individual H2Tnbpz molecules. A slight shortening of these features reflects a certain weakening of the O⋯N bonds, which are 2.90 Å for H2Tnbpz, 3.00 Å for (1) and 3.05 Å for (2). A decrease in the contribution of the O⋯O and O⋯N/N⋯O (pyrazole) contacts coincides with a rapid growth of contributions from the contacts with H atoms, which are only 14.1% for H2Tnbpz and 26.7% for H2Tnbpz·H2O, but are 38.5% for (1) and even 40.8% for (2). Similar to Li2(Tnbpz)(H2O)4 (Domasevitch et al., 2023), it may be postulated that the water hydrogen-bond donors have only a minor influence on the lone pair–π-hole interactions of the nitro groups, but rather they help to avoid less favourable nitro-O⋯O contacts. Unlike the NO2⋯NO2 bonds, O⋯H/H⋯O bonds for the hydrates are visualized by asymmetric fingerprint plots, clearly suggesting the prevalence of the interactions with water molecules. The upper spikes are most sharp and their lengths indicate particularly short separations of H⋯O = 1.65 Å, corresponding to the strongest N—H⋯OW hydrogen bonds. Lower spikes indicate OH2⋯O (nitro) bonds and they are perceptibly shorter (2.10–2.20 Å) and more diffuse. That the direct N—H⋯N hydrogen bonding between the pyrazole rings is irrelevant to the high hydrates coincides with the evolution of the upper spikes in the N⋯H/H⋯N plots. A very diffuse collection of points at large di + de seen for (1) completely disappears in the plot for the pentahydrate, leaving only a single lower sharp spike, which corresponds to OH2⋯N bonds with H⋯N = 2.00 Å.
Contacts | PITGEH | PITGIL | (1) | (2) | All⋯O,N | 80.7 | 69.3 | 56.7 | 54.5 | All⋯H | 14.1 | 26.7 | 38.5 | 40.8 | | | | | | O⋯H/H⋯O | 20.3 | 27.0 | 37.9 | 35.3 | N⋯H/H⋯N | 5.2 | 8.7 | 5.0 | 8.8 | C⋯H/H⋯C | 0.0 | 0.3 | 0.4 | 1.2 | N⋯O/O⋯N | 27.8 | 21.9 | 21.2 | 19.3 | C⋯O/O⋯C | 11.0 | 8.1 | 8.5 | 10.3 | C⋯N/N⋯C | 0.0 | 0.3 | 2.3 | 0.1 | O⋯O | 32.9 | 29.7 | 17.9 | 21.1 | N⋯N | 2.7 | 3.1 | 3.6 | 0.5 | H⋯H | 0.0 | 0.8 | 3.2 | 3.4 | | |
| Figure 9 2D fingerprint plots for the individual H2Tnbpz molecules, reflecting the subtle structural changes seen upon progressive hydration: (a) anhydrate PITGEH; (b) H2Tnbpz·H2O PITGIL; (c) H2Tnbpz·4H2O, (1); (d) H2Tnbpz·5H2O, (2), and delineated into the principal contributions of O⋯O, O⋯N/N⋯O, O⋯H/H⋯O and N⋯H/H⋯N contacts. For details on the other contributors to the surface areas, see Table 6. |
The intermolecular interaction energies were calculated using the CE B3LYP/6 31G(d,p) energy model in CrystalExplorer17 (Turner et al., 2017), where a set of symmetry-unique intermolecular paths were considered for the organic–organic and organic–aqua interactions for (1) and (2) (Table 7). From the plethora of intermolecular interactions, the highest energies (−49.3 to −56.1 kJ mol−1) correspond to strong N—H⋯O hydrogen bonds established by highly acidic pyrazole N—H donors and the primary contributor to the total energy is the electrostatic component (up to −83.7 kJ mol−1). These data are consistent with the energetics of azole–aqua interactions, which is Etot = −23.9 kJ mol−1 for relatively weak donor pyrazole itself, but −41.4 kJ mol−1 for the pair sustained by acidic pentazole and water (Chopra et al., 2018). In the case of hydrogen bonding to weakly basic pyrazole N-atom acceptors, Etot does not exceed −17.8 kJ mol−1. However, such interactions are superior to the O—H⋯O hydrogen bonds with the nitro O-atom acceptors (Table 6). Beyond the particular N—H⋯O hydrogen bonds, the strongest organic–aqua interactions are lone pair–π-hole bonds with nitro N-atom acceptors, with Etot approaching −26.4 kJ mol−1 in response to a higher contribution of the dispersion components. These values agree with the range of −8.5 to −33.6 kJ mol−1 reported for different nitro derivatives (Daszkiewicz, 2013).
Path | Typea | R (Å) | Eele | Epol | Edis | Erep | Etot | (1) | | | | | | | | Org⋯Orgix | A | 6.19 | −16.8 | −3.2 | −29.2 | 18.9 | −33.9 | Org⋯Orgi | B | 6.29 | −16.8 | −3.1 | −27.0 | 16.8 | −33.3 | Org⋯Orgviii | C | 8.30 | −2.9 | −1.2 | −14.1 | 8.6 | −10.9 | Org⋯Orgx | D | 6.86 | 1.3 | −4.1 | −28.1 | 19.2 | −14.3 | Org⋯Orgxi | E | 8.71 | −7.4 | −1.6 | −8.8 | 4.2 | −14.1 | Org⋯(O1W) | NH⋯OH2 | 5.54 | −82.2 | −20.1 | −7.9 | 85.0 | −56.1 | Org⋯(O4W)i | NH⋯OH2 | 5.59 | −62.8 | −16.1 | −6.9 | 56.8 | −49.3 | Org⋯(O2W)vi | OH⋯N | 5.90 | −28.2 | −4.7 | −5.9 | 33.4 | −17.8 | Org⋯(O4W)xii | OH⋯O | 5.69 | −12.0 | −1.6 | −5.9 | 11.1 | −12.1 | Org⋯(O3w)xii | OH⋯O | 5.09 | 0.4 | −0.8 | −4.1 | 3.0 | −1.8 | Org⋯(O1W)x | H2O⋯NO2 | 4.50 | −19.1 | −4.1 | −10.3 | 10.7 | −25.5 | | | | | | | | | (2) | | | | | | | | Org⋯Orgiv | A | 6.66 | −17.1 | −2.7 | −23.2 | 14.3 | −31.5 | Org⋯Orgx | B | 6.10 | −18.8 | −3.4 | −25.6 | 12.5 | −37.0 | Org⋯Orgviii | C | 6.33 | −14.7 | −2.9 | −29.6 | 22.3 | −29.8 | Org⋯Orgxi | D | 8.35 | −2.8 | −1.3 | −12.0 | 5.1 | −11.2 | Org⋯(O1W) | NH⋯OH2 | 5.56 | −83.7 | −20.6 | −6.7 | 87.3 | −55.7 | Org⋯(O4W)i | NH⋯OH2 | 5.61 | −68.5 | −17.5 | −8.1 | 64.1 | −52.9 | Org⋯(O3w) | OH⋯N | 5.72 | −22.3 | −5.0 | −7.2 | 41.8 | −7.8 | Org⋯(O2W)xii | OH⋯N | 5.88 | −11.9 | −4.3 | −6.5 | 24.1 | −6.6 | Org⋯(O4W)xiii | OH⋯O | 7.10 | 1.1 | −0.4 | −2.3 | 1.6 | −0.1 | Org⋯(O5W) | OH⋯O | 6.77 | −14.9 | −2.1 | −3.5 | 16.7 | −10.0 | Org⋯(O2W)v | OH⋯O | 6.21 | −4.6 | −0.8 | −3.6 | 3.7 | −6.3 | Org⋯(O1W)vi | H2O⋯NO2 | 4.43 | −19.2 | −3.6 | −11.7 | 10.8 | −26.4 | Org⋯(O5W)vii | H2O⋯π | 4.29 | −4.2 | −1.3 | −7.2 | 3.2 | −9.6 | Notes: (a) For details of the interaction modes in (1), see Fig. 10, and for details of the interaction modes in (2), see Fig. 11; R is the distance between the centroids of the interacting molecules. Symmetry codes for (1): (i) −x, −y + 1, −z + 1; (vi) x − , −y + , −z + 1; (viii) −x, −y, −z + 1; (ix) x, y − 1, z; (x) −x + , y − , z; (xi) −x, y − , −z + ; (xii) x, −y + , z − . Symmetry codes for (2): (i) x − 1, −y + , z + ; (iv) x, −y + , z − ; (v) −x + 1, −y, −z; (vi) −x + 1, −y, −z + 1; (vii) x, −y + , z + ; (viii) x, y, z − 1; (x) −x, −y, −z + 1; (xi) −x, −y, −z + 2; (xii) x − 1, y, z + 1; (xiii) x − 1, y, z. | |
It is not surprising that the Edis contributors are essentially larger (up to −29.6 kJ mol−1) for the different kinds of organic–organic stacks sustained by multiple lone pair–π-hole bonds. As a result, in every case, the total interaction energies are large and negative. This provides a realistic model for mutual stackings of H2Tnbpz, which could be rival for interactions with the solvent by weaker O—H⋯O hydrogen bonds or single H2O⋯NO2 bonds. This emphasis is important when involving a certain pre-nucleation effect for the crystallization of the present high hydrates (1) and (2). In the case of (1), the highest energies are associated with the formation of pairs according to Type A (Etot = −33.9 kJ mol−1), with double O⋯N (nitro) bonding, and Type B (Etot = −33.3 kJ −1), with two inversion-related O⋯N (pyrazole) interactions (Fig. 10). Three different pairs (Types A–C; Fig. 11) occur for (2), with very comparable interaction energies up to −37.0 kJ mol−1 (Table 6). The energetics of the individual O⋯N (nitro) and O⋯N (pyrazole) interactions are very similar, as is suggested by the values calculated for Types C and E in (1), and Type D in (2) (Etot are about −11 to −14 kJ mol−1). The latter values slightly exceed the energies of most O—H⋯O hydrogen bonds with nitro acceptors. Therefore, the extensive organic–aqua interactions in (1) and (2) concern predominantly the outer hydrophylic pyrazole N/NH regions, whereas two kinds of lone pair–π-hole bonds with nitro groups provide the main motifs for the packing of shape-complementary H2Tnbpz molecules. Also, the accumulation of additional organic–aqua interactions at the expense of the more favourable organic–organic bonding may be involved as a structure-destabilizing factor when comparing tetrahydrate (1) and very unstable pentahydrate (2).
| Figure 10 The five principal pathways (Types A–E) of the intermolecular bonding represented by different kinds of lone pair–π-hole interactions in the structure of (1). [Symmetry codes: (i) −x, −y + 1, −z + 1; (viii) −x, −y, −z + 1; (ix) x, y − 1, z; (x) −x + , y − , z; (xi) −x, y − , −z + .] |
| Figure 11 The NO2/NO2 and NO2/pyrazole N⋯O interactions in the structure of (2). The double interactions of Type D are more distal [N⋯O = 3.432 (2) Å] and weaker. [Symmetry codes: (iv) x, −y + , z − ; (viii) x, y, z − 1; (x) −x, −y, −z + 1; (xi) −x, −y, −z + 2.] |
In brief, the present systems provide an attractive paradigm for the crystallization of organic hydrates. The sequence of appearance of the various forms, moving from the exclusively unstable and elusive pentahydrate to the more stable tetrahydrate and finally to the stable monohydrate, is entirely consistent with Ostwald's rule of stages. The library of metastable high hydrates developed in the present study is important for the understanding of nucleation and crystal growth. In terms of Steed & Steed (2015), the metastable phases may be viewed as `fossil relics of the fastest growing crystal nuclei' and therefore they possibly reflect some preferable stacking configurations of the molecules and modes of their bonding with the solvent. Since the twisted molecules of 3,3′,5,5′-tetranitro-4,4′-bipyrazole, as well as the corresponding anionic species, perfectly fit one another like puzzle pieces and are prone to the generation of dense stacks, the postulated effects of pre-nucleation could be particularly significant in the present case. Therefore, this system may offer new insights not only into the area of organic hydrates, but also to the related problems of polymorphism and pseudopolymorphism.