Open access peer-reviewed chapter

Quantum Computational Chemistry: Modeling and Calculation of S-Block Metal Ion Complexes

Written By

Rakesh Kumar and Sangeeta Obrai

Submitted: 14 September 2019 Reviewed: 17 November 2019 Published: 11 March 2020

DOI: 10.5772/intechopen.90531

From the Edited Volume

Density Functional Theory Calculations

Edited by Sergio Ricardo De Lazaro, Luis Henrique Da Silveira Lacerda and Renan Augusto Pontes Ribeiro

Chapter metrics overview

783 Chapter Downloads

View Full Metrics

Abstract

The computational study of some s-block metal nitrophenolate complexes, [Ca(THEEN)(PIC)]+ (1), [Ca(THPEN)(H2O)2]2+ (2), Ba(THPEN)(PIC)2 (3) [Na(THPEN)]22+ (4), [Sr(THPEN)(H2O)2]22+ (5) and [Ba(THPEN)(H2O)2]22+ (6) (where THEEN (N,N,N′,N′-Tetrakis(2-hydroxyethyl)ethylenediamine) and THPEN (N,N,N′,N′-Tetrakis(2-hydroxypropyl)ethylenediamine) are tetrapodal ligands and PIC− is 2,4,6-trinitrophenolate anion), is presented here using density functional theory (DFT) in its hybrid form B3LYP. The geometries of the title complexes are described by the quantum-chemical approach using input coordinates obtained from the previously synthesized and X-ray characterized diffraction data of [Ca(THEEN)(PIC)](PIC), [Ca(THPEN)(H2O)2](PIC)2, Ba(THPEN)(PIC)2, [Na(THPEN)]2(PIC)2, [Sr(THPEN)(H2O)2]2(DNP)4 and [Ba(THPEN)(H2O)2]2(DNP)4 (where DNP is 3,5-dinitrophenolate). Only the primary coordination sphere of complexes (1–6) is optimized in the gaseous phase. Calculations of the energy gaps of frontier orbitals (HOMO-LUMO), 13C-NMR shifts and vibrational bands are carried out using B3LYP/6-31 g + (d,p)/LANL2DZ level of theory. The calculated geometric and spectral parameters reproduced the experimental data with a well agreement.

Keywords

  • DFT
  • s-block metal complexes
  • nitrophenolates
  • tetrapodal ligands

1. Introduction

The alkali and alkaline earth metal cations have inert gas electronic structures and are not expected to show any stereochemical requirements in their complex formation as do transition metal cations. They may be considered spherical even in the complex state. Their complexation is thus treated as recognition of spherical cations by organic ligands [1]. Depending upon the nature of the organic ligand and the anion, the metal ions can be separated as solvated ions, solvent separate, loose and tight ion pairs.

Alkali and alkaline earth metal ions form a large number of solid complexes with podands [2, 3, 4, 5, 6]. The podands are inherently flexible because the two ends of the molecule are not tied simultaneously. Polypodal ligands are acyclic multidentate ligands containing more than three arms. They form an unlimited family of structures which finally give rise to dendrimers. A predominant 1:1 complexation has also been observed in alkali and alkaline earth metal complexes of the tetra- and pentapodands. The stability constants of the terapodands are generally lower than those of the corresponding tripodands because of more severe steric hindrance to complexation [7]. Vögtle and coworkers have indicated that the ligands resembling tetrapodands are capable of forming 1:1 complexes with s-block metal ions [8, 9].

The s-block elements present a usual challenge in the molecular modeling, because the metal-ligand interactions in both cases are principally electrostatic. The types of alkali and alkaline earth metal complexes subjected to molecular modeling can be divided into five categories: crown ethers [10, 11, 12, 13, 14, 15, 16], cryptands [17, 18], spherands [19, 20], podands and other biologically important ligands, such as ionophores and cyclic antibiotics [21, 22, 23, 24] . The present work has been undertaken with the aim to computationally characterize the structure and nature of complexes of s-block metal ions with the tetrapodands THEEN and THPEN. Recently the computational studies of these tetrapodal ligands with Cu(II), Ag(I) and La(III) have been reported. Recently, synthesis, crystal structure and biological properties of [Co(edtp)Cl]·NO3·H2O complex was also determined, where edtp is N,N,N′,N′-Tetrakis(2-hydroxypropyl) ethylenediamine in which Co2+ ion is coordinated by the N,N′,O,O′,O″-pentadentate edtp ligand and a chloride to generate a distorted CoClN2O3 octahedron [25, 26, 27, 28].

Advertisement

2. Computational details

From the last 3 decades density functional theory has been the dominant method for the quantum mechanical simulation of periodic systems. In recent years it has also been adopted by quantum chemists and is now very widely used for the simulation of energy surfaces in molecules. The quantum-chemical calculations (DFT calculations) giving molecular geometries of minimum energies, molecular orbitals (HOMO-LUMO), 13C-NMR and vibrational spectra were performed using the Gaussian 09 [29]. Molecular orbitals were visualized using “Gauss view”. The method used was Becke’s three-parameter hybrid-exchange functional, the nonlocal correlation provided by the Lee, Yang and Parr expression, and the Vosko, Wilk, and Nuair 1980 local correlation functional (III) (B3LYP) [30, 31]. The 6-31 g + (d,p) basis set was used for C, N and O. The LANL2DZ basis set [32] and pseudopotentials of Hay and Wadt were used for Ca, Sr, Ba and Na metal atoms [33, 34]. DFT calculations were performed in the gaseous phase and the input coordinates were obtained from and then compared with crystal structure data of already reported complexes: [Ca(THEEN)(PIC)](PIC),[Ca(THPEN)(H2O)2](PIC)2, Ba(THPEN)(PIC)2, [Na(THPEN)]2(PIC)2, [Sr(THPEN)(H2O)2]2(DNP)4 and [Ba(THPEN)(H2O)2]2(DNP)4 (where DNP is 3,5-dinitrophenolate) [35]. The structural parameters were adjusted until an optimal agreement between calculated and experimental structure obtained throughout the entire range of available structures. HOMO-LUMO analyses and spectroscopic calculations were performed on the optimized geometries of the title complexes (16) using B3LYP/6-31 g + (d,p)/LANL2DZ level of theory.

