Substantial π-aromaticity in the anionic heavy-metal cluster [Th@Bi12]4

Actinide endohedral inter-metalloid clusters of the group 15 elements

Inter-metalloid clusters in Zintl chemistry have been extensively studied due to their unique electronic structures and potential applications. In this work, we explored a series of actinide endohedral inter-metalloid clusters of the group 15 elements [An@Bi12]4- (An = Th-U) and [An@Sb12]4- using density functional theory (DFT). [Th@Bi12]4- and [U@Bi12]4- exhibit Cs symmetry, while [Pa@Bi12]4- and [An@Sb12]4- (An = Th-U) have C1 structures. Bonding analyses such as bond order, molecular orbitals (MO) and quantum theory of atoms in molecules (QTAIM) show covalent An-Bi/An-Sb bonding in the clusters. All these clusters are highly stable according to the studied formation reactions and may be accessible experimentally. Compared with [An@Bi12]4-, [An@Sb12]4- possesses stronger bonding interactions, mainly arising from the higher electrostatic interaction energy. For clusters with the same group 15 elements, the bonding interactions increase gradually from Th to U, which is mainly determined by the covalent interactions of An-Bi/An-Sb bonding. This work is expected to provide potential avenues for the construction of robust inter-metalloid clusters of the group 15 elements.

Supramolecular trapping of a cationic all-metal σ-aromatic {Bi4} ring

Aromaticity in organic molecules is well defined, but its role in metal-only rings remains controversial. Here we introduce a supramolecular stabilization approach of a cationic {Bi4} rhomboid within the symmetric charge sphere of two bowl-shaped dianionic calix[4]pyrrolato indinates. Crystallographic and spectroscopic characterization, quantum chemical analysis and magnetically induced ring currents indicate σ-aromaticity in the formally tetracationic 16-valence electron [Bi4]4+ ring. Computational screening for other p-block elements identifies the planar rhomboid as the globally preferred structure for 16-valence electron four-atomic clusters. The aromatic [Bi4]4+ is isoelectronic to the [Al4]4-, a motif previously observed as antiaromatic in Li3[Al4]- in the gas phase. Thus, subtle factors such as charge isotropy seem to decide over aromaticity or antiaromaticity, advising for caution in debates based on the Hückel model-a concept valid for second-row elements but less deterministic for the heavier congeners.

Core-electron contributions to the magnetic response of molecules with heavy elements and their significance in aromaticity assessments

This study delves into the magnetic response of core electrons and their influence on the global magnetic response of planar and three-dimensional systems containing heavy elements, employing the removing valence electron (RVE) approximation. We also explore electronic aromaticity indices to understand the potential role of core electrons on electron delocalization in the absence of an external perturbation. The study reveals that core electrons significantly contribute to the overall magnetic response, especially to the magnetic shielding, affecting the interpretation of aromaticity. In contrast, the calculation of the electronic aromaticity indices suggests a negligible participation of the core electrons on electron delocalization. Despite their widespread use, the study emphasizes caution in labeling systems as strongly aromatic based solely on shielding function computations. It is noteworthy to emphasize the limitations associated with each aromaticity criterion; particularly in the context of magnetic shielding function calculations, the core-electron effect contamination is undeniable. Hence, the integration of various criteria becomes imperative for attaining a comprehensive understanding of magnetic responses within complex systems.

Synthesis of bismuthanyl-substituted monomeric triel hydrides

The syntheses and characterizations of the first bismuthanylborane monomers stabilized only by a donor in D·BH2Bi(SiMe3)2 (D = DMAP 1a, IDipp 1b, IMe41c; DMAP = 4-dimethylaminopyridine, IDipp = 1,3-bis(2,6-diisopropylphenyl)-imidazolin-2-ylidene, IMe4 = 1,3,4,5-tetramethylimidazol-2-ylidene) are presented. All compounds were synthesized by salt metathesis reactions between D·BH2I and KBi(SiMe3)2(THF)0.3 and represent some of the extremely rare compounds featuring a 2c-2e B-Bi bond in a molecular compound. The products display high sensitivity towards air and light and slowly decompose in solution even at -80 °C. By the reaction of IDipp·GaH2(SO3CF3) with KBi(SiMe3)2(THF)0.3, the synthesis of the first bismuthanylgallane IDipp·GaH2Bi(SiMe3)2 (2) stabilized only by a 2-electron donor was possible, as evident from single crystal X-ray structure determination, NMR spectroscopy and mass spectrometry. Computational studies shed light on the stability of the products and the electronic nature of the compounds.

Reactivity of Tetrel-Functionalized Heptaphosphane Clusters toward Azides

In this work, the reactivity of tetrel-functionalized phosphorus clusters toward organoazides is probed. Clusters (Me3Si)3P7 (1) and (Me3Ge)3P7 (2) were reacted with benzyl azide, phenyl azide, and 4-bromophenyl azide, and it was found that the [RN] (R = benzyl, phenyl, and 4-bromophenyl) unit from the azide inserted into the phosphorus-tetrel bonds on the cluster, accompanied by N2 elimination. Through control of the azide stoichiometry, the mono-, bis-, and tris-inserted products could be observed, consistent with these insertions proceeding in a stepwise manner. The bonding between the amine moieties and clusters was further investigated by computational chemistry, and the findings were consistent with the phosphorus cluster having undergone a formal oxidation. These insertion reactions are a convenient means of accessing Zintl clusters functionalized with exo-nitrogen-bonded moieties, which, to the best of our knowledge, were previously unknown.

