Structure Prediction and Mechanical Properties of Silicon Hexaboride on Ab Initio Level

Silicon borides represent very appealing industrial materials for research owing to their remarkable features, and, together with other boride and carbide-based materials, have very wide applications. Various Si–B phases have been investigated in the past, however a limited number of studies have been done on the pristine SiB6 compound. Structure prediction using a data mining ab initio approach has been performed in pure silicon hexaboride. Several novel structures, for which there are no previous experimental or theoretical data, have been discovered. Each of the structure candidates were locally optimized on the DFT level, employing the LDA-PZ and the GGA-PBE functional. Mechanical and elastic properties for each of the predicted and experimentally observed modifications have been investigated in great detail. In particular, the ductility/brittleness relationship, the character of the bonding, Young’s modulus E, bulk modulus B, and shear modulus K, including anisotropy, have been calculated and analyzed.


Introduction
Silicon borides are lightweight ceramic compounds formed between silicon and boron and are regarded as the most elusive refractory compounds [1]. Metal borides, especially silicon borides, are very appealing industrial materials with potentially very wide applications. Due to silicon boride's extraordinary features, it represents a very promising material for future research, and has excellent electrical conductivity [2], a high degree of hardness [2,3], a moderate melting point (2123 K), and a low specific gravity [2,3]. According to the literature, there are a few previously reported phases in the silicon boron system: SiB 3 , SiB 6 , and SiB n [4][5][6]. Silicon hexaboride was first reported in 1900 by Henry Moissan and Alfred Stock after briefly heating silicon and boron in a clay vessel, and was described as black, irregular crystals [7]. As it is known, one of the first reports of this structure described it as a cubic structure Pm-3m (no. 221) [8], but later confirmed it to be mechanically unstable [9]. Later studies presented SiB 6 as an orthorhombic structure [1,10], combined with interconnected icosahedra, icosihexahedra, and several isolated boron and silicon atoms [10]. Based on first-principle calculations and CALYPSO structure prediction methodology, two structure candidates, defined as mechanically stable, have been proposed: monoclinic P21/m and hexagonal R3-m structure types [9]. Furthermore, the first of these newly proposed structures, P21/m-SiB 6 , is found to be an indirect band gap of 0.41 eV, while R3m-SiB 6 -81, has a direct band gap of 1.654 eV. There is also an Due to its features, the silicon-boron system has been quite well-researched to examine its potential as a high-temperature material. Hence, it was suggested that SiB 6 exhibits very good thermoelectric material characteristics at high temperatures [2,24]. There have also been novel proposed synthesis methods for these surface-protected, oxidation-resistant semiconductor materials within the Si-B system, and these materials are considered to be very useful for various high temperature solar thermal, or solar electric applications [25]. Nowadays, studies are also focused on investigating the properties of SiB 6 compounds for their use in nanotechnology applications [25].

