Investigation of Structural and Electronic Properties of [Tris(Benzene-1,2-Dithiolato)M]3- (M = V, Cr, Mn, Fe and Co) Complexes: A Spectroscopic and Density Functional Theoretical Study


In this study, the first raw transition metals from V to Co complexes with benzene-1,2-dithiolate (L2-) ligand have been studied theoretically to elucidate the geometry, electronic structure and spectroscopic properties of the complexes. Density Functional Theory (DFT) and Time-Dependent Density Functional Theory (TD-DFT) methods have been used. The ground state geometries, binding energies, spectral properties (UV-vis), frontier molecular orbitals (FMOs) analysis, charge analysis and natural bond orbital (NBO) have been investigated. The geometrical parameters are in good agreement with the available experimental data. The metal-ligand binding energies are 1 order of magnitude larger than the physisorption energy of a benzene-1, 2-dthiolate molecule on a metallic surface. The electronic structures of the first raw transition metal series from V to Co have been elucidated by UV-vis spectroscopic using DFT calculations. In accordance with experiment the calculated electronic spectra of these tris complexes show bands at 522, 565, 559, 546 and 863 nm for V3+, Cr3+, Mn3+, Fe3+ and Co3+ respectively which are mainly attributed to ligand to metal charge transfer (LMCT) transitions. The electronic properties analysis shows that the highest occupied molecular orbital (HOMO) is mainly centered on metal coordinated sulfur atoms whereas the lowest unoccupied molecular orbital (LUMO) is mainly located on the metal surface. From calculation of intramolecular interactions and electron delocalization by natural bond orbital (NBO) analysis, the stability of the complexes was estimated. The NBO results showed significant charge transfer from sulfur to central metal ions in the complexes, as well as to the benzene. The calculated charges on metal ions are also reported at various charge schemes. The calculations show encouraging agreement with the available experimental data.

Share and Cite:

Matin, M. and Rahman, M. (2019) Investigation of Structural and Electronic Properties of [Tris(Benzene-1,2-Dithiolato)M]3- (M = V, Cr, Mn, Fe and Co) Complexes: A Spectroscopic and Density Functional Theoretical Study. Advances in Chemical Engineering and Science, 9, 317-332. doi: 10.4236/aces.2019.94023.

1. Introduction

A metalloligand can be defined as a complex that contains several potential donor groups that enable coordination to a variety of metal ions. Among the well-established metalloligands, interest in the chemistry of dithiolene complexes has increased tremendously in the past five decades since their initial popularity in the 1960s [1]. Early interest focused primarily on the structural geometries [2], redox properties [3] [4] [5] and magnetic properties [6] of this class of complexes, which arise from the noninnocent property [7] [8] of dithiolenes. Ward and McCleverty [7] have pointed out that the term noninnocent is applied properly when it is referred to a particular combination of the metal and the ligand rather than to redox-active ligands only.

By taking into account the dithiolene ligands, two main classes of molecules can be distinguished. Such as the non-benzenoid dithiolenes, where the p-delocalization is confined in the dithiolene core [9] and benzenoid systems, where the p-delocalization is extended to the aromatic ring [10]. Dithiolene ligands often are referred to as noninnocent when coordinated to transition metals [11] [12]. Let us turn our attention to the chemistry of tris (benzene-1,2-dithiolato) complexes. In contrast to their well-studied redox behavior, bond-making and bond-breaking reactivities of tris(benzene-1,2-dithiolato) complexes have been paid less attention [4]. The term “dithiolene” (L) will be used for two classes of ligands as shown in scheme 1(a) and 1(b) of Figure 1 [13], irrespective of their “true” oxidation level as monoanionic radical (L•)1−, or as dianionic, closed-shell (LRed)2− ligand or as neutral 1,2-dithioketone (LOx)0 [14]. Previously, Sproules et al. have calculated the geometry and electronic structures of the [V(L)3]z (z = 1+, 0, 1-, 2-, 3-, 4-) series by using density functional theory (DFT) and explored the stability of a trigonal prismatic structure over an octahedral one [15] [16]. No dianion containing ligand as shown in scheme 2(a) [13] in Figure 1 has been isolated previously. The most reduced, trianionic [Cr(dithiolene)3]3− form with an S = 3/2 ground state has been isolated only as [PPh4]3 [Cr(mnt)3] salt [17] [18] (as shown in scheme 2(b) [13] in Figure 1). Wieghardt et al. and others have successfully elucidated the electronic structures for the [Cr(L)3]z (z = 0, 1-, 2-, 3-) [19] [20] electron transfer series and found a Werner-type compound with dithiolates. The following dianionic manganese

Figure 1. Binding schemes of metals coordination complexes [13].

