Effect of Doping in Controlling the Structure ... - ACS Publications

Jun 29, 2017 - calculations. Their structural features can be understood on the basis of. Vegard,s law. The lattice parameters were found to correlate...
0 downloads 0 Views 2MB Size
Subscriber access provided by BOSTON UNIV

Article

Effect of Doping in Controlling Reactivity, Electronic and Optical Properties of Pristine and Ca(II) Intercalated Layered Silicene Tamiru Teshome, and Ayan Datta J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b03002 • Publication Date (Web): 29 Jun 2017 Downloaded from http://pubs.acs.org on July 5, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Effect of Doping in Controlling Structure, Reactivity and Electronic Properties of Pristine and Ca(II) Intercalated Layered Silicene Tamiru Teshome and Ayan Datta* Department of Spectroscopy, Indian Association for the Cultivation of Science, 2A and 2B Raja S. C. Mullick Road, Jadavpur – 700032, Kolkata, WestBengal, India. Email: [email protected]

Abstract

Hybrid two-dimensional (2D) materials composed of carbon and germanium doped silicene monolayers, Si1xCx,

Si1-xGex and Si1-2xCxGex (0 ≤ x ≤ 1) are investigated on the basis of first-principles calculations. Their

structural features can be understood on the basis of the Vegard’s Law. The lattice parameters are shown to correlate well with the arithmetic mean of the percentage doping of individual constituents. The electronic band gap show interesting bell-shaped behavior for Si1-xCx with respect to doping wherein Eg increases from 0 (for x = 0) to a maxima (for x = 0.5) and eventually again decrease to 0 (for x = 1). Clearly controlled carbon doping of silicene monolayer opens up the band-gaps and also provide a tool to harvest solar energy for semiconductor and photovoltaic applications. Apart from pristine monolayer silicene, the effect of Ca(II) intercalated multilayer silicene/germanene as available in the van der Waals mineral phases of CaSi2 and CaGe2 are also investigated. The stability and reaction energies for hydrogenation, fluorination and mixed hydrogenationfluorination for these Ca(II) intercalated structures are compared vis-a-vis pristine and doped monolayer silicene. Optical absorption calculations demonstrate that doping monolayers as well as Ca(II) intercalated multilayer silicene/germanene have strong photo-absorption in visible light region. The presence of Ca(II) stabilizes the silicene/germanene layers and also the Ca(II) is affecting the electronic structure of these 2D materials. 1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 30

Introduction Graphene, a 2D honeycomb structure of atomic carbon, has attracted attention due to its remarkable electronic and magnetic properties for potential advanced technology applications.1, 2,

3, 4

In connection with, there has

been an active interest to study the honeycomb lattices of other group-IV elements. Silicene, the all Si analogue of graphene, that has recently become extremely popular5, is a single silicon monolayer with buckled honeycomb structure6 and has useful chemical and physical properties7, 8, including ferromagnetism, 10

9,

halfmetallicity, 11quantum hall effect12 and giant magnetoresistance.

Two-Dimensional Honeycomb Structures of Silicon and Germanium13 have been studied. Unlike graphene, two-dimensional Si and Ge layers are buckled due to pseudo Jahn-Teller distortion.14 Recently, experimental synthesis of germanene on metal substrates has been achieved.15,

16

First-principles calculations on silicene,

graphene and germanene indicate that the π and π* bonds linearly cross at the Fermi level, reflecting zero band gap semiconducting character.17 This unfortunately makes both pristine silicene and germanene unsuitable in electronic devices with logic operation. To overcome this, extensive efforts have been devoted to realize controllable band gaps in these materials. Several strategies such as patterning into nanoribbons18 and using chemical functionalization19,

20

have been proposed to effectively open a band gap in silicene, graphene and

germanene. Silicane and germanane (fully hydrogenated silicene and germanene) are predicted to be wide band gap semiconductors, with energy gap values ranging between 2.9-4.0 eV for silicane and between 2.9-3.6 eV for germanane.21, 22 Similarly the band structure and optical properties of layered germanane (a van-der Waals solids) 23 has also been studied. Recently, we have shown that blue phosphorene (an indirect band gap system) can be readily converted into a direct band gap n-type semiconductor based on doping with CT molecules like TTF.24 Although intensive theoretical studies on mono- and multilayer silicene have been reported, isolated freestanding silicene has not been obtained experimentally. This is probably because of the inherent air instability of 2D silicon and high reactivity of Si-Si bonds. Even though the computed freestanding silicene and germanene structures are strongly buckled, it has been shown that this buckling can be significantly reduced by 2 ACS Paragon Plus Environment

Page 3 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

stacked functionalized and confined in bilayer graphene.

25, 26, 27

Recent synthesis of silicene28,

29

and

germanene30 have been a major breakthrough in this area because of their enhanced spin-orbit coupling31, 32 and higher electron correlations.33 However, the underlying metal substrate, which is necessary for the growth of these honeycomb lattices, strongly affects the Dirac electron behavior.34 As a result, there has an interest in layered Ca intercalated silicene/germanene system like CaSi2 and CaGe2 for which strong layer and substrate interactions are circumvented.35

Intercalated layered materials are of significant interest due to wide tunability of their conducting properties regulated by the interface states.36, 37, 38 Experimental progresses in recent times enables the synthesis of a wide range of layered materials with subsequent intercalations having precise control over the stoichiometry and number of layers.39 For multilayer silicene, strong Si-Si interactions (interlayer bonding) destroy the sp2 nature of the electrons.40 These limitations may be avoided by means of intercalation in a multilayer system wherein the substrate effects and interlayer sp3 bonding are absent41. Therefore the 2D sp2 electrons might be preserved. Although freestanding silicene probably cannot exist in nature, the material has been successfully deposited on metallic substrates such as ZrB2 (0001), Ir(111), and Ag(111).28, 42- 45 On the other hand, several semiconductors have been explored theoretically to overcome the strong interaction characteristic for metallic substrates. For example, Ding and Wang have reported that GaS nanosheets could preserve the linearly dispersing bands with a gap of 0.17 eV, but the lattice mismatch would amount to 7.5% .46 GaS preserves the Dirac cone of silicene but a large lattice mismatch of 7.5% is a major drawback. Although the lattice mismatch is smaller (2.3 %) for ZnS(0001) the substrate also perturbs the band structure.47 Bhattacharya and co-workers have investigated IIVII and III-VI semiconductor substrates such as AlAs(111) and ZnSe(111).48 They reported that the stability and electronic properties of these systems depend largely on the topmost atoms of the substrate. WSe2 would fit much better (lattice mismatch of 0.6%) but the silicene monolayer develops a large band gap of 0.3 eV.49 The excellent electronic properties and low thermal conductivity make silicene very promising for thermoelectric devices.50,

51

The structural and electronic properties of silicene on MgX2(X = Cl, Br and I) combines a 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 30

reasonably small lattice mismatch with an almost gapless Dirac cone and a Fermi velocity comparable to freestanding silicene as reported by Jiajie Zhu and Schwingenschögl.52 Recently, Noguchi et al. reported the existence of Dirac dispersion in the Ca(II) intercalated multilayer silicene system CaSi2 using angle resolved photo emission spectroscopy (ARPES). 53