Computational Methods
Our general approach to performing structure prediction and gaining insight into the structural stability of possible phases existing in the SiB 6 system is based on a data mining search. We performed data-mining-based explorations of the ICSD database [26,27] via resemblance to known crystallographic structures. According to some earlier studies, a data mining search could be also used as an additional method for some others, such as global optimization, for instance, because it has proven successful in finding additional possible modifications in some chemical systems. Additionally, many relevant structure candidates in a given chemical system exhibited a very similar crystal structure to some other compounds observed in another chemical system, even though there was no obvious chemical relationship between these two chemical systems [28]. We have used the well-known knowledge discovery in databases (KDD) process, which involves selection, preprocessing, transformation, and interpretation/evaluation (or post-processing), and it has been already used successfully in some previous studies [29][30][31]. All potential structure candidates that appear in the ICSD database have been extracted and subsequently submitted to local optimization at an ab initio level. Details about the KDD process and the combination of data mining with ab initio methods can be found elsewhere [32,33]. After the structure, candidates were identified using the data mining approach, and are submitted to density functional theory (DFT) calculations. Local optimizations (including the cell parameters and atom positions) were performed using the CRYSTAL17 code [34][35][36], which is based on linear combinations of atomic orbitals (LCAO). Structure optimizations were performed on the DFT level, employing two different functionals: the Generalized Gradient Approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) functional [37] and the Local Density Approximation (LDA) with the Perdew-Zunger (PZ) correlation functional [38] for comparison. An all-electron basis set based on Gaussian-type orbitals was employed; in particular, in the case of Silicon a [5s4p1d] basis set was used [39,40], and the [3s2p1d] basis set was used in the case of boron [41][42][43] (further basis sets information are given in the Appendix B). In each structural optimization, Fock/KS matrix mixing was set to 30%, and the tolerances for the convergence on energy were set to 10 −7 . K-point meshes of 8 × 8 × 8 Monkhorst-Pack scheme have been used. Chosen LCAO/Gaussian basis set type approach and DFT methods have been shown highly efficient and precise in our previous theoretical studies [44][45][46] and in comparison to the experimental data [30,47,48]. A computational strategy implemented in the CRYSTAL17 solid-state, quantum-chemical program has been performed for the accurate ab initio simulation of elastic and mechanical properties of crystalline materials [49]. Full elastic tensor has been generated using the keyword ELASTCON [50]. Elastic tensor analysis and visualization have been performed using ELATE code [51]. The symmetries of the analyzed structures were determined using the SFND [52] and RGS [53] algorithms implemented in the program KPLOT [54]. The structures were visualized using the Vesta3 program [55].

Structure Prediction of Silicon Hexaboride
The data mining searches were performed within the ICSD database, which included more than 250,000 crystal structures in the latest release, with more than 180,000 structures assigned to 9873 distinct structure types [26,27]. In order to find new structures in the SiB 6 , data mining-based searches to find all possible A 6 X structure types in the ICSD database [20,21] have been used. In particular, the data mining has resulted in the previously investigated silicon hexaboride structures: the c-SiB 6 (or CaB 6 ) type, the SiB 6 (Cmca) type, the SiB 6 (P21/m) type, the SiB 6 -81 (R3m) type, the SiB 6 (Pnnm) type, and the following AX 6 structure types, which together could be used in search of another chemical system: the PB 6 (or α-B 6 O) type, the β-B 6 O type the BaN 6 type, the BaSi 6 type, the HgN 6 type, the OsOF 5 type, the Al 6 Mn type, the Ga 6 Pu type, the MoCl 6 (P-3m1) type, the MoCl 6 (P-3c1) type, the Au 6 Sm type, the MnU 6 type, the PbN 6 (Pna21) type, the PbN 6 (Pcmn) type, the RbTe 6 type, the c-SF 6 type, the LT-SF 6 type, the SF 6 (C-1) type, the MoZn 6 type, the SrN 6 type, the TeOH 6 type, the CeCu 6 type, the WCl 6 (R-3H) type, the WCl 6 (P-3m1) type, the XeF 6 (P121/c1) type, the XeF 6 (C12/c1) type, the XeF 6 (Pc21n) type, the PrAu 6 type, the Cu 6 Nd type, and the UCl 6 type.
Nine additional structure candidates have been created using the Primitive Cell approach for Atom Exchange (PCAE) method [31,56], despite resulting in non-stoichiometric compounds. Since the above-mentioned prototypes are not common, data mining and PCAE based searches resulted in 44 structure candidates. A total number of structure candidates were finally reduced after performing full structural optimization at the ab initio level, and four final structures were distinguished as being the most energetically favorable options, regardless of the computational approach (Table 1). It is not surprising to find a large quantity of the energetically non-favorable structures, since many of the prototypes from the data mining-based search resulted in several non-equilibrium structures in the SiB 6 and other parent compounds, or in non-stoichiometric compounds, as previously observed [9,10,57,58]. Table 1. The total energy values (in E h ) and relative energies (compared to the global minimum α-SiB 6 structure taken as the zero of energy in E h ) of the SiB 6 modifications obtained from data-mining-based searches and local optimization on the GGA-PBE and LDA-PZ level of calculations.

