Open access peer-reviewed chapter

Reactivity Indexes and Structure of Fullerenes

Written By

Ernestina Mora Jiménez, Francisco J. Tenorio, David Alejandro Hernández-Velázquez, Jaime Gustavo Rodríguez-Zavala and Gregorio Guzmán-Ramírez

Submitted: 03 April 2017 Reviewed: 21 August 2017 Published: 20 December 2017

DOI: 10.5772/intechopen.70642

From the Edited Volume

Fullerenes and Relative Materials - Properties and Applications

Edited by Natalia V. Kamanina

Chapter metrics overview

1,103 Chapter Downloads

View Full Metrics


The discovery of fullerenes and their production in measurable quantities launched many studies about their reactivity and possible applications. Their peculiar structure opened possibilities for their study, initially replacing carbon atoms with alternative atoms. The surface also offers the possibility of attaching several species and the interior of their hollow structure represents a challenge because of the possibility of confining elements or molecules that may become less stable when attached to the exterior of the cage. These modifications may considerably affect both chemical and physical properties. In this chapter, we propose the encapsulation of 3–10 nitrogen atoms as aggregates inside the C70 cage. We also study the structures and reactivity indexes and the stabilization conferred as a result of being part of the fullerene. These aggregates are mainly of interest because of their possible application as energetic materials.


  • polynitrogen
  • endohedral fullerenes
  • Density Functional Theory
  • reactivity indexes
  • C70
  • energetic materials

1. Introduction

Over the course of time, carbon materials have become important components not only everyday aspects of life, but also vanguard research. The enormous diversity concerning their uses and applications is constantly increasing, representing an area of constant development. Whether as part of a compound or in pure form, carbon has always attracted attention. At least in its purest forms, its structural diversity is extremely attractive (Figure 1).

Figure 1.

Allotropic forms of carbon. From left to right and up to down, diamond, graphite, nanotube and fullerene.

Research on the formation and arrangement of carbon compounds that included long chains of this element in interstellar space lead to the discovery of fullerenes in 1985. Subsequently, they were produced in measurable quantities for use in experiments [1, 2, 3]. Even though their study may have evolved accidentally, insufficient effort has been dedicated to understanding their structure, properties and applications. The best known examples of fullerenes are C60 and C70, both of which consist of 12 pentagons, and only differ concerning the number of hexagons in their structure: 20 in the case of the first and 25 in the second. These differences result in the C60 structure, which is more like a soccer ball and C70, which is more like a rugby ball (Figure 2). While C60 is more often known as a buckyball, the rest of the hollow structures are inscribed to the general fullerene group.

Figure 2.

Comparison of C60 fullerene with a soccer ball and C70 with a rugby ball.

1.1. Types of fullerenes

Research into fullerenes has been extremely diverse, but can be categorized according to the structural aspect to be studied. Correspondingly, we have the following: exohedral fullerenes (with external ligands or molecules), endohedrals (with internal atoms or molecules) and heterofullerenes (with one or more carbon atoms in the cage replaced by heteroatoms).

Exohedral fullerenes have essentially resulted from efforts to “decorate” the surface of the fullerenes; mainly achieved by means of chemical functionalization. These types of fullerenes have proved to be extremely interesting, both for medical applications and for their possible applications in material sciences [4, 5, 6, 7, 8]. With respect to heterofullerenes, the substitution of carbon atoms by boron, nitrogen and silicon has been proposed; in the case of the first two, these have regularly been as -BN- pairs, thus ensuring that isoelectronic species are obtained, although this has not necessarily been the only purpose; in the case of silicon, the main motive has been to search for useful applications in electronics [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]. Last but not least, endohedral fullerenes exist because of the fact that they are hollow; so they have the potential application of confining or protecting molecules, thus acting as carriers or stabilizers. Particularly, those with metallic atoms inside them—endohedral metallofullerenes—have been attractive for chemistry, physics and interdisciplinary areas, such as materials and biological sciences [22, 23, 24, 25].

1.2. Polynitrogenated materials

Considerable interest in nitrogen compounds has emerged, especially those with a high content of this element, because of their particular properties and forms. Apart from the previously known molecular nitrogen (N2), little was known of species with higher nitrogen content, until the synthesis of the azide anion (N3) [26]. However, over time, new polynitrogenated species such as N3, N4, N5 and N5+ have been discovered or posited, whose synthesis has either represented a challenge or turned out to be surprisingly simple [27, 28, 29, 30]. In spite of this, until now, the species with the highest content of nitrogen produced in measurable quantities has contained 5 nitrogen atoms, meaning that in particular the cyclopentazole anion could be considered as a motivator for theoretical studies on its formation, as well as representing a building block towards more complex structures [31, 32, 33, 34, 35]. As species with higher nitrogen content have proved to be very unstable with respect to decomposition to N2, theoretical studies have outnumbered experimental ones. Of these, it has been determined that the preferred forms of N4 and N6 are acyclic, N8 must be formed from an azide-pentazole (N3–N5) and that more complex polynitrogens may exist such as N7, N12, N18, N20 and a great etcetera [36, 37, 38, 39, 40, 41, 42]. Considering larger structures, there is no doubt that polyhedral or fullerene-type structures have provided inspiration, as they are considered as alternative polynitrogenous arrangements with a greater number of nitrogen atoms [18, 19, 21, 41, 43, 44, 45, 46, 47]. In this sense, a hollow structure could in principle be useful for storing internal molecules, which is why it is common to find hollow carbon structures as candidates for confining polynitrogenous species. Currently, the most notable examples have been studies of N@C60, N@C70 or even N2@C60 [48, 49, 50, 51, 52], which have extended the use of the C60 cage as a way of confining up to 16 nitrogen atoms [53, 54]. Also recently, the encapsulation of a polymer nitrogen chain in a carbon nanotube was conceived, as in principle this could be stable up to room temperature [55]. This contrasts with other proposed polynitrogenated phases of nitrogen that have been inspired by the analogy to phosphorus and arsenic, all superior in energy to the “cubic gauche” (CG). These have been proposed for extreme conditions and some experimental evidence has even been found for them [56, 57, 58, 59, 60, 61]. Therefore, confinement represents an alternative for the stabilization of polynitrogenated species that are not necessarily related to high temperatures and pressure.

