Analysis of bonding motifs in unusual molecules I: planar hexacoordinated carbon atoms

The bonding structures of CO3Li3+ and CS3Li3+ are studied by means of oriented quasi-atomic orbitals (QUAOs) to assess the possibility of these molecules being planar hexacoordinated carbon (phC) systems. CH3Li and CO32- are employed as reference molecules. It is found that the introduction of Li+ ions into the molecular environment of carbonate has a greater effect on the orbital structure of the O atoms than it does on the C atom. Partial charges computed from QUAO populations imply repulsion between the positively charged C and Li atoms in CO3Li3+. Upon the transition from CO3Li3+ to CS3Li3+, the analysis reveals that the substitution of O atoms by S atoms inverts the polarity of the carbon-chalcogen σ bond. This is linked to the difference in s- and p-fractions of the QUAOs of C and S, as element electronegativities do not explain the observed polarity of the CSσ bond. Partial charges indicate that the larger electron population on the C atom in CS3Li3+ makes C-Li attraction possible. Upon comparison with the C-Li bond in methyllithium, it is found that the C-Li covalent interactions in CO3Li3+ and CS3Li3+ have about 14% and 6% of the strength of the C-Li covalent interaction in CH3Li, respectively. Consequently, it is concluded that only CS3Li3+ may be considered to be a phC system.

InnTl4-nH+ (n = 0∼4): Tetracoordinate Hydrogen in a Planar Fashion?

The recent report of planar tetracoordinate hydrogen (ptH) in In4H+ is very intriguing in planar hypercoordinate chemistry. Our high-level CCSD(T) calculations revealed that the proposed D4h-symmetric ptH In4H+ is a first-order saddle point with an imaginary frequency in the out-of-plane mode of the hydrogen atom. In fact, at the CCSD(T)/aug-cc-pV5Z/aug-cc-pV5Z-PP level, the C4v isomer, with the H atom located 0.70 Å above the In4 plane, is 0.5 kcal/mol more stable than the D4h isomer. However, given the small perturbation from planarity and essentially barrierless C4v ↔ D4h ↔ C4v transition, the vibrationally averaged structure can still be considered as a planar. Extending our exploration to the InnTl4-nH+ (n = 0-3) systems, we found all these ptH structures, except for In2Tl2H+, to be the putative global minimum. The single σ-delocalized interaction between the central hydrogen atom and InnTl4-n ligand rings proves pivotal in establishing planarity and aromaticity and conferring substantial stability upon these rule-breaking ptH species.

Planar Hexacoordinate Beryllium: Covalent Bonding Between s-block Metals

Achieving a planar hypercoordinate arrangement of s-block metals through covalent bonding with ligands is challenging due to the strong ionicity involved. Herein, we report the first case of a neutral binary global minimum containing a planar hexacoordinate beryllium atom. The central Be atom is coordinated by six active Be atoms, the latter in turn are enclosed by an equal number of more electronegative chlorine atoms in the periphery, forming a star-like phBe cluster (Be©Be6 Cl6 ). Importantly, the cluster exhibits dynamically stabilized stemming geometrically from the appropriate matching of metal-ligand size and electronically from adherence to the octet rule as well as possessing a 6σ/2π double aromaticity. Remarkably, energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) analysis reveals a significant covalent interaction between the ligand and the central metal beryllium atoms, a fact further supported by a large Wiberg bond index. This cluster is a promising synthetic as its excellent electronic, dynamic and thermodynamic stability.

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.

Searching for Systems with Planar Hexacoordinate Carbons

Here, we present evidence that the D2h M2C50/2+ (M = Li-K, Be-Ca, Al-In, and Zn) species comprises planar hexacoordinate carbon (phC) structures that exhibit four covalent and two electrostatic interactions. These findings have been made possible using evolutionary methods for exploring the potential energy surface (AUTOMATON program) and the Interacting Quantum Atoms (IQA) methodology, which support the observed bonding interactions. It is worth noting, however, that these structures are not the global minimum. Nonetheless, incorporating two cyclopentadienyl anion ligands (Cp) into the CaC52+ system has enhanced the relative stability of the phC isomer. Moreover, cycloparaphenylene ([8]CPP) provides system protection and kinetic stability. These results indicate that using appropriate ligands presents a promising approach for expanding the chemistry of phC species.