(IV) complexes have also been isolated: [PPh4]2 [Mn(mnt)3] and [N(nBu)4]2 [Mn(Cl4-bdt)3] both of which possess an S = 3/2 ground state [13] [18] [21] and only very recently (PNP)2 [Mn(bdt)3] (PNP+ = bis(triphenylphosphine)iminium) became the first structurally characterized tris(dithiolene)manganese compound [22]. They are isoelectronic with the corresponding trianions [Cr(dithiolene)3]3−. The corresponding tris(benzo-1,2-dithiolato)chromium complex [N(n-Bu)4] [CrIII(3,5LS,S)2(3,5LS,S)](S)1/2) has also been isolated: (3,5LS,S)2− represents the closed-shell dianion 3,5-di-tert-butylbenzene-1,2-dithiolate(2-), and (3,5LS,S)1− is its monoanionic ð radical complex. Milsmann et al. reported on a similar series of iron(III) complexes, where a redox-innocent, closed-shell cyclam ligand (1,4,8,11-tetraazacyclotetradecane) and a potentially redox-noninnocent dithiolene ligand such as toluene-3,4-dithiolate or maleonitrile dithiolate are coordinated to a low-spin ferric ion, affording the octahedral monocations [23] [24]. Only a few such complexes containing a cis-[FeIIIN4S2] core structure have been reported previously: [FeIII(L-N4Me2)-(bdt)]+(L-N4Me2 = N,N0-dimethyl-2,11-diaza [3.3]-(2,6)pyridinophane) [25], [FeIII(tripod)(bdt)]+ tripod = tris [(2-pyridyl)methyl]amine, and bis [(2-pyridyl)methyl] [(1-methylimidazole-2-yl)-methyl]amine [26] [(bdt)2− is unsubstituted benzene-1,2-dithiolate(2-)]. The former low-spin ferric compound exhibits an unusual thermal spin transition from an S = 1/2 low-spin ground state to an S = 3/2 intermediate-spin excited state [25]. [CoIII(trien)(3,5-DTBsq•)]Cl2 (S = 1/2) is one of the earliest examples (1976) [27], whereas [CrIII(tren)-(3,6-DTBsq•)](PF6)2 is probably the most thoroughly studied [28] [29]. Milsmann et al. employed the DFT to establish the presence of an S,S0-coordinated benzene-1,2-dithiolate(1-) π radical in the two species [CoIII(tren)(tmsdt•)]2+ and [CrIII(tren)(tmsdt•)]2+ [30].

The bis(dithiolate) metal and tris(dithiolate) metal complexes have been extensively studied [2] and only very recently, the spectroscopic methods like S K-edge X-ray absorption spectroscopy have been developed to the extent that the presence of a benezene-1,2-dithiolate(1-) π radical in a coordination compound can be established beyond doubt [31].

In the present study, we have performed systematically a theoretical investigation on the structures, binding energies, spectroscopic and electronic properties of [M(bdt)3]3− (M = V, Cr, Mn , Fe and Co) which contains bdt2− = benzene-1,2-dithiolate ligand (as shown in scheme 2(a) of Figure 1) in the gas phase using the first principles HF/DFT hybrid approach. The natural bond orbital (NBO) analysis on the tris complexes of benzene-1,2-ditholate (bdt2−) complexes including V3+ to Co3+ is reported. These findings help us understand the thermodynamic behavior of each system as a function of the quantum chemistry descriptors.

2. Computational Details

We have studied the [tris(bdt)M]3− complexes where M = V3+, Cr3+, Mn3+, Fe3+ and Co3+ with the bdt2− = benzene-1,2-dithiolate ligand (Figure 2). The dithiolene was taken to be oxidized (dithiolate). The compounds are theoretically examined by means of HF/DFT hybrid approach B3LYP [32] with the 6-311G(d,p) basis sets in gas phase. All the calculations were performed using Gaussian 09 program [33]. GaussView 5.0.8 was used for the visualization of the structures and simulated vibrational spectra. We have considered dithiolate as a weak-field ligand, the metal coordination complex was assumed to have high spin multiplicity. Geometry optimization was taken to be converged if the maximum atomic force was smaller than 0.00045 Hartree/Bohr and there was no imaginary frequency. No symmetry was imposed in all the calculations. After optimizing the geometry of the coordination complexes, the metal-ligand binding energies were calculated as [34] [35]

Figure 2. Metal binding to benzene-1,2-dithiolate, optimized geometry of the [Fe(bdt)3]3− complex.

Δ E = ( E complex E metal E ligand ) 3 (1)

where E complex , E metal , and E ligand are the energies of the tris(dithiolate) coordination complex, the metal ion and dithiolate ligand, respectively. Therefore, ΔE refers to the binding energy per ligand. The energy was calculated by optimizing the geometry of the complex and the ligand separately.

Electronically excited state calculation was carried out to compute the UV-Vis of the [M(dithiolate)3]3− complexes of V3+, Cr3+, Mn3+, Fe3+ and Co3+. The vertical excitation energies were obtained using the time-dependent DFT at the level of CAM-B3LYP [36] /6-311 + G(d,p) after the ground state optimization.

3. Results and Discussion

3.1. Structural Analysis

