Concentration-Diversified Magnetic and Electronic Properties of Halogen-Adsorbed Silicene

Diverse magnetic and electronic properties of halogen-adsorbed silicene are investigated by the first-principle theoretical framework, including the adatom-diversified geometric structures, atom-dominated energy bands, spatial spin density distributions, spatial charge density distributions and its variations, and orbital-projected density of states. Also, such physical quantities are sufficient to identify similar and different features in the double-side and single-side adsorptions. The former belongs to the concentration-depended finite gap semiconductors or p-type metals, while the latter display the valence energy bands with/without spin-splitting intersecting with the Fermi level. Both adsorption types show the halogen-related weakly dispersed bands at deep energies, the adatom-modified middle-energy σ bands, and the recovery of low-energy π bands during the decrease of the halogen concentrations. Such feature-rich band structures can be verified by the angle-resolved photoemission spectroscopy experiment.

fermion characteristic to charge carriers, which has been verified by angle-resolved photoemission spectroscopy and scanning tunneling spectroscopy measurements 35,36 . The possible high carrier mobility make silicene an ideal material, especially for the applications in field effect transistors 37 . Other fundamental properties of silicene have been revealed, such as a large energy gap opened by the spin-orbit coupling at Dirac point 38 , a quantum spin Hall effect observed in a low-temperature regime 39 , transition from a topological insulating phase to a band insulator generated by an electric field and electrically tunable energy gap 40,41 , and the emergence of a valley-polarized metal and anomalous quantum Hall effect 42 . Such potential properties identify silicene as a promising candidate to replace graphene not only due to the graphene-like features, but also its compatibility with silicon-based electronic devices.
Most intriguingly, the essential properties of silicene are extremely sensitive to chemical modifications. Among the chemical functionalizations, the atom dopings are one of the effective methods to dramatically change the fundamental properties. Hydrogen-adsorbed silicene exhibits semiconducting or metallic behavior that depend on the adatom concentrations 43 . The metalization of silicene, narrow gap semiconductor, and the semi-metallic or semiconducting behavior are revealed in alkali (Li, Na, K)-, alkali-earth (Be, Mg, Ca)-, and 3d transition metal atom (Ti, V, Cr, Mn, Fe, Co, Mo, and W)-adsorbed silicenes, respectively 44 . Besides, the chemical functionalization of silicene with boron (B), nitrogen (N), aluminum (Al), and phosphorus (P) adatoms have been explored in detail. Such systems present a metallic behavior with strongly bonded B, N, Al, and P adatoms, in which there exists the electron transfer from Si atoms to B, N, and P adatoms (p-type doping), and the opposite charge transfer (n-type doping) with weaker bonding for Al adatoms 45 . Furthermore, the electronic properties and optimal structure models of halogenated silicene have been explored by first-principles studies and low-temperature scanning tunneling microscopy measurements, respectively [46][47][48] . However, those studies mainly focused on a specific halogen concentration. As a result, the various halogenation effects on 2D silicene deserve a further look because of its promising applications, especially for electrode materials in lithium-ion battery.
In this work, the geometric, magnetic, and electronic properties of halogen-adsorbed silicene are fully investigated by the first-principle study. The double-side [ Fig. 1(a)] and single-side [ Fig. 1(b)] adsorptions are considered in calculations to evidence that these adsorption types can present the diverse magnetic and electronic properties, in which they are determined by the adatom-diversified geometric structures, atom-dominated energy bands, the spatial spin density distributions, the spatial charge density distributions and its difference, and the orbital-projected density of states (DOSs). Via such physical quantities, the similar and different characteristics in double-side and single-side adsorption cases are thoroughly identified, including the halogen-related weakly dispersed bands at deep energies, the adatom-modified middle-energy σ bands, the recovery of low-energy π bands during the decrease of the halogen concentrations; however, the breaking of the mirror symmetry in the single-side adsorptions leads to the critical discrepancy, i.e., there appear the valence energy bands with/without spin-splitting intersecting with the Fermi level as a result of p-type metallic behavior. The predicted results in the geometric structures, band structures, and DOSs can be verified by the scanning tunneling microscopy (STM), angle-resolved photoemission spectroscopy (ARPES), and scanning tunneling spectroscopy (STS) experiments, respectively.

