Effects of Hydrogen Passivation on Fullerene-Derived Si30C30 Clusters

Silicon carbide (SiC) is a technologically significant material. A recent report on the abundance of C60 fullerenes in interstellar space, along with the presence of SiC precursors, sparked interest in potentially similar SiC nanostructures. C60 fullerene was found experimentally and is an exceptionally stable form of carbon. As Si and C have similar valence electron properties, it has been long envisioned that Si and SiC could also form similar fullerene-type structures. In this paper, Si30C30 fullerene-derived clusters were studied from the first principle, starting from an ideal Si60 fullerene template with various arrangements of silicon and carbon atoms and then relaxing them without any symmetry constraint. Hydrogen passivation was considered as well to model the effect of ligands that may be available during chemical synthesis processes. We have found that, after passivation, the relative stability of different configurations of Si30C30 clusters changed compared to the unpassivated structures, while some structures collapsed. We have also noticed several Si–Si and Si–C double bonds in the unpassivated structures. Upon relaxation, some Si atoms lost their hydrogen, while carbon atoms capture those hydrogen atoms. Finally, we tested the endohedral doping of a transition metal, tungsten atom, to stabilize and magnetize Si30C30. With hydrogen passivation, the magnetic moment on W atom was enhanced. Overall, the effects of passivation on these fullerene structures are not very straightforward.


