Endohedral metalloborofullerenes (EMBFs) are novel boron analogues of the famous endohedral metallofullerenes (EMFs). Many EMBFs have been proposed by theoretical calculations thus far.However, in sharp contrast to EMFs, which trap most of the lanthanides with f electrons inside the cages, the corresponding lanthanide-based EMBFs have never been reported. In this work, the encapsulation of Eu and Gd in the B 38 and B 40 fullerenes was studied by means of density functional theory calculations. Our results revealed that Gd@B 38 ( 9 A), Eu@B 40 ( 8 B 2 ), and Gd@B 40 ( 7 A 00 ) all favor the endohedral configuration, and the electronic structures can be described as Gd 31 @B 32 38 , Eu 21 @B 22 40 , and Gd 31 @B 32 40 with jailed f electron spins. The large binding energies and sizable HOMO-LUMO gaps suggest that they may be achieved experimentally. They feature r and p double aromaticity, and their excellent stabilities were confirmed by the Born-Oppenheimer molecular dynamics simulations. Finally, the infrared and UV/vis spectra were simulated to assist experimental characterization.borofullerenes, density functional theory calculations, endohedral fullerenes, lanthanide element, spectra simulation 1 | I N TR ODU C TI ON Since the discovery of C 60 buckyball in 1985, [1] fullerenes and their various derivatives have attracted great attention due to their beautiful structures and unique physical and chemical properties. In particular, with the presence of metal elements in the cavity, endohedral metallofullerenes (EMFs) exhibit distinctive and tunable properties that differ dramatically from those of their parent cages. [2][3][4][5][6] To date, most of the lanthanides and even some actinide elements have been successfully encapsulated in the fullerene cages, with much attention paid to those with f electrons. [2][3][4][5][6][7] For example, gadolinium holds a half-filled 4f subshell (4f 7 ) with strong paramagnetic property, and the carbon cages are stable enough to prevent Gd from leakage under physiological conditions. Numerous researches have demonstrated that the Gd-based EMFs (e.g., Gd@C 82 ), after polyhydroxylation, are highly effective and nontoxic contrasting agents for magnetic resonance imaging. [8][9][10] Moreover, the isolated f electron spins on the metal could serve as a nanosized qubit in quantum computers and be applied in quantum information processing. [11] Besides the carbon fullerenes, the boron clusters of certain sizes can form cage-like structures as well, as exemplified by the recent fascinating discovery of the first borofullerene, D 2d -B -=0 40 (Figure 1). [12] After B 40 , more borofullerenes B n in neutral or ionic forms were observed in experiments or predicted by theory. [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29] Inspired by the facile preparation of EMFs, various endohedral derivatives of small borofullerenes were intensively proposed by computations as well. These complexes, named as endohedral metalloborofullerenes (EMBFs) [30][31][32] or endohedral meta...