Reaction energies due to the encapsulation of nitrogen structures for N
Abstract
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.
Keywords
- 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).
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.
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; N
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
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
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).
N | Anion | Neutral | Cation |
---|---|---|---|
For | 2N + N− + C70 → N3@C70− | 3N + C70 → N3@C70 | 2N + N+ + C70 → N3@C70+ |
−495.82 | −383.89 | −544.51 | |
N3− + C70 → N3@C70− | 1.5N2 + C70 → N3@C70 | N3+ + C70 → N3@C70+ | |
−39.13 | 101.29 | −106.97 | |
For | 3N + N− + C70 → N4@C70− | 4N + C70 → N4@C70 | 3N + N+ + C70 → N4@C70+ |
−678.69 | −606.42 | −764.55 | |
N2 + N2− + C70 → N4@C70− | 2N2 + C70 → N4@C70 | N2 + N2+ + C70 → N4@C70+ | |
−69.85 | 40.49 | −145.73 | |
For | 4N + N− + C70 → N5@C70− | 5N + C70 → N5@C70 | 4N + N+ + C70 → N5@C70+ |
5 N-a | −764.51 | −652.59 | −814.12 |
5 N-b | −745.26 | −646.04 | −804.31 |
5 N-a | N5− + C70 → N5@C70− | 2.5N2 + C70 → N5@C70 | 2N2 + N+ + C70 → N5@C70+ |
−0.13 | 156.04 | −167.21 | |
5 N-b | N2 + N3− + C70 → N5@C70− | N2 + N3 + C70 → N5@C70 | N2 + N3+ + C70 → N5@C70+ |
34.89 | 73.75 | −43.31 | |
For | 5N + N− + C70 → N6@C70− | 6N + C70 → N6@C70 | 5N + N+ + C70 → N6@C70+ |
−881.66 | −807.24 | −965.47 | |
2N2 + N2− + C70 → N6@C70− | 3N2 + C70 → N6@C70 | 2N2 + N2+ + C70 → N6@C70+ | |
50.62 | 163.12 | −23.19 | |
N5− + N + C70 → N6@C70− | N5− + N+ + C70 → N6@C70 | N3 + N3+ + C70 → N6@C70+ | |
−117.29 | −385.00 | −131.54 | |
N3− + N3+ + C70 → N6@C70 | |||
−255.14 |
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).
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.
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.
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.
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).
N | Anion (−) | Neutral (0) | Cation (+) |
---|---|---|---|
For | 6N + N− + C70 → N7@C70− | 7N + C70 → N7@C70 | 6N + N+ + C70 → N7@C70+ |
−979.94 | −866.37 | −1026.35 | |
N2 + N5− + C70 → N7@C70− | 3.5N2 + C70 → N7@C70 | 3N2 + N+ + C70 → N7@C70+ | |
107.89 | 265.72 | −55.99 | |
2N2 + N3−+C70 → N7@C70− | N5− + N2+ + C70 → N7@C70 | ||
123.67 | −148.76 | ||
For | 7N + N− + C70 → N8@C70− | 8N + C70 → N8@C70 | 7N + N+ + C70 → N8@C70+ |
−1615.18 | −1528.77 | −1688.71 | |
N3 + N3− + N2 → N8@C70− | 4N2 + C70 → N8@C70 | 3.5N2 + N++C70 → N8@C70+ | |
187.33 | 392.54 | 70.88 | |
N3 + N5−→ N8@C70− | N5− + N3+ + C70 → N8@C70 | ||
171.55 | −41.48 | ||
For | 8N + N− + C70 → N9@C70− | 9N + C70 → N9@C70 | 8N + N+ + C70 → N9@C70+ |
−1046.85 | −995.79 | −1149.58 | |
N5− + 2N2 + C70 → N9@C70− | 4.5N2 + C70 → N9@C70 | 4N2 + N+ + C70 → N9@C70+ | |
364.43 | 459.75 | 144.23 | |
3N2 + N3−+ C70 → N9@C70− | N3− + N3++ N3 + C70 → N9@C70 | ||
380.20 | −47.35 | ||
For | 9N + N− + C70 → N10@C70− | 10N + C70 → N10@C70 | 9N + N+ + C70 → N10@C70+ |
b | −1143.41 | −981.87 | −1192.45 |
a | −1150.42 | −996.12 | −1184.66 |
N3−+3.5N2 + C70 → N10@C70− | 5N2 + C70 → N10@C70 | 4N2 + N2++C70 → N10@C70+ | |
a | 438.36 | 621.14 | 396.74 |
N5− + 2.5N2 → N10@C70− | |||
422.59 |
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.
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.
I | A | μ | η | |
---|---|---|---|---|
N3@C70 | 7.51 | 4.49 | −6 | 1.51 |
N4@C70 | 7.61 | 2.77 | −5.19 | 2.42 |
N5@C70 | 7.47 | 4.49 | −5.98 | 1.49 |
N6@C70 | 7.61 | 2.87 | −5.24 | 2.37 |
N7@C70 | 7.54 | 4.56 | −6.05 | 1.49 |
N8@C70 | 7.54 | 3.39 | −5.46 | 2.08 |
N9@C70 | 7.81 | 1.85 | −4.83 | 2.98 |
N10@C70 | 6.3 | 6.33 | −6.31 | −0.02 |
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).
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 N
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.
Acknowledgments
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.
References
- 1.
Kroto HW, Heath JR, O’Brien SC, Curl RF, Smalley RE. C60: Buckminsterfullerene. Nature. 1985 Nov 14; 318 (6042):162-163 - 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.
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.
Prato M. [60]Fullerene chemistry for materials science applications. Journal of Materials Chemistry. 1997; 7 (7):1097-1109 - 5.
Da Ros T, Prato M. Medicinal chemistry with fullerenes and fullerene derivatives. Chemical Communications. 1999; 8 :663-669 - 6.
Wudl F. Fullerene materials. Journal of Materials Chemistry. 2002; 12 (7):1959-1963 - 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.
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.
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.
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.
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.
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.
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.
Harrison PG. Silicate cages: Precursors to new materials. Journal of Organometallic Chemistry. 1997; 542 (2):141-183 - 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.
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.
Baei MT. Remove of toxic pyridine from environmental systems by using B12N12 nano-cage. Superlattices and Microstructures. 2013; 58 :31-37 - 18.
Wang LJ, Zgierski MZ. Super-high energy-rich nitrogen cluster N60. Chemical Physics Letters. 2003; 376 :698-703 - 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.
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.
Bliznyuk AA, Shen M, Schaefer HF. The dodecahedral N20 molecule. Some theoretical predictions. Chemical Physics Letters. 1992; 198 (3):249-252 - 22.
Shinohara H. Endohedral metallofullerenes. Reports Progress in Physics. 2000; 63 (6):843 - 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.
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.
Popov AA, Yang S, Dunsch L. Endohedral Fullerenes. Chemical Reviews. 2013 Aug 14; 113 (8):5989-6113 - 26.
Curtius T. Ueber Stickstoffwasserstoffsäure (Azoimid) N3H. Berichte der Dtsch Chem Gesellschaft. 1890 Jul 1; 23 (2):3023-3033 - 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.
Cacace F, de Petris G, Troiani A. Experimental detection of Tetranitrogen. Science (80-). 2002 Jan 18; 295 (5554):-480 LP, 481 - 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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Balasubramanian K. Nuclear spin multiplets for the dodecahedral N20 cluster. Chemical Physics Letters. 1993; 202 (3):271-277 - 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.
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.
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.
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.
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.
Dinse K-P. EPR investigation of atoms in chemical traps. Physical Chemistry Chemical Physics. 2002; 4 (22):5442-5447 - 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.
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.
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.
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.
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.
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.
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.
Mailhiot C, Yang LH, McMahan AK. Polymeric nitrogen. Physical Review B. 1992; 46 (22):14419-14435 - 58.
Lewis SP, Cohen ML. High-pressure atomic phases of solid nitrogen. Physical Review B. 1992; 46 (17):11117-11120 - 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.
Eremets MI, Gavriliuk AG, Trojan IA, Dzivenko DA, Boehler R. Single-bonded cubic form of nitrogen. Nature Materials. 2004; 3 (8):558-563 - 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.
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.
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.
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.
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.
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.
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.
Becke AD. A new mixing of Hartree–Fock and local density-functional theories. The Journal of Chemical Physics. 1993; 98 (2):1372-1377 - 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.
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.
Parr RG, Pearson RG. Absolute hardness: Companion parameter to absolute electronegativity. Journal of the American Chemical Society. 1983 Dec 1; 105 (26):7512-7516