Plan
Comptes Rendus

Structural studies of bis(histidinato)nickel(II): Combined experimental and computational studies
Comptes Rendus. Chimie, Volume 20 (2017) no. 5, pp. 467-474.

Résumé

A distorted octahedral nickel(II) complex, [Ni(2-amino-3-(1H-imidazol-4-yl)propanoic acid)2] (1), has been synthesized by a solvothermal method and characterized by single-crystal X-ray diffraction. Geometry optimization in the gas phase and pyridine together with Hirshfeld surface and reduced density gradient analyses reveal that this complex shows different distortions from octahedral in the gas, liquid, and solid phases. The reason seems to be because of the presence of two intramolecular NH⋯O weak interactions in the gas phase and two sets of rather strong intermolecular NH⋯O and CH⋯O interactions in the solid phase. Time-dependent density functional theory (TD-DFT) calculations suggest that these different distortions result in different electronic absorption spectra.

Métadonnées
Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crci.2016.12.006
Mots clés : Nickel(II) complex, Histidine, X-ray structure, DFT calculations, TD-DFT calculations, QTAIM analysis
Alireza Abbasi 1 ; Barzin Safarkoopayeh 1 ; Nikoo Khosravi 1 ; Alireza Shayesteh 1

1 School of Chemistry, College of Science, University of Tehran, Tehran, Iran
@article{CRCHIM_2017__20_5_467_0,
     author = {Alireza Abbasi and Barzin Safarkoopayeh and Nikoo Khosravi and Alireza Shayesteh},
     title = {Structural studies of {bis(histidinato)nickel(II):} {Combined} experimental and computational studies},
     journal = {Comptes Rendus. Chimie},
     pages = {467--474},
     publisher = {Elsevier},
     volume = {20},
     number = {5},
     year = {2017},
     doi = {10.1016/j.crci.2016.12.006},
     language = {en},
}
TY  - JOUR
AU  - Alireza Abbasi
AU  - Barzin Safarkoopayeh
AU  - Nikoo Khosravi
AU  - Alireza Shayesteh
TI  - Structural studies of bis(histidinato)nickel(II): Combined experimental and computational studies
JO  - Comptes Rendus. Chimie
PY  - 2017
SP  - 467
EP  - 474
VL  - 20
IS  - 5
PB  - Elsevier
DO  - 10.1016/j.crci.2016.12.006
LA  - en
ID  - CRCHIM_2017__20_5_467_0
ER  - 
%0 Journal Article
%A Alireza Abbasi
%A Barzin Safarkoopayeh
%A Nikoo Khosravi
%A Alireza Shayesteh
%T Structural studies of bis(histidinato)nickel(II): Combined experimental and computational studies
%J Comptes Rendus. Chimie
%D 2017
%P 467-474
%V 20
%N 5
%I Elsevier
%R 10.1016/j.crci.2016.12.006
%G en
%F CRCHIM_2017__20_5_467_0
Alireza Abbasi; Barzin Safarkoopayeh; Nikoo Khosravi; Alireza Shayesteh. Structural studies of bis(histidinato)nickel(II): Combined experimental and computational studies. Comptes Rendus. Chimie, Volume 20 (2017) no. 5, pp. 467-474. doi : 10.1016/j.crci.2016.12.006. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2016.12.006/

Version originale du texte intégral

1 Introduction

In recent years, there has been an increasing interest in the syntheses of functional metal complexes having catalytic [1,2] and magnetic [3,4] properties or biological activities [5,6]. These complexes can be formed by using various multidentate ligands. Although less attention has been paid to Ni(II) complexes, they have various applications such as antibacterial activity, CC and CS cross-coupling, double-bond functionalization, and CO2 reduction reaction [7,8]. Herein, we report synthesis of a novel Ni(II) complex, using histidine (2-amino-3-(1H-imidazol-4-yl) propanoic acid) as the chelating agent. Histidine is an α-amino acid with an imidazole functional group, first isolated by German physician Albrecht Kossel [9]. The conjugate acid (protonated form) of the imidazole side chain in histidine has a pKa of approximately 6.0. For pH less than 6, the imidazole ring is mostly protonated [10], and the imidazole ring bears two NH bonds and has a positive charge. The positive charge is equally distributed between the two nitrogen atoms and can be represented with two equally important resonant structures. The imidazole ring of histidine is aromatic at all pH values [11]. It can have π stacking interactions [12], complicated by the positive charge [13]. The title compound was synthesized using the solvothermal method [14]. Solvothermal processes, using low or medium temperatures, can be considered as alternatives to the high temperature methods based in particular on melt growth or vapor growth processes. Furthermore, theoretical calculation can be used for structural elucidation [15–17]. In the present work, we report the newly synthesized bis(histidinato)nickel(II) complex, followed by TD-DFT, quantum theory of atoms in molecules (QTAIM), and reduced density gradient calculations.

