
==== Front
Nat Chem
Nat Chem
Nature Chemistry
1755-4330
1755-4349
Nature Publishing Group UK London

38760432
1530
10.1038/s41557-024-01530-z
Article
Supramolecular trapping of a cationic all-metal σ-aromatic {Bi4} ring
http://orcid.org/0000-0003-2375-7607
Yadav Ravi 12
http://orcid.org/0000-0001-8388-3096
Maiti Avijit 1
Schorpp Marcel 3
http://orcid.org/0000-0002-9107-4843
Graf Jürgen 4
http://orcid.org/0000-0001-5060-1689
Weigend Florian florian.weigend@chemie.uni-marburg.de

5
http://orcid.org/0000-0002-5253-419X
Greb Lutz greb@uni-heidelberg.de

1
1 https://ror.org/038t36y30 grid.7700.0 0000 0001 2190 4373 Anorganisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany
2 https://ror.org/01pe3t004 grid.462378.c 0000 0004 1764 2464 School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Thiruvananthapuram, India
3 https://ror.org/01eezs655 grid.7727.5 0000 0001 2190 5763 Universität Regensburg, Regensburg, Germany
4 https://ror.org/038t36y30 grid.7700.0 0000 0001 2190 4373 Organisch-Chemisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany
5 https://ror.org/01rdrb571 grid.10253.35 0000 0004 1936 9756 Fachbereich Chemie, Philipps-Universität Marburg, Marburg, Germany
17 5 2024
17 5 2024
2024
16 9 15231530
19 6 2023
5 4 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
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.

The synthesis of cationic all-metal aromatic systems without covalent functionalization remains an underexplored area in chemistry. Now a tetracationic [Bi4]4+ featuring all-metal σ-aromaticity has been stabilized through a supramolecular approach relying on dianionic calix[4]pyrrolato indiumbromide shells. This planar rhomboid represents the global minimum for 16 valence electron systems.

Subject terms

Organometallic chemistry
Molecular capsules
Chemical bonding
https://doi.org/10.13039/100010663 EC | EU Framework Programme for Research and Innovation H2020 | H2020 Priority Excellent Science | H2020 European Research Council (H2020 Excellent Science - European Research Council) 948708 Greb Lutz https://doi.org/10.13039/100005156 Alexander von Humboldt-Stiftung (Alexander von Humboldt Foundation) issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcMain

Seminal investigations by Hückel on organic π-conjugated molecules established the particular stability of planar rings with (4n + 2)π electrons1. Since then, aromaticity developed as a cornerstone in chemistry2,3. Over the past decades, the concept extended from π-electrons to delocalized σ-, δ- and φ-electrons in organic and inorganic compounds along different dimensions4–6. Aromaticity in pure metal compounds (all-metal aromaticity) sparked particular interest7,8, and numerous all-metal aromatics were predicted from quantum theory9. One of the simplest in silico examples is the σ-aromatic Li3+, the all-metal analogue of H3+ featuring two delocalized σ-electrons (Fig. 1a)10,11. Spectroscopic detections of the aromatic Al42− and antiaromatic Al44− in molecular beams represented seminal experimental work that indicated the generality of Hückel’s rule in all-metal ring structures (Fig. 1b)12,13. Cases such as the aromatic Au5Zn+ detected by mass spectrometry or the planar Ti3(CO)3 cluster observed in a solid argon matrix consolidated this perspective14,15. Aromatic metal clusters were also suggested in Zintl phases but appear restricted to anions16. For instance, the square [Bi4]2− has been regarded as an inorganic analogue of the prototypical π-aromatic [C4H4]2− (Fig. 1f)17. More recent density functional theory (DFT) computations challenged the aromaticity in [Bi4]2−, but revealed a small degree of antiaromatic character that contrasts the lighter [Pn4]2− congeners18. A case of π-aromaticity was determined in the heterometallic anion [Th@Bi12]4−, being only one example of the notable progress on bismuth clusters in recent years18,19.Fig. 1 Examples of aromatic rings.

a, Prototypical σ-aromatic cations. b, Gas phase observed aromatic and antiaromatic anions. c, Ligand-stabilized aromatic gallium ring. d,e Ligand-stabilized aromatic gold (d) and thorium rings (e). f, Aromatic Zintl anion. g, All-metal σ-aromatic formal [Bi4]4+ ring stabilized by supramolecular encapsulation.

Substantial insight on all-metal aromaticity was also gained with gallium-based three-, four- or five-membered ring systems with two π-electrons but that required stabilization by bulky terphenyl ligands (Fig. 1c)20–22. A recently reported thorium cluster extended σ-aromaticity towards actinides, while the electronic delocalization within the ring remained disputed (Fig. 1e)23–25. Aromatic systems in a cationic charge state are much less explored. Only mono- and di-cationic aromatic systems are known but have been restricted to compounds stabilized by large, covalent substituents26–30. For instance, a transition metal system with an aromatic Au3+ core coordinated by N-heterocyclic carbenes was presented (Fig. 1d)31. Unfortunately, any direct substitution perturbs the aromatic core’s electronic structure and impedes comparisons between experiment and the rich theoretical work on unperturbed ring systems9. Examples of cationic all-metal aromatic ring systems without covalent functionalization have not been reported. Notably, the highly electron-deficient nature of unsaturated cationic aromatics raised concerns about the chance of ever being isolated in the condensed phase32. Consequently, also the theoretical interest oriented more towards anions, while polycationic all-metal aromatics were rarely studied computationally8,9.

In this Article, we present a planar {Bi4} rhomboid trapped within two calix[4]pyrrolato indiumbromide shells (Fig. 1g). While supramolecular stabilization approaches proved effective for organic intermediates33–35, this strategy is far less developed for p-block elements36,37. The formally tetracationic [Bi4]4+ ring represents an isolable cationic all-metal aromatic ring without directly attached ligands. Interestingly, the [Bi4]4+ unit is isoelectronic to the gas-phase-observed Hückel-antiaromatic [Al4]4−.

Results and discussions

Synthesis and characterization

The reaction of the lithium salt of octaethylporphyrinogen [Li4·(thf)3EtCx] with BiCl3 in tetrahydrofurane (THF) resulted in the formation of [Li(thf)2EtCxBi] (1) in 80% yield (Fig. 2a). Single-crystal X-ray diffraction (scXRD) analysis of 1 revealed a Bi coordinated by four nitrogen atoms (Fig. 2b). The Li+ countercation is tetrahedrally coordinated by two thf, a nitrogen of one pyrrole ring and the α-C of the adjacent pyrrole ring. Despite the asymmetry in the solid state, the 1H-NMR spectrum of 1 in dichloromethane-d2 shows a singlet at δ = 6.28 ppm for the pyrrolyl β-hydrogen atoms, indicating the dynamic coordination of Li+ to the pyrrole rings. The Li+ can be replaced with non-coordinating PPh4+ by a metathesis reaction between 1 and PPh4Cl to obtain [PPh4][EtCxBi] (2), revealing a pseudo-C4v symmetric anion with the Bi adopting a non-VSEPR square pyramidal geometry (Fig. 2c; see Supplementary Section 2.3).Fig. 2 Synthesis and molecular structures of 1 and 2.

