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.

Global minimum and a heap of low-lying isomers with planar tetracoordinate carbon in the CAl3MgH2- system

A global minimum and a heap of low-lying isomers with planar tetracoordinate carbon (ptC) are identified in the CAl3MgH2- system by computational quantum chemical investigations. The nature of the chemical bonding in the global minimum ptC isomer is examined using the conceptual quantum chemical tools. The atoms in molecule (AIM) analysis reveals that the global minimum isomer possesses a ptC geometry. Additionally, the adaptive natural density partitioning (AdNDP), electron localization function (ELF), and nucleus-independent chemical shifts (NICS) analysis corroborate the presence of delocalization in the ptC isomer. The delocalization of electron density in the global minimum ptC isomer contributes to attaining structural stability. The results also suggest that the bridging hydrogen plays a crucial role in stabilizing the ptC system. Furthermore, the ab initio molecular dynamics study supports the structural stability of the ptC isomer.

Reconsidering the Structures of C2 Al4 - and C2 Al5. Chemistry 2023;29:e202301338. [PMID: 37498677 DOI: 10.1002/chem.202301338]

The study of C2 Al4 -/0 and C2 Al5 -/0 was conducted using anion photoelectron spectroscopy and quantum chemical computations. The present findings reveal that C2 Al4 - has a boat-like structure, with a single C2 unit surrounded by four aluminum atoms. In contrast, the neutral C2 Al4 species adopts a D2h planar structure with two planar tetracoordinate carbon (ptC) units, consistent with previous reports. Furthermore, the global minimum isomer of C2 Al5 - adopts a D3h symmetry, where the C2 unit interacts with five aluminum atoms. It was also found that a lower symmetry structure of C2 Al5 - , where all five aluminum atoms are located on the same side of the C2 unit, albeit slightly higher in energy compared to the D3h structure. These computations show that the D3h structure of C2 Al5 - is highly stable, exhibiting a large HOMO-LUMO gap.

3-D molecular stars with covalent axial bonding

In designing three-dimensional (3-D) molecular stars, it is very difficult to enhance the molecular rigidity through forming the covalent bonds between the axial and equatorial groups because corresponding axial groups will generally break the delocalized π bond over equatorial frameworks and thus break their star-like arrangement. In this work, exemplified by designing the 3-D stars Be₂ ©Be₅ E₅ ⁺ (E = Au, Cl, Br, I) with three delocalized σ bonds and delocalized π bond over the central Be₂ ©Be₅ moiety, we propose that the desired covalent bonding can be achieved by forming the delocalized σ bond(s) and delocalized π bond(s) simultaneously between the axial groups and equatorial framework. The covalency and rigidity of axial bonding can be demonstrated by the total Wiberg bond indices of 1.46-1.65 for axial Be atoms and ultrashort Be-Be distances of 1.834-1.841 Å, respectively. Beneficial also from the σ and π double aromaticity, these mono-cationic 3-D molecular stars are dynamically viable global energy minima with well-defined electronic structures, as reflected by wide HOMO-LUMO gaps (4.68-5.06 eV) and low electron affinities (4.70-4.82 eV), so they are the promising targets in the gas phase generation, mass-separation, and spectroscopic characterization.

Pentacoordinate Carbon Atoms in a Ferrocene Dication Derivative—[Fe(Si2-η5-C5H2)2]2+

Pentacoordinate carbon atoms are theoretically predicted here in a ferrocene dication derivative in the eclipsed-(1; C2v), gauche-(2; C2) and staggered-[Fe(Si2-η5-C5H2)2]2+(3; C2h) forms for the first time. Energetically, the relative energy gaps for 2 and 3 range from −3.06 to 16.74 and −2.78 to 40.34 kJ mol−1, respectively, when compared to the singlet electronic state of 1 at different levels. The planar tetracoordinate carbon (ptC) atom in the ligand Si2C5H2 becomes a pentacoordinate carbon upon complexation. The ligand with a ptC atom was predicted to be both a thermodynamically and kinetically stable molecule by some of us in our earlier theoretical works. Natural bond orbital and adaptive natural density partitioning analyses confirm the pentacoordinate nature of carbon in these three complexes (1–3). Although they are hypothetical at the moment, they support the idea of “hypercoordinate metallocenes” within organometallic chemistry. Moreover, ab initio molecular dynamics simulations carried out at 298 K temperature for 2000 fs suggest that these molecules are kinetically stable.

