Energy barriers of Be and B in passing through the C60 fullerene cage

A. V. Bibikov,A. V. Nikolaev,P. V. Borisyuk,E. V. Tkalya

2023-12-21

Abstract:We have studied the potential barriers for the penetration of atomic beryllium or boron inside the C60 fullerene by performing ab initio density functional theory (DFT) calculations with three variants for the exchange and correlation: B3LYP (hybrid functional), PW91 and PBE. Four principal trajectories to the inner part of C60 for the penetrating atom have been considered: through the center of six-member-carbon ring (hexagon), five-member-carbon ring (pentagon), and also through the center of the double C-C bond (D-bond) and the center of the single C-C bond (S-bond). Averaging over the three DFT variants yields the following barriers for beryllium penetrating inside a deformable fullerene: 3.2 eV (hexagon), 4.8 eV (S-bond), 5.3 eV (D-bond), 5.9~eV (pentagon). These barriers correspond to the slow and adiabatic penetration of Be, in contrast to the fast (non-adiabatic) penetration through the rigid cage of C60 resulting in 5.6 eV (hexagon), 16.3 eV (pentagon), 81.8 eV (S-bond) and 93.4 eV (D-bond). The potential barriers for the boron penetrating inside deformable/rigid C60 are: 3.7/105.4 eV (D-bond), 4.0/86.8 eV (S-bond), 4.7/7.8 eV (hexagon), 6.8/14.0 eV (pentagon). The potential barriers for Be and B escaping from the inner part of C$_{60}$ are higher by the value of 0.84 eV for Be and 0.81 eV for B. The considerable reduction of the potential barriers for the deformable fullerene is ascribed to the formation of the Be-C and B-C bonds. We discuss the difference between Be and B, compare three variants of DFT, and analyze the role of the dispersion interaction.

Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

### What problems does this paper attempt to solve?

This paper aims to study **the energy barriers for beryllium (Be) and boron (B) atoms to penetrate the C60 fullerene cage**. Specifically, the authors explored the energy barriers of these atoms entering or leaving the C60 cage along different paths through first - principles density - functional theory (DFT) calculations. The main focuses are as follows:

1. **Energy barriers of different paths**:
- Through four main paths: the six - membered carbon ring (hexagon), the five - membered carbon ring (pentagon), the double - bond center (D - bond), and the single - bond center (S - bond).

2. **Differences between rigid and deformable C60 cages**:
- The differences in the energy barriers for Be and B penetration in rigid (non - deformable) and deformable C60 cages were compared. The energy barriers in the rigid cage are significantly higher than those in the deformable cage.

3. **Influence of different DFT functionals**:
- Three different exchange - correlation functionals (B3LYP, PW91, and PBE) were used for calculations to evaluate their influence on the results.

4. **Role of dispersion interactions**:
- The differences in the results with and without van der Waals (vdW) interactions were analyzed.

5. **Differences between Be and B**:
- The differences in the behavior of Be and B during the penetration process were explored, especially their ability to form chemical bonds with carbon atoms.

6. **Potential applications**:
- The potential applications of these research results in the medical field were discussed, such as the design of new radiopharmaceuticals for boron neutron capture therapy (BNCT).

### Key formulas

- **Binding energy formula**:
\[
E_b = E(AC_{60}) - E(C_{60}) - E(A)
\]
where \( A \) represents a Be or B atom, \( E(A) \) is the energy of the Be or B atom, and \( E(C_{60}) \) is the energy of the unperturbed isolated C60 fullerene.

- **Relationship between escape energy barrier and penetration energy barrier**:
\[
\Delta E_0=\Delta E_{\infty}-E_c
\]
where \( \Delta E_0 \) is the energy barrier for escaping from the inside, \( \Delta E_{\infty} \) is the energy barrier for penetrating from the outside to the inside, and \( E_c \) is the binding energy when Be or B is located at the center of C60.

### Overview of results

- **For deformable C60**:
- The energy barriers for Be and B through the six - membered ring are the lowest, being 3.2 eV and 3.7 eV respectively.
- The energy barriers of the double - bond and single - bond paths are significantly reduced, especially the energy barrier for B through the double - bond path is only 3.0 eV.

- **For rigid C60**:
- The energy barriers increase significantly, especially the energy barriers through the double - bond and single - bond paths are as high as 81.8 eV and 93.4 eV (for Be), and 105.4 eV and 86.8 eV (for B) respectively.

These results indicate that **the displacement of carbon atoms and the formation of local chemical bonds significantly reduce the energy barriers for Be and B to penetrate the C60 cage**, thus providing a theoretical basis for the design of new radiopharmaceuticals.