In this article, we propose an effective method to improve the electronic properties of silicene/germanene doping with C atom at different concentration ratio increases and decreases the band gap of the monolayer. Alloying Si with C and Ge atoms in such 2D monolayers systems will be particularly interesting as it promises composition controlled modulation of electronic and optical properties. We report a systematic study of monolayer silicene/germanene along with their Ca(II) intercalated multilayer employing first-principles calculations within the framework of DFT. To the best of our knowledge, there has been no previous report for the reactivity of Ca(II) intercalated silicene/germanene towards hydrogenation/fluorination or mixed hydrofluorination and their optical properties. This work is very significant for the silicene field as it suggests alternative bilayer, tri-layer and multilayer compounds that expand ones imagination to go and synthesize them.

Computational Details In this work, all the calculations are performed within the framework of density functional theory (DFT) using the project-augmented plane-wave method, as implemented in Vienna ab-initio simulation package (VASP). 54, 55

The generalized gradient approximation of Perdew, Burke and Ernzerhof (GGA-PBE) is selected for the

exchange correlation potentials56, and the DFT is employed. The cutoff energy for the plane-wave basis is set to 500 eV and the structures are relaxed until the residual forces on the atoms become smaller than 0.01 eV/ Å. Brillouin zone integration are carried out with 5 x 5 x 1 k meshes for the structure optimization and 7 x 7 x 1 k meshes for optical spectra and band structure calculations. We analyze Dirac dispersion in several possible stacking sequences of AA and AB silicene/germanene intercalated compounds along with their bilayer, tri-layer sandwich and multilayer structures.57 The optical properties were computed utilizing the time-dependent (TD) 4 ACS Paragon Plus Environment

Page 5 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

DFT scheme for the inclusion of electron-hole interactions58 within the random phase approximation (RPA) approach for which the local field effects are included at Hartree level. Only inter-band transitions are taken into account. For a transitionally invariant system, the Fourier transform of the frequency dependent symmetric dielectric matrix in the RPA is given by.59 

ε, q, ω = δ, q, ω − |||

χ q, ω | ,

(1)

Where G and G’ are reciprocal lattice vectors and q stands for the Bloch vector of the incident wave. The matrix χ q, ω is the irreducible polarizability matrix in case of the independent particle derived by Adler and Wiser60, 61 in the context of the self-consistent field approach. Optical properties of any system are in general calculated with the help of a frequency dependent complex dielectric function ɛ ω = ɛ ω + iɛ ω . The imaginary part is determined by a summation over empty states using the equation: ɛ ω =

