Exploring the Potential Energy Surface of Trimetallic Deltahedral Zintl Ions: Lowest-Energy [Sn6Ge2Bi]3- and [(Sn6Ge2Bi)2]4- Structures

Exploration of the potential energy surface in mixed Zintl clusters applying an automatic Johnson polyhedra generator: the case of arachno E6M24- (E = Si, Ge, Sn; M = Sb, Bi)

A new algorithm called Automatic Johnson Cluster Generator (AJCG) is presented, which, as its name indicates, allows the definition of the desired Johnson polyhedron to subsequently carry out all the possible permutations between the atoms that form this polyhedron. This new algorithm allows the exhaustive study of the structures' potential energy surface (PES). In addition, the AJCG algorithm is helpful for the study of three-dimensional compounds such as boranes or Zintl clusters and their structural derivatives with two or more different atoms. The automatic filling of vertices is particularly useful in mixed compounds because of the possibility of taking into account all possible configurations in the structure. As a test system, we investigated the arachno-type E6M24- (E = Si, Ge, Sn; M = Sb, Bi) structure which has eight vertices and complies with Wade-Mingos rules. Initially, we defined a bipyramidal structure (10 vertices), and filled the vertices with the atoms in all possible configurations. Since the selected system has eight atoms, the two remaining vertices were filled with pseudo atoms to complete the structure. After re-optimizing the initial population generated with AJCG, a large number of isomers with energy below 10 kcal mol-1 are identified. These results show that the most stable isomers possess homonuclear M-M bonds, except Sn6Bi24-. Although the overall putative minima differ at the PBE0-D3 and DLPNO-CCSD(T) levels, they are always competitive minima. In addition to using high-precision methodologies to correctly study relative energies, applying solvent effects in highly charged systems becomes mandatory. The aromatic character of these studied systems was demonstrated qualitatively with two- and three-dimensional mapping and quantitatively by calculating the value of the z-component of the induced magnetic field at the cage center, including scalar and spin-orbit correction for relativistic effects. The compounds studied have a high degree of aromaticity, which allows us to establish that despite structural modifications (i.e., from closo to arachno), the aromaticity is preserved.

Functionalized nona-silicide [Si9R3] Zintl clusters: a new class of superhalogens

Superatoms, due to their various applications in redox and materials chemistry, have been a major topic of study in cluster science. Superhalogens constitute a special class of superatoms that mimic the chemistry of halogens and serve as building blocks of novel materials such as super and hyper salts, perovskite-based solar cells, solid-state electrolytes, and ferroelectric materials. These applications have led to a constant search for new class of superhalogens. In this study, using density functional theory, we show that recently synthesized [Si₉{Si (^tBu)₂H}₃] and [Si₉{Si (TMS)₃}₃] Zintl clusters not only behave like halogens but also when functionalized with suitable ligands exhibit superhalogen characteristics. Frontier molecular orbital (FMO) analyses give insights into the electron-accepting nature of the Zintl clusters. Additional bonding techniques such as energy density at the bond critical point (BCP) and adaptive natural density partitioning (AdNDP) gives complementary information about the nature of bonding in Si₉-based Zintl clusters. The potential of these Zintl clusters in the synthesis of new electrolytes in Li-ion batteries is also investigated.

Why an integrated approach between search algorithms and chemical intuition is necessary?

Though search algorithms are appropriate tools for identifying low-energy isomers, fixing several constraints seems to be a fundamental prerequisite to successfully running any structural search program. This causes some potential setbacks as far as identifying all possible isomers, close to the lowest-energy isomer, for any elemental composition. The number of explored candidates, the choice of method, basis set, and availability of CPU time needed to analyze the various initial test structures become necessary restrictions in resolving the issues of structural isomerism reasonably. While one could arrive at new structures through chemical intuition, reproducing or achieving those exact same structures requires increasing the number of variables in any given program, which causes further constraints in exploring the potential energy surface in a reasonable amount of time. Thus, it is emphasized here that an integrated approach between search algorithms and chemical intuition is necessary by taking the C₁₂O₂Mg₂ system as an example. Our initial search through the AUTOMATON program yielded 1450 different geometries. However, through chemical intuition, we found eighteen new geometries within 40.0 kcal mol^-1 at the PBE0-D3/def2-TZVP level. These results indirectly emphasize that an integrated approach between search algorithms and chemical intuition is necessary to further our knowledge in chemical space for any given elemental composition.