[Bi10{RuPPh3}3]-: Paramagnetic 13-Vertex Polybismuthide Heteroanion

Ru(PPh3)3Cl2 reacts with Binn- from an ethylenediamine (en) solution of K5Bi4 to yield a new architype of 13-vertex [Bi10{RuPPh3}3]- (1) composed of unprecedented incomplete cuboidal Bi73- and triangular Bi33- held together by {RuPPh3}2+.

On the aromaticity of actinide compounds

The chemistry of actinides has flourished since the late 2010s with the synthesis of new actinide complexes and clusters. On the theoretical side, a range of tools is available for the characterization of these heavy element-containing compounds, but discrepancies in the assessment of aromaticity using different tools have led to controversies. In this Perspective, we examine the origin of controversies relating to the aromaticity of metallic compounds, with a focus on actinides. The aromaticity of actinides is important, not because these molecules are numerous or have a special role in catalysis or reactivity, but because this topic pushes theories of aromaticity to their limits. Owing to its reference independence, the magnetic criterion of aromaticity has been the most popular choice for the characterization of the aromaticity of metallic compounds, including actinide compounds. Through examination of several case studies, we show why this criterion might be misleading for metallic species and explain how findings relating to actinide compounds could reshape theories of aromaticity, not just for actinides but perhaps also for well-known hydrocarbons.

Recent advances in atomic cluster synthesis: a perspective from chemical elements

Despite its potential significance, "cluster chemistry" remains a somewhat marginalized topic within the chemistry field. However, atomic clusters with their unusual and unique structures and properties represent a novel material group situated between molecules and nanoparticles or solid matter, judging from both scientific standpoints and historical backgrounds. Surveying an entire material group, including all substances that can be regarded as a cluster, is essential for establishing cluster chemistry as a more prominent chemistry field. This review aims to provide a comprehensive understanding by categorizing, summarizing, and reviewing clusters, focusing on their constituent elements in the periodic table. However, because numerous disparate synthetic processes have been individually developed to date, their straightforward and uniform classification is a challenging task. As such, comprehensively reviewing this field from a chemical composition viewpoint presents significant obstacles. It should be therefore noted that despite adopting a synthetic method-based classification in this review, the discussions presented herein could entail inaccuracies. Nevertheless, this unorthodox viewpoint unfolds a new scientific perspective which accentuates the common ground between different development processes by emphasizing the lack of a definitive border between their synthetic methods and material groups, thus opening new avenues for cementing cluster chemistry as an attractive chemistry field.

[GeRu6(CO)18HI]: A Germanium-Centered Ruthenium Carbonyl Cluster with Aromatic Ring Current

The carbonyl cluster compound [GeRu6(CO)18HI] is unique in regard to its structure and bonding with a GeRu6 cluster core, a planar GeRu4HI unit, extensive multi-center bonding, and an aromatic ring current similar to benzene (9-10 nA T-1). The open-shell cluster core is a Ge-centered five-membered Ru4(Ru2) ring with CO ligands and an additional H and I atom, each bridging two Ru atoms on opposite sides of the cluster core. The compound is prepared at 130 °C in a weakly-coordinating ionic liquid.

Multiconfigurational actinide nitrides assisted by double Möbius aromaticity

Understanding the bonding nature between actinides and main-group elements remains a key challenge in actinide chemistry due to the involvement of f orbitals. Herein, we propose a unique "aromaticity-assisted multiconfiguration" (AAM) model to elucidate the bonding nature in actinide nitrides (An2N2, An = Ac, Th, Pa, U). Each planar four-membered An2N2 with equivalent An-N bonds possesses four delocalized π electrons and four delocalized σ electrons, forming a new family of double Möbius aromaticity that contributes to the molecular stability. The unprecedented aromaticity further supports actinide nitrides to exhibit multiconfigurational characters, where the unpaired electrons (2, 4 or 6 in naked Th2N2, Pa2N2 or U2N2, respectively) either are spin-free and localized on metal centres or form metal-ligand bonds. High-level multiconfigurational computations confirm an open-shell singlet ground state for actinide nitrides, with small energy gaps to high spin states. This is consistent with the antiferromagnetic nature observed experimentally in uranium nitrides. The novel AAM bonding model can be authenticated in both experimentally identified compounds containing a U2N2 motif and other theoretically modelled An2N2 clusters and is thus expected to be a general chemical bonding pattern between actinides and main-group elements.

Effect of La3+ on the Formation of Endohedral Zintl Clusters Featuring In/Bi Shells

Investigating the interactions of f-block metal ions with p-block polyanions in multinary cluster compounds is becoming increasingly attractive but remains a challenge in terms of both the synthetic approach and the control of the structures that are formed during the syntheses. So far, two types of reactions were dominant for the formation of corresponding clusters: the reaction of binary anions of p-block elements in 1,2-diamino-ethane (en) solutions or the reaction of organobismuth compounds with corresponding f-block metal complexes in THF. Herein, we report the synthesis of [La@In2Bi11]4- (1) and its doubly μ-Bi-bridged analogue in the doubly [K(crypt-222)]+-coordinated {[K(crypt-222)]2[La@In2Bi11](μ-Bi)2[La@In2Bi11]}4- (2) as their [K(crypt-222)]+ salts [K(crypt-222)]41 and [K(crypt-222)]42, respectively, achieved by reactions of [InMes3] and [La(C5Me4H)3] (Mes = mesityl, C5Me4H = tetramethylcyclopentadienyl) with K10Ga3Bi6.65/crypt-222 (crypt-222 = 4,7,13,16,21,24-hexaoxa-1,10-diazabicyclo[8.8.8]hexacosane) in en. In the absence of [La(C5Me4H)3], the otherwise unchanged reactions afford the anions [Bi6(InMes)(InMes2)]3- (3) and [Mes3In-InMes3]2- (4) instead, which can be isolated as their [K(crypt-222)]+ salts [K(crypt-222)]33 and [K(crypt-222)]24·tol (tol = toluene), respectively. The {Bi6} fragment observed in anion 3 is assumed to be one of the key intermediates not only toward the formation of 1 and 2 but also on the way to more general bismuth rich compounds.