INTRODUCTION
Due to the similar valence electronic structures of Si and C, there are reports on silicon fullerene being similar to C 60 fullerene. However, it is well known that the hollow-cage-type clusters by Si are not energetically preferred configurations, as compact structures are usually preferred due to sp 3 bonding. In some earlier studies, the Si 60 clusters were simulated with the same symmetry as C 60 . The results showed significant distortions and that Si 60 was not nearly as stable as C 60 (Li et al., 2000;Menon, 2001;Sun et al., 2003;Wang et al., 2006;Bainglass et al., 2017). A recent study showed that, for Si 60 , a lower energy structure could be achieved by a triangular-shaped structure with three Si 20 units at the vertices (Bainglass et al., 2017) rather than by a hollow cage structure. Cage-like clusters are desirable because the hollow of the cage can be used for endohedral doping. Several endohedral dopants in the Si cages were considered to stabilize fullerene-like structures and to produce materials with tailored properties (Hiura et al., 2001;Kumar and Kawazoe, 2001). Some of these endohedral transition-metal-doped nanostructures can exhibit magnetic properties as well.
Besides that, it has been proposed that the spin property of the doped atoms inside the cage can be utilized in spintronics applications (Zhu and Khanna, 2012;Wang and Liu, 2016).
For similar reasons, as mentioned above, it is tempting to construct SiC fullerene structures. Silicon carbide is a technologically important material, for example, in power electronics (Gorai et al., 2019) and for device applications in extreme conditions (Cheung, 2006;Wijesundara and Azevedo, 2011). In addition, a recent astrophysical finding of the abundance of C 60 fullerenes in interstellar space also sparked interest in SiC clusters. SiC nanoclusters formed near dead stars were proposed as the precursor for the formation of these interstellar C 60 fullerenes under certain conditions (Candian, 2019). It is possible, for example, for some other intermediate "deformed fullerene" Si m C n (m + n = 60) structures to be formed before turning into C 60 or other carbon-based clusters.
Silicon carbide clusters were observed experimentally as well as studied theoretically (Flores and Largo-Cabrerizo, 1987;Nakajima et al., 1995;Fye and Jarrold, 1997;Pellarin et al., 1997Pellarin et al., , 1999Yamamoto and Asaoka, 2001;Ray and Huda, 2006;Li et al., 2011;Huda and Ray, 2012;Song et al., 2013;Anafcheh and Ghafouri, 2014;Huda, 2014;Yong et al., 2014). The theoretical studies by Srinivasan et al. (2005Srinivasan et al. ( , 2006 on Si 60 C 2n (n = 1, 2), Si 40 C 20 , Si 36 C 24 , and Si 60 C 20 showed that the SiC clusters were more stable than the bare Si 60 cage and further reported that the stability of the cluster depended on the carbon atoms' configurations inside or on the surface of the cages. It also reported that the Si-C interactions significantly affected the stability of the cages as compared to the C-C interactions. A similar conclusion was reached for these sizes of SiC clusters by  by noting, in an alternate manner, that the higher number of Si atoms away from the Si-C interface increases the instability of the structure. On the other hand, studies on Si 12 C 12 clusters showed that the stability of these structures was due to the maximizing strong C-C bonds in the segregated carbon region, minimizing the weaker Si-Si bonds in the segregated silicon regions and the Si-C bonds which are at the interface (Duan and Burggraf, 2015).
Smaller SiC cage-type structures were reported as well; for example, a highly stable SiC cluster consists of an eight-atom silicon cube with a carbon dimer inside it (Huda and Ray, 2004a). In another theoretical study, various stable isomers of Si 20 cages were reported, where carbon atoms were put inside these cages; it was reported that ionic interactions play an important role in the stabilization and the symmetrization of these Si 20 C n cages (Huda and Ray, 2004b). On the other hand, starting from a set of smaller C 20 fullerene structures, Si atoms were put systematically on the C 20 fullerene surfaces. These set of Si m C n (m + n = 20) clusters were then studied computationally (Huda and Ray, 2008). It was found that, even though SiC fullerene could not be formed starting from C-fullerene, the process provided some interesting thermodynamically stable structural motifs. In fact, both studies provided several useful metastable cage structures.
It is worth noting that it is customary to identify the ground-state structure for a given size of cluster with the first principle density functional theory (DFT) calculations. However, metastable structures can be stabilized during synthesis via kinetic processes, such as stabilization by surface energy or by surface passivations. For example, even though diamond is the ground-state structure for the bulk Si, ultra-small Si nanocrystals can be stabilized in wurtzite form, especially by hydrogen passivation on its surface (Mayfield and Huda, 2013). One such example for Si 30 C 30 clusters can be found in this work as well, where a less stable isomer becomes more stable due to passivation. On the other hand, the synthesis route can influence the favorability of one structure upon the other (Abreu et al., 2015). In addition, the relative thermodynamic stability of a structure obtained by DFT is at 0 K. In contrast, at a higher temperature, a highly stable structure can be dynamically unstable, and in comparison, a relatively less stable structure can become dynamically more stable (Scipioni et al., 2011).
In this paper, we have considered Si 30 C 30 clusters with various arrangements of silicon and carbon atoms, starting from a Si 60 fullerene structure with icosahedral symmetry and then relaxing those configurations without any symmetry constraint. These structures were relaxed without and with hydrogen passivation, where only Si atoms were passivated. The goal of this study is to understand the effect of hydrogen passivation on SiC fullerene's stability. Note that the aim was not to predict the ground-state structure of the similar-sized SiC clusters. Previous studies (Wang et al., 2006) showed that Si 60 H 60 is energetically stable with high symmetry due to the strong interaction between Si and H. On the other hand, it has also been demonstrated that Si 60 H 60 is less stable than Si 60 because the H atoms, while encouraging sp 3 bonds with Si, take away the sp 2 -bond-like feature from the Si 60 Si-Si bond network (Bainglass et al., 2017). In addition, the literature review above shows inconsistent reports regarding the relative importance of C-C, Si-C, and Si-Si bonds on the stability of the SiC cage-type structures. These questions will be addressed in this study, with emphasis on how passivation impacts these bonds. Finally, a representative transition metal, W, was incorporated in the most stable Si 30 C 30 cluster to study its impact on cluster stability along with hydrogen passivation. Brief comments will be made regarding the magnetic properties of the Si 30 C 30 structure within the scope of this study with the transition metal doped in it. For the transition metal, only the tungsten atom is considered. The rationale for selecting the W atom was that W has six valence electrons, and it has almost the same electronegativity as the Si atoms, with first ionization potential of about 8 eV. So, it is expected that the W atom will retain its spin moment inside the cage. In fact, this was also one of the first atoms that were encapsulated in Si cages experimentally (Hiura et al., 2001).

COMPUTATIONAL METHOD
The DFT (Kohn et al., 1996), with B3LYP hybrid functional (Becke, 1993) as implemented in the Gaussian 09 suite of codes (McLean and Chandler, 1980), was employed to study the Si 60 , Si 30 C 30 , and Si 30 C 30 H 30 fullerenes and fullerenederived structures. The atomic wave functions were populated by the atomic centered all-electron 6-311G basis set. Binding energy (E b ), highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), and several structural features were calculated. The binding energy (E b ) values were calculated with respect to the infinitely separated atomic limit, namely by: where E(si), E(C), and E(Si m C n ) and are total energies of the Si and C atoms and the Si m C n clusters in their respective ground-state spin multiplicities; m and n are the numbers of Si and C atoms in the Si m C n clusters, and for this study, m = n = 30. A negative number will imply a thermodynamically stable configuration; the more negative a binding energy is, the more stable a structure is. For hydrogen-passivated structures, the modified formulas will be discussed in a later section. The basis set and hybrid functional used here have been used in several earlier studies, including in one of our studies (Bainglass et al., 2017) (Chase et al., 1982;Park et al., 2006), respectively. These numbers show very good agreements between our calculated and experimental geometric features. Therefore, the relaxed structures in this study, both passivated and unpassivated, represent reliable configurations to study their structural and electronic properties. Finally, for W-doped Si 30 C 30, we used the LANL2DZ basis set (Chiodo et al., 2006) as available in Gaussian 09 package, keeping the B3LYP functional. Such a method for W is justified in previous publications (Sun et al., 2002). The binding energy was calculated in a similar approach as in Eq. (1): here l is the number of W atoms inside the Si 30 C 30 cluster; for the present study, we considered l = 1 and 2. Tungsten atomic energy, E(W), is in its lowest energy spin-state. For Si 30 C 30 W l, we have only considered the lowest energy structure of Si 30 C 30 obtained in this study and its passivated counterpart.

RESULTS AND DISCUSSION
Si 30 C 30 Fullerene

Structural Analysis
The well-known C 60 fullerene cage has 20 hexagons and 12 pentagons with icosahedral symmetry. Figure 1 shows a relaxed Si 60 structure with the same symmetry. This is the only structure in this study that was relaxed with symmetry constraints. As a result, the final structure remained a fullerene structure; in FIGURE 1 | Si 60 fullerene relaxed with symmetry constraint. This structure was used as a template to construct the other Si 30 C 30 configurations.
general, it has been reported that a fully relaxed structure of Si 60 is usually a distorted one (Li et al., 2000;Menon, 2001;Wang et al., 2006). The average bond length of Si 60 in Figure 1 is 2.03 Å. The purpose of symmetry-constrained relaxation was that the relaxed fullerene Si 60 was used as a template to construct the Si 30 C 30 structures to be studied in this paper. As mentioned in the introduction, earlier silicon carbide clusters based on Si 60 fullerene resulted in highly distorted structures. Among those, some SiC structures maintained cagelike features. For the present study, to build Si 30 C 30 , the first step was to take the symmetry-restricted relaxed Si 60 fullerene, and then 30 carbon atoms replaced 30 Si atoms at different possible configurations. Due to differences in the electronic properties of Si-Si and C-C bonds, SiC fullerene-like structures do not have high symmetry (Srinivasan et al., 2006;Huda et al., 2009). The relaxed Si 30 -C 30 configurations that have been investigated in this study are shown in Figure 2. Note that all these Si 30 C 30 structures, although started from an ideal Si 60 fullerene structure, relaxed completely without any symmetry constraint. Table 1 shows the binding energies, HOMO-LUMO gaps, and bond lengths for these Si 30 C 30 structures.
As the initial geometry was a fullerene structure, in the structure the Si atoms will have unsaturated 3-fold coordinated bonds, where the Si atoms are supposed to make π-bonds and the bond lengths are too large to Si-p orbital's side-wise overlap. To passivate these unsaturated bonds, hydrogen atoms were added to each Si atom in the unrelaxed Si 30 C 30 fullerenes initially and then optimized completely (Si 30 C 30 H 30 ); the initial set of structures was the same as for those in Figure 2, but here with passivating H atoms. The relaxed set of passivated structures is shown in Figure 3. After relaxation, in some structures, a few Si atoms lost passivating H atoms, and these H atoms then re-boned with nearby C atoms. Hence, these H atoms also act as an indicator to investigate the local bonding coordination and configurations. The binding energies, the HOMO-LUMO gaps, and the bond lengths for Si 30 C 30 H 30 are given in Table 2.

Nomenclature
In the following discussion, each structure in Figure 2 will be assigned a name to facilitate further discussions. In Figure 3, the initial structures were the same as in Figure 2, but with added H passivations, so the names in Figure 3 have an H prefix with each TABLE 1 | The following table shows the binding energies (E b ), highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gaps, the shortest Si-Si, Si-C bond lengths, and the average bond lengths of Si-C, C-C, and Si-Si bonds for different configurations of Si 30 C 30 structures. of them. The structures in Figures 2, 3 are arranged in order of their relative stability. In the remaining discussion of the paper, these names will be referred to when discussing these structures. The first one in Figure 2A is the C 30 -bowl structure: the left side of this structure has 30 carbon atoms making a bowllike configuration, and 30 Si atoms are on the right of it. The initial Si 60 fullerene was segregated with 30 C and 30 Si atoms and then relaxed entirely without any symmetry constraints. The 30 carbon atoms relaxed into almost half of the C 60 fullerene sphere. The C 30 -bowl structure has Si-C bonds only at the two segregated regions' boundary. There are 10 Si-C bonds with an average of 1.88 Å.
In Figure 2B, the structure C 15 -Si 30 -C 15 is presented, where 15 C atoms make caps in the left and the right ends of the fullerene cage and the 30 silicon atoms are in the middle. In Figure 2C, a similar model was constructed, but this time the 15 silicon atoms are located in the right and the left end, and the 30 atoms of carbons are in the middle of fullerene, Si 15 -C 30 -Si 15 . However, because of the shape of the relaxed structure, we named it as a tubular-Si 30 C 30 structure. Both structures have 18 Si-C bonds with the same average bond lengths of 1.86 Å.
The next structure in Figure 2D is named as 6-penta, where C atoms substituted six pentagons in Si 60 , and the rest were Si atoms. In this structure, there are Si and C pentagons, and the others are Si-C hexagons. There are 18 Si-C bonds with an average bond length of 1.84 Å.
In Figure 2E, the structure is named 4-Hexa-C, where we select the four hexagon faces in the diagonally opposite sides of Si 60 fullerene sphere (front and back in the viewing direction), replaced by carbon atoms, and three carbon atoms at the top and bottom in the pentagon faces of that fullerene. In Figure 2F, a complementary structure of Figure 2E was constructed, where the positions of Si and C were swapped. We call it 4-Hexa-Si. The final relaxed structures of these two configurations are significantly different, as seen in Figures 2E,F. However, the total Si-C bonds in these two structures are the same, 30, with equal average bond lengths of 1.89 Å.
In Figures 2G,H, we model the structures to minimize the number of C-C and Si-Si bonds and to maximize the number of Si-C bonds. In the structure, as shown in Figure 2G, the Si and the C atoms were placed alternately on Si 60 fullerene as much as possible, however, there remained some Si-Si and C-C bonds in the structure. The relaxed structure was somewhat a spherical one; therefore, we named it as sphere. The total number of Si-C bonds is 56 for this structure with an average bond length of 1.88 Å. The last structure that we tried in this group is shown in Figure 2H, where the structure starts with a pentagon face by silicon atoms (the frontal pentagon). The next neighbor atoms were replaced by carbon atoms and, after that, the next neighbors were kept as silicon atoms and went on alternating until the structure ends with the pentagon face of carbon atoms (the back pentagon). So, in this structure, two diagonally opposite pentagons are made up of Si and C atoms, respectively. This structure also preserves its spherical nature, and we named it sphere-2. The total number of Si-C bonds in this structure is 60, the highest in all the Si 30 C 30 configurations studied in this paper. The average Si-C bond length here is 1.86 Å.
The hydrogen-passivated structures are presented in Figure 3. They were constructed in the same way as described above, with added hydrogen on the Si atoms; hence, these passivated Si 30 C 30 H 30 structures are named with an H-prefix.

Electronic Properties
The relaxed C 30 -bowl structure in Figure 2A shows the distortion in the Si 30 side, as expected. The reason is the tendency for Si atoms to form sp 3 bonds with each other, as compared to the favorability of sp 2 in C-C bonds. On the C-side, the C 30 bowl remained almost a part of the C 60 fullerene, except at the boundary where the Si-C bonds formed. The average bond length of C-C is 1.39 Å, and the average Si-C bond length is 1.88 Å (presented in Table 1). The experimental value of the S = C double bond length was reported to be 1.77 Å (Bravo-Zhivotovskii et al., 2008), which is close to the shortest S-C bond length of 1.81 Å in this structure, as listed in Table 1. The bond length of this structure compares well with previously published values for a similar structure also calculated by DFT (Matsubara and Massobrio, 2005). The binding energy of this structure was −6.214 eV/atom, which is the most stable one in the unpassivated group of structures studied here. This structure is more stable than the previously studied similar Si and C segregated Si 40 C 20 and Si 36 C 24 fullerene-derived structures with binding energies of −4.748 and −4.977 eV/atom, respectively (Srinivasan et al., 2005). A trend that emerges from these three segregated silicon carbide structures is that, as C atoms are increasing in the carbon side, the clusters' stabilities become higher. All these three structures have dipole moments of about 1.4 D, so the orbital charge distribution is not spherical, indicating a polar nature of the structure. A similar conclusion was reached by , where they also found this Si 30 C 30 configuration as the one with the lowest energy. Figure 2B shows the relaxed C 15 -Si 30 -C 15 structure, where two groups of 15 carbon atoms are placed on the right and the left side of the fullerene. These two C 15 caps consist of two hexagons and two pentagons, and they almost became flat upon relaxation, not like the C 15 bowls. As a result, the 30 Si atoms are spread out almost in a bi-layer structure with two C 15 caps in the middle of each layer, with a bi-layer separation distance of 4.55 Å for C-C distance and, around the edge, the shortest Si-Si distance was 2.38 Å. The shortest Si-C bond length is 1.80 Å; the average of the C-C bond length in a layer is also 1.44 Å. This structure is the second TABLE 2 | The following table shows the binding energies (E b ), the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gaps, the shortest Si-Si, Si-C bond lengths, and average C-C, Si-C, and Si-Si bond lengths of the passivated structure Si 30 C 30 H 30 .
Type of passivated structure (Si 30 C 30 H 30 ) stable fullerene-derived structure with E b of −6.174 eV/atom and HOMO-LUMO gap of 1.171 eV.
Due to the strong C-C bonds and weak Si-Si bonds in the tubular-Si 30 C 30 , the 30 carbon atoms tightened the fullerene around the middle, and the silicon atoms that located on the caps extended a little bit without breaking the Si-Si bond, and the structure looks like a closed tube. The E b is −6.124 eV/atom, with −0.050 eV higher than C 15 Si 30 C 15 . However, the Si and the C atom distribution is more symmetric than the previous structures, as is evident by the lower dipole moment, with a reduction of about 50% and with a less polar structure.
The relaxed 6-penta is shown in Figure 2D, and it has a binding energy of −6.088 eV/atom. The average C-C bond length in this structure is 1.37 Å, the lowest among all Si 30 C 30 structures considered in Table 1, whereas the average Si-C bond length was also lower at 1.84 Å. This structure has a similarity with the previous structure of Figure 2C tubular-Si 30 C 30 , the main difference being that 6-penta has 30 C atoms as pentagons around the middle. Although lower in stability, 6-penta is higher in symmetry and with a very low dipole moment of 0.13 D.
The 4-Hexa-C structure shows the distorted open structure because some Si-Si bonds broke while relaxing the structure ( Figure 2E). The binding energy is 0.144 eV/atom higher than the most stable structure in Table 1, and it has a smaller value of the HOMO-LUMO gap of 0.998 eV. Its complementary structure is the 4-Hexa-Si structure shown in Figure 2F, where the Si atoms were switched with C atoms in 4-Hexa-C. The 4-Hexa-Si has a higher binding energy of −6.056 eV/atom; however, the binding energy difference between these two structures is only 0.014 eV/atom. It also has a higher HOMO-LUMO gap of 1.503 eV. It also does not show distorted open structure as was found in 4-Hexa-C structure with a smaller average Si-Si bond length of 2.27 Å, which is 0.12 Å lower than that of 4-Hexa-C. The average Si-C bond length for these two structures is the same, 1.89 Å; it is the largest Si-C bond length of all relaxed structures in Table 1. These two structures have zero dipole moments, implying that the effective negative and positive charge centers coincide at the center of these two cages. This also signifies an overall "spherical" nature of orbital charge distributions for these two structures.
The structure named sphere in Figure 2G is 12.060 eV higher in energy compared to that of the most stable C 30 -bowl, which agrees well with the total energy difference value between two similar structures obtained by Scipioni et al. (2011). In our calculation, this sphere structure has the largest HOMO-LUMO gaps of 1.605 eV. The last structure is shown in Figure 2H, sphere-2. It has the highest E b of −5.975 eV/atom, hence the least stable and with the lowest HOMO-LUMO gap of 0.972 eV, as presented in Table 1. These last two almost spherical Si 30 C 30 structures have a binding energy difference of 0.038 eV/atom. After relaxing these two configurations, we found that, even though the apparent spherical nature of fullerene was not significantly distorted like the other configurations in this study, the orbital charge distributions were significantly distorted. The dipole moments were 2.15 and 3.20 D, respectively, for sphere and sphere-2 structures.
All the structures discussed above have no higher spin state, that is, they have a spin multiplicity of one. The range of binding energy within the scope of this study is −6.089 ± 0.080 eV/atom, which implies that the fullerenederived Si 30 C 30 structures overall are very stable. The Si-C Frontiers in Materials | www.frontiersin.org bond lengths remained consistent across the structures, with an average value of 1.868 ± 0.017 Å. In addition, the binding energy trend follows the increasing number of C-C bonds in the structures. Figure 4 represents such a plot, which shows that the stability of the structures increases as the number of C-C bond increases in the Si 30 C 30 structures considered for this study. Except for the two least stable structures on the top-left corner of the plot, all the structures follow the trend. For these two structures, sphere and sphere-2, they have the least number of C-C bonds and the highest number of Si-C bonds; these are the most polar structures. A competition between the covalent and the ionic bonding changes the trend for these two structures.
So far, we have found that C 30 -bowl configuration, with segregated Si and C atoms, is the most favorable one for Si 30 C 30 cage structure within this study. We now ask the question on how the pre-relaxed input precursor structures impact the relative stability of this configuration. To investigate, we considered C 30 and Si 30 "half-fullerene" as two bowl structures separately, as shown in Figures 5A,B, respectively. The binding energy for C 30 is −7.669 eV/atom, and it has HOMO-LUMO gap of 1.118 eV; the Si 30 structure is relatively less stable with a higher binding energy of −3.722 eV/atom and with a slightly higher HOMO-LUMO gap of 1.139 eV, all of which are listed in Table 1. Then, we put these two Si 30 and C 30 half-fullerene at a distance of 5 Å apart, as shown in Figure 6. After full relaxation, these two "half-fullerenes" came together and made the bond between Si and C atoms at the periphery of these two bowl structures, with a final relaxed distorted Si 30 C 30 configuration. This configuration, named as Si 30 -C 30 -bowl, has a binding energy of −6.130 eV/atom and HOMO-LUMO gap of 1.046 eV. The main difference between the Si 30 -C 30 -bowl and the C 30 -bowl structures in Figure 2A is the initial structure: before relaxation, the initial structure for C 30 -bowl was constructed from a relaxed Si 60 fullerene. The C 30 -bowl has −0.084 eV/atom more stable than the Si 30 -C 30 -bowl structure. In Si 30 -C 30 -bowl, the carbon part of the structure is also more distorted than the C 30 -bowl. The two approaches for obtaining a similar segregated Si 30 C 30 configuration shows that (i) due to different initial conditions, even the stronger part of the structure, here the C 30 part, can be susceptible to distortions and (ii) initial conditions also determine the final stability of the clusters. Hence, in selecting the precursors while synthesizing these clusters, initial conditions need to be optimized.

The Hydrogen-Passivated Fullerene
Under passivation, the fullerene structures relaxed differently compared to their unpassivated counterparts. We note here again that the initial structures were the same as for the unpassivated ones, as presented in Figure 2; the difference is that, for the passivated initial fullerene configuration, each Si atom was bonded with an additional hydrogen atom. Upon relaxation, some of the passivated structures have shown interesting behaviors, such as some hydrogen atoms breaking bonds with the Si atoms and then re-bonding with nearby C atoms. This phenomenon resulted in significant structural changes as compared to their unpassivated counterparts. Hence, the relative stability order of the structures changed due to H-passivation, as can be found from comparing the binding energies in Tables 1, 2.
To estimate the stabilities of the passivated structures, we calculated the binding energies of Si 30 C 30 H 30 in the following way: First, we calculated the energy of the carbon-hydrogen bond [E(C−H) b ] by removing a hydrogen atom from the surface of a passivated structure that was bonded with a carbon atom and calculating the energy of the remaining structure without relaxing. Then, this energy and the energy of a hydrogen atom 1 2 E(H 2 ) was subtracted from the total energy of the   Figure 1A. However, the carbon side looks more distorted due to the different initial conditions. equations: where E(. . .) refers to the total energy of the system indicated in parenthesis; also b = 30 and d = 30; hence, b + d= 60. The results of E b and HOMO-LUMO gaps of the relaxed passivated structures are listed in Table 2. The H-4-Hexa-C structure in Figure 3A shows the most stable structure with hydrogen passivation, with E b of −6.439 eV/atom and with a significantly lower HOMO-LUMO gap of 0.375 eV. The higher stability is mainly because the hydrogen-passivated Si atoms have 4-fold coordination and created a favorable sp 3 -like bonding environment. Among all the hydrogen-passivated structures in this study, this is the one of the two structures where all the 30 Si atoms are 4-fold coordinated and the 30 C atoms are 3fold coordinated. Hence, the two C hexagons on each side of the structure have adopted a flatter configuration due to C-C sp 2 bonding, and hydrogenated silicon atoms occupy the curved edge sides with sp 3 bonds. Therefore, the initial Si 30 C 30 sphere became a squashed structure by carbon hexagons from both sides. The stability of H-4-Hexa-C increased significantly compared to its unpassivated counterpart 4-Hexa-C ( Figure 2E). Even though stability increased by −0.369 eV/atom due to hydrogen passivation, the dipole moment increased by 0.16 D compared to its unpassivated structures.
The second most stable passivated structure was found to be the H-C 30 -bowl structure (Figure 3B), with a binding energy of −6.249 eV/atom. The passivation helped to reduce the distortion on the silicon side as compared to the unpassivated C 30 -bowl structure. The relaxed H-C 30 -bowl has four H atoms migrated to make C-H bonds. These four hydrogen atoms made bonds with the neighboring carbon atoms, and the Si-C-H bonds made a sp 3 bonding environment. The value of the HOMO-LUMO gap of the H-C 30 H-C 30 -bowl structure is 1.108 eV. Passivation for this structure also increased the dipole moment by 0.97 D, with a total dipole moment of 2.64 D. The higher dipole moment in a cluster usually implies a relatively higher contribution of ionic properties in overall bonding.
The H-tubular-Si 30 C 30 structure ( Figure 2C) shows that all H atoms are still bonded with Si atoms after relaxation, and the final structure looks like a dumbbell-shaped structure, as the 30 carbon atoms in the middle slightly tighten the structure. The binding energy is −5.880 eV/atom, and the HOMO-LUMO gap is 1.328 eV, with a significantly short Si-Si average bond length of 2.29 Å. In the following discussion, we will show, by analyzing the HOMO plots, that some of these bonds may qualify for Si-Si double bonds. In this structure, hydrogen passivation reduces the dipole moment by 0.10 D when compared to its unpassivated counterpart. Figure 2D shows the H-C 15 -Si 30 -C 15 structure, where the C 15 units squeezed the structure from the right and the left sides of the fullerene. Some of the carbon atoms from these two C 15units has significant distortions after relaxation; as a result, 12 hydrogen atoms broke bonds with silicon atoms and created new bonds with unsaturated carbon atoms, with seven hydrogen atoms on one side and five on the other. Higher C-H bonds show here a higher order of distortions. The number of Si-C bonds in the unpassivated structure ( Figure 2B) was 18; after the passivation, the number of Si-C bonds reduced to 14. This shows significant reconfigurations within the structure after hydrogen passivation; as a result, the passivated structure is less stable by 0.979 eV/atom as compared to the unpassivated one (C 15 -Si 30 -C 15 ). The passivated structure has HOMO-LUMO gap of 1.537 eV, and the dipole moment of this structure is 3.81 D.
The H-6-penta structure in Figure 3E has an easily recognizable symmetry; there is a vertical mirror plane, and the two sides are then related by the mirror reflection and a 180 • rotation, which results in an inversion symmetry. Due to the inversion symmetry, this is the only hydrogen-passivated structure with a dipole moment of 0.00 D. The silicon atoms around the middle tightened the passivated relaxed structure. The mid-section Si atoms lose their hydrogen atoms to the carbon atoms, as those Si atoms are 4-fold coordinated. Each side of the H-6-penta structure has three pentagon faces of carbon atoms, which bonded with three hydrogen atoms. This structure loses four Si-C bonds upon passivation, as does the previous structure. However, here some Si atoms, which lost hydrogen, bonded with each other to make the midsection 4-fold coordination as seen in Figure 3E. This structure has binding energy of −5.632 eV/atom and HOMO-LUMO gap of 0.699 eV. Figure 3F presents the H-sphere2 structure, which distorted significantly after relaxation; its distortions are quite evident from the fact that 19 hydrogen atoms rebounded with carbon atoms after relaxation. Therefore, the bonding features have been significantly changed as well throughout the surface of the structure. In addition, it has the highest dipole moment among all the Si 30 C 30 structures considered here, about 6.00 D. Note that sphere2 (Figure 3H) has the highest dipole moment of 3.20 D in their respective unpassivated structure group as well. The H-sphere2 structure has binding energy of −5.631 eV/atom and HOMO-LUMO gap of 1.343 eV.
For H-sphere in Figure 3G, the stability decreased by 0.922 eV/atom, and the HOMO-LUMO gap increased by 0.192 eV compared to the most stable structure H-4-Hexa-C. No significant distortion was found in this structure after relaxing; in fact, it remained like a spherical cage. In addition, the average Si-C bond length, 1.74 Å, listed in Table 2, is smaller than the value of Si = C double bond length reported in earlier literature (Matsubara and Massobrio, 2005;Bravo-Zhivotovskii et al., 2008), which was 1.77 Å. Hence, the H-sphere structure has some Si = C double bonds in it. Even though 16 hydrogen atoms migrated to carbon atoms from Si atoms, the structure remains mostly spherical.
The least stable H-4-Hexa-Si structure is presented in Figure 3H; it is a distorted passivated structure. In this structure, 21 hydrogen atoms broke bonds with silicon atoms and rebonded with carbon atoms. There were 30 Si-C bonds in 4-Hexa-Si structure, whereas H-4-Hexa-Si has 23 of them after passivation. The unpassivated 4-Hexa-Si also has zero dipole moment with inversion symmetry, while with passivation, it gained a moment of 5.42 D. These are indications of the extent of distortions. The binding energy is −5.483 eV/atom and the HOMO-LUMO gap is 1.668 eV. Its stability is the lowest and the gap is the highest among all the Si 30 C 30 structures considered here.
A note regarding the migration of hydrogen atoms from the silicon atoms to the carbon atoms: The starting template of Si 30 C 30 was Si 60 fullerene, where all the Si atoms are 3fold coordinated, and all bond lengths, involving both Si and C, are equal to Si 60 fullerene's Si-Si bond length. Silicon atoms do not prefer 3-fold coordinated bonds among themselves. So, to provide their favorable sp 3 bonding environment, each Si atom in Si 30 C 30 was passivated with one hydrogen atom. On the other hand, upon relaxation, the C 30 part of Si 30 C 30 remains 3fold coordinated and creates almost a local environment of a half-fullerene of C 60 . Now, first, usually Si-Si and Si-C bonds are too large to allow for double bonds. However, due to the configurational effect in some structures, upon relaxation, some of the Si-C and Si-Si distances shorten enough to make double bonds. When a double bond is made involving a Si atom, that Si atom may not need the previously attached hydrogen atom. As a result, Si-H bonds became weaker during relaxation. Second, when a carbon atom near a Si atom favors sp 3 bonds due to a local distortion in structure, it gets the hydrogen from the Si atom as a C-H bond is much favorable than a Si-H bond. In this case, a carbon atom with a hydrogen atom became 4-fold coordinated. In both scenarios, hydrogen detached from these Si and then reattached to a nearby sp 3 -like C atom. In general, these types of hydrogen-passivated Si 30 C 30 structures are relatively more distorted from the ideal fullerene shape.

HOMO-LUMO Properties
In general, silicon prefers sp 3 hybridization, and forcing it to form sp 2 bonds leads to destabilization. The carbon can accommodate both sp 3 as well as sp 2 hybridization due to shorter C-C bond lengths (Huda and Ray, 2004a;Huda et al., 2009;De Padova et al., 2011). However, there were reports of sp 2like hybridization in silicon carbide, such as in SiC multi-wall nanotubes (Huda et al., 2009;De Padova et al., 2011). Due to local chemical configurations in Si 30 C 30 fullerenes, the molecular orbital calculations may establish the sp 2 -like π bonds that we hypothesized in earlier sections. The orbital plots can verify such features. Figures 7, 8 show the HOMO and LUMO plots for all the structures, without and with hydrogen passivation, respectively. In the following paragraphs, we discuss a few selected HOMO and LUMO plots from the figures.
The HOMO plot of the C 30 -bowl structure in Figure 7A shows the π-bond-like feature on the silicon side at the leftmost part of the figure where the Si-Si bond length is 2.31 Å. This bond length is comparable to the previously published Si-Si double bonds, which were in the range of 2.15 to 2.23 Å (Wang et al., 2008;Fischer and Power, 2010;Kira, 2012;Gau et al., 2016). On the other hand, for passivated H-C 30 -bowl structure in Figure 8B, there is no sp 2 hybridization on the silicon side, as all Si atoms, except four Si atoms, are passivated with hydrogen and have 4-fold sp 3 -like coordination, however, the π bond does show up on the carbon side of the structure. In Figure 8C, for the H-tubular-Si 30 C 30 structure, both the HOMO and the LUMO appear on the carbon atoms in the middle of the structure. Si atoms do contribute significantly to these frontier orbitals because the hydrogen-passivated silicon atoms created an almost perfect sp 3 -like bonding environment with bonding angles of about 104 • , whereas the ideal sp 3 bonding angles are 109.5 • . Figure 7B shows in the C 15 -Si 30 -C 15 structure that HOMO is almost symmetrically distributed on both C 15 units; except for a pair of carbon atoms on both faces, these two carbon  atoms are doubly bonded on their respected faces. Most of the Si atoms are not contributing to HOMO. Si atoms from one of the most distorted sides of the structure dominate the LUMO. Interestingly, H-C 15 -Si 30 -C 15 in Figure 8D shows the HOMO in only one of the C 15 caps, particularly on the side where more hydrogen atoms are broken away from Si atoms and re-bonded with C-atoms. On the other hand, LUMO, as it is for its unpassivated counterpart, is on the side of the structure that is more distorted with Si atoms without hydrogen atoms. The highest contribution for LUMO came from a pair of Si and C atoms that are separated by 3.66 Å.
Highest occupied molecular orbital and LUMO are mainly distributed on the Si caps in the 6-penta structure in Figure 7D. Unlike some other structures, no π-bond-like feature was found on the HOMO plot. H-6-penta in Figure 8E has main contributions to HOMO and LUMO, which are from the Si atoms that are 3-fold coordinated and without hydrogen atoms and from a few carbon atoms. Therefore, bondingantibonding splitting that created HOMO and LUMO occurs almost on the same set of atoms. The eight 4-fold coordinated Si atoms in the mid-section of the structure remained tightly bound, with no contribution to either HOMO or LUMO. H-sphere2 has the shortest Si-Si bond, 2.24 Å, which is significantly close to experimental Si-Si double bond length values, as quoted earlier (Wang et al., 2008;Fischer and Power, 2010;Kira, 2012;Gau et al., 2016). The HOMO plot in Figure 8F reflects the π bond between these two silicon atoms as well.
Regarding the double bonds, we make some remarks here: when a bond is well formed and saturated, the orbitals' energy for that bond becomes lower in energy, and it may not contribute to HOMO. In the unpassivated Si 30 C 30 structures, atoms that contribute to HOMO are from the part of the structure that is prone to distortion or are not still in their preferred bonding state. For example, in the C 30 -bowl structure, both HOMO and LUMO have major contributions from the unpassivated Si atoms. Once these Si atoms are hydrogenated, in the H-C 30bowl structure, the Si part of the structure does not contribute to these frontier orbitals anymore; these Si atoms tend to have sp 3 -like bonding with the help of passivating agents such as hydrogen. Therefore, some Si-Si double bonds showing up in HOMO in the unpassivated structures are prone to be reactive with chemical stimulations because, with passivation, these bonds are more likely oriented in standard sp 3 hybridization. For example, in both the most stable C 30 -bowl and lest stable sphere-2 structures, there were indications for Si-Si double bonds from the HOMO plots; then, with hydrogen passivation, these double bond features were lost. In general, in any chemical synthesis of these structures, there will possibly be passivating ligands available in the synthesis environment, and hence the detection of Si-Si double bonds in these larger structures is unlikely.

W Doping in Si 30 C 30
For the endohedral doping of tungsten atoms, the C 30 -bowl configuration is considered. The reason for selecting the C 30bowl is that, firstly, without passivation, it was the most stable Si 30 C 30 configuration, and it kept a cage-type feature; secondly, its passivation increased its stability (by 0.035 eV/atom) and improved its overall cage feature. Figure 9A shows the relaxed W-C 30 -bowl structure. The initial structure was constructed by placing W atom at the C 30 -bowl sphere's center, and then the structure was relaxed without any symmetry constraint. Figure 9B shows a similarly, obtained structure; this is a W-dimer inside the cage, W 2 -C 30 -bowl. The W 2 -C 30 -bowl was found to be magnetically not very interesting, as the W-dimer in the cage prefers not to retain spin as much magnetic moments as in W-C 30 -bowl-based structures. So, the following discussion is based on a single W atom in Si 30 C 30 , W-C 30 -bowl structure ( Figure 9A).
Several spin-multiplicities have been considered for the W-C 30 -bowl structure. From the binding energy calculations, the spin-triplet state was found to be 28 meV/atom more stable than the singlet state. In all spin-states considered, including the FIGURE 9 | The figures show the C 30 -bowl structure (as in Figure 2A) with tungsten atom(s) in the center of the structures: (A) with one W atom, C 30 -bowl-W and (B) with two W atoms, W 2 -C 30 -bowl. quintet state, after relaxation, the W atom moved toward the Si side. The C 30 side remained almost unaffected with fullerenelike structure. In the triplet state, there is only one W-Si bond, with a bond length of 2.57 Å. On the other hand, in the singlet state, the W atom is 3-fold coordinated, with an average bond length of 2.53 Å. For the triplet spin state structure, the W atom does possess almost all the spin, namely, 1.85 m B , out of a total magnetic moment of 2 m B . In contrast, in the hydrogenpassivated H-W-C 30 -bowl triplet structure, where hydrogen atoms passivate Si atoms, the spin-magnetic moment on the W atom is 1.99 m B , and the W atom remains almost at the center of the cage without any bond. A similar effect was found for the structure with an even higher spin-multiplicity of five states; 3.99 m B was found on the W atom out of the total 4 m B and sitting at the center of the H-W-C 30 -bowl structure. Retention of spinmagnetic moment on the endohedral transition metal atom is much higher in Si 30 C 30 cages than what was reported for smaller Si cages earlier (Guo et al., 2008).
In addition, the W atom increased the binding of the C 30bowl structure by 84 meV/atom and decreased the dipole moment by 59%; this implies that the stability of the overall structure improved by W atom doping. On the other hand, for the passivated H-C 30 -bowl structure, the HOMO-LUMO gap increased by 0.715 eV, implying increased electronic stability compared to the undoped structure. In short, two important findings with endohedral W atom doping in conjunction with hydrogen passivations are: (i) higher spin moment is retained at the W atom sitting at the center of the cage and (ii) the cage stability increases.

CONCLUSION
Starting from Si 60 fullerene template, several configurations of Si 30 C 30 clusters were studied first without any passivating agent. Then, all 3-fold coordinated 30 Si atoms were passivated by hydrogen atoms, and the effects of the passivation of the structures have been analyzed. The first noticeable effect of passivation is the change of the relative stability of the Si 30 C 30 isomers; the order of relative stability of the unpassivated Si 30 C 30 isomers changes when they are passivated. Second, the position of HOMO shifted significantly with H-passivation. Third, depending on the local chemical environment, such as coordination order, neighboring elements, and bond angles, some C-H bonds were found after relaxation, even though no carbon atom was bonded with H before relaxing the structures. Fourth, hydrogen passivation has significant effects on the endohedral metal dopant, such as it helps keep the metal atom almost at the center of the Si 30 C 30 cage and retain a very high spin-magnetic moment on the dopant metal atom. So, inside the hydrogenated Si 30 C 30 cage, the spin state of the W atom can be preserved from the outside chemical perturbations; this can have a potential application in spintronics applications.
Finally, enhancing the stability of SiC fullerene structures by surface passivation is not a very straightforward phenomenon. The integrity of the structure and the cohesion between the atoms within the structure upon passivation depend on the local chemical environment on the fullerene surface. Even though all Si atoms in both Si 60 and Si 30 C 30 are 3-fold coordinated, none of the hydrogenated Si 30 C 30 structures could maintain their fullerene sphere as was in Si 60 H 60 , albeit among the Si 30 C 30 structures considered here, a few isomers could maintain sufficiently large cage-type structures. Note that in the case of cage structures, there are no bulk atoms, and all are surface atoms. When bonding mismatch weakens the bonds on the surface atoms, even during passivations, distortions and even collapse of the cages are possible, as there is no support from inside; several of the Si 30 C 30 structures showed the behavior in this study. In addition, having a prior knowledge for the need of passivation for a particular set of atoms for a large cage structure can be a daunting task. In this study, we have found that, in several cases, some hydrogen atoms that initially bonded with 3-fold Si atoms were later moved to carbon atoms due to structural relaxations. The migration of hydrogen atoms from silicon to carbon atoms causes significant distortions in the cage structures; the higher the number of hydrogen atoms moved to the carbon atoms, the higher distortions were found in the Si 30 C 30 cages. The near-spherical one is the one where the sp 3 -type bonding was distributed almost equally over the structure, including the carbon part. For example, in the H-sphere structure in Figure 3G, even though it is one of the least stable isomers, it maintains a near-spherical form. A symmetrization of the bonding over the surface of the structure is needed to maintain a large fullerene-like spherical structure.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/supplementary material.

AUTHOR CONTRIBUTIONS
HA performed the computations. HA and NA performed the analysis. HA and MH mainly wrote the manuscript. NA helped in organizing the references. MH planned the project and supervised the work. All authors contributed to the article and approved the submitted version.