4!  e 1 lim  ) 2ω+ δ E,/ − E./ − ω × 〈34+ + 56& 738+ 〉〈34+ + 5:& 738+ 〉 ∗ Ω &→ ( ,.../

where ω, the frequency of electromagnetic (EM) radiation in energy unit is Ω represents the volume of the primitive cell, q is the electron momentum operator, c and v are the conduction band and valance state, respectively,ω+ , is the k point weight E,/ , E./ and 34+ , 38+ are the eigenvalues and wave function at the k point, respectively and 56 , 5: are unit vectors for the three Cartesian direction. The phonon calculations are carried out using the PHONOPY code combined with density functional perturbation theory (DFPT) method in VASP.

Results and Discussion A. Structures and Stability The necessary and sufficient condition for dynamical stability of a crystal at low temperature is phonon stability. We have examined dynamical stability of silicene, germanene, graphene, and CaSi2 and CaGe2 phonon dispersion along the high symmetry lines ᴦ-K-M-ᴦ as shown in Fig. 1 (a-e), respectively. The absence of soft modes within the entire Brillouin zone clearly indicating that silicene, germanene, graphene, CaSi2 and CaGe2 are indeed stable. 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 30

We considered 2D Si-C and Si-Ge hybrids with different Si concentrations. First, we produced a 3 x 3 unit cell of silicene, in which silicon atoms are gradually substituted by carbon/germanium atoms to obtain different Si/C and Si/Ge ratios. After carefully considering the possible structures of all of these configurations for Si1-xCx, Ge1-xCx and Si1-xGex including isomeric forms (different structures with the same composition), the energies are compared to obtain the most stable configuration. The suitability of our level of theory might be gauged from the reproducible band gap of SiC and Dirac cones for pristine monolayer silicene and germanene. Therefore, we employed our method to investigate the effect of doping and optical properties of the 2D monolayers.

The well-known Vegard’s law (Eqn. 2) was utilized to understand the qualitative hybrid structures. According to Vegard’s law, unit cell parameters should vary linearly with composition for a continuous and constant temperature substitution in a solid solution for which the substituted atoms/ions are randomly distributed. The lattice parameters are controlled by the relative size of the atoms or species exchanged.62, 63 Table 1: Doping ratio (x), lattice constant “a” (Å), Energy gap Eg (eV) and optical excitation energy Ee (eV) for Si1-xCx,Si1-xGex and Ge1-xCx. doping Ratio(x)

C-doping silicene

Ge-doping silicene

C-doping germanene

a

Eg

Ee

a

Eg

Ee

a

Eg

Ee

0

3.86

0

4.0

3.86

0

4.0

4.06

0

4.0

0.111

3.69

0.32

4.0

3.88

0.02

3.75

3.77

0.27

1.95

0.222

3.52

0.51

5.14

3.90

0.02

3.06

3.65

0.36

2.01

0.333

3.15

1.99

5.31

3.92

0.04

3.70

3.53

0.49

2.5

0.5

3.09

2.55

5.35

3.95

0.05

1.69

3.26

2.07

3.5

0.556

3.03

1.38

2.63

3.96

0.03

0.86

3.14

1.02

3.65

0.778

2.73

0.65

1.75

4.01

0.03

0.75

2.69

0.80

2.1

1.0

2.46

0

4.0

4.06

0

0.61

2.46

0

4.0

The simplest mathematical expression for Vegard’s law for a binary solid solution A-B is: 6 ACS Paragon Plus Environment

Page 7 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

< = a> 1 − ? + a@ ?

(2)

where x = ?A is the mole fraction of component B a = a> and a = a@ are the lattice parameters of pure components A and B, respectively. The doping silicene sheets were generated without prior bias or preference for any structure by replacing the Si atoms by C/Ge and C-Ge. The concentration of individual C and Ge doping were varied from 0 ≤ x ≤ 1 where x is the relative composition of one type of atom with respect to the total number of atoms in the unit cell. Positions of the doping have also been varied for the different concentration of substitution doping to analyses the effect of doping. Note that, for multiple configurations at the same doping concentration, the most stable structure was chosen (supporting information Fig. S1and S2). Graphene is a pure planar monolayer as shown in Fig. 2 (I(i)) while silicene and germanene undergo slight out-of-plane buckling with buckling height, dH = 0.45 Å and 0.69 Å, respectively as consequence of pseudo Jahn-Teller distortion as discussed above. The optimized structures of 2D Si1-xCx, Si1-xGex, Ge1-xCx and Si1-2xCxGex and doping at different ratio are listed in Fig. 2 (I-IV). The isotropic silicene-graphene (Si0.5C0.5) hybrid is stable in a planar hexagonal lattice with Si-C bond length of 1.78 Å. With the increases of the substitution ratio-of C, the buckling of the planar structure also decreases. In Fig. 2(II and III), the optimized structures of 2D Si1-xGex and Ge1-xCx monolayers at different doping ratios are shown. The Ge0.5C0.5 (GeC) monolayer in planar or lower buckled structure obtained germanene doping carbon due to the symmetry destroyed and electronegativity difference between the C and Ge atoms; hence, GeC has a polar structure. In contrast, with the increases of Ge doping silicene, the buckling height increases in between silicene-germanene (dH=0.45-0.69 Å). As shown that in Fig. 2(IV), for Si1-2xCxGex we have C and Ge co-doping silicene as obtained by increasing the concentration ratios from x = 0.06 to x = 0.5.

The variation of the lattice constants and average Si-X (X = C, Ge) bond lengths for Si1-xCx, Si1-xGex and Ge1xCx

with respect to doping are shown in Fig. 3((b), (d) and (f)), respectively. The x=0 (un-doping) substitution

ratios correspond to pristine silicene sheet, while 1(100%) substitution ratio graphene (germanene) monolayer. For Si1-xCx and Ge1-xCx, increasing x leads to shrinkage of the lattice constants as shown in Fig. 3(a) and Fig. 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 30

3(e) which is concomitant with increasing the fraction of the shorter Si-C and Ge-C bonds at the expense of the longer Si-Si and Ge-Ge bonds. Thus, evidences from reduced average bond-length in the hybrid 2D monolayers (Si1-xCx and Ge1-xCx). The lattice constant decreases linearly from 3.86 Å for pure silicene (x=0) to 2.47 Å for pure graphene (x=1) and lattice constant decreases from germanene 4.06 Å to graphene 2.47 Å as shown in Fig. 3(a) and Fig. 3(e), respectively. Similar trend is also observed for Si1-xGex for which sequential substitution of Si by Ge leads to an expanded lattice constant in Fig. 3(c) which again correlates well with an increasing fraction of the longer Si-Ge bonds within a shorter Si-Si silicene monolayer in Fig. 3(d). Here again, the lattice constant monotonically increases from 3.86 Å from pure silicene to 4.06 Å for pure germanene (x=1). Within the hybrid layers, the bond lengths of Si-C and Si-Ge are in the ranges of 1.75 Å - 1.84 Å and 2.34 Å - 2.40 Å which are in the fraction of x=0.111 - 0.778. Clearly, the effect of doping carbon and germanium on silicene element is specific as C doping reduces the lattice constant while it is vice-versa for Ge. Nevertheless, for both the cases, an increment or decrement of lattice constants is correlated linearly with doping concentration in harmony with the Vegard’s law.

We further investigate the bilayer sandwiched compound silicene-Ca-silicene. Although this compound has not yet been synthesized experimentally, it represents an intercalated bilayer exfoliated from bulk CaSi2. In fact recently, based on molecular dynamics simulations, we have shown that liquid-phase exfoliation provides an excellent synthetic route to produce mono to few layer 2D materials 64. The atomic structures for several of Ca intercalated silicene/germanene multilayer where in buckled silicene/germanene layers and trigonal Ca layers are alternatively stacked (supporting information Fig. S5-S9). Here, we choose to present the electronic properties of bilayer, tri-layer sandwich together with multilayer CaSi2 and CaGe2 with AA and AB stacking sequences. For further confirmation of momentum shifts of Dirac cones, we show the two-dimensional dispersion in both shifted the ᴦ-k-M plane and A-H-L plane of the Brillouin zone. This happens actually due to the enhanced hopping within the silicene/germanene layers, which further increases in the presence of intercalated Ca atoms. Furthermore, A sublattice points of one silicene/germanene layer and B sublattice points 8 ACS Paragon Plus Environment

Page 9 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of the adjacent silicene/germanene layer in CaSi2 and CaGe2 come closer with enhanced buckling, breaking the in-plane sublattice symmetry and making the interlayer hopping within different sublattice points asymmetric. This is evident from the appearance of the Dirac cone on high symmetric point in the case of AA stacked CaSi2 and CaGe2 systems, in spite of the fact that, the hopping does not break the sublattice symmetry. The AB staking does not show the shift of the Dirac cone. This can be attributed to the in-plane sublattice symmetry breaking. We further investigate the bilayer sandwich compound of Si2-Ca-Si2 and Ge2-Ca-Ge2. That geometry also indicates the enhanced coupling induced interaction of individual Dirac dispersions that leads to the shift of Dirac cone. This makes the next-next-neighbor electron hopping within A sublattice points and within B sublattice points different. Because of the buckled geometry of silicene/germanene layers, only one kind of sublattice points in each layer is close proximity to the intercalated Ca atoms.

The reactivity of freestanding silicene/germanene towards hydrogenation, fluorination and mixed hydrogenation-fluorination (H-silicene-F and H-germanene-F) was studied by computing their binding energies. The binding energies, Eb are calculated as: 

EB = C ED − EEFB − EGH

(3)

Where ET is the total energy of system, Esub is the energy of substrate, ESi is the energy of silicene and N is the number of Si, Ge, F and H atoms in the unit cells. The binding energy of H-Silicene-F and H-Germanene-F are -4.01 eV and -4.74 eV per HF, respectively. The band gap of freestanding H-silicene-F and H-germanene-F are 1.65 eV and 0.60 eV, respectively. Interestingly, on mixed functionalization of the top-layer of germanene (Hgermanene-F) in CaGe2 multilayer, the Ge-monolayer becomes planar. This is in marked contrast to freestanding H-germanene-F (hydrofluro-germanene) which is strongly buckled. To understand this effect, the total charge associated with each atom and the zero flux surfaces’ defining the Bader volumes was computed as: ∆ρ = ρD − ρEFB − ρGH

(4)

where ρD refers to total charge, ρEFB is the charge of substrate and ρSi, Ge= silicane, germanane, fluorosilicene, fluorogermanene, half hydrogenated and half fluorinated silicene/germanene charge densities, respectively. The 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 30

charge density differences as displayed in Fig. 4(a) highlight that, the charge transfer occurs mainly within the silicene lower lying layer. Fig. 4(b) shows that the fluorine atoms accept 1.09 electrons while the Ca atoms donate 1.97 electrons. Fig. 4(c) indicates the mixed hydrogenated-fluorinated structure of CaSi2 multilayer wherein a charge density difference of 1.24 and 0.96 electrons exist for F and H, respectively. As shown in Fig. 5(a) for hydrogenation of the topmost germanene sheet in CaGe2, hydrogen accepts 0.76 electrons and which are donated by the Ge-atoms. The topmost germanane layer and the CaGe2 substrate show weak van der Waals interaction between the two systems with the binding energy -2.03 eV per H atoms. Such phase-separated state was reported for graphene and graphane to maximize/minimize the boundary length between graphenegraphane.65 Fig. 5(b) refers to the fluorinated CaGe2 multilayer which shows a charge density difference of 1.11, 1.09 and 1.97 electrons for F, Ge and Ca atoms, respectively. Extremely, the lower F atoms interacts with Ca atoms which results in almost doubling of the binding energy (-4.31 eV) vis-a-vis hydrogenated CaGe2. In case of mixed hydrogenated-fluorinated germanene, the Ge-monolayer becomes planar due to the strong binding of F-atoms with Ca(II) ions (-6.47 eV) as shown in Fig. 5(c).

B. Band gaps The band gaps of 2D Si1-xCx and Si1-xGex monolayers are related not only with the substitution ratio of C and Ge but also on the detailed Si-C and Si-Ge order. As observed in Fig. 6(a) for Si1-xCx at a given composition, several possible structures exist and they have different band gaps. This leads to a rich variation in the band gaps when the composition of the hybrid is changed. This is also shown in Fig. 6(a) for which different structure at the same composition creates different gaps and hence, a spread in the band-gap versus composition plots. The band gap first increases with increasing x and reaches a peak of 2.55 eV at x = 0.5 (Si:C = 1:1). Beyond this, the band gap starts to decrease and eventually collapses to 0 for x=0. Interestingly, the band gap of Si1-xCx shows a non-monotonic bell-shaped relationship with doping. For less than 50% doping of C in silicene, the band gap increases with increase in doping then it reaches a maxima at x = 0.5 and eventually again decays to zero for 100% doping. The same band gap trend is observed in Ge1-xCx in Fig. 6(b). In other words, 10 ACS Paragon Plus Environment

Page 11 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

intermediate doping connects the two zero-gap semiconductors pair namely silicene-graphene and germanenegraphene forms new gap opened semiconductors. Our studies of the electronic properties results agree with previous result reported.66 Due to the difference in electronegativity between Si and C atoms (silicon: 1.90 and carbon: 2.55 at the Pauling’s scale), the valence electrons of Si tend to transfer towards the nearest C atoms for 2D Si0.5C0.5 sheet creating charge-transfer states and electron localizations. This is evident from the partial charge density of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) for Si0.5C0.5 (supporting Information Fig. S4). Such localized electron distribution, consequently opens up a band gap of 2.55 eV in Si0.5C0.5 which is also in agreement with the previous reports.67, 68 The SiC and GeC monolayer structure have large band gaps as shown in Fig. 6(a) and 6(b) at x=0.5. The valance band width of the Si1-xCx and Ge1-xCx decreases with the increase in the electronegativity difference between Si/Ge and the C atom. A different situation arises in the case of Si1-xGex for which the band gap remains constant (zero band gap) with increasing x as shown in Table 1.