Ternary aromatic and anti-aromatic clusters derived from the hypho species [Sn2Sb5]3

Heterometallic clusters have attracted broad interests in the synthetic chemistry due to their various coordination modes and potential applications in heterogeneous catalysis. Here we report the synthesis, experimental, and theoretical characterizations of four ternary clusters ([M₂(CO)₆Sn₂Sb₅]^3- (M = Cr, Mo), and [(MSn₂Sb₅)₂]^4-, (M = Cu, Ag)) in the process of capturing the hypho- [Sn₂Sb₅]^3- in ethylenediamine (en) solution. We show that the coordination of the binary anion to transition-metal ions or fragments provides additional stabilization due to the formation of locally σ-aromatic units, producing a spherical aromatic shielding region in the cages. While in the case of [Mo₂(CO)₆Sn₂Sb₅]^3- stabilization arises from locally σ-aromatic three-centre and five-centre two-electron bonds, aromatic islands in [(AgSn₂Sb₅)₂]^4- and [(CuSn₂Sb₅)₂]^4- render them globally antiaromatic. This work describes the coordination chemistry of the versatile building block [Sn₂Sb₅]^3-, thus providing conceptual advances in the field of metal-metal bonding in clusters.

Structure, Stability, and Spectroscopic Properties of Small Acetonitrile Cation Clusters

Ionization and fragmentation pathways induced by ionizing agents are key to understanding the formation of complex molecules in astrophysical environments. Acetonitrile (CH₃CN), the simplest organic nitrile, is an important molecule present in the interstellar medium. In this work, DFT and MP2 calculations were performed in order to obtain the low energy structures of the most relevant cations formed from electron-stimulated ion desorption of CH₃CN ices. Selected reaction pathways and spectroscopic properties were also calculated. Our results indicate that the most stable acetonitrile cation structure is CH₂CNH⁺ and that hydrogenation can occur successively without isomerization steps until its complete saturation. Moreover, the stability of distinct cluster families formed from the interaction of acetonitrile with small fragments, such as CH_n⁺, C₂H_n⁺, and CH_nCNH⁺, is discussed in terms of their respective binding energies. Some of these molecular clusters are stabilized by hydrogen bonds, leading to species whose infrared features are characterized by a strong redshift of the N-H stretching mode. Finally, the rotational spectra of CH₃CN and protonated acetonitrile, CH₃CNH⁺, were simulated using distinct computational protocols based on DFT, MP2, and CCSD(T) considering centrifugal distortion, vibrational-rotational coupling, and vibrational anharmonicity corrections. By adopting an empirical scaling procedure for calculating spectroscopic parameters, we were able to estimate the rotational frequencies of CH₃CNH⁺ with an expected average error below 1 MHz for J values up to 10.

Searching for double σ- and π-aromaticity in borazine derivatives

Inspired by the double-aromatic (σ and π) C₆H₃⁺, C₆I₆²⁺, and C₆(SePh)₆²⁺ ring-shaped compounds, herein we theoretically study their borazine derivative analogues. The systems studied are the cation and dications with formulas B₃N₃H₃⁺, B₃N₃Br₆²⁺, B₃N₃I₆²⁺, B₃N₃(SeH)₆²⁺, and B₃N₃(TeH)₆²⁺. Our DFT calculations indicate that the ring-shaped planar structures of B₃N₃H₃⁺, B₃N₃I₆²⁺, and B₃N₃(TeH)₆²⁺ are more stable in the singlet state, while those of B₃N₃Br₆²⁺ and B₃N₃(SeH)₆²⁺ prefer the triplet state. Besides, exploration of the potential energy surface shows that the ring-shaped structure is the putative global minimum only for B₃N₃I₆²⁺. According to chemical bonding analysis, B₃N₃H₃⁺, B₃N₃I₆²⁺, and B₃N₃(TeH)₆²⁺ have σ and π delocalized bonds. The number of delocalized σ/π electrons is 2/6 for the first, and 10/6 for the second and third, similar to what their carbon analogs exhibit. Finally, the analysis of the magnetically induced current density allows B₃N₃H₃⁺, B₃N₃I₆²⁺, and B₃N₃(TeH)₆²⁺ to be classified as strongly σ aromatic, and poorly π aromatic compounds.

Evolutionary algorithms, Born–Oppenheimer molecular dynamics and the magnetic criteria of aromaticity have been used to evaluate the stability and σ–π aromaticity of borazine derivatives in order to expand the family of double aromatics systems.