As our work group has studied fullerenes and their reactivity [7, 16, 62, 63, 64, 65] and it has been proposed that C60 might be an ideal candidate for trapping nitrogen by a polymer [53, 54], the intention now is to show that C70 might also be an alternative way of storing polynitrogenous species, contributing to the development of new materials in the future.


2. Methodology

Initiating with considerable structural diversity, we calculated and modeled endohedral fullerenes; Nn@C70 (n = 3–10). For this, we used the C70 structure that corresponds to the isomer from this composition and complies with the isolated pentagon rule, IPR. The Avogadro visualizer was used as an auxiliary [66]. Given the considerable number of possible isomers, we determined the minimum energy structures in two stages. Initially, we obtained the geometry optimizations for all molecules, using the PM6 method [67]. To ensure that the global minimum for each composition has been identified, the use of search algorithms, such as those inspired by genetic algorithms, is essential, but this is currently beyond the scope of this type of system. However, the diversity of proven structures inspires confidence in our determination of the most important and representative species. Subsequently, both the geometry and the electronic structure were refined for the lower energy isomers of each composition. Likewise, the calculation of vibrational frequencies was undertaken in order to corroborate that the stationary points located on the potential energy surface correspond to a minimum (NImag = 0). All this was undertaken using the hybrid functional B3LYP with the base set 6-311G [68, 69]. As the structures of the polynitrogen species in free state may differ according to charge, this factor was also evaluated, determining the structures for the isomers neutral, cation and anion. All calculations were performed using the Gaussian09 program [70].

We also calculated ionization energies (IE) and electron affinities (EA). These were therefore calculated to reveal the following energy differences: IE = Ecation − Eneutral, EA = Eneutral − Eanion.

Chemical potential was computed by conceptual Density Functional Theory approximation. For an N electron system with an external potential v(r) and total energy E, electronegativity is defined as the partial energy derivative to the number of electrons at constant potential and then by the definition of Mulliken as the mean of IE and EA and the negative of electronegativity is the molecular chemical potential, μ:


Chemical hardness was calculated as defined by Parr and Pearson [71], differentiating the chemical potential to the number of electrons, also at constant energy potential:


Energies were also obtained for the stabilization reactions for polymer species within the C70 structure, applying the formula: ∆E = ΣE(products) − ΣE(reagents) for two possible schemes. In the first, the stabilization of isolated nitrogen atoms could be analyzed through the reaction nN + C70 → Nn@C70. In the second scheme, the reagents are substituted with nitrogen as found in standard state (N2) or detected in experiments (N3, N5, etc).


3. Results

It was previously reported that the methodology consisted of an optimization at the PM6 level. As the refinement and electronic structure were performed at the B3LYP/6-311G level, they are presented exclusively in the following.

For the majority of the structures obtained, neutral and charged species correspond to the same arrangement. In cases where it is not observed, it will be mentioned in due course.

3.1. Nn@C70 structures and reaction energies

N3. The minimum energy structure corresponds to a linear and centered N3, which shows the same structure in charged and neutral systems, see Figure 3. There is also found a second stable structure that may closely resemble the association between the N2 molecule and an “isolated” nitrogen atom. The difference in the relative energies exceeds 40 kcal/mol, leading us to assume that the linear isomer encapsulated in C70 would predominate, if it was ever produced. Distances between nitrogen atoms (1.20 Å) suggest a weaker bond than that observed for N2 (1.11 Å).

Figure 3.

Structure of Nn@C70, where n = 3–10, endohedral fullerenes, calculated using the B3LYP/6-311G method.

Reaction energies were calculated based on various assumptions. First, it is considering only the stabilization due to the encapsulation of nitrogen atoms as reagents. Second, considering structures that can be found at room temperature as starting materials (Table 1).

For n = 32N + N + C70 → N3@C703N + C70 → N3@C702N + N+ + C70 → N3@C70+
N3 + C70 → N3@C701.5N2 + C70 → N3@C70N3+ + C70 → N3@C70+
For n = 43N + N + C70 → N4@C704N + C70 → N4@C703N + N+ + C70 → N4@C70+
N2 + N2 + C70 → N4@C702N2 + C70 → N4@C70N2 + N2+ + C70 → N4@C70+
For n = 54N + N + C70 → N5@C705N + C70 → N5@C704N + N+ + C70 → N5@C70+
5 N-a−764.51−652.59−814.12
5 N-b−745.26−646.04−804.31
5 N-aN5 + C70 → N5@C702.5N2 + C70 → N5@C702N2 + N+ + C70 → N5@C70+
5 N-bN2 + N3 + C70 → N5@C70N2 + N3 + C70 → N5@C70N2 + N3+ + C70 → N5@C70+
For n = 65N + N + C70 → N6@C706N + C70 → N6@C705N + N+ + C70 → N6@C70+
2N2 + N2 + C70 → N6@C703N2 + C70 → N6@C702N2 + N2+ + C70 → N6@C70+
N5 + N + C70 → N6@C70N5 + N+ + C70 → N6@C70N3 + N3+ + C70 → N6@C70+
N3 + N3+ + C70 → N6@C70