Structural Evolution of Carbon-Doped Aluminum Clusters AlnC- (n = 6-15): Anion Photoelectron Spectroscopy and Theoretical Calculations

Carbon-doped aluminum cluster anions, Al_nC^- (n = 6-15), were generated by laser vaporization and investigated by mass-selected anion photoelectron spectroscopy. The geometric structures of Al_nC^- (n = 6-15) anions were determined by the comparison of theoretical calculations with the experimental results. It is found that the most stable structure of Al₆C^- is a carbon endohedral triangular prism. The Al₇C^- anion is a magic cluster with high stability. The structures of Al_7-9C^- can be viewed as the additional aluminum atoms attached around the triangular prism Al₆C^-. Two isomers of Al₁₀C^- have been detected in the experiments. The most stable one has a planar tetracoordinate carbon structure. The second one derives from Al₉C^- with the carbon atom located in a pentagonal bipyramid. The Al₁₁C^- anion has a bilayer structure composed of one planar tetracoordinate carbon and one aluminum-centered hexagon, in which the major interactions between two layers are multicenter bonds. The structures of Al_12-14C^- can be viewed as evolving from Al₁₁C^- by adding aluminum atoms to interact with the carbon atom. In Al₁₅C^-, the carbon atom stays at the surface with a tetracoordinate structure, and an icosahedral Al₁₃ unit can be identified as a part of the geometric structure of Al₁₅C^-.

Binary mono-anions with unprecedented anti-aromatic planar tetracoordinate carbon and nitrogen atoms

Global minima (CBe₄^- and NBe₄^-) with planar tetracoordinate carbon and nitrogen atoms are stabilized by multicentric bonds. Their unusual anti-aromaticity is due to the strong localized Be-Be bonds preventing the full delocalization of the electrons in the CBe₄/NBe₄ skeleton. These binary mono-anions are expected to be experimentally realized in the gas phase.

From Molecules with a Planar Tetracoordinate Carbon to an Astronomically Known C5H2 Carbene

Ethynylcyclopropenylidene (2), an isomer of C₅H₂, is a known molecule in the laboratory and has recently been identified in Taurus Molecular Cloud-1 (TMC-1). Using high-level coupled-cluster methods up to the CCSDT(Q)/CBS level of theory, it is shown that two isomers of C₅H₂ with a planar tetracoordinate carbon (ptC) atom, (SP-4)-spiro[2.2]pent-1,4-dien-1,4-diyl (11) and (SP-4)-spiro[2.2]pent-1,4-dien-1,5-diyl (13), serve as the reactive intermediates for the formation of 2. Here, a theoretical connection has been established between molecules containing ptC atoms (11 and 13) and a molecule (2) that is present nearly 430 light years away, thus providing evidence for the existence of ptC species in the interstellar medium. The reaction pathways connecting the transition states and the reactants and products have been confirmed by intrinsic reaction coordinate calculations at the CCSDT(Q)/CBS//B3LYP-D3BJ/cc-pVTZ level. While isomer 11 is non-polar (μ = 0), isomers 2 and 13 are polar, with dipole moment values of 3.52 and 5.17 Debye at the CCSD(T)/cc-pVTZ level. Therefore, 13 is also a suitable candidate for both laboratory and radioastronomical studies.

CSiGaAl2 -/0 and CGeGaAl2 -/0 having planar tetracoordinate carbon atoms in their global minimum energy structures

Density functional theory (DFT) is used to explore the structure, stability, and bonding in CSiGaAl₂ ^-/0 and CGeGaAl₂ ^-/0 systems having planar tetracoordinate carbon (ptC). The neutral systems have 17 valence electrons and the mono-anionic systems have 18 valence electrons. The ab initio molecular dynamics simulations for 2000 fs time at two different temperatures (300 and 500 K) supported the kinetic stability of the systems. From the natural bond orbital (NBO) analysis it is shown that there is a strong electron donation from the ligand atoms to the ptC atom. We have used Li⁺ ion for the neutralization of the mono-anionic systems and more interestingly it does not disrupt the planar structure. The most preferable site for binding of Li⁺ ion is along the AlAl bond in both of the mono-anionic systems. All the systems in this work have both σ and π aromaticity which is predicted from the computations of nucleus independent chemical shift (NICS). Although the anionic species obey the 18 valence electronic rule, the neutral systems break the rule with 17 valence electrons. However, both sets of systems are stable in the planar form. The bonding analysis of the systems includes molecular orbital, adaptive natural density partitioning (AdNDP), quantum theory of atoms in molecules (QTAIM), electron localization function (ELF) basin, and aromaticity analyses. The energy decomposition analysis (EDA) determines the interaction of Li⁺ ion with CSiGaAl₂ ^- and CGeGaAl₂ ^- in Li@SiGaAl₂ and Li@GeGaAl₂ , respectively.