a–c, Synthesis of 1 and 2 (a) and molecular structure of 1 (b) and 2 (c) in solid state. dcm, dichloromethane; r.t., room temperature. Hydrogen atoms and PPh4+ (for 2) are omitted for clarity.

During our investigation of 2 as a transmetallation agent, we observed an unusual reactivity with InBr3. The reaction of 2 with InBr3 in dichloromethane resulted in a rapid colour change from light orange to dark red. After standing for 3 days at room temperature, dark-coloured crystals of 3 developed (Fig. 3f). The synthesis is reproducible, yielding up to 28% of isolated 3. scXRD analysis showed a planar rhomboid {Bi4} ring encapsulated by two indiumbromide-calix[4]pyrrolates, [(Bi4)(EtCxInBr)2] (Fig. 3b). The In–Br (2.5355(4) Å) bond lengths in 3 are in line with literature-known bromidoindinates In–Br (2.5541(4)–2.5615(4) Å)38. The In–N (average 2.240(3) Å) are in the range of pentacoordinated indium centres, for example, (bis)amidinate indiumbromide (2.1727(18)–2.2503(19) Å)39. Each indiumbromide-calix[4]pyrrolate can be understood as dianionic, rendering the four-membered Bi ring formally tetracationic, with a formal oxidation state of each Bi atom as +1. The Bi–Bi bond lengths of the rhomboid are 3.0223(3) Å and 3.0238(3) Å, which are in the range of literature-known Bi(I)–Bi(I) bonds (2.970 to 3.038 Å)40,41.Fig. 3 Synthesis and molecular structures of 3 and 4.

a, Synthesis of compound 3 and 4; *undefined stoichiometries of Bi metal formed. dcm, dichloromethane. b, The molecular structure of 3 in solid state. Hydrogen atoms are omitted for clarity. The {Bi4} ring is disordered over two positions. Two bismuth atoms are symmetry generated (inversion in the centre of the Bi4 rhomboid 1 − x, 1 − y and 1 − z; Supplementary Fig. 44). c, The molecular structure of 4 in solid state. Hydrogen atoms are omitted for clarity. Cyclopropyl ring is highlighted in red. d, Top view of the {Bi4} ring in 3 along with the Bi–Bi bond lengths. e, Side view of the {Bi4} ring in 3 showing the planarity. f, A microscopic photograph of the single crystals of 3.

The Bi1–Bi1′ distance (3.3772(5) Å) is substantially elongated, indicating only weak transannular interaction. The closest Bi–Cpyrrole distance in 3 is 2.982 Å, which is longer than reported Bi–C covalent bonds (around 2.26 Å)42–46. Despite the Bi–C distance in 3 being shorter than the sum of van der Waals radii (3.770 Å), these parameters indicate minor covalent interaction and justify the formal tetracationic classification as a [Bi4]4+ ring. This interpretation is further supported by bond analysis tools, as described in the Quantum chemical investigations section.

