DOI:
10.1039/D3RA06603G
(Paper)
RSC Adv., 2023,
13, 34224-34229
The curious case of the crystalline tri-thorium cluster: cyclic delocalization without aromatic stabilization?
Received
27th September 2023
, Accepted 17th November 2023
First published on 22nd November 2023
Abstract
Actinides have been known to form extremely weak homonuclear bonds with their d-type orbitals, and one should therefore expect the superposition of cyclic resonance forms containing such bonds to bring rather marginal aromatic stabilization to the system, if any. It is for this very reason that the discovery of the cyclically delocalized Th3 σ-bonding in the crystalline cluster isolated by Liddle and co-workers has sparked such vigorous discussion on the actual role of molecular aromaticity on the periphery of the periodic table. It has recently been argued that the tri-thorium ring at the heart of the cluster features considerable aromatic stabilization energy comparable to the heterocyclic π-aromatic rings such as thiophene. However, previous investigations involved highly ionized model clusters like Th3Cl64+ or Th310+ in which aromatic stabilization associated with the cyclic delocalization of charge is dramatically exaggerated. In this work we investigate the model tri-thorium clusters at different geometries and ionization states to show that cyclic delocalization of electrons in the isolated crystalline cluster may be associated with rather marginal σ-aromatic stabilization energy thus strongly suggesting its non-aromatic character.
Introduction
When Kekulé first pictured the cyclic structure of benzene,1 hardly anyone could have imagined that electron delocalization and aromaticity may play a central role in an incredible variety of inorganic rings and clusters containing s-, d- and even f-block metal atoms.2 Very recently, Liddle and co-workers isolated under normal experimental conditions the crystalline actinide cluster [{Th(η8-C8H8)(μ3-Cl)2}3{K(THF)2}2]∞ containing at its heart the Th3 ring with a pair of uniformly distributed electrons.3 This discovery has been hailed for extending the range of aromaticity to the record seventh row of the periodic table. Also, it has been suggested that aromaticity can be an important factor in the design of stable complexes with actinide–actinide bonds, which have previously been reported in the literature to be weak and localized.4 However, the follow-up investigations revealed that both experimental and computational data (especially results of aromaticity assessment based on the magnetic criteria)5 may not provide solid proof for the existence of the aromatic Th3 bonding6,7 and the unique multicenter (delocalized) charge-shift bonding (ThCl2)3 was shown to be the actual factor determining high symmetry (D3h) and exceptional thermodynamical stability of the crystalline actinide cluster.7 Lin and Mo provided further computational data supporting the charge-shift bonding in the actinide cluster.8 Moreover, the authors found that the overall charge transfer between Th and Cl atoms is effectively ten times more stabilizing than the electron delocalization in the Th3 ring.8 Lin and Mo linked the latter with the stabilizing effect resulting from the periodic boundary conditions (PBC) for the electron distribution,9 and made the conclusion that the Th3 bonding is “truly delocalized and σ-aromatic”.8 Similar conclusion was recently drawn by Tomeček et al.,10 who utilized the entire arsenal of different sophisticated computational methods to demonstrate that Th3 contains cyclically delocalized electrons and should be regarded σ-aromatic. In both works,8,10 however, the simplified ionized model cluster Th3Cl64+ was investigated, in which electron binding energy associated with the highest occupied molecular orbital (HOMO) is incomparably higher in magnitude than in the neutral crystalline cluster. This leads to dramatic exaggeration of σ-bonding between thorium atoms (a direct manifestation of which is contraction of the Th–Th bond length by about 0.3 Å) and makes the conclusions by Lin and Mo,8 and Tomeček et al.,10 not directly transferable to the Th3 core found in the experimentally isolated actinide cluster.3
In this work we show how separation of charge in the ionized model tri-thorium clusters investigated previously by different authors affects aromaticity predictions. Furthermore, we demonstrate that the tri-thorium σ-bonding in the cyclic isomer brings rather marginal extra stabilization compared to the linear one, which suggests non-aromatic character in accordance with the IUPAC recommendation.11
The covalent resonance energy issue in Th3
To quantitatively assess aromatic stabilization in the Th3 core of the isolated crystalline cluster, Lin and Mo utilized the highly ionized model Th310+ derived from the Th3Cl64+ cage at its equilibrium geometry (with the Th–Th bond length of 3.684 Å, thus being much shorter from the experimentally measured 3.991 Å).8 In both highly charged models, the highest-occupied molecular orbitals (HOMO) visually resemble the corresponding HOMO in the original model cluster investigated by Liddle and co-workers.3 The authors used a sophisticated computational method called the block-localized wavefunction (BLW) within the framework of the density functional theory (DFT)12 to quantitatively assess the covalent resonance stabilization energy associated with the delocalization of σ-bond in the Th310+ ring, ΔEcovRE (Fig. 1a).8 Next, by comparison of ΔEcovRE for cyclic (D3h) and linear (D∞h) isomers of Th310+ the authors found the extra cyclic resonance energy (ECRE) to be equal 18.7 kcal mol−1 (Fig. 1b). This value falls between the corresponding values for the archetypical σ-aromatic H3+ (31.9 kcal mol−1) and non-aromatic Li3+ (0.2 kcal mol−1), which prompted Lin and Mo to make the conclusion that the tri-thorium ring is considerably σ-aromatic; for comparison, the reported aromatic stabilization energy in thiophene is about 18.6 kcal mol−1.13
![image file: d3ra06603g-f1.tif](/image/article/2023/RA/d3ra06603g/d3ra06603g-f1.gif) |
| Fig. 1 (a) Definition of the ECRE for Th310+ proposed by Lin and Mo. (b) The results of the BLW-DFT based calculations of ECRE with the corresponding HOMO isosurfaces (visualized at τ = 0.07e). Numbers in brackets represent degree of charge localization based on the molecular orbital coefficients. The ECRE calculations were performed using the same methods and software as used in the original work by Lin and Mo.8 | |
Unfortunately, the definition of ECRE (Fig. 1a) by Lin and Mo, and the resulting prediction of aromatic stabilization energy suffer from several important issues. First, the authors assume that electronic structures of cyclic (D3h) and linear (D∞h) isomers of Th310+ are predominated by the resonance of covalent forms. However, even a cursory look at the shapes of HOMOs in both isomers reveals that ionic forms prevail over the covalent ones with about 84–91% of charge localized on the central Th atoms. Consequently, the calculated resonance stabilization energy of 21.22 kcal mol−1 represents the effect of charge polarization rather than σ-bonding delocalization, which makes it not directly comparable with the covalent resonance energy of the cyclic isomer (39.88 kcal mol−1). Thus, the energy gap of 18.7 kcal mol−1 calculated by Lin and Mo cannot be strictly interpreted as the ECRE. Second, this particular result cannot be strictly interpreted as the σ-aromatic stabilization energy since it comes down with dramatic and systematic overestimation of the PBC effect.9 This is because the total lengths of the linear (L = 2 × RTh–Th) and cyclic (L = 3 × RTh–Th) isomer of Th310+ are different (as well as the total number of Th–Th bonds), which prevents direct assessment of the effect of PBC on their relative stability.9 Lin and Mo, being aware of this potential issue, also compared the resonance energies in the cyclic (D3h) form of Th310+ with the linear (D∞h) form of Th414+, and found the ECRE to be equal 28.8 kcal mol−1.8 But, despite having the same number of thorium–thorium bonds, the internal potential from nuclei experienced by delocalized electrons in both cases is fundamentally different making the calculated covalent resonance stabilization energies incomparable.
The charge separation issue in Th3
Aside from the issues associated with the definition of ECRE, the aromaticity predictions made for such extremely ionized model clusters like Th3Cl64+ and Th310+ may be not directly transferable to the tri-thorium core in the isolated crystalline cluster. This is because separation of charge between the Th3Cl6 cage and the cycloocta-1,3,5,7-tetraene (COT) ligands affects energy levels in the linear and cyclic isomers to a different extent, and, as will be shown in the following, the magnitude of the extra cyclic delocalization energy increases proportionally with the overall charge of the system. Moreover, the calculated covalent bond order14 between each Th atom and the COT ligand is equal to 0.96, thus marking sharing of the 7s2 electrons (Th) with the π-system of COT as an important component of chemical bonding along with strong electrovalent interactions.
To illustrate to what extent separation of charge affects the tri-thorium bonding in the core of the crystalline cluster, we performed the symmetry-constrained relativistic calculations of the 3c2e-orbital energy levels in Th310+, Th3Cl64+, Th3, and the original model by investigated by Liddle and co-workers (Fig. 2).3 It should be explained at this point that in the neutral cluster Th3 the interference of the 3c2e σ-type HOMO−2 (−3.19 eV) with a pair of the degenerate HOMO and HOMO−1 (−2.79 eV) does not affect the resulting distribution of charge (D3h symmetry), and, more importantly, it has no effect on the electron binding in the tri-thorium core. Indeed, a unitary transformation of the Kohn–Sham orbitals HOMO−2, HOMO−1, and HOMO into the natural bond orbitals (NBO)15 representation leads a subset of 3 doubly-occupied and fully localized (2c2e) σ-type bonds with the NBO energies equal to the original HOMO−2 energy level (−3.19 eV). This fact, as well as the observations reported below, justify the use of HOMO−2 in the neutral Th3 as an interesting alternative to the charged models of the tri-thorium bonding investigated by Lin and Mo.8
![image file: d3ra06603g-f2.tif](/image/article/2023/RA/d3ra06603g/d3ra06603g-f2.gif) |
| Fig. 2 Isosurfaces (τ = 0.07e) of the 3c2e molecular orbitals in selected charged and neutral model tri-thorium clusters with the corresponding total one-electron energies as well as the energy contributions from 6d and 7s orbitals (Th). Method: ωB97X/jorge-TZP-DKH with the Douglas–Kroll–Hess 2nd-order Hamiltonian. | |
The results of relativistic calculations clearly show that energy levels associated with the three-center two-electron Th3 σ-bonding in charged clusters are essentially different from the corresponding MO energies in the neutral ones. Moreover, orbitals in charged systems are much more inhomogeneous due to significant contribution of the atomic orbitals 7s (especially in Th310+). Also, despite visual similarity between HOMOs in Th3Cl64+ and the model by Liddle and co-workers,3 the overlap of three 6d orbitals in the former is associated with about six times higher electron binding energy. This results in significantly shorter Th–Th bond lengths (3.684 Å) compared to the average value from the experiment (3.991 Å) and calculations by Liddle and co-workers (4.039 Å).3 Thus, it is clear that high ionization in the cationic Th3 models result in dramatic exaggeration of bonding strength, which make them hardly comparable with the experimentally validated model of the crystalline cluster.3 Consequently, the results of theoretical investigations involving Th310+ and Th3Cl64+ cannot be used to justify σ-aromaticity in the isolated actinide cluster.
Aromatic stabilization from first principles
Although, to the best of our knowledge, there is no reliable and universal method to quantitatively assess the effect of aromatic stabilization in all-metal clusters, the three-center two-electron bonding is a particularly special case in which stabilizing effects of the periodic boundary conditions on the electron distribution can be derived directly from the MO energy diagrams for linear and cyclic isomers. This is mainly because within the 1-electron approximation the effect of PBC on such 3c2e bonding is directly proportional to the corresponding effect on the 3c2e-type MO itself (Fig. 3). Accordingly, by comparison of the MO energy levels of the equal-length linear (C∞v) and cyclic (D3h) forms of Th3, Th34+, and Th310+, it can easily be found that only in the cases of ionic clusters the isomers with cyclically delocalized charge have orbital energies lower than the C∞v isomers (with the electrons shared only by two thorium atoms) (Fig. 4). This means that only the charged model clusters can be considered effectively stabilized by the PBC, and hence feature aromatic stabilization.9 Contrariwise, in the neutral model cluster Th3 the corresponding 3c2e orbital in the cyclic isomer clearly makes no profits from the PBC as it has the same orbital energy as its linear counterpart with much more localized charge. It should be emphasized that the orbital energy differences between isomers C∞v and D3h (0.00 eV, −3.68 eV, and −11.37 eV) correlate linearly (R2 = 0.993) with the corresponding total charges of the model tri-thorium clusters (0, +4, +10), while lack of correlation is observed for the corresponding orbital energy differences between isomers D∞h and D3h. This is a direct manifestation of the foregoing incompatibility of these isomers in the context of assessment of the PBC effects due to different total lengths, L.9 Nevertheless, the fact that localized isomer C∞v has much lower orbital energy than the contracted isomer D∞h with more extended charge delocalization clearly shows that the 6d atomic orbitals are reluctant to σ-conjugation at this distance. This, in turn, is in full agreement with the conclusions from our previous study.7
![image file: d3ra06603g-f3.tif](/image/article/2023/RA/d3ra06603g/d3ra06603g-f3.gif) |
| Fig. 3 The correlation of the changes in the one-electron HOMO energy and the total energy from the self-consistent field calculations with the BLW constraints in the model cluster Th310+. The calculations were performed using the same methods and software as used in the original work by Lin and Mo.8 | |
![image file: d3ra06603g-f4.tif](/image/article/2023/RA/d3ra06603g/d3ra06603g-f4.gif) |
| Fig. 4 Comparison of the absolute and relative KS-MO energy levels in selected cyclic and acyclic isomers of the simplified neutral and charged tri-thorium clusters at the experimental geometry (RTh–Th = 3.991 Å). Method: ωB97X/jorge-TZP-DKH with the Douglas–Kroll–Hess 2nd-order relativistic Hamiltonian. | |
It should also be noted that at distances shorter than 3.991 Å the overlap of 6d orbitals in the neutral tri-thorium core would be more effective, and weak but noticeable σ-aromatic stabilization could be theoretically possible (Fig. 5). However, one must realize that even a small contraction of the Th–Th bond length in the isolated crystalline actinide cluster implies size reduction of the entire Th3Cl6 cage and the resulting enhancement of the Pauli repulsion between lone pairs of chlorine atoms crowed in the relatively small volume. This has been argued to have a negative impact on thermodynamic stability of the entire cluster.7 Thus, the strong multicenter charge-shift bonding in the Th3Cl6 cage seems to stretch the small ring of weakly bonded thorium atoms to the extent at which the extra cyclic stabilization is no more effective.
![image file: d3ra06603g-f5.tif](/image/article/2023/RA/d3ra06603g/d3ra06603g-f5.gif) |
| Fig. 5 The stabilizing and destabilizing effects of the PBC on the 3c2e σ-bonding in the neutral model cluster Th3 at different bond lengths, RTh–Th. Method: ωB97X/jorge-TZP-DKH with the Douglas–Kroll 2nd-order relativistic Hamiltonian (numbers in brackets refer to the corresponding relativistic calculations including the spin–orbit coupling). | |
Conclusions
The crystalline tri-thorium cluster represents an interesting and important case in which cyclic delocalization of electrons inside the all-metal core is not associated with aromatic stabilization, but is dictated by high symmetry of the network of surrounding Th–Cl bonds.7 In fact, the charge transfer between thorium and chlorine atoms makes the multicenter charge-shift bonding in the entire Th3Cl6 cage particularly strong.8 But, at the same time it impairs the direct σ-bonding between thorium atoms to the extent at which the Th3 ring exhibits no extra cyclic stabilization compared to its linear and more localized isomer (C∞v). Thus, in accordance with the very first sentence of the recommendation by IUPAC,11
‘Aromaticity is the concept of spatial and electronic structure of cyclic molecular systems displaying the effects of cyclic electron delocalization which provide for their enhanced thermodynamic stability relative to acyclic structural analogues (…)’
the tri-thorium core in the isolated crystalline cluster is non-aromatic. Similar conclusion was made by Foroutan-Nejad and co-workers who investigated the magnetic-response properties of the model tri-thorium cluster at the fully relativistic level of the theory.6 Interestingly, it has been demonstrated that the oxidation of the original cluster to its counterpart without the characteristic 3c2e-type HOMO noticeably enhances the charge-shift bonding in the Th3Cl6 cage (shorter Th–Cl bonds), which almost doubles the HOMO–LUMO gap.7 This means that the cyclic delocalization of electrons in the Th3 ring may affect thermodynamic stability of the entire complex thus putting into question the core–shell syngenetic model proposed by Lin and Mo.8 It should be noted that the proposed methodology to assess the ASE directly from the MO energy levels is limited to the three-membered rings investigated at the one-determinant (HF or DFT) theory level provided that the corresponding 3c2e-type MO does not effectively interfere with other orbitals in the MO spectrum (which is the case in the neutral Th3). Moreover, the predicted vanishing ASE may still be subject to some uncertainty since the proposed simplified model clusters (to a limited extent) suffer from the charge separation issue (the Th3 core in the isolated cluster is not neutral).3,10 Nevertheless, the presented results clearly show that aromatic stabilization in Th3 increases linearly with the total charge of the system, and therefore the results of investigations involving highly ionized model clusters cannot be used to justify σ-aromatic character of the isolated crystalline cluster.
Methodology
The extra cyclic resonance energy calculations were performed using the same methods and software as utilized in the original investigation by Lin and Mo.8 The input files (as well as details on implementation of the BLW-DFT code) were kindly provided by Dr Lin., and they are available from the author on request. All the molecular geometries of the model tri-thorium clusters were taken from the original article by Liddle and co-workers3 and the recently published paper by Lin and Mo.8 The relativistic calculations were performed with Gaussian G16 (rev. C01)16 and the implemented therein DKH2-relativistic Hamiltonian;17 for the neutral Th3 model cluster the spin–orbit coupling (SOC) effects were included in the calculations at different RTh–Th via the keyword DKHSO, but the results shown rather marginal effect on the calculated energy gaps (Fig. 5). The exchange-correlation functional ωB97X (ref. 18) by the Head-Gordon group was used to get reliable description of the one-electron energy levels, especially the frontiers ones, which are of crucial importance for the presented study. In this context it should be noticed that ωB97X has demonstrated that performs excellent in reproduction of the experimental ionization potentials for the aromatic systems of different types within the framework of the DFT-Koopman's theorem.19 In all relativistic calculations the all-electron triple-zeta valence basis set developed by Jorge and co-workers was used.20 It should be noted that other all-electron basis sets available for Th were found to suffer from significant Rydberg contributions to the ground-state density matrix by inclusion of orbitals with larger azimuthal quantum numbers (8s, 6f, etc.), which is unphysical and leads to exaggeration of bonding density between thorium atoms; a more comprehensive examination of the role of the choice of the basis set and the exchange-correlation functional on the description of bonding in actinide clusters will be the subject of a separate study. The SuperFine grid and very tight convergence criteria were used in the self-consistent field (SCF) calculations. To calculate the covalent bond orders and estimate the degree of charge localization in all the model Th3 clusters the linear coefficients of natural atomic orbitals were used from the NBO7 program (the keyword NAOMO in the $NBO section).15 The original code of the BLW-DFT method has been modified to facilitate direct calculations of the PBC effects; the modified BLW-DFT code is available on request from the author.
Conflicts of interest
There are no conflicts to declare.
Acknowledgements
The author expresses his special thanks to Prof. Xuhui Lin for his assistance with performing all the BLW-DFT calculations. The research was supported by the National Science Centre, Poland (2021/42/E/ST4/00332). The computer facilities and support were provided by the Polish high-performance computing infrastructure PLGrid (HPC Centers: ACK Cyfronet AGH) within computational grant no. PLG/2022/015950.
References
- A. Robinson, Nature, 2010, 465, 36 CrossRef CAS.
- J. M. Mercero, A. I. Boldyrev, G. Merino and J. M. Ugalde, Chem. Soc. Rev., 2015, 44, 6519 RSC.
- J. T. Boronski, J. A. Seed, D. Hunger, A. W. Woodward, J. van Slageren, A. J. Wooles, L. S. Natrajan, N. Kaltsoyannis and S. T. Liddle, Nature, 2021, 598, 72 CrossRef CAS PubMed; J. T. Boronski, J. A. Seed, D. Hunger, A. W. Woodward, J. van Slageren, A. J. Wooles, L. S. Natrajan, N. Kaltsoyannis and S. T. Liddle, Nature, 2022, 603, E21 CrossRef PubMed.
- G. Cavigliasso and N. Kaltsoyannis, Inorg. Chem., 2006, 45, 6828 CrossRef CAS PubMed; B. O. Roos, P.-Å. Malmqvist and L. Gagliardi, J. Am. Chem. Soc., 2006, 128, 17000 CrossRef PubMed; T. Cheisson, K. D. Kersey, N. Mahieu, A. McSkimming, M. R. Gau, P. J. Carroll and E. J. Schelter, J. Am. Chem. Soc., 2019, 141, 9185 CrossRef PubMed.
- R. Gershoni-Poranne and A. Stanger, Chem. Soc. Rev., 2015, 44, 6597 RSC.
- A. B. J. Cuyacot and C. Foroutan-Nejad, Nature, 2022, 603, E18 CrossRef PubMed.
- D. W. Szczepanik, Angew. Chem., Int. Ed., 2022, 61, e202204337 CrossRef CAS PubMed.
- X. Lin and Y. Mo, Angew. Chem., Int. Ed., 2022, 61, e202209658 CrossRef CAS PubMed.
- M. Solà and F. M. Bickelhaupt, J. Chem. Educ., 2022, 99, 3497 CrossRef PubMed.
- J. Tomeček, S. T. Liddle and N. Kaltsoyannis, ChemPhysChem, 2023, 24, e202300366 CrossRef PubMed.
- V. I. Minkin, Pure Appl. Chem., 1999, 71, 1919 CrossRef CAS.
- Y. Mo, L. Song and Y. Lin, J. Phys. Chem. A, 2007, 111, 8291 CrossRef CAS PubMed.
- K. Najmidin, A. Kerim, P. Abdirishit, H. Kalam and T. Tawar, J. Mol. Model., 2013, 19, 3529 CrossRef CAS PubMed.
- D. W. Szczepanik and J. Mrozek, J. Theor. Comput. Chem., 2011, 10, 471 CrossRef CAS; D. W. Szczepanik and J. Mrozek, J. Math. Chem., 2013, 51, 1619 CrossRef.
- E. D. Glendening, C. R. Landis and F. Weinhold, J. Comput. Chem., 2019, 40, 2234 CrossRef CAS PubMed.
- M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery Jr., J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman, and D. J. Fox, Gaussian 16, Revision C.01, Gaussian, Inc., Wallingford CT, 2016 Search PubMed.
- M. Douglas and N. M. Kroll, Ann. Phys., 1974, 82, 89 Search PubMed; B. A. Hess, Phys. Rev. A: At., Mol., Opt. Phys., 1985, 32, 756 CrossRef CAS PubMed; B. A. Hess, Phys. Rev. A: At., Mol., Opt. Phys., 1986, 33, 3742 CrossRef PubMed.
- J. D. Chai and M. Head-Gordon, Phys. Chem. Chem. Phys., 2008, 10, 6615 RSC.
- D. W. Szczepanik, M. Solà, M. Andrzejak, B. Pawełek, J. Dominikowska, M. Kukułka, K. Dyduch, T. M. Krygowski and H. Szatylowicz, J. Comput. Chem., 2017, 38, 1640 CrossRef CAS PubMed.
- A. Z. de Oliveira, C. T. Campos, F. E. Jorge, I. B. Ferreira and P. A. Fantin, Comput. Theor. Chem., 2018, 1135, 28 CrossRef CAS.
|
This journal is © The Royal Society of Chemistry 2023 |