Obviously, the energy bands of the silicene, graphene and germanene without doping present the expected zerogap character, with π and π* bands both cross the Fermi level at the K-point, as displayed in Fig. 7(I(a-c)), respectively. The Dirac point is found at the symmetric K point in the hexagonal Brillouin zone. Therefore, it is expected that doping will mainly affect the occupation of the π and π* bands. Fig. 7(II(a)) and Fig. 7 (II(b)) shows the band structures of CaSi2 and CaGe2, respectively. This can be understood from the analyses along the high

symmetry

points

of

the

first

Brillouin

zone.

It

is

evident

that

unlike

freestanding

silicene/germanene/graphene which shows a Dirac cone at K-point at the Fermi energy as shown Fig. 7(I (a-c)) respectively the same cannot be detected for CaSi2 and CaGe2. According to our calculations, the two Dirac cones at K and H lie at 2.11 eV and 2.71 eV for CaSi2, respectively below the Fermi energy. In case of CaGe2 Dirac cones found at 2.06 eV at K-point and 2.92 eV at H-point below the Fermi energy. As a result, a large band gap of 4.29 eV and 4.05 eV open up at the K and H points respectively. As a consequence of the electron transfer from Ca atoms to the silicene/germanene layers, the conical dispersion is shifted below the Fermi 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 30

energy. These observations clearly indicate that the buckling silicene/germanene layers plays an important role in introducing sub-lattice symmetry breaking and consequent asymmetric interlayer hopping within different sublattice points. That makes the individual Dirac dispersions of individual buckled layers interact with each other, thus shifting the Dirac cone from the high symmetric points. Their similar behavior have been observed experimentally and theoretically.35, 42, 69

Having established the favorable structures of 2D Si1-xCx, Si1-xGex, Ge1-xCx and Si1-2xCxGex sheets, we proceed to explore their novel properties. We found that the band gaps could be direct or indirect, which also depends on the detailed configurations. For further understanding the change in the band gaps, we plotted the energies band structures for 2D Si1-xCx, Si1-xGex, Ge1-xCx and Si1-2xCxGex monolayer systems, respectively (supporting information Fig.S3 (I-IV)). There are indirect band gaps for x=0.111, x=0.222, and x=0.333 in C doping silicene (Si1-xCx) and for x=0.111, x=333 and x=0.778 in C doping germanene (Ge1-xCx). Similarly, for codoping C-Ge silicene (Si1-2xCxGex) at x=0.278, x=0.333, 0.389 and 0.444 have indirect band gaps. Their indirect and direct band structures are illustrated in (supporting information Table S3). Overall, the energy and type of the band gaps are closely related with the concentration ratio of C atoms and the Si-C and Ge-C order. For small doping (x=0.111) in Si1-xCx, the valence band maximum (VBM) and conduction band minima (CBM) almost meet at the Fermi energy at ᴦ-point. As doping increases, both VB and CB flatten out considerably to create the largest band gap for Si0.5C0.5 in the Si1-xCx family. On further increase in the mole fraction of C atoms in silicene, the CBM starts to come closer to the Fermi-energy and becomes almost a semi-metal at x=0.889. The same properties observed for Ge1-xCx monolayer. In case of Si1-xGex for all doping concentrations, VBM and CBM are always close to the Fermi level which results in direct and negligible band gap (Eg= 0.02 eV) at the ᴦpoint. When C and Ge are co-doped within silicene, the band-gaps increase as we increase the fraction of C-Ge. For x=0.06 the VB and CB touch the Fermi level creating semimetalic state. As we increase the doping ratio of x (C-Ge) the band gap starts to open and both VB and CB are far from the Fermi level. Since Ge-doping alone does not affect the band-gap in silicene, the remarkable opening of the band gaps on co-doping C arises entirely 12 ACS Paragon Plus Environment

Page 13 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

from carbon atoms. As seen from Table 1, Si1-xCx and Ge1-xCx shows the discussed bell-shaped behavior, Si1xGex

shows zero band gaps and Si1-2xCxGex shows an energy gap ranging from 0.04 eV to 2.07 eV depending

upon the doping and their ratio. Generally, silicene/germanene monolayer doping with C atom concentration open up the band gap. This can be understood from the band gap flexible by increasing or decreasing the concentration ratio. Thus, evidences are the monolayer silicene and germanene doping with C atom at different concentration ratio increases (from x=0.111 to x=0.5) and decreases (from x=0.5 to x=1) the band gap of the monolayer. On the other hand, hydrogenation and fluorination freestanding silicene/germanene obtained open up band gap. In opposition to, silicene doping with Ge atom and silicene/germanene intercalated Ca layers functionalized with H, F and mixed H and F together with bilayer sandwiched and multilayer does not open up the band gap. The energy shift of Dirac dispersion below the Fermi energy in CaSi2 and CaGe2 in AA stacking sequence systems are expected due to electron transfer from Ca atom to silicene/germanene layers (supporting information S5).

C. Optical properties To explain origin of the peaks, contributions of specific transitions between the bands along the high-symmetry ᴦ-K-M-ᴦ direction as displayed in Fig. 7(I (a-c)) were calculated. In the optical spectra, one can expect only π to π* and σ to σ* transitions as allowed if the light is polarized parallel to the silicene/germanene/doping-silicene monolayer. In contrast, only π to σ* and σ to π* transitions are allowed if the light is polarized perpendicular to the sheets. The absorption spectrum of pure silicene, graphene and germanene are shown in Fig. 7(I (d-f)). We found that silicene due to the low-buckled structure exhibits two prominent absorption peaks: the π to π* at photon energy 1.65 eV and 4 eV for E⊥c and 6.5 eV for the polarization of light E∥c as shown in Fig. 7(I(d)). The absorption spectra of graphene the π to π* transition along the K-M direction was found from two peaks at 2.0 eV and 4.0 eV for the polarization of light E⊥c. In the energy range from 7 eV to 11 eV the imaginary part of dielectric function is negligible in the polarization of light E∥c. The main feature in the region from 12 to 24 eV is an intense resonance at 14 eV which arises from σ to σ* transitions around the M points in both ᴦ-M and 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 30