Table 1.

Reaction energies due to the encapsulation of nitrogen structures for Nn@C70, where n = 3–6, obtained at B3LYP/6-311G level of theory.

The 3N + C70 → N3@C70 reaction was used to evaluate the role of the C70 box for the stabilization (or not) of 3 isolated nitrogen atoms, compared to placing them within the cage in polymer form as N3. Apparently reactions are energetically favored. This energy difference is also highest for the cation, followed by the anion and finally the neutral (Table 1).

We evaluated energy differences for the reactions presented, while attempting to evaluate the role of the box, when initial material consisted of existing structures (such as N3 and N2, rather than isolated nitrogen atoms). Evidently, as is stated in Table 1, the cation has lower ∆E energy. This molecule would be the most viable option at the time of the reaction, as the carbon cage would effectively stabilize the nitrogen polymer. Similarly, fullerene C70 is able to stabilize and cage the anion polynitrogen, as it also has a negative ΔE value. Contrastingly, in the case of the neutral structure, this is not predicted as a favorable reaction. In spite of this, stabilization of N2 + N in the C70 cage is favorable (∆E = −60.44 kcal/mol).

N4. We analyzed structures consisting of two pairs of N2 molecules, Figure 3, in almost parallel positions (with an approximate distance of 2.67 Å between the two N2 pairs). An arrangement very similar to the neutral one was observed for the charged systems. Unstable structures for charged systems were also found, as in essence they could almost be considered as 4 atoms with no obvious interacting link.

The stabilization of arrangements of 4 nitrogen atoms in the cage was evaluated in terms of reaction energy, results for which are shown in Table 1. Evidently, the encapsulation would be energetically favored. Similar to that observed in the previous isomer, it is predicted that the interaction with the cationic system would be greater than with the anionic or with the neutral one.

Analyzing encapsulation reactions initiating with different materials (such as N2, neutral and charged), evidently the neutral one is not predicted to be energetically favorable. This would lead to a structure which could be very interesting for possible high energy density material (HEDM) applications, a field where metastable compounds and single nitrogen-nitrogen bonds are encouraged. However, the distance between nitrogen atoms (1.11 Å) is characteristic for a triple bond, which discards (at least in principle) this structure as a candidate. For charged systems, the reaction energy is favorable. This gives an idea of the role played by the C70 cage in general, as it stabilizes systems with higher charge density.

N5. The geometries for N5@C70 nitrogen polymers are shown in Figure 3. Captions 5 N-a and 5 N-b symbolize the optimized structures, noting that 5 N-a represents a cycle of 5 nitrogen atoms that is quite similar to cyclopentazole; a system found only as an anion [29, 31]. Distances between nitrogen atoms (~1.34 Å) for the cycle suggest a double-like bond. Contrastingly, 5 N-b isomers essentially consist of two systems; N2 and N3, separated by approximately 2.4 Å. We attempted to optimize a structure similar to N5+ that was predicted theoretically; however, our efforts were in vain, as all the structures presented fragmentation. Distances between nitrogen atoms are, for bended N3 fragment, 1.19 Å, whereas for N2 dimer 1.11 Å.

Geometry acquired by 5 N-a molecule represents the most stable. The energy of neutral 5 N-b is 6.5 kcal/mol greater than 5 N-a, which may be considerable, but would not be as difficult to achieve under experimental conditions. Likewise, the energy difference between the anion and cation systems exceeds 10 kcal/mol, with the cyclic isomers proving the most stable.

The formation energies of the different N5@C70 isomers were calculated by assuming that nitrogen atoms were the material at initiation, indicating at all times that this reaction would be favored. The highest value corresponded to cationic species, followed by anionic species and finally neutral species. Changing the reagents for other species reveals interesting dilemmas. Assuming formation with the pentazolate anion would imply a slightly favored reaction, but when substitute for N2 + N3, these systems become energetically unfavorable. Although this approach might discourage synthesis, it could also act as a motivator, as it would be considered as metastable and a candidate for high energy material, HEDMs.

N6. For N6@C70, the most favored structure corresponds to the arrangement of three N2 molecules, separated 2.4 Å each of the other (Figure 3). For previous reports of N6@C60 [53, 54], the lower energy structure corresponds to the pot-shaped hexagon, in contrast to that presented here. We suggest that this difference is due to the larger volume of C70 as compared to C60, which allows a bigger dispersion within the cage. Cyclic structure with boat-shaped conformation is found 100 kcal/mol greater in energy than reported here and, therefore, is not considered as a possible candidate.

When the reactants are represented by isolated nitrogen atoms, formation energy shows an energetically favorable reaction. However, by changing the reactants for N2, the prediction becomes highly unfavorable. We undertook analysis by changing the reagents, finding that when we proposed charged species as reactants, negative formation energies appeared. We attribute this mainly to the stabilization of charges by the cage.

N7. Of isomers identified in this system, lowest energy consisted of an N5–N2 partnership, similar to that found previously for C60 [53, 54]. Other systems involving N3 with 2N2 or even rings of N7 proved to be stable but considerably higher in energy (Figure 3). Bond distances for nitrogen-nitrogen in the pentagon are 1.31–1.33 Å long, which indicate a double bond. The second energy isomer in order (with N3 and 2N2 fragments) is 28 kcal/mol higher in energy, whereas a distorted heptagon N7 cycle is about 90 kcal/mol higher. This is an indication that the pentagon-N2 association would be the only isomer found in a hypothetical experiment.