Result and Discussion
Geometric structure. The ab initio calculations are performed in order to explore the geometric, magnetic, and electronic properties of halogen-adsorbed silicene under the various concentrations. Two typical adsorption types of fully halogenated silicene are chosen for this systematic study. The double-side adsorption, in which both sides of the silicene plane are thoroughly covered by halogen adatoms, as shown in Fig. 1(a). The halogen adatoms are uniformly distributed so that each adatom is separated by its six nearest adatoms with the same distance [ Fig. 1(c)], and the back-side adatoms (green balls) are set in the middle of three front-side adatoms (red balls). For the single-side adsorption, the adatoms only cover a single side of the silicene plane [ Fig. 1(b)], i.e., all the back-side adatoms are removed so that the concentration is reduced to half of the double-side adsorption. The typical concentrations vary from 100%, 12,5%, and 5,6%, where the ratios between halogen adatoms and Si atoms correspond to 2:2, 1:8 and 1:18, respectively [ Fig. 1(c)]; more cases are shown in Tables 1 and 2. The binding energy is calculated by E b = (E T − E P − n E A )/n, in which the E T , E P , and E A are, respectively, the energies of the fully halogen adsorbed silicene, pristine silicene, and isolated halogen adatoms; n represents the number of halogen adatoms. The lower binding energy corresponds to a higher stability. The most optimal geometric symmetry is used in the calculations. Via the detailed examinations, the most optimal adsorption position is situated at the top site (on the top of the upper silicon atom), as compared with the hollow site (above the center of the hexagonal silicon ring), valley site (on the top of the lower silicon atom), and bridge site (on the top of the Si-Si bond), irrespective of any doping case 45 . Among the halogen adatoms, the magnitude of binding energy declines with the increase of the atomic number [Tables 1 and 2], i.e., the F-adsorbed systems achieve the lowest binding energy − . 5 3 eV-−4.7 eV. Furthermore, the halogen-absorbed silicene [ Tables 1 and 2] possesses a higher geometric stability than halogen-adsorbed graphene 49 because of its highly reactive buckled surface. The shortest halogen-Si bond length is revealed in the F-adsorbed system .1 63 Å-1.64 Å, due to its smallest atomic number among halogen atoms. The stable halogen-Si bondings are formed by the transfer of electrons from Si atoms to halogen adatoms (discussed in the charge distribution section), and they are almost insensitive to the adatom concentration. The charge transfer renders the nearest Si-Si bonds weakened, i.e., the nearest Si-Si bond lengths are lengthened (2.31 Å-2.33 Å) as compared with the pristine Si-Si bond length (2.25 Å). The Si-Si bonds nearest to the halogen adatom are termed as the nearest Si-Si bonds. Besides, the halogen-Si-Si bond angles increase and buckling structures become obvious when the adatom concentration declines [Tables 1 and 2]. It is notable that in a pristine silicene, the planar Si-Si bondings are formed by Si-(3 s, 3p x , 3p y ) orbitals, while the π-π overlap of the parallel 3p z orbitals of nearest Si atoms reduces in strength, and thus Si atoms have a mix of sp 2 /sp 3 hybridization; therefore, the adatom adsorption on silicene could be easier to tune than purely sp 2 hybridized graphene. When the halogen adatoms are bonded with Si atoms, the Si atoms are changed from a mix of sp 2 /sp 3 hybridization to sp 3 hybridization. The chemical bonding scheme is determined by the hybridization between the Si-3p z orbitals and halogen-(p x , p y , and p z ) orbitals, and the weak sp 3 hybridization of four orbitals (3s, 3p x , 3p y , and 3p z ). The double-side and single-side adsorptions present similar geometric distortions; however, they might lead to different electronic and magnetic properties, mainly owing to the broken mirror symmetry [ Fig. 1(b)].
For experimental verification, STM is a powerful instrument for imaging surfaces at the atomic level. This tool can directly detect the surface structures in real-space under the atomic resolution, including the atomic lattice, the very short bond lengths, the planar or buckled structures, the achiral or chiral edges, the surface adsorptions, and the atom substitutions. Up to date, high-resolution STM observations have successfully verified the atomic structures of graphenes 50 , few-layer graphenes 51 , hydrogenation of graphenes 52 , fluorinated graphenes on copper 53 , 1D graphene nanoribbons 54 , silicene layers grown on Ag(111) 30 , silicene nanoribbons on Ag(110) 55 , hydrogenation of silicene films grown on Ag(111) 56 , chlorination of monolayer silicenes 48 , and multilayer silicenes 57 . Apparently, the predicted results in the buckled silicene after halogen adsorptions, including the top-site positions of halogen adatoms, halogen-Si & Si-Si bond lengths, buckling structures, and bond angles could be verified by high-resolution STM measurements, being very useful in confirming the multi-orbital hybridizations in the chemical bonds.  Fig. 2) owing to the extended π bondings of Si-3p z orbitals, and a negligible energy gap of 1 meV. The low-lying linear energy bands turn into parabolic dispersions with a saddle point at the M point, as shown in Fig. 3(a). These energy bands are mainly dominated by 3p z orbitals of the two nearest Si atoms, regarding as π bands. Moreover, the weak sp 3 hybridization illustrates an obvious separation of the π bands and σ bands within ±2 eV. The occupied and unoccupied σ parabolic bands are respectively initiated at −1 eV and 1.2 eV, indicated by the black and red rectangles in Fig. 3(a). The highest occupied bands display a double degeneration at the Γ point, while the lowest unoccupied bands are non-degenerate at the M point. The former bands are mainly contributed by the (3p x and 3p y ) orbitals, and they are hybridized with the 3p z orbitals at lower energy of −2.5 eV. The third occupied σ band initiates at −1 eV, dominating by the 3 s orbitals, and gradually becomes parabolically dispersed. This band becomes partially flat at −3.2 eV near the Γ point, in which there exist an orbital hybridization between the 3 s orbitals and 3p z orbitals (discussed in DOSs section). These fundamental features are dramatically modified after halogen adsorptions.
The double-side adsorption cases exhibit semiconducting or metallic behaviors, depending on distinct concentrations. The fully fluorinated, chlorinated, brominated, and iodinated silicene present a direct energy gap of 0.47 eV, 1.16 eV, 1.09 eV, and 0.51 eV, as shown in Fig. 3   www.nature.com/scientificreports www.nature.com/scientificreports/ section) and the weakened nearest Si-Si bond strength (discussed in the spatial charge distribution section). Among halogen adatoms, F adatom-related weakly dispersed bands dominate at the deepest energies due to its shortest F-Si bond length [blue circles in Fig. 3(b)]. This band dispersion becomes more obvious and gradually dominates at higher energies with an increase of its atomic numbers, mainly owing to the higher halogen-Si bond lengths [blue circles Fig. 3(e-h)].
As the adatom concentration declines, the halogen-Si bonds reduce in number, and rather compete with the effect of the weak sp 3 hybridization to determine the energy bands. Quite different from the 100% system, the 25% system presents an indirect energy gap of 0.75 eV, which is determined by HOS and LUS at the Γ and M points [ Fig. 3(c) and Table 1]. The critical concentration is revealed at 11% [ Fig. 3(d)], exhibiting the p-type metallic behavior, in which there exist valence energy bands intersecting with the Fermi level, so that the unoccupied   Table 2. Representation similar to that in Table 1, but shown for single-side adsorptions. www.nature.com/scientificreports www.nature.com/scientificreports/ valence states between the Fermi level and the top of valence bands belong to free holes. Also, the energy bands near the Fermi level with a weaker dispersion are mainly contributed by the non-passivated Si-3p z orbitals [ Fig. 3(c,d)], revealing that the low-lying energy bands gradually changed from the σ bands to the π bands. Besides, the middle-energy σ band of Si-(3p x and 3p y ) orbitals in pristine silicene [ Fig. 3(a)] is located − .  www.nature.com/scientificreports www.nature.com/scientificreports/ non-passivated Si atoms, in which the latter atoms mainly contribute to the bands in the 11% system because of its higher concentration. Notably, the direct energy gap in the 100% system becomes indirect when the concentration reduces up to 25% [ Table 1], and any halogen concentration lower than the critical concentration of 11% becomes the p-type metals. This remains true for most of the halogen adatoms. These revealed band features are sufficient to comprehend the concentration-dependent electronic properties.
In the single-side adsorptions, the destruction of the mirror symmetry [ Fig. 1(b)] creates the valence energy bands intersecting with the Fermi level, in which the Fermi level is situated at the valence band, so the unoccupied valence states between the Fermi level and the half-occupied valence bands belong to free holes [ Fig. 4]. As a result, these systems can also be regarded as the p-type metals. The full single-side adsorption cases [ Fig. 4(a,d,f-h)] possess valence parabolic bands intersecting with the Fermi level, which are only related to the three non-passivated Si atoms [gray circles in Fig. 1(b)] nearest to the passivated Si atoms [red circles in Fig. 1(b)]. The parabolic dispersion change to the partially flat bands when the concentration declines, as shown in Fig. 4(b/e,c) for 12.5% and 5.6% systems, respectively. Furthermore, the π bondings between non-passivated Si atoms and passivated Si atoms [ Fig. 1(b)] is seriously suppressed owing to the high charge transfer from the latter Si atoms and the halogen adatoms. This affirms that the non-passivated Si atoms can not form the extended π electronic states via the interaction with the passivated Si atoms. The clear discrepancy between the single-side and double-side adsorption cases is the valence energy bands intersecting with the Fermi level. Also, the halogen-induced weakly dispersed bands in the double-side adsorption cases [ Fig. 3] reduce in number due to its decreased concentration [ Fig. 4].
As compared with the pristine case [ Fig. 3(a)], the Dirac cone structure is seriously distorted in the fully fluorinated case [ Fig. 3(b)], and thoroughly destroyed in other cases, mainly owing to the seriously destruction or termination of the π bondings. Also, there exist other similar features, such as the recovery of low-energy π bands in the low-concentration systems, the adatom-modified middle-energy σ bands, and the halogen-induced weakly dispersed bands at deep energies. When the concentration decreases, both 5.6% single-side [ Fig. 4(c)] and 11% double-side [ Fig. 3(d)] systems show the low-lying π bands with large energy widths. These low-energy π bands are mainly contributed by the non-passivated Si-3p z orbitals, indicating that the π bondings are recovered among the Si atoms. The middle-energy σ bands (−0.2 eV-−3.5 eV) [Figs 3(b-d) and 4(a-c)] are very similar to each other in terms of their band width, initiated energy, and orbital contribution, mainly owing to the hardly affected σ bondings. In addition, the halogen-related weakly dispersed bands are confined at similar energies in both single-side and double-side systems since the stable halogen-Si bondings are insensitive with the halogen concentrations (discussed in geometric structure section). That is to conclude that the broken geometric symmetries hardly have an impact on the aforementioned similar features.
To date, ARPES is an unique tool to verify the band structures. This experimental equipment can directly measure the direction, speed, and scattering process of valence electrons in the sample being studied. The high-resolution ARPES experiments have successfully verified the feature-rich band structures of graphenes on Ir(111) 58 , few-layer graphenes 59 , halogenated graphenes 60 , 1D graphene nanoribbons 61 , monolayer silicenes 35 , hydrogenated silicenes 62 , and multilayer silicene nanoribbons 63 . Apparently, similar ARPES measurements can be utilized to examine the feature-rich band structures after halogen adsorptions, including the destruction or recovery of low-energy π bands, the halogen-induced weakly dispersed bands at deep energies, and the valence energy bands intersecting with the Fermi level. Spatial spin distribution. The highly geometric symmetry in the double-side adsorption cases [ Fig. 1(a)] thoroughly destroys the spatial spin orientation, so that their energy bands are degenerate. When the geometric symmetry no longer exists [ Fig. 1(b)], the spatial spin density distributions can be performed in the single-side adsorption cases, leading to the spin-splitting energy bands. Only except for the 50% fluorinated system [ Fig. 4(a)] exhibits the degenerate parabolic band. The other 50% systems have the spin-splitting parabolic bands, which can clearly observe near the Fermi level, as shown by red (spin up) and black (spin down) curves in Fig. 4(d,f-h). These spin-splitting parabolic bands are mainly contributed by the non-passivated Si atoms, verified by the spatial spin density distributions near the non-passivated Si atoms (pure spin up configuration) [ Fig. 5(a-c)]. However, the spin-splitting bands become weaker in the 12.5% systems [ Fig. 4(b,e)]. This indicates that the spin-up and spin-down states can co-exist in its spin density distributions under the decrease of the halogen concentration (spin up and spin down mixed configuration), as shown by red (spin up state) and blue (spin down state) balls in Fig. 5(d-f). The difference in the spin up and spin down states directly determines the strength of the net magnetic moments, i.e., 0.72 μ B , 0.89 μ B , and 0.79 μ B in the 50% systems and 0.64 μ B , 0.67 μ B , and 0.45 μ B in the 12.5% systems [ Table 2]. Notably, the spin density distributions are fully absence at the critical concentration of 5.6% (vanishing the net magnetic moment in Table 2), thus its energy bands are highly degenerate [ Fig. 4(c)]. Furthermore, any halogen concentration beyond the critical concentration results in the absence of the spin density distributions. This remains true for all halogen adatoms. The single-side adsorption-induced magnetic properties on the silicene surface can be verified by the spin-polarized STM.