2 Experimental section

2.1 Materials

All chemicals were purchased from commercial sources and used as received.

Compound 1 was prepared from the mixture of NiCO3·xH2O (0.375 mmol), histidine (0.75 mmol), 4,4′-bipyridine (0.375 mmol), and methanol/H2O (2 mL with a ratio of 10:1). The suspension was stirred for 15 min and then transferred to a 23 mL Teflon-lined stainless steel autoclave followed by heating at 150 °C for 48 h. After slow cooling to room temperature, the product was filtered, washed with methanol, and dried in air. Elemental analysis calcd (%) for NiC12H16N6O4 + 1 H2O: C, 37.44%; H, 4.71%; N, 21.83%; O, 20.87%; found: C, 37.40%; H, 4.75%; N, 21.78%; O, 20.90%.

2.2 Crystal structure determination and refinement

The data were collected on Marresearch 345 dtb diffractometer using Mo Kα radiation (λ = 0.71073 Å) at 295 K. The structure was solved by direct methods using SHELXS-97 and the obtained model refined with SHELXL-97 [18]. All non-hydrogen atoms were refined anisotropically. Hydrogen atoms were added at ideal positions and refined using a riding model. The data were deposited in Cambridge Crystallographic Data Centre, with deposition number CCDC-979056. The Oak Ridge thermal ellipsoid plot (ORTEP) view of 1 and its selected crystallographic and experimental details are summarized in Fig. 1 and Table 1, respectively.

Fig. 1

ORTEP diagram of 1. Thermal ellipsoids are drawn at the 50% probability level.

Table 1

Crystal and structure refinement data for bis(histidinato)nickel(II) (1).

Table 1
Compound1
Empirical formulaC12H16N6NiO4
Molecular weight367.00
Crystal systemTetragonal
Space groupP43212
Temperature (K)293(2)
a = b (Å)7.3157(10)
c (Å)28.870(6)
α = β = γ (°)90.00
Cell volume (Å3)1545.1(4)
Z4
ρ (g cm−3)1.578
μ (mm−1)1.285
Total reflections5255
Unique reflections1465
Observed reflections [F2 > 2σ(F2)]1292
Data/restraints/parameters1465/0/121
Goodness-of-fit on F21.068
R [F2 > 2σ(F2)] (R1, wR2)0.0719, 0.1611
R (all data) (R1, wR2)0.0830, 0.1735

2.3 Thermal gravimetric analysis

The thermal properties of 1 were studied by Thermogravimetric analysis/Differential thermal analysis (TGA/DTA) in the temperature range 30–600 °C (Fig. 2). It showed that the complex is stable until 355 °C and starts to decompose at about 371 °C. The thermogram exhibits two distinct weight losses at 141 and 390 °C. The first weight loss at 141 °C could be ascribed to the loss of surface water. The second weight loss at 390 °C indicates loss of histidine ligands and decomposition of the complex.

Fig. 2

TGA/DTA characterization.

3 Computational details

To get insight into the electronic structure of 1, density functional calculations were carried out with the ORCA 3.0.3 program [19]. We used the PBE0 functional and the def2-TZVP (triple zeta valence polarized) basis sets to optimize the geometries in the gas phase and pyridine using the COnductor-like Screening MOdel (COSMO) solvation model [19]. Vibrational frequency calculations were carried on the optimized structures in solvent and gas phase, to confirm that both structures were minima. The electronic structure of the title compound was investigated by QTAIM and by reduced density gradient analysis using Multiwfn 3.7 program [20].

4 Results and discussion

4.1 Molecular structure