Advertisement

3. Geometrical structures

3.1 Results and discussion

Complexes (16) were successfully modeled by using the input coordinates of crystal data. Tables 1 and 2 represent comparison of calculated and experimental M-Ligand bond lengths (Å) of complexes (1–3) and (46), respectively. The picrates and dinitrophenolates that are excluded from the primary coordination spheres of metal atoms in the crystallographic determinations, are not optimized in the computed structures. Tables 3 and 4 represent comparison of calculated and experimental torsion angles of ligand in complexes (1–3) and (46) respectively.

Complex 1 (M = Ca)Complex 2 (M = Ca)Complex 3 (M = Ba)
Theo.Exp.Dev.Theo.Exp.Dev.Theo.Exp.Dev.
M-O12.4502.4100.0402.3892.3410.0482.7202.722−0.002
M-O22.4502.3800.0702.4952.498−0.0032.8072.812−0.006
M-O32.4102.4800.0702.7532.7530.000
M-O42.4102.3700.0402.8122.816−0.004
M-O52.6812.687−0.006
M-O63.1273.135−0.008
M-O122.3102.300.012.7282.735−0.007
M-O132.4702.730.26
M-O182.9772.990−0.013
M-O1W2.4402.442−0.002
M-O2W
M-N12.8242.5910.2332.6012.603−0.0023.0383.0420.004
M-N22.7382.6580.083.0263.032−0.006

Table 1.

Comparison of experimental and calculated M-ligand bond lengths (Å) of complexes (1–3).

Complex 4 (M = Na)Complex 5 (M = Sr)Complex 6 (M = Ba)
Theo.Exp.Dev.Theo.Exp.Dev.Theo.Exp.Dev.
M-O12.4122.416−0.0042.6172.628−0.0112.7362.743−0.007
M-O22.3932.396−0.0032.6112.618−0.0072.7632.767−0.004
M-O32.5052.508−0.0032.5062.516−0.0102.7562.762−0.006
M-O42.6282.632−0.0042.6182.626−0.0082.6582.668−0.010
M-O1W2.7012.711−0.0102.8802.884−0.004
M-O2W2.6992.705−0.0062.9882.995−0.007
M-O2WA2.7262.732−0.006
M-N12.8352.842−0.0073.0093.008−0.001
M-N22.8492.857−0.0083.0103.0100.000

Table 2.

Comparison of experimental and calculated M-ligand bond lengths (Å) of complexes (46).

Name of theComplex (1)Complex (2)Complex (3)
AtomsTheo.Exp.Dev.Theo.Exp.Dev.Theo.Exp.Dev.
O1-C1-C2-N1−55.1−55.050.0557.958.00.1−32.3−32.4−0.1
C1-C2-N1-C5165.1165.10.082.882.70.1161.9162.00.1
C2-N1-C5-C6−159.0−159.00.0−165.5−165.50.0−157.6−157.8−0.2
N1-C5-C6-N265.465.30.160.460.50.165.966.10.2
C5-C6-N2-C8−159.1−159.00.1−159.3−159.40.1
C6-N2-C8-C789.189.10.0163.6163.70.1
N2-C8-C7-O360.960.90.0−34.1−34.2−0.1
O2-C3-C4-N136.536.70.255.555.60.146.846.9−0.1
C3-C4-N1-C5−121.6−121.7−0.1−161.0−161.00.0−132.6−132.7−0.1
C4-N1-C5-C677.277.20.072.672.60.083.783.80.1
C5-C6-N2-C1079.479.40.079.079.10.1
C6-N2-C10-C9−151.8−151.80.0−134.9−134.90.0
N2-C10-C9-O461.361.40.152.152.20.1

Table 3.

Comparison of calculated and experimental torsion angles (°) of ligand in the complexes (1–3).

Name of theComplex (4)Complex (5)Complex (6)
AtomsTheo.Exp.Dev.Theo.Exp.Dev.Theo.Exp.Dev.
Complex−55.6−55.7−0.132.332.30.0−43.4−43.50.1
O1-C1-C2-N1152.7152.60.1−143.9−143.90.0162.3162.40.1
C1-C2-N1-C5−80.2−80.20.086.186.10.0−92.7−92.80.1
C2-N1-C5-C6−64.0−63.9−0.143.643.60.0−50.1−50.10.0
N1-C5-C6-N2165.5165.40.1−146.4−146.40.0156.5156.70.2
C5-C6-N2-C8−85.7−85.70.0156155.90.1−89.6−89.90.3
C6-N2-C8-C7−60.3−60.4−0.1−44.7−44.70.0−30.4−30.40.0
N2-C8-C7-O3−59.7−59.8−0.1−17.4−17.40.047.848.10.3
O2-C3-C4-N1−91.0−91.1−0.1112.7112.70.0−135.1−135.20.1
C3-C4-N1-C5156.0156.00.0−159.1−159.20.1151.9151.90.0
C4-N1-C5-C6−72.4−72.5−0.197.497.40.0−92.1−92.40.3
C5-C6-N2-C10156.6156.60.0−151.8−151.80.0143.4143.90.5
C6-N2-C10-C9−59.2−59.3−0.141.641.70.15.14.90.3
N2-C10-C9-O4−55.6−55.7−0.132.332.30.0−43.4−43.50.1

Table 4.

Comparison of experimental and calculated and torsion angles (°) of ligand in the complexes (4–6).

