Defect induced electronic structure of uranofullerene.

The interaction between the inner atoms/cluster and the outer fullerene cage is the source of various novel properties of endohedral metallofullerenes. Herein, we introduce an adatom-type spin polarization defect on the surface of a typical endohedral stable U(2)@C(60) to predict the associated structure and electronic properties of U(2)@C(61) based on the density functional theory method. We found that defect induces obvious changes in the electronic structure of this metallofullerene. More interestingly, the ground state of U(2)@C(61) is nonet spin in contrast to the septet of U(2)@C(60). Electronic structure analysis shows that the inner U atoms and the C ad-atom on the surface of the cage contribute together to this spin state, which is brought about by a ferromagnetic coupling between the spin of the unpaired electrons of the U atoms and the C ad-atom. This discovery may provide a possible approach to adapt the electronic structure properties of endohedral metallofullerenes.

Endohedral metallofullerenes (EMFs) are novel fullerene derivatives that can encapsulate metal atom or even clusters in their inner space12. Previous studies show that the inner cluster and/or the outer cage of the unique structure preferentially adopt the stable closed-shell electronic configuration34, such as La2@C725, La2@C8067, Sc2@C668, Sc2C2@C689, Sc3N@Cn (n = 68, 78, 80)10111213, Ti2C2@C7814. There is a significant charge transfer between the metal atoms and the fullerene1516. In EMFs, the internal metal atoms always donate electrons to the fullerene cage and carry considerable positive charges. This can result in an enhanced chemical reactivity as compared to the hollow fullerene12. EMFs have strong ability of donate and accept electrons, implying that they can be easily oxidized and reduced17. The stability of the whole system can be altered by this oxidation or reduction. The physical properties of EMFs strongly depend on the properties of the inner metal atom. Such a series of fantastic M@Cn structures have drawn extensive research interests in recent years. It is well known that defects can be easily formed in a large variety of carbon nanostructures, such as fullerenes, nanotubes and graphene. This defect can also bring about a significant change in the electronic structure as well as many other properties181920. But the natures of defective EMFs system are relatively unexplored till date. This work is aimed at gaining an understanding in this area through the study of an EMF structure with a typical spin polarized defective fullerene, which will be potentially beneficial in the design or synthesis of novel materials in several areas including catalysis and medicine2122.

EMFs contain U, Th, Pu and other actinide elements are very important complexes. Actinide metallofullerene materials not only have the thermal stability of the fullerene-like, but also have the physicochemical properties of the actinide elements, and therefore it is foreseeable that they will have important potential application in advanced fuel technology. So, studies on EMFs containing U atom have been widely reported. Guo et al. had prepared uranium fullerenes U@C28 and U@C60 experimentally for the first time. C28 is usually unstable. The U atom wrapped into C28 stabilizes the fullerene cage23. Subsequently, Diener et al. synthesized U@C74 and U@C8224. Zhao and his co-workers successfully synthesized and isolated a new type of metallofullerene containing metal atoms of transuranic elements for the first time, proving that actinide metal atoms can stabilize the otherwise unstable C82 fullerene25, and explained the structural stability of U@C82 through theoretical calculation of density functional theory (DFT)26. The uranium dimer compound is a rare compound whose presence is evidenced only in experimental spectra27. Therefore, its theoretical study is scarce. Gagliardi and Roos have theoretically studied the U2 molecule using high-level ab initio method and found that this molecule contains three electron-pair bonds and four one-electron bonds, corresponding to a unique quintuple bond, and the U2 molecule has a septet ground state28. Theorists have also predicted the possibility of fullerene encapsulating two U atoms based on quantum first principles method. Usually, the exohedral dimetal fullerenes are more stable than the endohedral structures, e.g. Cr2, Mo2 and W229. However, the U2 case is different. Its endohedral complex is significantly more stable than the exohedral and U2 lies at the center of the C60 fullerene29. The U-U metal bond tends to dissociate and the U atoms interact with the inner wall of the cage when the neutral U2 molecule is encapsulated into carbon cage; but due to the special size of the C60 fullerene, the U-U bond is weakened but not completely broken in the U2@C60 EMF. The results from a DFT study using the GGA-PBE exchange correlation functional shows that the dimetalloendofullerene, U2@C60, has an unprecedented U-U multiple bond consisting solely of sixfold ferromagnetically coupled one-electron-two-center bonds, which are dominated by the uranium 5f atomic orbitals, and it is septet in ground state30. This work opened up a new chemical field of the interaction between the actinide metal-metal bond and the fullerene, and provided the basis for many other studies. The size of the cavity can affect the properties of EMFs. If we enlarge the size to, say, C70 and C84, the metal-metal bond will be completely dissociated but the whole systems still have a septet ground state29.