Spatial charge distribution. The spatial charge density distribution (ρ) can illustrate the strength of the
Si-Si and halogen-Si bondings, as shown in Fig. 6(a-g). ρ displays a strong σ bonding between two Si atoms in the pristine silicene [ Fig. 6(a)] or two non-passivated Si atoms in the low-concentration systems [ Fig. 6(f,g)], as presented by the enclosed black rectangle. The strength of such σ bonding becomes weaker when Si atoms are bonded to halogen adatoms, as indicated by enclosed gray rectangle in Fig. 6(b-g), resulting in the extended Si-Si bond length. This accounts for the similarity in the adatom-modified middle-energy σ bands between the single-side and double-side adsorption cases. Away from the horizontal line between two Si atoms, ρ is lowered; however, its charge distribution is extended, evidencing for the π-electronic states, as shown by enclosed red rectangle in Fig. 6(a,f,g). Nevertheless, hardly extended states, corresponding to the π-charge-depleted region, are www.nature.com/scientificreports www.nature.com/scientificreports/ perceived between the passivated and non-passivated Si atoms (enclosed purple rectangle in Fig. 6(d-g)). This charge distribution is similar to isolated Si atoms and evidences for the destruction of π bondings. Furthermore, the charge density of halogen-Si bondings is much higher than that of Si-Si bondings, as indicated by the dashed black rectangle in Fig. 6(b-g). Thus, such stable halogen-Si bondings are insensitive to the halogen concentrations. This illustrates that the halogen-Si bonding-induced weakly dispersed bands are limited to similar energies in both single-side and double-side systems. In order to further understand the detailed charge transfers among all orbitals, the charge density difference (Δρ) is presented in Fig. 6(h-m). Δρ is created by subtracting www.nature.com/scientificreports www.nature.com/scientificreports/ the charge density of the pristine silicene and isolated halogen adatoms from that of the composite system. All the halogen-adsorbed systems exhibit a high charge transfer from Si atoms to halogen adatoms to form the stable halogen-Si bondings (dashed black rectangle in Fig. 6(h-m)). Also, such clear variation of charge density accounts for a weaker bond strength of the nearest Si-Si bonds around the halogen adatoms.
Diverse density of states. The orbital-projected density of states (DOSs) is used to verify the band structures and comprehend the orbital hybridizations in the chemical bonds. For pristine silicene [ Fig. 7(a)], the low-energy DOSs exhibits the vanished value at E F = 0, V-shaped structure, and two symmetric logarithmic divergent peaks at ~± E 1 c v , eV, are induced by Si-3p z orbitals (purple solid curve in Fig. 7(a)). They respectively originate from the negligible-energy gap, the Dirac cone structure, and the saddle points of π bands at the M point [ Fig. 3(a)]. The middle-energy DOSs shows a shoulder structure at −1 eV and a prominent peak at −2.6 eV, resulting from the maximum point and the saddle point of σ bands, respectively. They are mainly contributed by Si-(3p x and 3p y ) orbitals (green and blue solid curves in Fig. 7(a)). Moreover, the hybridization between Si-3s (red solid curve in Fig. 7(a)) and Si-3p z orbitals creates a prominent peak at −3.2 eV, coming from the partially flat  Fig. 3(h)]. The middle-energy DOSs displays a shoulder structure at −0.2 eV and the prominent peak at −1.5 eV, are dominated by Si-(3p x and 3p y ) orbitals (green and blue solid curves in Fig. 7(b)), are very similar to that of pristine silicene [ Fig. 7(a)], accounting for the adatom-modified middle-energy σ band. Such modified middle-energy DOSs structure is also observed in Fig. 7(e-h). Furthermore, the extra strong peaks at deep energies are dominated by halogen-(2p x and 2p y ) orbitals www.nature.com/scientificreports www.nature.com/scientificreports/ [ Fig. 7(b)], standing for the halogen-related weakly dispersed bands at deep energies [ Fig. 3(b)]. These strong peaks become weaker and dominate at higher energies as their atomic numbers increasing [ Fig. 7(e-h)]. When the halogen concentration declines, three low-energy prominent peaks appear at −0.6 eV, 0.8 eV, and 1.4 eV in the 25% system [ Fig. 7(c)]. The bandwidth of these prominent peaks about 0.6 eV, and their states are related to Si-3p z orbitals (purple solid curve). Also, their middle-energy DOSs are contributed by Si-(3p x and 3p y ) orbitals; however, the single peak structure in the full adsorption case [ Fig. 7(b)] separate into several sub-peaks owing to an enhancement of sub-bands. At the critical concentration of 11%, most of the low-energy DOSs within ±1 eV are mainly dominated by the non-passivated Si-3p z orbitals [ Fig. 7(d)], illustrating a π bandwidth similar to the www.nature.com/scientificreports www.nature.com/scientificreports/ pristine case [ Fig. 7(a)]. This evidences the recovery of low-energy π bands in low-concentration systems. Also, these two low-concentration systems show the halogen-dominated extra strong peaks at deep energies similar to the 100% system [ Fig. 7(b)]. However, the intensity of these deep-energy extra strong peaks is roughly proportional to the halogen concentration.
The fundamental features of DOSs in the single-side and double-side adsorption cases are the strong peaks, which are closely associated with the Si-and halogen-related bondings. The similarities are identified between the single-side cases and its corresponding double-side adsorptions, including the extended low-energy π band widths [Figs 7(d) and 8(c)], the three Si-3p z orbital-related prominent peaks within −1 eV to 2 eV [Figs 7(c) and 8(b)], www.nature.com/scientificreports www.nature.com/scientificreports/ the Si-(3p x and 3p y ) orbital-dominated single-peak structure [Figs 7(b) and 8(a)], and the halogen-related extra strong peaks at deep energies Figs 7(b) and 8(a); 7(c) and 8(b); 7(d) and 8(c)]. These extra strong peaks are confined at the same energy range; however, their intensities are gradually declined as the halogen concentration reduced. The similar DOSs features is also found in other halogen-adsorption cases [Figs 7(e-h) and 8(d-h)]. On the other hand, the breaking of the mirror symmetry in the single-side adsorption cases leads to a critical discrepancy, i.e., there exist the finite value of DOSs at the Fermi level as a result of the metallic behaviors [ Fig. 8(a-h)]. Also, the asymmetric prominent peaks near the Fermi level [ Fig. 8(b,d-h)] provide a very useful information in confirming the spin-splitting energy bands [ Fig. 3].
On the experimental side, the STS measurements, where the tunneling differential conductance (dI/dV) is proportional to DOSs, can provide sufficient information in examining the special features in DOSs. High-resolution STS measurements are able to distinguish the semiconducting, semi-metallic and metallic behavior. Moreover, they can be used to identify the close relations between the electronic energy spectra and the orbital hybridizations of the chemical bonds. To date, such experimental measurements have been successfully used to confirm the electronic band structure near the Fermi level and the dimension-diversified van Hove singularities in monolayer graphene 64 and few-layer graphene systems [65][66][67] , 1D graphene nanoribbons 68,69 , adatom-adsorbed graphene 70 , monolayer silicene 71 , and hydrogenated silicene 72 . Obviously, the theoretical predictions on the halogen-diversified DOSs can be thoroughly affirmed by the STS experiments, including the vanishing or existing of the finite value of DOSs at the Fermi level, the Si-or halogen-dominated strong peaks, and the halogen-modified σ shoulder structures. Furthermore, spin-polarized STS is available for verifying the asymmetric peaks near the Fermi level in DOSs.