σ-Aromatic MAl6S6 (M = Ni, Pd, Pt) Stars Containing Planar Hexacoordinate Transition Metals

Hypercoordinate transition-metal species are mainly dominated by the 18-valence-electron (18ve) counting. Herein, we report ternary MAl₆S₆ (M = Ni, Pd, Pt) clusters with the planar hexacoordinate metal (phM) centers, which feature 16ve counting instead of the classic 18ve rule. These global-minimum clusters are established via unbiased global searches, followed by PBE0 and single-point CCSD(T) calculations. The phM MAl₆ units are stabilized by six peripheral bridging S atoms in these star-like species. Chemical bonding analyses reveal that there are 10 delocalized electrons around the phM center, which can render the aromaticity according to the (4n + 2) Hückel rule. It is worth noting that adding an (or two) electron(s) to its π-type lowest unoccupied molecular orbital (LUMO) will make the system unstable.

Two-dimensional Be2Al and Be2Ga monolayer: anti-van't Hoff/Le Bel planar hexacoordinate bonding and superconductivity

Because of the electron deficiency of boron, a triangular network with planar hexacoordination is the most common structural and bonding property for isolated boron clusters and two-dimensional (2D) boron sheets. However, this network is a rule-breaking structure and bonding case for all other main-group elements. Herein, the Be₂M (M = Al and Ga) 2D monolayer with P6/mmm space group was found to be the lowest-energy structure with planar hexacoordinate Be/Al/Ga motifs. More interestingly, Be₂Al and Be₂Ga were observed to be intrinsic phonon-mediated superconductors with a superconducting critical temperature (T_c) of 5.9 and 3.6 K, respectively, where compressive strain could further enhance their T_c. The high thermochemical and kinetic stability of Be₂M make a promising candidate for experimental realization, considering its high cohesive energy, absence of soft phonon modes, and good resistance to high temperature. Moreover, the feasibility of directly growing Be₂M on the electride Ca₂N substrate was further demonstrated, where its intriguing electronic and superconducting properties were well maintained in comparison with the freestanding monolayer. The Be₂M monolayer with rule-breaking planar hexacoordinate motifs firmly pushes the ultimate connection of the "anti-van't Hoff/Le Bel" structure with promising physical properties.

Structure and Bonding in Planar Hypercoordinate Carbon Compounds

The term hypercoordination refers to the extent of the coordination of an element by its normal value. In the hypercoordination sphere, the element can achieve planar and/or non-planar molecular shape. Hence, planar hypercoordinate carbon species violate two structural rules: (i) The highest coordination number of carbon is four and (ii) the tetrahedral orientation by the connected elements and/or groups. The unusual planar orientations are mostly stabilized by the electronic interactions of the central atom with the surrounding ligands. In this review article, we will talk about the current progress in the theoretical prediction of viable planar hypercoordinate carbon compounds. Primary knowledge of the planar hypercoordinate chemistry will lead to its forthcoming expansion. Experimental and theoretical interests in planar tetracoordinate carbon (ptC), planar pentacoordinate carbon (ppC), and planar hexacoordinate carbon (phC) are continued. The proposed electronic and mechanical strategies are helpful for the designing of the ptC compounds. Moreover, the 18-valence electron rule can guide the design of new ptC clusters computationally as well as experimentally. However, the counting of 18-valence electrons is not a requisite condition to contain a ptC in a cluster. Furthermore, this ptC idea is expanded to the probability of a greater coordination number of carbon in planar orientations. Unfortunately, until now, there are no such logical approaches to designing ppC, phC, or higher-coordinate carbon molecules/ions. There exist a few global minimum structures of phC clusters identified computationally, but none have been detected experimentally. All planar hypercoordinate carbon species in the global minima may be feasible in the gas phase.