Carbon nanostructures occurring in nature mostly contain defects, which cause changes in the electronic structure of carbon based nanomaterials and have been widely reported. However, fullerene is a necessary ingredient in EMFs. If the cage surface of the EMF which encapsulates U2 molecule within it contains a defect, the nature of such a structure is expected to be much different. So, we have considered an adatom defect by introducing a carbon ad-atom over the [6, 6] bond on the typical U2@C60 system and analyzed the features of this defective EMF U2@C61.

Results

In order to examine the reliability of the exchange correlation (XC) functional and basis sets, we calculated the isolated U2 molecule, the C60 fullerene and the EMF, U2@C60. All the results qualitatively agree with the previous reports. Furthermore, the results showed that defective C61 has a triplet ground state. The adatom contributes significantly to the spin. The details for the calculations have been presented in the supplementary information Part 1.3.

In C61 fullerene, we have considered six isomers with the two uranium atoms at different positions. The full geometry optimization results showed that the total energy of the whole system is progressively reduced as a uranium atom gradually moves closer to the defect. Eventually, we found the structure c in all isomers to be the most stable. Structure c can be characterized as one U atom located below the defect and another U atom facing a six-ring as shown in Figure 1. The details for searching the lowest energy structure are provided in the supplementary information Part 2.

Figure 1

The equilibrium structure (Structure c) with the lowest energy of all the U2@C61 isomers.

To investigate the possible mechanism responsible for the structural stability of c, we have proceeded systematically. First of all, a qualitative explanation is provided as follows. The single U atom is reported to tend to adsorb on a hexagonal face from the inside when it is encapsulated in a carbon cage26. The neutral U2 molecule is found to show a tendency to dissociate and also adsorb on the hexagonal faces in C60 fullerene2930. Interestingly, the adatom-type defect can enlarge the local size of the C60. We have calculated and compared the size of the C61 and C60. The diameter of C61 is about 0.31 Å longer than the original C60 (Ih) along the direction of the defect. Moreover, the defect dissociates the original [6, 6] bond of the C60. Consequently, the distance between these two C atoms formerly forming the [6, 6] bond in the original C60 increases by 0.79 Å. This causes the inner U atom to sense the outer unsaturated ad-atom. Based on the discussion above, there is no doubt that both the size effect and the impact of the adatom directionally guide the motion of one of the two U atoms to the defect area. As a result, the whole system is energetically favored to form the unique structure c as depicted in Figure 1. The distance of the two U atoms in C61 fullerene (2.80 Å of BP86, 2.88 Å of PBE0, see Table 1) is larger than the U2@C60 case which has been calculated at the same level (2.70 Å of BP86, see supplementary information Part 1.3), including all the U2@C60 results reported earlier (2.72 Å of PBE30, 2.73 Å of BP8629). This also confirms the conclusions of Gagliardi etc.29, that is, if the space around U2 in U2@Cn becomes larger, the two U atoms will tend to dissociate and interact with the inner wall of the cage29. It should be noted that, one of the adsorption sites of the two U atoms in C61 fullerene is different from that of U2@C60, the C atoms in the defective area form a 10-membered ring-like structure and a U atom shifts below the center of this 10-membered ring. This position is different from the case of U2@C60 in which the two U atoms face the centers of two hexagonal rings. In this way, we have arrived at an interesting ‘dissociation-adsorption’ equilibrium state structure of neutral U2 molecule encapsulated in defective C61 fullerene. In Part 1.3 of the supplementary information, we have explained all the details of size comparison between C60 and C61 fullerenes.

Table 1

Full geometry optimization results of U2@C61 (structure c)