Reaction energies show that the stabilization of individual nitrogen atoms is favored, whereas if the reagents are exchanged for others found experimentally such as N5 + N2, this would not be the case. Again this presents a possible other use (as HEDM) for this type of molecule. We also present alternative reactions with energetically favored schemes, which we essentially attribute to the stabilization of charges in the fullerene cage (Table 2).

Nn@C70Anion (−)Neutral (0)Cation (+)
For n = 76N + N + C70 → N7@C707N + C70 → N7@C706N + N+ + C70 → N7@C70+
N2 + N5 + C70 → N7@C703.5N2 + C70 → N7@C703N2 + N+ + C70 → N7@C70+
2N2 + N3+C70 → N7@C70N5 + N2+ + C70 → N7@C70
For n = 87N + N + C70 → N8@C708N + C70 → N8@C707N + N+ + C70 → N8@C70+
N3 + N3 + N2 → N8@C704N2 + C70 → N8@C703.5N2 + N++C70 → N8@C70+
N3 + N5→ N8@C70N5 + N3+ + C70 → N8@C70
For n = 98N + N + C70 → N9@C709N + C70 → N9@C708N + N+ + C70 → N9@C70+
N5 + 2N2 + C70 → N9@C704.5N2 + C70 → N9@C704N2 + N+ + C70 → N9@C70+
3N2 + N3+ C70 → N9@C70N3 + N3++ N3 + C70 → N9@C70
For n = 109N + N + C70 → N10@C7010N + C70 → N10@C709N + N+ + C70 → N10@C70+
N3+3.5N2 + C70 → N10@C705N2 + C70 → N10@C704N2 + N2++C70 → N10@C70+
N5 + 2.5N2 → N10@C70

Table 2.

Reaction energies due to the encapsulation of nitrogen structures for Nn@C70, where n = 7–10, obtained at B3LYP/6-311G level of theory.

N8. The most stable isomer for this composition consists of the assembly between a hexagon and a dimer as two separate units. Bond distances for N─N are around 1.33–1.34 Å, more related to a double bond. Likewise, isomers that involve the formation of a distorted—but unstable—octagon are found with considerably higher energies, making it possible to ensure that this isomer would predominate in an experiment. The nitrogen hexagon shows distances that suggest double bonds, ranging from 1.31 to 1.35 Å, whereas for the unstable octagon they range from 1.26 to 1.39 Å.

Similarly, the reaction energy from the stabilization of nitrogen atoms shows that this would be a favored case, starting with the cation, followed by the anion and finally the neutral. Changing the reactants it is found that metastable species would be formed with molecules present at laboratory conditions, such as N2 and N5. This alternate reaction scheme opens up the possibility to produce HEDM candidates.

N9. The most favored structure for this encapsulation consists of the assembly formed from an N5 ring and two dimers parallel to this. Nitrogen-nitrogen bonds for the cycle are very similar to the corresponding for similar cycles (1.3–1.33 Å) and the distance to N2 is about 2.29 Å. Qualitatively, it is similar to that found for C60, although with considerably greater distances [53, 54]. Stable structures involving ensembles formed with N3 and 3N2 species are 18.78 kcal/mol higher in energy. The fact that it is preferred a more complex structure instead simple ones, indicates us to suggest a more determinant role of the C70 cage in the stabilization of the structure. The negative values of the reaction energy reflect the stabilization of the polynitrogenated species and evidently the reaction is more favored when the system is positively charged, after which follows the negative system and finally the neutral system. The energy difference revealing a positive energy reflects an unfavorable and endergonic reaction, although probably with metastable products. Moreover, by exploring other reaction paths to obtain the isomers of 9 nitrogen atoms, favorable reactions are obtained.

N10. As the amount of nitrogen atoms in the fullerene cage increases, a greater variety of nitrogen polymer structures can be perceived. For this system, the neutral minimum corresponds to a distorted cycle of 10 atoms, which contrasts with that reported previously for C60 [53, 54]. However, the second energy structure for charged systems (with an energy difference not exceeding 8 kcal/mol compared to the most stable) corresponds to that reported for C60. Bonding distances between nitrogen atoms are 1.23–1.43 Å for the lowest energy isomer and 1.24–1.48 Å for the second. The formation energy starting from isolated nitrogen atoms shows negative values, suggesting stabilization on the part of the box for N10, which is highest for the cation, followed by the anion and finally the neutral. The formation reaction that initiates from a variety of reagents, such as N3, N5 and N2, shows positive values, which suggest metastable structures.

3.2. Global reactivity indexes

We obtained the following chemical reactivity indexes: ionization potential, electron affinity, as well as electronic hardness and chemical potential and the results are summarized in Table 3. For ionization potentials, apparently fullerenes encapsulating 3–8 nitrogen atoms fall within a range of 7.4–7.6 eV. Notably, when the fullerenes contain 4 and 6 nitrogen atoms, the force with which the nitrogen endohedral fullerenes retain the electrons is slightly greater than when encapsulating between 3 and 5 nitrogen atoms. When encapsulating between 7 and 8 nitrogen atoms, this force tends to be the same, but with 9 nitrogen atoms within the fullerene, this energy increases to 7.81 eV, implying that more energy is required in order to remove an electron. Contrastingly, when there are 10 nitrogen atoms within the fullerene network, it is easier to remove an electron because the ionization potential decreases. The even/odd effect observed in the first clusters may be explained by the resistance of systems with closed shell (even number of nitrogens) to be ionized, in contrast to those with an open shell (odd number of nitrogens) that contrarily would present a greater tendency to lose an electron, resulting in a closed shell.


