Skeletal Rearrangements of the C240 Fullerene: Efficient Topological Descriptors for Monitoring Stone–Wales Transformations

Stone–Wales rearrangements of the fullerene surface are an uncharted field in theoretical chemistry. Here, we study them on the example of the giant icosahedral fullerene C240 to demonstrate the complex chemical mechanisms emerging on its carbon skeleton. The Stone–Wales transformations of C240 can produce the defected isomers containing heptagons, extra pentagons and other unordinary rings. Their formations have been described in terms of (i) quantum-chemically calculated energetic, molecular, and geometric parameters; and (ii) topological indices. We have found the correlations between the quantities from the two sets that point out the role of long-range topological defects in governing the formation and the chemical reactivity of fullerene molecules.


Introduction
Most chemical studies on fullerenes are focused on the reactions and compounds of C60 and C70, the two most abundant fullerenes. Nevertheless, the diversity of the sizes and shapes within the fullerenes' family is much wider. Indeed, there are many fullerene structures, smaller [1] and larger [2], which are produced in smaller amounts. The chemistry of such small and large fullerenes has growing experimental and theoretical interests.
We pay attention to the giant C240 fullerene with icosahedral symmetry. Its structure is the next after C60 in the series of icosahedral Goldberg polyhedra [3][4][5] having 60a 2 vertices (where a is the integer). In the fundamental aspect, such chemical structures are very close to their mathematical models, so direct bridges between the chemistry and mathematics could be built, which means the possible correlations between the structural and topological (mathematical) indices and molecular (chemical) properties of the molecules [4][5][6]. It is noteworthy that the C240 structure obtained with quantum chemical calculations corresponds to the idealized one, i.e. the highest symmetry (Ih) allowed by the topology of the molecule (there is a possibility of the Jahn-Teller symmetry reduction in the case of other icosahedral fullerenes, e.g. C20 and C80) [6]. In practical aspects, such giant cages as C240 attract as they have a large size and may be classified as nano-objects, but they are still molecules. Due to their size, giant fullerenes are considered as nano-vessels for gas storage [7] and the closed structure makes C240 a molecular Faraday cage screening the guest atoms/molecules from the external electric fields [8]. Note that C240 (Ih) has not been yet synthesized nor isolated. Under arcdischarge fullerene synthesis, it should be produced in extremely minute quantities (it is also potentially insoluble compound due to its size; this makes an obstacle of its extraction) [9,10]. However, its complexes like C60@C240, C240@C560, and C60@C240@C560 have been obtained and visualized with atom-force microscopy [11][12][13][14]. As follows from these studies, the C240 molecule is most likely icosahedral. The thermodynamic [5,15,16] and dielectric properties [8,[17][18][19] of C240 and other giant fullerenes have been studied with relevant computational techniques.
Giant fullerenes relate to the processes of classical fullerene formation [20]. They are proposedly formed from the C2 gas and convert to C60 and C70 [21]. The presence of a large cavity (much larger than in C60) makes giant fullerenes chemically unstable [21] and sensitive to various transformations. The Stone-Wales isomerization is one of such topological processes that may occur on the fullerene surface [22,23]. Since 1986 [24], it is well-known that Stone-Wales (SW) rotations have a theoretical role in connecting the isomers of a given Cn fullerene.
Under the action of the SW5|6 operator (Figure 1a), the 1812 isomers of C60 are grouped in 13 sets. The larger one includes 1709 cages connected to C60 (Ih) through one or more SW5|6 transformations whereas 31 isomers remain separated from the rest of the isomers [24]. Later, the introduction of the so-called generalized Stone-Wales transformations (Figure 1b) overcame this limitation offering a method of generating the complete C60 isomeric space starting from just one C60 isomers [25]. The instrumental role of chemical graphs in studying fullerene rearrangements was immediately evident [25]. The description of fullerene molecule Cn as a planar graph made by n nodes (carbon atoms), = 3 /2 edges (carbon-carbon bonds) forming f5 = 12 pentagons and f6 = n/2 -10 hexagons is a convenient and powerful tool that is also suitable for computerized algorithms. Chemical graph represents the topological skeleton of the molecule and allows a clear picture of the Stone-Wales isomerization mechanism. In Figure 1a, the SW5|6 transformation (the so-called pyracylene rearrangement) that is capable of bridging Cn isomers with different symmetries is shown in the direct (bottom) and dual graph (top). In the rotation of the central bond is obtained by cutting the two graph edges indicated by the red arrows and then by connecting, pairwise, the green and the orange nodes. This topological transformation is perfectly isomeric since the numbers of nodes, edges, and faces are conserved. The generation of the SW5|6 rotation became even simpler when occurring on the dual fullerene graph composed by = 2 + 2 dual nodes, which correspond to the faces of the direct Cn fullerene and = 3 /2 edges. Both fullerene graphs, direct and dual , have the same number of edges . The top part of Figure 1a illustrates the pyracylene transformation in the dual graph . In this case the SW5|6 rearrangement is obtained just by rotating the central edge (red arrow) evidencing a general feature of the SW transformations the dual space :  In the dual space, a SW rotation consists in the reversible rotation of just one dual bond. This straightforward mechanism allows the generation of extended topological defects, so-called generalized Stone-Wales transformations (gSW) that modify the fullerene surfaces in a peculiar fashion. Figure 1b shows the typical gSW involving two 5|6 pairs (shown in Figure 1a) and η pairs of the internal faces (white circles) included between those 5|6. The white circles represent indifferently hexagons or pentagons. The surface transformation of the dual fullerene ( Figure 1b) corresponds therefore to the gSWη=3 of size η = 3 whereas the generalized rotation gSWη = 0 of order η = 0 represents instead to the pyracylene rearrangement SW5|6. Operators like gSW depend on integer parameter η ( Figure 1b) and constitute only one set of the infinite class S of reversible non-local topological rearrangements. The extension (or the size) of the gSW defect is expressed by integer η.  The gSW size η corresponds to the number of pairs of internal faces (white circles) included between the two 5|6 pairs. All SW operators represent full-isomeric rearrangements of the fullerene surface, which therefore preserve the numbers of direct nodes n, dual nodes , bonds B and faces. They also preserve the type of internal faces included in a dislocation dipole. Figure 1 gives nice examples of these topological modifications. The generalized Stone-Wales reversible rotation gSW of order η = 3 in the dual space. In both cases, blue (yellow) elements represent hexagonal (pentagonal) rings; white circles may be hexagonal or pentagonal faces indifferently and they are preserved before/after the gSW rotation.
In S another interesting group of operators is represented by the specific transformations that generate and then propagate 5|7 pairs in a hexagonal network [26,27]. In a regular hexagonal mesh, this original mechanism produces a peculiar linear extended defect, the so-called SW wave (SWw) with variable length η. Figure 2 describes in the dual graph the sequence of specific SW operators that create (SW6|6) and move (SW6|7) that wave in a graphenic lattice. The first rotation SW6|6 ( Figure  2a) of the arrowed dual edge creates the 5|7 double pairs (Figure 2b). The second rotation SW6|7 separates them by the insertion of a couple of hexagons (Figure 2c). At this point, the propagation of the 5|7 pair has begun and successive SW6|7 transformations propagate it in the lattice producing the topological SW wave. This topological defect SWw is known in the literature also under the name of dislocation dipole, whose length η corresponds to the minimal number of hexagon pairs between the 5|7 pairs. Figure 2 shows the dislocation dipoles (SWw) with η = 0, 1, and 2.
Theoretically, the existence of extended topological defects such as SWw with η ≥ 3 has been documented by molecular dynamics and ab initio studies [28,29] reporting about the presence of 5|7 dislocations in pristine graphene and at the grain boundaries of polycrystalline graphenic lattices. In the latter case, the authors state that these linear defects cannot be annealed by local reorganization of the honeycomb mesh. Specifically, work [28] affirmed that "neither an isolated pentagon, or a heptagon, or their pair 5|7" cannot be corrected by lattice rearrangements -an erroneous conclusion not considering that the SWw mechanism allows the 5|7 dislocations to be annealed by moving the 5|7 pair backward by applying the SW transformations ( Figure 2) in the reverse direction.
Electronic properties of honeycomb planar systems indicate that single heptagon-pentagon dislocations are stable defects whereas 5|7 adjacent pairs are dynamically unstable [29]. The energy Ef of the formation of the 5|7 double pairs (Figure 2b) determined by the extended Hückel [30] and molecular mechanical simulations [31] depends on multiple structural factors: (i) the size of the honeycomb flat lattice around the defect, with Ef that assumes the minimum value Ef ≅ 6 eV for "hexagonal layers" made by 9,000 carbon atoms [31]; (ii) the presence a cylindrical curvature in the honeycomb mesh that facilitates furthermore the creation of the dislocation dipole by lowering the energy barrier to Ef ≅ 3 eV for carbon nanotubes. Formation energy of the 5|7 double pair drops even more, down to Ef ≅ 2 eV, when interstitial defects or ad-atoms are present in the networks [32,33]. Pentagon-heptagon pairs have been moreover detected by high-resolution TEM studies [34][35][36][37] whereby energetic particles, such as electrons and ions, generate 5|7 pairs in graphite layers or singlewalled CNTs because of atom displacements.
We consider that both appearing elements favor the creation of dislocation dipoles, i.e. curvature and presence of hexagonal regions of variable size, are specific characteristics of fullerenes' morphology. Therefore, the aim of this work is addressed, mainly from a topological corner, to the following basal theoretical questions:  In which ways is the fullerene surface modified by SWw topological defects?  Is the creation of SWw defects energetically favored? In the following, the first attempt for answering the above points is made by taking into account the C240 (Ih) molecule as a case study. The reasons for this choice are the high symmetry of the molecule and the previous of well-documented ab initio [16] and topological studies [27] about its stability. The next paragraph is devoted to the description of the mechanisms for the formation of dislocation dipoles on the C240 (Ih) surface; then the energetic and topological considerations complete this work.

SW Waves on the C240 (Ih) Fullerene
The theoretical evidences of the stability of 5|7 defects in various hexagonal systems [28][29][30][31][32][33][34][35][36] suggest the possibility to create extended defects even on the surface of a fullerene large enough to have graphene-like zones tiled only with hexagons. This section provides a concise introduction to the topological wave-like mechanism for the creation and the diffusion (or annihilation) of the pentagon-heptagon linear defect on the fullerene surface. The present investigation, combining extended topological rearrangements of the carbon networks and specific modifications of fullerene topology, sheds new light on possible mechanism that could be a base of fullerene formation and chemical properties. Figure 3 shows the icosahedral C240 cage with f6 = 110 hexagonal rings and the symmetryindependent rearrangements of type SW6|6, which may take place on its surface (images sourced from the web [38]). Due to the icosahedral structure of the C240 molecule, one has in fact just three symmetry-independent ways A, B, and C of selecting quartets of hexagons admitting the SW6|6 rotation. The three regions A, B, and C are shown in Figure 3a prior to applying the SW operator according to the mechanisms illustrated above ( Figure 2). The SW6|6 operator rotates the dual bonds pictured in black in each diamond in Figure 3a. structure generated by applying the SW6|6 rotation to quartet A in (a). The double pair 5|7 and SWw with length η = 0 are shown. Yellow (green) polygons represent pentagonal (heptagonal) rings. The pair 6|6 is ready for the SW6|7 rotation and swap place with the nearby pair 5|7.
Our density functional theory (DFT) calculations show that, among the three isomers of the C240 (Ih) fullerene generated by the SW6|6 operator applied to the quartets A, B, and C (Figure 3a), the energy effect favors the rotation of the quartet A. Indeed, the calculated ΔE values for modes A, B, and C equal 147.9, 323.5, and 347.4 kJ/mol (see Appendix A for the details of the calculations). Note that the structures obtained via the rotations of quartets A and B have the CS symmetry whereas the structure derived by the quartet C rearrangement has no symmetry. The configuration of the molecule after the SW6|7 rotation of quartet A is given in Figure 3b, which evidences the new 5|7 double pair. At this point, the SW wave with η = 0 has been created on the surface of the new isomer C with f5 = 14, f6 = 106, and f7 = 7. The two extra pentagons and the two heptagons are labeled with the number of respective edges in Figure 3b. According to the propagation mechanism (Figure 2c,d), the second SW6|7 operator rotates the dual bond between the heptagon and the nearby hexagon (the arrowed dual bonds in Figure 2b,c), moving the 5|7 pair along the fullerene network. At each step η, the 5|7 pair swaps its place with a 6|6 one (see also Figure 3b); after η iterations of the SW6|7 rearrangement, the C isomer is formed. Figure 4 shows the isomer with η = 2. The described propagation process is applied to the isolated 5|7 dislocation monopole as well as to the 5|7 double pair arising from the SW6|6 rearrangement. It is worth mentioning that similar topological tools are used in other disciplines such as biology whereby wave-like diffusion processes and icosahedral patterns of viruses are normal mechanisms to model biochemical processes [39][40][41][42].
From the pure topological point of view, it is noteworthy that each fullerene network in the current study could be considered as the result of an instantaneous transformation generated by a single non-local generalized SW rotation. This kind of transformations, previously proposed in [25], constitutes a potentially infinite class of global rearrangements aiming to generate the entire isomeric space of a given Cn fullerene starting from a limited number (even only one) of inequivalent cages. fullerene. It exhibits the extended dislocation dipole generated by applying two SW6|7 rotation to the configuration shown in Figure 3b. The topological defect consists of two pairs of 5|7 and two pairs of 6|6 in between (the faces are numbered according to the number of carbon atoms in the ring). The molecule is ready for another SW6|7 rearrangement involving the lower pair 5|7 and two neighboring hexagons (marked with yellow diamond) to produce the C fullerene (not represented).
The final part of this section briefly introduces the topological invariants and how they measure the relative stability of fullerene isomers C modified by SW waves. Graph-theoretical results are presented and compared with energy values of the systems coming from quantum chemical simulations in the next sections. According to the topological modeling (TM) methods [43], chemical structures with n atoms and B bonds are schematically described as simple graphs G with n nodes and B edges. Geometrical, structural, and other physicochemical properties (maximal molecular symmetry being in the opinion of the authors the prominent, and somehow surprising, one) are embedded in the graph topology as the study of distance-based invariants evidenced some decades ago in pioneering studies [44,45]. In this context, distances dij are integer numbers computed by counting the number of edges connecting two atoms i and j with the shortest path in G. The distances are graph invariants (i.e., they do not depend on the labeling of the nodes) and generate a vast parade of distance-based graph descriptors. Graph diameter is defined as M = max{dij}.
TM techniques mainly make use of two invariants W (1) and (3) to measure graph compactness and roundness, respectively. Assuming dii = 0, the Wiener index W is derived as the semi-sum of all distances in G: counting the contributions arising from each node in the graph: Intuitively, Equation (2) admits the minimum ( ) and maximum ( ) values for central (peripheral) nodes in G. The quotient between the two extremals: called extreme topological roundness denotes the topological symmetry of the graph: the lower (3) is the more symmetric G is (with the general constraint ≥ 1). We assign the role of the topological potential of the system described by G to the Wiener index and extreme topological roundness. Indeed, both quantities obey the minimum principle. Therefore, the TM model assumes that similar carbon systems arrange their structures to minimize W and . Heuristically, this approach is confirmed by the fact that icosahedral fullerenes C60 and C240 minimize both indices in the respective isomeric sets. Recent works [46][47][48] on C66 and C84 molecules confirm that stable fullerenes combine high topological symmetry (low ) with high compactness (low W). This fast and elegant computational approach is applicable to the evolution of the SWw-defected fullerene molecules C . The electronic configurations of the atoms constituting the C240 isomers with defective structures in Figure 3; Figure 4 keep the sp 2 -hybridization state. On the other hand, surface reconstructions induced by SWw cause complex modifications of "the electronic charge density and this varies the bond lengths within and around" [49], with effects also on the local curvature of the molecular cage.
We have performed the preliminary density functional theory (DFT) computations to monitor the topology effect on the geometrical parameters and molecular properties of the defected C240 cages. The DFT computations have been performed with standard procedures; the details are collected in Appendix A. We have computationally studied SWw on the C240 surface up to η = 6 ( Figure 5). The SW6|6 operator implies replacing six original hexagons with the emergence of four "defected" polygons (two 5|7 pairs). From the point of structural chemistry, we approximate SWw as a movement of one "defected" pair from the other due to similar operations. Thus, each C with η = 0-3 contains 12 original pentagons, 2 heptagons and 2 pentagons, static and migrating, generated by the wave (the remaining faces being hexagons). In Figure 5, static and dynamic defects are colored differently for clarity. The C structure contains one migrating tetragon resulting from involving an original pentagon by SWw. In further structures, this tetragon converts to a new static pentagon. In other words, the wave shifts the position of the original pentagon in two operations. Note that the structures with η = 3 and η = 5 contain pentagon-pentagon fusions. These structural features are presented in Table 1.  Local curvature (k) is a key parameter we pay attention as it strongly correlates with the stability [1,2,50] and reactivity [51][52][53] of the fullerene species. The SWw transformations occur locally on the C240 surface. However, the calculated curvatures indicate the change in the shape throughout the whole surface (Equation (A2), Appendix A). To demonstrate it, we have chosen the curvature in the pentagon regions, both original and emerged upon SWw. As follows from Table 1, all pentagons have the same curvature in original C240 but SWw makes the curvatures more diversified. Accordingly, some of pentagons become more flattened as the others become more curved. Herewith, the structures with a tetragon ( C ) and fused pentagons ( C and C ) manifest the highest curvatures in the set. Sphericity (Ψ) is another parameter indicating the changes in the shape (Equation (A1), Appendix A). The calculated sphericities decrease from the starting value 0.9829 as the C240 system loses its initial icosahedral symmetry. Note that original C240 (Ih) and C have very close values because the structure with η = 0 still has a low symmetry (CS). The further changes are more pronounced. Nevertheless, Ψ does not become lower than 0.97 though SWw perturbs the shape. For comparison, Ψ may be significantly reduced (Ψ < 0.9) in highly functionalized fullerene compounds (e.g., C60 halogenides [54]). Sphericity (with associated quantity of volume, V) does not change monotonically with SWw propagation, i.e. there are no (and should not be) correlations between Ψ (or V) and η. However, we demonstrate below the correlations between these geometric parameters with topological ones.
Pristine and almost spherical C240 (Ih) has n = 240 direct nodes and B = 360, three symmetrydistinct sets of atoms (orbits) having multiplicities 60, 60, and 120, and 240 atoms with the same = = 2312. This last condition implies W = 277,400 and = 1, indicating that the C240 (Ih) fullerene is a maximal roundness molecule even in presence of non-symmetry equivalent atoms (a special condition featured by tetrahedral C40 too [49]). We would like to remember here that the complete individuation of all families of non-transitive fullerenes with = 1 is still an open task.

Topological Modeling
The first step of the SW wave on the C240 (Ih) fullerene graph consists of a SW6|6 bond flip that produces the 5|7 double pair (Figure 3b). This surface transformation decreases the lattice compactness to W = 277,122 derivable from the direct computation of the graph chemical distances dij and Equation (1). In our approximated model, this -0.11% variation of the Wiener index represents the gain in topological compactness allowed by the creation of the SW defect in the molecule C .
The subsequent SW6|7 rotation (η = 1) splits the 5|7 pairs by inserting in between one pair 6|6 according to the topological mechanism described in (Figure 2a,b). This decreases by another -0.14% the value of the topological descriptor W = 276,730 for newly formed C . This "gaining" behavior holds for each propagation steps in the interval η = [0, 3] of SWw over the fullerene surface, increasing the topological compactness of the system.
The reduction of the topological index W follows a quasi-parabolic trend ( Figure 6). The curve evidences that the creation of an extended topological defect up to η = 3 is favored; the topology of the fullerene network allows the topological diffusion of the 5|7 pair over a quite extended region. The SWw propagation is then impeded by a small +0.02% increment of the topological index W having a local peak for η = 4. The Wiener index curve in Figure 6 evidences, and this result is reported here for the first time for the fullerene molecules, that the 5|7 defects are able to migrate on the fullerene surface to form extended linear defects. This is a specific behavior that matches similar findings previously reported for graphene planar layers [26] and curved carbon nanotubes and nanotori [33]. The present TM simulations show that creation of SW waves is allowed also on quasispherical surfaces tiled by hexagons. Based on this outcome, we may therefore propose the following important conjecture: the formation and propagation of the SW waves are allowed on surfaces tiled with hexagons independently of the surface genus.
This effect has a pure topological root and strongly correlates with the connectivity properties of the pentagon-heptagon pairs embedded in the hexagonal mesh and the "edge" effect induced by the presence of pristine pentagons on the surface. We note in fact here that the small barrier effect encountered by the 5|7 pair at η = 4 ( Figure 6) is an effect generated by SWw colliding with one of the original pentagons of the C240 (Ih) fullerene. A different situation is encountered in the variations of the topological roundness (Equation (3)) of the system upon SWw diffusion, see Figure 6. The generation of the 5|7 double pairs reduces the topological symmetry and moves the curve upward after the first SW6|6 rotation, with a topological increase = 1.0241 for C compared with the unitary values of the original C240 (Ih). The successive SW6|7 flip starts the dipole propagation still increasing the topological invariant . The diffusion of the wave is opposed by the action of the invariant , which shows a barrier effect confirmed at the successive steps in Figure 5. At η = 2 and 3, the quantity show an almost flat curve in the same region where W shows a local minimum ( Figure 6). This result is quite interesting since it confirms that system, when the dislocation dipole has size η = 2 and 3, stops the obstacle to the creation of defective molecules (C represented in Figure 4). Topological invariants are listed in Table 2 and represented in Figure 6. We have found the inverse correlations between the extremal roundness (Table 2), volumes and sphericities of the C240 cages (Figure 7a). The found correlations are the same as in the case of the C84 isomeric series (made up with the isomers obeying the isolated-pentagon rule) [55] and similar to the findings of comprehensive work [56] on all possible isomers of C60. The geometry and topology of the molecules obtain direct linkage in the relations of Ψ and . Indeed, Ψ → 1-for sphere-like species whereas → 1+ in the case of structures with "equilibrated connectivity". These terms, in the case of the fullerenes, coincide with the notion that is manifested with the found correlations. As for the correlation between the volume and extremal roundness of the C240 cages, it follows from the interdependence of the Ψ and V values (Figure 7b).

DFT Simulations on Energetic and Molecular Parameters of the Defected C240 Cages
We have assessed the possibility of SWw transformations on the C240 cage from the thermodynamic point of view. For this purpose, we have operated with the DFT-computed total energies of the original and defected cages (see Appendix A for the details of the DFT study). As expected, the energies of all SWw-generated C240 cages are much higher as compared with the pristine one (Table 3). However, the energy effects of the stepwise transitions between C shows that some of the transitions may be favorable. Indeed, the transitions are endothermic up to C → C and exothermic for the two next steps ( Table 3). The exothermicity is explainable in terms of the structural chemistry of fullerenes. The structure C , a participant of the two exothermic steps, has two adjacent pentagons. This is the instability factor making the conversion C → C favorable. As for η = 4, the precursor of C , it has the tetragon-heptagon junction, and exothermicity C → C indicates that such an introduction of tetragon to the C240 cage makes the structure less favorable as compared with one containing the fused pentagons.
The HOMO-LUMO gap is another indicator of the fullerene stability [55] (HOMO and LUMO are the highest occupied and lowest unoccupied molecular orbitals, respectively). It is calculated as the difference in the HOMO and LUMO energies. In our isomeric C240 series, only the original structure manifests the largest gap value that indicates its stability. Other (defected) cages, as follows from the orbital energies, should be kinetically unstable. We have tried to collate the obtained HOMO-LUMO gaps with the results of the TM modeling and found a correlation between the gap and extreme topological roundness values (Figure 8). Indeed, both values, being different in their origin, should reflect the stability of the C240 cages. However, the correlation is weak (R 2 = 0.861) and requires further investigation, which includes a larger number of the defected fullerene structures.

Discussion
The results of the TM simulations are summarized in Table 2 and Figure 6. The analysis of the evolution of the topological compactness W of the isomers C ( Figure 6) allows observing how the generation η = 1 of the SW dipole causes an immediate decrease in the network invariant W = 277,122 (Equation (1)). In our topological model, this negative variation in the Wiener index represents the topological gain in the graph compactness in respect to the pristine C240 (Ih) fullerene. The creation of the SW defect is therefore able to make the fullerene cage more compact. The successive SW6|7 rotations further augment the topological compactness of the structure, the reduction in W showing an almost parabolic trend for the first propagation steps η = 1, 2, 3. This characteristic behavior of the W index, which favors the expansion of the SW wave, has a pure topological root and it relates to the connectivity properties of the hexagonal regions partially covering the fullerene surface. This result, in fact, matches with previously reported behaviors of SW waves when they propagate in pure honeycomb lattices [26]. The topological diffusion of the 5|7 pair with creating extended dislocation dipoles is stopped at η = 4 and W reaching a local maximum (W = 276,607). Here, the SW wave "collapses" with one of the pentagons of the pristine icosahedral fullerene with the creation of a 4ring of carbon atoms. Figure 6 shows instead that the topological order gets diminished by the creation and the propagation of the SW wave. The increase in the extreme topological roundness (Equation (3)) represents, in fact, an evident "topological response" of the cage to the diffusion of the pentagonheptagon pairs over the fullerene surface. The evolution of demonstrates a characteristic plateau (η = 3, 4 with ≅1.05 and η = 4, 5, 6 with ≅1.06) where the topological symmetry of the cages remains pretty uncanged, opening the door for the presence of extended dislocation dipoles embedded in the fullerene mesh.
In summary, the main outcomes of the TM simulations state that isomers C have greater compactness (lower W) than the parent C240 (Ih), favoring the creation and the propagation of the SWw defects over the fullerene surface. In these molecules, however, at each step of the SWw propagation, the index increases (lowering their extreme topological roundness) and obstructs the diffusion of the SWw itself; this barrier effect is reduced in the case of the large η values.
The DFT-estimated energy parameters, in general, indicate endothermicity of the transformations towards the defected fullerene isomers. Thus, such SW processes may take place under high-energy (nonequilibrium) conditions (plasma, interstellar environments etc. [20,21]). The computed HOMO-LUMO gaps also reveal the kinetic instability of the defected cages (in contrast to the parent one having the largest gap). We do not deeply discuss the HOMO-LUMO vs. plot ( Figure 8). There values seem correlated, but the justification of the trend requires extending the correlation field.
The title defect introduced to the honeycomb carbon-containing system is called a dislocation dipole. In the case of the fullerenes, we associate it with nonzero dipole moments (Table 3) of the defected fullerene cages (indeed, it is impossible to make such attribution in the case of introducing a dislocation dipole to the infinite systems, like graphene [57], as their dipole moment is not defined). The influence of the defects on the polarizability of fullerenes has been studied in a few works concluding the increase with the topological defects (see in review [58]). In further studies, we will investigate this important property in the context of SW waves.
As usually considered, the title defect also generates negative-curvature regions on the fullerene surface. In our estimates deduced from the pyramidality angles (Equation (A2)), we have found the atoms with decreased curvatures of the surroundings but which are still positive. Thus, the presence of the pair 5|7 is not the only condition of the negatively curved surface. Additionally, upon SWw, more curved regions also emerge, especially with tetragon inclusion or pentagon-pentagon fusions. Curvature is an important parameter for assessing the reactivity of fullerenes and their derivatives [1,[50][51][52][53]. We think that our estimates may be useful for assessing chemical properties of the giant fullerenes, hollow [4][5][6]50] or filled (as a part of carbon nano-onions [59,60]).
The calculated curvature values correspond to the local changes in the C240 structures upon SWw. Sphericity is another parameter reflecting the general change in the shape of the molecule. As expected, it goes down with introducing the defects. More important here, sphericity, a geometric structural parameter, is correlated with extreme roundness, a topological structural parameter. This indicated the deducibility of the structural and energetic parameters from the topology of the molecules.

Conclusions and Prospective
In the present work, topological, energetic and structural modifications induced on the C240 fullerene by the presence of a propagating expanding SWw have been investigated for the first time. Although cage compactness is enhanced for the dislocation-dipole isomers, their energies increase with the size η of the SW dislocation. This means a relatively low probability of the formation of these defective fullerenes. However, the existence of these defective isomers, especially those molecules with large η, which combine low W with non-increasing , may relate to special thermodynamic conditions like during fullerene growth or in interstellar media.
As found, the extreme topological roundness correlate with volumes and sphericities of the C240 cages. Such correlations are general for fullerene structures (as they have also been found for other isomeric series). Considering these regularities, we consider the efficiency of the presented topological approach for structural studies on fullerenes and deducing their molecular properties.
These results are perfectly in line with the topological determination of the relative stability of smaller fullerenes. The ability of W in selecting the most compact cage as the physically stable one is documented by the non-trivial case of the C60 fullerene [44] that, among its 1812 non-isomorphic isomers, minimizes the Wiener index in correspondence to the C60 (Ih) with W = 8340 and = 1. In more complex C66 case [46], the stability of the whole set of 4478 isomers has been ranked with maximizing two graph properties, i.e. topological compactness (low W) and topological roundness (low ). As concluded, the most stable C66 isomers minimize both graph invariants and, hence, may be found in the local minima of the (W, ) plane. We recently found that this purely topological approach works also for the C84 fullerene. Its most stable forms have minimal graph invariants W and [47,48]. The latter result also evidences the computational convenience of the TM simulations that rank just in a matter of minutes the topological stability of 51,592 distinct C84 isomers. Thus, we highlight that the TM approximated method is generally applicable to the formation and stabilization mechanisms of complex sp 2 carbon structures, like large fullerenes. As an additional advance, topological studies allow quick sieving the possible structures before detailed but time-consuming ab initio computations.

Acknowledgments:
The friendly support and suggestions of Franco Cataldo (Actinium Chemical Research Institute, Italy) and Mihai V. Putz (West University of Timisoara, Romania) are warmly recognized. The work was performed under the theme "Novel theoretical approaches and software for modeling complex chemical processes and compounds with tunable physicochemical properties" (registration number AAAA-A19-119022290011-6, Institute of Petrochemistry and Catalysis of RAS).

Conflicts of Interest:
The authors declare no conflict of interest.
We used the Cartesian coordinates of the atoms for calculating the volumes, sphericities, and curvatures of the fullerene cages. The term volume of the fullerene cage means the volume of the polyhedron made up with nuclei of the atoms of the cage (i.e. nuclear volume). The algorithm for the volume calculations has been previously presented (in general, it implies the triangulation of the fullerene cage, i.e. its partition over the disjoint pyramid primitives having one common vertex in the center of mass of the cage) [64]. The sphericities (Ψ) of the C240 cages are calculated as previously described [47,54]: where V and A are the volume and the surface area of the fullerene cage. Sphericity is the ratio of the surface area of a sphere, having the same volume as the given particle, to the surface area of the particle.
The local curvature of each site (atom) is calculated using the optimized geometry of a fullerene as [63] k = 2sin θP/a, where a is the average distance between the site and its three neighbors, and θP is the pyramidality angle of the site. The local curvature of polygons are the average values of the constituting atoms.