System	Method	El. state	Total energy (H)	RU-U (Å)
U2@C61 (c)	BP86	7A′	−3265.91611	2.77
		9A′	−3265.91672	2.80
		11A′	−3265.89768	2.79
	PBE0	7A″	−3262.79316	2.84
		9A′	−3262.79947	2.88
		11A″	−3262.78450	2.71

Table 1 shows the full geometry optimization results of U2@C61 (structure c) calculated at BP86/ANO_ECP60MWB~3-21G level and also with the PBE0 data as a test. From the energy relations of different electronic state we can see both the two XC functionals show that the nonet state is the most stable state and there is no imaginary frequency in the frequency calculation. Hence, structure c is identified as the ground state structure and we express this electronic state as U2@C61 (c, 9A′).

Since infrared (IR) and Raman spectroscopy have been widely used to characterize metallofullerenes1, we showed the IR and Raman spectrums of U2@C61 (c, 9A′) predicted at BP86/ANO_ECP60MWB~3-21G level in this paper. These results are expected to guide experimental research in this area. We first examined the infrared and Raman spectrums of the C60 and U2@C60 and the results are consistent with the previous studies (see supplementary information Part 1.2). Also, we performed the same calculation for defective C61 fullerene. The IR-active vibrational mode in C61 is more complex than that in C60. Unique peaks due to the ad-atom appear in the IR spectrum at the frequencies of 1257.5 and 1353.6 cm−1. These vibrational modes correspond to the asymmetric and symmetric stretching modes of the two C-C bonds of the ad-atom respectively (see supplementary information Part 1.2).

Figure 2 shows the IR and Raman spectra of U2@C61 (c, 9A′) respectively which are predicted at the BP86/ANO_ECP60MWB~3-21G level. We examined all the vibrational mode related to the U atoms and found that the symmetric U-U stretching mode is IR-inactive but Raman-active. So, we can find this signal in the Raman spectrum at the frequency of 127.4 cm−1. This mode can suggest that there is still interaction between the two U atoms in C61 fullerene, the special space size of C61 may not cause them to dissociate completely, although the distance between the two U atoms was enlarged in C61. The asymmetric U-U stretching mode is IR-active and its weak intensity is predicted at a frequency of 142.0 cm−1. This mode explains the interaction between the U atoms and the inner wall of the cage. The unique IR-active peaks arising due to the adatom lie at the frequencies of 1327.5 and 1341.6 cm−1. The vibrational modes still correspond to the asymmetric and symmetric stretching modes of the two C-C bond of the adatom respectively.

Figure 2

Simulated IR and Raman spectrums of U2@C61 (c, 9A′) at BP86/ANO_ECP60MWB~3-21G level.

Furthermore, we compared the vibrational frequencies of U atoms between U2@C60 and U2@C61 systems. From Table 2 we can clearly see that the relative frequency positions in the two EMF systems are reversed. The frequency of symmetric U-U stretching mode in U2@C61 becomes significantly smaller than that in the U2@C60. This interesting phenomenon illustrate that the bond between the two U atoms in C61 is very weak.

Table 2

The frequencies of symmetric and asymmetric U-U stretching modes in U2@C60 and U2@C61 systems. All of these results were calculated at BP86/ANO_ECP60MWB~3-21G level

System	Symmetric U-U stretching	Asymmetric U-U stretching
U2@C60	162.9	157.2
U2@C61 (c, 9A′)	127.4	142.0

The key factor that governs the stability of EMFs is the magnitude of electronic charge transfer between the inner atom/cluster and the cage. Therefore, it is necessary to analyze its charge distribution. We have named the adatom as C61th, the U atom close to the defect as U1st, and the other U atom as U2nd. We have particularly focused on the charge distribution on these three atoms and the relevant data is listed in Table 3. Usually, metal atoms in fullerene, including the U atom, donate electrons to the cage and attain a considerable positive charge. Owing to the different equilibrium positions of the two U atoms in the EMF U2@C61, U1st and U2nd also carry different amount of positive charges. In Table 3, we observe that the U1st is more positively charged than the U2nd. It means that U1st which lies below the defect area transfers more electronic charge to the fullerene cage than the U2nd which faces a six-member ring in the C61 fullerene. This result can be explained as follows. The U1st adsorbs on a 10-membered ring-like structure in contrast to the U2nd which faces a hexagonal structure. U1st interacts with more C atoms than the U2nd. As a result, it donates more electronic charge. The C atoms of the cage gather electronic charge from the inner U atoms and thus get negatively charged except the adatom. Both the two XC-functionals show that the C ad-atom transfers a small amount of electronic charge to the cage and acquires positive charge like the inner U atoms. The ten C atoms of the 10-membered ring mentioned above are comparatively more negatively charged (about −0.297) than other C atoms on the cage, as they have gained electronic charge from both the inner U atom and the outer C adatom.

