Palash J. Thakuria,
Kangkan Sarmah,
Siddharth K. Purkayastha,
Amit Das and
Ankur K. Guha
*
Advanced Computational Chemistry Centre, Cotton University, Panbazar, Guwahati, Assam 781001, India. E-mail: ankurkantiguha@gmail.com
First published on 2nd July 2025
Planar pentacoordinate sulfur atoms (ppS) are rare, with only one example experimentally verified to date. Herein, we report a planar pentacoordinate sulfur atom (ppS) as the global minimum [Mg5S6]2− cluster. Despite dianionic nature, it is stable towards spontaneous electron detachment making it a suitable candidate for experimental detection. The global minimum of this cluster features perfect D5h symmetry. Detailed electronic structure calculations reveal the presence of significant covalency in the Mg–S bonds.
Extending the coordination of the carbon atom in a plane from four to five4 and six has also been predicted.22,23 Planar pentacoordinate atoms (ppA) other than carbon have also been predicted to be GEMs.25,35–39 This included planar pentacoordinate nitrogen (ppN),25 planar pentacoordinate hydrogen (ppH),35 planar pentacoordinate beryllium (ppBe),36 planar pentacoordinate s-block metals (ppM)37 and planar pentacoordinate halogens (ppX).38 Very recently, planar pentacoordinate oxygen (ppO)39 has been predicted to be the GEM in an experimentally observed [Be5O6]2− cluster.40
Interestingly, the [Be5O6]2− cluster is the first example of ppO in an experimentally detected and computationally predicted cluster, although a high energy linear isomer was computationally predicted earlier.41 Recently, Merino et al. have predicted 35 global minimum geometries containing planar tetracoordinate oxygen atoms (ptO) stabilized by group 13 elements.42 Despite sulfur being in the same group of the periodic table, planar tetra or pentacoordination of sulfur is very rare. Herein, we propose a planar pentacoordinate sulfur atom (ppS) as the GEM of the [Mg5S6]2− cluster. To the best of our knowledge, only one experimental report (by Müller and Henkel) of ppS in the complex [Ni5S(StBu)5]− is there in the literature,43 although some examples of computationally designed ppS in the hydrometal complexes, Ag5H5S2+ and Au5H5S2+ are there in the literature.44 But it is unclear whether they are GEMs or not3 and hence, their experimental detection is still a matter of question.
The potential energy surface of the [Mg5S6]2− cluster was explored using the TPSS-D3(BJ)/def2-SVP level45–48 in combination with the ABCluster code.49,50 Low energy isomers were then re-optimized at the TPSS-D3(BJ)/aug-cc-pVTZ level. Harmonic frequency calculations reveal that all the structures are local minima. Low energy isomers were then energetically refined using single point calculations at the CCSD(T)/aug-cc-pVTZ51 level on top of TPSS-D3(BJ) optimized geometries. All these calculations were performed using the Gaussian 16 suite of programs.52 Electronic structure was analyzed using the adaptive natural density partitioning (AdNDP)53 scheme implemented in the Multiwfn program code.54
Fig. 1 shows the optimized GEM geometry of the [Mg5S6]2− cluster, which adopts D5h symmetry in the singlet ground state. Low lying isomers are shown in Fig. S1 (ESI†). All the structures have shown reliable T1 diagnostic values (<0.02)55 suggesting negligible multi-reference character. The closest isomer lies 10.1 kcal mol−1 higher in energy at the CCSD(T)/aug-cc-pVTZ//TPSS-D3(BJ)/aug-cc-pVTZ level. The lowest energy triplet state lies 58.3 kcal mol−1 higher in energy. The Mg–Scentral (Scentral means the central S atom) bond length is 2.60 Å with Mayer bond order of 0.28, while the terminal Mg–S distance is 2.33 Å with Mayer bond order of 0.67. The calculated Mayer bond order of 0.28 for the Mg–Scentral bond indicates significant covalency. The Mg–Mg distance is 3.06 Å with Mayer bond order of 0.14. The calculated Hirshfeld charge at the central S atom is −0.35 e, while Mg carries positive charge of 0.28 e. This indicates that significant electrostatic interaction is also present in the GEM of the [Mg5S6]2− cluster.
![]() | ||
Fig. 1 TPSS-D3(BJ)/aug-cc-pVTZ optimized global minimum geometry of the [Mg5S6]2− cluster. Bond lengths are in Å and Hirshfeld charges (in red font) are in electrons. |
The electronic stability of the dianion was assessed through the calculation of its vertical detachment energy (VDE) using the outer valence Green's function (OVGF) method at the OVGF/aug-cc-pVTZ level.56 The calculated first VDE is found to be positive (2.48 eV with a pole strength of 0.87) indicating that the dianion is stable towards spontaneous electron detachment. Similar positive VDE was also computed for the ppO [Be5O6]2− cluster,39 an experimentally detected dianion.40 Furthermore, the electronic structure was analyzed using the adaptive natural density partitioning scheme (AdNDP) at the TPSS-D3(BJ)/def2-TZVP level. This analysis provides a description of the electron distribution in nc-ne bonds and is a very important tool to analyse the electronic structure having unusual chemical bonds. As shown in Fig. 2, AdNDP recovers five 1c-2e S lone pairs, ten 2c-2e Mg–S σ bonds and five 3c-2e S–Mg–S π bonds (set A). All these orbitals have significant occupation numbers (1.95–1.98 |e|). In set B, one 1c-2e lone pair orbital at the central S atom was recovered with an occupation number of 1.90 |e|. In set C, three 1c-2e lone pairs at the central sulfur atom were recovered. However, their occupation numbers are significantly less (1.76–1.83 |e|). Note that the difference of 0.23–0.16 |e| between the occupancies of the lone pairs and 6c-2e bonds is not negligible. The lower occupation numbers in set C indicate that it may not be an adequate description of the electronic structure. An alternative scheme is also shown (set D) where three 6c-2e σ bonds are recovered with significant occupation number (1.99 |e|). Therefore, the alternative scheme (set D) is more appropriate description of the electronic structure. This is also in tune with the computed Mayer bond order of 0.28 for the Mg–Scentral bond, which indicated covalency. Moreover, the delocalization index – an atoms in molecules (AIM)57 – based indicator of covalency computed for the Mg–Scentral bond is significant (0.49), further supporting significant covalency.
To investigate the nature of interaction, the GEM structure has been further analysed using natural orbital for chemical valence extended to transition state (ETS-NOCV)58 using Multiwfn program code at the TPSS-D3(BJ)/def2-TZVP level. For this, the S2− dianion and neutral Mg5S5 fragments have been considered as interacting fragments. Fig. 3 shows the deformation densities corresponding to ΔEorb(n) along with the eigen values |νn| of charge transfer (from yellow to green). Significant orbital interactions have been observed where charge transfer taking place from the central S atom to Mg atoms has been observed (Fig. 3), suggesting covalency in the Mg–S bonds.
![]() | ||
Fig. 3 Deformation densities and orbital interaction energies (kcal mol−1) along with eigen values |ν| of charge transfer. Charge transfer is from yellow to green. |
In summary, the present quantum chemical establishes a planar pentacoordinate sulfur atom as the global minimum of the [Mg5S6]2− cluster. Despite being dianionic in nature, it is stable towards spontaneous electron detachment as revealed by its positive vertical detachment energy. Detailed electronic structure study reveals significant covalency in the Mg–S bonds. The presence of three 6c-2e σ bonds renders stability to the global minimum. We feel that owing to the binary nature of the cluster, its experimental detection is quite likely.
Footnote |
† Electronic supplementary information (ESI) available: It contains complete citation of ref. 51, Fig. S1 and Cartesian coordinates of the optimized geometries. See DOI: https://doi.org/10.1039/d5cc03041b |
This journal is © The Royal Society of Chemistry 2025 |