3.1.1 [Ca(THEEN)(PIC)]+ (1)

The coordination number of Ca(II) ion is eight with distorted cube geometry in the optimized geometry of cationic complex (1) (Figure 1). THEEN is interacting with Ca(II) ion through all the six potential donor atoms. The seventh and eighth coordination sites of Ca(II) are occupied by picrate anion through phenolic oxygen and one of the o-nitro oxygen. The observed and calculated positions of calcium and donor atoms are in agreement. A comparison of bond lengths and bond angles provided a maximum tolerance of 0.27 Å and 16.33°Å, respectively (Table 1 and Table S1; Figure S1a,b) of computed title complex (1) and crystallographically determined complex [Ca(THEEN)(PIC)](PIC). The torsion angles of the ligand THEEN in theoretically determined and experimental complex are also in well consistency with each other (Table 3). The HOMO-LUMO analysis has indicated that there is maximum distribution of HOMO over the carbon atoms and to a lesser extent on the oxygen of the coordinated picrate. LUMO is mainly distributed or concentrated on all the atoms of coordinated picrate except p-nitro group. Neither HOMO nor LUMO distribution is found either on calcium ion or on any ligand atom. The HOMO-LUMO gap (ΔE) is found to be 0.852 eV (Figure S1c).

Figure 1.

(a) Optimized geometric structure of [Ca(THEEN)(PIC)]+ (1) (b) distorted cube geometry.

3.1.2 [Ca(THPEN)(H2O)2]2+ (2)

Ca(II) is eight coordinated in monomeric cationic complex with a distorted square-antiprismatic geometry in complex (2) as is observed in crystal structure of [Ca(THPEN)(H2O)2] (PIC)2. H2O (Figure 2). THPEN is acting as hexadentate ligand towards the metal ion. Two remaining sites around the metal ion are occupied by two water molecules. The observed and calculated positions of calcium and donor atoms are in agreement. All strain energy minimized structures reproduced the observed X-ray structures with almost no deviation in M-L bond length and torsion angle of the ligand (Table 3). The maximum tolerance of L-M-L bond angle is 15.8° (Tables 1 and Table S1, Figure S2a,b). Complex 2 is displaying the main distribution of HOMO as well as LUMO only on water molecules. There is slight distribution of HOMO and LUMO on hydroxyl oxygen of THPEN ligand also with HOMO-LUMO gap (HLG) of 0.419 eV (Figure S2c).

Figure 2.

(a) Optimized geometric structure of [Ca(THPEN)(H2O)2]2+ (2) (b) distorted square-antiprismatic geometry.

3.1.3 Ba(THPEN)(PIC)2 (3)

Ba(II) is ten-coordinate in its monomeric complex (3) (Figure 3). The ligand THPEN coordinates through both its nitrogen atoms and all the four hydroxyl oxygen. The two picrate ions are also directly coordinated to Ba2+ ion through the phenolic oxygen and one oxygen of the o-nitro group. Since both the picrates are directly interacting with the cation, the complex is termed as tight ion-paired complex as is found in the crystallographically determined complex. The observed and calculated positions of barium and donor atoms are in agreement. There is negligible deviation found in bond length, bond angle (0.11°) and torsion angle of ligand THPEN for calculated title complex (3) and experimental Ba(THPEN)(PIC)2 (Tables 1 and 3 and S1, Figure S3a,b). An analysis of HOMO-LUMO has illustrated that HOMO is mainly distributed on the ligand THPEN with a small distribution on coordinated picrates. In contrast to this LUMO is mainly distributed on coordinated picrates with no distribution over the ligand or metal atom. The complex is showing HLG of 0.118 eV (Figure S3c).

Figure 3.

(a) Optimized geometric structure of Ba(THPEN)(PIC)2 (3) (b) bicapped square-antiprismatic geometry.

3.1.4 [Na(THPEN)]22+ (4)

Figure 4a shows the optimized structure of the cationic complex of sodium (4). This structure is centrosymmetric dimer. The coordination number around each Na+ ion is seven. Each THPEN ligand coordinates to Na(I) ion in a heptadenate fashion (Figure 4). In other words, one hydroxyl group of the THPEN ligand acts as a bridge between two Na+ ions. The geometry of the complex is distorted monocapped octahedron. The Na…Na non bridging distance is 3.429 Å. Almost no deviation has been observed in bond length (M-L), torsion angle of ligand THPEN and bond angles (L-M-L) of computed title complex (4) and crystallographically determined complex [Na(THPEN)]2(PIC)2 (Tables 2 and 4 and Table S1, Figures S4a,b). The HOMO-LUMO study has revealed that in the title complex (4) HOMO is mainly concentrated on bridged hydroxyl oxygens and sodium metal but to a smaller extent on the other coordinated hydroxyl oxygens and amine nitrogens. LUMO is mainly concentrated on bridged hydroxyl oxygens and to a smaller extent on other coordinated oxygen atoms. The HOMO-LUMO gap is 0.261 eV (Figure S4c).

Figure 4.

(a) Optimized geometric structure of [Na(THPEN)]22+(4) (b) distorted monocapped octahedron geometry.

3.1.5 [Sr(THPEN)(H2O)2]22+ (5)

Sr(II) is nine-coordinated in its dimeric complex (5) (Figure 5). This cationic complex is having tricapped trigonal-prismatic geometry as is found in crystallographically determined complex [Sr(THPEN)(H2O)2]2(DNP)4. Each Sr2+ ion in the complex is coordinated by six potential donor sites of the ligand THPEN and four water molecules. Out of four water molecules, two are bridging and the third one is non-bridging. Sr…Sr non bridging distance is 4.346 Å indicating the existence of van der Waals contact. All strain energy minimized structures reproduced the observed X-ray structures to a maximum tolerance of 9.07° bond angle (Table 2 and Table S1, Figure S5a,b). Almost no deviation of M-L bond length and torsion angle of ligand THPEN has been observed for the title complex (5) and crystallographically determined complex (Tables 2 and 4). The HOMO-LUMO analysis has shown that the complex (5) is having maximum distribution of HOMO and LUMO on bridged coordinated water molecules and there is HLG of 0.0225 eV (Figure S5c).