Table 3

Mulliken charge distribution of U2@C61 (c, 9A′)

System	Method	Atom	Charge (e)
U2@C61 (c, 9A′)	BP86	C61st	0.122
		U1st	2.383
		U2nd	2.316
	PBE0	C61st	0.181
		U1st	2.491
		U2nd	2.410

Electron density map shows the charge distribution in the system more clearly. In Figure 3, the region between the U atoms and the inner wall of the cage shows apparent electron density, suggesting obvious bonding interaction in these areas. The electron density in the space between the U1st and C61st is very low. We have also given the electron density contour line map of U2@C61 (c, 9A′) as Figure S6 in the supplementary information in Part 1.3.

Figure 3

The electron density color-filled map of U2@C61 (c, 9A′).

Molecular orbital (MO) analysis also helps us to understand the electronic structure. The multiplicity of the ground state of U2@C61 is M = 9, which corresponds to eight unpaired electrons. So, we are concerned with the characteristics of LUMOα,HOMOα,HOMOα-1 to HOMOα-7, which are singly-occupied frontier molecular orbitals of the U2@C61(c, 9A′). It can be seen from Figure 4 that the two U atoms provide significant contribution to all these orbitals. Due to impact of adatom-type defect on the fullerene, the shape of the orbitals localized on the U atoms has been deformed and no longer has any symmetry of U2@C6030. Furthermore, we can obtain the natural atomic orbitals (NAO) from the density matrix. These NAOs can be understood as deformed atomic orbitals in the molecular structure. Molecular orbitals can be expanded by these NAOs, which are orthogonal to each other. So, we can calculate the sum of the square of the corresponding orbital coefficients as the contribution to these molecular orbitals. As the 5f and 6d atomic orbitals of U atom are not fully occupied, they usually have a significant contribution to bonding. The 2s and 2p orbitals are the valence orbitals of C atoms. These atomic orbitals are active in intermolecular interaction. Table 4 shows the contributions in the form of a percentage.

Figure 4

The frontier molecular orbitals of the U2@C61 (c, 9A′) at the isovalue of 0.03.

Table 4

The percentages of atomic orbitals active in intermolecular interaction

U2@C61 (c, 9A′)	U1st		U2nd		C61st
	5f	6d	5f	6d	2s	2p
LUMOα	50.59%	0.03%	25.35%	0.16%	0.03%	0.57%
HOMOα	24.97%	0.63%	46.62%	0.24%	0.09%	1.09%
HOMOα-1	40.01%	1.17%	22.54%	0.18%	1.15%	6.30%
HOMOα-2	29.94%	0.31%	40.72%	0.23%	0	0
HOMOα-3	19.52%	0.52%	26.16%	0.02%	0.65%	2.77%
HOMOα-4	12.84%	0.01%	21.55%	0.06%	0.12%	32.08%
HOMOα-5	8.67%	0	8.08%	0.07%	0.11%	43.12%
HOMOα-6	9.96%	0.01%	15.20%	0	0	0.06%
HOMOα-7	25.27%	0.12%	29.89%	0.03%	0.01%	2.89%

According to the shape of the MOs in Figure 4 and corresponding to the data in Table 4, all these MOs reflect the bonding between the two U atoms. The HOMOα-4 and HOMOα-5 reflect the contribution from the 2p electrons of the ad-atom. The 5f electrons of U atoms have a significant contribution to these frontier MOs but the contribution of 6d shell is very small. These results show that the eight singly-occupied electrons mainly come from the 5f of U atoms and the 2p of the ad-atom.