Geometry optimization was started from the single crystal X-ray experimental atomic positions and the optimized geometric parameters (bond distances and bond angles) are given in Table 2. All of the bond lengths are slightly overestimated compared with the experimental values. By comparing the bond lengths, maximum deviations of 0.052 and 0.010 Å for Ni2O8 were found for the optimized structure in the gas phase and pyridine, respectively. All bond lengths are in good agreement with their corresponding values from X-ray structural analysis, but the angles around central Ni(II) show rather high deviations. Maximum deviations of 11.35° and 6.08° for the O8aNi2O8 angle, 6.61° and 1.12° for N13aNi2N5, and 9.59° and 2.83° for N13aNi2N13 were obtained for the optimized structure in the gas phase and pyridine, respectively. Further investigations revealed that this deviation originated from two sets of weak and moderate O⋯H intramolecular and intermolecular interactions present in the gas and solid phases, respectively. The presence of two intramolecular NH⋯O weak interactions in the gas phase is responsible for increasing the N13aNi2N5 angle and decreasing the N13aNi2N13 angle (Fig. 3). These interactions do not exist in the solid phase because of the participation of respected oxygen atoms in a stronger intermolecular hydrogen bond with neighboring molecules. By optimization of the structure in pyridine solvent, the N13aNi2N5 and N13aNi2N13 angles get closer to the experimental values. However, the O8aNi2O8 angle still shows rather high deviations and we discuss the reasons for this discrepancy in the following sections.

Table 2

Selected bond lengths and angles around central nickel atom in 1.

Table 2
ExperimentalOptimized (gas)Optimized (pyridine)
Bonda (Å)
Ni2O82.0602.0262.070
Ni2N52.0992.1512.109
Ni2N132.0952.1152.104
Angle (degree)
O8aNi2O883.9795.3290.05
N13Ni2O880.7679.8979.34
N5Ni2O887.2886.1787.63
N13Ni2N586.3287.0987.68
N13aNi2N5103.17107.79104.29
N13aNi2N13165.24155.65162.41

a Because of P43212 space group contrary bonds are equivalent.

Fig. 3

The presence of two intramolecular weak interactions in the gas phase, which do not exist in solid phase.

4.2 Hirshfeld surfaces

To gain further insight into the intermolecular interactions of molecular crystals, Hirshfeld surface [21–23] and the related two-dimensional (2D)-fingerprint plot [24–26] were calculated using the Crystal Explorer program [27]. Hirshfeld surface analysis reveals that the dominant interactions in the solid phase are two sets of bifurcated intermolecular NH⋯O and CH⋯O bonds between imidazole ring from four of the neighboring molecules and carboxyl groups of 1 (Fig. 4). The related data of these hydrogen bonds are summarized in Table 3. The presence of these hydrogen bonds has two effects on the solid phase geometry: the intramolecular NH⋯O weak interactions disappear and the O8aNi2O8 angle becomes 83.97°.

Fig. 4

The presence of two sets of bifurcated intermolecular NH⋯O and CH⋯O bonds between imidazole ring from four of the neighboring molecules and carboxyl groups of 1.

Table 3

Hydrogen bonds specifications in the molecular crystal.

Table 3
H⋯A (Å)D⋯A (Å)D–H⋯A (degree)
N21bH⋯O11.9942.845169.50
N21bH⋯O82.5963.062115.09
N13bH⋯O12.2843.171167.41
C22cH⋯O12.1253.043169.14

The Hirshfeld surfaces were mapped with dnorm for 1 and the corresponding fingerprint plots are shown in Fig. 5. The dominant interaction in the molecular crystal is O⋯H and can be seen in the Hirshfeld surface as bright red areas. Two-dimensional fingerprint plots complement these surfaces, quantitatively summarizing the nature and type of intermolecular contacts experienced by molecules in the crystal. The O⋯H/H⋯O intermolecular interactions appear as distinct spikes in the 2D plot. The fingerprint plots can highlight atom pair close contacts and enable one to separate different interaction types, which overlap in the full fingerprint. The sum of O⋯H/H⋯O interactions comprises 36.5% of the Hirshfeld surfaces. Furthermore, Hirshfeld surfaces show no sign of the presence of the intramolecular NH⋯O weak interactions in the solid phase, although it is present in the gas phase.

Fig. 5