After full optimization, the structural parameters such as bond distances and angles of the [tris(bdt)M]3− (M = V3+, Cr3+, Mn3+, Fe3+ and Co3+) complexes are calculated. The calculated geometrical parameters as depicted in Figure 2 and Table 1. The metal and the sulfur atoms are 2.51 Å apart for V, 2.48 Å for Cr, 2.55 Å for Mn, 2.54 Å for Fe and 2.60 Å for Co. The DFT and X-ray analysis of the tris-dithiolate complexes reported V-S distances of 2.337 Å and 2.350 Å respectively [15], for Cr-S distances were 2.283 Å and 2.299 Å [19]. The experimental average Cr-S distance in the monoanion of the tris complexes at 2.299 ± 0.003 Å, is shorter by 0.086 Å than the calculated one. The experimental average Cr-S distance in monoanion of the tris complex is 0.175 Å shorter than that for our DFT calculation of the [tris(bdt)Cr]3− complex. The Mn-S distances are 2.383 Å and 2.347 Å respectively in the [MnIV(Cl2-bdt)3] [Net4]2. CH2Cl2 complex and 2.331 Å for the [Mn(bdt)3]2− complex [13]. The DFT calculated [tris(bdt)Mn]3− complex is 2.55 (±0.17) and 0.17 Å higher than experiment one [13]. The calculated Fe-S bond distance of the tris complex is 2.54 Å which is 0.29 Å higher than the experimental value of the complex of [FeIII(cyclam)(tdt)] (PF6) where the Fe-S distance is 2.25Å [23]. The experimental value deviates due to present

Table 1. Geometric parameters of the tris(benzene-1,2-dithiolato) complexes of various transition metal ions. The metal-S bond lengths (dMS’s), C-S bond lengths (dCS’s) and the S-M-S bending angles (θSMS’s), the C-C-S bending angles (θCCS’s) and Metal-S-C bending angles (θMSC’s) are shown for trivalent metal ions ranging from V3+ to Co3+. The average values are listed with the standard deviations in parentheses.

bulky substitute group cyclam. The theoretically calculated average Co-S bond lengths is 2.60 Å for the tris complex whereas experimentally measured value was found 2.25 Å [30] of the complex [CoIII(tren)(bdt)](PF6)1/2(Cl)1/2∙H2O. This deviation was due to present bulky groups tren and bdt in the complex [CoIII(tren)(bdt)](PF6)1/2(Cl)1/2∙H2O.

The present S-M-S angles were significantly larger than the previous X-ray crystal structures. The average angle is nearly octahedral angle 90˚. The angles formed by the metal ion and the two sulfur atoms of the dithiolate ligands are within 90.30˚ to 90.80˚ on an average indicating almost octahedral geometry in the coordination. The average S-Cr-S angle of the complex [N(n-Bu)4] [CrIII(3,5LS,S)2(3,5LS,S)] is 84.22˚ [19] whereas theoretically calculated value of the present tris complex of [Cr(bdt)3]3− is 90.80˚. The experimental S-Mn-S angles of the [MnIV(Cl2-bdt)3] [NEt4]2∙CH2Cl2 and [Mn(S,S-C6H4)3]2− complexes was 85.2˚ and 92.28˚ respectively [13] The present DFT calculation estimates 90.30˚ (±6.66) which are within the experimental range. The theoretical S-Fe-S angles in our calculation was 90.47˚, in good agreement with experimentally measured value (90.72˚) in the [(n-Bu)4N]2 [Fe(bdt)2]2 complex [37]. The DFT calculated average bond angle S-Co-S was found 90.51˚. Figure 2 shows the average M-S bond distance’s dMS’s and the angles of S-M-S triplets, θSMS’s. The fluctuation in dMS was significant. All the, θSMS’s however fluctuate significantly from their average values.

3.2. Binding Energy

In this study we have investigated the metal-ligand binding energies, ΔE, of the trivalent metal ions considered and are presented in Figure 3 and Table 2. The ΔE values have greatly exceeded the physiorption energy of benzene-1,2-dithiol adsorbed on metallic surface. Lee et al. reported that a single benzene-1,2-dithiol molecule binds to a gold (111) surface with a binding energy of 6.0 kcal/mol [38]. A similar energy was reported in a DFT study of the adsorption of benzene-1, 2-dithiol to silver (111) surface (3.30 kcal/mol) [38]. The present ΔE

Figure 3. Metal-ligand binding energies of the coordination complexes of various metal ions for tris [M(bdt)3]3− (M = V to Co) complexes.

Table 2. The calculated metal-ligand binding energies ΔE, zero point corrected binding energies, ΔEZPEs, thermal energies ΔEtots, enthalpies ΔHs and Gibbs free energies ΔGs. The metal-ligand binding energies ΔEs were corrected by using the thermodynamic condition, 298.15 K and 1 atm. All energies are in units of kcal/mol.

values (467.19 - 497.20 kcal/mol) are an order of magnitude larger than these. Interestingly, the previous study also reported that trivalent 1st row transition metal coordination produced a 1 (one) order of magnitude increase in the cross-linking strength [39]. Figure 3 shows that the overall ΔE increases with increasing nuclear charge due to the increased electrostatic interaction with increasing nuclear charge. The decreasing ΔE with changing Cr to Mn arises from the Jahn-Teller distortion of the metal complex. The dip in ΔE observed for Fe is ascribed to the d5 configuration of Fe (III). In this case, as all of the d orbitals are occupied, the ligand-to-metal charge transfer is inefficient resulting in a reduced ΔE. The peak at Cr can be understood by noting that the ligand field stabilization should be largest for the d3 configurations of a high-spin complex. We have calculated the binding energies by including the zero-point energies (ZPEs), thermal energies, enthalpies and Gibbs free energies as summarized in Table 2. The vibrational, thermal, and entropic contributions to ΔE values turned out to be small, presumably due to the covalent nature of the metal ligand binding. Regardless of the metal ion, the ZPE-, thermal-energy-, enthalpy-, Gibbs free energy corrected ΔE value were all within 6% of the uncorrected binding energies.