[K2(Bi@Pd12@Bi20)]4-: An Endohedral Inorganic Fullerene with Spherical Aromaticity

Inorganic fullerene clusters have attracted widespread attention due to their highly symmetrical geometric structures and intrinsic electronic properties. However, cage-like clusters composed of heavy metal elements with high symmetry are rarely reported, and their synthesis is also highly challenging. In this study, we present the synthesis of a [K2(Bi@Pd12@Bi20)]4- cluster that incorporates a {Bi20} cage with pseudo-Ih symmetry, making it the largest main group metal cluster compound composed of the bismuth element to date. Magnetic characterization and theoretical calculations suggest that the spin state of the overall cluster is a quartet. Quantum chemical calculations reveal that the [Bi20]3- cluster has a similar electronic configuration to C606- and the [Bi@Pd12@Bi20]6- cluster exhibits a unique open-shell aromatic character.

f-Element Zintl Chemistry: Actinide-Mediated Dehydrocoupling of H2Sb1- Affords the Trithorium and Triuranium Undeca-Antimontriide Zintl Clusters [{An(TrenTIPS)}3(μ3-Sb11)] (An = Th, U; TrenTIPS = {N(CH2CH2NSiiPr3)3}3-)

Reaction of the cesium antimonide complex [Cs(18C6)2][SbH2] (1, 18C6 = 18-crown-6 ether) with the triamidoamine actinide separated ion pairs [An(TrenTIPS)(L)][BPh4] (TrenTIPS = {N(CH2CH2NSiiPr3)3}3-; An/L = Th/DME (2Th); U/THF (2U)) affords the triactinide undeca-antimontriide Zintl clusters [{An(TrenTIPS)}3(μ3-Sb11)] (An = Th (3Th), U (3U)) by dehydrocoupling. Clusters 3Th and 3U provide two new examples of the Sb113- Zintl trianion and are unprecedented examples of molecular Sb113- being coordinated to anything since all previous reports featured isolated Sb113- Zintl trianions in separated ion quadruple formulations with noncoordinating cations. Quantum chemical calculations describe dominant ionic An-Sb interactions in 3Th and 3U, though the data suggest that the latter exhibits slightly more covalent An-Sb linkages than the former. Complexes 3Th and 3U have been characterized by single crystal X-ray diffraction, NMR, IR, and UV/vis/NIR spectroscopies, elemental analysis, and quantum chemical calculations.

Thorium(IV)-antimony complexes exhibiting single, double, and triple polar covalent metal-metal bonds

There is continued burgeoning interest in metal-metal multiple bonding to further our understanding of chemical bonding across the periodic table. However, although polar covalent metal-metal multiple bonding is well known for the d and p blocks, it is relatively underdeveloped for actinides. Homometallic examples are found in spectroscopic or fullerene-confined species, and heterometallic variants exhibiting a polar covalent σ bond supplemented by up to two dative π bonds are more prevalent. Hence, securing polar covalent actinide double and triple metal-metal bonds under normal experimental conditions has been a fundamental target. Here we exploit the protonolysis and dehydrocoupling chemistry of the parent dihydrogen-antimonide anion, to report one-, two- and three-fold thorium-antimony bonds, thus introducing polar covalent actinide-metal multiple bonding under normal experimental conditions between some of the heaviest ions in the periodic table with little or no bulky-substituent protection at the antimony centre. This provides fundamental insights into heavy element multiple bonding, in particular the tension between orbital-energy-driven and overlap-driven covalency for the actinides in a relativistic regime.

Consistent Analytical Second Derivatives of the Kohn-Sham DFT Energy in the Framework of the Conductor-Like Screening Model through Gaussian Charge Distributions

The use of implicit solvation models such as the conductor-like screening model (COSMO) in quantum chemical calculations is very common, as both a rough estimate of solvation effects as well as a general tool for stabilizing ionic molecular structures. In order to generate a smooth potential energy surface as well as consistent gradients, it is necessary to apply the Gaussian charge model (GCM) for the COSMO charges. This work introduces an efficient implementation for consistent analytical second derivatives of the electronic energy with COSMO-GCM in the framework of the Kohn-Sham density functional theory. This is used to investigate the infrared spectroscopy of amino acids in aqueous solution, where the impact of pH on the molecular structure and vibrational spectra is examined. Furthermore, the structure and stability of selected all-metal aromatic cluster ions are assessed.

[Ag(Sn9-Sn9)]5- and [(η4-Sn9)Ag(η1-Sn9)]7-, as aggregates of spherical aromatic building blocks. Persistence of aromaticity upon cluster gathering

