Revisiting the Structure and Bonding in Li5H6- and the Exploration of Reactivity: Planar Pentacoordinate Hydrogen

Recently, Guha and co-workers (Sarmah, K.; Kalita, A.; Purkayastha, S.; Guha, A. K. Pushing The Extreme of Multicentre Bonding: Planar Pentacoordinate Hydride. Angew. Chem. Int. Ed. 2024, e202318741) reported a highly intriguing bonding motif: planar pentacoordinate hydrogen (ppH) in Li5H6-, featuring C2v symmetry in the singlet state with two distinct H-Li (center-ring) bond distances. We herein revisited the potential energy surface of Li5H6- by using a target-oriented genetic algorithm. Our investigation revealed that the lowest-energy structure of Li5H6- exhibits a ppH configuration with very high D5h symmetry and a 1A1' electronic state. We did not find any electronic effect like Jahn-Teller distortion that could be responsible for lowering its symmetry. Moreover, our calculations demonstrated significant differences in the relative energies of other low-lying isomers. An energetically very competitive planar tetracoordinate hydrogen (ptH) isomer is also located, but it corresponds to a very shallow minimum on the potential energy surface depending on the used level of theory. Chemical bonding analyses, including AdNDP and EDA-NOCV, uncover that the optimal Lewis structure for Li5H6- involves H- ions stabilized by the Li5H5 crown. Surprisingly, despite the dominance of electrostatic interactions, the contribution from covalent bonding is also significant between ppH and the Li5H5 moiety, derived from H-(1s) → Li5H5 σ donation. Magnetically induced current density analysis revealed that due to minimal orbital overlap and the highly polar nature of the H-Li covalent interaction, the ppH exhibits local diatropic ring currents around the H centers, which fails to result in a global aromatic ring current. The coordination of Li5H6- with Lewis acids, BH3 and BMe3, instantly converts the ppH configuration to (quasi) ptH. These Lewis acid-bound ptH complexes show high electronic stability and high thermochemical stability against dissociation and, therefore, will be ideal candidates for the experimental realization.