3.3. Spectroscopic Data

In order to assign the electronic absorption bands of the complexes, TD-DFT calculations have been carried out on these complexes in gas phase optimized geometries at the CAM-B3LYP [34] /6-311 + G(d,p) level of theory. TD-DFT is a useful method for studying excitation energies, and its application has increased in the recent years. The electronic absorption spectra of complexes formed by V, Cr, Mn, Fe and Co have been calculated. The calculated electronic absorption spectra of the tris complexes formed by V to Co are presented in Figure 4 and Table 3. Transition metal complexes generally show three types of electronic excitation bands which cover a wide wavelength range: d-d (crystal-field) transitions (300 - 1500 nm); metal-to-ligand charge-transfer (MLCT) and ligand-to-metal charge transfer (LMCT) transitions (200 - 500 nm) [40]. The LMCT transitions localized on the ligands commonly known as intra-ligand charge transfer (ILCT) transitions, which regularly occur in the ultraviolet region. ILCT results

Figure 4. UV-vis spectrum of the tris [M(bdt)3]3− (M = V to Co) complexes obtained from the present calculation.

Table 3. Electronic spectra of [tris(bdt)M]3− (M = V, Cr, Mn, Fe and Co) complexes.

from n → π and n → π* transitions and is affected by the type of coordination. The calculated spectra of the model compound [V(bdt)3]3− are 431, 448, 522 and 691 respectively whereas in [PPh4]2 [V(bdt)3] where vanadium has (IV) oxidation state showed the absorption maxima at 878, 681, 553, 426 and 327 respectively [15]. The calculated spectra of the model compound [Cr(bdt)3]3− are 442 nm and two d-d transitions with relatively low intensities in the visible region at 565 and 687 nm, as typical for chromium (III) octahedral complex. The electronic structure of this complex is best described as [CrIII(3,5LS,S)3]3− with a chromium (III) central ion (d3, S = 3/2) and three closed-shell dithiolate-(2-) ligands [19]. The electronic spectra of these tris(dithiolene)chromium complexes resemble closely those reported previously for the corresponding tris(dioxolene) chromium complexes [19].

For the trianion manganese [tris(bdt)Mn]3− complex, the absorption spectra were obtained at 453, 488 and 559 nm. These are in the range of experimentally measured absorption spectra (363, 508, 582 and 741 nm) of the [Mn(S,S-C6H4)3]2− complex (where Mn is d3) in CH3CN solution [22]. Our DFT calculated absorption spectra of the tris complex of iron [Fe(bdt)3]3− was found at 394, 546 and 650 nm respectively. The experimentally measured values of the complex of toluene-3,4-Dithiolate(tdt), [FeIII(cyclam)(tdt)](PF6) complex, the monocation, shows absorption bands above 400 nm with maximum at 580 nm and 825 nm respectively [23]. The calculated absorption spectra of the Co tris complex [tris(bdt)Co]3− was found at 687, 731 and 863 nm which are overestimated compared to the measured spectra of the [Co(III)(tren)(bdt)](PF6)1/2(Cl)1/2∙H2O complex at 463 nm and 646 nm respectively [30]. This discrepancies due to present bulky groups tren and bdt of the complex [Co(III)(tren)(bdt)](PF6)1/2(Cl)1/2∙H2O.

3.4. Electronic Properties

The electronic structure parameters are very important for understanding the molecular interactions with other species. We have calculated the electronic energy gap Eg, electron affinity A, electronegativity χ, chemical hardness η and chemical softness S for all the considered complexes. The results are presented in Table 4. The energetics of frontier molecular orbitals (FMOs), i.e., the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) energies are calculated. HOMO, which can be considered of as the outermost orbital containing electrons, tends to give these electrons and acts as an electron donor. On the other hand LUMO can be considered as the innermost orbital containing free places and accepts electrons. The difference between HOMO and LUMO provides the electronic gap Eg as depicted in Figure 5. Our investigation shows that Eg for the considered complexes increases in the following direction Mn3+ > V3+ > Cr3+ > Co3+ > Fe3+. The electron affinity and electro negativity values also follow the same trend (Table 4). From MOs plots (Figure 5), it can be concluded that the HOMO is mainly centered on metal coordinated sulfur atoms (which acts as a donor atoms) whereas LUMO is mainly located on the metal surface.

The Energy gap characterizes the molecular chemical stability, and it is a critical parameter in determining molecular electrical transport properties because it is a measure of electron conductivity. The energy gap between HOMO and LUMO (Eg) also plays important role in predicting the polarizability of a molecule. Smaller the energy gap, the more polarizable the molecule is. The chemical stability of molecule can be also studied by the calculation of the chemical hardness η and chemical softness S as shown in Table 4. Large energy gap is an indication of hard molecule and small energy gap is the sign of soft molecule. Therefore, among the [tris(bdt)M]3− (M = V, Cr, Mn, Fe and Co) complexes, the

Table 4. The electronic parameters calculated for the tris(benzene-1,2-dithiolato) complexes of considered transition metals.

