Open Access Article
Damir A.
Safin
*a,
Maria G.
Babashkina
a,
Michael
Bolte
b,
Mariusz P.
Mitoraj
*c and
Axel
Klein
*a
aInstitut für Anorganische Chemie, Universität zu Köln, Greinstrasse 6, D-50939 Köln, Germany. E-mail: damir.a.safin@gmail.com; axel.klein@uni-koeln.de; Fax: +49 221 4705196; Tel: +49 221 4702913
bInstitut für Anorganische Chemie J.-W.-Goethe-Universität, Frankfurt/Main, Germany
cDepartment of Theoretical Chemistry, Faculty of Chemistry, Jagiellonian University, R. Ingardena 3, 30-060 Cracow, Poland. E-mail: mitoraj@chemia.uj.edu.pl
First published on 7th July 2015
Reaction of the in situ deprotonated N-thiophosphorylated thiourea 2-MeO(O)CC6H4NHC(S)NHP(S)(OiPr)2 (HL) with MCl2 (M = ZnII, CdII) in aqueous ethanol leads to complexes of the formula [ML2]. Both compounds crystallise in the triclinic space group P
with Z = 2 and the metal cations are found in a tetrahedral S2S′2 coordination environment formed by the C–S and P–S sulfur atoms. The crystal structures reveal intramolecular N–H⋯O
C hydrogen bonds formed within the 2-MeO(O)CC6H4NH fragments. Both structures are further stabilised by intermolecular π⋯π stacking interactions, which are more efficient in [CdL2]. Here, a pronounced dimeric intermolecular aggregate is observed which goes along with a pronounced distortion of the chelate [(S)CNP(S)]− backbone of the ligand upon coordination to CdII as well as a significantly distorted coordination tetrahedron CdS2S′2. The aggregation is also reflected in the positive electrospray ionisation (ESI) mass spectrum of the CdII complex, which exhibits peaks for the dimeric cations [Cd2L3]+, [Cd2L4 + H]+ and [Cd2L4 + Na]+, while for the ZnII analogue only monomeric species were observed. Quantum chemical ETS-NOCV (ADF) calculations confirm the higher stability of dimers in [CdL2] compared with [ZnL2]. The π⋯π stacking interactions are prodominantly due to dispersion contributions, though the electrostatic and orbital interaction components are also important. QTAIM (ADF) type calculations additionally quantify the covalent and non-covalent interactions in the momomers.
In this contribution we describe the synthesis of new ZnII and CdII complexes with the N-thiophosphorylated thiourea 2-MeO(O)CC6H4NHC(S)NHP(S)(OiPr)2 (HL).9 We also describe a complete structural investigation of the obtained complexes [ZnL2] and [CdL2] both in solution and solid state together with their thermal properties. The experimental results were supported by detailed quantum chemical calculations.
The IR spectra of complexes are very similar and contain a band at about 560 cm−1 representing the P–S group of the anionic form L− with a delocalised 6 π electron S–C–N–P–S system (bond order ∼1.5).10 This band is shifted by ∼90 cm−1 to low frequencies compared with that in the spectrum of the parent HL.9 Further bands at 1530 and 1695 cm−1 correspond to the conjugated SCN fragment and C
O group, respectively. In addition, there is a broad intense absorption arising from the POC group at about 970–990 cm−1. Also, the characteristic band for the arylNH group is found at about 3215 cm−1.
The 31P{1H} NMR spectra of [ZnL2] and [CdL2] in CDCl3 exhibit a unique signal at 55.7 and 55.9 ppm, respectively, which indicates the exclusive presence of 1,5-S,S′-coordinated ligands in the complexes.9–11 The 1H NMR spectra of complexes in CDCl3 each contain one set of signals. The signals of the isopropyl CH3 protons are observed at about 1.40 ppm, while the aryl CH3 protons are at 3.90 ppm. The isopropyl CH(O) protons appear as a doublet of septets at 4.81 ppm with the characteristic coupling constants 3JPOCH = 10.4 Hz and 3JH,H = 6.1 Hz. The phenylene protons are observed as four multiplet peaks at about 7.05, 7.45, 8.00 and 8.65 ppm with the typical coupling constants 3JH,H ranging from 7.9 to 8.3 Hz. The spectra also contain a doublet for the arylNH protons, observed at 11.25 ppm with a coupling constant 4JPNCNH of 7.5 Hz. The latter is exclusively observed when the structure of the H–N–C–N–P fragment meets the so-called “W-criterion” (marked by red in Scheme 1).12 The low-field shift of the signal for the arylNH protons is due to the formation of intramolecular hydrogen bonds of the type arylN–H⋯O
C.
Crystals of the complexes were obtained by slow evaporation of the solvent from CH2Cl2–n-hexane solutions. Both structures represent spirocyclic chelates and were refined in the triclinic space group P
, each containing one independent molecule in the unit cell. In the crystal, both Δ and Λ enantiomers of the complexes are present due to the inversion centre in the non-chiral space group P
(Fig. S1 in the ESI†). Such pairs of enantiomers can be expected for nonplanar bis-chelate complexes.13 Both metal cations are found in a tetrahedral S2S′2 coordination environment formed by the C–S and P–S sulfur atoms (Fig. 1). The six-membered M–S–C–N–P–S metallocycles have an asymmetric boat form. The values of the endocyclic S–M–S angles are about 108.9° (Table 1) in the structure of [ZnL2] and very close to that of the ideal tetrahedron (109.5°). The same angles in the structure of [CdL2] deviate significantly with 101.0(1) and 104.3(1)°. The same trend is found for the exocyclic S(C)–M–S(C), S(P)–M–S(P) and S(C)–M–S(P) angles: the values fall in the range of 108.6–110.9° in the structure of the ZnII complex, while the corresponding values of about 116.3, 112.0, 107.5 and 115.6 were observed in the CdII analog. The two values for the N–P–S, Zn–S–C and Zn–S–P angles are very similar in [ZnL2] representing a quite symmetric surrounding. This is not the case for [CdL2] where the corresponding pairs show significant deviations (Table 1).
| [ZnL2] | [CdL2] | |
|---|---|---|
C S |
1.743(3), 1.747(2) | 1.770(4), 1.778(4) |
P S |
1.991(1), 1.991(1) | 1.972(1), 1.992(2) |
| P–N | 1.602(2), 1.603(2) | 1.583(3), 1.604(3) |
| C–N(C) | 1.356(3), 1.361(3) | 1.342(5), 1.351(4) |
| C–N(P) | 1.304(3), 1.305(3) | 1.291(5), 1.343(5) |
| M–S(C) | 2.3087(8), 2.3137(7) | 2.433(1), 2.513(1) |
| M–S(P) | 2.3365(7), 2.3505(7) | 2.531(1), 2.621(1) |
| S–C–N(C) | 110.5(2), 111.9(2) | 111.3(3), 111.9(3) |
| S–C–N(P) | 128.4(2), 129.2(2) | 128.4(2), 128.5(3) |
| N–C–N | 119.7(2), 120.3(2) | 119.7(3), 120.2(3) |
| N–P–S | 117.4(1), 117.5(1) | 112.4(1), 119.5(1) |
| C–N–P | 127.3(2), 128.4(2) | 127.2(2), 128.5(3) |
| M–S–C | 106.1(1), 108.9(1) | 97.1(1), 101.2(1) |
| M–S–P | 95.16(3), 97.34(3) | 93.63(5), 97.71(6) |
| S–M–Sendo | 108.89(2), 108.94(3) | 101.00(5), 104.32(4) |
| S(C)–M–S(C)exo | 110.77(3) | 116.32(4) |
| S(P)–M–S(P)exo | 108.71(3) | 112.04(5) |
| S(C)–M–S(P)exo | 108.62(3), 110.85(3) | 107.49(4), 115.62(4) |
The lengthening of the C–S and P–S and shortening of the C–N(P) and P–N bonds in both complexes, compared with the values for the parent ligand HL,9 are in agreement with the IR data. The relative higher symmetry of the coordination surrounding of the metal atoms in [ZnL2] compared with [CdL2] is also visible in the M–S(C) and M–S(P) bonds. While the each two values for the M–S(C) and M–S(P) bonds in [ZnL2] are very similar and range from 2.31 to 2.35 Å, the two M–S(C) and M–S(P) bonds in the structure of [CdL2] are rather dissimilar and range from 2.43 to 2.62 Å.
The arylNH protons are involved in intramolecular hydrogen bonds of the type arylN–H⋯O
C (Fig. 1 and Table 2) in line with the structure concluded from NMR (Scheme 1). A closer inspection of the crystal structures revealed further H⋯X hydrogen bonds (Table S1 in the ESI†), however, based on established criteria14 we consider them not to be determining the crystal or molecular structures.
| D–H⋯A | d(D–H) | d(H⋯A) | d(D⋯A) | ∠(DHA) | |
|---|---|---|---|---|---|
| [ZnL2] | N(2)–H(2)⋯O(3) | 0.79(3) | 2.02(3) | 2.677(3) | 142(3) |
| N(4)–H(4)⋯O(7) | 0.83(4) | 1.91(4) | 2.637(3) | 146(3) | |
| [CdL2] | N(2)–H(2)⋯O(3) | 0.89(4) | 1.94(4) | 2.680(4) | 140(4) |
| N(4)–H(4)⋯O(7) | 0.73(4) | 2.00(4) | 2.620(5) | 143(4) |
Both structures are further stabilised by an intermolecular π⋯π stacking interaction leading to dimers for both [ZnL2] and [CdL2] (Fig. 2, Table 3). The π⋯π stacking seems to be much more efficient in [CdL2] leading to an almost coplanar arrangement (interplanar angle α = 0.03°) and a quite short interplanar distance of 3.756(3) Å with a displacement angle β of 24.12°.15 For [ZnL2] a second, only slightly longer, π⋯π stacking interaction (Table 3) finally leads to a 1D polymeric chain (Fig. 3), which is not observed for the Cd derivative. This marked difference seems to be caused by the pronounced distortion of the chelate [(S)CNP(S)]− fragment of the ligands upon coordination in the structure of [CdL2] as well as the significantly distorted coordination tetrahedron CdS2S′2 (Fig. 4) discussed already above. Thus, the nature and especially the size of the metal cation (ZnII: 74 pm vs. CdII: 92 pm)16 drive the supramolecular aggregation of molecules in the structures of [ZnL2] and [CdL2]. A similar influence of ZnIIvs. CdII was found for the formation of supramolecular coordination complexes of the N-thiophosphorylated 2,5-dithiobiurea [NHC(S)NHP(S)(OiPr)2]2.5e The dinuclear mesocate structure was formed upon reacting with ZnII, while the tetranuclear nanoscaled aggregate was isolated in the reaction with CdII. Furthermore, the idea to use the different sizes of these two d10 configured so-called “spherical ions”17 to control the supramolecular aggregation has been worked out recently.17,18 However, in most cases the observed supramolecular aggregate is the product of several interactions and carefuly ligand design is necessary.
![]() | ||
| Fig. 3 1D polymeric chain formed via π⋯π stacking interactions of the π⋯π stacked dimers in the structure of [ZnL2]. Hydrogen atoms except NH were omitted for clarity. | ||
![]() | ||
| Fig. 4 Molecule overlay of [ZnL2] (red) and [CdL2] (black). Hydrogen atoms were omitted for clarity. | ||
| Cg(I) | Cg(J) | Cg–Cg (Å) | α (°) | β (°) | |
|---|---|---|---|---|---|
| a Symmetry codes: #1 1 − x, 1 − y, 1 − z; #2 −x, −y, −z. Cg(3): C(11)–C(12)–C(13)–C(14)–C(15)–C(16), Cg(4): C(31)–C(32)–C(33)–C(34)–C(35)–C(36). b Symmetry codes: #1 1 − x, 1 − y, 1 − z. Cg(4): C(31)–C(32)–C(33)–C(34)–C(35)–C(36). | |||||
| [ZnL2] | Cg(3) | Cg(3)#1 | 4.1083(15) | 0.02 | 34.75 |
| Cg(4) | Cg(4)#2 | 4.112(2) | 0.03 | 32.77 | |
| [CdL2] | Cg(4) | Cg(4)#1 | 3.756(3) | 0.03 | 24.12 |
Frequently, CdII exhibits higher coordination numbers than ZnII thus changing the crystal structure.19 It is also quite common that going from ZnII to CdII leads to a higher distortion of the same coordination polyhedron.5c,20 In contrast to this, the two here presented compounds represent a rare example in which the increasing size of the central metal ion well-nigh “switches on” a very strong π⋯π stacking interaction leading to a marked change in the supramolecular interactions and crystal structure; e.g., the cell volume of 2006.2(2) Å3 for [ZnL2] is reduced to 1949.9(7) Å3 for [CdL2].
The bulk samples of [ZnL2] and [CdL2] were studied by means of X-ray powder diffraction analysis (Fig. 5). The experimental X-ray powder patterns are in agreement with the calculated powder patterns obtained from a single crystal X-ray analysis, showing that the bulk materials of [ZnL2] and [CdL2] are free from phase impurities.
![]() | ||
| Fig. 5 Calculated (black) and experimental (red) X-ray powder diffraction patterns of [ZnL2] (bottom) and [CdL2] (top). | ||
The formation of a rather stable dimeric aggregate via efficient intermolecular π⋯π stacking interactions is also reflected in the positive electrospray ionisation (ESI) mass spectrum of the CdII complex, which exhibits intense peaks for the dimeric cations [Cd2L3]+, [Cd2L4 + H]+ and [Cd2L4 + Na]+ (Fig. 6). For the ZnII analogue only monomeric species were found in the positive ion ESI mass spectrum.
![]() | ||
| Fig. 6 The positive ion ESI-mass spectrum of [CdL2] (top). Calculated (black) and experimental (red) patterns for the [Cd2L3]+ (middle) and [Cd2L4 + H]+ (bottom) cations. | ||
The thermal properties of [ZnL2] and [CdL2] in air atmosphere were studied by means of TG analyses in order to determine their respective stabilities (Fig. 7). The compound [ZnL2] is stable up to about 100 °C and decomposed in two steps. [CdL2] is stable up to 140 °C and decomposed in three steps with the third step to be poorly defined. The observed final residues of 11.6 and 16.4% are in excellent agreement with the calculated 11.54 and 16.21% for ZnS and CdS, respectively. The formation of these sulfides, each exhibiting a hexagonal form (wurtzite for ZnS and greenockite for CdS, respectively), was proved based on the powder X-ray diffraction analysis with full correspondence with data calculated from a corresponding single crystal analysis (Fig. 8).21
![]() | ||
| Fig. 8 Calculated (black) and experimental (red) X-ray powder diffraction patterns of ZnS (bottom) and CdS (top) obtained from annealing of [ZnL2] and [CdL2], respectively, at 750 °C. | ||
Thus, the complexes [ZnL2] and [CdL2], being easily obtained, and air/mosture stable might be very suitable single source precursors for the formation of ZnS and CdS, respectively.
In order to shed some light on the nature of bonding between the monomers in [ZnL2] and [CdL2], we have applied the charge and energy decomposition method (ETS-NOCV)22 as implemented in the ADF2012.01 program suite.23 We have used DFT/BLYP-D3 when performing the ETS-NOCV calculations as it was shown that such computational details provides satisfactory results for the noncovalent interactions as compared with the accurate CCSD(T) results.24 The coordinates of the dimers were extracted from the crystal structures.
The total interaction energy (ΔEtotal) in the dimer of [CdL2] is significantly more pronounced compared to that of [ZnL2] (Table 4). These results are in line with the ESI-mass spectra (Fig. 6), which also suggest a higher stability of the dimer of [CdL2]. Decomposition of the π⋯π stacking interaction energies shows that the dispersion contribution (ΔEdispersion) is the most important for the overall stabilisation and it contributes as much as −31.6 kcal mol−1 for [CdL2] and −40.3 kcal mol−1 for [ZnL2] (Table 4). Quantitatively less important is the electrostatic term (ΔEelstat), which appeared to be quite similar for both dimers and of about −10.0 kcal mol−1 (Table 4). Finally, the least contribution is the orbital interaction (ΔEorb), which is approximately twice more important for the cadmium-containing system. In order to shed further light on the ΔEorb term we have plotted the overall deformation density (Δρorb) upon formation of dimers (Fig. 9). The formation of both dimers leads rather to internal polarisations within the monomers, no typical charge transfer between the stacking rings is observed. It is noticeable that changes in the electron density are not only within the stacking rings but they also cover further regions including the OiPr groups as well as the sulphur atoms. The domination of the dispersion contribution in such stacking interactions is in line with the literature.25
| [ZnL2] | [CdL2] | |
|---|---|---|
| ΔEtotal = ΔEelstat + ΔEPauli + ΔEorb + ΔEdispersion. | ||
| ΔEelstat | −9.7 | −10.6 |
| ΔEPauli | 20.8 | 26.0 |
| ΔEorb | −4.6 | −8.8 |
| ΔEdispersion | −31.6 | −40.3 |
| ΔEtotal | −25.2 | −33.7 |
We have further performed the geometry optimisation of the monomers of [ZnL2] and [CdL2] in the gas phase based on DFT/BLYP-D3/TZP (Fig. 10), and the calculated parameters are in qualitative agreement with the experimental values (Tables 1 and 2 and Table S2 in the ESI†).
![]() | ||
| Fig. 10 The optimised structures with the selected bond lengths (Å) for monomers of [ZnL2] (bottom) and [CdL2] (top) based on ADF/DFT/BLYP-D3/TZP. | ||
We have finally performed a preliminary study of bonding in the monomers of [ZnL2] and [CdL2] based on the Quantum Theory of Atoms in Molecules (QTAIM)26 method as implemented in the ADF program.23 It was found that in both cases not only the strongest coordinating bonds Zn–S and Cd–S are observed, but also less important secondary intramolecular non-covalent interactions of the types O⋯H–N, N⋯H–C and S⋯H–C. From density values at the bond critical points as well as distances one can further infer that the strength of these interactions decreases in the following order O⋯H–N > N⋯H–C > S⋯H–C. Finally, given the fact that the bond critical points implies that stabilisation exists from the electronic exchange channel between atoms, further studies are required to more deeply describe each contribution of intramolecular close contacts. Neverthless, we believe that these preliminary QTAIM based results indentify important factors that might influence the overall stability of the [ZnL2] and [CdL2] monomers (Fig. 11).
C hydrogen bonds favour the all-trans arrangement of the ligands H–N–C–N–P fragment (the so-called “W-criterion”). In the solid, the structure was elucidated by single crystal X-ray diffraction analysis, revealing that both compounds crystallised isostructural in the triclinic space group P
, each containing one independent molecule in the unit cell. The metal cations display a tetrahedral S2S′2 coordination environment formed by the C–S and P–S sulfur atoms and strong intramolecular N–H⋯O
C hydrogen bonds were observed within the 2-MeO(O)CC6H4NH fragments. Both structures are further stabilised by intermolecular π⋯π stacking interactions, which are more efficient in [CdL2] leading to isolated stongly connected dimers. In contrast to this, a second, weaker π⋯π stacking is observed in the Zn analogue which leads to a 1D polymeric structure in the solid. This difference between the two complexes is also reflected in the positive electrospray ionisation (ESI) mass spectrum of the CdII complex, which exhibits peaks for the dimeric cations [Cd2L3]+, [Cd2L4 + H]+ and [Cd2L4 + Na]+. Only peaks for the monomeric species were found in the ESI mass spectrum of the ZnII analogue. Theoretical calculations confirm a higher stability of [CdL2] compared with [ZnL2]. Furthermore, the π⋯π stacking interactions are prodominantly due to dispersion contributions, though the electrostatic and orbital interaction components are also important. Thus, there is strong evidence from experiment and theory that formation of the dimeric intermolecular aggregate in the structure of [CdL2] is driven by the pronounced distortion of the chelate [(S)CNP(S)]− backbones of ligands upon coordination to CdII as well as by a significantly distorted coordination tetrahedron CdS2S′2. Thus, the nature and most of all the size of the metal cation (ZnIIvs. CdII) drives the supramolecular aggregation of molecules in the structures of [ZnL2] and [CdL2].
:
4 mixture of CH2Cl2 and n-hexane.
S), 970, 987 (POC), 1531 (SCN), 1694 (C
O), 3211 (NH). 1H NMR δ (ppm): 1.37 (d, 3JH,H = 6.2 Hz, 24H, CH3, iPr), 3.92 (s, 6H, CH3, Me), 4.81 (d. sept, 3JPOCH = 10.4 Hz, 3JH,H = 6.1 Hz, 4H, OCH), 7.04 (d. t, 3JH,H = 8.2 Hz, 4JH,H = 1.1 Hz, 2H, p-H, C6H4), 7.44 (d. t, 3JH,H = 8.1 Hz, 4JH,H = 1.6 Hz, 2H, m-H, C6H4), 7.99 (d. d, 3JH,H = 8.0 Hz, 4JH,H = 1.5 Hz, 2H, o-H, C6H4), 8.65 (d, 3JH,H = 8.3 Hz, 2H, m-H, C6H4), 11.26 (d, 4JPNCNH = 7.4 Hz, 2H, arylNH). 31P{1H} NMR δ (ppm): 55.7. ES-MS positive ion, m/z (%): 454.2 (18.4) [ZnL]+, 866.4 (100) [ZnL2 + Na]+. ES-MS negative ion, m/z (%): 389.2 (100) [L]−, 1232.8 (37.2) [ZnL3]−. Anal. Calc. for C30H44N4O8P2S4Zn (844.27): C 42.68, H 5.25, N 6.64. Found: C 42.61, H 5.22, N 6.71%.
S), 970. 982 (POC), 1530 (SCN), 1695 (C
O), 3215 (NH). 1H NMR δ (ppm): 1.38 (d, 3JH,H = 6.2 Hz, 24H, CH3, iPr), 3.93 (s, 6H, CH3, Me), 4.81 (d. sept, 3JPOCH = 10.3 Hz, 3JH,H = 6.2 Hz, 4H, OCH), 7.04 (br. t, 3JH,H = 7.9 Hz, 2H, p-H, C6H4), 7.43 (d. t, 3JH,H = 8.0 Hz, 4JH,H = 1.4 Hz, 2H, m-H, C6H4), 7.98 (d. d, 3JH,H = 8.1 Hz, 4JH,H = 1.5 Hz, 2H, o-H, C6H4), 8.62 (d, 3JH,H = 8.2 Hz, 2H, m-H, C6H4), 11.28 (d, 4JPNCNH = 7.8 Hz, 2H, arylNH). 31P{1H} NMR δ (ppm): 55.9. ES-MS positive ion, m/z (%): 892.7 (63.4) [CdL2 + H]+, 1394.1 (100) [Cd2L3]+, 1783.2 (27.1) [Cd2L4 + H]+, 1806.7 (7.1) [Cd2L4 + Na]+. ES-MS negative ion, m/z (%): 389.1 (100) [L]−. Anal. Calc. for C30H44CdN4O8P2S4 (891.30): C 40.43, H 4.98, N 6.29. Found: C 40.52, H 5.06, N 6.37%.
, a = 10.4280(5), b = 11.4053(6), c = 17.3530(9) Å, α = 78.001(4), β = 84.932(4), γ = 85.131(4)°, V = 2006.19(18) Å3, Z = 2, ρ = 1.398 g cm−3, μ(Mo-Kα) = 0.949 mm−1, reflections: 32
272 collected, 7511 unique, Rint = 0.0510, R1(all) = 0.0454, wR2(all) = 0.1055.
, a = 9.7368(19), b = 10.046(2), c = 21.882(4) Å, α = 97.42(3), β = 99.37(3), γ = 109.62(3)°, V = 1949.9(7) Å3, Z = 2, ρ = 1.518 g cm−3, μ(Mo-Kα) = 0.907 mm−1, reflections: 18
402 collected, 7261 unique, Rint = 0.0853, R1(all) = 0.0539, wR2(all) = 0.1060.
Footnotes |
| † Electronic supplementary information (ESI) available: Fig. S1, showing the crystal packing of the complexes, Table S1 with selected structural data, and Table S2 with a comparison of calculated and experimental bond parameters. CCDC reference numbers 1062050 ([ZnL2]) and 1062051 ([CdL2]). For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c5dt02189h |
| ‡ Dedicated to Professor F. Ekkehardt Hahn on the Occasion of his 60th Birthday. |
| This journal is © The Royal Society of Chemistry 2015 |