Figure 5.

(a) Optimized geometric structure of [Sr(THPEN)(H2O)2]2(DNP)4 (b) Trigonal- prismatic geometry.

3.1.6 [Ba(THPEN)(H2O)2]22+ (6)

Ba(II) is ten-coordinate in the cationic title complex (6) as is found in the crystallographic determined complex [Ba(THPEN)(H2O)2](DNP)4. The geometry around Ba(II) is bicapped cubic. Each Ba(II) in dimer is interacting with THPEN through all its six potential donor sites and four water molecules. The latter are bridging in nature. The existence of van der Waals contact between non-bridging Ba…Ba is indicated by their larger distance (4.196 Å). The observed and calculated positions of the metal and donor atoms are in agreement as almost negligible deviation of M-L bond length and torsion angle of THPEN. A deviation of 0.14° of bond angle L-M-L has been observed for complex (6) [Ba(THPEN)(H2O)2](DNP)4 (Tables 2 and 4 and Table S1; Figures S6a,b and S7). The region of distribution of HOMO and LUMO is only over two of the bridged water molecules in complex (6) with a very small HOMO-LUMO gap (ΔE = 0.0375 eV) indicating the soft nature of complex (Figure 6).

Figure 6.

(a) Optimized geometric structure of [Ba(THPEN)(H2O)2]22+ (6) (b) Bicapped cubic geometry.

3.2 HOMO-LUMO analyses

Smaller is the HUMO-LUMO gap (HLG) softer is the complex [36, 37]. The frontier orbitals HOMO and LUMO are very important parameters for chemical reaction and take part in chemical stability [38, 39, 40]. It is predicted from the HOMO-LUMO gaps that the title complexes are soft as is obvious from their smaller HUMO-LUMO energy gaps (HLG) relative to the similar reported complexes of copper, silver and lanthanoid [25, 26, 27] (Table 5). It has been observed in the present computational study that the dinitrophenolate complexes are softer than trinitrophenolate and among the latter, [Ca(THEEN)(PIC)]+ (1) is displaying least softness. It is pertinent to mention here that the complex (1) is having THEEN ligand whereas THPEN is ligand in rest of the complexes (26).

HLG (eV)Reported complexesHLG (eV)Reported complexesHLG (eV)
[Ca(THEEN)(PIC)]+ (1)0.852[Cu(THEEN)(H2O)](PIC)23.537[La(THEEN)(PIC) (H2O)2] (PIC)2.2H2O3.428
[Ca(THPEN)(H2O)2]2+ (2)0.419[Cu(THPEN)](PIC)2.C3H8O3.467[La(TEAH3)(H2O)2] (PIC)33.673
Ba(THPEN)(PIC)2 (3)0.118[Cu(TEAH3)(PIC)](PIC).H2O3.619
[Na(THPEN)]22+ (4)0.261[Ag (THEEN)]2 (PIC)22.530
[Sr(THPEN)(H2O)2]22+ (5)0.0225[Ag (THPEN)]2 (PIC)22.640
[Ba(THPEN)(H2O)2]22+ (6)0.0375[Ag(TEAH3)2] (PIC)1.061

Table 5.

Comparison of HOMO-LUMO energy gaps of complexes (16) with earlier reported complexes [25, 26, 27].

3.3 Spectral data

Nuclear magnetic resonance spectra (NMR) and infrared (IR) spectroscopy can be useful for studying the coordination of various ligating sites. The 13C-NMR spectra were predicted for complexes (1–6) using DFT/B3LYP/6-31G** method and the spectral data was compared with experimental data reported in literature [35]. The computed NMR spectral data is fairly in agreement with experimental data (Tables 6 and 7). Small deviations are due to the fact that H-bonding interactions or any type of lattice interactions are not modeled in theoretically computed structures. The terminal methyl groups in theoretically predicted complexes (25) are displaying quiet high upfield shifts relative to the experimentally obtained due to more free movements in the gaseous phase than in solution or solid phase.

Assignments (δ)[Ca(THEEN)(PIC)]+ (1)Ba(THPEN)(PIC)2 (3)
Theo.Exp.Theo.Exp.
▬CH3**4.2219.61
▬CH34.4919.89
▬NCH252.3155.0351.7355.34
▬NCH243.3555.20
▬OCH2, ▬OCH59.0357.4750.7563.00
▬ArCH113.63124.33117.85124.36
p-ArCN134.75123.60122.06123.68
o-ArCN140.68140.78140.98140.78
▬ArCO161.57160.20150.21160.20

Table 6.

Comparison of calculated and experimental 13C-NMR spectral data for complexes (1 and 3).

* Group absent.


Assignment (δ)[Ca(THPEN)(H2O)2]2+ (2)[Na(THPEN)]22+ (4)[Sr(THPEN)(H2O)2]22+ (5)[Ba(THPEN)(H2O)2]22+ (6)
Theo.Exp.Theo.Exp.Theo.Exp.Theo.Exp.
▬CH311.3518.794.0120.303.7718.3127.0218.35
▬CH311.3819.025.6320.363.8518.5927.6918.63
▬NCH239.6050.9943.0550.4243.6048.3161.6959.80
▬NCH239.6151.7644.0552.9057.8159.6861.6960.59
▬OCH**50.7355.7341.8560.6260.3760.66
▬OCH69.2461.7457.9055.9551.9161.0860.3862.53

Table 7.