Total Energy
Relative Energy The most relevant structure candidates predicted for the SiB 6 compound are α-SiB 6, β-SiB 6 , γ-SiB 6 , and δ-SiB 6 modifications (Table 1). Full structural data of all favorable candidates are given in Table 2 for calculations with the PBE functional, and those computed with the LDA functional are shown in Table A1. Table 3 presents structure details and unit cell parameters of the four predicted candidates chosen for further theoretical analysis in comparison to the previous experimental and theoretical results where available. Besides the experimentally known cubic SiB 6 , the orthorhombic phase in space group Cmce (no. 64) has been recently predicted [17], and the results obtained within this study are in good agreement with those for both known structures on the GGA-PBE and LDA-PZ levels of calculation. We noted that, in both previous theoretical calculations involving SiB6 structures DFT (GGA-PBE) and CASTEP, code has been used based on the robust methods of a plane-wave (PW) basis set and pseudopotentials (PPs) (Table 3 [9,17]). Our LCAO-GGA-PBE calculations concur with these PW/PPs-GGA-PBE data, as we expect to agree with possible PAW/PPs-GGA-PBE and FP/APW + LO-PBE calculations [59]. Moreover, we have predicted two additional structures denoted as αand δ-SiB 6 modifications, for which there are no previous data, but the results of structural relaxation agree between two levels of calculations (GGA-PBE and LDA-PZ). Table 2. Full structural details (modifications, space group, unit cell parameters, and atomic positions) for predicted SiB 6 modifications obtained from data-mining-based searches and local optimization on the GGA-PBE level. The data mining ab initio method resulted in four final structure candidates in the SiB 6 system. The energetically most favorable modification after local optimization is denoted as α-SiB 6 -type and appears in hexagonal space group R-3mH (no. 166) with unit cell parameters of a = 6.164 Å and c = 12.079 Å on the GGA-PBE level of calculation. The α-SiB 6 phase is visualized in Figure 1, while full structural data are presented in Table 2 for calculation with the PBE functional, and in Table A1 when computed with the LDA functional, respectively.