We can directly find out the source of the contribution of the spin state from the spin density analysis. The ground state spin density analysis of U2@C61 (c, 9A′) is shown in Figure 5. The ground state of the neutral U2 molecule is known to have a septet spin state. The six electrons of the same spin are contributed by the 5f orbital of the two U atoms. The ground state of C61 fullerene cages is triplet due to the existence of the defect, and the unpaired electron comes mainly from the 2p of the C adatom (spin = 1.82, supplementary information Part 1.3). We can clearly see from Figure 5 that the nonet ground state of U2@C61 mainly gets the spin contribution from both the U atoms and the adatoms. The unpaired electrons on the U atoms and the C ad-atom are of the same spin (the color of them are all blue). Calculated results show that the spin density of U1st is 3.09, U2nd of 3.30 and the adatom of 1.51. And also, there is a small amount of spin contribution arising from the other C atoms of the fullerene cage. So, the reason of the nonet spin state is that there is a ferromagnetic coupling between the spin of the unpaired electrons of the U atoms and the adatom.

Figure 5

The spin density distributions of the U2@C61 (c, 9A′).

Discussion

In summary, we have constituted an EMF of U2@C61 structure through a spin polarization adatom-type defect on the surface of the typical dimetalloendofullerene U2@C60. DFT theoretical calculations have shown that the defect could expand the local space within the C60 fullerene. This expanded local space and the electronic impact of the ad-atom directionally ‘attract’ one of the two U atoms to the defect area. The bond length in U2 molecule elongates significantly in C61 fullerenes as compared to that of U2@C60. From the vibrational spectra and molecular orbital analysis, bonding interactions are still observed between the two U atoms in C61 fullerene, signifying that the expansion of the local space in the C60 cavity caused by the ad-atom defect does not dissociate the U2 molecule completely. Moreover, this EMF has a unique electronic state unlike the U2@C60. The electronic ground state of U2@C61 has a nonet spin multiplicity with a ferromagnetic coupling between the spin of the unpaired electrons of the U atoms and the adatom.

Because of its unique electronic structure, EMFs have very important promising applications in functional nanomaterials, biomedicines, etc. Our study found that we can further modify the equilibrium structure and electronic structure of EMFs by the introduction of spin polarization defect. This work will contribute to a profound understanding of the complex interplay between the inner metal atoms and outer surface of the fullerene cage in different defect-induced environments or conditions. It can pave possible new ways to regulate the electronic structure of the EMFs to realize its practical applications.

Methods

The relativistic effects and electron correlation must be taken into the consideration in a systematic investigation of actinide-containing systems. A significant advantage of DFT calculations is that accounting for the electron correlation does not increase the CPU time considerably. Therefore, DFT has been widely used in the calculation of actinides31323334. Relativistic pseudopotential method is able to incorporate relativistic effects with a small increase in computational time3536. The impact of the pseudopotential core size of U on chemical calculations based on DFT has been reported3738. Previous studies show that pure DFT methods perform better than hybrid DFT or many-body perturbation methods in describing metal diatomics394041, and therefore pure DFT, such as BP86 and PBE, was used in U2@C60 studies2930. Based on the above, the gradient-corrected BP8642 XC-functional was employed in this paper, and the hybrid functional PBE04344 have also been used for a quality check. Based on testing results and analysis of the structure and charge population characteristics, 3-21G basis set45 was chosen for the carbon atoms4647484950. Similarly, 3-21G basis set for C atoms of fullerene has been applied in a previous lanthanide EMFs study51. Scalar relativistic effects were applied to the uranium atoms. The (14s13p10d8f6g)/[6s6p5d4f3g] named ANO52 valence basis set and ECP60MWB5253 pseudopotential basis set with 60 core-electrons were used. All the calculations have been performed in Gaussian 09 program54.

Author Contributions

X.D. performed most of the numerical simulations. Z.G.W. and X.D. carried out most of the theoretical analysis. C.C., W.Z. and M.S.X. carried out some numerical simulations and theoretical analysis. Z.G.W., P.H. and R.Q.Z. contributed most of the ideas and wrote the paper. All authors discussed the results and commented on the manuscript.