Comparison of calculated and experimental 13C-NMR spectral data for complexes (2, 4–6).

* Group absent.


The computed IR spectral peaks that appear in the range of 1370–600 cm−1 are fairly in agreement with the experimental data (Tables 8 and 9). The absorption peaks due to the presence of hydroxyl groups were observed only for complexes (1) and (6) in the computed IR spectra. It is pertinent to mention here that both of these complexes possess the picrate anion in their coordination sphere. The less extent of H-bonding is reported in the crystallographic description of these complexes [35]. The theoretical absorption band appears at 3223 and 3235 cm−1 for both complexes (1) and (6) whereas the experimental band is reported at 3300 cm−1 for both of them [35].

Assignments (cm−1)[Ca(THEEN)(PIC)]+ (1)[Ca(THPEN)(H2O)2]2+ (2)Ba(THPEN)(PIC)2 (3)
Theo.Exp.Theo.Exp.Theo.Exp.
ν (NO2)1363.561360 m1320.811360 vs1384.201370
δ (〓CH)696.80700 m784.16790 m800800

Table 8.

Comparison of calculated and experimental IR spectral data for complexes (1–3).

Assignments (cm−1)[Na(THPEN)]22+ (4)[Sr(THPEN)(H2O)2]22+ (5)[Ba(THPEN)(H2O)2]22+ (6)
Theo.Exp.Theo.Exp.Theo.Exp.
ν (NO2)13661371.2113301329.2413401329.43
δ (〓CH)790784.54760762.52790787.17

Table 9.

Comparison of calculated and experimental IR spectral data for complexes (4–6).

Advertisement

4. Conclusions

The coordination number of the title s-block complexes is varying from 7 to 10 in the present work. As the size of the metal increases, coordination number of central metal ion also increases. Out of the six complexes presented, three are monomeric (1–3) and three are dimeric (4–6). The longer distances between two M…M distances in the dinuclear complexes indicate the existence of van der Waals contacts between the two s-block metal ions. All the title complexes are cationic except Ba(THPEN)(PIC)2 (3). The latter is tight ion-paired complex. All strain energy minimized structures obtained using quantum-chemical approach reproduced the observed X-ray structures with geometric parameters in well agreement. HOMO-LUMO studies suggest the softness of the title s-block complexes relative to the similar already reported copper, silver and lanthanoid complexes. The theoretical spectral data (13C-NMR and IR) computed using DFT and experimental data is fairly in agreement with each other. The accuracy of the results predicts that the DFT studies performed using B3LYP/6-31 g + (d,p)/LANL2DZ level of theory is the appropriate quantum-chemical method for reproducing the experimental results for the title s-block complexes. This quantum-chemical approach has potential for molecular modeling of other s-block complexes and exploring their chemistry.

Small deviations in geometric as well as spectral parameters may be attributed to the lack of H-bonding and packing interactions within lattice which were not modeled during the computational study of the entitled s-block complexes. Moreover, the quantum-chemical approach of DFT studies has been carried out in the gaseous phase whereas the already reported experimental crystal and IR spectral data is in the solid phase while 13C-NMR spectral data is in the solution phase.

Advertisement

Appendix

Figure S1.

Plot showing the deviations of theoretical and experimental (a) bond lengths (Å) and (b) bond angles (°) for the complex (1) and (c) (HOMO-LUMO) of the complex (1) with energy gap.

Figure S2.

Plot showing the deviations of theoretical and experimental (a) bond lengths (Å) and (b) bond angles (°) for the complex (2) (c) (HOMO-LUMO) of the complex (2) with energy gap.

Figure S3.

Plot showing the deviations of theoretical and experimental (a) bond lengths (Å) and (b) bond angles (°) for the complex (3) (c) (HOMO-LUMO) of the complex (3) with energy gap.

Figure S4.

Plot showing the deviations of theoretical and experimental (a) bond lengths (Å) and (b) bond angles (°) for the complex (4) (c) (HOMO-LUMO) of the complex (4) with energy gap.

Figure S5.

Plot showing the deviations of theoretical and experimental (a) bond lengths (Å) and (b) bond angles (°) for the complex (5) (c) (HOMO-LUMO) of the complex (5) with energy gap.

Figure S6.

Plot showing the deviations of theoretical and experimental (a) bond lengths (Å) and (b) bond angles (°) for the complex (6) (c) (HOMO-LUMO) of the complex (6) with energy gap.

Figure S7.

R = H, THEEN; R = CH3, THPEN.