Formation of cluster-based materials requires a fundamental understanding of the resulting cluster aggregation processes. The Sn94- Zintl-ion structure can be viewed as a building block featuring a spherical aromatic species, leading to a cluster gathering upon oxidative coupling and/or mediated by transition metals. Here, we evaluate the spherical aromatic properties of [Sn9-Sn9]6-, [Ag(Sn9-Sn9)]5- and [(η4-Sn9)Ag(η1-Sn9)]7-, as aggregates of two Sn9 building units held together via oxidative coupling and mediated by a Ag(I) transition metal center. Our results from magnetic criteria of aromaticity show that the inherent spherical aromatic characteristics of the parent Sn94- cluster are persistent in the overall aggregate where the enabled shielding cones ascribed to each Sn9 unit are able to interplay between them, leading to an overlap of the shielding regions. Hence, the two approaches for bringing cluster units together are able to retain the inherent spherical aromatic features for each Sn9 unit, leading to a cluster-based dimer where the parent properties remain. Thus, further cluster-based materials can be envisaged from aggregation upon oxidative coupling and/or mediated by transition metals, where the constituent building blocks retain their initial features, useful to guide the formation of more complex cluster-based aggregates.

f-Element heavy pnictogen chemistry

The coordination and organometallic chemistry of the f-elements, that is group 3, lanthanide, and actinide ions, supported by nitrogen ligands, e.g. amides, imides, and nitrides, has become well developed over many decades. In contrast, the corresponding f-element chemisty with the heavier pnictogen analogues phosphorus, arsenic, antimony, and bismuth has remained significantly underdeveloped, due largely to a lack of suitable synthetic methodologies and also the inherent hard(f-element)-soft(heavier pnictogen) acid-base mismatch, but has begun to flourish in recent years. Here, we review complexes containing chemical bonds between the f-elements and heavy pnictogens from phosphorus to bismuth that spans five decades of endeavour. We focus on complexes whose identity has been unambiguously established by structural authentication by single-crystal X-ray diffraction with respect to their synthesis, characterisation, bonding, and reactivity, in order to provide a representative overview of this burgeoning area. By highlighting that much has been achieved but that there is still much to do this review aims to inspire, focus and guide future efforts in this area.

Exact two-component theory becoming an efficient tool for NMR shieldings and shifts with spin-orbit coupling

We present a gauge-origin invariant exact two-component (X2C) approach within a modern density functional framework, supporting meta-generalized gradient approximations such as TPSS and range-separated hybrid functionals such as CAM-B3LYP. The complete exchange-correlation kernel is applied, including the direct contribution of the field-dependent basis functions and the reorthonormalization contribution from the perturbed overlap matrix. Additionally, the finite nucleus model is available for the electron-nucleus potential and the vector potential throughout. Efficiency is ensured by the diagonal local approximation to the unitary decoupling transformation in X2C as well as the (multipole-accelerated) resolution of the identity approximation for the Coulomb term (MARI-J, RI-J) and the seminumerical exchange approximation. Errors introduced by these approximations are assessed and found to be clearly negligible. The applicability of our implementation to large-scale calculations is demonstrated for a tin pincer-type system as well as low-valent tin and lead complexes. Here, the calculation of the Sn nuclear magnetic resonance shifts for the pincer-type ligand with about 2400 basis functions requires less than 1 h for hybrid density functionals. Further, the impact of spin-orbit coupling on the nucleus-independent chemical shifts and the corresponding ring currents of all-metal aromatic systems is studied.

On-Surface Synthesis and Real-Space Visualization of Aromatic P3 N3

On-surface synthesis is at the verge of emerging as the method of choice for the generation and visualization of unstable or unconventional molecules, which could not be obtained via traditional synthetic methods. A case in point is the on-surface synthesis of the structurally elusive cyclotriphosphazene (P3 N3 ), an inorganic aromatic analogue of benzene. Here, we report the preparation of this fleetingly existing species on Cu(111) and Au(111) surfaces at 5.2 K through molecular manipulation with unprecedented precision, i.e., voltage pulse-induced sextuple dechlorination of an ultra-small (about 6 Å) hexachlorophosphazene P3 N3 Cl6 precursor by the tip of a scanning probe microscope. Real-space atomic-level imaging of cyclotriphosphazene reveals its planar D3h -symmetric ring structure. Furthermore, this demasking strategy has been expanded to generate cyclotriphosphazene from a hexaazide precursor P3 N21 via a different stimulation method (photolysis) for complementary measurements by matrix isolation infrared and ultraviolet spectroscopy.

Stabilization of Pu(IV) in PuBr4(OPCy3)2 and Comparisons with Structurally Similar ThX4(OPR3)2 (R = Cy, Ph) Molecules