CB6Al0/+: Planar hexacoordinate boron (phB) in the global minimum structure

Herein, we report for the first time the presence of a planar hexacoordinate boron (phB) atom in the global minimum energy structure of a neutral cluster system. The potential energy surface (PES) has been explored for CB₆Al^0/+/- systems using density functional theory (DFT). The global minima of CB₆Al (1a) and CB₆Al⁺ (1b) contain a phB center. However, the global minimum of CB₆Al^- (1c) does not have a phB atom. The CCSD(T)/aug-cc-pVTZ level of theory has been applied to compute the relative energies of the low-lying isomers with respect to the 1a and 1b structures of CB₆Al and CB₆Al⁺ systems, respectively. The exploration of the PES of CB₆^0/+/- systems indicates that the global minima do not contain a phB atom. However, the incorporation of an aluminium (Al) atom into the CB₆ moiety produces structures containing a phB center in the CB₆Al^0/+ systems. Hence, the Al metal has an important role in attaining a planar geometry having a hexacoordinate boron center. The dynamical stability of CB₆Al (1a) and CB₆Al⁺ (1b) was confirmed from the atom-centered density matrix propagation (ADMP) simulation over 20 ps of time at temperatures of 300 K and 400 K. The natural charge computations showed that the charges on the phB are almost zero in both systems. The 1a structure has σ/π-dual aromaticity as predicted from the nucleus independent chemical shift (NICS) values and the gauge-including magnetically induced ring current (GIMIC).

Infrared intensities of [Formula: see text]: a true challenge for DFT methods

Absolute infrared intensities of [Formula: see text] were evaluated with a great variety of DFT and ab initio methods and basis sets. It is shown that the intensities calculated by different levels of theory may not agree with each other even in the qualitative (weak/strong) sense. Geometrical parameters, as well as net atomic charges evaluated from multiple partition schemes, did not vary as much as the intensities and thus cannot explain the tremendous differences found for the latter. As there are no experimental estimates for the intensities to guide the theoretical evaluation, deciding on the best level of theory is reduced to comparisons between the different DFT methods and QCISD or CCSD, believed to be the best theoretical estimates in the set. The differences found among the various DFT methods suggest the development of new methods, instead of converging to a focal point, is rather diverging.

Ternary CE2Ba2 (E = As, Sb) Clusters: New Pentaatomic Planar Tetracoordinate Carbon Species with 18 Valence Electrons

18-valence-electron (ve) rule is one important guide for us to design planar tetracoordinate carbon (ptC) species. Using the "polarization of ligands" strategy, the new pentaatomic ptC species CE₂Ba₂ (E = As, Sb) with 18 ve are designed in this work. Computer structural searches and high-level calculations reveal that the ptC CE₂Ba₂ (E = As, Sb) species are global minima (GMs) on the potential energy surfaces, whose C center is coordinated by the interspaced E and Ba atoms. CE₂Ba₂ (E = As, Sb) are also kinetically stable. Chemical bonding analyses reveal that the ptC core is stabilized by two localized C-E σ bonds, one delocalized five-center two-electron (5c-2e) σ bond and one delocalized 5c-2e π bond. One π and three σ bonds collectively conform to the 8-electron counting, which determines the stability of ptC CE₂Ba₂ (E = As, Sb) species. Interestingly, the delocalized 2π and 2σ electrons render the ptC systems π/σ double aromaticity. Additional 10 electrons contribute to peripheral lone pairs of E and E-Ba bonding.

Persistent Planar Tetracoordinate Carbon in Global Minima Structures of Silicon-Carbon Clusters

Recently, we reported a series of global minima whose structures consist of carbon rings decorated with heavier group 14 elements. Interestingly, these structures feature planar tetracoordinate carbons (ptCs) and result from the replacement of five or six protons (H+) from the cyclopentadienyl anion (C5H5−) or the pentalene dianion (C8H62−) by three or four E2+ dications (E = Si–Pb), respectively. The silicon derivatives of these series are the Si3C5 and Si4C8 clusters. Here we show that ptC persists in some clusters with an equivalent number of C and Si atoms, i.e., Si5C5, Si8C8, and Si9C9. In all these species, the ptC is embedded in a pentagonal C5 ring and participates in a three-center, two-electron (3c-2e) Si-ptC-Si σ-bond. Furthermore, these clusters are π-aromatic species according to chemical bonding analysis and magnetic criteria.