M-ᴦ direction. For germanene, the π to π* absorption spectra was found to be near-infrared at 0.56 eV along the K-M direction when the light is polarized along E⊥c. The σ to σ* transitions are around the ᴦ-M and M-ᴦ points at energy 3.85 eV in the near UV when the polarization of light E∥c as shown in Fig. 7(I(e)). The excitonic effects in 2D silicene and germanene are significant and stronger than that in graphene. Those are due to the fact that the Si-Si and Ge-Ge bond are weaker than the C-C bond. The results found were in agreement with previous calculations.70

In Fig. 7(II (c and d)) we list all bound excited states in the absorption spectra of CaSi2 and CaGe2, respectively. The top of σ-band is situated closer to Fermi energy than the crossing point of π and π* bands at the K and H points. This suggests that electrons from Ca atoms are pre-dominantly donated to π* band in CaSi2 and CaGe2 with a maximum absorption frequency of 3.15 eV and 2.75 eV, respectively. The first exciton peak of CaSi2 for E⊥c was found to be located at 0.5 eV, the second exciton peak lies at 1.25 eV for E⊥c and the third exciton peak exists at 2 eV. The most intense peak in the absorption spectrum exists at 2.75 eV for CaSi2 which begins in the near-infrared region and gradually increases to the visible frequency as shown in Fig. 7(II(c)). More interestingly, the first exciton is observed at 1.85 eV for E∥c in CaGe2 and an intense peak exists at photon energy of 2.85 eV as shown in Fig. 7(II(d)). The optical absorption spectra of freestanding silicene and germanene functionalized with H, F and mixed H-F were also calculated (supporting information Fig. S12 and Fig. S13). The optical absorption shows maxima at 3.82 eV for E⊥c and 7.35 eV for E∥c in silicane. Vertical transitions from the top of the valance band to the bottom of the conduction band close to ᴦ-point contribute towards the absorption spectrum. For germanene, the optical absorption spectra begins in the near infrared region (0.95 eV) and is observed in visible region at photon energy 3.36 eV for E⊥c and peaks at 7.46 eV for E∥c. In case fluorine functionalized silicene, the optical absorption spectrum peaks a 4.00 eV and 8.75 eV for E⊥c and E∥c, respectively. For F-germanene, the optical absorption spectrum shows maxima 3.04 eV and 7.12 eV for E⊥c and E∥c, respectively. For H-Silicene-F and H-Germanene-F, the maximum exciton peak lies around 3.32 eV and 3.01 eV for E⊥c, respectively. In general, for freestanding silicene functionalized with H 14 ACS Paragon Plus Environment

Page 15 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

and F, the absorption occurs mostly in the visible region while that for germanene it is in the infrared for E∥c. However, when E∥c both silicene/germanene in their pristine together with functionalized forms exhibit UV absorption.

The typical optical spectra of C, Ge and C-Ge doping silicene are shown in Fig. 8(I-IV). With increase of the C doping, the intensity of the main absorption peak (E⊥c) increases in the range of 2.50-5.50 eV and then decrease after x=0.556-0.778 as shown in Fig. 8(I(a-f)). For polarization along E∥c as we increases the Cdoping concentration the intensity of the absorption peak increases and at concentration x=0.778 this peak disappears. As shown in Fig. 8(I-III(a-f)) when the light is polarized E∥c the photon energy increases in the UV region for silicene doped with C/Ge and germanene doped with C. In other cases of mixed C-Ge doped system as shown in Fig. 8(IV(a-f)), with increases of doping concentration, the intense peak vary when the light is polarized E⊥c and for E∥c, the absorption peak is found at 8 eV for x=0.06 as shown in Fig. 8(IV(a)). Therefore, in the presence of C-Ge atoms in the system gives rise in-plane low energy π excitation peak due to polarization along the 2D-plane. For C-doping systems in both silicene and germanene the absorption is redshifted with respect to pristine silicene and graphene. For E∥c polarization the absorption spectrum is blueshifted compared to the pristine systems irrespective of doping.

Conclusion In summary, we have performed a comprehensive structural search on the two-dimensional Si1-xCx, Si1-xGex, Ge1-xCx and Si1-2xCxGex sheet with 0 ≤ x ≤ 1 within density functional theory. The structures of the doping systems are linearly correlated with the structures of their individual atomic constituents as predicted by the Vegard’s law. All the in-plane hybrid structures are semiconductors with widely tunable band gaps on doping. The band gaps increase from 0 eV at x=0 or 1 to 2.55 eV at x=0.5 for 2D Si1-xCx. Furthermore, we obtain the band gap 2.07 eV for monolayer Ge1-xCx at x=0.5. The band gaps for C-Ge co-doping systems are controlled by C atoms as Ge doping alone is insufficient to open up the band gap. Optical properties of 2D Si1-xCx, Si1-xGex, 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 30

Ge1-xCx and Si1-2xCxGex hybrids indicates that with different C, Ge and C-Ge composition ratios the optical responses can be remarkably different and yet absorb visible light for specific compositions. The optical spectrum of silicene/graphene/germanene/doping-silicene monolayer sheets have been calculated for light polarization parallel and perpendicular to the plane. When C-doping in silicene is small (x=0.111 to x=0.5), the main absorption peak occurs in the range of 2.50-5.50 eV and almost vanishes at higher doping. For Ge-doping silicene, the intensity of the main absorption peak decreases in the visible region from 3.75-0.75 eV. Amongst the four systems studied here (three binary and one ternary), Si1-xCx is expected to be most useful to experimentalist as one might be able to access them by sequential replacement of Si by C or vice-versa and observe a monotonic band gap dependence till x ≤ 0.5. It is important to note that for the same stoichiometric ratio of Si/C, there can be multiple low-energy configurations which would result in a spread across the bellshaped curve of the band-gap with respect to doping. Presence, of such configurations would be all the more important under ambient experimental conditions. Nevertheless, it would be worthwhile to experimentally explore the predicted bell-shaped band gap behavior for C doping silicene/germanene and optical properties of such low-dimensional materials. Clearly controlled C doping of silicene and germanene monolayer opens up the band-gaps and also provide a tool to harvest solar energy for semiconductor and photovoltaic applications. Interestingly, apart from pristine monolayer silicene, Ca(II) intercalated multilayer silicene/germanene to the electronic properties and stability of their hydrogenated and fluorinated as available in the van der Waals mineral phases of CaSi2 and CaGe2 are also investigated. Based on phonon spectra, we showed that silicene, germanene, graphene, CaSi2 and CaGe2 are stable as it does not exhibit imaginary phonon mode. Ca(II) intercalated multilayer silicene/germanene towards hydrogenation, fluorination and mixed hydrogenationfluorination are studied. Optical absorption calculations demonstrate that doping monolayers along with Ca(II) intercalated multilayer silicene/germanene have strong photo-absorption in visible light region. The presence of Ca(II) stabilizes the silicene/germanene layers and also the Ca(II) is affecting the electronic structure of these 2D materials.

16 ACS Paragon Plus Environment

Page 17 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ASSOCIATED CONTENT Supporting Information. AUTHOR INFORMATION Corresponding Author E-mail: [email protected]. Tel.: +91-33-24734971. Notes The authors declare no competing financial interest.