Hirshfeld surfaces mapped with dnorm for 1. (a) Fingerprint plots, (b) hole, and (c) O⋯H/H⋯O highlighted.

4.3 QTAIM and reduced density gradient

The QTAIM approach proposed by Bader is the most refined method for the analysis of the topology of the electron density, which is repeatedly used for revealing and enlightening the bonding situation between the atoms [28,29]. Here, we used this method to gain insight into binding characteristics and to confirm the presence of NH⋯O weak interactions in the gas phase structure. The data are shown in Table 4. Although QTAIM shows no critical point for the intramolecular NH⋯O weak interactions in the gas phase, the delocalization index for this bond in the gas phase shows some degree of electron sharing for a weak interaction between hydrogen and nitrogen. The delocalization index values in the gas phase are twice those in the solvent (Table 5), which confirms the presence of a weak interaction in the gas phase. Recently, it has been suggested that for weak intramolecular hydrogen bonds, the QTAIM analysis may show no bond critical points [30], and in these cases one might use the reduced density gradient analysis [31]. This analysis confirmed the presence of intramolecular NH⋯O weak interactions in the gas phase. Gradient isosurfaces for the NH⋯O weak interaction are shown in Fig. 6. The isosurface represents the weak interaction area clearly, although this interaction is categorized as weak. Repeating this analysis on the solvent-optimized structure revealed that such interaction is not present in solvent.

Table 4

Bond critical point electron densities (ρ), their Laplacian (∇2ρ), and ellipticities (ɛ).

Table 4
ρ2ρɛ
Bond(pyridine solvent)
Ni2O80.0630.3090.015
Ni2N130.0700.2940.072
Ni2N50.0660.3030.037
O8C300.362−0.5540.091
C30C160.251−0.6470.096
C16N130.255−0.6270.029
C16C410.245−0.6100.025
C30O10.392−0.4060.104
C41C200.259−0.6870.058
C20N50.311−0.9070.169
C20C400.331−0.9650.359
C22N50.360−1.1170.290
C22N210.329−0.8850.263
C40N210.305−0.7560.199
Bond(gas phase)
Ni2O80.0710.3460.005
Ni2N130.0690.2880.090
Ni2N50.0610.2830.050
O8C300.356−0.5770.101
C30C160.244−0.6100.095
C16N130.254−0.6230.030
C16C410.246−0.6160.031
C30O10.406−0.3110.130
C41C200.259−0.6890.050
C20N50.311−0.9120.158
C20C400.331−0.9710.358
C22N50.365−1.2000.302
C22N210.323−0.8540.264
C40N210.301−0.7250.209
Table 5

Calculated delocalization index for intermolecular hydrogen bonds.

Table 5
Bondδ(A⋯B)Gas phaseδ(A⋯B)Pyridine
N13aH⋯O80.02170.0117
N13H⋯O8a0.02170.0117
Fig. 6

Gradient isosurfaces (s = 0.35) showing weak intramolecular interaction present in the gas phase.

4.4 Spectroscopic studies

Because of insolubility of the titled compound in common solvents, solid state ultraviolet–visible (UV–vis) spectrum was recorded in the 200–700 nm range, and the electronic absorption spectrum was calculated using the gas phase and pyridine optimized geometries with TD-DFT method. We used the B3P86 functional and the def2-TZVP basis sets with C2 symmetry. The same calculations were carried out using the solid-state geometry without optimization (see Fig. 7).

Fig. 7

Calculated electronic absorption spectra using the gas phase, pyridine, and solid-state geometries with TD-DFT method using B3P86 functional.

The calculated electronic absorption spectra of 1 in the solid and gas phases and pyridine together with solid state UV–vis spectrum reveal that different distortions from octahedral geometry in the gas phase, pyridine, and the solid phase for 1 result in different electronic absorption spectrum in each phase. Calculation of electronic absorption spectra using solid-state geometry without optimization results in a spectrum, which is in good agreement with experimental data (Fig. 8). Repeating the same calculations with PBE0 functional gives similar result as B3P86.

Fig. 8

Solid-state UV–vis spectrum and calculated spectrum using B386 applying def2-TZVP basis set 200–800 nm.