Two-dimensional Be2Au Monolayer with Planar Hexacoordinate s-Block Metal Atoms: A Superconducting Global Minimum Dirac Material with Two Perfect Dirac Node-Loops

Using a starlike Be₆Au₇⁻ cluster as a building block and following the bottom-up strategy, an intriguing two-dimensional (2D) binary s-block metal Be₂Au monolayer with a P6/mmm space group was theoretically designed. Both the Be₆Au₇⁻ cluster and the 2D monolayer are global minima featuring rule-breaking planar hexacoordinate motifs (anti-van't Hoff/Le Bel arrangement), and their high stabilities are attributed to good electron delocalization and electronic-stabilization-induced steric force. Strikingly, the Be₂Au monolayer is a rare Dirac material with two perfect Dirac node-loops in the band structure and is a phonon-mediated superconductor with a critical temperature of 4.0 K. The critical temperature can be enhanced up to 11.0 K by applying compressive strain at only 1.6%. This study not only identifies a new binary s-block metal 2D material, namely Be₂Au, which features planar hexacoordination, and a candidate superconducting material for further explorations, but also provides a new strategy to construct 2D materials with novel chemical bonding.

A topological superconductor, named Be₂Au monolayer, containing planar hexacoordinate s-block metal (Be and Au) atoms was theoretically designed by rationally assembling related clusters.

Bare and ligand protected planar hexacoordinate silicon in SiSb3M3+ (M = Ca, Sr, Ba) clusters

The global minimum of SiSb3M3+ (M = Ca, Sr, Ba) is a D3h symmetric structure containing an elusive planar hexacoordinate silicon (phSi) atom. Most importantly, the phSi core remains intact in ligand protected environment as well.

Planar Hypercoordinate Carbons in Alkali Metal Decorated CE3 2- and CE2 2- Dianions

After exploring the potential energy surfaces of M_m CE₂ ^p (E=S-Te, M=Li-Cs, m=2, 3 and p=m-2) and M_n CE₃ ^q (E=S-Te, M=Li-Cs, n=1, 2, q=n-2) combinations, we introduce 38 new global minima containing a planar hypercoordinate carbon atom (24 with a planar tetracoordinate carbon and 14 with a planar pentacoordinate carbon). These exotic clusters result from the decoration of V-shaped CE₂ ^2- and Y-shaped CE₃ ^2- dianions, respectively, with alkali counterions. All these 38 systems fulfill the geometrical and electronic criteria to be considered as true planar hypercoordinate carbon systems. Chemical bonding analyses indicate that carbon is covalently bonded to chalcogens and ionically connected to alkali metals.

Planar hexacoordinate gallium

We report the first planar hexacoordinate gallium (phGa) center in the global minimum of the GaBe6Au6 + cluster which has a star-like D 6h geometry with 1A1g electronic state, possessing a central gallium atom encompassed by a Be6 hexagon and each Be-Be edge is further capped by an Au atom. The electronic delocalization resulting in double aromaticity (both σ and π) provides electronic stability in the planar form of the GaBe6Au6 + cluster. The high kinetic stability of the title cluster is also understood by Born-Oppenheimer molecular dynamics simulations. The energy decomposition analysis in combination with the 'natural orbitals for chemical valence' theory reveals that the bonding in the GaBe6Au6 + cluster is best expressed as the doublet Ga atom with 4s24p⊥ 1 electronic configuration forming an electron-sharing π bond with the doublet Be6Au6 + moiety followed by Ga(s)→[Be6Au6 +] σ-backdonation and two sets of Ga(p‖)←[Be6Au6 +] σ-donations.

σ-Aromaticity Planar Pentacoordinate Beryllium Atoms

Six-valence-electron planar pentacoordinate beryllium (ppBe) is explored herein as a global minimum, which is only constructed by s-block metals in BeM₅⁺ (M = Cu, Ag, Au). The bonding in ppBe can be regarded as the excited-stated Be with a 2p_x¹2p_y¹ electronic configuration, forming electron sharing with doublet M₅⁺ motifs followed by two sets of Be(p_∥) → [M₅⁺] σ donations and one Be(s) ← [M₅⁺] σ back-donation. Thus, the σ aromaticity originating from three delocalized σ orbitals gives rise to the whole stability of the high D_5h-symmetry ppBe and strongly enriches s-block planar hypercoordinate bonding.

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.

Analysis of Local and Global Aromaticity in Si3C5 and Si4C8 Clusters. Aromatic Species Containing Planar Tetracoordinate Carbon

The minimum energy structures of the Si3C5 and Si4C8 clusters are planar and contain planar tetracoordinate carbons (ptCs). These species have been classified, qualitatively, as global (π) and local (σ) aromatics according to the adaptive natural density partitioning (AdNDP) method, which is an orbital localization method. This work evaluates these species’ aromaticity, focusing on confirming and quantifying their global and local aromatic character. For this purpose, we use an orbital localization method based on the partitioning of the molecular space according to the topology of the electronic localization function (LOC-ELF). In addition, the magnetically induced current density is analyzed. The LOC-ELF-based analysis coincides with the AdNDP study (double aromaticity, global, and local). Moreover, the current density analysis detects global and local ring currents. The strength of the global and local current circuit is significant, involving 4n + 2 π- and σ-electrons, respectively. The latter implicates the Si-ptC-Si fragment, which would be related to the 3c-2e σ-bond detected by the orbital localization methods in this fragment.

Planar Hexacoordinate Carbons: Half Covalent, Half Ionic

Herein, the first global minima containing a planar hexacoordinate carbon (phC) atom are reported. The fifteen structures belong to the CE₃ M₃ ⁺ (E=S-Te and M=Li-Cs) series and satisfy both geometric and electronic criteria to be considered as a true phC. The design strategy consisted of replacing oxygen in the D_3h  CO₃ Li₃ ⁺ structure with heavy and less electronegative chalcogens, inducing a negative charge on the C atom and an attractive electrostatic interaction between C and the alkali-metal cations. The chemical bonding analyses indicate that carbon is covalently bonded to three chalcogens and ionically connected to the three alkali metals.

CAl4Mg0/−: Global Minima with a Planar Tetracoordinate Carbon Atom

Isomers of CAl4Mg and CAl4Mg− have been theoretically characterized for the first time. The most stable isomer for both the neutral and anion contain a planar tetracoordinate carbon (ptC) atom. Unlike the isovalent CAl4Be case, which contains a planar pentacoordinate carbon atom as the global minimum geometry, replacing beryllium with magnesium makes the ptC isomer the global minimum due to increased ionic radii of magnesium. However, it is relatively easier to conduct experimental studies for CAl4Mg0/− as beryllium is toxic. While the neutral molecule containing the ptC atom follows the 18 valence electron rule, the anion breaks the rule with 19 valence electrons. The electron affinity of CAl4Mg is in the range of 1.96–2.05 eV. Both the global minima exhibit π/σ double aromaticity. Ab initio molecular dynamics simulations were carried out for both the global minima at 298 K for 10 ps to confirm their kinetic stability.

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.

Computational design of species with ultrashort Be–Be distances using planar hexacoordinate carbon structures as the templates

Replacing the planar hexacoordinate carbon in CX₃M₃⁺ species with the Be₂ moiety leads to isoelectronic species with ultrashort Be–Be distances.

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.

Viable aromatic BenHn stars enclosing a planar hypercoordinate boron or late transition metal

Similar to B_n rings, star-like Be_nH_n rings can serve as the n-electron σ-donors for designing species with planar hypercoordinate atom.

CB3E2q (q = ±1): a family of “hyparene” analogues with a planar pentacoordinate carbon

Highlighted elements are the feasible E atoms of CB₃E₂^q template ppC structures, in which CB₃Mg₂⁻ is suitable for photoelectron spectroscopy.

A two-dimensional TiB4 monolayer exhibits planar octacoordinate Ti

At present, the concept of planar hypercoordination in chemistry meets the fast development of two-dimensional (2D) nanomaterials, leading to considerable interest in searching for 2D materials with planar hypercoordinate atoms. In this work, by means of the swarm-intelligence structure search method and first-principles calculations, we predict a hitherto unknown 2D TiB₄ monolayer with a planar octacoordinate Ti moiety, in which each Ti atom binds to eight B atoms with equal distances in a perfect plane, and has the highest coordination of Ti known for 2D materials thus far. Systematic ab initio calculations demonstrate the superior thermodynamic and dynamic stabilities of the predicted TiB₄ monolayer, indicating the high feasibility for experimental synthesis. The stabilization of this perfect planar structure originates from the geometric fit between the atomic radius of Ti and the size of the 8-membered B ring, as well as the electron transfer from Ti to B atoms which compensates the electron deficiency of the full sp² hybridized B network. Motivated by the unforeseen geometry of the TiB₄ monolayer, a series of other 2D transition metal borides (TMB₄, TM = V, Cr, Mo, W and Os) with quasi-planar octacoordinate TM atoms are further designed and discussed. The present work provides a useful roadmap for the discovery of 2D hypercoordinate materials.

Prediction of a new BeC monolayer with perfectly planar tetracoordinate carbons

Recently, there has been a growing interest in exploring planar hypercoordinate carbons in two-dimensional nanostructures. However, atomic monolayers with ideal planar hypercoordinate carbon are quite rare due to the challenge in stabilizing the exotic motifs. We predicted a global minimum two-dimensional BeC monolayer using the global particle-swarm optimization method. Each carbon binds peripheral four atoms in the same plane, forming a perfectly planar tetracoordinate carbon moiety. The cohesive energy, phonon-spectrum and mechanical stability criteria confirm the stability of monolayer BeC. In addition, the BeC monolayer has a large in-plane stiffness (145.54 N m^-1) and thermo-dynamical stability (up to 2000 K). Furthermore, BeC is an indirect semiconductor with a band gap of 1.01 eV and possesses exceptionally high carrier mobilities (∼10⁴ cm² V^-1 s^-1), rendering it suitable for application in electronics and photoelectronics.

Testing the effectiveness of the isoelectronic substitution principle through the transformation of aromatic osmathiophene derivatives into their inorganic analogues

The isoelectronic substitution principle has been employed in the [OsCl₂(SC₃H₃)(PH₃)₂]⁺ complex, generating a new set of complexes, where C has been replaced by B⁻ and/or N⁺.

Linear, planar, and tubular molecular structures constructed by double planar tetracoordinate carbon D2hC2(BeH)4 species via hydrogen-bridged -BeH2Be- bonds

This computational study identifies the rhombic D2hC2 (BeH)4 (2a) to be a species featuring double planar tetracoordinate carbons (ptCs). Aromaticity and the peripheral BeBeBeBe bonding around CC core contribute to the stabilization of the ptC structure. Although the ptC structure is not a global minimum, its high kinetic stability and its distinct feature of having a bonded C2 core from having two separated carbon atoms in the global minimum and other low-lying minima could make the ptC structure to be preferred if the carbon source is dominated by C2 species. The electron deficiency of the BeH group allows the ptC species to serve as building blocks to construct large/nanostructures, such as linear chains, planar sheets, and tubes, via intermolecular hydrogen-bridged bonds (HBBs). Formation of one HBB bond releases more than 30.0 kcal/mol of energy, implying the highly exothermic formation processes and the possibility to synthesize these nano-size structures.

Computational design of organometallic oligomers featuring 1,3-metal-carbon bonding and planar tetracoordinate carbon atoms

Density functional theory computations (B3LYP) have been used to explore the chemistry of titanium-aromatic carbon "edge complexes" with 1,3-metal-carbon (1,3-MC) bonding between Ti and planar tetracoordinate Cβ . The titanium-coordinated, end-capping chlorides are replaced with OH or SH groups to afford two series of difunctional monomers that can undergo condensation to form oxide- and sulfide-bridged oligomers. The sulfide-linked oligomers have less molecular strain and are more exergonic than the corresponding oxide-linked oligomers. The HOMO-LUMO gap of the oligomers varies with their composition and decreases with growing oligomer chain. This theoretical study is intended to enrich 1,3-MC bonding and planar tetracoordinate carbon chemistry and provide interesting ideas to experimentalists. Organometallic complexes with the TiE2 (E = OH and SH) decoration on the edge of aromatic hydrocarbons have been computationally designed, which feature 1,3-metal-carbon (1,3-MC) bonding between titanium and planar tetracoordinate β-carbon. Condensation of these difunctional monomers by eliminating small molecules (H2O and H2S) produce chain-like oligomers. The HOMO-LUMO gaps of the oligomers decreases with growing oligomer chain, a trend that suggests possible semiconductor properties for oligomers with longer chains.

Planar tetracoordinate carbons with a double bond in CAl3E clusters

The potential energy surfaces of a series of clusters with the formula CAl3E (E = P, As, Sb, Bi) are systematically explored using density functional theory and high level ab initio calculations. The global minimum structure of these clusters contains a planar tetracoordinate carbon atom. The presence of a C=E double bond is supported by the Wiberg bond indices, the adaptive natural density partitioning analysis, and the magnetic response. Our results show that these planar tetracoordinate carbon clusters are both thermodynamically and kinetically viable species in the gas phase.

Exploring the geometrical structures of X©BnHnm [(X, m) = (B, +1), (C, +2) for n = 5; (X, m) = (Be, 0), (B, +1) for n = 6] by an electronic method

The (quasi-)planar wheel-type structures can be obtained by adding electrons.

Planar pentacoordinate carbons in CBe54− derivatives

Ab initio computations show that the global minimum structure of the CBe₅Li_nⁿ⁻⁴ clusters (n = 1 to 5) contains a planar pentacoordinate carbon atom.

D(3h) [A-CE3-A]⁻ (E = Al and Ga, A = Si, Ge, Sn, and Pb): a new class of hexatomic mono-anionic species with trigonal bipyramidal carbon

The non-classical trigonal bipyramidal carbon (TBPC) arrangement generally exists as transition states (TSs) in nucleophilic bimolecular substitution (S(N)2) reactions. Nevertheless, chemists have been curious about whether such a carbon bonding could be stable in equilibrium structures for decades. As the TBPC arrangement was normally realized as cationic species theoretically and experimentally, only one anionic example ([At-C(CN)3-At](-)) was computationally devised. Herein, we report the design of a new class of anionic TBPC species by using the strategy similar to that for stabilizing the non-classical planar hypercoordinate carbon. When electron deficient Al and Ga were used as the equatorial ligands, eight D(3h) [A-CE3-A](-) (E = Al and Ga, A = Si, Ge, Sn, and Pb) TBPC structures were found to be the energy minima rather than TSs at both the B3LYP and MP2 levels. Remarkably, the energetic results at the CCSD(T) optimization level further identify [Ge-CAl3-Ge](-) and [Sn-CGa3-Sn](-) even to be the global minima and [Si-CAl3-Si](-) and [Ge-CGa3-Ge](-) to be the local minima, only slightly higher than their global minima. The electronic structure analyses reveal that the substantial ionic C-E bonding, the peripheral E-A covalent bonding, and the axial mc-2e (multi center-two electrons) bonding play roles in stabilizing these TBPC structures. The structural simplicity and the high thermodynamic stability suggest that some of these species may be generated and captured in the gas phase. Furthermore, as mono-anionic species, their first vertical detachment energies are differentiable from those of their nearest isomers, which would facilitate their characterization via experiments such as the negative ion photoelectron spectroscopy.