Acknowledgements This work was supported by The World Academy of Sciences (TWAS) - for the advancement of science in developing countries (FR number: 3240280472) and Indian Association for the Cultivation of Science (IACS) is gratefully acknowledged. AD thanks DST, INSA and BRNS for partial funding.

References 1. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, A. A. Firsov,Two-dimensional Gas of Massless Dirac Fermions in Graphene, Nature 2005, 438, 197-200. 2. C. Berger, Z. Song, X. Li, X.Wu, N. Brown, C. Naud, D. Mayou, T. Li, J. Hass, A. N. Marchenkov, E. H. Conrad, P. N. First, W. A. de Heer, Electronic Confinement and Coherence in Patterned Epitaxial Graphene, Science 2006, 312, 1191-1196. 3. L. Zhang, Y. Zhang, J. Camacho, M. Khodas, I. Zaliznyak, The Experimental Observation of Quantum Hall Effect of l=3 Chiral Quasiparticles in Trilayer Graphene, Nature Phys. 2011, 7, 953-957. 4.

M. J. Allen, V. C. Tung, R. B. Kaner, Honeycomb Carbon: A Review of Graphene, Chem. Rev. 2010, 110, 132-145 .

5. G. Gao, N. W. Ashcroft, R. Hoffmann, The Unusual and The Expected in The Si/C Phase Diagram, J. Am. Chem. Soc. 2013, 135, 11651-11656. 6. D. Jose, A. Datta, Structures and Chemical Properties of Silicene: Unlike Graphene, Acc. Chem. Res. 2014, 47, 593-602. 7. T. K. Mandal, D. Jose, A. Nijamudheen, Structures and Electronic Properties of Heavier Congeners of Disk-Like Molecules: (Si, Ge) Sulflower and (Si, Ge) Olympicene, A. Datta, J. Phys. Chem. C 2014, 118, 12115 -12120. 8.

D. Jose, A. Nijamudheen, A. Datta, Tip Enhanced Raman spectroscopy (TERS) as a Probe for The Buckling Distortion in Silicene, Phys. Chem. Chem. Phys. 2013, 15, 8700-8704.

9. X.Q.Wang, H.D.Li, J. T. Wang, Induced Ferromagnetism in one-side Semihydrogenated Silicene and Germanene, Phys. Chem. Chem. Phys. 2012, 14, 3031-3036. 17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 30

10. C. Zhang, S. Yan, First-Principles Study of Ferromagnetism in Two-Dimensional Silicene with Hydrogenation, J. Phys. Chem. C 2012, 116, 4163-4166. 11. F. Zheng, C. Zhang, The Electronic and Magnetic Properties of Functionalized Silicene: a First-Principles Study, Nano. Res. Lett. 2012, 7, 422-5. 12. C. C. Liu, W. Feng, Y. Yao, Quantum Spin Hall Effect in Silicene and Two-Dimensional Germanium, Phys. Rev. Lett. 2011, 107, 076802-076806. 13. S. Cahangirov, M. Topsakal, E. Akturk, H. Sahin, S. Ciraci, Two- and One-Dimensional Honeycomb Structures of Silicon and Germanium, Phys. Rev. Lett. 2009, 102, 236804-4. 14. A. Nijamudheen, R. Bhattacharjee, S. Choudhury, Ayan Datta, Electronic and Chemical Properties of Germanene: The Crucial Role of Buckling, J. Phys. Chem. C 2015, 119, 3802-3809. 15. M. E. Davila, L. Xian, S. Cahangirov, A. Rubio, and G. Le. Lay, Germanene: a Novel Two-dimensional Germanium Allotrope akin to Graphene and Silicene, New J. Phys. 2014, 16, 095002-10. 16. G. R. Bhimanapati, Z. Lin, V. Meunier, Y. Jung, J. Cha, S. Das, Di Xiao, Y. Son, M. S. Strano, V. R. Cooper, L. Liang, S. G. Louie, E. Ringe, Wu Zhou, S. S. Kim, R. R. Naik, B. G. Sumpter, H. Terrones, F. Xia, Y. Wang, J. Zhu, D. Akinwande, N. Alem, J. A. Schuller, R. E. Schaak, M. Terrones, J. A. Robinson, Recent Advances in Two-Dimensional Materials Beyond Graphene, ACS Nano 2015, 9, 11509-11539. 17. L. Matthes, O. Pulci, F. Bechstedt, Massive Dirac Quasiparticles in the Optical Absorbance of Graphene, Silicene, Germanene, and Tinene, J. Phys.: Condens. Matter 2013, 25, 395305-395312. 18. Y. Ding, Y. Wang, Density Functional Theory Study of The Silicene-like SiX and XSi3 (X = B, C, N, Al, P) Honeycomb Lattices: The Various Buckled Structures and Versatile Electronic Properties, J. Phys. Chem. C 2013, 117, 18266-18278. 19. Z. Ni, Qi. Liu, K. Tang, J. Zheng, J. Zhou, R. Qin, Z. Gao, D. Yu, J. Lu, Tunable Band gap in Silicene and Germanene, Nano Lett. 2012, 12, 113-118. 20. H. Sahin, S. Cahangirov, M. Topsakal, E. Bekaroglu, E. Akturk, R. T. Senger, S. Ciraci, Monolayer Honeycomb Structures of Group-IV Elements and III-V Binary Compounds: First-Principles Calculations, Phys. Rev. B 2009, 80, 155453-155465. 21. L. C. Lew Yan Voon, E. Sandberg, R. S. Aga, A. A. Farajian, Hydrogen Compounds of Group-IV Nanosheets, Appl. Phys. Lett. 2010, 97, 163114-163123. 22. M. Houssa, E.Scalise, K. Sankaran, G. Pourtois, V. V. Afanas’ev, A. Stesmans, Electronic Properties of Hydrogenated Silicene and germanene, Appl. Phys. Lett. 2011, 98, 223107-223110. 23. N. Gao, W. T. Zheng, Q. Jinag, Density Functional Theory Calculations for Two-Dimensional Silicene with Halogen Functionalization, Phys. Chem. Chem. Phys. 2012, 14, 257-261.

18 ACS Paragon Plus Environment

Page 19 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