Conclusion
The geometric, magnetic, and electronic properties of halogen-adsorbed silicene are studied by the ab initio calculations. The rich physical and chemical properties are diversified by the halogen concentrations. Specifically, the stable halogen-Si bondings are the critical factor in affecting the geometric structures, spatial spin density distributions, spatial charge density distributions and its variations, atom-dominated energy bands, and DOSs. The transfer of electrons from Si atoms to halogen adatoms leads to variations in the charge density distributions in the chemical bonds, directly resulting in the drastically modulation of the magnetic and electronic properties. The rich band structures after halogen adsorptions are highlighted by the destruction or recovery of low-lying π bands, the halogen-induced weakly dispersed bands at deep energies, and the valence energy bands with/without spin-splitting intersecting with the Fermi level. The first feature is determined by the competition between the halogen-Si bondings and the weak sp 3 hybridization. The energy gap determined by the σ bands is indirect in the full double-side systems. As the concentration declines, the low-lying π bands gradually recover, rendering the gap size smaller with an indirect behavior in the 25% system, and even though becoming the p-type metals at the critical concentration of 11%. The second feature, mainly owing to the stable halogen-Si bonding insensitive with the halogen concentrations, is the similar characteristic of both single-side and double-side adsorptions. The third feature, originating from the broken mirror symmetry, is a distinct characteristic in the single-side adsorptions. Also, the spatial spin density distributions only reveal in the specific single-side systems, resulting in the spin-splitting energy bands or the asymmetric prominent peaks near the Fermi level in DOSs. The above-mentioned band properties are fully confirmed by their DOSs, including the vacant region centered at the Fermi level, the recovery of π band-related low-energy prominent peaks, the halogen-related extra strong peaks at deep energies, and the finite value of DOSs at the Fermi level. The feature-rich band structures and the diverse DOSs can be verified by the ARPES and STS measurements, respectively. The concentration-tuned diverse magnetic and electronic properties of halogen-adsorbed silicene are very potential for applications in both nanoelectronics and spintronics.

Method
The diverse magnetic and electronic properties of halogen-adsorbed silicene are investigated by the spin-polarized density functional theory (DFT) implemented in VASP (Vienna ab initio simulation package) 73 . The exchange and correlation energies, which come from the many-particle electron-electron interactions, are evaluated by the Perdew-Burke-Ernzerhof (PBE) functional under the generalized gradient approximation 74 . Furthermore, the projector-augmented wave (PAW) pseudopotentials can characterize the electron-ion interactions. The plane-waves basis set with a kinetic energy cutoff of 500 eV is utilized to calculate wave functions and state energies. A vacuum space of 10 Å is inserted between periodic planes to avoid the interaction of two planes. The first Brillouin zone is sampled by 12 × 12 × 1 and 100 × 100 × 1 k-point meshes within the Monkhorst-Pack scheme for geometric optimizations and electronic structure calculations, respectively. Such points are sufficient for obtaining reliable DOSs. The convergence for energy is set to be 10 −5 eV between two consecutive steps, and the maximum Hellmann-Feynman force acting on each atom is less than 0.01 eV/Å during the ionic relaxations.