Bond distances (Å)TheoreticalExperimentalDev.Bond angles (°)TheoreticalExperimentalDev.
Complex (1)
Ca-N12.6002.6000.000N1-Ca-N1A71.5471.600.06
Ca-O12.3872.3880.001N1-Ca-O2A81.4165.5215.89
Ca-O22.4952.4960.001N1-Ca-O1W119.08119.090.01
Ca-O1W2.4392.4390.000N1-Ca-O1WA143.89143.890.00
O1-Ca-O2102.93102.950.02
O1-Ca-N167.8067.790.01
O1-Ca-O1W79.8279.850.03
O1-Ca-N1A132.22132.220.00
O1-Ca-O2A84.3784.360.01
O1-Ca-O1A159.11159.100.01
O1-Ca-O1WA83.5283.540.02
O2-Ca-N165.5065.520.02
O2-Ca-O1W74.1174.100.01
O2-Ca-N1A81.4281.400.02
O2-Ca-O2A139.33139.890.56
O2A-Ca-O1W145.88145.900.02
Complex (2)
Ca-O12.452.410.04O1-Ca-O294.43102.087.65
Ca-O22.452.380.07O1-Ca-O3172.00175.213.21
Ca-O32.412.480.07O1-Ca-O485.5674.3511.21
Ca-O42.412.370.04O1-Ca-O12107.80104.183.62
Ca-O122.312.300.01O1-Ca-O1361.9164.852.94
Ca-O132.472.730.26O1-Ca-N165.9368.502.57
Ca-N12.742.590.15O1-Ca-N2104.33115.4911.16
Ca-N22.822.650.17O2-Ca-O390.9979.2711.72
O2-Ca-O4180.00168.6511.35
O2-Ca-O12108.06107.200.86
O2-Ca-O1361.3869.928.54
O2-Ca-N164.2068.434.23
O2-Ca-N2114.39105.728.67
O3-Ca-O489.01105.2416.23
O3-Ca-O1275.9571.064.89
O3-Ca-O13126.03110.6915.34
O3-Ca-N1111.54116.184.64
O3-Ca-N267.9768.260.29
O4-Ca-O1271.9384.2012.27
O4-Ca-O13118.62116.911.71
O4-Ca-N1115.80100.3815.42
O4-Ca-N265.6167.361.75
O12-Ca-O1371.3262.199.13
O12-Ca-N1168.56169.550.99
O12-Ca-N2123.47120.542.93
O13-Ca-N197.26107.4610.2
O13-Ca-N2163.73175.6011.87
N1-Ca-N267.9569.921.97
Complex (3)
Ba-O12.7202.7200.000O1-Ba-O286.9587.000.05
Ba-O22.8072.8010.006O1-Ba-O3172.87173.000.13
Ba-O32.7532.7500.003O1-Ba-O488.5588.500.05
Ba-O42.8122.8120.000O1-Ba-O672.9072.900.00
Ba-O52.6812.6830.002O1-Ba-O1261.7961.800.01
Ba-O63.1273.1280.001O1-Ba-O18109.63109.700.07
Ba-O122.7282.7290.001O1-Ba-N158.1258.100.02
Ba-O182.9772.9790.002O1-Ba-N2114.85114.800.05
Ba-N13.0383.0380.000O2-Ba-O4129.37129.400.03
Ba-N23.0263.0280.002O2-Ba-O664.7964.800.01
O2-Ba-O18161.62160.600.02
O2-Ba-N279.3679.400.03
O2-Ba-N157.0357.120.05
O3-Ba-O292.0892.100.02
O3-Ba-O486.6186.700.09
O3-Ba-O6113.02113.100.08
O3-Ba-O1872.7072.800.10
O3-Ba-N1115.63115.600.03
O3-Ba-N258.0658.200.14
O4-Ba-O6156.82156.800.02
O4-Ba-O1863.1363.200.07
O4-Ba-N178.2078.200.00
O4-Ba-N257.4357.400.03
O5-Ba-O1123.12123.200.08
O5-Ba-O285.9485.900.04
O5-Ba-O363.8063.900.10
O5-Ba-O4135.92135.900.02
O5-Ba-O653.2653.300.03
O5-Ba-O1284.084.000.00
O5-Ba-O1876.6776.900.13
O5-Ba-N1142.90142.900.00
O5-Ba-N2118.97118.900.07
O12-Ba-O2133.15133.200.05
O12-Ba-O3123.02123.000.02
O12-Ba-O486.4586.400.05
O12-Ba-O672.6972.700.01
O12-Ba-O1854.0554.000.05
O12-Ba-N1117.95117.960.01
O12-Ba-N2143.78143.800.02
O18-Ba-O6109.54109.600.06
O18-Ba-N1140.27140.400.13
O18-Ba-N2101.49101.500.01
N2-Ba-N161.2761.300.03
N2-Ba-O6143.21143.200.01
N1-Ba-O6102.39102.400.01
Complex (4)
Na-N12.5522.5540.002O1-Na-O299.5699.500.06
Na-N22.5652.5660.001O1-Na-O386.5986.600.01
Na-O12.4122.4120.000O1-Na-O4174.34174.300.04
Na-O22.3932.3930.000O1-Na-O4A90.1890230.05
Na-O32.5052.5050.000O1-Na-N169.8269.820.00
Na-O42.6282.6290.001O1-Na-N2109.13109.100.03
Na-O4A2.4422.4430.001O2-Na-O3164.10164.130.03
Na-NaA3.4293.4300.001O2-Na-O478.2678.310.05
O2-Na-O4A88.5688.570.01
O2-Na-N170.6970.700.01
O2-Na-N2121.68121.710.03
O3-Na-O496.8696.870.01
O3-Na-N1125.19125.170.02
O3-Na-N268.9368.910.02
O3-Na-O4A76.6976.720.03
O4-Na-O4A94.9695.000.03
O4-Na-N168.2568.250.00
O4-Na-N268.2568.250.00
N1-Na-N273.3873.370.01
N1-Na-O4A147.48147.510.03
N2-Na-O4A138.85138.870.02
Complex(5)
Sr-O12.6172.6180.001O3-Sr-O2151.28158.707.42
Sr-O22.6112.6100.001O3-Sr-O173.2477.203.96
Sr-O32.5062.5050.001O2-Sr-O1113.98114.300.32
Sr-O42.6182.6180.000O3-Sr-O487.9188.200.29
Sr-O1W2.7012.7020.001O2-Sr-O477.8973.504.39
Sr-O2W2.6992.6990.000O1-Sr-O4158.36151.606.76
Sr-O2WA2.7262.6990.000O3-Sr-N1102.34113.6011.26
Sr-N12.8352.8350.000O4-Sr-N1113.34102.5010.84
Sr-N22.8492.8490.000O3-Sr-N261.3864.403.02
O2-Sr-N289.9297.107.18
O1-Sr-N296.8090.26.60
O4-Sr-N264.1361.602.53
O3-Sr-O2W106.07106.400.33
O2-Sr-O2W70.5877.106.52
O1-Sr-O2W129.35140.5011.15
O4-Sr-O2W70.8067.902.90
O2-Sr-O1W79.5265.4014.12
O1-Sr-O1W65.1179.8014.69
O4-Sr-O1W136.42127.409.02
O3’-Sr-O1W126.14127.901.76
O2W-Sr-O1W70.7571.000.25
O2W-Sr-N2128.11128.400.29
O1W-Sr-N2152.06152.400.34
N1-Sr-N264.7765.000.23
O2W-Sr-N1139.03139.300.27
O1W-Sr-N187.5287.800.28
Complex (6)
Ba-O12.7362.7380.002O4-Ba-O1162.70162.900.20
Ba-O22.7632.7610.002O4-Ba-O392.4592.700.25
Ba-O32.7562.7570.001O1-Ba-O379.6980.000.31
Ba-O42.6582.6570.001O4-Ba-O284.2584.600.35
Ba-O1W2.8802.8800.000O1-Ba-O294.2294.500.28
Ba-O2W2.9882.990-0.002O3-Ba-O2148.19148.400.21
Ba-N13.0093.0020.007O4-Ba-N1107.19107.500.31
Ba-N23.0103.0040.006O1-Ba-N158.1958.400.21
O3-Ba-N191.9092.000.01
O2-Ba-N159.4359.600.17
O4-Ba-N259.9561.101.15
O1-Ba-N2103.03103.200.17
O3-Ba-N258.6858.800.12
O2-Ba-N293.293.400.20
N1-Ba-N261.8462.000.16
O4-Ba-O1W126.17126.500.33
O1-Ba-O1W66.2966.500.23
O3-Ba-O1W72.5372.800.27
O2-Ba-O1W133.50133.800.30
O1-Ba-O2W73.4373.700.27
O3-Ba-O2W131.18131.500.32
O2-Ba-O2W74.7475.000.26
O2W-Ba-O1WA59.3659.500.14
N1-Ba-O1WA137.03137.200.17
N1-Ba-O1W124.19124.190.00
N1-Ba-O2W106.68106.900.22
N2-Ba-O1WA123.73124.000.27
N2-Ba-O1W131.21131.500.29
N2-Ba-O2W166.91167.000.10
N2-Ba-O2WA100.84100.870.03