Electronic absorption spectrum of 1 in solid state (Fig. 8) shows several bands roughly around 240, 290, 350, and 540 nm. By investigating the percentage weights of the molecular orbitals involved in electronic excitations, we found that in the band around 500–600 nm, there are contributions from three types of electronic excitations, that is, metal-to-ligand charge transfer, intraligand transitions, and metal d–d transitions. Furthermore, TD-DFT calculation results suggest that band around 500–580 mainly related to the transitions from HOMO-3 and HOMO-4 to LUMO, the band around 350 mainly related to the transitions from HOMO-12 and HOMO-11 to LUMO+1, and the band around 240 and 290 mainly related to the transition from HOMO and HOMO-1 to LUMO+2. The free histidine ligand has no absorption band at wavelength longer than 245 nm [32], which is consistent with our TD-DFT calculations on free ligand.

4.5 Frontier molecular orbitals

The frontier molecular orbitals (HOMO and LUMO) can determine molecular reactivity and its photophysical properties. The isodensity plots for the HOMOs and LUMOs orbitals are shown in Fig. 9. It can be seen that HOMO is mainly localized on the metal atom and on the 2-amino-propanoic acid part of the histidine ligand, whereas the LUMO is most on the imidazole ring and to a lesser extent on the metal. Because the first electron transfer occurs from the HOMO to LUMO, it can be inferred that the first low energy electronic transition is combined ligand to ligand charge transfer (LLCT) from organic moiety and d–d metal electron transition. By using HOMO and LUMO energy values for the molecule, electronegativity and chemical hardness and softness were estimated using well-known relations [33–36] and are reported in Table 6 along with the dipole moment of 1.

Fig. 9

The isodensity plots for HOMOs and LUMOs.

Table 6

Orbital energies, chemical hardness, chemical softness, electronegativity, and dipole moment of 1.

Table 6
PBE0/def2-TZVPGas phase (ɛ = 1)Pyridine (ɛ = 12.5)Water (ɛ = 80.40)
EHOMO (eV)−4.85−5.26−5.30
ELUMO (eV)−0.70−0.22−0.15
ΔEHOMO–LUMO4.145.045.14
Chemical hardness (η)2.072.522.57
Chemical softness (S)0.240.200.19
Electronegativity (χ)2.782.742.73
Dipole moment (Debye)16.023.324.6

5 Conclusions

A new bis(histidinato)nickel(II) has been synthesized, using solvothermal crystal growth. The resulting complex has distorted octahedral geometry. The title compound shows different octahedral distortions in the gas, liquid, and solid phases using geometry optimization in the gas phase and pyridine together with Hirshfeld surface and reduced density gradient analyses. This can be related to the presence of two weak intramolecular NH⋯O weak interactions in the gas phase and two sets of rather strong intermolecular NH⋯O and CH⋯O interactions in the solid phase. Although intramolecular NH⋯O weak interactions in the gas phase are not confirmed by QTAIM analysis, the reduced density gradient and delocalization indices approve this weak interaction. Furthermore, TD-DFT calculations suggest that these different distortions result in different electronic absorption spectrum in each phase.

Acknowledgments

The authors appreciate the University of Tehran for financial support.


Bibliographie

[1] A.K. Sengupta; O.P. Bajaj; U. Chandra J. Indian Chem. Soc., 55 (1978), p. 962

[2] H. Singh; L.D.S. Yadav; B.K.J. Battacharya J. Indian Chem. Soc., 56 (1979), p. 1013

[3] S. Giri; H. Singh; L.D.S. Yadav; R.K. Khare J. Indian Chem. Soc., 55 (1978), p. 168

[4] A. Cansiz; S. Servi; M. Koparir; M. Altintas; M. Digrak J. Chem. Soc., 23 (2001), p. 237

[5] I. Kostova Curr. Med. Chem., 13 (2006), pp. 1085-1107

[6] T. Ishiyama; J. Takagi; J.F. Hartwig; N. Miyaura Angew. Chem., Int. Ed., 41 (2002), p. 3056

[7] J. Breitenfeld; R. Scopelliti; X. Hu Organometallics, 31 (2012), p. 2128

[8] A.M.A. Alaghaz; M.E. Zayed; S.A. Alharbi; R.A. Ammar; A. Elhenawy J. Mol. Struct., 1084 (2015), p. 352