Modification and Structure Type Space Group Cell Parameters Position of Atoms
From the literature data, it appears that B12 icosahedra, as the single primary unit of α-rhombohedral boron, is a fundamental structural element for most of the Si-B compounds [9,60]. This is the case in the α-SiB 6 phase, with the PB 6 structure type, as in the case of the α-B 6 O compound [46,58,61,62]. Boron atoms form B12 icosahedra with interatomic distances from 1.73 Å to 2.95 Å, as shown in Figure 1a. Moreover, silicon atoms in the second coordination polyhedra (CP) form corner-connected tetrahedra with Si-Si distances 1 × 2.448 Å and 3 × 3.893 Å (Figure 1b). On the other hand, when analyzing Si-B distances, the silicon atom is surrounded by only three boron atoms in the first CP (3 × 2.005 Å), and with 30 boron atoms in the second CP, creating a complex polyhedra, shown in Figure 1c.
A second, energetically favorable structure candidate, found through the data mining approach, is referred to as a β-SiB 6 -type of modification. It crystallizes in space group Cmce (no. 64) with unit cell parameters a = 5.894, b = 11.184, and c = 8.420 Å (GGA), and it has five non-equivalent Si1, B1, B2, B3, and B4 atoms in the structure, for which full structural data are given in Table 2. This orthorhombic structure is consisted of three layers of boron icosahedra within the unit cell, with layers of silicon atoms between them. As in the previous α-SiB 6 -phase, boron atoms within this β-modification form four distinct B12 icosahedra with atom-atom distances from 1.73 Å to 2.87 Å, and the structure is  Figure 2a. Besides, in the second coordination polyhedra (CP), silicon atoms are in four-fold coordination with interatomic distances of 2 × 3.638 Å and 2 × 3.977 Å, respectively (Figure 2b). Upon examination of Si-B bonding, it appears that, within the first coordination polyhedra (CP), the silicon atom is surrounded by four boron atoms, thus forming a tetrahedra with interatomic distances of 1 × 2.021 Å-B, 1 × 2.064 Å-B, and 2 × 2.085 Å-B, which is visualized in Figure 2c.
From the literature data, it appears that B12 icosahedra, as the single primary unit of α-rhombohedral boron, is a fundamental structural element for most of the Si-B compounds [9,60]. This is the case in the α-SiB6 phase, with the PB6 structure type, as in the case of the α-B6O compound [46,58,61,62]. Boron atoms form B12 icosahedra with interatomic distances from 1.73 Å to 2.95 Å , as shown in Figure 1a. Moreover, silicon atoms in the second coordination polyhedra (CP) form corner-connected tetrahedra with Si-Si distances 1 × 2.448 Å and 3 × 3.893 Å (Figure 1b). On the other hand, when analyzing Si-B distances, the silicon atom is surrounded by only three boron atoms in the first CP (3 × 2.005 Å ), and with 30 boron atoms in the second CP, creating a complex polyhedra, shown in Figure 1c. A second, energetically favorable structure candidate, found through the data mining approach, is referred to as a β-SiB6-type of modification. It crystallizes in space group Cmce (no. 64) with unit cell parameters a = 5.894, b = 11.184, and c = 8.420 Å (GGA), and it has five non-equivalent Si1, B1, B2, B3, and B4 atoms in the structure, for which full structural data are given in Table 2. This orthorhombic structure is consisted of three layers of boron icosahedra within the unit cell, with layers of silicon atoms between them. As in the previous α-SiB6-phase, boron atoms within this β-modification form four distinct B12 icosahedra with atom-atom distances from 1.73 Å to 2.87 Å , and the structure is visualized in Figure 2a. Besides, in the second coordination polyhedra (CP), silicon atoms are in fourfold coordination with interatomic distances of 2 × 3.638 Å and 2 × 3.977 Å , respectively ( Figure 2b). Upon examination of Si-B bonding, it appears that, within the first coordination polyhedra (CP), the silicon atom is surrounded by four boron atoms, thus forming a tetrahedra with interatomic distances of 1 × 2.021 Å -B, 1 × 2.064 Å -B, and 2 × 2.085 Å -B, which is visualized in Figure 2c. Maybe the most investigated structure in the silicon hexaboride, both experimentally and theoretically, is found as a cubic phase, denoted as γ-SiB 6 -type. It appears in space group Pm-3m (no. 221), while full structural data are given in Tables 1, 2 and A1. Previous reports show unit cell parameters of a = 4.130 Å, which concur very well with our GGA (a = 4.161 Å) and LDA (a = 4.109 Å) results (Table 1). First reports of the SiB 6 compound described it as a cubic phase and indicated that it is isomorphous with CaB 6 [8]. Hence, within this γ-modification, boron atoms have five-fold coordination, with interatomic distances of 1 × 1.663 Å and 4 × 1.766 Å, and the structure is visualized in Figure 3a. Additionally, in the second coordination polyhedra (CP), silicon atoms form octahedra with atom-atom distances of 6 × 4.161 Å, which is visualized in Figure 3b. Furthermore, analyzing Si-B distances of silicon atoms were found to be surrounded by 24 boron atoms, with the distance between the atoms being 24 × 3.057 Å-B, as shown in Figure 3c. In addition, it has been confirmed that there is ionic bonding in the isomorphous CaB 6 phase between the boron group and the Ca atom [63].
The last modification of four yielded from our searches was the rhombohedral structure, denoted as δ-SiB 6 -type, that crystallizes in space group P3m1 (no. 156) with unit cell parameters a = 3.503 and c = 6.407 Å (GGA). Besides silicon, there are two different boron atoms, B1 and B2, and full structural data are listed in Table 2 (LDA in Table A1). Within δ-SiB 6 -type of the structure represented with hexagonal axes, boron atoms are in six-fold coordination and form two distinct octahedra with atom-atom distances (B1 2 × 1.736 Å-B, 2 × 1.766 Å-B, 2 × 1.952 Å-B, B2 2 × 1.749 Å-B, 2 × 1.754 Å-B, 2 × 1.952 Å-B). Within δ-modification, six-fold coordinated silicon atoms with interatomic distances of 6 × 3.503 Å form a layered-like structure, with boron octahedra situated between these two layers 7 of 17 ( Figure 4b). In addition, when examining Si-B distances, silicon atoms were surrounded by six boron atoms, with a distance between atoms of 6 × 2.479 Å-B, which is visualized in Figure 4c. Maybe the most investigated structure in the silicon hexaboride, both experimentally and theoretically, is found as a cubic phase, denoted as γ-SiB6-type. It appears in space group Pm-3m (no. 221), while full structural data are given in Tables 1, 2 and A1. Previous reports show unit cell parameters of a = 4.130 Å , which concur very well with our GGA (a = 4.161 Å ) and LDA (a = 4.109 Å ) results (Table 1). First reports of the SiB6 compound described it as a cubic phase and indicated that it is isomorphous with CaB6 [8]. Hence, within this γ-modification, boron atoms have five-fold coordination, with interatomic distances of 1 × 1.663 Å and 4 × 1.766 Å , and the structure is visualized in Figure 3a. Additionally, in the second coordination polyhedra (CP), silicon atoms form octahedra with atom-atom distances of 6 × 4.161 Å , which is visualized in Figure 3b. Furthermore, analyzing Si-B distances of silicon atoms were found to be surrounded by 24 boron atoms, with the distance between the atoms being 24 × 3.057 Å -B, as shown in Figure 3c. In addition, it has been confirmed that there is ionic bonding in the isomorphous CaB6 phase between the boron group and the Ca atom [63]. The last modification of four yielded from our searches was the rhombohedral structure, denoted as δ-SiB6-type, that crystallizes in space group P3m1 (no. 156) with unit cell parameters a = 3.503 and c = 6.407 Å (GGA). Besides silicon, there are two different boron atoms, B1 and B2, and full structural data are listed in Table 2 (LDA in Table A1). Within δ-SiB6-type of the structure represented with hexagonal axes, boron atoms are in six-fold coordination and form two distinct octahedra with atom-atom distances (B1 2 × 1.736 Å -

Elastic and Mechanical Properties of SiB6
Surprisingly, there is a limited number of investigations on the elastic and mechanical properties of SiB6, both in theory and in experiment. Experimental work is mostly devoted to non-stochiometric SiB6 or doped compounds [3,[64][65][66][67], while there are only a few recent theoretical studies on silicon hexaboride [9,11,17]. In this study, the elastic constants Cij for the most relevant silicon hexaboride modifications (α-SiB6, β-SiB6, and γ-SiB6) have been calculated using a GGA-PBE and LDA-PZ approach and were compared to previous theoretical data where found. Cubic γ-SiB6 modification has only three independent elastic constants, namely C11, C12, and C44, and they are given in Tables 4 and A2. These calculated elastic constants are in good agreement with available theoretical data from Ref. [9] ( Table 4).

Elastic and Mechanical Properties of SiB 6
Surprisingly, there is a limited number of investigations on the elastic and mechanical properties of SiB 6 , both in theory and in experiment. Experimental work is mostly devoted to non-stochiometric SiB 6 or doped compounds [3,[64][65][66][67], while there are only a few recent theoretical studies on silicon hexaboride [9,11,17]. In this study, the elastic constants Cij for the most relevant silicon hexaboride modifications (α-SiB 6 , β-SiB 6 , and γ-SiB 6 ) have been calculated using a GGA-PBE and LDA-PZ approach and were compared to previous theoretical data where found. Cubic γ-SiB 6 modification has only three independent elastic constants, namely C 11 , C 12 , and C 44 , and they are given in Tables 4 and A2. These calculated elastic constants are in good agreement with available theoretical data from Ref. [9] (Table 4).
Using the elastic constants, one can calculate the mechanical stability of the corresponding modifications using the mechanical stability criteria [68,69]. There are three conditions for cubic crystals: The results of the elastic constants (Tables 4 and A2) are not satisfying one of the mechanical stability criteria, since the C 44 constant is negative, indicating mechanical instability in the cubic γ-SiB 6 structure. This is in agreement with previous calculations, where γ-phase is found to be mechanically unstable [9].
In the case of the α-SiB 6 phase, these are the first reports of the elastic constants, calculated using LDA (Table 4) and GGA (Table A2). There are four conditions for hexagonal phases [68]. C 11 > |C 12 |; 2C 13 2 < C 33 (C 11 + C 12 ); C 44 > 0; C 66 > 0, and C 66 constant has been calculated with C 66 = 1/2(C 11 − C 12 ) [49,50]. The α-SiB 6 modification, calculated in a hexagonal setting, appears mechanically stable, and was calculated  Finally, we move to the β-SiB 6 phase with a lower orthorhombic symmetry and a larger number of independent elastic constants. They were all calculated using both the LDA and GGA approach, and concur with previous theoretical data [17] (Tables 4 and A2). There are six conditions for orthorhombic crystal system [68]: The calculated β-SiB 6 modification appears mechanically stable regardless of the DFT method applied, which is in agreement with previous calculations, where it has been found as mechanically and dynamically stable [17].
In this study, we have calculated bulk modulus B, shear modulus K, Young's modulus E, Poisson's ratio v, and Pugh's criterion B/K, for three SiB 6 modifications using LDA and GGA approximations. The value of the bulk modulus has been predicted between 146.44-153.43 GPa using GGA and 154.59-169.21 GPa using LDA for all three SiB 6 phases ( Table 5). This is in good agreement with other theoretical studies on the Si-B compound, where the bulk modulus was calculated in the range between 118 and 183 GPa [9,11,17]. Shear modulus has been predicted in a wide range from 22.21 to 93.88 GPa using GGA, and 23.63-92.09 GPa using the LDA method. This is in reasonable agreement, since literature data on other silicon borides show an even wider span of calculated B values 11,70]  Apart from the elastic moduli (B, K, and E) shown above, Poisson's ratio is another important mechanical property. Poisson's ratio, ν, is the negative ratio of the lateral or transverse strain to the axial strain in tensile loading and is thus interrelated with Young's modulus E [71,72]. Our DFT calculations show ν in the range between 0.24 and 0.43, regardless of the computational approach (Table 5), and are in a reasonable agreement with the values of 0.17-0.35 reported for various silicon borides [9,11,70]. Moreover, Poisson's ratio provides information about the ductility/brittleness of the materials. If the v value is smaller than 0.26, the material will have brittle behavior; otherwise, it is ductile. It appears that β-SiB 6 has a brittle character, while αand γ-phase appear to be ductile, after using both GGA and LDA methods (Table 5).
This can be further investigated with Pugh's forecast material delay/brittle empirical criterion (B/K) ( Table 5) [73]. The critical value of Pugh's criterion, which separates ductile and brittle materials, is around or higher than 1.75. If higher than this value, the material behaves in a ductile manner, otherwise, the material behaves in a brittle manner [74]. Again, according to the B/K relationship, the β-SiB 6 phase has a brittle character, while αand γ-modifications show ductile behavior regardless of the computational approach (Table 5).
In order to illustrate Young's and shear modulus anisotropy, we have plotted the 3D anisotropic surface figures of the K and E modulus under the spherical coordinates for the α-phase ( Figure 5). The content of anisotropy depends on the deviation degree from the spherical shape. The degree of deviation between the sphere and the surfaces which we obtained suggests a high degree of elastic anisotropy, especially in the shear modulus (Figure 5b), while Young's modulus shows the smallest deviation from the sphere in the xy plane ( Figure 5a). Moreover, 3D contour plots of the anisotropic surface figures of Young's and shear modulus for the β-SiB 6 modification have been performed ( Figure 6). In contrast to the α-phase, the β-SiB 6 shows a lower degree of elastic anisotropy in Young's modulus, since it is more spherical in all three planes (Figure 6a), while shear modulus is mostly spherical in the xy plane (Figure 6b). In addition, visualization and anisotropic analysis of the linear compressibility and Poisson's ratio in 3D for both αand β-phases have been presented in Appendix A. A smaller deviation degree from the spherical shape has been found for the linear compressibility, while Poisson's ratio largely deviates, indicating a high degree of elastic anisotropy (Figures A1 and A2).

Conclusions
Structure prediction and mechanical properties investigations of silicon hexaboride on an ab initio level have been performed. Data-mining-based searches over the ICSD database combined with the PCAE method produced 44 structure candidates, which, after full structural optimization using two DFT methods (LDA and GGA), have been reduced to four final SiB6 structures, regardless of the computational approach. Two novel structures are denoted as the α-SiB6-type, appearing in the hexagonal space group R-3mH (no. 166), and as δ-SiB6-type that crystallizes in the space group P3m1 (no. 156) have been predicted, for which there are no previous experimental or theoretical data. Our DFT calculations on the experimentally known cubic γ-SiB6 and recently proposed orthorhombic β-SiB6 phase are in very good agreement with previous findings. Elastic and mechanical properties of the predicted structures were investigated in the next phase. There is a limited number of such investigations and, in this study, the elastic constants for the most relevant silicon hexaboride modifications have been calculated using the GGA-PBE and LDA-PZ approaches, and were compared to previous theoretical data where found.
Calculated elastic constants are in good agreement with available theoretical data and show α-SiB6 and β-SiB6 as mechanically stable. Besides, we have calculated bulk

Conclusions
Structure prediction and mechanical properties investigations of silicon hexaboride on an ab initio level have been performed. Data-mining-based searches over the ICSD database combined with the PCAE method produced 44 structure candidates, which, after full structural optimization using two DFT methods (LDA and GGA), have been reduced to four final SiB6 structures, regardless of the computational approach. Two novel structures are denoted as the α-SiB6-type, appearing in the hexagonal space group R-3mH (no. 166), and as δ-SiB6-type that crystallizes in the space group P3m1 (no. 156) have been predicted, for which there are no previous experimental or theoretical data. Our DFT calculations on the experimentally known cubic γ-SiB6 and recently proposed orthorhombic β-SiB6 phase are in very good agreement with previous findings. Elastic and mechanical properties of the predicted structures were investigated in the next phase. There is a limited number of such investigations and, in this study, the elastic constants for the most relevant silicon hexaboride modifications have been calculated using the GGA-PBE and LDA-PZ approaches, and were compared to previous theoretical data where found.
Calculated elastic constants are in good agreement with available theoretical data and show α-SiB6 and β-SiB6 as mechanically stable. Besides, we have calculated bulk

Conclusions
Structure prediction and mechanical properties investigations of silicon hexaboride on an ab initio level have been performed. Data-mining-based searches over the ICSD database combined with the PCAE method produced 44 structure candidates, which, after full structural optimization using two DFT methods (LDA and GGA), have been reduced to four final SiB 6 structures, regardless of the computational approach. Two novel structures are denoted as the α-SiB 6 -type, appearing in the hexagonal space group R-3mH (no. 166), and as δ-SiB 6 -type that crystallizes in the space group P3m1 (no. 156) have been predicted, for which there are no previous experimental or theoretical data. Our DFT calculations on the experimentally known cubic γ-SiB 6 and recently proposed orthorhombic β-SiB 6 phase are in very good agreement with previous findings. Elastic and mechanical properties of the predicted structures were investigated in the next phase. There is a limited number of such investigations and, in this study, the elastic constants for the most relevant silicon hexaboride modifications have been calculated using the GGA-PBE and LDA-PZ approaches, and were compared to previous theoretical data where found.
Calculated elastic constants are in good agreement with available theoretical data and show α-SiB 6 and β-SiB 6 as mechanically stable. Besides, we have calculated bulk modulus B, shear modulus K, Young's modulus E, Poisson's ratio v, and Pugh's criterion B/K for various SiB 6 modifications using LDA and GGA approximations, and our results concur with other theoretical studies on the related Si-B compound. From the calculated Poisson's ratio and Pugh's criterion (B/K), it appears that β-SiB 6 will have brittle character, while αand γ-phase appear to be ductile using both GGA and LDA methods. In addition, we have plotted the 3D anisotropic surface figures of K and E modulus under the spherical coordinates for the αand β-phase. We believe that our results could potentially have a great impact on the industrial and technological applications of silicon boride-based materials.  Data Availability Statement: Additional data may be requested from the first author.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Table A1. Full structural details (modifications, space group, unit cell parameters, and atomic positions) for SiB 6 modifications obtained from data-mining-based searches and local optimization on the LDA-PZ level.  Table A2. Calculated elastic constants Cij (GPa) for various SiB 6 modifications using LDA approximation and compared to previous calculations.