Figure 5. Frontier MOs of the tris(benzene-1,2-dithiolato) complexes. The H, C and S atoms are shown as white, gray, and yellow spheres, respectively. The orbital lobes shown in green and red represent the opposite phases.

[tris(benzene-1,2-dithiolato)Fe]3− complex is the hardest, followed by those of Co3+, Cr3+, V3+ and Mn3+ respectively. The soft molecules are more polarizable than the hard ones as they need small energy for excitation.

3.5. Atomic Charges

The atomic charge on the metal ion has been calculated and summarized in Table 5. Four different schemes were used to calculate the charges: Natural Population Analysis (NPA) [41], MK [42] (Merz-Singh-Kollman), CHelpG [43] (Charges from ELectrostatic Potentials using a Grid based method), and CHelp [44] methods. The atomic charges vary according to the schemes. For example, according to NPA method, the charge is 1.513 for Fe, but the charge is 2.038 if the MK scheme was used. The atomic charge decreased with increasing nuclear charge. The present drift of the atomic charge deviated considerably from a uniform decrease with increasing nuclear charge. The uniform decrease was attributed to the increased covalent nature of the metal-sulfur bond (hence, the charge transfer increased from the ligand to metal). For the present complexes, however, the atomic charge was a minimum for Co3+, regardless of the charge scheme. The charge also decreased as the metal atom was changed from Cr to Mn and Fe, according to the MK, CHelpG and CHelp schemes.

3.6. NBO Analysis

In the natural bond orbital (NBO) [45] analysis, electronic wave function is explained in terms of occupied Lewis and unoccupied Lewis localized orbitals. The delocalization of electron density(ED) between occupied Lewis-type and unoccupied non-Lewis NBOs correspond to stabilizing donor-acceptor interactions that contribute predominantly to the stabilization of the entire molecular system. The strength of donor-acceptor interactions, E i j ( 2 ) , are evaluated by second-order perturbation theory. For each donor (i) and acceptor (j) in the complexes, the stabilization energy or second-order perturbation energy, E i j ( 2 ) associated with the delocalization from i j was estimated using Equation (2).

Δ E i j ( 2 ) = q ( F ^ i j ) 2 ε j ε i (2)

In Equation (2), q is the donor orbital occupancy, ε i and ε j are diagonal elements (orbital energies) of donor and acceptor NBOs respectively and F ^ i j is the off-diagonal NBO Fock matrix element. Values of this energy are proportional to the intensities of NBO interactions or to the extent of intramolecular charge transfer (ICT) within a molecular entity. The greater the electron donating tendency from donor to acceptor NBOs is, the larger the E i j ( 2 ) values and the more intense the interaction between the electron donors and the electron acceptors are. Table 6 summarizes the E i j ( 2 ) values of the tris(benzene-1,2-dithiol) complexes of trivalent V, Cr, Mn, Fe and Co complexes in gas phase for the important NBO interactions. The most important interaction energies of these complexes are due to interactions between the lone pair electron of the S atom(LP0) and antibonding orbital of the M3+( LP M * ) (M3+ = V, Cr, Mn, Fe and Co). According to Table 6, E i j ( 2 ) values decrease with a decrease in the ionic sizes. Therefore, donor-acceptor interaction energies are the greatest values in the tris(benzene-1,2-dithiol) complex of trivalent V ion except Mn. The donor-acceptor interaction value of Mn complex is higher due to Jahn-Teller

Table 5. Atomic charges on the metal atoms obtained using various charge schemes, NPA, MK, CHelpG and CHelp schemes.

aNatural Population Analysis (NPA) [41]. bMK(Merz-Singh-Kollman) [42]. cCHelpG (CHarges from ELectrostatic Potentials using a Grid based method) [43]. dCHelp methods to fit the electrostatic potential method [44].

Table 6. Significant donor-acceptor interaction energies of [M(benzene-1,2-dithiolato)3]3− (M = V3+, Cr3+, Mn3+, Fe3+ and Co3+) complexes.

effect. The strongest interaction in V-complex is LPS12 → LPV, in Cr-complex is LPS37 → LPCr, in Mn-complex is LPS12 → LPMn. In Fe-complex the strongest interaction is LPS37 → LPFe and in Co-complex is LPS36 → LPCo. The results of NBO analysis reflect generally charge transfer from lone pair orbitals located on the donor atoms S to the central metal ions. NBO analysis provides the most accurate possible “natural Lewis structure”.

4. Conclusion

The structural, energetic, spectroscopic and electronic properties of [M(benzene-1,2-dithiolato)3]3− (M = V3+, Cr3+, Mn3+, Fe3+ and Co3+) complexes have been investigated theoretically using quantum chemical calculations based on HF/DFT hybrid approach. The calculated geometrical parameters of the tris (benzene-1,2-dithiolato) complexes are in good agreement with the experimental measurement values. The binding energies of the tris complexes range from 467.19 - 497.20 kcal∙mol−1 which are 1 (one) order higher than physically adsorbed on metallic surface. TD-DFT calculations have been successively employed to simulate the electronic spectra of the V, Cr, Mn, Fe and Co complexes and facilitated a transition assignment. With the electronic structure of this tris-complexes, we have begun to make meaningful comparisons to experimental spectra for investigation into many tris(dithiolato) complexes of early transition metals. The investigation of electronic parameters including frontier molecular orbital energies suggests that Fe containing complexes should have the largest band gap and therefore being the hardest among all. The results of NBO analysis reflect charge transfer from lone pair orbitals located on the donor atoms to the central metal ions. The present metal-ligand binding energies, structures, stability and atomic charges of the metal-benzene-1,2-ditholato (bdt2−) complexes will serve as a cornerstone for such modeling using molecular dynamics or Monte Carlo simulations.