The pursuit of a trivalent plutonium halide phosphine oxide compound, e.g., "PuBr3(OPR)3," instead led to the isolation of the tetravalent trans-PuIVBr4(OPCy3)2, PuBr/Cy, compound by spontaneous oxidation of PuIII. The donating nature of phosphine oxides has allowed the isolation and characterization of PuBr/Cy by crystallographic, multinuclear NMR, solid state, and solution phase UV-vis-NIR spectroscopic techniques. The presence of a putative plutonyl(VI) complex formulated as "trans-PuVIO2Br2(OPCy3)2" was also observed spectroscopically and tentatively by single-crystal X-ray diffraction as a cocrystal of PuBr/Cy. A series of trans-ThX4(OPCy3)2 (X = Cl, ThCl/Cy; Br, ThBr/Cy; I, ThI/Cy) complexes were synthesized for comparison to PuBr/Cy. The triphenylphosphine oxide, OPPh3, complexes, trans-AnI4(OPPh3)2 (An = Th, ThI/Ph; U, UI/Ph), were also synthesized for comparison, completing the series trans-UX4(OPPh3)2 (X = Cl, Br, I), UX/Ph. To enable the synthesis of ThI/Cy and ThI/Ph, a new nonaqueous thorium iodide starting material, ThI4(Et2O)2, was synthesized. The syntheses of organic solvent soluble ThI4L2 (L = Et2O, OPCy3, and OPPh3) are the first examples of crystallographically characterized neutral thorium tetraiodide materials beyond binary ThI4. To show the viability of ThI4(Et2O)2 as a starting material for organothorium chemistry, (C5Me4H)3ThI was synthesized and crystallographically characterized.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

The s-p Nonhybrid Nature Causes Adaptive Superatomic States of Bismuth Clusters

We report a joint experimental and theoretical study on the stability and reactivity of Bin + (n=5-33) clusters. The alternating odd-even effect on the reaction rates of Bin + clusters with NO is observed, and Bi7 + finds the most inertness. First-principles calculation results reveal that the lowest energy structures of Bi6-9 + exhibit quasi-spherical geometry pertaining to the jellium shell model; however, the Bin + (n≥10) clusters adopt assembly structures. The prominent stability of Bi7 + is associated with its highly symmetric structure and superatomic states with a magic number of 34e closed shell. For the first time, we demonstrate that the unique s-p nonhybrid feature in bismuth rationalizes the stability of Bi6-9 + clusters within the jellium model, by filling the 6s electrons into the superatomic orbitals (forming "s-band"). Interestingly, the stability of 18e "s-band" coincides with the compact structure for Bin + at n≤9 but assembly structures for n≥10, showing an accommodation of the s electrons to the geometric structure. The atomic p-orbitals also allow to form superatomic orbitals at higher energy levels, contributing to the preferable structures of tridentate binding units. We illustrate the s-p nonhybrid nature accommodates the structure and superatomic states of bismuth clusters.

Insight into the formation of bismuth-tungsten carbonyl clusters

Multimetallic clusters play a key role as models to doped metals, as candidates to new types of superatomic catalysts and as precursors to new multimetallic solids. Understanding formation pathways is an essential and necessary step forward in the development of cluster synthesis and research, yet remains considerably lacking owing to difficulty in identification of intermediates and the ill-defined nature of common starting materials. Here we show progress in this regard by investigating the reactivity of an intermetallic solid of nominal composition 'K₅Ga₂Bi₄' with [W(cod)(CO)₄] upon extraction with ethane-1,2-diamine (en) and 4,7,13,16,21,24-hexaoxa-1,10-diazabicyclo[8.8.8]hexacosane (crypt-222). Several polybismuthide intermediates and by-products were identified along the reaction pathway, ultimately forming the new polybismuthide salt [K(crypt-222)]₃[µ:η³-Bi₃{W(CO)₃}₂]∙en∙tol. DFT calculations revealed plausible reaction schemes for the transformations taking place in the reaction mixture providing insight into the complex reactivity of 'K₅Ga₂Bi₄' on the basis of in situ generation of Bi₂^2-.

Planar Four-Membered Diboron Actinide Compound with Double Möbius Aromaticity

The Möbius rule predicts that a planar four-membered metallacycle can be aromatic with four mobile electrons, but such a simple ring has escaped recognition because it usually favors Hückel anti-aromaticity. Here, we report that a quasi-square four-membered actinide compound (Pa₂B₂) is doubly Möbius aromatic. Chemical bonding analyses reveal that this diboron protactinium molecule has four delocalized π electrons in addition to four delocalized σ electrons, satisfying the 4n Möbius rule for both σ and π components. Energetically, the block-localized wavefunction method, which is the simplest variant of ab initio valence bond theory, shows that the delocalization energy for the π and σ electrons reaches up to 65.0 and 72.3 kcal/mol, respectively, while the extra cyclic resonance energy (ECRE) amounts to 45 kcal/mol. The large positive ECRE values strongly confirm the unprecedented double Möbius aromaticity in Pa₂B₂. We anticipate that this new type of aromatic molecule can enrich the concept of Möbius aromaticity and open a new avenue for actinide compounds.

Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy

The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.

φ-Aromaticity in prismatic {Bi6}-based clusters

The occurrence of aromaticity in organic molecules is widely accepted, but its occurrence in purely metallic systems is less widespread. Molecules comprising only metal atoms (M) are known to be able to exhibit aromatic behaviour, sustaining ring currents inside an external magnetic field along M-M connection axes (σ-aromaticity) or above and below the plane (π-aromaticity) for cyclic or cage-type compounds. However, all-metal compounds provide an extension of the electrons' mobility also in other directions. Here, we show that regular {Bi₆} prisms exhibit a non-localizable molecular orbital of f-type symmetry and generate a strong ring current that leads to a behaviour referred to as φ-aromaticity. The experimentally observed heterometallic cluster [{CpRu}₃Bi₆]^-, based on a regular prismatic {Bi₆} unit, displays aromatic behaviour; according to quantum chemical calculations, the corresponding hypothetical Bi₆^2- prism shows a similar behaviour. By contrast, [{(cod)Ir}₃Bi₆] features a distorted Bi₆ moiety that inhibits φ-aromaticity.

