Open Access Article
Md Gishana,
Puspendu Middyaa,
Michael G. B. Drewb,
Antonio Frontera
*c and
Shouvik Chattopadhyay
*a
aDepartment of Chemistry, Jadavpur University, Kolkata 700032, India. E-mail: shouvik.chattopadhyay@jadavpuruniversity.in
bSchool of Chemistry, The University of Reading, P.O. Box 224, Whiteknights, Reading RG6 6AD, UK
cDepartment of Chemistry, Universitat de les Illes Balears, Crta de Valldemossa km 7.5, 07122 Palma de Mallorca, Baleares, Spain. E-mail: toni.frontera@uib.es
First published on 27th September 2024
In this article, we present the synthesis and characterization of three zinc(II) complexes, [ZnII(HL1)2] (1), [ZnII(HL2)2]·2H2O (2) and [ZnII(HL3)2] (3), with three tridentate Schiff base ligands, H2L1, H2L2, and H2L3. The structures of the complexes were confirmed by single-crystal X-ray diffraction analysis. DFT calculations were performed to gain insights into the self-assembly of the complexes in their solid-state structures. Complex 1 exhibits dual halogen-bonding interactions (Br⋯Br and Br⋯O) in its solid-state structure, which have been thoroughly investigated through molecular electrostatic potential (MEP) surface calculations, alongside QTAIM and NCIPlot analyses. Furthermore, complex 2 features a fascinating hydrogen-bonding network involving lattice water molecules, which serves to link the [ZnII(HL2)2] units into a one-dimensional supramolecular polymer. This network has been meticulously examined using QTAIM and NCIplot analyses, allowing for an estimation of the hydrogen bond strengths. The significance of H-bonds and CH⋯π interactions in complex 3 was investigated, as these interactions are crucial for the formation of infinite 1D chains in the solid state.
Moreover, the exploration of supramolecular interactions, e.g., hydrogen bonding, π⋯π stacking, cation⋯π, C–H⋯π, lone pair⋯π and anion⋯π interactions, in the solid-state structure of many complexes is an interesting research area in structural biology and crystal engineering.38–47 Supramolecular chemists and theoretical chemists are also interested in the solid-state structure of transition and non-transition metal complexes to explore σ-hole interactions.48–50 Usually, the σ-hole indicates a region of lower electron density compared to its neighbouring regions. It is produced by the anisotropic distribution of the electron density along the σ-skeleton.51 In this article, we shall concentrate on the halogen-bonding interaction, which is an example of σ-hole interaction between the σ-hole created on a halogen atom in a molecular entity and an electron-rich region in another, or the same, molecular entity.52,53 Other types of σ-hole interactions are also known to crystal engineers, such as spodium bonding (for group 12), tetrel bonding (for group 14), chalcogen bonding (for group 16) etc.54–57
In the present work, three tridentate N,N,O-donor Schiff bases (H2L1, H2L2 and H2L3) were used to form three zinc complexes, [ZnII(HL1)2] (1), [ZnII(HL2)2]·2H2O (2) and [ZnII(HL3)2] (3) H2L1 = 2-(3-hydroxypropylimino)methyl-4,6-dibromo phenol; H2L2 = 2-(1-hydroxybutan-2-ylimino)methyl-6-ethoxyphenol; H2L3 = 2-(1-hydroxy-2-methylpropan-2-ylimino)methyl-4-bromophenol). Two types of halogen-bonding interactions (characterized by Br⋯Br and Br⋯O contacts) were observed in the solid-state structures of complex 1, as confirmed by DFT calculations. Furthermore, an interesting R44(8) supramolecular ring was observed in the solid-state structure of complex 2. This provides a high degree of stabilization and explains the basis of co-crystallization of complex 2 with two water molecules. The optimization and QTAIM analysis elucidated the strength and orientation of these H-bonds. For complex 3, the relative importance of H-bonds and CH⋯π interactions was investigated, as they determine the formation of infinite 1D chains in the solid state.
Yield: 472 mg (∼64%) based on zinc(II). Anal. calc. for C20H20Br4N2O4Zn (FW: 737.37): C, 32.58; H, 2.73; N, 3.80. Found: C, 32.60; H, 2.71; N, 3.81%. FT-IR (KBr, cm−1): 3576–3177 (νO–H), 3058–2841 (νC–H), 1620 (νCN). UV-vis, λmax (nm), [εmax (L mol−1 cm−1)] (MeOH): 376 (1.132 × 104).
Yield: 369.16 mg (∼65%) based on zinc(II). Anal. calc. for C26H34N2O8Zn (FW: 567.94): C, 54.98; H, 6.03; N, 4.93. Found: C, 55.01; H, 6.15; N, 4.92%. FT-IR (KBr, cm−1): 3477–3215 (νOH2, νOH), 2974–2854 (νC–H), 1610 (νCN). UV-vis, λmax (nm), [εmax (L mol−1 cm−1)] (MeOH): 375 (1.124 × 104).
Yield: 431 mg (∼71%) based on zinc(II). Anal. calc. for C22H26Br2N2O4Zn (FW: 607.64): C, 43.49; H, 4.31; N, 4.61. Found: C, 43.50; H, 4.32; N, 4.62%. FT-IR (KBr, cm−1): 3549–3107 (νOH), 2958–2926 (νC–H), 1622 (νC
N). UV-vis, λmax (nm), [εmax (L mol−1 cm−1)] (MeOH): 371 (1.048 × 104).
| Complex | 1 | 2 | 3 |
|---|---|---|---|
| a (ξ)R1 = Σ||Fo| − |Fe||/Σ|Fo|; (ψ) wR2 = Σw(|Fo|2 − |Fe|2)2/Σw|Fo|2)1/2. | |||
| Formula | C20H20Br4N2O4Zn | C26H34N2O8Zn | C22H26Br2N2O4Zn |
| Formula weight | 737.37 | 567.94 | 607.64 |
| Temperature (K) | 273 | 273 | 273 |
| Crystal system | Monoclinic | Orthorhombic | Monoclinic |
| Space group | C2/c | A b a 2 | P21/n |
| a (Å) | 21.561(1) | 21.494(2) | 12.5138(5) |
| b (Å) | 4.595(2) | 13.8469(13) | 11.9667(4) |
| c (Å) | 25.047(1) | 9.6418(1) | 16.7858(7) |
| β | 110.368(1) | (90) | 108.173(1) |
| Z | 4 | 4 | 4 |
| V (Å3) | 2326.44(17) | 2869.6(5) | 2388.27(16) |
| dcalc (g cm−3) | 2.105 | 1.315 | 1.690 |
| μ (mm−1) | 7.956 | 0.903 | 4.405 |
| F(000) | 1424.0 | 1192.0 | 1216.0 |
| Total reflections | 48 902 |
24 328 |
33 325 |
| Unique reflections | 3637 | 2572 | 5297 |
| Observed data [I > 2σ(I)] | 4428 | 2091 | 3785 |
| No. of parameters | 147 | 188 | 281 |
| R(int) | 0.0396 | 0.0494 | 0.0681 |
| ξR1, ψwR2 (all data) | 0.0541, 0.1008 | 0.0641, 0.1437 | 0.0923, 0.2249 |
| ξR1, ψwR2 [I > 2σ(I)] | 0.0396, 0.0915 | 0.0494, 0.1314 | 0.0681, 0.2004 |
| CCDC no. | 2346042 | 2346043 | 2346044 |
| Bonds | Complex 1 | Complex 2 | Complex 3 |
|---|---|---|---|
| Zn(1)–O(1) | 1.933(2) | 1.912(6) | 2.054(4) |
| Zn(1)–N(1) | 2.009(2) | 2.020(7) | 2.073(5) |
| Zn(1)–O(2) | — | — | 2.288(4) |
| Zn(1)–O(3) | — | — | 1.955(5) |
| Zn(1)–O(4) | — | — | 2.355(6) |
| Zn(1)–N(2) | — | — | 2.093(5) |
| Angles | Complex 1 | Complex 2 | Complex 3 |
|---|---|---|---|
| a Symmetry transformation = i = 1 − x, y, 1.5 − z, ii = 1 − x, 1 − y, z. | |||
| O(1)–Zn(1)–N(1) | 94.95(8) | 95.10(2) | 89.80(15) |
| O(1)–Zn(1)–O(1)i | 123.27(8) | — | — |
| O(1)–Zn(1)–O(1)ii | — | 117.60(2) | — |
| O(1)–Zn(1)–N(1)i | 111.87(8) | — | — |
| O(1)–Zn(1)–N(1)ii | — | 112.80(2) | — |
| N(1)–Zn(1)–N(1)i | 122.27(9) | — | — |
| N(1)–Zn(1)–N(1)ii | — | 125.30(3) | — |
| O(1)–Zn(1)–O(2) | — | — | 158.81(16) |
| O(1)–Zn(1)–O(3) | — | — | 100.46(17) |
| O(1)–Zn(1)–O(4) | — | — | 87.35(18) |
| O(1)–Zn(1)–N(2) | — | — | 99.10(17) |
| O(1)–Zn(1)–N(1) | — | — | 89.80(15) |
| O(2)–Zn(1)–O(3) | — | — | 96.87(18) |
| O(2)–Zn(1)–O(4) | — | — | 78.35(19) |
| O(2)–Zn(1)–N(1) | — | — | 74.31(16) |
| O(2)–Zn(1)–N(2) | — | — | 92.34(19) |
| O(3)–Zn(1)–O(4) | — | — | 166.10(2) |
| O(3)–Zn(1)–N(1) | — | — | 103.38(18) |
| O(3)–Zn(1)–N(2) | — | — | 92.02(18) |
| O(4)–Zn(1)–N(1) | — | — | 88.00(2) |
| O(3)–Zn(1)–N(2) | — | — | 92.02(18) |
| O(4)–Zn(1)–N(2) | — | — | 75.30(2) |
| N(1)–Zn(1)–N(2) | — | — | 160.56(17) |
:
1 molar ratio in methanol/acetonitrile according to literature methods.68–77 These Schiff bases were not isolated and purified; the methanol/acetonitrile solution of the Schiff bases were used in situ for the synthesis of Zn(II) complexes, [Zn(HL1)2] (1), [Zn(HL2)2]·2H2O (2) and [Zn(HL3)2] (3), by adding zinc(II) chloride and zinc(II) acetate dehydrate with constant stirring for the formation of the complexes 1, 2 and 3 respectively (in a 1
:
1 molar ratio), as shown in Scheme 1. Additionally, in complex 2, the cocrystallized water molecules are assumed to come from the precursor zinc(II) acetate dihydrate salt used during the course of the reaction. The synthetic routes to complexes 1, 2 and 3 are shown in Schemes 1–3, respectively.
N) stretching vibration appears at 1610–1622 cm−1.78,79 Bands in the region of 3058–2841 cm−1 may be assigned as C–H stretching vibrations.80,81 The appearance of a broad band around 3576–3107 cm−1 indicates the presence of O–H stretching vibrations82 of the hydroxyl groups of the complexes and the coordinated water molecule in the case of complex 2. The IR spectra of the complexes are shown in Fig. S1–S3 (ESI).†
The electronic spectra of complexes 1–3 were recorded in methanol. There is only one band in the UV-vis spectra of all the complexes. The band at 376 nm (in 1), 375 (in 2) and 371 (in 3) may be attributed as π → π* or n → π* transition.83–85 There is no band corresponding to d–d electronic transitions, as expected for zinc(II) complexes with d10 electronic configuration.86 The UV spectra of complexes 1, 2 and 3 are shown in Fig. S4–S6 (ESI),† respectively.
Upon excitation at 380 nm in methanol solution, complex 3 shows strong emission of around 461 nm.87 This emission may be tentatively attributed to the intra-ligand transitions modified by metal coordination. The fluorescence spectrum of the complex is shown in Fig. S7 (ESI).†
![]() | ||
| Fig. 1 The perspective view of complex 1. i = Symmetry transformation = 1 − x, y, 1.5 − z. Only the major component of the pendant disordered –CH2OH group is shown. | ||
The coordination polyhedron around the zinc(II) centre is best described as a distorted tetrahedral with a ZnN2O2 chromophore, furnished by two imine nitrogen atoms [N(1) and its symmetry-related counterpart, N(1)i] and two phenolate oxygen atoms [O(1) and its symmetry-related counterpart, O(1)i] of the tridentate Schiff base ligand (i = symmetry transformation = 1 − x, y, 1.5 − z). Two well-known geometries for a tetra-coordinated metal ion are square planar and tetrahedral. The geometry of the tetra-coordinated zinc(II) centre is distorted tetrahedral, as confirmed by calculating the τ4 index, which is defined as τ4 = {360 − (α + β)}/141, where α and β are the two largest ligand–metal–ligand angles in the coordination sphere.88 The τ4 index is equal to 0 for a perfectly square planar geometry. Similarly, the τ4 index is equal to 1 for a regular tetrahedral geometry. In the present case, the τ4 index is 0.81, which confirms that the geometry of the zinc(II) centre is distorted tetrahedral. The unsaturated six-membered chelate ring [Zn(1)–N(1)–C(7)–C(6)–C(1)–O(1)] presents an envelope conformation (Fig. S8, ESI†) with puckering parameters q = 0.1575(16) Å; θ = 55.6(7)°; ϕ = 5.3(10)°.89,90
In complex 1, there are Br⋯Br and Br⋯O halogen-bonding interactions involving neighboring molecules. The bromine atom Br(1) of the Schiff base unit forms halogen bonds (HaB) (Fig. 2) with the bromine atom Br(1)iii of the neighboring molecule (symmetry transformation = iii = 1.5 − x, −1/2 + y, 1.5 − z) to form a supramolecular dimer. On the other hand, the oxygen atom, O(2A), forms a HaB with the bromine atom, Br(2)iv, of the neighbouring molecule (symmetry transformation = iv = 1 − x, −y, 1 − z) to form a supramolecular dimer. The details of HaB interactions are presented in Table 4.
![]() | ||
| Fig. 2 The perspective views of the HaB interactions of complex 1. Hydrogen atoms are omitted for clarity. Symmetry transformation = iii = 1.5 − x, −1/2 + y, 1.5 − z, iv = 1 − x, −y, 1 − z. | ||
| Complex | Br⋯Br/O | Br⋯Br/O (Å) | Symmetry transformation |
|---|---|---|---|
| 1 | Br(1)⋯Br(1)iii | 3.7487(5) | iii = 1.5 − x, −1/2 + y, 1.5 − z |
| Br(2)iv⋯O(2A) | 3.768(6) | iv = 1 − x, −y, 1 − z |
In complex 1, three C–H⋯Br hydrogen-bonding interactions with the neighbouring molecules are present. The hydrogen atom H(9A) attached to the carbon atom C(9) of the Schiff base unit forms a hydrogen bond with the bromine atom Br(1)iii of the neighbouring molecule to form a supramolecular dimer (Fig. 3a) (symmetry transformation = iii = 1.5 − x, −1/2 + y, 1.5 − z). Another two hydrogen atoms H(5) and H(7), which are attached to the carbon atoms C(5) and C(7), respectively, of the Schiff base, form hydrogen bonds with the bromine atom Br(2)iv of the neighbouring molecule (symmetry transformation = iv = 1 − x, −y, 1 − z) to form a supramolecular dimer (Fig. 3b). The geometrical parameters are presented in Table 5.
| Complex | C–H⋯X | C–H (Å) | H⋯X (Å) | C–H⋯X (Å) | ∠C–H⋯X (°) | Symmetry transformation |
|---|---|---|---|---|---|---|
| 1 (X= Br) | C(7)–H(7)⋯Br(2)iv | 0.930 | 3.10 | 4.026(2) | 171.5 | iv = 1 − x, −y, 1 − z |
| C(5)–H(5)⋯Br(2)iv | 0.930 | 3.34 | 4.215(3) | 156.8 | ||
| C(9)–H(9A)⋯Br(1)iii | 0.970 | 3.53 | 4.135(3) | 122.3 | iii = 1.5 − x, −1/2 + y, 1.5 − z | |
| 3 (X = O) | O(2)–H(2O)⋯O(1)iii | 0.929 | 1.81 | 2.668(6) | 151.5 |
The zinc(II) center's coordination polyhedron can be best understood as a deformed tetrahedral with a ZnN2O2 chromophore, comprised by two imine nitrogen atoms [N(1) and its symmetry-related counterpart, N(1)ii] and two phenolate oxygen atoms [O(1) and its symmetry-related counterpart, O(1)ii] of the tridentate Schiff base ligand (ii = symmetry transformation = 1 − x, 1 − y, z). The zinc(II) center's deformed tetrahedral geometry is confirmed by the τ4 index of 0.83. The unsaturated six-membered chelate ring [Zn(1)–N(1)–C(9)–C(7)–C(8)–O(1)] presents an envelope conformation (Fig. S9, ESI†) with puckering parameters q = 0.182(5) Å; θ = 111(2)°; ϕ = 179(2)°.89,90
In complex 3, there is one O–H⋯O hydrogen-bonding interaction present with the symmetry-related neighbouring molecule. The hydrogen atom H(2O) attached to the oxygen atom O(2) of the Schiff base ligand (HL3)− forms a hydrogen bond with the oxygen atom O(1)iii of the symmetry-related neighbouring molecule, thus forming a one-dimensional supramolecular chain (Fig. 6) of complex 1 (symmetry transformation iii = 1.5 − x, −1/2 + y, 1.5 − z). The hydrogen bond distance and angle are presented in Table 5.
![]() | ||
| Fig. 6 Perspective view of the hydrogen-bonding interactions of complex 3. Only the relevant hydrogen atoms are shown, for clarity. Symmetry transformation iii = 1.5 − x, −1/2 + y, 1.5 − z. | ||
In complex 3, two C–H⋯π interactions are present with the symmetry-related neighbouring molecules. The aromatic ring R(1), containing carbon atoms C(12)–C(13)–C(14)–C(15)–C(16)–C(17), is involved in two C–H⋯π interactions with two hydrogen atoms, H(11A)iii and H(3)v, which are attached to carbon atoms C(11)iii and C(3)v, respectively, of the Schiff base ligand (HL3)− of the symmetry-related neighbouring molecules to form a supramolecular chain (Fig. 7). Symmetry transformation iii = 1.5 − x, −1/2 + y, 1.5 − z and v = 1.5 − x, −1/2 + y, 1.5 − z. The geometric parameters of C–H⋯π interactions are presented in Table 6.
| Complex | C–H⋯π | H⋯π (Å) | C–H⋯π (Å) | ∠C–H⋯Cg (°) | Symmetry transformation |
|---|---|---|---|---|---|
| 3 | C(11)iii–H(11A)iii⋯Cg(1) | 2.66 | 3.56 | 153.6 | iii = 1.5 − x, −1/2 + y, 1.5 − z |
| C(3)v–H(3)v⋯Cg(1) | 3.07 | 3.49 | 109.8 | v = 1.5 − x, −1/2 + y, 1.5 − z |
![]() | ||
| Fig. 8 Fingerprint plot: different contacts contributed to the total Hirshfeld surface area of complex 1. | ||
![]() | ||
| Fig. 9 Fingerprint plot: different contacts contributed to the total Hirshfeld surface area of complex 2. | ||
![]() | ||
| Fig. 10 Fingerprint plot: different contacts contributed to the total Hirshfeld surface area of complex 3. | ||
In complex 1, the O⋯H/H⋯O interactions comprise 16.7% of the Hirshfeld surfaces of the complex. This O⋯H/H⋯O interaction appears as two distinct spikes in the 2D fingerprint plot. The O⋯H interaction is represented by the lower spike (di = 1.1 and de = 0.75 Å), and the H⋯O interaction is represented by the upper spike (di = 0.75 and de = 1.1 Å), which can be viewed as bright-red spots on the dnorm surface. Similarly, the Br⋯H/H⋯Br and C⋯H/H⋯C interactions comprise 30.1% and 9.7% of the Hirshfeld surfaces of the complex. The Br⋯H interaction is represented by the lower spike (di = 1.8 and de = 1.1 Å), and the H⋯Br interaction is represented by the upper spike (di = 1.1 and de = 1.8 Å), which can be viewed as bright-red spots on the dnorm surface. The C⋯H interaction is represented by the lower spike (di = 2.2 and de = 1.5 Å), and the H⋯C interaction is represented by the upper spike (di = 1.5 and de = 2.2 Å), which can be viewed as bright-red spots on the dnorm surface. The Br⋯Br interactions comprise 3.6% and Br⋯O/O⋯Br interactions comprise 1.3% of the Hirshfeld surfaces of complex 1.
In complex 2, the O⋯H/H⋯O interaction comprises 17.6% of the Hirshfeld surfaces of the complex. This O⋯H/H⋯O interaction appears as two distinct spikes in the 2D fingerprint plot. The O⋯H interaction is represented by the lower spike (di = 1.45 and de = 1.15 Å), and the H⋯O interaction is represented by the upper spike (di = 1.15 and de = 1.45 Å), which can be viewed as bright-red spots on the dnorm surface. Similarly, the C⋯H/H⋯C interactions comprise 16.5% of the Hirshfeld surfaces of the complex. The C⋯H interaction is represented by the lower spike (di = 1.7 and de = 1.0 Å), and the H⋯C interaction is represented by the upper spike (di = 1.0 and de = 1.7 Å), which can be viewed as bright-red spots on the dnorm surface.
In complex 3, the O⋯H/H⋯O interactions comprise 6.2% of the Hirshfeld surfaces of the complex. This O⋯H/H⋯O interaction appears as two distinct spikes in the 2D fingerprint plots. The O⋯H interaction is represented by lower spike (di = 1.15 and de = 0.7 Å), and the H⋯O interaction is represented by upper spike (di = 0.7 and de = 1.15 Å), which can be viewed as bright-red spots on the dnorm surface. Similarly, the Br⋯H/H⋯Br and C⋯H/H⋯C interactions comprise 24.4% and 20.8% of the Hirshfeld surfaces of the complex. The Br⋯H interaction is represented by the lower spike (di = 1.85 and de = 1.1 Å), and the H⋯Br interaction is represented by the upper spike (di = 1.1 and de = 1.85 Å), which can be viewed as bright-red spots on the dnorm surface. The C⋯H interaction is represented by the lower spike (di = 1.6 and de = 1.1 Å), and the H⋯C interaction is represented by the upper spike (di = 1.1 and de = 1.6 Å), which can be viewed as bright-red spots on the dnorm surface.
![]() | ||
| Fig. 11 Partial views of the X-ray structures of complex 1 (a and b), 2 (c) and 3 (d). Distances in Å. H-atoms omitted for clarity (in complex 2, the water H atoms could not be located). | ||
Complex 2, lacking halogen atoms, demonstrates the significant impact of co-crystallization with two water molecules on the packing. Specifically, these water molecules, in synergy with the alcohol groups and the ethoxy–phenoxy O atoms of the Schiff-base complex, form a complex H-bonding network, although the hydrogen atoms could not be located because of disorder. However, it is clear from Fig. 11 that the hydrogen atoms form a central ring of hydrogen bonds that could run either clockwise or anticlockwise. This network is pivotal in the propagation of complex 2 into a 1D supramolecular polymer. The characteristics and significance of the H-bonding network are detailed through QTAIM and NCIplot analyses.
Complex 3, despite having a Br atom, does not engage in halogen bonding. Instead, it forms infinite 1D chains propagated by H-bonding interactions involving the phenoxido atoms as acceptor and the hydroxyl groups as donor (see Fig. 11d). The characteristics and energetic significance of such 1D chain have also been analyzed through QTAIM and NCIPlot analyses.
First, we calculated the molecular electrostatic potential (MEP) surfaces for complexes 1–3, presented in Fig. 12. For complex 1, the MEP's maximum value is found at the hydrogen atom of the hydroxyl groups (+50 kcal mol−1), while its minimum is at the phenoxy oxygen atoms (−33 kcal mol−1). In Fig. 12a, the MEP surface surrounding the bromine atoms, depicted using a reduced scale (±10 kcal mol−1), reveals the anisotropy around the bromine atom. This illustrates a σ-hole along the extension of the C–Br bond (+9.4 kcal mol−1) and a surrounding area of negative potential (−10.0 kcal mol−1). Notably, the σ-hole cone angle measures 35°, indicating that the C–Br⋯O interaction observed in complex 1 (Fig. 11b) can be classified as a σ-hole interaction. However, this interaction is anticipated to be relatively weak because it occurs in a region near the bromine atom where the MEP is approximately zero (electroneutral). In complex 2, the MEP's highest value is at the hydroxyl group (+55 kcal mol−1), with the lowest in the area influenced by both ethoxy and phenoxy oxygen atoms. Additionally, the oxygen atom of the hydroxyl group, serving as an intramolecular hydrogen bond donor, also displays a significantly negative MEP value (−39 kcal mol−1), highlighting its potential role as H-bond acceptor. Two different views of the MEP surface of complex 3 are presented in Fig. 12c. The MEP maximum is located at the hydroxyl group (+59 kcal mol−1), highlighting its potential role as an H-bond donor. This value is significantly more positive than that of complex 2, due to the coordination of the hydroxyl's oxygen atom to the Zn metal center, which increases the acidity of the hydroxyl proton. The minimum is located at the coordinated phenoxide oxygen atoms (−47 kcal mol−1), indicating their potential role as H-bond acceptors. The MEP is also negative at the belt of the bromine atom (−20 kcal mol−1) and over the aromatic ring (−17 kcal mol−1). Additionally, the MEP is positive at the hydrogen atoms of the alkyl chain connecting the imino nitrogen atoms of the Schiff-base ligand.
![]() | ||
| Fig. 12 MEP surfaces of complexes 1 (a), 2 (b) and 3 (c) using the 0.001 a.u. isovalue. The energies at selected points are given in kcal mol−1. | ||
Fig. 13a displays the QTAIM/NCIplot analysis for the Br⋯Br dimer of complex 1, extracted from the assembly depicted in Fig. 11a. The halogen-bonding interaction is confirmed by QTAIM analysis through a bond critical point (displayed as a fuchsia sphere) and a dashed bond path connecting the two Br atoms. The QTAIM analysis also uncovered a CH⋯Br interaction involving an aliphatic CH bond and several hydrophobic Csp2–H⋯H–Csp3 interactions. Small, green reduced density gradient (RDG) isosurfaces, aligned with the bond critical points, further elucidate the halogen-bonding (HaB), hydrogen-bonding (HB), and hydrophobic interactions, demonstrating their attractive nature (λ2 < 0). The interaction energy for this dimer is modest, at −4.9 kcal mol−1. To discern the contribution of halogen bonding, we calculated a mutated dimer where the propanol arm is substituted with a hydrogen atom, isolating the HaB interaction. This modification reduces the interaction energy to −1.6 kcal mol−1, indicating the HaB's contribution. This small interaction is consistent with the long Br⋯Br experimental distance and the small MEP value at the σ-hole.
Fig. 13b presents the combined QTAIM/NCIplot analysis of the self-assembled Br⋯O dimer. Here, bond critical points (BCPs) and bond paths link the Br and O atoms, highlighting the Br⋯O interaction. Additional CH⋯Br and Csp2H⋯HCsp2 interactions are identified, characterized by their respective BCPs, bond paths, and green RDG isosurfaces. The dimerization energy is −7.8 kcal mol−1, attributed to all identified interactions. To assess the HaBs' contribution specifically, a mutated dimer was analyzed, in which the propanol arms (indicated by small arrows) were replaced by hydrogen atoms, thus eliminating the HaBs. This alteration leads to an interaction energy of −4.5 kcal mol−1. The difference in dimerization energies (−3.3 kcal mol−1) is ascribed to the energy of the HaBs, with each Br⋯O contact contributing −1.65 kcal mol−1.
The hydrogen-bonding network of complex 2 was meticulously examined. The positions of the hydrogen atoms were determined through optimizations using the PBE0-D3/def2-SVP level of theory, while the coordinates of the heavier atoms were kept fixed. These positions of necessity break the two-fold symmetry of the crystal structure and indicate that the crystal structure must be disordered to accommodate two possible hydrogen bond patterns, clockwise as shown in Fig. 14, or anticlockwise. This disorder may well account for the hydrogen atoms not being located in the crystal structure. Subsequent single-point calculations were carried out at the PBE0-D3/def2-TZVP level. The energy associated with each hydrogen bond was approximated by the potential energy density (Vg) at the BCP identifying each H-bond. These findings, illustrated in Fig. 14, reveal an intermolecular OH⋯O hydrogen bond between the central water molecules and a sequential orientation of H-bonds, forming a R44(8) supramolecular ring. This includes the intramolecular OH⋯O hydrogen bond between the hydroxyl groups on the pendant arms. This configuration does not align with the Aba2 point group observed for this complex, suggesting a potential disorder in the hydrogen atoms' positioning, allowing for both clockwise and anticlockwise orientations of the R44(8) supramolecular ring. Furthermore, the water molecules form two bifurcated hydrogen bonds with the oxygen atoms of the ethoxy and phenoxy groups. Consequently, a total of eight strong hydrogen bonds are established, contributing to a significant stabilization energy of −44.2 kcal mol−1, which accounts solely for the intermolecular hydrogen bonds. This substantial interaction energy correlates with the dark-blue RDG isosurfaces characterizing the R44(8) supramolecular ring, with HB energies ranging from −8.3 to −8.9 kcal mol−1. The bifurcated hydrogen bonds, though weaker, exhibit energies between −3.5 and −6.3 kcal mol−1. This analysis confirms the importance of the H-bonding network dictating the solid-state architecture of complex 2.
![]() | ||
| Fig. 14 Combined QTAIM and NCIplot analyses (see Theoretical methods for settings) of the H-bonding network in complex 2. The energies of the HBs are indicated. | ||
Finally, a dimer extracted from the infinite 1D chain presented in Fig. 11c was analyzed. The QTAIM and NCIplot analyses support the strong nature of the OH⋯O bond (see Fig. 15), as also anticipated by the MEP surface analysis. The RDG isosurface characterizing the H-bond is strongly blue, and the calculated energy using the Vg predictor is −10.4 kcal mol−1, stronger than that determined for the H-bonding network of complex 3. The total dimerization energy is −15.4 kcal mol−1, primarily due to the OH⋯O H-bond. The QTAIM analysis also reveals the existence of an ancillary CH⋯O H-bond that accounts for −0.9 kcal mol−1. Interestingly, the QTAIM/NCIplot analysis also reveals the presence of a CH⋯π interaction characterized by a bond critical point (BCP) and bond path linking the aliphatic hydrogen atom to one carbon atom of the ring. The formation of this contact is consistent with the MEP surface analysis (Fig. 12c), which showed negative MEP over the π-system and positive MEP at the aliphatic hydrogen atoms. The π-nature of the interaction is clearly revealed by the NCIplot analysis, which shows a green RDG isosurface embracing most of the aromatic ring. The contribution of the CH⋯π interactions can be roughly estimated by the difference between the total contribution of the H-bonds (−11.3 kcal mol−1) and the dimerization energy (−15.4 kcal mol−1), resulting in −4.1 kcal mol−1.
![]() | ||
| Fig. 15 Combined QTAIM and NCIplot analyses (see Theoretical methods for settings) of a selected dimer of complex 3. The energies of the HBs are indicated in red. | ||
The implications of the results are not restricted to the development of the single crystals of different zinc complexes only. The results may point to the wider importance of understanding halogen and hydrogen bonding in the context of numerous physical and chemical phenomena in catalysis, materials science, ion transport and recognition etc. We trust that these findings will captivate not only the crystal engineering community but also researchers in supramolecular chemistry and theoreticians.
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 2346042–2346044. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4ra06217e |
| This journal is © The Royal Society of Chemistry 2024 |