24. C. Chowdhury, S. Jahiruddin, A. Datta, Pseudo- JahnTeller Distortion in Two-Dimensional Phosphorus: Origin of Black and Blue Phases of Phosphorene and Band Gap Modulation by Molecular Charge Transfer, J. Phys. Chem. Lett. 2016, 7, 1288-1297. 25. G. R. Berdiyorov, M. Neek-Amal, F.M. Peeters, Adri C. T. van Duin, Stabized Silicene within Bilayer Graphene: A Proposal Based on Molecular Dynamics and Density-Functional Tight-binding Calculation, Phys. Rev. B 2014, 89, 024107-6 26. Pablo A. Denis, Stacking Functionalized Silicene: a Powerful System to Adjust the Electronic Structure of Silicene, Phys. Chem. chem. Phys. 2015, 17, 5393-5402. 27. J. Zhu, A. Chroneos, UdoSchwingenschögl, Silicene/ Germanene on MgX2 (X = Cl, Br, and I) for Li-ion Battery Applications, Nanoscale 2016, 8, 7272-7277. 28. P. Vogt, P. De Padova, C. Quaresima, J. Avila, E. Frantzeskakis, M. C. Asensio, A. Resta, B. Ealet, G. Le Lay, Silicene: Compelling Experimental Evidence for Graphene like Two-Dimensional Silicon, Phys. Rev. Lett. 2012, 108, 155501-155506. 29. P. De Padova, O. Kubo, B. Olivieri, Multilayer Silicene Nanoribbons, C. Quaresima, T. Nakayama, M. Aono, G. Le Lay, Nano Lett. 2012, 12, 5500-5503. 30. Dávila, M. E. Le Lay G. Few Layer Epitaxial Germanene: a Novel Two-Dimensional Dirac Material, Sci. Rep. 2016, 6, 20714-6: doi: 10.1038/srep20714. 31. C. C. Liu, W. Feng, Y. Yao, Quantum Spin Hall Effect in Silicene and Two-Dimensional Germanium, Phys. Rev. Lett. 2011, 107, 076802-076806. 32. M. Ezawa, Valley-Polarized Metals and Quantum Anomalous Hall Effect in Silicene, Phys. Rev. Lett. 2012, 109, 055502-055507. 33. S. Dutta, K. Wakabayashi, Spin and Charge Excitations in Zigzag Honeycomb Nanoribbons: Effect of Many Body Correlations, Jpn. J. Appl. Phys. 2014, 53, 06JD01- 4. 34. S. Cahangirov, M. Audiffred, P. Tang, A. Iacomino, W. Duan, G. Merino, A. Rubio, Electronic Structure of Silicene on Ag(111): Strong Hybridization Effects, Phys. Rev. B 2013, 88, 035432-035438. 35. S. M. Pratik, A. Nijamudheen, A. Datta, Topochemical Transformations of CaX2 (X=C, Si, Ge) to Form Free-Standing Two-Dimensional Materials, Chem. Eur. J. 2015, 21, 18454-18460. 36. B. Huang, H. Deng, H. Lee, M. Yoon, B. G. Sumpter, F. Liu, S. C. Smith, Su-Huai Wei, Exceptional Optoelectronic Properties of Hydrogenated Bilayer Silicene, Phy. Rev. X 2014, 4, 021029-12. 37. M. S. Dresselhaus, G. Dresselhaus, Intercalation Compounds of Graphite, Adv. Phys. 2002, 51, 1-186. 38. S.Ichinokuta, K. Sugawara, A. Takayama, T. Takahashi, Shuji Hasegawa, Superconducting CalciumIntercalated Bilayer Graphene, ACS Nano 2016, 10, 2761-2765.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 30

39. K. Kanetani, K. Sugawara, T. Sato, R. Shimizu, K. Iwaya, T. Hitosugi, and T. Takahashi, Ca Intercalated Bilayer Graphene as a Thinnest Limit of Superconducting C6Ca, Proc. Natl. Acad. Sci. 2012, 109, 1961019613. 40. H. Liu, N. Han, J. Zhao, Band Gap Opening in Bilayer Silicene by Alkali Metal Intercalation, J. Phys.: Condens. Matter 2014, 26, 475303-475310. 41. R. Quhe, Y. Yuan, J. Zheng, Y. Wang, Z. Ni, J. Shi, D. Yu, J. Yang, J. Lu, Does The Dirac Cone Exist in Silicene on Metal Substrates?, Sci. Rep., 5476; DOI:10.1038/srep05476 2014, 4. 42. Fleurence, A.; Friedlein, R.; Ozaki, T.; Kawai, H.; Wang, Y.; Yamada-Takamura, Y., Experimental Evidence for Epitaxial Silicene on Diboride Thin Films, Phys. Rev. Lett. 2012, 108, 245501-245506. 43. B. Feng, Z. Ding, Sh. Meng, Y. Yao, X. He, P. Cheng, L. Chen, K. Wu, Evidence of Silicene in Honeycomb Structures of Silicon on Ag(111), Nano Lett. 2012, 12, 3507-3511. 44. Tsoutsou, D.; Xenogiannopoulou, E.; Golias, E.; Tsipas, P.; Dimoulas, A., Evidence for Hybrid Surface Metallic Band in (4 x 4) Silicene on Ag(111), Appl. Phys. Lett. 2013, 103, 231604-231609. 45. H. Jamgotchian, Y. Colignon, N. Hamzaoui, B. Ealet, J. Y. Hoarau, B. Aufray, J. P. Biberian, Growth of Silicone layers on Ag(111): Unexpected Effect of The Substrate Temperature, J. Phys.: Condens. Matter 2012, 24, 172001-172009. 46. Yi. Ding, Y. Wang, Electronic Structures of Silicene/GaS Heterosheets, Appl. Phys. Lett. 2013, 03, 043114-043118. 47. Houssa, M.; van den Broek, B.; Scalise, E.; Pourtois, G.; Afanasev, V. V.; Stesmans, A., An Electric Field Tunable Energy Band Gap at Silicene/(0001) ZnS Interfaces, Phys. Chem. Chem. Phys. 2013, 15, 37023705. 48. A. Bhattacharya, S. Bhattacharya, G. P. Das, Exploring Semiconductor Substrates for Silicene Epitaxy, Appl. Phys. Lett. 2013, 103, 123113-123117. 49. J. Zhu, U. Schwingenschlögl, Stability and Electronic Properties of Silicene on WSe2 J. Mater. Chem. C 2015, 3, 3946-3953. 50. K. Yang, S. Cahangirov, A. Cantarero, A. Rubio, R. DAgosta, Thermoelectric Properties of Atomically Thin Silicene and Germanene Nanostructures, Phys. Rev. B 2014, 89, 125403-125416. 51. K. Zberecki, M. Wierzbicki, J. Barnas, R. Swirkowicz, Thermoelectric Effects in Silicene Nanoribbons, Phys. Rev. B 2013, 88, 115404. 52. J. Zhu, UdoSchwingenschögl, Structural and Electronic Properties of Silicene on MgX2 (X = Cl, Br, and I), ACS Appl. Mater. Interface 2014, 6, 11675-11681. 53. E. Noguchi, K. Sugawara, R. Yaokawa, T. Hitosugi, H. Nakano, T. Takahashi, Direct Observation of Dirac Cone in Multilayer Silicene Intercalation Compound CaSi2, Adv. Mater. 2015, 27, 856-860. 20 ACS Paragon Plus Environment

Page 21 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