[9] J.D. Kopple; M.E. Swendseid J. Clin. Invest., 55 (1975), p. 881

[10] F.A. Carey, McGraw-Hill (2000), p. 868

[11] M. Agnieszka; K.W. Janina; K. Katarzyna J. Mol. Struct., 655 (2003), p. 397

[12] L. Wang; Na. Sun; T. Simon; Z. Xuejun; D.R. Benson Biochemistry, 45 (2006), p. 13750

[13] R.H. Blessing; E.L. McGandy J. Am. Chem. Soc., 94 (1972), p. 4034

[14] R.A. Laudise; R.L. Barns IEEE Trans. Ultrason. Ferr., 35 (1988), p. 277

[15] F. Senn; L. Helm; A. Borel; C.A. Daul C. R. Chim., 15 (2012), p. 250

[16] A. Elkechai; A. Boucekkine; L. Belkhiri; D. Hauchard; C. Clappe; M. Ephritikhine C. R. Chim., 13 (2010), p. 860

[17] D. Farmanzadeh; T. Abdollahi; M. Jafary C. R. Chim., 19 (2016), p. 579

[18] G.M. Sheldrick Acta Crystallogr., A, 64 (2008), p. 112

[19] F. Neese Comput. Mol. Sci., 2 (2012), p. 73

[20] T. Lu; F. Chen J. Comput. Chem., 33 (2012), p. 580

[21] M.A. Spackman; D. Jayatilaka Cryst. Eng. Commun., 11 (2009), p. 19

[22] F.L. Hirshfeld Theor. Chim. Acta, 44 (1977), p. 129

[23] H.F. Clausen; M.S. Chevallier; M.A. Spackman; B.B. Iversen New J. Chem., 34 (2010), p. 193

[24] A.L. Rohl; M. Moret; W. Kaminsky; K. Claborn; J.J. McKinnon; B. Kahr Cryst. Growth Des., 8 (2008), p. 4517

[25] A. Parkin; G. Barr; W. Dong; C.J. Gilmore; D. Jayatilaka; J.J. McKinnon; M.A. Spackman; C.C. Wilson Cryst. Eng. Commun., 9 (2007), p. 648

[26] M.A. Spackman; J.J. McKinnon Cryst. Eng. Commun., 4 (2002), p. 378

[27] S.K. Wolff; D.J. Grimwood; J.J. McKinnon; D. Jayatilaka; M.A. Spackman CrystalExplorer 2.0, University of Western Australia, Perth, Australia, 2007

[28] R.F.W. Bader Atoms in Molecules a Quantum Theory, Oxford University Press, Oxford, UK, 1990

[29] C.F. Matta; R.J. Boyd The Quantum Theory of Atoms in Molecules. From Solid State to DNA and Drug Design, Wiley-VCH, 2007

[30] J.R. Lane; J. Contreras-García; J.P. Piquemal; B.J. Miller; H.G. Kjaergaard J. Chem. Theory Comput., 9 (2013), p. 3263

[31] E.R. Johnson; S. Keinan; P. Mori-Sánchez; J. Contreras-García; A.J. Cohen; W. Yang J. Am. Chem. Soc., 9 (2010), p. 6498

[32] J.M. Dixon; M. Taniguchi; J.S. Lindsey Photochem. Photobiol., 81 (2005), pp. 212-221

[33] Z. Demircioglu; Ç. Albayrak; O. Büyükgüngör J. Mol. Struct., 1065–1066 (2014), p. 210

[34] N. Özdemir; S. Dayan; O. Dayan; M. Dinçer; N.Ö. Kalaycıoglu Mol. Phys., 111 (2013), p. 707

[35] H. Tanak; A.A. Agar; O. Büyükgüngör J. Mol. Struct., 1048 (2013), p. 41

[36] H. Pir; N. Günay; D. Avcı; Y. Atalay Spectrochim. Acta, Part A, 96 (2012), p. 916


Commentaires - Politique


Ces articles pourraient vous intéresser

Synthesis, crystal structure, and density functional theory study of a zinc(II) complex containing terpyridine and pyridine-2,6-dicarboxylic acid ligands: Analysis of the interactions with amoxicillin

Adedibu C. Tella; Joshua A. Obaleye; Margaret D. Olawale; ...

C. R. Chim (2019)