Planar pentacoordinate carbon in a sulphur-surrounded boron wheel: the global minimum of CB5S5. Chem Commun (Camb) 2022;58:2552-2555. [PMID: 35103735 DOI: 10.1039/d1cc07313c]

We report a σ + π double aromatic CB₅S₅⁺ cluster, the first global minimum unusually having a planar hypercoordinate carbon inside a boron wheel. Five peripheral sulfur atoms stabilize the carbon-centered boron wheel by weakening the electron deficiency of the boron atoms through strong S → B π back-bonding.

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.

Sulphur-Bridged BAl5S5+ with 17 Counting Electrons: A Regular Planar Pentacoordinate Boron System

At present, most of the reported planar pentacoordinate clusters are similar to the isoelectronic substitution of CAl₅⁺, with 18 counting electrons. Meanwhile, the regular planar pentacoordinate boron systems are rarely reported. Hereby, a sulphur-bridged BAl₅S₅⁺ system with a five-pointed star configuration and 17 counting electrons is identified at the global energy minimum through the particle-swarm optimization method, based on the previous recognition on bridged sulphur as the peripheral tactics to the stable planar tetracoordinate carbon and boron. Its outstanding stability has been demonstrated by thermodynamic analysis at 900 K, electronic properties and chemical bonding analysis. This study provides adequately theoretical basis and referable data for its experimental capture and testing.

Photoelectron spectroscopy and theoretical study of AlnC5-/0 (n = 1-5) clusters: structural evolution, relative stability of star-like clusters, and planar tetracoordinate carbon structures

The AlnC5- (n = 1-5) clusters were detected in the gas-phase and were investigated via mass-selected anion photoelectron spectroscopy. The structures of AlnC5-/0 (n = 1-5) were explored by theoretical calculations. It is found that the structures of AlC5-/0 and Al2C5-/0 are linear while those of Al3C5-/0, Al4C5-/0, and Al5C5-/0 are two-dimensional. The most stable structures of AlC5-/0 and Al2C5-/0 are linear with the Al atoms attached to the ends of C5 chain. The most stable structures of Al3C5-/0 can be viewed as three Al atoms interacting with a nonlinear C5 chain. The most stable structure of Al4C5- anion is a planar structure composed of a C2 unit, a C3 unit, and two Al2 units, while that of the neutral Al4C5 cluster has four Al atoms connected to different positions of a distorted C5 chain. The global minimum structures of Al5C5-/0 are planar structures composed of an Al4C quadrilateral, two C2 groups, and an Al atom connected to two C2 groups. Planar tetracoordinate carbon (ptC) has been identified in the structures of both anionic and neutral Al5C5. It is worth mentioning that the star-like structure of Al5C5- is slightly higher in energy than the ground state structure. The comparison of theoretical calculations with the experimental spectra indicates the star-like structure of Al5C5- may also appear in our experiments.

Anion Photoelectron Spectroscopy and Theoretical Studies of Al4C6-/0: Global Minimum Triangle-Shaped Structures and Hexacoordinated Aluminum

Anion photoelectron spectroscopy and theoretical calculations were used to investigate the structural and bonding properties of Al₄C₆^-/0 clusters. The vertical detachment energy of Al₄C₆^- was measured to be 3.36 ± 0.08 eV. The structure of the Al₄C₆^- anion is confirmed to be a bowl-shaped distorted triangle with an Al atom at the center and three Al atoms at the vertices. The global minimum isomer of neutral Al₄C₆ has a planar triangle-shaped structure with D_3h symmetry. Both anionic and neutral Al₄C₆ have a hexacoordinated Al atom surrounded by three C≡C groups. Compared with the structure of neutral Al₄C₆, the structure of Al₄C₆^- is distorted due to the addition of the excess electron. The molecular orbital analysis shows that the singly occupied molecular orbital of Al₄C₆^- mainly locates on one side of the triangle plane and the neutral Al₄C₆ has a large highest occupied molecular orbital and lowest unoccupied molecular orbital gap. Theoretical calculations indicate that neutral Al₄C₆ has some aromaticity.

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.

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⁺.

