Next Article in Journal
Multi-Component Syntheses of Spiro[furan-2,3′-indoline]-3-carboxylate Derivatives Using Ionic Liquid Catalysts
Previous Article in Journal
Stereoselective [4+3]-Cycloaddition of 2-Amino-β-nitrostyrenes with Azaoxyallyl Cations to Access Functionalized 1,4-Benzodiazepin-3-ones
Previous Article in Special Issue
Nanocomposites Based on Spin-Crossover Nanoparticles and Silica-Coated Gold Nanorods: A Nonlinear Optical Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 29
Issue 6
10.3390/molecules29061222
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Insights into ThB40: Stability, Electronic Structure, and Interaction

by
Yutian Li
1,*,
Yingying Wang
1,
Zhanrong Zhou
1,
Yang Gao
1,
Yiming Chen
1,
Guoqing Zhang
1 and
Chao Ma
2
1
Xi’an Research Institute of High Technology, Xi’an 710025, China
2
Foundation Department, Engineering University of PAP, Xi’an 710086, China
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(6), 1222; https://doi.org/10.3390/molecules29061222
Submission received: 4 February 2024 / Revised: 1 March 2024 / Accepted: 5 March 2024 / Published: 8 March 2024
(This article belongs to the Special Issue Inorganic Young Investigators' Contributions Collection)
Browse Figures
Versions Notes

Abstract

:
The interaction between nonmetal and metal atoms has attracted great interest in the development of organometallic compounds and their promising applications. In this study, we explored the interaction between boron and thorium atoms, based on the stable B40Th coordination compound, by employing density functional theory calculations. We elucidated the stability and geometries of the B40Th coordination compound and revealed the electron transfer from the metal atom Th to B40, which is evidenced by the natural bond orbital calculations. This electron transfer is attributed to the electron-withdrawing character of the boron atom and results in clear electrostatic interaction. Additionally, bond critical analysis and bond order calculations show obvious covalent characters between the metal and nonmetal atoms. The IR spectrum was simulated to give detailed information to identify this targeted compound in future experiments. This study is expected to enhance the understanding of metal–nonmetal interactions and provides useful information for constructing new organometallic compounds based on actinium metal atoms.

1. Introduction

The interaction between metal and nonmetal atoms, particularly actinium metal atoms, is complicated due to the complex electronic structure of metal atoms [1,2,3]. Understanding these metal–nonmetal interactions is crucial for the development of organometallic compounds [4,5,6], which have a wide range of promising applications in energy, environment, and information fields. This has attracted much attention and prompted further research into metal–nonmetal interactions to prepare new organometallic structures.
In the last few decades, endohedral metallofullerenes have been regarded as the ideal models to study the interaction between encaged metal atoms and nonmetal atoms, including actinium metal atoms [7,8,9,10,11,12,13]. The discovery of the first lanthanum-based endohedral metallofullerene in 1985 with nuclear-shell geometries provided a good environment to further study the interaction between metal and carbon atoms, because the spherical electronic structure of fullerene is resistant to external interaction [14]. Since then, various interatomic interactions have been studied, including metal–metal and metal–nonmetal interactions, in which the types of atoms cover a wide range of atoms, including carbon, oxygen, nitrogen, main group metals, transition metals, lanthanide metal (La-Lu), and actinium metal (Ac-Lr) [15,16,17,18,19]. Particularly, the actinium–carbon interaction is focused on according to the successful isolation of actinium-based metallofullerenes, including U@C2n (76 ≤ 2n ≤ 86) [20,21], UN@C82 [22], Th@C2n (76 ≤ 2n ≤ 86) [23,24,25], ThC2@C82 [25], U2@C80 [26], and Th2@C80 [27]. However, there has been limited research on the boron–actinium interaction.
Boron chemistry has rapidly developed in the last few decades, including research on elemental boron and boron-based organometallic structure [28,29,30,31,32]. Thus, detailed studies on metal–boron interactions are desired. It is known to all that boron atoms possess similar electronic structures to neighboring carbon atoms, and boron-based nanomaterials have also been developed with similar structures to carbon nanomaterials [33,34,35,36]. A previous study was successfully conducted on the borophene [37], which is a 2D-graphene like structure. An all-boron cage-shaped fullerene nanostructure was also recognized, with similar geometry to fullerenes, in 2014 by Wang and Li [38,39]. This study provided an opportunity to study the nonmetal–boron and metal–boron interactions based on metalloborospherenes. For example, Li et al. theoretically studied the Be-B and Zn-B interactions based on their core-shell octahedral structure [40], and the interaction between boron atoms and transition metal atoms (Sc, Ti, and V atoms) based on transition metal-centered endohedral seashell-like metalloborospherenes [41]. Clearly, previous reports have successfully conducted theoretical studies on the interaction between boron and main group atoms, as well as transition metal atoms [41]. However, there is scarce research on actinide–boron interactions to date.
In this study, we aimed to investigate the interaction between thorium and boron atoms based on metalloborospherenes by using density functional theory calculations (DFT). Electronic structures are studied by using natural bond orbital (NBO) calculations and frontier molecular orbital analysis. Additionally, metal–nonmetal interactions are studied through bond critical point (BCP) analysis and bond order calculations [42]. The corresponding IR spectrum is simulated in theory to provide useful information for geometrical recognition in future experiments. This research will contribute to a deeper understanding of metal–nonmetal interactions and may lead to the development of new organometallic compounds.

2. Results and Discussion

2.1. Stability of ThB40

The B40 isomers, D2d-B40 and Cs-B40, identified as potential candidates for ThB40 compounds following screening with the AM1 semiempirical molecular orbital method [43], were further optimized on the TPSSh/6-31G* theoretical level [44,45]. As shown in Figure 1, the relative energy of D2d-B40 was found to be 15.2 kcal/mol lower than that of Cs-B40, indicating the enhanced stability of the spherical-like D2d-B40 compared to Cs-B40. The lowest frequencies recorded were 168.6 cm−1 and 49.0 cm-1 for D2d-B40 and Cs-B40, respectively, indicating the structural stability of both isomers on the potential energy surface, with the potential well of the spherical borospherene D2d-B40 being notably deeper than that of the bowl-like Cs-B40. The lowest frequency of D2d-B40 and Cs-B40 also indicates the greater stability of the former isomer. Previous reports on the stability of these two B40 isomers align with the outcomes of the present theoretical calculations [39,46]. Furthermore, the electronic spin states of the two stable B40 isomers have been considered, confirming the singlet ground state (Table S1), as defined by the energy difference between singlet and triplet states. The energy of D2d-B40 and Cs-B40 in singlet state is 40.8 kcal/mol and 7.8 kcal/mol lower than that in triplet state, respectively, indicating the singlet-ground electron spin state. The lower energy difference for Cs-B40 is attributed to the likely plane geometry, ensuring the share of electrons.
The selected candidate isomers were utilized to form ThB40 compounds, and a total of twelve ThB40 isomers were considered, as shown in Figure S1. The relative energy for all of these ThB40 isomers is shown in Table 1. The findings show that encaged Th into D2d-B40 results in the same geometries with the lowest relative energy. Notably, the bowl-like Th@Cs-B40 possesses a higher relative energy of 64.2 kcal/mol, indicating the thermodynamic stability of Th-based metalloborospherenes. In order to further confirm the theoretical level, we conducted the optimization of Th@Cs-B40 and Th@D2d-B40, as shown in Table S2, and all of the results indicate that Th@D2d-B40 is more stable than Th@Cs-B40 based on the relative energy, which is in line with the results on the TPSSh/6-31G*~SDD theoretical level. The stable isomer, with the lowest relative energy, possesses an energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of 0.80 eV, indicating its kinetic stability in a room-temperature environment.
In order to further confirm the stability of ThB40 isomers, the Boltzmann distribution of three2-Th@D2d-B40 and hexa1-Th@Cs-B40 was carried out at 1000 K. The results indicate that the concentration of three2-Th@D2d-B40 exceeds 90%. All of these results strongly support the suitability of D2d-B40 as a host for encapsulating the Th atom within its hollow cavity. The following discussion is solely based on the Th-based endohedral metallic borospherene denoted as Th@D2d-B40, which serves as an ideal model to further study the interaction between Th and B atoms.

2.2. Geometries and Electronic Structures

The nature of molecules is determined by their geometry, and the geometry of stable Th@D2d-B40 is shown in Figure 2. The Th atom is situated just off the inner center of the borospherene D2d-B40 and lies within a plane of Cs symmetry. The distance between Th and B atoms ranges from 2.94 Å to 3.50 Å, which is larger than the Th-C distances ranging from 2.8 Å to 2.5 Å in Th@C76, Th@C80, Th@C82, and Th@C86 [23,27,47,48]. This phenomenon is attributed to the strong coordination of the Th atom with either a pentagon or hexagon in Th-based endohedral metallofullerenes, despite the larger carbon cage. This suggests the weaker interaction of Th-B in endohedral metalloborospherenes compared to Th-C in endohedral metallofullerenes, a point that will be further explored in the following discussions.
The singlet and triplet Th@D2d-B40 structures were optimized to determine their spin ground state. The energy of the triplet Th@D2d-B40 is 1.5 kcal/mol higher than that of the singlet structure on the TPSSh/6-31G*~SDD level. To further confirm the spin ground state, the energy difference between singlet and triplet states was calculated on the B3LYP, BP86, PBE0, and HSE06 functionals with 6-31G*~SDD, and all of the results in Table S3 indicate the lower relative energy of the singlet state. Under the consideration of the valence electronic structure of the Th atom (6d27s2), we inferred a four-electron transfer from Th atom to D2d-B40. NBO calculations (shown in Table 2) were carried out to explore the electron transfer and clarify the electronic structures of Th@D2d-B40. It is evident that the electrons of 6d and 7s orbitals of the Th atom have been lost, and a four-electron transfer is accepted by the borospherene D2d-B40. Consequently, the electronic structure of the endohedral metallic borospherene can be described as Th4+@[D2d-B40]4−.
As shown in Figure 3, the frontier molecular orbitals have been mapped to further confirm the anticipated four-electron transfer from the Th atom to the B40 cage. Notably, two orbitals, LUMO + 1 and LUMO + 2, exhibit degeneracy. The LUMO + 1 and LUMO of the optimized D2d-B40 transition to become the HOMO and HOMO-1 of the optimized Th@D2d-B40, also providing further evidence of the anticipated four-electron transfer from the Th atom to the borospherene D2d-B40 cage. Additionally, orbital overlap is observed in the HOMO and HOMO-1, indicating electron backdonation from the outer cage to the inner metal atom and suggesting a covalent interaction between Th and B atoms. This suggests that electron transition and chemical functionalization will occur in the outer borospherene cage.

2.3. Interactions between Th and B40

Due to the four-electron transfer, there is a clear electrostatic interaction between Th4+ and B404−. On the other hand, according to the NBO analysis and the mapping of frontier molecular orbitals, there is a clear electron backdonation from borospherene D2d-B40 to the Th atom, suggesting the presence of covalent interaction characteristics. Additionally, Mayer bond orders have been calculated to characterize the covalent interactions between them, as shown in Table 3 [49]. The Mayer bond order ranges from about 0.27 to about 0.34 between the Th atom and B atoms in the outer borospherene D2d-B40 cage, with the values of Mayer bond order positively correlating to the distance between Th and B atoms. These calculations indicate the presence of covalent characteristics between the inner metal atom Th and the outer borospherene D2d-B40 cage.
To further study the interaction between Th and B40, BCP analysis (Figure 4) was carried out [50,51], focusing specifically on the BCPs identified between the Th and B40 cages. The detailed parameters of these BCPs are shown in Table S4. In Figure 4, there is the presence of BCPs between Th and the outer borospherene D2d-B40 cage, with the average parameters of these BCPs calculated. Generally, the closer the |VBCP|/ρBCP ratio is to 1, the greater the covalent characteristics present, with VBCP and ρBCP denoting the potential energy density and electronic density of the BCP, respectively. The ρBCP value is larger than 0. Additionally, the average value of HBCP/ρBCP is negative, in which HBCP is the energy density of the BCP. Accordingly, there is a combination of electrostatic interactions and covalent interactions. In detail, the electrostatic interaction arises from the electron transfer between the inner metal atoms and the outer borospherene cage, while the covalent interaction is attributed to the backdonation from the outer borospherene cage to the inner Th atom. This suggests that the intra-interaction in Th-based endohedral metallic borospherene resembles the ionic model described in endohedral metallofullerens. Clearly, in Th-based endohedral metallic borospherene compounds, the bonding form differs from traditional metal-based coordinated compounds, in which the coordination interaction stems from the electron transfer from organic ligands to the central metal atom, and the electron backdonation occurs from the central metal atom to the organic ligand.

2.4. Simulated IR Spectra of D2d-B40 and Th@D2d-B40

Furthermore, the IR spectrum was simulated to provide more information to be able to identify Th@B40 in future experiments. The theoretically simulated IR spectra of D2d-B40 and Th@D2d-B40, as shown in Figure 5 at the TPSSh/6-31G* theoretical level, reveal several significant absorption peaks. Clearly, there are two absorption peaks at round 66 and 79 cm−1 for Th@D2d-B40 compared with D2d-B40. These two absorption peaks correspond to the vibrations which involve the metal atom Th. The maps of displacement vectors, as shown in Figure S3, provide further insight into these vibrations. Around 200 cm−1, two absorption peaks are observed, corresponding to the distortion vibration (at about ~187 cm−1) and telescopic vibration (at about ~202 cm−1). Additionally, clear B-B bonding stretch vibrations are evident at about 1200~1400 cm−1, while the absorption peaks in the range of about 400~800 cm−1 indicate shearing vibrations of the outer cages. Clearly, the presence of the Th atom not only leads to new featured absorption peaks, but also plays important roles in the vibrations of D2d-B40 because of the intra-molecular interaction. These featured absorption peaks in the simulated IR spectrum hold valuable implications for future experiments involving Th-based endohedral metallic borospherenes, offering crucial insights into their structural and vibration properties.

3. Calculation Methods

All of the density functional theory calculations in this work were performed by using Gaussian 16 software [52]. The semiempirical molecular orbital (AM1) method was used to preliminarily screen the stability of borospherene B40, and several isomers were selected with relative low energy. These selected isomers were further optimized on the TPSSh/6-31G* theoretical level [44,45,53]. The frequency calculations were also carried out on the same theoretical level, free from imaginary frequency for the optimized geometries, indicating the local energy minimum point on the potential energy surface for each optimized isomer. Then, geometry optimizations were carried out for the ThB40 isomers on the theoretical levels of TPSSh/3-21G~SDD and TPSSh/6-31G*~SDD [45,54]. Frequency calculations of optimized ThB40 were conducted on the TPSSh/6-31G*~SDD theoretical level, also free from imaginary frequency, indicating the potential stability of ThB40. In order to confirm the spin ground state of the most stable Th@B40, we considered the single and triplet states at the TPSSh/6-31G*~SDD theoretical level. The Boltzmann distribution was carried out at 1000 K according to the equation e − Δ E i / R T ∑ j e − Δ E j / R T , in which ΔE is the relative energy and T is the absolute temperature. Frontier molecular orbital and NBO calculations were carried out to analyze the electronic structures [55], and BCP analysis and Mayer bond order (MBO) calculations were used to analyze the bonding characteristics and the interaction characteristics between actinide Th and B40 of the most stable Th@B40 [49]. The BCP analysis and Mayer bond order calculations are based on the Multiwfn 3.8 software [50].

4. Conclusions

In conclusion, density functional theory results have revealed that borospherene can serve as a stable carrier for the actinide Th, as demonstrated by Th@B40, in which an ionic interaction is uncovered in combination with covalent characteristics between Th and B atoms. The most stable isomer, Th@D2d-B40, with a singlet ground state, has been identified. The HOMO-LUMO gap (0.804 eV), defined as the energy difference between the energy of the HOMO and LUMO, and the energy levels of HOMO (−4.18 eV) and LUMO (−3.37 eV) suggest that redox processes can easily occur on the outer borospherene cage. Natural bond orbital analysis has revealed that there is a formal four-electron transfer from the inner Th atom to the outer borospherene B40 cage, with backdonation occurring from the outer borospherene cage to the inner metal atom Th. Furthermore, according to the frontier molecular orbital analysis of the hollow borospherene cage and Th-based endohedral metallic borospherene, the electronic configuration has been confirmed as Th4+@[D2d-B40]4−. The NBO and frontier molecular orbital analysis have indicated a significant electrostatic interaction between the inner metal Th atom and the outer borospherene cage. Additionally, the characters of covalent interaction have also been confirmed by the calculation of the MBO between the Th atom and B atom, with the MBO being determined by their distances. The BCP analysis also suggests ionic interaction between the inner metal atom Th and the outer borospherene cage, in combination with covalent characteristics. The simulated IR spectrum has provided valuable information for geometrical identification in future experiments of Th-based endohedral metallic borospherenes. This is expected to enhance our understanding of the interaction between nonmetal atoms and metal atoms, particularly regarding the interaction between B atom and f-metal atoms, and to provide a certain type of guiding information for the research and application of borospherenes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29061222/s1, Figure S1: Geometries of potential isomers of ThB40 based on D2d-B40 and Cs-B40 isomers, Table S1: relative energy of D2d-B40 and Cs-B40 in different spin multiplicities, Tables S2 and S3: relative energy of Th@D2d-B40 and Th@Cs-B40 at different functionals, and relative energy of Th@D2d-B40 in different spin multiplicities. Figures S2 and S3, Table S4: Geometry of Th@D2d-B40 with atom labels, parameters of bond critical points, featured displacement vectors of Th@D2d-B40, and coordinates of optimized geometries. Coordinates of D2d-B40, Cs-B40, Th@D2d-B40, and ThCs-B40.

Author Contributions

Y.L. conducted all of the DFT computations, Y.W., Z.Z., Y.C., Y.G. and C.M. provided valuable discussion, Y.L. wrote the manuscript, Y.L. and G.Z. revised the manuscript, and Y.L. conceived the project idea and supervised the project. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge support from the Natural Science Foundation of Shaanxi Province (Grant No. 2023-JC-QN-0082, 2023-JC-QN-0034, and 2022JM-341), the Youth Fund of PLA Rocket Force University of Engineering (Grant No. 2021QN-B030), and the National Natural Science Foundation of China (Grant No. 12304317).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Conflicts of Interest

The authors declare no competing financial interests.

References

  1. Luo, M.; Xu, Y.E.; Song, Y.X. Ab initio study on nonmetal and nonmagnetic metal atoms doped arsenene. JETP Lett. 2017, 106, 434–439. [Google Scholar] [CrossRef]
  2. Li, M.Y.; Cui, J.B.; Zhao, Y.X.; Zhao, P.; Li, Q.Z.; Zhao, X. Unexpected diverseness on electronic density and bonding behaviours for Sc2X@C2v(63751)-C86 and Sc2X@C1(63755)-C86 (X = S and 0). Chem. Phys. Lett. 2018, 707, 93–100. [Google Scholar] [CrossRef]
  3. Ifthikar, J.; Shahib, I.I.; Jawad, A.; Gendy, E.A.; Wang, S.Q.; Wu, B.B.; Chen, Z.Q.; Chen, Z.L. The excursion covered for the elimination of chromate by exploring the coordination mechanisms between chromium species and various functional groups. Coord. Chem. Rev. 2021, 437, 213868. [Google Scholar] [CrossRef]
  4. Dimic, D.; Eichhorn, T.; Milenkovic, D.; Kaluderovic, G.N. Synthesis, Structural, and Quantum Chemical Analysis of Neutral and Cationic Ruthenium(II) Complexes with Nicotinate-Polyethylene Glycol Ester Ligands. Inorganics 2023, 11, 460. [Google Scholar] [CrossRef]
  5. Mato, M.; Cornella, J. Bismuth in Radical Chemistry and Catalysis. Angew. Chem. Int. Ed. 2023, 63, e202315046. [Google Scholar] [CrossRef] [PubMed]
  6. Weinhold, F.; Glendening, E.D. Natural resonance-theoretic conceptions of extreme electronic delocalization in soft materials. Phys. Chem. Chem. Phys. 2024, 26, 2815–2822. [Google Scholar] [CrossRef] [PubMed]
  7. Hirsch, A.; Chen, Z.; Jiao, H. Spherical Aromaticity in Ih Symmetrical Fullerenes The 2(N+1)2 Rule. Angew. Chem. Int. Ed. 2000, 39, 3915–3917. [Google Scholar] [CrossRef]
  8. Popov, A.A.; Dunsch, L. Bonding in endohedral metallofullerenes as studied by quantum theory of atoms in molecules. Chemistry 2009, 15, 9707–9729. [Google Scholar] [CrossRef]
  9. Canchaya, J.G.; Wang, Y.; Alcami, M.; Martin, F.; Busnengo, H.F. Study of the interaction between short alkanethiols from ab initio calculations. Phys. Chem. Chem. Phys. 2010, 12, 7555–7565. [Google Scholar] [CrossRef]
  10. Rodriguez-Fortea, A.; Balch, A.L.; Poblet, J.M. Endohedral metallofullerenes: A unique host-guest association. Chem. Soc. Rev. 2011, 40, 3551–3563. [Google Scholar] [CrossRef]
  11. Yang, S.; Liu, F.; Chen, C.; Jiao, M.; Wei, T. Fullerenes encaging metal clusters—Clusterfullerenes. Chem. Commun. 2011, 47, 11822–11839. [Google Scholar] [CrossRef]
  12. Lu, X.; Feng, L.; Akasaka, T.; Nagase, S. Current status and future developments of endohedral metallofullerenes. Chem. Soc. Rev. 2012, 41, 7723–7760. [Google Scholar] [CrossRef]
  13. Bao, L.; Peng, P.; Lu, X. Bonding inside and outside Fullerene Cages. Acc. Chem. Res. 2018, 51, 810–815. [Google Scholar] [CrossRef] [PubMed]
  14. Heath, J.; O’Brien, S.; Zhang, Q.; Liu, Y.; Curl, R.; Tittel, F.; Smalley, R. Lanthanum complexes of spheroidal carbon shells. J. Am. Chem. Soc. 1985, 107, 7779–7780. [Google Scholar] [CrossRef]
  15. Alvarez, M.M.; Gillan, E.G.; Holczer, K.; Kaner, R.B.; Min, K.S.; Whetten, R.L. Lanthanum carbide (La2C80): A soluble dimetallofullerene. J. Phys. Chem. 1991, 95, 10561–10563. [Google Scholar] [CrossRef]
  16. Chai, Y.; Guo, T.; Jin, C.; Haufler, R.E.; Chibante, L.F.; Fure, J.; Wang, L.; Alford, J.M.; Smalley, R.E. Fullerenes with metals inside. J. Phys. Chem. 1991, 95, 7564–7568. [Google Scholar] [CrossRef]
  17. Moreno-Vicente, A.; Roselló, Y.; Chen, N.; Echegoyen, L.; Dunk, P.W.; Rodríguez-Fortea, A.; de Graaf, C.; Poblet, J.M. Are U–U Bonds Inside Fullerenes Really Unwilling Bonds? J. Am. Chem. Soc. 2023, 145, 6710–6718. [Google Scholar] [CrossRef] [PubMed]
  18. Li, M.Y.; Luo, X.L.; Zhao, Y.X.; Zhang, W.X.; Yuan, K.; Zhao, X. Metal Atoms (Li, Na, and K) Tuning the Configuration of Pyrrole for the Selective Recognition of C60. Inorg. Chem. 2023, 62, 4618–4624. [Google Scholar] [CrossRef]
  19. Jiang, H.; Yu, X.; Guo, M.; Yao, Y.-R.; Meng, Q.; Echegoyen, L.; Autschbach, J.; Chen, N. USc2C2 and USc2NC Clusters with U–C Triple Bond Character Stabilized Inside Fullerene Cages. J. Am. Chem. Soc. 2023, 145, 5645–5654. [Google Scholar] [CrossRef]
  20. Cai, W.; Abella, L.; Zhuang, J.; Zhang, X.; Feng, L.; Wang, Y.; Morales-Martinez, R.; Esper, R.; Boero, M.; Metta-Magana, A.; et al. Synthesis and Characterization of Non-Isolated-Pentagon-Rule Actinide Endohedral Metallofullerenes U@ C1(17418)-C76, U@C1(28324)-C80, and Th@C1(28324)-C80: Low-Symmetry Cage Selection Directed by a Tetravalent Ion. J. Am. Chem. Soc. 2018, 140, 18039–18050. [Google Scholar] [CrossRef]
  21. Zhao, Y.X.; Li, M.Y.; Zhao, P.; Ehara, M.; Zhao, X. New Insight into U@C80: Missing U@D3(31921)-C80 and Nuanced Enantiomers of U@C1(28324)-C80. Inorg. Chem. 2019, 58, 14159–14166. [Google Scholar] [CrossRef]
  22. Liu, F.; Spree, L. Molecular spinning top: Visualizing the dynamics of M3N@C80 with variable temperature single crystal X-ray diffraction. Chem. Commun. 2019, 55, 13000–13003. [Google Scholar] [CrossRef]
  23. Yan, Y.; Morales-Martinez, R.; Zhuang, J.; Yao, Y.R.; Li, X.; Poblet, J.M.; Rodriguez-Fortea, A.; Chen, N. Th@D5h(6)-C80: A highly symmetric fullerene cage stabilized by a single metal ion. Chem. Commun. 2021, 57, 6624–6627. [Google Scholar] [CrossRef]
  24. Jin, M.; Zhuang, J.; Wang, Y.; Yang, W.; Liu, X.; Chen, N. Th@Td(19151)-C76: A Highly Symmetric Fullerene Cage Stabilized by a Tetravalent Actinide Metal Ion. Inorg. Chem. 2019, 58, 16722–16726. [Google Scholar] [CrossRef]
  25. Shen, Y.; Yu, X.; Meng, Q.; Yao, Y.-R.; Autschbach, J.; Chen, N. ThC2@C82 versus Th@C84: Unexpected formation of triangular thorium carbide cluster inside fullerenes. Chem. Sci. 2023, 13, 12980–12986. [Google Scholar] [CrossRef]
  26. Zhang, X.; Wang, Y.; Morales-Martínez, R.; Zhong, J.; de Graaf, C.; Rodríguez-Fortea, A.; Poblet, J.M.; Echegoyen, L.; Feng, L.; Chen, N. U2@Ih(7)-C80: Crystallographic Characterization of a Long-Sought Dimetallic Actinide Endohedral Fullerene. J. Am. Chem. Soc. 2018, 140, 3907–3915. [Google Scholar] [CrossRef]
  27. Zhuang, J.; Morales-Martínez, R.; Zhang, J.; Wang, Y.; Yao, Y.R.; Pei, C.; Rodríguez-Fortea, A.; Wang, S.; Echegoyen, L.; de Graaf, C.; et al. Characterization of a strong covalent Th3+-Th3+ bond inside an Ih(7)-C80 fullerene cage. Nat. Commun. 2021, 12, 2372. [Google Scholar] [CrossRef]
  28. Zhao, Y.F.; Lusk, M.T.; Dillon, A.C.; Heben, M.J.; Zhang, S.B. Boron-based organometallic nanostructures: Hydrogen storage properties and structure stability. Nano Lett. 2008, 8, 157–161. [Google Scholar] [CrossRef]
  29. Ramachandrant, P.V.; Burghardt, T.E. Recent developments in the chiral synthesis of homoallylic amines via organoboranes. Pure Appl. Chem. 2006, 78, 1397–1406. [Google Scholar] [CrossRef]
  30. Bigmore, H.R.; Lawrence, S.C.; Mountford, P.; Tredget, C.S. Coordination, organometallic and related chemistry of tris(pyrazolyl) methane ligands. Dalton Transact. 2005, 1, 635–651. [Google Scholar] [CrossRef]
  31. Xue, K.; Li, H.Y.; Pan, L.; Li, C.J.; Zhang, X.W.; Zou, J.J. Preparation and performance characterization of functionalized boron-based energetic-microcapsules with uniform size. Chem. Eng. J. 2023, 469, 143917. [Google Scholar] [CrossRef]
  32. Ali, F.; Hosmane, N.S.; Zhu, Y.H. Boron Chemistry for Medical Applications. Molecules 2020, 25, 828. [Google Scholar] [CrossRef]
  33. Chen, W.H.; Liu, J.Y.; Sun, W.M.; He, H.M.; Yu, S.S.; Li, Y.; Li, Z.R. Metalloborospherenes with the stabilized classical fullerene-like borospherene B36 as electric field manipulated second-order nonlinear optical switches. New J. Chem. 2022, 46, 22246–22255. [Google Scholar] [CrossRef]
  34. Soliman, K.A.; Aal, S.A. Ti, Ni, and Cu decorated borospherene as potential molecular sensor for phosgene. Mat. Sci. Semicon. Proc. 2022, 144, 106574. [Google Scholar] [CrossRef]
  35. Pei, L.; Zhang, L.J.; Li, D.Z. Theoretical study on exohedral complexes C6H6TMB40 (TM = Sc-Ni). Phys. Chem. Chem. Phys. 2022, 24, 21794–21799. [Google Scholar] [CrossRef]
  36. Kaur, H.; Kaur, J.; Kumar, R. Comparative study of symmetrical and asymmetrical B40 molecular junctions. J. Comput. Electron 2022, 21, 599–607. [Google Scholar] [CrossRef]
  37. Chinnalagu, D.K.; Murugesan, B.; Arumugam, M.; Chinniah, K.; Ganesan, S.; Cai, Y.R.; Mahalingam, S. Fabrication of 2D-Borophene nanosheets anchored S, N-mesoporous carbon nanocomposite (SNC-Bp//SNC-Bp) symmetric device for high-performance supercapacitor application. J. Energy Storage 2023, 74, 109328. [Google Scholar] [CrossRef]
  38. Wang, L.S. Photoelectron spectroscopy of size-selected boron clusters: From planar structures to borophenes and borospherenes. Int. Rev. Phys. Chem. 2016, 35, 69–142. [Google Scholar] [CrossRef]
  39. Zhai, H.J.; Zhao, Y.F.; Li, W.L.; Chen, Q.; Bai, H.; Hu, H.S.; Piazza, Z.A.; Tian, W.J.; Lu, H.G.; Wu, Y.B.; et al. Observation of an all-boron fullerene. Nat. Chem. 2014, 6, 727–731. [Google Scholar] [CrossRef]
  40. Yan, Q.Q.; Zhao, X.; Zhang, T.; Li, S.D. Perfect Core-Shell Octahedral B@B38+, Be@B38, and Zn@B38 with an Octa-Coordinate Center as Superatoms Following the Octet Rule. ChemPhysChem 2023, 24, e202200947. [Google Scholar] [CrossRef]
  41. Zhang, T.; Zhang, M.; Lu, X.-Q.; Yan, Q.-Q.; Zhao, X.-N.; Li, S.-D. Sc@B28−, Ti@B28, V@B28+, and V@B292−: Spherically Aromatic Endohedral Seashell-like Metallo-Borospherenes. Molecules 2023, 28, 3892. [Google Scholar] [CrossRef]
  42. Müller, M.; Hansen, A.; Grimme, S. An atom-in-molecule adaptive polarized valence single-ζ atomic orbital basis for electronic structure calculations. J. Chem. Phys. 2023, 159, 164108. [Google Scholar] [CrossRef]
  43. Ishihara, M.; Hatano, H.; Kawase, M.; Sakagami, H. Estimation of Relationship Between the Structure of 1,2,3,4-Tetrahydroisoquinoline Derivatives Determined by a Semiempirical Molecular-Orbital Method and their Cytotoxicity. Anticancer Res. 2009, 29, 2265–2271. [Google Scholar]
  44. Conradie, M.M.; Conradie, J.; Ghosh, A. Capturing the spin state diversity of iron(III)-aryl porphyrins OLYP is better than TPSSh. J. Inorg. Biochem. 2011, 105, 84–91. [Google Scholar] [CrossRef]
  45. Jensen, K.P. Bioinorganic Chemistry Modeled with the TPSSh Density Functional. Inorg. Chem. 2008, 47, 10357–10365. [Google Scholar] [CrossRef]
  46. Fa, W.; Chen, S.; Pande, S.; Zeng, X.C. Stability of Metal-Encapsulating Boron Fullerene B40. J. Phys. Chem. A 2015, 119, 11208–11214. [Google Scholar] [CrossRef]
  47. Zhao, P.; Ehara, M. Theoretical Insights into Monometallofullerene Th@C76: Strong Covalent Interaction between Thorium and the Carbon Cage. Inorg. Chem. 2018, 57, 2961–2964. [Google Scholar] [CrossRef]
  48. Muñoz-Castro, A.; King, R.B. Th@C86, Th@C82, Th@C80, and Th@C76: Role of thorium encapsulation in determining spherical aromatic and bonding properties on medium-sized endohedral metallofullerenes. Phys. Chem. Chem. Phys. 2020, 22, 23920–23928. [Google Scholar] [CrossRef]
  49. Bridgeman, A.J.; Cavigliasso, G.; Ireland, L.R.; Rothery, J. The Mayer bond order as a tool in inorganic chemistry. J. Chem. Soc. Dalton 2001, 1, 2095–2108. [Google Scholar] [CrossRef]
  50. Lu, T.; Chen, F.W. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  51. Cao, H.; Han, X.Y.; Wu, H.Y. Efficient atom-molecule conversion in Bose-Einstein condensates based on nonlinear dressed-state scheme. Chaos Soliton Fract. 2023, 174, 113882. [Google Scholar] [CrossRef]
  52. Frisch, M.; Trucks, G.; Schlegel, H.; Scuseria, G.; Robb, M.; Cheeseman, J.; Scalmani, G.; Barone, V.; Petersson, G.; Nakatsuji, H. Gaussian 16, Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2019. [Google Scholar]
  53. Prascher, B.P.; Woon, D.E.; Peterson, K.A.; Dunning, T.H.; Wilson, A.K. Gaussian basis sets for use in correlated molecular calculations. VII. Valence, core-valence, and scalar relativistic basis sets for Li, Be, Na, and Mg. Theor. Chem. Acc. 2011, 128, 69–82. [Google Scholar] [CrossRef]
  54. Woon, D.E.; Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 1993, 98, 1358–1371. [Google Scholar] [CrossRef]
  55. Glendening, E.D.; Landis, C.R.; Weinhold, F. Natural bond orbital methods. WIR Comput. Mol. Sci. 2012, 2, 1–42. [Google Scholar] [CrossRef]
Molecules 29 01222 g001
Figure 1. Geometries of D2d-B40 and Cs-B40 optimized on the TPSSh/6-31G* theoretical level, including their relative energy in kcal/mol and lowest frequency.
Figure 1. Geometries of D2d-B40 and Cs-B40 optimized on the TPSSh/6-31G* theoretical level, including their relative energy in kcal/mol and lowest frequency.
Molecules 29 01222 g001
Molecules 29 01222 g002
Figure 2. Geometries of optimized Th@D2d-B40 on the TPSSh/6-31G*~SDD theoretical level, including top view and side view. The Th and B atoms are colored in blue and pink, respectively.
Figure 2. Geometries of optimized Th@D2d-B40 on the TPSSh/6-31G*~SDD theoretical level, including top view and side view. The Th and B atoms are colored in blue and pink, respectively.
Molecules 29 01222 g002
Molecules 29 01222 g003
Figure 3. Frontier molecular orbitals of optimized D2d-B40 and Th@D2d-B40, including the energy level in eV.
Figure 3. Frontier molecular orbitals of optimized D2d-B40 and Th@D2d-B40, including the energy level in eV.
Molecules 29 01222 g003
Molecules 29 01222 g004
Figure 4. Searching bond critical points of optimized Th@D2d-B40, including the average parameters for bond critical points between inner metal atom Th and outer borospherene D2d-B40 cage.
Figure 4. Searching bond critical points of optimized Th@D2d-B40, including the average parameters for bond critical points between inner metal atom Th and outer borospherene D2d-B40 cage.
Molecules 29 01222 g004
Molecules 29 01222 g005
Figure 5. Simulated IR spectra of (a) D2d-B40 and (b) Th@D2d-B40 at the TPSSh/6-31G* theoretical level. The broadening function is selected as Lorentzian function, and the width is set as 12 cm−1 at half maximum.
Figure 5. Simulated IR spectra of (a) D2d-B40 and (b) Th@D2d-B40 at the TPSSh/6-31G* theoretical level. The broadening function is selected as Lorentzian function, and the width is set as 12 cm−1 at half maximum.
Molecules 29 01222 g005
Table 1. Relative energies ΔE of ThB40 isomers optimized on the TPSSh/6-31G*~SDD theoretical level and their HOMO-LUMO gaps.
Table 1. Relative energies ΔE of ThB40 isomers optimized on the TPSSh/6-31G*~SDD theoretical level and their HOMO-LUMO gaps.
IsomersΔE (kcal/mol)Gap (eV)IsomersΔE (kcal/mol)Gap (eV)
three2-Th@D2d-B400.00.80hexa1-Th@Cs-B4064.20.86
hepta-Th@D2d-B400.00.80Th@Cs-B4066.61.0
three-Th@D2d-B400.00.80three-Th@Cs-B4080.90.87
two-Th@D2d-B400.00.80hexa2-Th@Cs-B4091.70.97
center-Th@D2d-B400.00.80---
hexa-Th@D2d-B400.00.80---
hexa2-Th@D2d-B4084.20.98---
hexa3-Th@D2d-B4091.61.12---
Table 2. Natural population analysis for optimized Th@D2d-B40, including Th and B atoms in which the targeted B atom is closest to Th atom, on the TPSSh/6-31G*~SDD theoretical level.
Table 2. Natural population analysis for optimized Th@D2d-B40, including Th and B atoms in which the targeted B atom is closest to Th atom, on the TPSSh/6-31G*~SDD theoretical level.
AtomsPopulations
Th7s0.015f0.376d0.147p0.528s0.19
B2s0.532p2.453s0.013p0.02
Table 3. Mayer bond order between Th atom and several B atoms close to the metal atom (Figure S2), and their bond length in Å in optimized Th@D2d-B40 on the TPSSh/6-31G*~SDD theoretical level.
Table 3. Mayer bond order between Th atom and several B atoms close to the metal atom (Figure S2), and their bond length in Å in optimized Th@D2d-B40 on the TPSSh/6-31G*~SDD theoretical level.
BondsMayer Bond OrderBond Length
Th-B70.3422.94
Th-B100.2693.26
Th-B170.3422.94
Th-B320.2902.97
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Y.; Wang, Y.; Zhou, Z.; Gao, Y.; Chen, Y.; Zhang, G.; Ma, C. Insights into ThB40: Stability, Electronic Structure, and Interaction. Molecules 2024, 29, 1222. https://doi.org/10.3390/molecules29061222

AMA Style

Li Y, Wang Y, Zhou Z, Gao Y, Chen Y, Zhang G, Ma C. Insights into ThB40: Stability, Electronic Structure, and Interaction. Molecules. 2024; 29(6):1222. https://doi.org/10.3390/molecules29061222

Chicago/Turabian Style

Li, Yutian, Yingying Wang, Zhanrong Zhou, Yang Gao, Yiming Chen, Guoqing Zhang, and Chao Ma. 2024. "Insights into ThB40: Stability, Electronic Structure, and Interaction" Molecules 29, no. 6: 1222. https://doi.org/10.3390/molecules29061222

APA Style

Li, Y., Wang, Y., Zhou, Z., Gao, Y., Chen, Y., Zhang, G., & Ma, C. (2024). Insights into ThB40: Stability, Electronic Structure, and Interaction. Molecules, 29(6), 1222. https://doi.org/10.3390/molecules29061222

Article Metrics

Back to TopTop