[M@Sn14-xSbx]q- (M = La, Ce, or U; x = 6-8; q = 3, 4): Interaction of 4f or 5f Metal Ions with 5p Metal Atoms in Intermetalloid Clusters

The impact of 4f metal ions Ln³⁺ (Ln = La or Ce) versus 5f metal ions Uⁿ⁺ (n = 3 or 4) on the compositions and distribution of 5p metal atoms in the cluster shells of endohedral species [M@Sn_14-xSb_x]^q- (M = La, Ce, or U; x = 6-8; q = 3, 4) was studied by means of combined experimental and quantum chemical investigations. While all known f-block metal ion-centered endohedral clusters possessed combinations of larger main group metal atoms so far (Sn/Bi or Pb/Bi), resulting in mixtures of 13- and 14-atom cages, the 14-atom cages reported herein comprise exclusively Sn and Sb atoms and therefore are challenged in accommodating the large 4f and 5f ions. We show that the clusters form in reactions of (Sn₂Sb₂)^2- anions with [Ln(C₅Me₄H)₃] or [U(C₅Me₄H)₃Cl], and that salts of [La@Sn₆Sb₈]^3-, [La@Sn₇Sb₇]^4-, [U@Sn₈Sb₆]^4-, and [U@Sn₇Sb₇]^3- can be isolated from them. The assignment of Sn versus Sb in the encapsulating cage follows a simple rule. Different central atoms cause only slight differences in this regard and with respect to distortions of the cluster shells. The reactions also yielded the salt of the new binary anion (Sn₄Sb₄)^2- that was recently predicted by quantum chemical studies.

Actinide Pnictinidene Chemistry: A Terminal Thorium Parent-Arsinidene Complex Stabilised by a Super-Bulky Triamidoamine Ligand

We report the direct synthesis of the terminal pnictidenes [An(TrenTCHS )(PnH)][M(2,2,2-cryptand)] (TrenTCHS ={N(CH2 CH2 NSiCy3 )3 }3- ; An/Pn/M=Th/P/Na 5, Th/As/K 6, U/P/Na 7, U/As/K 8) and their conversion to the pnictides [An(TrenTCHS )(PnH2 )] (An/Pn=Th/P 9, Th/As 10, U/P 11, U/As 12). Use of the super-bulky TrenTCHS ligand was essential to accessing complete families, and 6 is an unprecedented example of a terminal thorium-arsinidene complex and only the second structurally authenticated parent terminal arsinidene complex of any metal. Comparison of the terminal Th=AsH unit of 6 to the bridging ThAs(H)K linkage in structurally analogous [Th(TrenTIPS ){μ-As(H)K(15-crown-5)}] (TrenTIPS ={N(CH2 CH2 NSiPri 3 )3 }3- ) reveals a stronger Th-As bond in the former compared to the latter, and a large response overall to the nature of the Th=AsH bonding upon removal of the electrostatically-bound K-ion; the σ-bond changes little but the π-bond is significantly perturbed.

Reactive Solubilization of Heterometallic Clusters by Treatment of (TrBi3 )2- Anions (Tr=Ga, In, Tl) with [Mn{N(SiMe3 )2 }2 ]

Lowering the charge of Zintl anions by (element-)organic substituents allows their use as sources of (semi)metal nanostructures in common organic solvents, as realized for group 15 anions or Ge9 4- and Sn9 4- . We developed a new strategy for other anions, using low-coordinate 3d metal complexes as electrophiles. [K(crypt-222)]+ salts of (TrBi3 )2- anions dissolved in situ in Et2 O and/or THF when reacted with [Mn(hmds)2 ]. Work-up afforded soluble [K(crypt-222)]+ salts of [{(hmds)2 Mn}2 (TlBi3 )]2- (in 1), [{(hmds)2 Mn}2 (Bi2 )]2- (in 2), and [{(hmds)Mn}4 (Bi2 )2 ]2- (in 3) (crypt-222=4,7,13,16,21,24-hexaoxa-1,10-diazabicyclo[8.8.8]hexacosane; Tr=Ga, In, Tl; hmds=N(SiMe3 )2 ), representing rare cases of Zintl clusters with open-shell metal atoms. 1 comprises the first coordination compound of the (TlBi3 )2- anion, 2 features a diamond-shaped {Pn2 M2 } unit, and 3 is a mixed-valent MnI /MnII compound. The uncommon electronic structures in 1-3 and magnetic coupling were studied by comprehensive DFT calculations.

Between Elemental Match and Mismatch: From K12 Ge3.5 Sb6 to Salts of (Ge2 Sb2 )2- , (Ge4 Sb12 )4- , and (Ge4 Sb14 )4

The solid mixture "K2 GeSb" was shown to comprise single-crystalline K12 Ge3.5 Sb6 (1), a double salt of K5 [GeSb3 ] with carbonate-like [GeSb3 ]5- anions, and the metallic Zintl phase K2 Ge1.5 . Extraction of 1 with ethane-1,2-diamine in the presence of crypt-222 afforded [K(crypt-222)]+ salts of several novel binary Zintl anions: (Ge2 Sb2 )2- (in 2), (Ge4 Sb12 )4- (in 3), and in the presence of [AuMePPh3 ] also (Ge4 Sb14 )4- (in 4). The anion in 2 represents a predicted, yet heretofore missing pseudo-tetrahedral anion. 4 comprises a cluster analogous to (Ge4 Bi14 )4- and (Ga2 Bi16 )4- , and thus one of the most Sb-rich binary p-block anions. The unprecedented cluster topology in 3 can be viewed as a defect-version of the one in 4 upon following a "dead end" of cluster growth. The findings indicate that Ge and Sb atoms are at the border of a well-matching and a mismatch elemental combination. We discuss the syntheses, the geometric structures, and the electronic structures of the new compounds.