The degree of π electron delocalization and the formation of 3D-extensible sandwich structures

DFT B3LYP/6-31G(d) calculations were performed to examine the feasibility of graphene-like C42H18 and starbenzene C6(BeH)6 (SBz) polymers as ligands of 3D-extensible sandwich compounds (3D-ESCs) with uninterrupted sandwich arrays. The results revealed that sandwich compounds with three or more C42H18 ligands were not feasible. The possible reason may be the localization of π electrons on certain C6 hexagons due to π-metal interactions, which makes the whole ligand lose its electronic structure basis (higher degree of π electron delocalization) to maintain the planar structure. For comparison, with the aid of benzene (Bz) molecules, the SBz polymers can be feasible ligands for designing 3D-ESCs because the C-Be interactions in individual SBz are largely ionic, which will deter the π electrons on one C6 ring from connecting to those on neighbouring C6 rings. This means that high degree of π electron delocalization is not necessary for maintaining the planarity of SBz polymers. Such a locally delocalized π electron structure is desirable for the ligands of 3D-ESCs. Remarkably, the formation of a sandwich compound with SBz is thermodynamically more favourable than that found for bis(Bz)chromium. The assembly of 3D-ESCs is largely exothermic, which will facilitate future experimental synthesis. The different variation trends on the HOMO-LUMO gaps in different directions (relative to the sandwich axes) suggest that they can be developed to form directional conductors or semiconductors, which may be useful in the production of electronic devices.

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.

How can carbon favor planar multi-coordination in boron-based clusters? Global structures of CBxEy2−(E = Al, Ga, x + y = 4)

CB₂E₂Mg (E = Al, Ga) designed in the present study represents the first successful design of a boron-based planar penta-coordinate carbon (ppC) structures as the global minima.

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.

AlxC Monolayer Sheets: Two-Dimensional Networks with Planar Tetracoordinate Carbon and Potential Applications as Donor Materials in Solar Cell

We perform a global search of the most stable structures of 2D stoichiometric AlxC (x = 1/3, 1, 2, and 3) monolayer sheets. In the most stable 2D planar AlC network, every carbon atom is tetracoordinated. In addition to the structure of AlC, structures of the most stable Al2C and Al3C monolayer sheets are also predicted for the first time. AlC and Al2C monolayers are semiconducting, while Al3C monolayer is metallic. In particular, Al2C monolayer possesses a bandgap of 1.05 eV (based on HSE06 calculation), a value suitable for photovoltaic applications. Moreover, three Al2C/WSe2, Al2C/MoTe2, and AlC/ZnO van der Waals heterobilayers are investigated, and their power conversion efficiencies are estimated to be in the range of 12-18%. The near-perfect match in lattice constants between the Al2C monolayer and PdO (100) surface suggests strong likelihood of experimental realization of the Al2C monolayer on the PdO (100) substrate.

Re-examination of the C(6)Li(6) structure: to be, or not to be symmetric

The potential energy surface of C6 Li6 was re-examined and a new non-symmetric global minimum was found. The new structure can be described as three C2 (2-) fragments strongly aggregated through lithium bridges. At high temperatures, fluxionality is perceived instead of dissociation. At 600 and 900 K, the BOMD simulations show that the lithium mobility is high, indicating that the cluster behaves in a liquid-like manner (BOMD=Born-Oppenheimer molecular dynamics).

Heavy periodane

The potential energy surface of the hypothetical NaMgAlSiPSCl system (heavy periodane) is exhaustively analyzed via the gradient embedded genetic algorithm (GEGA) in combination with density functional theory (DFT) computations. The electronegativity differences among the elements in both the second and third rows of the periodic table indicate that low-energy heavy periodane structures are obtained when highly electronegative and electropositive elements are bound together, but the global minimum of the heavy periodane system is completely different to its second-row analog (LiBeBCNOF).