Authors gratefully thank the Ministry of Science and Technology, Bangladesh for financial support of this study. We also express our sincere thanks to CARS, DU for providing support for this research.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Nomura, M., Cauchy, T. and Fourmigué, M. (2010) Radical Cpni(Dithiolene) and Cpni(Diselenolene) Complexes: Synthetic Routes and Molecular Properties. Coordination Chemistry Reviews, 254, 1406-1418.
[2] Beswick, C.L., Schulman, J.M. and Stiefel, E.I. (2004) Structures and Structural Trends in Homoleptic Dithiolene Complexes. In: Stiefel, E.I., Ed., Dithiolene Chemistry: Synthesis, Properties, and Applications, Volume 52, John Wiley & Sons, Inc., Hoboken, 52-55.
[3] Kirk, M.L., McNaughton, R.L. and Helton, M.E. (2004) The Electronic Structure and Spectroscopy of Metallo-Dithiolene Complexes. In: Stiefel, E.I., Ed., Dithiolene Chemistry: Synthesis, Properties, and Applications, Volume 52, John Wiley & Sons, Inc., Hoboken, 111-212.
[4] Wang, K. (2004) Electrochemical and Chemical Reactivity of Dithiolene Complexes. In: Stiefel, E.I., Ed., Dithiolene Chemistry: Synthesis, Properties, and Applications, Volume 52, John Wiley & Sons, Inc., Hoboken, 267-314.
[5] Sellmann, D. and Sutter, J. (2004) Dithiolenes in More Complex Ligands. ChemInform, 35, 585-681.
[6] Cassoux, P., Valade, L., Kobayashi, H., Kobayashi, A., Clark, R.A. and Underhill, A.E. (1991) Molecular Metals and Superconductors Derived from Metal Complexes of 1,3-Dithiol-2-Thione-4,5-Dithiolate (Dmit). Coordination Chemistry Reviews, 110, 115-160.
[7] Ward, M.D. and McCleverty, J.A. (2002) Non-Innocent Behaviour in Mononuclear and Polynuclear Complexes: Consequences for Redox and Electronic Spectroscopic Properties. Journal of the Chemical Society, Dalton Transactions, 3, 275-288.
[8] Fernández, L., Fernando Pérez-Pla, F., Tuñón, I. and Llopis, E. (2016) DFT Study on the Interaction of Tris(benzene-1,2-dithiolato)molybdenum Complex with Water. A Hydrolysis Mechanism Involving a Feasible Seven-Coordinate Aquomolybdenum Intermediate. The Journal of Physical Chemistry A, 120, 9636-9646.
[9] Lim, B.S., Fomitchev, D.V. and Holm, R.H. (2001) Nickel Dithiolenes Revisited: Structures and Electron Distribution from Density Functional Theory for the Three-Member Electron-Transfer Series [Ni(S2c2me2)2]0,1-,2. Inorganic Chemistry, 40, 4257-4262.
[10] Sellmann, D., Binder, H., Häußinger, D., Heinemann, F.W. and Sutter, J. (2000) Transition Metal Complexes with Sulfur Ligands: Part Cxliv. Square Planar Nickel Complexes with Nis4 Cores in Three Different Oxidation States: Synthesis, X-Ray Structural and Spectroscopic Studies. Inorganica Chimica Acta, 300-302, 829-836.
[11] Eisenberg, R. and Gray, H.B. (2011) Noninnocence in Metal Complexes: A Dithiolene Dawn. Inorganic Chemistry, 50, 9741-9751.
[12] Periyasamy, G., Burton, N.A., Hillier, I.H., Vincent, M.A., Disley, H., McMaster, J. and Garner, C.D. (2007) The Dithiolene Ligand-“Innocent” or “Non-Innocent”? A Theoretical and Experimental Study of Some Cobalt-Dithiolene Complexes. Faraday Discussion, 135, 469-488.
[13] Banerjeee, P., Sproules, S., Weyhermüller, T., George, S.D. and Wieghardt, K. (2009) Electronic Structure of the [Tris(dithiolene)chromium]z (z = 0, 1−, 2−, 3−) Electron Transfer Series and Their Manganese(IV) Analogues. An X-Ray Absorption Spectroscopic and Density Functional Theoretical Study. Inorganic Chemistry, 48, 5829-5847.
[14] McCleverty, J.A. (2007) Metal 1,2-Dithiolene and Related Complexes. In: Cotton, F.A., Ed., Progress in Inorganic Chemistry, Volume 10, John Wiley & Sons, Inc., Hoboken, 49-221.
[15] Sproules, S., Weyhermüller, T., DeBeer, S. and Wieghardt, K. (2010) Six-Membered Electron Transfer Series [V(Dithiolene)3]Z (Z = 1+, 0, 1−, 2−, 3−, 4−). An X-Ray Absorption Spectroscopic and Density Functional Theoretical Study. Inorganic Chemistry, 49, 5241-5261.
[16] Graham, M.J., Krzyaniak, M.D., Wasielewski, M.R. and Freedman, D.E. (2017) Probing Nuclear Spin Effects on Electronic Spin Coherence via EPR Measurements of Vanadium (IV) Complexes. Inorganic Chemistry, 56, 8106-8113.
[17] Davison, A., Edelstein, N., Holm, R.H. and Maki, A.H. (1964) Synthetic and Electron Spin Resonance Studies of Six-Coordinate Complexes Related by Electron-Transfer Reactions. Journal of the American Chemical Society, 86, 2799-2805.
[18] McCleverty, J.A., Locke, J., Wharton, E.J. and Gerloch, M. (1968) Transition-Metal Dithiolene Complexes. Part IV. Trisdicyano-1, 2-Dithiolene Complexes of First-Row Metals, and of Molybdenum and Tungsten. Journal of the Chemical Society A, 816-823.
[19] Kapre, R.R., Bothe, E., Weyhermüller, T., DeBeer George, S., Muresan, N. and Wieghardt, K. (2007) Electronic Structures of Tris(Dioxolene)Chromium and Tris(Dithiolene)Chromium Complexes of the Electron-Transfer Series [Cr(Dioxolene)3]Z and [Cr(Dithiolene)3]Z (Z = 0, 1−, 2−, 3−). A Combined Experimental and Density Functional Theoretical Study. Inorganic Chemistry, 46, 7827-7839.
[20] Sproules, S., Benedito, F.L., Bill. E., Weyhermüller, T., DeBeer George, S. and Wieghardt, K. (2009) Characterization and Electronic Structures of Five Members of the Electron Transfer Series [Re(Benzene-1,2-Dithiolato)3]Z (Z = 1+, 0, 1−, 2−, 3−): A Spectroscopic and Density Functional Theoretical Study. Inorganic Chemistry, 48, 10926-10941.
[21] Wharton, E.J. and McCleverty, J.A. (1969) Transition-Metal Dithiolenes. Part VIII. Metal Complexes of Tris-Tetrachlorobenzene-1, 2-Dithiol and Related Compounds. Journal of the Chemical Society A, 2258-2266.
[22] Lin, C.-H., Chen, C.-G., Tsai, M.-L., Lee, G.-H. and Liaw, W.-F. (2008) Monoanionic {Mn(No)}5 and Dianionic {Mn(No)}6 Thiolatonitrosylmanganese Complexes: [(No)Mn(L)2]− and [(No)Mn(L)2]2− (Lh2 = 1,2-Benzenedithiol and Toluene-3,4-Dithiol). Inorganic Chemistry, 47, 11435-11443.
[23] Milsmann, C., Patra, G.K., Bill, E., Weyhermüller, T., DeBeer George, S. and Wieghardt, K. (2009) Octahedral Monodithiolene Complexes of Iron: Characterization of S,S’-Coordinated Dithiolate(1−) Π Radical Monoanions: A Spectroscopic and Density Functional Theoretical Investigation. Inorganic Chemistry, 48, 7430-7445.
[24] Castillo, O., Delgada, E. Hernández, D., Hernández, E. and Martìn, A. (2019) Synthesis and Crystal Structures of Ion-Pairs Based on Anionic Irondithiolenes and Alkylammonium as Counteraction. Journal of Molecular Structure, 1196, 323-331.
[25] Koch, W.O., Schünemann, V., Gerdan, M., Trautwein, A.X. and Krüger, H.-J. (1998) Evidence for an Unusual Thermally Induced Low-Spin (S = 1/2) Intermediate-Spin (S = 3/2) Transition in a Six-Coordinate Iron(III) Complex: Structure and Electronic Properties of a (1,2-Benzenedithiolato)Iron(III) Complex Containing N,N’-Dimethyl-2,11-Diaza [3.3](2,6)Pyridinophane as Ligand. Chemistry—A European Journal, 4, 686-691.<686::AID-CHEM686>3.0.CO;2-P
[26] Merkel, M., Pascaly, M., Wieting, M., Duda, M. and Rompel, A. (2003) Synthesis and Characterization of the First (1,2-Benzenedithiolato)iron(III) Complexes with Tripodal Ligands. Journal of Inorganic and General Chemistry, 629, 2216-2221.
[27] Wicklund, P.A., Beckmann, L.S. and Brown, D.G. (1976) Preparation and Properties of a Stable Semiquinone Complex. Inorganic Chemistry, 15, 1996-1997.
[28] Wheeler, D.E. and McCusker, J.K. (1998) Electron Exchange and the Photophysics of Metal-Quinone Complexes. 1. Synthesis and Spectroscopy of Chromium-Quinone Dyads. Inorganic Chemistry, 37, 2296-2307.
[29] Rodriguez, J.H., Wheeler, D.E. and McCusker, J.K. (1998) Density Functional Studies of a Heisenberg Spin Coupled Chromium-Semiquinone Complex and Its Chromium-Catechol Analog. Journal of the American Chemical Society, 120, 12051-12068.
[30] Milsmann, C., Bothe, E., Bill, E., Weyhermüller, T. and Wieghardt, K. (2009) Octahedral Monodithiolene Complexes of Cobalt(III) and Chromium(III). Spectroscopic and Density Functional Theoretical Characterization of S,S’-Coordinated Benzene-1,2-Dithiolate(1−) Π Radicals. Inorganic Chemistry, 48, 6211-6221.
[31] Sarangi, R., DeBeer George, S., Rudd, D.J., Szilagyi, R.K., Ribas, X., Rovira, C., Almeida, M., Hodgson, K.O., Hedman, B. and Solomon, E.I. (2007) Sulfur K-Edge X-Ray Absorption Spectroscopy as a Probe of Ligand-Metal Bond Covalency: Metal vs Ligand Oxidation in Copper and Nickel Dithiolene Complexes. Journal of the American Chemical Society, 129, 2316-2326.
[32] Becke, A.D. (1993) Density-Functional Thermochemistry. III. The Role of Exact Exchange. Journal of Chemical Physics, 98, 5648-5652.
[33] Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G.A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H.P., Izmaylov, A.F., Bloino, J., Zheng, G., Sonnenberg, J.L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery Jr., J.A., Peralta, J.E., Ogliaro, F., Bearpark, M., Heyd, J.J., Brothers, E., Kudin, K.N., Staroverov, V.N., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J.C., Iyengar, S.S., Tomasi, J., Cossi, M., Rega, N., Millam, J.M., Klene, M., Knox, J.E., Cross, J.B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R.E., Yazyev, O., Austin, A.J., Cammi, R., Pomelli, C., Ochterski, J.W., Martin, R.L., Morokuma, K., Zakrzewski, V.G., Voth, G.A., Salvador, P., Dannenberg, J.J., Dapprich, S., Daniels, A.D., Farkas, O., Foresman, J.B., Ortiz, J.V., Cioslowski, J. and Fox, D.J. (2010) Gaussian 09, Revision B.01. Gaussian Inc., Wallingford.
[34] Myradalyyev, S., Limpanuparb, T., Wang, X. and Hirao, H. (2013) Comparative Computational Analysis of Binding Energies between Several Divalent First-Row Transition Metals (Cr2+, Mn2+, Fe2+, Co2+, Ni2+, and Cu2+) and Ligands (Porphine, Corrin, and Tmc). Polyhedron, 52, 96-101.
[35] Matin, M.A., Islam, M.M., Bredow, T. and Aziz, M.A. (2017) The Effects of Oxidation States, Spin States and Solvents on Molecular Structure, Stability and Spectroscopic Properties of Fe-Catechol Complexes: A Theoretical Study. Advances in Chemical Engineering and Science, 7, 137-153.
[36] Yanai, T., Tew, D.P. and Handy, N.C. (2004) A New Hybrid Exchange-Correlation Functional Using the Coulomb-Attenuating Method (Cam-B3lyp). Chemical Physics Letter, 393, 51-57.
[37] Cerebra, A.C., Simão, D., Santos, I.C., Machado, A., Pereira, L.C.J., Waerenborgh, J.C., Henriques, R.T. and Almeida, M. (2008) (N-Bu4n) [Fe(Cbdt)2]: Synthesis, Crystal Structure and Magnetic Characterisation of a New Fe III Bis-Dithiolene Complex. Inorganica Chimica Acta, 361, 3836-3841.
[38] Lee, Y.J., Jeon, I.C., Paik, W.-K. and Kim, K. (1996) Self-Assembly of 1,2-Benzenedithiol on Gold and Silver: Fourier Transform Infrared Spectroscopy and Quartz Crystal Microbalance Study. Langmuir, 12, 5830-5837.
[39] Matin, M.A., Chitumalla, R.K., Lim, M., Gao, X. and Jang, J.K. (2015) Density Functional Theory Study on the Cross-Linking of Mussel Adhesive Proteins. Journal of Physical Chemistry B, 119, 5496-5504.
[40] Mohammed, L.S., Hamza, I.S., AL-Deen, F.R.M. and Muhyedeen, B.R.J. (2014) DFT and MP2 Study of Pd(II) and Ni(II) PhCN, DMSO and Dithiooxamide Complexes—Part II: Theoretical. Journal of Applicable Chemistry, 3, 2102-2122.
[41] Reed, A.E., Weinstock, R.B. and Weinhold, F. (1985) Natural Population Analysis. Journal of Chemical Physics, 83, 735-746.
[42] Singh, U.C. and Kollman, P.A. (1984) An Approach to Computing Electrostatic Charges for Molecules. Journal of Computational Chemistry, 5, 129-145.
[43] Breneman, C.M. and Wiberg, K.B. (1990) Determining Atom-Centered Monopoles from Molecular Electrostatic Potentials. The Need for High Sampling Density in Formamide Conformational Analysis. Journal of Computational Chemistry, 11, 361-373.
[44] Chirlian, L.E. and Francl, M.M. (1987) Atomic Charges Derived from Electrostatic Potentials: A Detailed Study. Journal of Computational Chemistry, 8, 894-905.
[45] Weinhold, F. and Landis, C.R. (2001) Natural Bond Orbitals and Extensions of Localized Bonding Concepts. Chemistry Education Research and Practice, 2, 91-104.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.