On the Bonding Nature in the Crystalline Tri-Thorium Cluster: Core-Shell Syngenetic σ-Aromaticity

A unique thorium-thorium bond was observed in the crystalline tri-thorium cluster [{Th(η8 -C8 H8 )(μ3 -Cl)2 }3 {K(THF)2 }2 ]∞ , though the claim of σ-aromaticity for Th3 bond has been questioned. Herein, a new type of core-shell syngenetic bonding model is proposed to describe the stability of this tri-thorium cluster. The model involves a 3c-2e bond in the Th3 core and a multicentered (ThCl2 )3 charge-shift bond with 12 electrons scattering along the outer shell. To differentiate the strengths of the 3c-2e bond and the charge-shift bond, the block-localized wavefunction (BLW) method which falls into the ab initio valence bond (VB) theory is employed to construct a strictly core/shell localized state and its contributing covalent resonance structure for the Th3 core bond. By comparing with the σ-aromatic H3 + and nonaromatic Li3 + , the computed resonance energies and extra cyclic resonance energies confirm that this Th3 core bond is truly delocalized and σ-aromatic.

A Planar Five-Membered Aromatic Ring Stabilized by Only Two π-Electrons

Many chemicals known today are partially or fully aromatic, since a ring framework experiences additional stabilization through the delocalization of π-electrons. While aromatic rings with equal numbers of π-electrons and ring atoms such as benzene are particularly stable, those with the minimally required two π-electrons are very rare and yet remain limited to three- and four-membered rings if not stabilized in the coordination sphere of heavy metals. Here we report the facile synthesis of a dipotassium cyclopentagallene, a unique example of a five-membered aromatic ring stabilized by only two π-electrons. Single-crystal X-ray diffraction revealed a planar Ga5 ring with almost equal gallium-gallium bond lengths, which together with computational and spectroscopic data confirm its aromatic character. Our results prove that aromatic stabilization goes far beyond what has previously been assumed as minimum π-electron count in a five-atom ring fragment.

Paramagnetic NMR Shielding Tensors Based on Scalar Exact Two-Component and Spin-Orbit Perturbation Theory

The temperature-dependent Fermi-contact and pseudocontact terms are important contributions to the paramagnetic NMR shielding tensor. Herein, we augment the scalar-relativistic (local) exact two-component (X2C) framework with spin-orbit perturbation theory including the screened nuclear spin-orbit correction for the EPR hyperfine coupling and g tensor to compute these temperature-dependent terms. The accuracy of this perturbative ansatz is assessed with the self-consistent spin-orbit two-component and four-component treatments serving as reference. This shows that the Fermi-contact and pseudocontact interaction is sufficiently described for paramagnetic NMR shifts; however, larger deviations are found for the EPR spectra and the principle components of the EPR properties of heavy elements. The impact of the perturbative treatment is further compared to that of the density functional approximation and the basis set. Large-scale calculations are routinely possible with the multipole-accelerated resolution of the identity approximation and the seminumerical exchange approximation, as shown for [CeTi₆O₃(OiPr)₉(salicylate)₆].

Communication: Impact of the current density on paramagnetic NMR properties

Meta-generalized gradient approximations (meta-GGAs) and local hybrid functionals generally depend on the kinetic energy density τ. For magnetic properties, this necessitates generalizations to ensure gauge invariance. In most implementations, τ is generalized by incorporating the external magnetic field. However, this introduces artifacts in the response of the density matrix and does not satisfy the iso-orbital constraint. Here, we extend previous approaches based on the current density to paramagnetic NMR shieldings and EPR g-tensors. The impact is assessed for main-group compounds and transition-metal complexes considering 25 density functional approximations. It is shown that the current density leads to substantial improvements-especially for the popular Minnesota and SCAN functional families. Thus, we strongly recommend to use the current density generalized τ in paramagnetic NMR and EPR calculations with meta-GGAs.

Inorganic Ferrocene Analogue [Fe(P4)2]2

Inorganic metallocene derivatives containing only cyclo-P_n ligands have been targeted for more than 20 years, but their syntheses have never been achieved by pursuing the conventional route of using P₄ phosphorus except for the generation of [Ti(η⁵-P₅)₂]^2-. Herein, we report a facile one-step method for the synthesis of the homoleptic iron complex [Fe(P₄)₂]^2- by the Zintl-phase-type precursor KP. ³¹P NMR analyses indicate that upon dissolving the KP phase in ethylenediamine P₄^2- was generated only in the presence of 2,2,2-crypt. The amounts of cation-sequestering agents, the type of iron precursor, and their consuming ratio have a decisive impact on the yield of [Fe(P₄)₂]^2-. Both the Fe^II and the Fe^III precursors can oxidize P₄^2- to give a concomitant product [(P₇)Fe(P₄)]^3-, which can be partially inhibited by the addition of potassium to produce relatively pure crystalline [K(2,2,2-crypt)]₂[Fe(P₄)₂].

Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance

We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C₅Me₅)₂Y(μ-S)₂Mo(μ-S)₂Y(C₅Me₅)₂]^-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C₆Cl₅)₄]^- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.