54. J. P. Perdew, K. Burke, M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 1996, 77, 3865-3869. 55. G. Onida, L. Reining, A. Rubio, Electronic Excitations: Density-Functional Versus Many-Body GreensFunction Approaches, Rev. Mod. Phys. 2002, 74, 601-659. 56. B. Adolph, J. Furthmüller, F. Bechstedt, Optical Properties of Semiconductors using Projector-Augmented Waves, Phys. Rev. B 2001, 63, 125108-125116. 57. A.V. Rozhkov, A.O. Sboychakov, A.L Rakhmanov, Franco Nori, Electronic Properties of Graphene-Based Bilayer Systems, Phys. Reports, 2016, 648, 1-104. 58. M. Gajdos, K. Hummer, G. Kresse, J. Furthmuller, F. Bechstedt, Linear Optical Properties in The Projector Augmented Wave Methodology, Phys. Rev. B 2006, 73, 045112-045121. 59. D. Jana, C.-L. Sun, L.-C. Chen, K.-H. Chen, Effect of Chemical Doping of Boron and Nitrogen on the Electronic, Optical and Electrochemical Properties of Carbon Nanotubes, Prog. Mater. Sci. 2013, 58, 565635. 60. S. Brodersen, D. Lukas, W. Schattke, Calculation of The Dielectric Function in a Local Representation, Phys. Rev. B 2002, 6, 085111-8. 61. N. Wiser, Dielectric Constant with Local Field Effects Included, Phys. Rev. B 1963, 129, 62-9. 62. A. R. Denton, N. W. Ashcroft, Vegard’s law, Phys. Rev. A 1991, 43, 3161-3164. 63. S. Mukherjee, A. Nag, V. Kocevski, P. K. Santra, M. Balasubramanian, S. Chattopadhyay, T. Shibata, F. Schaefers, J. Rusz, C. Gerard, O. Eriksson, C. U. Segre, D. D. Sarma, Microscopic Description of The Evolution of The Local Structure and an Evaluation of The Chemical Pressure Concept in a Solid Solution, Phys. Rev. B 2014, 89, 224105-11. 64. T. K. Mukhopadhyay, A. Datta, Deciphering the Role of Solvents in the Liquid Phase Exfoliation of Hexagonal Boron Nitride: A Molecular Dynamics Simulation Study, J. Phys. Chem. C. 2017, 121, 811822. 65. A.L Rakhmanov, A.V. Rozhkov, A.V. Rozhkov, Franco Nori, Phase Separation of Hydrogen Atoms adsorbed on Graphene and The Smoothness of The Graphene-Graphane Interface, Phys. Rev. B, 2012, 85, 035408-6. 66. Z. Shi, Z. Zhang, A. Kutana, B. I. Yakobson, Predicting Two-Dimensional Silicon Carbide Monolayers, ACS Nano 2015, 9, 9802-9809. 67. G. Liu, M. S. Wu, C. Y. Ouyang, B. Xu, Structural and Electronic Evolution from SiC sheet to Silicene, Int. J. Mod. Phys. B 2013, 27, 1350188-9.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 30

68. Tie-Yu Lü, Xia-Xia Liao, Hui-Qiong Wang, Jin-Cheng Zheng, Tuning the indirect-direct band gap transition of SiC, GeC and SnC Monolayer in a Graphene-like Honeycomb Structure by Strain Engineering: a Quasiparticle GW Study, J. mater. Chem., 2012, 22, 10062-10068. 69. D. Jose, A. Datta, Understanding of The Buckling Distortions in Silicene, J. Phys. Chem. C. 2012, 116, 24639-24648. 70. W. Wei, Y. Dai, B. Huang, T. Jacob, Many-Body Effects in Silicene, Silicane, Germanene and Germanane, Phys. Chem. Chem. Phys. 2013, 15, 8789-8794.

22 ACS Paragon Plus Environment

Page 23 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1: Phonon dispersion of (a) Silicene, (b) Germanene, (c) Graphene, (d) CaSi2 and (e) CaGe2. ᴦ (0, 0, 0), K (2/3, 1/3, 0), M (½, 0, 0) refer to special points in the first Brillouin zone in reciprocal space.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 30

Figure 2: Lowest energy structures of the most stable structures of: (I) 2D Si1-xCx at different doped ratio of C atoms (a) silicene (x=0.0), (b) x=0.111, (c) x=0.222, (d) x=0.333, (e) x=0.5, (f) x=0.556, (g) x=0.778) (h) x=0.889 and (i) graphene (x=1.0). (II) 2D Si1-xGex at different doped ratio of Ge atoms (a) silicene (x=0.0), (b) x=0.111, (c) x=0.222, (d) x=0.333, (e) x=0.5, (f) x=0.556, (g) x=0.778) (h) x=0.889 and (i) germanene (x=1.0). (III) 2D Ge1-xCx at different doped ratio of C atoms (a) germanene (x=0.0), (b) x=0.111, (c) x=0.222, (d) x=0.333, (e) x=0.5, (f) x=0.556, (g) x=0.778) (h) x=0.889 and (i) graphene (x=1.0). (IV) 2D Si1-2xCxGex at different co-doped ratio of C and Ge atoms (a) silicene (x=0.0), (b) x=0.056, (c) x=0.111, (d) x=0.167, (e) x=0.278, (f) x=0.333, (g) x=0.389, (h) x=0.444 and (i) x=0.5. Blue spheres represent Si atoms, dark gray spheres represent Ge atoms and brown spheres represent C atoms.

24 ACS Paragon Plus Environment

Page 25 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3: Lattice parameters and bond lengths (both in Å) for (a) Si1-xCx, (b) Si-C, (c) Si1-xGex, (d) Si-Ge (e) Ge1-xCx and (f) Ge-C at different fractions of C and Ge.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 30

Figure 4: Charge density difference for (a) Hydrogenation of CaSi2, (b) Fluorination of CaSi2 and (c) Mixed hydrogenation-fluorination of CaSi2. Note that all the reactions occur on the topmost silicene layers of CaSi2. Ca, Si, F and H atoms are shown in green, blue, gray and white respectively. Pink/yellow represents charge accumulation/depletion, the isosurfaces refer to isovalues=4.85 x10-3 electrons/bohr3.

Figure 5: Charge density difference for (a) Hydrogenation of CaGe2, (b) Fluorination of CaGe2 and (c) Mixed hydrogenation-fluorination of CaGe2. Note that all the reactions occur on the topmost germanene layers of CaGe2. Ca, Ge, F and H atoms are shown in green, violet, gray and white respectively. Pink/yellow represents charge accumulation/depletion, the isosurfaces refer to isovalues=4.85 x10-3 electrons/bohr3.

26 ACS Paragon Plus Environment

Page 27 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 6: Calculated band gaps in (a) Si1-xCx and (b) Ge1-xCx sheets as the function of doping ratio of C atoms. Note that the spread in the band gaps for intermediate doping arises due to multiple configurations/structures at the same composition. The solid lines connect the band gaps for the most stable configuration at a given doping concentration.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 30

(I) Silicene-Graphene-Germanene

(II) CaSi2 and CaGe2

Figure 7: (I) Band structure of pure (a) silicene, (b) graphene and (c) germanene and the imaginary part of dielectric function of optical spectra (d) silicene, (e) graphene and (f) germanene. (II) The band structures along the high symmetric points of (a) the intercalated compound CaSi2 with alternatively stacked silicene and Ca in AA stacking sequence, (b) CaGe2 with alternatively stacked germanene and Ca in AA stacking sequence, (c) Optical absorption spectrum of CaSi2 and (d) CaGe2 for light polarization E⊥c and E∥c to the surface plane shown in blue and red respectively. 28 ACS Paragon Plus Environment

Page 29 of 30

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(II) Si1-xGex

(I) Si1-xCx

(III) Ge1-xCx

(IV) Si1-2xCxGex

Figure 8: Imaginary part of dielectric function of 2D for E⊥c and E∥c in (I) 2D Si1-xCx at different doping ratio of C atoms (a) x=0.111, (b) x=0.222, (c) x=0.333, (d) x=0.5, (e) x=0.556 and (f) x=0.778). (II) 2D Si1-xGex at different doping ratio of Ge atoms (a) x=0.111, (b) x=0.222, (c) x=0.333, (d) x=0.5, (e) x=0.556 and (f) x= 0.778). (III) Ge1-xCx at different doping ratio of C atoms (a) x=0.111, (b) x=0.222, (c) x=0.333, (d) x=0.5, (e) x = 0.556 and (f) x=0.778). (IV) 2D Si1-2xCxGex at different co-doping ratio of C and Ge atoms (a) x=0.06, (b) x = 0.111, (c) x=0.167, (d) x=0.25, (e) x=0.278 and (f) x=0.389 monolayer system.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table of Contents

30 ACS Paragon Plus Environment

Page 30 of 30