Breaking the Hoff/Le Bel rule by an electron-compensation strategy: the global energy minimum of NGa4S4. Phys Chem Chem Phys 2024;26:3907-3911. [PMID: 38230710 DOI: 10.1039/d3cp05290g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

In tetracoordinate chemistry, there is an attractive scientific problem of how to make the planar configuration more stable than the tetrahedral configuration. For tetracoordinate nitrogen, the abundant studies indicate that the planar tetracoordinate nitrogen (ptN) is far less stable than the tetrahedral tetracoordinate nitrogen (ttN). Herein, we introduced four S atoms to the unstable ptN-NGa4+ and stable ttN-NGa4+ by following an electron-compensation strategy. Surprisingly, ptN-NGa4S4+ is more stable than ttN-NGa4S4+. Thermodynamically, ptN-NGa4S4+ is the global energy minimum, which is 46.7 kcal mol-1 lower in energy than ttN-NGa4S4+. Dynamically, the BOMD simulations indicated that ptN-NGa4S4+ has excellent dynamic stability at 4, 298, 500 and 1000 K, but the ttN-NGa4S4+ is isomerized at 1000 K. Electronically, the HOMO-LUMO gap of ptN-NGa4S4+ (6.91 eV) is much wider than that of ttN-NGa4S4+ (5.25 eV). Moreover, AdNDP analyses showed that the eight 2c-2e Ga-S σ-bonds eliminated the 4s2 lone pair/4s2 lone pair repulsion between the four Ga atoms and provided a strong spatial protection for ptN-NGa4S4+; and that the four 3c-2e Ga-S-Ga π back-bonds could compensate electrons for Ga, weakening the electron-deficiency of Ga. Simultaneously, the double 6σ/2π aromaticity further enhanced the stability of ptN-NGa4S4+. Thus, as the dynamically stable global energy minimum displaying double aromaticity, ptN-NGa4S4+ will be more promising than ttN-NGa4S4+ in gas phase generation.

Planar pentacoordinate s-block metals

The presence of a delocalized π-bond is often considered an essential criterion for achieving planar hypercoordination. Herein, we show that σ-delocalization could be sufficient to make the planar configuration the most stable isomer in a series of planar pentacoordinate s-block metals. High-level ab initio computations reveal that the global minimum of a series of interalkali and interalkali-alkaline earth clusters (LiNa5, Li5Mg+, Na5Mg+, K5Ca+, CaRb5+, Rb5Sr+, and SrCs5+) adopts a singlet D5h structure with a planar pentacoordinate lithium or alkaline earth metal (AE = Mg, Ca, Sr). These clusters are unusual combinations to stabilize a planar pentacoordinate atom, as all their constituents are electropositive. Despite the absence of π-electrons, Hückel's rule is fulfilled by the six σ-electrons. Furthermore, the systems exhibit a diatropic ring current in response to an external magnetic field and a strong magnetic shielding, so they might be classified as σ-aromatic. Therefore, multicenter σ-bonds and the resulting σ-delocalization stabilize these clusters, even though they lack π-aromaticity.

Planar Pentacoordinate Zinc Group Elements Stabilized by Multicentric Bonds

Planar pentacoordinate zinc group elements, (M = Zn, Cd, Hg) were computationally found to be at a global minimum in Li₅M⁺ clusters. The stability of these clusters is due to the presence of multicentric bonds. The central element (Zn, Cd, Hg) in each cluster features a negative oxidation state owing to the in-plane electron donation by the Li₅⁺ framework. A similar global minimum planar pentacoordinate structure is found in Na₅Zn⁺ and Na₅Cd⁺ clusters.

BAl4Mg−/0/+: Global Minima with a Planar Tetracoordinate or Hypercoordinate Boron Atom

We have explored the chemical space of BAl4Mg−/0/+ for the first time and theoretically characterized several isomers with interesting bonding patterns. We have used chemical intuition and a cluster building method based on the tabu-search algorithm implemented in the Python program for aggregation and reaction (PyAR) to obtain the maximum number of possible stationary points. The global minimum geometries for the anion (1a) and cation (1c) contain a planar tetracoordinate boron (ptB) atom, whereas the global minimum geometry for the neutral (1n) exhibits a planar pentacoordinate boron (ppB) atom. The low-lying isomers of the anion (2a) and cation (3c) also contain a ppB atom. The low-lying isomer of the neutral (2n) exhibits a ptB atom. Ab initio molecular dynamics simulations carried out at 298 K for 2000 fs suggest that all isomers are kinetically stable, except the cation 3c. Simulations carried out at low temperatures (100 and 200 K) for 2000 fs predict that even 3c is kinetically stable, which contains a ppB atom. Various bonding analyses (NBO, AdNDP, AIM, etc.) are carried out for these six different geometries of BAl4Mg−/0/+ to understand the bonding patterns. Based on these results, we conclude that ptB/ppB scenarios are prevalent in these systems. Compared to the carbon counter-part, CAl4Mg−, here the anion (BAl4Mg−) obeys the 18 valence electron rule, as B has one electron fewer than C. However, the neutral and cation species break the rule with 17 and 16 valence electrons, respectively. The electron affinity (EA) of BAl4Mg is slightly higher (2.15 eV) than the electron affinity of CAl4Mg (2.05 eV). Based on the EA value, it is believed that these molecules can be identified in the gas phase. All the ptB/ppB isomers exhibit π/σ double aromaticity. Energy decomposition analysis predicts that the interaction between BAl4−/0/+ and Mg is ionic in all these six systems.

Spontaneous bond dissociation cascades induced by Ben clusters (n = 2,4)

High-level single and multireference ab initio calculations show that the Be₄ cluster behaves as a very efficient Lewis acid when interacting with conventional Lewis bases such as ammonia, water or hydrogen fluoride, to the point that the corresponding acid-base interaction triggers a sequential dissociation of all the bonds of the Lewis base. Notably, this behavior is already found for the simplest beryllium cluster, the Be₂ dimer. However, whereas for Be₂ the first dissociation process involves a low activation barrier which is above the reactants, for Be₄ all the bond dissociation processes involve barriers below the entrance channel leading to a cascade of successive exothermic processes, which end up spontaneously in a global minimum in which the bonding patterns of both the base and the Lewis acid are completely destroyed. Indeed, the global minimum, in all cases, is stabilized by three-center Be-H-Be bonds and covalent interactions between the Be atoms and the basic center of the base, which replace the initial metallic bond stabilizing the Be₄ cluster. As a consequence, in the global minimum the basic atoms (N, O and F) behave as hyper-coordinated centers. Also importantly, the Be₄ cluster and its complexes present RHF-UHF instabilities (not reported before for Be₄), which require the use of multireference methods to correctly describe them.

Planar nitrogen compared with planar carbon in a series of compounds: the effect of replacing with nitrogen the potentially planar tetracoordinate carbon of a series of polycyclic alkanes

Previous work by Radom and Rasmussen had been extended by the author (E. Lewars. 2019. Can. J. Chem. 97(2): 154) to identify a smaller (C17H16 cf. C23H24) polycyclic alkane with a potentially planar tetracoordinate “central” carbon. Here, the series of alkanes in that latter study were computationally modified by replacing the central carbon with isoelectronic nitrogen, N+. A series of six ammonium ions was found by ab initio calculations to potentially all, with one exception, harbour a planar tetracoordinate nitrogen in a potential energy relative minimum at the MP2/6-31G* and MP2/(large basis) levels, as indicated by the absence of an imaginary vibrational frequency. This is in contrast to the hydrocarbons, where only two such structures were indicated to be relative minima.

Nitrogen versus carbon in planar pentacoordinate environments supported by Be5Hn rings

Nitrogen is a better fit for Be₅H_n rings, both geometrically and electronically, than carbon, leading to the viable planar pentacoordinate nitrogen species.

Planar Pentacoordinate Nitrogen in a Pseudo-Double-Aromatic NBe5H4+ Cluster

High-level quantum-chemical calculations have been used to predict a cationic ternary NBe₅H₄⁺ cluster containing a planar pentacoordinate nitrogen atom. The proposed cluster has pseudo dual aromaticity and is kinetically and thermodynamically very stable.

Avoided spin coupling: an unexpected σ-σ diradical in global planar pentacoordinate carbon

We present a global planar pentacoordinate carbon (ppC) featuring a hitherto unreported σ-σ diradical characteristic. Using the multi-reference approach combined with the CCSD(T)/aug-cc-pVTZ method, the ppC C₃Li₃^- was found to be an intriguing triplet ground state, in which the unpaired density is mostly located at three Li ligands. Chemical bonding analysis reveals that the 2p_zπ electrons of C₃Li₃^- are fully located at the C₃ ring formed by C-C multiple bonds, in contrast to the perfect 2p_zπ-delocalization found in the well-known ppCs.

A sixteen-valence-electron carbon-group 13 family with global penta-atomic planar tetracoordinate carbon: an ionic strategy

The design of planar tetracoordinate carbon (ptC) has always been a challenge due to its unique bonding mode that necessitates the perfect balance between the carbon center and surrounding ligands both electronically and mechanically. A unique type of 18-valence-electron (18ve) template, i.e., CAl42-, has been found to be very effective in designing various novel 18ve-species upon skeletal substitution. In this work, we showed that though ptC is not the global structure for the parent 16ve-CAl4, suitable skeletal substitution can allow for a series of global minimum ptC species. Theoretical calculations at the level of CCSD(T)/def2-QZVP//B3LYP/def2-QZVP for 35 carbon-group 13 systems with 16-ve, i.e., CXaYbZcKd (X, Y, Z, K = Al/Ga/In/Tl; 0 ≤ a, b, c, d ≤ 4, a + b + c + d = 4), showed that 9 systems (CAl3Tl, CGa3Tl, CGa2Tl2, CAl2GaTl, CAl2InTl, CGa2InTl, CAlGa2Tl, CGa2InTl and CAlGaInTl) possess global minimum ptC and 2 systems (CAl3In and CAl2Tl2) have quasi-GM ptC. Except for CAl3Tl and CAl3In, all the ptCs were predicted for the first time. All these stable ptC structures have the same skeleton and can be described as the same ionic sub-structure, i.e., [A-]B+. This study not only enriches 16ve-ptC, but also directly demonstrates that utilizing an ionic strategy, non-ptC CAl4 also can be used as a template to extend the ptC family.

Mono-silicon isoelectronic replacement in CAl4 : van't hoff/le bel carbon or not?

In cluster studies, the isoelectronic replacement strategy has been successfully used to introduce new elements into a known structure while maintaining the desired topology. The well-known penta-atomic 18 valence electron (ve) species C Al 4 2 - and its Al^- /Si or Al/Si⁺ isoelectronically replaced clusters CAl₃ Si^- , CAl₂ Si₂ , C AlSi 3 - , and C Si 4 2 + , all possess the same anti-van't Hoff/Le Bel skeletons, that is, nontraditional planar tetracoordinate carbon (ptC) structure. In this article, however, we found that such isoelectronic replacement between Si and Al does not work for the 16ve-CAl₄ with the traditional van't Hoff/Le Bel tetrahedral carbon (thC) and its isoelectronic derivatives CAl₃ X (X = Ga/In/Tl). At the level of CCSD(T)/def2-QZVP//B3LYP/def2-QZVP, none of the global minima of the 16ve mono-Si-containing clusters CAl₂ SiX⁺ (X = Al/Ga/In/Tl) maintains thC as the parent CAl₄ does. Instead, X = Al/Ga globally favors an unusual ptC structure that has one long C─X distance yet with significant bond index value, and X = In/Tl prefers the planar tricoordinate carbon. The frustrated formation of thC in these clusters is ascribed to the CSi bonding that prefers a planar fashion. Inclusion of chloride ion would further stabilize the ptC of CAl₂ SiAl⁺ and CAl₂ SiGa⁺ . The unexpectedly disclosed CAl₂ SiAl⁺ and CAl₂ SiGa⁺ represent the first type of 16ve-cationic ptCs with multiple bonds. © 2019 Wiley Periodicals, Inc.

CAl3X (X = B/Al/Ga/In/Tl) with 16 valence electrons: can planar tetracoordinate carbon be stable?

As a perpetual chemical curiosity, planar tetracoordinate carbon (ptC) that violates the traditional tetrahedral carbon (thC) has made enormous achievements. In particular, the 18-valence-electron (18ve) counting rule has been found to be very effective in predicting ptC structures, as in CX₄^2- (X = Al/Ga/In/Tl). By contrast, the corresponding neutral CX₄ with 16ve each takes the thC form like methane. Herein, we report a mono-substituted neutral 16ve-CAl₃X (X = Al/Ga/In/Tl). Our theoretical results showed that the competition between thC and ptC can be well tuned upon variation of X, and for X = In and Tl, the ptC structure becomes isoenergetic to and even more stable than thC, respectively. Thus, a low-lying ptC can be achieved in the 16ve-CAl₃X set without acquiring additional electrons. This unintuitive result can be ascribed to the increased energetic preference of the ionic sub-structure [CAl₃^-]X⁺ from X = Al to Tl. We thus predict the first penta-atomic ptC species with 16ve, and the ionic strategy presented in this work is expected to promote novel designs of ptC molecules.

Planar pentacoordinate carbon in CGa5+ derivatives

A single isoelectronic substitution by a smaller sized element, Be, is enough to stabilize a CGa₅⁺ cluster with a planar pentacoordinate carbon at the centre.

A template for a planar tetracoordinate heavier group 14 atom: a global study of C2Si2X(q) (X = C, Si, Ge, Sn, Pb; q = +1, 0, -1)

Through a global isomeric study, we computationally identified the first structural template C2Si2X that could encompass a planar tetracoordinate X for all the heavier group 14 elements X in the 0, +1 or -1 charge state. We thus significantly expanded the traditional 16/17/18ve rules to 19/20/21ve for ptX.