Introduction of plumbole to f-element chemistry

Herein, we present the synthesis and characterization of heteroleptic lanthanide complexes bearing a dianionic η5-plumbole ligand in their coordination sphere. The reaction proceeds via a salt elimination reaction between the dilithioplumbole ([Li(thf)]2[1,4-bis-tert-butyl-dimethylsilyl-2,3-bis-phenyl-plumbolyl] = [Li2(thf)2(η5-LPb)]) and specifically designed [Ln(η8-COTTIPS)BH4] precursors (Ln = lanthanide, La, Ce, Sm, Er; COTTIPS = 1,4-bis-triisopropylsilyl-cyclooctatetraenyl), that are capable of stabilizing a planar plumbole moiety in the coordination sphere of different trivalent lanthanide ions. In-depth ab initio calculations show that the aromaticity of the dianionic plumbole is retained upon coordination. Electron delocalization occurs from the plumbole HOMO to an orbital of mainly d-character at the lanthanide ion. The magnetic properties of the erbium congener were investigated in detail, leading to the observation of magnetic hysteresis up to 5 K (200 Oe s-1), an unequivocal proof for single molecule magnet behavior in this system. The magnetic behavior of the erbium species can be modulated by manipulating the position of the lithium cation in the complex, which directly influences the bonding metrics in the central [(η5-LPb)Er(η8-COTTIPS)]- fragment. This allowed us to assess a fundamental magneto-structural correlation in an otherwise identical inner coordination sphere.

NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation

We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.

Electronic structure and bonding in endohedral Zintl clusters

Endohedral Zintl clusters-multi-metallic anionic molecules in which a d-block or f-block metal atom is enclosed by p-block (semi)metal atoms-are very topical in contemporary inorganic chemistry. Not only do they provide insight into the embryonic states of intermetallic compounds and show promise in catalytic applications, they also shed light on the nature of chemical bonding between metal atoms. Over the past two decades, a plethora of endohedral Zintl clusters have been synthesized, revealing a fascinating diversity of molecular architectures. Many different perspectives on the bonding in them have emerged in the literature, sometimes complementary and sometimes conflicting, and there has been no concerted effort to classify the entire family based on a small number of unifying principles. A closer look, however, reveals distinct patterns in structure and bonding that reflect the extent to which valence electrons are shared between the endohedral atom and the cluster shell. We show that there is a much more uniform relationship between the total valence electron count and the structure and bonding patterns of these clusters than previously anticipated. All of the p-block (semi)metal shells can be placed on a ladder of total valence electron count that ranges between 4n+2 (closo deltahedra), 5n (closed, three-bonded polyhedra) and 6n (crown-like structures). Although some structural isomerism can occur for a given electron count, the presence of a central metal cation imposes a preference for rather regular and approximately spherical structures which maximise electrostatic interactions between the metal and the shell. In cases where the endohedral metal has relatively accessible valence electrons (from the d or f shells), it can also contribute its valence electrons to the total electron count of the cluster shell, raising the effective electron count and often altering the structural preferences. The electronic situation in any given cluster is considered from different perspectives, some more physical and some more chemical, in a way that highlights the important point that, in the end, they explain the same situation. This article provides a unifying perspective of bonding that captures the structural diversity across this diverse family of multimetallic clusters.

Tetrahedral [Sb(AuMe)4 ]3- Occurring in Multimetallic Cluster Syntheses: About the Structure-Directing Role of Methyl Groups

The anion of [K(crypt-222)]3 [Sb(AuMe)4 ]⋅py (1; crypt-222=4,7,13,16,21,24-hexaoxa-1,10-diazabicyclo[8.8.8]hexacosane; py=pyridine) represents a rare example of a homoleptic heavy p-block metal atom being surrounded by four free-standing transition metal complex fragments, and the third example for a corresponding Sb compound. In contrast to all reported complexes of this type, the transition metal atoms possess twofold coordination only, hence the complex as a whole does not exhibit significant steric shielding or further linkage of the metal atoms. This is reflected in a high flexibility, as confirmed by slight deviations from a tetrahedral coordination of the Sb atom in the crystal and soft vibrational modes. An alternative pyramidal conformer, observed for a related arsenic compound with terminal phosphine ligands, is apparently disfavored owing to electron correlation effects. The compound is formed in a reaction that in another solvent or at other reactant concentrations yields salts of ternary cluster anions. By a combined experimental and theoretical study of different reaction conditions and previously unidentified side-products, we provide insight into multimetallic cluster synthesis reactions.

Paramagnetic NMR Shielding Tensors and Ring Currents: Efficient Implementation and Application to Heavy Element Compounds

We present an efficient implementation of paramagnetic NMR shielding tensors and shifts in a nonrelativistic and scalar-relativistic density functional theory framework. For the latter, we make use of the scalar exact two-component Hamiltonian in its local approximation, and generally we apply the well established (multipole-accelerated) resolution of the identity approximation and the seminumerical exchange approximation. The perturbed density matrix of a paramagnetic NMR shielding calculation is further used to study the magnetically induced current density and ring currents of open-shell systems as illustrated for [U@Bi₁₂]^3-. [U@Bi₁₂]^3- features delocalized highest occupied molecular orbitals and sustains a net diatropic ring current of ca. 18 nA/T through the Bi₁₂ torus similar to the all-metal aromatic heavy-element cluster [Th@Bi₁₂]^4-.