Table 3.

Global reactivity indexes for Nn@C70, where n = 7–10, obtained at B3LYP/6-311G level of theory.

The odd/even behavior of endohedral fullerenes, on gaining an electron and forming a negative ion, is also reflected in electron affinity; as apparently when C70 fullerenes encapsulate 3, 5 and 7 nitrogen atoms, the energy released is greater than when 4, 6 and 8 nitrogen atoms are caged. Notably, the energy of 8 nitrogen atoms increases to a lesser degree than 4 and 6 atoms of nitrogen. When encapsulating 9 and 10 nitrogen atoms, the sequence tends to change, because higher energy release corresponds to the 10 nitrogen structure and the lower energy release corresponds to the 9 nitrogen structure.

In terms of electronic chemical potential, endohedral fullerenes: of 4, 6 and 8 N atoms present the highest values for chemical potential, so they tend to accept electrons. Fullerenes of 3, 5 and 7 atoms have the lowest value, so they can donate electrons more easily. However, on reaching 9 and 10 atoms, the sequence changes because the isomer with 9 nitrogen atoms is the one that has the highest value of chemical potential than all other values, whereas the one with 10 nitrogen atoms has the lowest value of all.

With respect to overall hardness, there is a sequence ranging from 3 to 7 nitrogen atoms where the highest hardness values are the even numbers. These, therefore, have the greatest resistance for modifying their electronic density. The lowest hardness values correspond to the odd numbers of nitrogen atoms. At 8 nitrogen atoms, the sequence changes direction: now the odd number of nitrogen atoms is the one that presents the greatest value for global hardness, besides being the molecule that is least reactive in terms of hardness. Considering the case of N10, this structure has the lowest hardness value. Consequently, it can modify its electronic density more easily, so that it represents the most reactive molecule. Previous results for C60 encapsulation [53] showed a reduced HOMO-LUMO gap, as the number of nitrogens increased: 2, 4, 6, 8, 10. This is somewhat similar to the results for hardness, as one of the approximations to this parameter is precisely the HOMO-LUMO gap (Figure 4).

Figure 4.

Frontier orbitals for Nn@C70, where n = 3–10, endohedral fullerenes, calculated using the B3LYP/6-311G method.

3.3. Analysis of frontier orbitals

Previous studies have mentioned that from 2 to 10 nitrogen atoms, there would be a charge transfer from the C60 cage to the nitrogen polymer, explained as resulting from the greater ionization and electronegativity potential of N compared to C. For more than 10 nitrogen atoms, there would be a reverse transfer from the nitrogen to the carbon atoms. This was attributed to reduced space availability and overlapping of orbitals [53]. In this work, the frontier orbitals of the most stable structures were analyzed, showing that there is an appreciable contribution from frontier orbitals to the nitrogen atoms. This implies that there must be areas susceptible to receiving charge density. It is especially interesting that when the number of N atoms within the fullerene increases, the contribution of Nitrogen-centered orbitals for LUMO is even greater. Particularly, in the case of N8 and higher systems, a very considerable contribution of HOMO is observed, and a sure indication that possible interactions with the cage are initiating. This may also relate to the suggestion that greater interaction between the cage and the Nn takes place, as the amount of nitrogen atoms increases. The reaction energy may also result from this, particularly as it is apparent that from N3 to N7 there is a more or less constant increase of 150 kcal/mol in stabilization energy. This reaction energy changes with N8 and higher, suggesting greater participation on the part of the cage and its interaction with Nn.


4. Conclusions

The C70 fullerene is viable for storage and stabilization of nitrogen aggregates of at least 3–10 atoms, without presenting any concurrent structural deformation in the carbon network or obvious interaction (bond-like) between the nitrogen atoms and the carbon of the fullerene.

Evidently, the formation energy from the isolated nitrogen atoms is favored in all cases and from N3 to N8 there is a progressive increase, indicating the role of the cage as an encapsulator and stabilizer of polynitrogenous structures. When we propose the formation of endohedral fullerenes, initiating with materials that can be found under laboratory conditions, we predict that metastable structures will form, a fact that may be of interest in the potential use of these materials as HEDMs. The stabilizing contribution of the cage becomes more evident, as the number of nitrogen atoms within it increases, manifested in a considerable contribution on the part of frontier orbitals as potential charge stabilizers.

As C70 fullerene is a large molecule of approximately 1 nm diameter, nitrogen, polymers of 3–7 nitrogen atoms prefer to remain as nitrogen molecules or as azides. However, from 8 atoms onwards, the fullerene cage begins to have greater interaction with the nitrogen polymers, reflected in the fact that these begin to compact and create more complex nitrogen polymers (in the form of rings).

According to indexes of global chemical reactivity, endohedral fullerenes present an odd/even behavior that corresponds to that expected for open layer species (odd number of electrons: lower ionization potential and hardness, greater electronic affinity and chemical potential) and closed layer (even number of electrons: greater ionization potential and hardness, lower electronic affinity and chemical potential). However, for 8, 9 and 10 encapsulated nitrogen atoms, there is a change in behavior. This coincides with the change in formation energies and the analysis of frontier orbitals, reflecting greater participation of the cage in fullerene behavior. The fact that smaller fullerenes have been able to encapsulate more atoms is an indication that the point of saturation in this structure has not yet been reached, providing an incentive to find new structures with greater complexity, as the degree of encapsulation progresses. These would thus constitute candidates for energetic materials and studies are ongoing.