Compound 3 is insoluble in any solvent and could not be characterized in solution, but solid-state 13C- and 15N-MAS-NMR (MAS, magic-angle-spinning), infra-red and Raman spectroscopy are in line with the expected signatures and with the DFT-computed spectroscopic features. Elemental analysis of the isolated crystals of complex 3 gave carbon values consistently 8–10% lower than expected. We surmised the formation of bismuth metal covering the crystalline material as a source for this deviation. Indeed, powder XRD (pXRD) of the isolated solid showed the diffraction pattern of 3 (matching peaks with the calculated pattern; Supplementary Fig. 11), beside a species matching with metallic bismuth (Supplementary Fig. 12; http://www.crystallography.net/cod/9008576.html). pXRD analysis of the solid residue forming along with the crystals of 3 showed a diffraction pattern corresponding to pure bismuth metal only (Supplementary Fig. 13).

The planar, rhomboidal structure of [Bi4]4+ is surprising and contrasts with previous literature reports of four-membered Bi(I) systems. Indeed, the four-membered R4Bi4 structures exhibit a pronounced buckling in a butterfly shape40,41. In one of the first experimental reports on cationic bismuth clusters, the presence of a [Bi4]4+ was suggested from spectral studies on bismuth in liquid BiCl3, but the solid phase was consistently found to contain the [Bi9]5+ ion in a lower oxidation state47. As discussed in a later section, DFT computations identified the planar rhomboid as the global minimum structure also for the parent, free [Bi4]4+. However, before treating the electronic structure, the formation mechanism shall be discussed.

Since the deposition of Bi metal was observed during the process, we assumed that [Bi4]4+ might represent an intermediate in metal formation that becomes stabilized as the side product 3 in the presence of indiumbromide-calix[4]pyrrolates. But what serves as the reducing agent? The redox chemistry of calix[4]pyrrolates has been of interest within the context of transition metal complexes48–50 and recently been extended by us to complexes of antimony51. The oxidation by 2e− results in a so-called Δ-form, in which a cyclopropyl ring forms upon dearomatization of two pyrrole rings (Fig. 4a). Indeed, characterization of the supernatant of the reaction of 2 with InBr3 (Fig. 3a) revealed the formation of [ΔEtCxInBr] (4) and PPh4InBr4 as stoichiometric by-products by scXRD (Fig. 3c) and NMR spectroscopy (Supplementary Figs. 21–24; see Supplementary Section 1.2.3 for details). Hence, the formation of 4 revealed the calix[4]pyrrolato ligand as a source of the electrons for the reduction of Bi(III) over Bi(I) to elemental bismuth. The exact sequence of steps and inner-sphere versus outer-sphere electron transfer is impossible to distinguish, but the following is one plausible pathway. Transmetallation of Bi(III) with In(III) is furnishing [EtCxInBr]2− and [BiBr3−n]n+ (Fig. 4b). Dianionic [EtCxInBr]2− reduces [BiBr3−n]n+ to [BiIBr1−n]mn+ along with the formation of the two-electron oxidized [ΔEtCxInBr] (4).Fig. 4 Underlying redox chemistry, formation mechanism, control experiments and derivatives of 3.

a, Redox chemistry of the calix[4]pyrrole ligand. dcm, dichloromethane. b, Proposed reaction pathway for the formation of 3 and 4. c, In situ generation of compound 7 and reaction with BiCl3 for the synthesis of 3 and 4 to support the proposed mechanism. d, Reaction of 2 with AlCl3 and GaCl3. e, The molecular structure of the cation in 8Al in solid state. Hydrogen atoms and anion are omitted for clarity. *Bi metal is formed due to the overreduction of Bi(I) species via other reaction channels.

The released bromide ions are captured by InBr3, as confirmed by the stoichiometric formation of PPh4InBr4. Alternatively, InBr3 might activate BiBr3 during the reduction process by bromide abstraction. Thus, from the perspective of the reducing agent [EtCxInBr]2−, the formation of [ΔEtCxInBr] 4 might also be described as a Lewis acid-assisted oxidative C–C coupling. Ultimately, the Bi(I)-cationic cluster forms and gets trapped within remaining bowl-shaped dianionic [EtCxInBr]2−. This reaction illustrates the multi-purpose nature of calix[4]pyrrolato ligand as both a mild reductant and supramolecular host.

To check whether the transmetallation product [EtCxInBr]2− can reduce Bi(III) to Bi(I), we generated [PPh4]2[EtCxInBr] 7 in situ by an alternative synthesis and reacted with BiCl3 (Fig. 4c and Supplementary Section 1.2.4). Indeed, the formation of 3, 4, PPh4Cl and bismuth metal was observed, supporting the above-suggested pathway. Interestingly, [PPh4]2[EtCxInBr] 7 is thermally unstable and decomposes within 4 h in solution and within 5 days in solid state at room temperature. On the other hand, after encapsulating [Bi4]4+ in compound 3, [EtCxInBr]2− is stable at room temperature for months. Based on the redox potentials of 2 and 7 determined by cyclic voltammetry (Supplementary Section 1.2.6), 2 cannot be excluded to serve as an alternative reducing agent towards deliberated Bi(III). However, 2 did appear unreactive towards BiCl3, corroborating the idea that transmetallation is required to initiate the reaction sequence.

To confirm the generality of the formation of 3 and to obtain potential derivatives thereof, 2 was reacted with the lighter group 13 halides AlCl3 and GaCl3. Indeed, also in these cases, the formation of an encapsulated [Bi4]4+ species was observed, but the supramolecular cavity, this time composed of two monoanionic [EtCxBi]− originating from starting materials 2, with two additional AlCl<sub4</sub− or GaCl4− counteranions, respectively (8Al and 8Ga; Fig. 4d). At the same time, the corresponding transmetallated, oxidized products [ΔEtCxECl] 4Al and 4Ga and PPh4ECl4 (E = Al, Ga) were observed by NMR spectroscopy and high-resolution mass spectrometry (HR-MS) in the supernatant. The exclusive detection of oxidized 4Al and 4Ga (instead of corresponding oxidized Bi-containing compounds) further supports that the transmetallated dianions are favoured as reducing agents compared with 2. It is interesting to note that literature-known pentacoordinated calix[4]aluminates and gallates that are isostructural to the putatively involved dianions [EtCxAlX]2− and [EtCxGaX]2− do not feature the bowl-shaped structure of [EtCxInBr]2− but a ruffled conformation52,53 and are thus not suitable to trap [Bi4]4+. Corresponding trends were confirmed by the computation of related thermodynamics (Supplementary Section 3.6). Hence, supramolecular trapping requires a bowl-shaped species, which in this case represents [EtCxBi]−. Of note, the bond lengths inside the {Bi4} unit of a crystallographically characterized salt of a salt containing 8Al are similar to that in 3 (see comparison in Supplementary Section 2.3), despite the different charge state of the supramolecular host (−1 compared with −2 in 3). This observation indicates that the electronic structure in {Bi4} does not strongly depend on the trapping sphere and supports the above-chosen description and discussion as [Bi4]4+.

Quantum chemical investigations

For the analysis of the electronic structure, density functional calculations (Perdew–Burke–Ernzerhof (PBE) functional, basis sets of triple zeta valence quality) were performed for 3 and M4 rings (for details and references, see Supplementary Information). For all charged species, the conductor-like screening model was used. Natural population analysis on the experimental structure of 3 yields a charge of +0.69 |e| for Bi1 and Bi3 (Fig. 3) and +0.45 |e| for Bi2 and Bi4, and thus a total charge of +2.3 |e| for the Bi4 ring. Since these numbers are not observables but indicate directions, we perceive the idealized description as [Bi4]4+[L2−]2 with four Bi(I) as reasonable. From calculation and analysis of localized molecular orbitals (MOs), it is evident that the difference between ideal and calculated charge comes from a partial delocalization of eight pi-type bonds to empty orbitals of the Bi44+ ring (Supplementary Section 3.5). For further corroboration, a topological investigation of the density for the [Bi4]4+[L2−]2 interaction was performed and compared with that of the pentamethyl cyclopentadiene anion (Cp*) with Bi+, indicating a higher ionic character for the title compound due to the total density and the ratio of its curvatures at the bond critical points (BCPs). Moreover, energy decomposition analysis yields a predominantly ionic nature in bonding (80%, as discussed in Supplementary Section 3.5). Further support for this assignment is the structural parameters of the computed global minimum for the [Bi4]4+ without any ligand sphere, yielding a rhombus with 298/313/508 pm for Bi1–Bi2/Bi2–Bi4/Bi1–Bi3 that agrees with the encapsulated ring structure in 3 (301/332/502 pm for the computed and 302/337/502 pm for the scXRD structure).

For simplicity, we start the discussion of with the electronic structure of the isolated [Bi4]4+ rhombus. The four highest occupied MOs are shown in Fig. 5 (left).Fig. 5 Highest occupied MOs in Bi44+ and corresponding counterparts in 3.

Contours are drawn at ±0.04 a.u.

The π-type highest occupied molecular orbital (HOMO) is the bonding linear combination of the 6p orbitals perpendicular to the plane, transforming like b1u for the idealized D2h-symmetric structure. The following three orbitals are σ-type linear combinations of in-plane 6p orbitals, transforming like b2u, b1g and ag. The ag orbital is bonding between Bi2 and Bi4, and the b2u orbital is antibonding between Bi1 and Bi3, rationalizing the rhomboid distortion. The four MOs arising from the Bi-6s orbitals are much lower in energy with E − EHOMO amounting to −7.9/−8.5/−9.7/−11.2 eV.

Four similar MOs are found for the encapsulated system in 3, with HOMO and HOMO-1 being energetically interchanged (Fig. 5, right). These four MOs are intrinsically delocalized. If one applies a localization procedure to [Bi4]4+—we have chosen that of Pipek and Mezey as an example—one obtains an orbital that is practically identical to the HOMO (Supplementary Fig. 52). Consequently, the other three MOs result in three orbitals representing three-centre bonds distributed over the five Bi–Bi contacts. In such a situation, it is likely that the presence of a magnetic field will induce a ring current, which is an indication of aromaticity. Calculations of densities of induced magnetic ring currents from the response of the density to an external magnetic field perpendicular to the Bi4 ring (obtained with TURBOMOLE’s NMR module) via Biot–Savart’s law were done with the GIMIC tool (Fig. 6)54,55. The integration plane was chosen orthogonal to one of the four equivalent Bi–Bi bonds. Integration starts at the ring centre (green line in Fig. 6a, middle), and extends from 10 a.u. below to 10 a.u. above the ring plane. This way, a ring current of 9.1 nA T−1 is obtained for the isolated [Bi4]4+ (twice the area under the black curve in the current profile shown in Fig. 6e). The contribution of the π-type HOMO to the current can be estimated from the calculation of Bi46+, where this orbital is unoccupied. For this system, one gets 6.5 nA T−1; thus, the π-contribution to the current in [Bi4]4+ is roughly one quarter only.Fig. 6 Induced ring-current analysis.

a, Side and top view of absolute values of current densities in [Bi4]4+ (the green line indicates the orientation of the integration plane). b, Side view of absolute values of current densities in 3. c, Side view of absolute values of current densities in [EtCxInBr]24− (the vertical line indicates the z direction, a horizontal line is drawn at 2 bohr. d, The difference between the current densities in b and c. Contours are drawn at 0.025 a.u, blue (red) contours indicate paratropic/diatropic areas, that is, areas with currents flowing counter-clockwise or clockwise. e, Current (j) profiles for a–d along the z direction defined above.

For comparison, the current in benzene, a typical π-aromat, is calculated to 11.4 nA T−1, and when removing the π-electrons (MOs e1g and a2u) it vanishes (−0.1 nA T−1). Hence, electronic delocalization in the planar rhomboid of [Bi4]4+ is predominant in the σ-framework of the ring.

Does the situation change in the encapsulated ring in compound 3? The absolute values of the current density of 3 are shown in Fig. 6b, that of the dianionic ligands only in Fig. 6c and that of the difference of both in Fig. 6d. Evidently, the ring current is maintained in the presence of the ligands, albeit somewhat extended to them. Integration over a plane extending from 2 a.u. below to 2 a.u. above the ring plane, which is a plausible choice according to Fig. 6c (green horizontal line), yields a current strength of 14.1 nA T−1 for the entire system. The contribution from the ligands in this region is 1.6 nA T−1. The difference of both, 13.2 nA T−1, corresponds to the current arising from the current density of the Bi4 ring in the presence of the ligands and is even higher than that of the isolated [Bi4]4+. The increase of currents and current densities in the presence of ligands is also reflected by nucleus-independent shifts (see Supplementary Information for details). For the bare [Bi4]4+ unit, they amount to ‒14.5 ppm at the ring centre and to ‒15.0 ppm at 1 a.u. above, whereas for the entire system 3, one finds ‒34.2 and ‒31.6 ppm. Overall, the magnetic analyses indicate aromaticity in [Bi4]4+, within both the isolated and supramolecular trapped [Bi4]4+ unit. Indeed, in the experimental 13C-CP-MAS NMR (CP, cross-polarization) of 3, unusually down-field shifted signals (>180 ppm) were observed. These shifts were reproduced by DFT computations and assigned to the pyrrole α-carbon atoms closest to the bismuth ring. The unusual shifts disappear upon replacing the [Bi4]4+ with 4 Xe atoms (Supplementary Section 3.7). Hence, this observation might be a spectroscopic indication of all-metal aromaticity, but further experiments and derivatives are required for a conclusive statement.

To inspect the generality of the rhomboid structural motif, we compared [Bi4]4+ with isoelectronic Pb4 and [Tl4]4− and with their lighter homologues from rows 3p to 5p. All results hold when including spin–orbit coupling and changing the functional to PBE0. Details are documented in Supplementary Information; we summarize the main results here. For all 12 species, the global minimum is the rhombic ring with an electronic singlet state. Stationary points with quadrangular shape (triplets) or rectangular shape (singlets) were calculated as saddle points. Further local minima are tetrahedra (singlets), disfavoured by 62 kJ mol−1 [Tl4]4− to 317 kJ mol−1 ([P4]4+), and butterfly structures (triplets), disfavoured by 17 kJ mol−1 [Tl4]4− to 136 kJ mol−1 ([As4]4+). For the rhombic arrangements of [Al4]4−, [Ga4]4−, [In4]4− and [Tl4]4−, which show comparably small HOMO–LUMO gaps (where LUMO is the lowest unoccupied molecular orbital) for the singlet state, additional triplet state minima are found, but were disfavoured by 11 kJ mol−1 to 14 kJ mol−1. Ring currents for the rhombic structures range from 5.2 nA T−1 ([Tl4]4−) to 14.8 nA T−1 ([Al4]4−), reflecting moderate to substantial aromatic character. Hence, the aromatic rhombic structure is favoured for all 16-valence electron (VE) M4 species, including Al44−. This latter ion has been observed in the gas phase by mass spectrometry as Li3Al4− (ref. 13). But why has it been described as a rectangular antiaromat? Indeed, the rectangular topology of [Al4]4− becomes a minimum in the anisotropic sphere of four Li+ and the current is overall paratropic, −5.0 nA T−1, reflecting weak antiaromaticity. The latter comparison offers an interesting conclusion: within an anisotropic Coulomb field (for example, in Li3[Al4]−), the antiaromatic rectangular is a minimum for 16VE M4 but an aromatic rhomboid minimum within a spherical (more isotropic) charge sphere (for example, as Al44− in a dielectric continuum, or as [Bi4]4+ observed in isoelectronic 3).

Conclusions

Here, we describe the isolation and characterization of a cationic all-metal σ-aromatic ring, the formally tetracationic {Bi4}. The highly electron-deficient unit is trapped as compound 3 in a highly symmetric charge field of a cavity spanned by two π-electron-rich calix[4]pyrrolato units. The similarity of the bowl-shaped, dianionic calix[4]pyrrolato indinate to the versatile host corannulene is to be noted56. The reduction of a Bi(III) precursor to the Bi(I) containing 3 occurs along with the formation of elemental bismuth, and it can be speculated whether [Bi4]4+ represents a rare snapshot during metal formation57,58. Interestingly, the reduction cascade is initiated by adding a Lewis acid, while the reducing equivalents originate from the electron-rich ligand—contrasting with common strategies based on reducing agents. The planar rhomboid of [Bi4]4+ is identified as the preferred structure for all isoelectronic four-atomic 16VE p-block element systems, providing rich opportunities to extend experimental and theoretical directions in metal-cluster design. Further, this study indicates that changes in the charge distribution around isoelectronic ionic metal rings can influence the preference over aromatic (rhomboid 16VE systems in symmetric charge sphere, for example, 3) versus antiaromatic (rectangular in asymmetric charge sphere, for example, Li3[Al4]−) minima. These subtle effects might provide new opportunities to control electron delocalization or charge carrier mobility with implications for materials design. From a more fundamental perspective, these insights question if excessively deterministic discussions on aromaticity, non-aromaticity or antiaromaticity are always reasonable, or if the projection of Hückel’s model finds its limitation if leaving from the second-row elements towards the heavier analogues.

Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41557-024-01530-z.

Supplementary information

Supplementary Information Supplementary Figs. 1–55, discussion and Tables 1–13.

Supplementary Data 1 Crystallographic data for compound 1; CCDC reference no. 2264827.

Supplementary Data 2 Crystallographic data for compound 2; CCDC reference no. 2264830.

Supplementary Data 3 Crystallographic data for compound 3; CCDC reference no. 2264829.

Supplementary Data 4 Crystallographic data for compound 4; CCDC reference no. 2264828.

Supplementary Data 5 Crystallographic data for compound 8_AlOAl; CCDC reference no. 2334621.

Supplementary Data 6 Cartesian coordinates of all DFT computed structures.

Supplementary information

The online version contains supplementary material available at 10.1038/s41557-024-01530-z.

Acknowledgements

Financial support for this project was provided by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no 948708). We acknowledge support by the state of Baden-Württemberg through bwHPC and the German Research Foundation (DFG) through grant no. INST 40/575-1 FUGG (JUSTUS 2 cluster). We acknowledge D. Hartmann, H. Ruppert and D. Roth for collection of scXRD data, and M. Schmitt for advice in scXRD structure refinement. J. Zaumseil and S. Lindenthal are acknowledged for Raman spectroscopy measurements, D. Arian and J. Gross for mass spectrometric measurements and J. Klingeler for SQUID measurements.

Author contributions

R.Y. and L.G. conceived the project. R.Y. and A.M. conducted all experimental work and analysed the data. Structure solution and refinement for scXRD analysis was carried out by M.S., and solid-state MAS NMR measurements by J.G. Quantum chemical computations were performed and analysed by F.W. The manuscript was written by all authors. L.G. directed the project.

Peer review

Peer review information

Nature Chemistry thanks Alison Edwards, Moumita Majumdar and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Data availability

All data generated or analysed during this study are included in this manuscript and its supplementary information files. The structures of compounds 1, 2, 3, 4 and 8AlOAl were determined by scXRD. Crystallographic data for the structures reported have been deposited at the Cambridge Crystallographic Data Centre, under deposition numbers CCDC 2264827 (1), 2264830 (2), 2264829 (3), 2264828 (4) and 2334621 (8AlOAl). Copies of the data can be obtained free of charge via https://www.ccdc.cam.ac.uk/structures/. Details of the quantum chemical calculations can be found in Supplementary Information together with the relevant references to the methods used and information on their availability. The optimized Cartesian coordinates of all compounds mentioned explicitly in the main text or Supplementary Information are provided in a supplementary data file.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Ravi Yadav, Avijit Maiti.
==== Refs
References

1. Hückel E Quantentheoretische Beiträge zum Benzolproblem Z. Physik 1931 70 204 286 10.1007/BF01339530
Hückel, E. Quantentheoretische Beiträge zum Benzolproblem. Z. Physik 70, 204–286 (1931).10.1007/BF01339530
2. Schleyer PVR Introduction: aromaticity Chem. Rev. 2001 101 1115 1118 10.1021/cr0103221 11749368
Schleyer, P. V. R. Introduction: aromaticity. Chem. Rev. 101, 1115–1118 (2001).11749368 10.1021/cr0103221
3. Solà M Aromaticity rules Nat. Chem. 2022 14 585 590 10.1038/s41557-022-00961-w 35668210
Solà, M. Aromaticity rules. Nat. Chem. 14, 585–590 (2022).35668210 10.1038/s41557-022-00961-w
4. Bühl M Hirsch A Spherical aromaticity of fullerenes Chem. Rev. 2001 101 1153 1184 10.1021/cr990332q 11710216
Bühl, M. & Hirsch, A. Spherical aromaticity of fullerenes. Chem. Rev. 101, 1153–1184 (2001).11710216 10.1021/cr990332q
5. King RB Three-dimensional aromaticity in polyhedral boranes and related molecules Chem. Rev. 2001 101 1119 1152 10.1021/cr000442t 11710215
King, R. B. Three-dimensional aromaticity in polyhedral boranes and related molecules. Chem. Rev. 101, 1119–1152 (2001).11710215 10.1021/cr000442t
6. Minkin VI Minyaev RM Cyclic aromatic systems with hypervalent centers Chem. Rev. 2001 101 1247 1266 10.1021/cr990358h 11710220
Minkin, V. I. & Minyaev, R. M. Cyclic aromatic systems with hypervalent centers. Chem. Rev. 101, 1247–1266 (2001).11710220 10.1021/cr990358h
7. Boldyrev AI Wang L-S All-metal aromaticity and antiaromaticity Chem. Rev. 2005 105 3716 3757 10.1021/cr030091t 16218565
Boldyrev, A. I. & Wang, L.-S. All-metal aromaticity and antiaromaticity. Chem. Rev. 105, 3716–3757 (2005).16218565 10.1021/cr030091t
8. Popov IA Starikova AA Steglenko DV Boldyrev AI Usefulness of the σ-aromaticity and σ-antiaromaticity concepts for clusters and solid-state compounds Chem. Eur. J. 2018 24 292 305 10.1002/chem.201702035 28593736
Popov, I. A., Starikova, A. A., Steglenko, D. V. & Boldyrev, A. I. Usefulness of the σ-aromaticity and σ-antiaromaticity concepts for clusters and solid-state compounds. Chem. Eur. J. 24, 292–305 (2018).28593736 10.1002/chem.201702035
9. Mercero JM Boldyrev AI Merino G Ugalde JM Recent developments and future prospects of all-metal aromatic compounds Chem. Soc. Rev. 2015 44 6519 6534 10.1039/C5CS00341E 26058731
Mercero, J. M., Boldyrev, A. I., Merino, G. & Ugalde, J. M. Recent developments and future prospects of all-metal aromatic compounds. Chem. Soc. Rev. 44, 6519–6534 (2015).26058731 10.1039/C5CS00341E
10. Alexandrova AN Boldyrev AI σ-Aromaticity and σ-antiaromaticity in alkali metal and alkaline earth metal small clusters J. Phys. Chem. A 2003 107 554 560 10.1021/jp027008a
Alexandrova, A. N. & Boldyrev, A. I. σ-Aromaticity and σ-antiaromaticity in alkali metal and alkaline earth metal small clusters. J. Phys. Chem. A 107, 554–560 (2003).10.1021/jp027008a
11. Havenith RWA De Proft F Fowler PW Geerlings P σ-Aromaticity in H3+ and Li3+: insights from ring-current maps Chem. Phys. Lett. 2005 407 391 396 10.1016/j.cplett.2005.03.099
Havenith, R. W. A., De Proft, F., Fowler, P. W. & Geerlings, P. σ-Aromaticity in H3 + and Li3 +: insights from ring-current maps. Chem. Phys. Lett. 407, 391–396 (2005).10.1016/j.cplett.2005.03.099
12. Li X Kuznetsov AE Zhang H-F Boldyrev AI Wang L-S Observation of all-metal aromatic molecules Science 2001 291 859 861 10.1126/science.291.5505.859 11157162
Li, X., Kuznetsov, A. E., Zhang, H.-F., Boldyrev, A. I. & Wang, L.-S. Observation of all-metal aromatic molecules. Science 291, 859–861 (2001).11157162 10.1126/science.291.5505.859
13. Kuznetsov AE All-metal antiaromatic molecule: rectangular Al44− in the Li3Al4− anion Science 2003 300 622 625 10.1126/science.1082477 12714740
Kuznetsov, A. E. et al. All-metal antiaromatic molecule: rectangular Al4 4− in the Li3Al4− anion. Science 300, 622–625 (2003).12714740 10.1126/science.1082477
14. Tanaka H Neukermans S Janssens E Silverans RE Lievens P σ aromaticity of the bimetallic Au5Zn+ cluster J. Am. Chem. Soc. 2003 125 2862 2863 10.1021/ja029157c 12617638
Tanaka, H., Neukermans, S., Janssens, E., Silverans, R. E. & Lievens, P. σ aromaticity of the bimetallic Au5Zn+ cluster. J. Am. Chem. Soc. 125, 2862–2863 (2003).12617638 10.1021/ja029157c
15. Xu Q Jiang L Tsumori N cyclo-Ti3[η2(μ2-C,O)]3: a side-on-bonded polycarbonyl titanium cluster with potentially antiaromatic character Angew. Chem. Int. Ed. 2005 44 4338 4342 10.1002/anie.200500364
Xu, Q., Jiang, L. & Tsumori, N. cyclo-Ti3[η2(μ2-C,O)]3: a side-on-bonded polycarbonyl titanium cluster with potentially antiaromatic character. Angew. Chem. Int. Ed. 44, 4338–4342 (2005).10.1002/anie.200500364
16. Liu C Popov IA Chen Z Boldyrev AI Sun Z-M Aromaticity and antiaromaticity in Zintl clusters Chem. Eur. J. 2018 24 14583 14597 10.1002/chem.201801715 29774973
Liu, C., Popov, I. A., Chen, Z., Boldyrev, A. I. & Sun, Z.-M. Aromaticity and antiaromaticity in Zintl clusters. Chem. Eur. J. 24, 14583–14597 (2018).29774973 10.1002/chem.201801715
17. Cisar A Corbett JD Polybismuth anions. Synthesis and crystal structure of a salt of the tetrabismuthide(2-) ion, Bi42−. A basis for the interpretation of the structure of some complex intermetallic phases Inorg. Chem. 1977 16 2482 2487 10.1021/ic50176a011
Cisar, A. & Corbett, J. D. Polybismuth anions. Synthesis and crystal structure of a salt of the tetrabismuthide(2-) ion, Bi4 2−. A basis for the interpretation of the structure of some complex intermetallic phases. Inorg. Chem. 16, 2482–2487 (1977).10.1021/ic50176a011
18. Eulenstein AR Substantial π-aromaticity in the anionic heavy-metal cluster [Th@Bi12]4− Nat. Chem. 2021 13 149 155 10.1038/s41557-020-00592-z 33288891
Eulenstein, A. R. et al. Substantial π-aromaticity in the anionic heavy-metal cluster [Th@Bi12]4−. Nat. Chem. 13, 149–155 (2021).33288891 10.1038/s41557-020-00592-z
19. Heine J Peerless B Dehnen S Lichtenberg C Charge makes a difference: molecular ionic bismuth compounds Angew. Chem. Int. Ed. 2023 62 e202218771 10.1002/anie.202218771
Heine, J., Peerless, B., Dehnen, S. & Lichtenberg, C. Charge makes a difference: molecular ionic bismuth compounds. Angew. Chem. Int. Ed. 62, e202218771 (2023).10.1002/anie.202218771
20. Li X-W Pennington WT Robinson GH Metallic system with aromatic character. Synthesis and molecular structure of Na2[[(2,4,6-Me3C6H2)2C6H3]Ga]3 the first cyclogallane J. Am. Chem. Soc. 1995 117 7578 7579 10.1021/ja00133a045
Li, X.-W., Pennington, W. T. & Robinson, G. H. Metallic system with aromatic character. Synthesis and molecular structure of Na2[[(2,4,6-Me3C6H2)2C6H3]Ga]3 the first cyclogallane. J. Am. Chem. Soc. 117, 7578–7579 (1995).10.1021/ja00133a045
21. Twamley, B. & Power, P. P. Synthesis of the square-planar gallium species K2[Ga4(C6H3-2,6-Trip2)2] (Trip=C6H2-2,4,6-iPr3): the role of aryl–alkali metal ion interactions in the structure of gallium clusters. Angew. Chem. Int. Ed. 39, 3500–3503 (2000).
22. Kysliak O A planar five-membered aromatic ring stabilized by only two pi-electrons Angew. Chem. Int. Ed. 2022 61 e202206963 10.1002/anie.202206963
Kysliak, O. et al. A planar five-membered aromatic ring stabilized by only two pi-electrons. Angew. Chem. Int. Ed. 61, e202206963 (2022).10.1002/anie.202206963
23. Boronski JT A crystalline tri-thorium cluster with σ-aromatic metal–metal bonding Nature 2021 598 72 75 10.1038/s41586-021-03888-3 34425584
Boronski, J. T. et al. A crystalline tri-thorium cluster with σ-aromatic metal–metal bonding. Nature 598, 72–75 (2021).34425584 10.1038/s41586-021-03888-3
24. Cuyacot BJR Foroutan-Nejad C [{Th(C8H8)Cl2}3]2− is stable but not aromatic Nature 2022 603 E18 E20 10.1038/s41586-021-04319-z 35322245
Cuyacot, B. J. R. & Foroutan-Nejad, C. [{Th(C8H8)Cl2}3]2− is stable but not aromatic. Nature 603, E18–E20 (2022).35322245 10.1038/s41586-021-04319-z
25. Lin X Mo Y On the bonding nature in the crystalline tri-thorium cluster: core–shell syngenetic sigma-aromaticity Angew. Chem. Int. Ed. 2022 61 e202209658 10.1002/anie.202209658
Lin, X. & Mo, Y. On the bonding nature in the crystalline tri-thorium cluster: core–shell syngenetic sigma-aromaticity. Angew. Chem. Int. Ed. 61, e202209658 (2022).10.1002/anie.202209658
26. Gleiter R Structure and bonding in cyclic sulfur–nitrogen compounds—molecular orbital considerations Angew. Chem. Int. Ed. 1981 20 444 452 10.1002/anie.198104441
Gleiter, R. Structure and bonding in cyclic sulfur–nitrogen compounds—molecular orbital considerations. Angew. Chem. Int. Ed. 20, 444–452 (1981).10.1002/anie.198104441
27. Sun X Simler T Yadav R Köppe R Roesky PW A stable aromatic tetrasilacyclobutadiene dication J. Am. Chem. Soc. 2019 141 14987 14990 10.1021/jacs.9b08018 31507176
Sun, X., Simler, T., Yadav, R., Köppe, R. & Roesky, P. W. A stable aromatic tetrasilacyclobutadiene dication. J. Am. Chem. Soc. 141, 14987–14990 (2019).31507176 10.1021/jacs.9b08018
28. Freitag K The σ-aromatic clusters [Zn3]+ and [Zn2Cu]: embryonic brass Angew. Chem. Int. Ed. 2015 54 4370 4374 10.1002/anie.201410737
Freitag, K. et al. The σ-aromatic clusters [Zn3]+ and [Zn2Cu]: embryonic brass. Angew. Chem. Int. Ed. 54, 4370–4374 (2015).10.1002/anie.201410737
29. Jin L Trinuclear gold clusters supported by cyclic (alkyl)(amino)carbene ligands: mimics for gold heterogeneous catalysts Angew. Chem. Int. Ed. 2014 53 9059 9063 10.1002/anie.201404665
Jin, L. et al. Trinuclear gold clusters supported by cyclic (alkyl)(amino)carbene ligands: mimics for gold heterogeneous catalysts. Angew. Chem. Int. Ed. 53, 9059–9063 (2014).10.1002/anie.201404665
30. Blanchard S Synthesis of triangular tripalladium cations as noble-metal analogues of the cyclopropenyl cation Angew. Chem. Int. Ed. 2014 53 1987 1991 10.1002/anie.201310204
Blanchard, S. et al. Synthesis of triangular tripalladium cations as noble-metal analogues of the cyclopropenyl cation. Angew. Chem. Int. Ed. 53, 1987–1991 (2014).10.1002/anie.201310204
31. Robilotto TJ Bacsa J Gray TG Sadighi JP Synthesis of a trigold monocation: an isolobal analogue of [H3]+ Angew. Chem. Int. Ed. 2012 51 12077 12080 10.1002/anie.201206712
Robilotto, T. J., Bacsa, J., Gray, T. G. & Sadighi, J. P. Synthesis of a trigold monocation: an isolobal analogue of [H3]+. Angew. Chem. Int. Ed. 51, 12077–12080 (2012).10.1002/anie.201206712
32. Hoffmann R The many guises of aromaticity Am. Sci. 2015 103 18 22 10.1511/2015.112.18
Hoffmann, R. The many guises of aromaticity. Am. Sci. 103, 18–22 (2015).10.1511/2015.112.18
33. Kaphan DM Levin MD Bergman RG Raymond KN Toste FD A supramolecular microenvironment strategy for transition metal catalysis Science 2015 350 1235 1238 10.1126/science.aad3087 26785485
Kaphan, D. M., Levin, M. D., Bergman, R. G., Raymond, K. N. & Toste, F. D. A supramolecular microenvironment strategy for transition metal catalysis. Science 350, 1235–1238 (2015).26785485 10.1126/science.aad3087
34. Wang K Jordan J Hu XY Wang L Supramolecular strategies for controlling reactivity within confined nanospace Angew. Chem. Int. Ed. 2020 59 13712 13721 10.1002/anie.202000045
Wang, K., Jordan, J., Hu, X. Y. & Wang, L. Supramolecular strategies for controlling reactivity within confined nanospace. Angew. Chem. Int. Ed. 59, 13712–13721 (2020).10.1002/anie.202000045
35. Yu Y Yang JM Rebek J Molecules in confined spaces: reactivities and possibilities in cavitands Chem 2020 6 1265 1274 10.1016/j.chempr.2020.04.014
Yu, Y., Yang, J. M. & Rebek, J. Molecules in confined spaces: reactivities and possibilities in cavitands. Chem 6, 1265–1274 (2020).10.1016/j.chempr.2020.04.014
36. Rupar PA Staroverov VN Baines KM A cryptand-encapsulated germanium(II) dication Science 2008 322 1360 1363 10.1126/science.1163033 19039131
Rupar, P. A., Staroverov, V. N. & Baines, K. M. A cryptand-encapsulated germanium(II) dication. Science 322, 1360–1363 (2008).19039131 10.1126/science.1163033
37. Mal P Breiner B Rissanen K Nitschke JR White phosphorus is air-stable within a self-assembled tetrahedral capsule Science 2009 324 1697 1699 10.1126/science.1175313 19556504
Mal, P., Breiner, B., Rissanen, K. & Nitschke, J. R. White phosphorus is air-stable within a self-assembled tetrahedral capsule. Science 324, 1697–1699 (2009).19556504 10.1126/science.1175313
38. Schwamm RJ Anker MD Lein M Coles MP Fitchett CM Indyllithium and the indyl anion [InL]−: heavy analogues of N-heterocyclic carbenes Angew. Chem. Int. Ed. 2018 57 5885 5887 10.1002/anie.201802444
Schwamm, R. J., Anker, M. D., Lein, M., Coles, M. P. & Fitchett, C. M. Indyllithium and the indyl anion [InL]−: heavy analogues of N-heterocyclic carbenes. Angew. Chem. Int. Ed. 57, 5885–5887 (2018).10.1002/anie.201802444
39. Banerjee S Amidinate based indium(III) monohalides and β-diketiminate stabilized In(II)–In(II) bond: synthesis, crystal structure, and computational study Dalton Trans. 2020 49 14231 14236 10.1039/D0DT03161E 33026029
Banerjee, S. et al. Amidinate based indium(III) monohalides and β-diketiminate stabilized In(II)–In(II) bond: synthesis, crystal structure, and computational study. Dalton Trans. 49, 14231–14236 (2020).33026029 10.1039/D0DT03161E
40. Balazs, L., Breunig Hans, J., Lork, E. & Silvestru, C. Low-valent organobismuth compounds with intramolecular coordination: cyclo-R3Bi3, cyclo-R4Bi4, RBi[W(CO)5]2, and R4Bi2 [R=2-(Me2NCH2)C6H4]. Eur. J. Inorg. Chem. 2003, 1361–1365 (2003).
41. Balázs L Breunig HJ Organometallic compounds with Sb–Sb or Bi–Bi bonds Coord. Chem. Rev. 2004 248 603 621 10.1016/j.ccr.2003.12.001
Balázs, L. & Breunig, H. J. Organometallic compounds with Sb–Sb or Bi–Bi bonds. Coord. Chem. Rev. 248, 603–621 (2004).10.1016/j.ccr.2003.12.001
42. Breunig HJ Rösler R Lork E The first organobismuth rings: (RBi)3 and (RBi)4, R=(Me3Si)2CH Angew. Chem. Int. Ed. 1998 37 3175 3177 10.1002/(SICI)1521-3773(19981204)37:22<3175::AID-ANIE3175>3.0.CO;2-M
Breunig, H. J., Rösler, R. & Lork, E. The first organobismuth rings: (RBi)3 and (RBi)4, R=(Me3Si)2CH. Angew. Chem. Int. Ed. 37, 3175–3177 (1998).10.1002/(SICI)1521-3773(19981204)37:22<3175::AID-ANIE3175>3.0.CO;2-M
43. Breunig HJ Organometallic compounds with homonuclear bonds between bismuth atoms, 70 years after Paneth’ report on the violet dimethyl bismuth compound Z. Anorg. Allg. Chem. 2005 631 621 631 10.1002/zaac.200400476
Breunig, H. J. Organometallic compounds with homonuclear bonds between bismuth atoms, 70 years after Paneth’ report on the violet dimethyl bismuth compound. Z. Anorg. Allg. Chem. 631, 621–631 (2005).10.1002/zaac.200400476
44. Pyykkö P Atsumi M Molecular double-bond covalent radii for elements Li–E112 Chem. Eur. J. 2009 15 12770 12779 10.1002/chem.200901472 19856342
Pyykkö, P. & Atsumi, M. Molecular double-bond covalent radii for elements Li–E112. Chem. Eur. J. 15, 12770–12779 (2009).19856342 10.1002/chem.200901472
45. Benjamin SL Karagiannidis L Levason W Reid G Rogers MC Hybrid dibismuthines and distibines: preparation and properties of antimony and bismuth oxygen, sulfur, and nitrogen donor ligands Organometallics 2011 30 895 904 10.1021/om1010148
Benjamin, S. L., Karagiannidis, L., Levason, W., Reid, G. & Rogers, M. C. Hybrid dibismuthines and distibines: preparation and properties of antimony and bismuth oxygen, sulfur, and nitrogen donor ligands. Organometallics 30, 895–904 (2011).10.1021/om1010148
46. Groom CR Bruno IJ Lightfoot MP Ward SC The Cambridge Structural Database Acta Crystallogr. Sect. B 2016 72 171 179 10.1107/S2052520616003954
Groom, C. R., Bruno, I. J., Lightfoot, M. P. & Ward, S. C. The Cambridge Structural Database. Acta Crystallogr. Sect. B 72, 171–179 (2016).10.1107/S2052520616003954
47. Hershaft A Corbett JD The crystal structure of bismuth subchloride. Identification of the ion Bi95+ Inorg. Chem. 1963 2 979 985 10.1021/ic50009a025
Hershaft, A. & Corbett, J. D. The crystal structure of bismuth subchloride. Identification of the ion Bi9 5+. Inorg. Chem. 2, 979–985 (1963).10.1021/ic50009a025
48. De Angelis S Solari E Floriani C Chiesi-Villa A Rizzoli C Oxidation of metal-meso-octaethylporphyrinogen complexes leading to novel oxidized forms of porphyrinogen other than porphyrins. 1. The redox chemistry of nickel(II)- and copper(II)-meso-octaethylporphyrinogen complexes occurring with the formation and cleavage of a cyclopropane unit J. Am. Chem. Soc. 1994 116 5691 5701 10.1021/ja00092a021
De Angelis, S., Solari, E., Floriani, C., Chiesi-Villa, A. & Rizzoli, C. Oxidation of metal-meso-octaethylporphyrinogen complexes leading to novel oxidized forms of porphyrinogen other than porphyrins. 1. The redox chemistry of nickel(II)- and copper(II)-meso-octaethylporphyrinogen complexes occurring with the formation and cleavage of a cyclopropane unit. J. Am. Chem. Soc. 116, 5691–5701 (1994).10.1021/ja00092a021
49. Bachmann J Nocera DG Multielectron chemistry of zinc porphyrinogen: a ligand-based platform for two-electron mixed valency J. Am. Chem. Soc. 2004 126 2829 2837 10.1021/ja039617h 14995200
Bachmann, J. & Nocera, D. G. Multielectron chemistry of zinc porphyrinogen: a ligand-based platform for two-electron mixed valency. J. Am. Chem. Soc. 126, 2829–2837 (2004).14995200 10.1021/ja039617h
50. Bachmann J Nocera DG Multielectron redox chemistry of iron porphyrinogens J. Am. Chem. Soc. 2005 127 4730 4743 10.1021/ja043132r 15796540
Bachmann, J. & Nocera, D. G. Multielectron redox chemistry of iron porphyrinogens. J. Am. Chem. Soc. 127, 4730–4743 (2005).15796540 10.1021/ja043132r
51. Schorpp M Yadav R Roth D Greb L Calix[4]pyrrolato stibenium: Lewis superacidity by antimony(III)–antimony(V) electromerism Angew. Chem. Int. Ed. 2022 61 e202207963 10.1002/anie.202207963
Schorpp, M., Yadav, R., Roth, D. & Greb, L. Calix[4]pyrrolato stibenium: Lewis superacidity by antimony(III)–antimony(V) electromerism. Angew. Chem. Int. Ed. 61, e202207963 (2022).10.1002/anie.202207963
52. Ebner F Wadepohl H Greb L Calix[4]pyrrole aluminate: a planar tetracoordinate aluminum(III) anion and its unusual Lewis acidity J. Am. Chem. Soc. 2019 141 18009 18012 10.1021/jacs.9b10628 31665601
Ebner, F., Wadepohl, H. & Greb, L. Calix[4]pyrrole aluminate: a planar tetracoordinate aluminum(III) anion and its unusual Lewis acidity. J. Am. Chem. Soc. 141, 18009–18012 (2019).31665601 10.1021/jacs.9b10628
53. Sigmund LM Engels E Richert N Greb L Calix[4]pyrrolato gallate: square planar-coordinated gallium(iii) and its metal–ligand cooperative reactivity with CO2 and alcohols Chem. Sci. 2022 13 11215 11220 10.1039/D2SC03054C 36320463
Sigmund, L. M., Engels, E., Richert, N. & Greb, L. Calix[4]pyrrolato gallate: square planar-coordinated gallium(iii) and its metal–ligand cooperative reactivity with CO2 and alcohols. Chem. Sci. 13, 11215–11220 (2022).36320463 10.1039/D2SC03054C
54. Jusélius J Sundholm D Gauss J Calculation of current densities using gauge-including atomic orbitals J. Chem. Phys. 2004 121 3952 3963 10.1063/1.1773136 15332941
Jusélius, J., Sundholm, D. & Gauss, J. Calculation of current densities using gauge-including atomic orbitals. J. Chem. Phys. 121, 3952–3963 (2004).15332941 10.1063/1.1773136
55. Sundholm D Dimitrova M Berger RJF Current density and molecular magnetic properties Chem. Commun. 2021 57 12362 12378 10.1039/D1CC03350F
Sundholm, D., Dimitrova, M. & Berger, R. J. F. Current density and molecular magnetic properties. Chem. Commun. 57, 12362–12378 (2021).10.1039/D1CC03350F
56. Nestoros E Stuparu MC Corannulene: a molecular bowl of carbon with multifaceted properties and diverse applications Chem. Commun. 2018 54 6503 6519 10.1039/C8CC02179A
Nestoros, E. & Stuparu, M. C. Corannulene: a molecular bowl of carbon with multifaceted properties and diverse applications. Chem. Commun. 54, 6503–6519 (2018).10.1039/C8CC02179A
57. Ecker A Weckert E Schnöckel H Synthesis and structural characterization of an AI77 cluster Nature 1997 387 379 381 10.1038/387379a0
Ecker, A., Weckert, E. & Schnöckel, H. Synthesis and structural characterization of an AI77 cluster. Nature 387, 379–381 (1997).10.1038/387379a0
58. Schnöckel H Structures and properties of metalloid Al and Ga clusters open our eyes to the diversity and complexity of fundamental chemical and physical processes during formation and dissolution of metals Chem. Rev. 2010 110 4125 4163 10.1021/cr900375g 20540559
Schnöckel, H. Structures and properties of metalloid Al and Ga clusters open our eyes to the diversity and complexity of fundamental chemical and physical processes during formation and dissolution of metals. Chem. Rev. 110, 4125–4163 (2010).20540559 10.1021/cr900375g