Table S1.

Comparison of selected experimental and calculated geometric parameters bond lengths (Å) and bond angles (°) for complexes (1–6).

References

  1. 1. Lehn JM. Structure and bonding. Alkali Metal Complexes with Organic Ligands. In: Dunitz JD, Hemmirch P, Ibers JA, Jorgensen CK, Neilands JB, Reinan D, Williams RJP, editors. Vol. 16, 7. New York: Springer-Verlag; 1973
  2. 2. Tümmler B, Maas G, Weber E, Wehner W, Vögtle F. Noncyclic crown-type polyethers, pyridinophane cryptands, and their alkali metal ion complexes: Synthesis, complex stability, and kinetics. Journal of the American Chemical Society. 1977;99:4688
  3. 3. Cram DJ, Cram JM. Design of complexes between synthetic hosts and organic guests. Accounts of Chemical Research. 1978;11:8
  4. 4. Cram DJ. Preorganization - from solvents to spherands. Angewandte Chemie (International Ed. in English). 1986;25:1039
  5. 5. Poonia NS, Bajaj AV. Coordination chemistry of alkali and alkaline earth cations. Chemical Reviews. 1979;79:145
  6. 6. Bajaj AV, Poonia NS. Comprehensive coordination chemistry of alkali and alkaline earth cations with macrocyclic multidentates: Latest position. Coordination Chemistry Reviews. 1988;87:55
  7. 7. Vögtle F, Müller WM, Buhleier E, Weber W. Synthesis and selectivity of novel four armed noncyclic neutral ligands. Chemische Berichte. 1979;112:899
  8. 8. Vögtle F, Weber E. Multidentate acyclic neutral ligands and their complexation. Angewandte Chemie (International Ed. in English). 1979;18:753
  9. 9. Vögtle F, Weber E. Crystalline 1: 1 alkali metal complexes of noncyclic neutral ligands. Tetrahedron Letters. 1975;16:2415
  10. 10. Hay BP, Rustad JR. Structural criteria for the rational design of selective ligands: Extension of the MM3 force field to aliphatic ether complexes of the alkali and alkaline earth cations. Journal of the American Chemical Society. 1994;116:6316
  11. 11. Wipff G, Weiner P, Kollman P. A molecular-mechanics study of 18-crown-6 and its alkali complexes: An analysis of structural flexibility, ligand specificity, and the macrocyclic effect. Journal of the American Chemical Society. 1982;104:3249
  12. 12. Grootenhuis PDJ, Kollman PA. Molecular mechanics and dynamics studies of crown ether-cation interactions: Free energy calculations on the cation selectivity of dibenzo-18-crown-6 and dibenzo-30-crown-10. Journal of the American Chemical Society. 1989;111:2152
  13. 13. Pretsch E, Badertscher M, Welti M, Maruizumi T, Morf WE, Simon W. Design features of ionophores for ion selective electrodes. Pure and Applied Chemistry. 1988;60:567
  14. 14. Badertscher M, Welti M, Pretsch E, Maruizumi T, Ha T–K, Simon W. Combined application of pair potentials and the MM2 force field for the modeling of ionophores. Journal of Computational Chemistry. 1990;11:819
  15. 15. Varnek AA, Glebov AS, Petrukhin OM, Kolycheva NV, Ozerov RP. Selectivity of reactions of 18-crown-6 complexation with alkali-metals. Koordinatsionnaya Khimiya (Russian Journal of Coordination Chemistry). 1989;15:600
  16. 16. Burns JH, Kessler RM. Structural and molecular mechanics studies of bis(dibutylphosphato)aquastrontium-18crown-6 and analogous alkaline-earth-metal complexes. Inorganic Chemistry. 1987;26:1370
  17. 17. Wipff G, Kollman P. Molecular mechanical calculations on a macrocyclic receptor: The 222 cryptand and its alkali complexes. Nouveau Journal de Chimie. 1985;9:457
  18. 18. Damu KV, Hancock RD, wade PW, Boeyens JCA, Billing DG, Dobson SM. Control of metal-ion size-based selectivity through the structure of the oxygen-donor pendant groups on lariat ethers. A crystallographic and thermodynamic study. Journal of the Chemical Society, Dalton Transactions. 1991:293
  19. 19. Kollman P, Wipff G, Singh UC. Molecular mechanical studies of inclusion of alkali cations into anisole spherands. Journal of the American Chemical Society. 1985;107:2212
  20. 20. Maye PV, Venanzi CA. Host-guest preorganization and complementarity: A molecular mechanics and molecular dynamics study of cation complexes of a cyclic urea-anisole spherand. Journal of Computational Chemistry. 1991;12:994
  21. 21. Lifson S, Felder CS, Shanzer A. Computational design of polylactone macrocyclic ionophores. Journal of the American Chemical Society. 1983;105:3866
  22. 22. Lifson S, Felder CS, Shanzer A. Design, synthesis, and conformation of ion carriers. Bopolymers. 1983;22:409
  23. 23. Lifson S, Felder CS, Shanzer A. Internal and external alkali ion complexes of enniatin B: An empirical force field analysis. Biochemistry. 1984;23:2577
  24. 24. Lifson S, Felder CS, Shanzer A. Enniatin B and valinomycin as ion carriers: An empirical force field analysis. Journal of Biomolecular Structure & Dynamics. 1984;2:641
  25. 25. Kumar R, Obrai S, Kaur A, Hundal G, Meehnian H, Jana AK. Synthesis, crystal structure determination and antimicrobial studies of copper(II) picrate complexes with N, N, N′′′-tetrakis(2-hydroxyethyl/propyl) ethylenediamine and tris (2-hydroxyethyl)amine. Polyhedron. 2013;56:55
  26. 26. Kumar R, Obrai S, Kaur A, Hundal MS, Meehnian H, Jana AK. Synthesis, crystal structure investigation, DFT analyses and antimicrobial studies of silver(I) complexes with N, N, N′′′-tetrakis (2-hydroxyethyl/propyl) ethylenediamine and tris(2-hydroxyethyl)amine. New Journal of Chemistry. 2014;38:1186
  27. 27. Kumar R, Obrai S, Jassal AK, Hundal MS. Supramolecular architectures of N,N,N′′,N′′-tetrakis-(2-hydroxyethyl)ethylenediamine and tris(2-hydroxyethyl)amine with La(III) picrate. RSC Advances. 2014;4:59248
  28. 28. Akhtar MN, Harrison WTA, Shahid M, Khan IU, Ejaz JI. Synthesis, crystal structure and biological activity of a cobalt(II) complex of N,N,N′′-tetrakis(2-hydroxypropyl) ethylenediamine. Transition Metal Chemistry. 2016;41:325-330
  29. 29. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Gaussian 03, (Revision B.04). Wallingford, Connecticut: Gaussian, Inc.; 2004
  30. 30. Becke AD. Density-functional thermochemistry. III. The role of exact exchange. The Journal of Chemical Physics. 1993;98:5648
  31. 31. Lee C, Yang W, Parr RG. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Physical Review B. 1988;37:785
  32. 32. Patersson GA, Al-Laham MA. A complete basis set model chemistry. II. Open-shell systems and the total energies of the first-row atoms. The Journal of Chemical Physics. 1991;94:6081
  33. 33. Hay PJ, Wadt WR. Ab initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg. The Journal of Chemical Physics. 1985;82:270
  34. 34. Wadt WR, Hay PJ. Ab initio effective core potentials for molecular calculations. Potentials for main group elements Na to Bi. The Journal of Chemical Physics. 1985;82:284
  35. 35. Hundal G, Hundal MS, Obrai S, Poonia NS, Kumar S. Metal complexes of tetrapodal ligands: synthesis, spectroscopic and thermal studies, and X-ray crystal structure studies of Na(I), Ca(II), Sr(II), and Ba(II) complexes of tetrapodal ligands N,N,N′,N′-tetrakis(2-hydroxypropyl)ethylenediamine and N,N,N′,N′-tetrakis(2-hydroxyethyl)ethylenediamine. Inorganic Chemistry. 2002;41:2077
  36. 36. Mariappan G, Sundaraganesan N. Spectral and structural studies of the anti-cancer drug Flutamide by density functional theoretical method. Spectrochimica Acta Part A. 2014;117:604
  37. 37. Fleming I. Frontier Orbitals and Organic Chemical Reactions. New York: John Wiley and Sons; 1976
  38. 38. Balachandran V, Lalitha S, Rajeswari S, Rastogi VK. Theoretical investigations on the molecular structure, vibrational spectra, thermodynamics, HOMO-LUMO, NBO analyses and paramagnetic susceptibility properties of p-(p-hydroxyphenoxy)benzoic acid. Spectrochimica Acta Part A. 2014;121:575
  39. 39. Raif K, Sinem O, Azizoglu A. Experimental and computational study on [2,6-bis(3,5-dimethyl-N-pyrazolyl)pyridine]-(dithiocyanato)mercury(II). Polyhedron. 2007;26:5069
  40. 40. Uesugi Y, Mizuno M, Shimojima A, Takahashi H. Resonance Raman and ab initio MO calculation studies of the structures and vibrational assignments of the T1 state and the anion radical of coumarin and its isotopically substituted analogues. The Journal of Physical Chemistry. 1997;101:268

Written By

Rakesh Kumar and Sangeeta Obrai

Submitted: 14 September 2019 Reviewed: 17 November 2019 Published: 11 March 2020