The authors like to thank PRODEP (formerly PROMEP) for support provided through the 103.5/13/6900 office. FJTR like to thank the University of Guadalajara for authorizing sabbatical leave. DAHV like to thank CONACYT for support provided by the program “Apoyos para la Incorporación de Investigadores Vinculada a la Consolidación Institucional de Grupos de Investigación y/o Fortalecimiento del Posgrado Nacional.” CONACYT through Project 52,827 is also acknowledged.


  1. 1. Kroto HW, Heath JR, O’Brien SC, Curl RF, Smalley RE. C60: Buckminsterfullerene. Nature. 1985 Nov 14;318(6042):162-163
  2. 2. Kratschmer W, Lamb LD, Fostiropoulos K, Huffman DR. Solid C60: A new form of carbon. Nature. 1990 Sep 27;347(6291):354-358
  3. 3. Wragg JL, Chamberlain JE, White HW, Kratschmer W, Huffman DR. Scanning tunnelling microscopy of solid C60/C70. Nature. 1990 Dec 13;348(6302):623-624
  4. 4. Prato M. [60]Fullerene chemistry for materials science applications. Journal of Materials Chemistry. 1997;7(7):1097-1109
  5. 5. Da Ros T, Prato M. Medicinal chemistry with fullerenes and fullerene derivatives. Chemical Communications. 1999;8:663-669
  6. 6. Wudl F. Fullerene materials. Journal of Materials Chemistry. 2002;12(7):1959-1963
  7. 7. Rodriguez-Zavala JG, Tenorio FJ, Samaniego C, Mendez-Barrientos CI, Pena-Lecona FG, Munoz-Maciel J, et al. Theoretical study on the sequential hydroxylation of C-82 fullerene based on Fukui function. Molecular Physics. 2011;109(14):1771-1783
  8. 8. Yan W, Seifermann SM, Pierrat P, Bräse S. Synthesis of highly functionalized C60 fullerene derivatives and their applications in material and life sciences. Organic & Biomolecular Chemistry. 2015;13(1):25-54
  9. 9. Oku T, Nishiwaki A, Narita I. Formation and atomic structures of BnNn (n = 24-60) clusters studied by mass spectrometry, high-resolution electron microscopy and molecular orbital calculations. Physica B: Condensed Matter. 2004 Aug 15;351(1-2):184-190
  10. 10. Golberg D, Bando Y, Stéphan O, Kurashima K. Octahedral boron nitride fullerenes formed by electron beam irradiation. Applied Physics Letters. 1998 Oct 21;73(17):2441-2443
  11. 11. Oku T, Narita I, Nishiwaki A. Formation and structures of B36N36 and Y@B36N36 clusters studied by high-resolution electron microscopy and mass spectrometry. Journal of Physics and Chemistry of Solids. 2004 Mar;65(2-3):369-372
  12. 12. Blase X, Rubio A, Louie S, Cohen M. Stability and Band Gap Constancy of Boron Nitride Nanotubes. EPL (Europhysics Letters). 2007;28(5):335-340
  13. 13. Nishiwaki A, Oku T, Suganuma K. Atomic and electronic structures of endohedral B36N36 clusters with doping elements studied by molecular orbital calculations. Physica B: Condensed Matter. 2004 Jun 15;349(1-4):254-259
  14. 14. Harrison PG. Silicate cages: Precursors to new materials. Journal of Organometallic Chemistry. 1997;542(2):141-183
  15. 15. Bernal Texca FG, Chigo-Anota E, Tepech Carrillo L, Castro M. A DFT study of the electronic and magnetic properties of C36Si24 fullerenes. Computational & Theoretical Chemistry. 2017 [cited 19 May 2017];1103:1-10
  16. 16. Tenorio FJ, Robles J. On the stability and reactivity of C-Si heterofullerenes. International Journal of Quantum Chemistry. 2000 Oct;80(2):220-226
  17. 17. Baei MT. Remove of toxic pyridine from environmental systems by using B12N12 nano-cage. Superlattices and Microstructures. 2013;58:31-37
  18. 18. Wang LJ, Zgierski MZ. Super-high energy-rich nitrogen cluster N60. Chemical Physics Letters. 2003;376:698-703
  19. 19. Wright JS, McKay DJ, DiLabio GA. Dodecahedral molecular nitrogen (N20) and related structures. Journal of Molecular Structure: THEOCHEM. 1998;424(1):47-55
  20. 20. Su B, Feng X, Guo X, Li N. Polynitrogen clusters encapsulated inside B24N24 fullerene-like nanocages: Nanoscale high energy materials studied by density functional theory. Inorganica Chim Acta. 2017;456:128-135
  21. 21. Bliznyuk AA, Shen M, Schaefer HF. The dodecahedral N20 molecule. Some theoretical predictions. Chemical Physics Letters. 1992;198(3):249-252
  22. 22. Shinohara H. Endohedral metallofullerenes. Reports Progress in Physics. 2000;63(6):843
  23. 23. Cerón MR, Maffeis V, Stevenson S, Echegoyen L. Endohedral fullerenes: Synthesis, isolation, mono- and bis-functionalization. Inorganica Chim Acta. 2017. DOI:10.1016/j.ica.2017.03.040. In press
  24. 24. Yeh W-Y. Coordination and reactivity of functionalized fullerenes, open-cage fullerenes, and endohedral metallofullerenes by organometallic complexes. Journal of Organometallic Chemistry. 2015;784:13-23
  25. 25. Popov AA, Yang S, Dunsch L. Endohedral Fullerenes. Chemical Reviews. 2013 Aug 14;113(8):5989-6113
  26. 26. Curtius T. Ueber Stickstoffwasserstoffsäure (Azoimid) N3H. Berichte der Dtsch Chem Gesellschaft. 1890 Jul 1;23(2):3023-3033
  27. 27. Thrush BA. The detection of free radicals in the high intensity photolysis of hydrogen Azide. Proceedings of the Royal Society of London Series A Mathematical Physical Science. 1956 Apr 10;235(1200):143 LP-143147
  28. 28. Cacace F, de Petris G, Troiani A. Experimental detection of Tetranitrogen. Science (80-). 2002 Jan 18;295(5554):-480 LP, 481
  29. 29. Vij A, Pavlovich JG, Wilson WW, Vij V, Christe KO. Experimental detection of the Pentaazacyclopentadienide (Pentazolate) anion, cyclo-N5−. Angewandte Chemie International Edition. 2002 Aug 16;41(16):3051-3054
  30. 30. Christe KO, Wilson WW, Sheehy JA, Boatz JA. N5+: A novel Homoleptic Polynitrogen ion as a high energy density material. Angewandte Chemie International Edition. 1999 Jul 12;38(13-14):2004-2009
  31. 31. Bazanov B, Geiger U, Carmieli R, Grinstein D, Welner S, Haas Y. Detection of Cyclo-N5− in THF solution. Angewandte Chemie International Edition. 2016 Oct 10;55(42):13233-13235
  32. 32. Zhang C, Yang C, Hu B, Yu C, Zheng Z, Sun C. A symmetric co(N5)2(H2O)4·4 H2O high-nitrogen compound formed by cobalt(II) cation trapping of a cyclo-N5− anion. Angewandte Chemie International Edition. 2017 Apr 10;56(16):4512-4514
  33. 33. Geiger U, Haas Y. Preparation of the Cyclopentazole anion in the bulk: A computational study. The Journal of Physical Chemistry. B. 2016 Jul 7;120(26):6208-6214
  34. 34. Choi C, Yoo H-W, Goh EM, Cho SG, Jung Y. Ti(N5)4 as a potential nitrogen-rich stable high-energy density material. The Journal of Physical Chemistry. A. 2016 Jun 23;120(24):4249-4255
  35. 35. Cui J, Zhang Y, Zhao F, Yang J, Shen G, Xu Y. HB (N5)3M (M = Li, Na, K, and Rb): A new kind of pentazolides as HEDMs. Progress in Natural Science. 2009;19(1):41-45
  36. 36. Wang X, Ren Y, Shuai M-B, Wong N-B, Li W-K, Tian A-M. Structure and stability of new N7 isomers. Journal of Molecular Structure: THEOCHEM. 2001;538(1):145-156
  37. 37. Gagliard L, Evangelisti S, Roos BO, Widmark P-O. A theoretical study of ten N8 isomers. Journal of Molecular Structure: THEOCHEM. 1998;428(1):1-8
  38. 38. Hirshberg B, Gerber RB. Decomposition mechanisms and dynamics of N6: Bond orders and partial charges along classical trajectories. Chemical Physics Letters. 2012;531:46-51
  39. 39. Klapötke TM. Ab initio calculations of the open-chain N6 diazide molecule. Journal of Molecular Structure: THEOCHEM. 2000 Mar 17;499(1-3):99-104
  40. 40. Klapötke TM, Harcourt RD. The interconversion of N12 to N8 and two equivalents of N2. Journal of Molecular Structure: THEOCHEM. 31 May 2001;541(1-3):237-242
  41. 41. Gu JD, Chen KX, Jiang HL, Chen JZ, Ji RY, Ren Y, Tian AM. N18: A computational investigation. Journal of Molecular Structure: THEOCHEM. 1998 Mar 23;428(1-3):183-188
  42. 42. Sturdivant SE, Nelson FA, Strout DL. Trends in stability for N18 cages. The Journal of Physical Chemistry. A. 2004 Aug 1;108(34):7087-7090
  43. 43. Balasubramanian K. Nuclear spin multiplets for the dodecahedral N20 cluster. Chemical Physics Letters. 1993;202(3):271-277
  44. 44. Ha T-K, Suleimenov O, Nguyen MT. A quantum chemical study of three isomers of N20. Chemical Physics Letters. 1999;315(5):327-334
  45. 45. Liu F-L, Yang F, Zhang L-X. Computational note on C10N30: Two intact dodecahedral C5N15 cages linked by five carbon–carbon single bonds. Journal of Molecular Structure: THEOCHEM. 2010;950(1-3):98-99
  46. 46. Zhou H, Beuve M, Yang F, Wong N-B, Li W-K. Theoretical investigation on the cylinder-shaped N66 cage. Computational & Theoretical Chemistry. 2013 Feb 1;1005:68-74
  47. 47. Zhou H, Wong N-B, Tian A. Theoretical study on the cylinder-shaped N78 cage. Journal of Molecular Graphics & Modelling. 2006 Dec;25(4):578-583
  48. 48. Knapp C, Dinse K-P, Pietzak B, Waiblinger M, Weidinger A. Fourier transform EPR study of N@C60 in solution. Chemical Physics Letters. 1997;272(5):433-437
  49. 49. Dinse K-P. EPR investigation of atoms in chemical traps. Physical Chemistry Chemical Physics. 2002;4(22):5442-5447
  50. 50. Meyer C, Harneit W, Lips K, Weidinger A, Jakes P, Dinse KP. Alignment of the Endohedral fullerenes N@C60 and N@C70 in a liquid-crystal matrix. Physical Review A. 2002;65:61201
  51. 51. Dietel E, Hirsch A, Pietzak B, Waiblinger M, Lips K, Weidinger A, et al. Atomic nitrogen encapsulated in fullerenes: Effects of cage variations. Journal of the American Chemical Society. 1999;121:2432
  52. 52. Franco L, Ceola S, Corvaja C, Bolzonella S, Harneit W, Maggini M. Synthesis and magnetic properties of N@C60 derivatives. Chemical Physics Letters. 2006;422(1):100-105
  53. 53. Sharma H, Garg I, Dharamvir K, Jindal VK. Structure of polynitrogen clusters encapsulated in C60: A density functional study. Journal of Physical Chemistry C. 20 May 2010;114(19):9153-9160
  54. 54. Sharma H, Garg I, Dharamvir K, Jindal VK. Structural, electronic, and vibrational properties of C60−nNn (n = 1−12). The Journal of Physical Chemistry. A. 2009 Aug 6;113(31):9002-9013
  55. 55. Zhong Y, Jaidann M, Zhang Y, Zhang G, Liu H, Ioan Ionescu M, et al. Synthesis of high nitrogen doping of carbon nanotubes and modeling the stabilization of filled DAATO@CNTs (10,10) for nanoenergetic materials. Journal of Physics and Chemistry of Solids. 2010;71(2):134-139
  56. 56. Mitáš L, Martin RM. Quantum Monte Carlo of nitrogen: Atom, dimer, atomic, and molecular solids. Physical Review Letters. 1994 Apr 11;72(15):2438-2441
  57. 57. Mailhiot C, Yang LH, McMahan AK. Polymeric nitrogen. Physical Review B. 1992;46(22):14419-14435
  58. 58. Lewis SP, Cohen ML. High-pressure atomic phases of solid nitrogen. Physical Review B. 1992;46(17):11117-11120
  59. 59. Martin RM, Needs RJ. Theoretical study of the molecular-to-nonmolecular transformation of nitrogen at high pressures. Physical Review B. 1986;34(8):5082-5092
  60. 60. Eremets MI, Gavriliuk AG, Trojan IA, Dzivenko DA, Boehler R. Single-bonded cubic form of nitrogen. Nature Materials. 2004;3(8):558-563
  61. 61. Mattson WD, Sanchez-Portal D, Chiesa S, Martin RM. Prediction of New Phases of Nitrogen at High Pressure from First-Principles Simulations. Physical Review Letters. 2004;93(12):125501
  62. 62. Rodríguez-Zavala JG, Barajas-Barraza RE, Padilla-Osuna I, Guirado-López RA. Hydration behaviour of polyhydroxylated fullerenes. Journal of Physics B: Atomic, Molecular and Optical Physics. 2011;44(20):205104
  63. 63. Rodríguez-Zavala JG, Padilla-Osuna I, Guirado-López RA. Modifying the endohedral la-ion location in both neutral and positively charged polyhydroxylated metallofullerenes. Physical Review B. 2008 Oct 24;78(15):155426
  64. 64. Rodríguez-Zavala JG, Guirado-López RA. Stability of highly OH-covered C60 fullerenes: Role of Coadsorbed O impurities and of the charge state of the cage in the formation of carbon-opened structures. The Journal of Physical Chemistry. A. 2006 Aug 1;110(30):9459-9468
  65. 65. Rodríguez-Zavala JG, Guirado-López RA. Structure and energetics of polyhydroxylated carbon fullerenes. Physical Review B. 2004 Feb 23;69(7):75411
  66. 66. Hanwell MD, Curtis DE, Lonie DC, Vandermeersch T, Zurek E, Hutchison GR. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. Journal of Cheminformatics. 2012;4(1):17
  67. 67. Stewart JJP. Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements. Journal of Molecular Modeling. 2007;13(12):1173-1213
  68. 68. Becke AD. A new mixing of Hartree–Fock and local density-functional theories. The Journal of Chemical Physics. 1993;98(2):1372-1377
  69. 69. Krishnan R, Binkley JS, Seeger R, Pople JA. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. The Journal of Chemical Physics. 1980 Jan 1;72(1):650-654
  70. 70. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA, Nakatsuji H, Caricato M, Li X, Hratchian HP, Izmaylov AF, Bloino J, Zheng G, Sonnenberg JL, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Vreven T, Montgomery JA Jr, Peralta JE, Ogliaro F, Bearpark M, Heyd JJ, Brothers E, Kudin KN, Staroverov VN, Kobayashi R, Normand J, Raghavachari K, Rendell A, Burant JC, Iyengar SS, Tomasi J, Cossi M, Rega N, Millam JM, Klene M, Knox JE, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Martin RL, Morokuma K, Zakrzewski VG, Voth GA, Salvador P, Dannenberg JJ, Dapprich S, Daniels AD, Farkas Ö, Foresman JB, Ortiz JV, Cioslowski J, Fox DJ. Gaussian, Inc., Wallingford CT, 2009. Gaussian 09, Revision C.01
  71. 71. Parr RG, Pearson RG. Absolute hardness: Companion parameter to absolute electronegativity. Journal of the American Chemical Society. 1983 Dec 1;105(26):7512-7516

Written By

Ernestina Mora Jiménez, Francisco J. Tenorio, David Alejandro Hernández-Velázquez, Jaime Gustavo Rodríguez-Zavala and Gregorio Guzmán-Ramírez

Submitted: 03 April 2017 Reviewed: 21 August 2017 Published: 20 December 2017