arxiv: v1 [cond-mat.mes-hall] 5 Feb 2018

Transcription

1 Semimetal behavior of bilayer stanene I. Evazzade a,, M. R. Roknabadi a, T. Morshedloo b, M. Modarresi a, Y. Mogulkoc c,d, H. Nematifar a a Department of Physics, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran b Nano Structured Coatings Institute of Yazd Payame Noor University, , Yazd, Iran c Department of Engineering Physics, Ankara University, 6, Tandogan, Ankara, Turkey d Institute for Molecules and Materials, Radboud University, Heyendaalseweg 35, 655AJ Nijmegen,The Netherlands arxiv:8.35v [cond-mat.mes-hall] 5 Feb 8 Abstract Stanene is a two-dimensional (D) buckled honeycomb structure which has been studied recently owing to its promising electronic properties for potential electronic and spintronic applications in nanodevices. In this article we present a first-principles study of electronic properties of fluorinated bilayer stanene. The effect of tensile strain, intrinsic spin-orbit and van der Waals interactions are considered within the framework of density functional theory. The electronic band structure shows a very small overlap between valence and conduction bands at theγpoint which is a characteristic of semimetal in fluorinated bilayer stanene. A relatively high value of tensile strain is needed to open an energy band gap in the electronic band structure and the parity analysis reveals that the strained nanostructure is a trivial insulator. According to our results, despite the monolayer fluorinated stanene, the bilayer one is not an appropriate candidate for topological insulator. Keywords: bilayer stanene, semimetal behavior, topological insulators, first-principles calculations. Introduction Recently theoretical investigation of materials based on group IV elements attract considerable attention [,, 3, 4, 5]. Among them, Stanene the two-dimensional (D) tin as a topological insulators has been studied thoroughly because of being ideal candidates to quantum spin Hall (QSH) effect [6, 7, 8, 9,, ]. The topological state was observed experimentally in HgTe quantum wells [], three dimensional structures e.g., Bi x S b x [3], S nte [4], Bi Te 3 [5], TiBiS e [6] and predicted theoretically for two dimensional nanostructures [7, 8, 9,,,, ]. These materials are associated with topological invariant which is different from ordinary insulators having topologically protected gapless edge states that are robust against time-reversal-invariant weak disorders [, 3]. In spite of the fact that D materials have been on the focus of a large number of scientific researches, not all of them are fully utilized as QSH insulators. Graphene has been predicted to be a QSH insulator by Kane and Mele at low temperatures with a too small energy gap ( 3 mev) which is opened by spin-orbit coupling at the Dirac point []. The functionality of QSH insulators is having an experimentally accessible gap at room temperature that is due to large spin-orbit coupling strength [4]. In the pursuit of large-gap QSH insulators, many novel structures have been predicted [, 7, 8, 9, 4, 4, 5] including slightly buckled honeycomb lattice of tin atoms[6], also known as stanene, which recently has been fabricated by molecular beam epitaxy in laboratory [7]. Due to the relatively large atomic mass of tin atoms the intrinsic spin-orbit interaction opens an energy gap in the electronic structure of stanene which is around 7 mev [4, 8]. Theoretical calculations predict a larger gap for halogenated monolayer [4] and ethynyle-derivative [4] stanene, as well as chemically decorated plumbene [5], that suggest these structures for application in nanodevices based on topological insulators. Indeed, it was shown that oxygen adsorption may cause topological behavior to trivial phase transition in stanene nanoribbons [5]. Corresponding auther address: i.evazzade@mail.um.ac.ir (I. Evazzade) Preprint submitted to Physica E February 6, 8

2 Figure : (a) Top and (b) side view of atomic structure of FBSn From a realistic point of view, the synthesis of an exact monolayer of D structure on a substrate is an experimental challenge for experimental scientists. In the case of stanene, the latest attempt by molecular beam epitaxy leads to few layer stanene [7]. The main question arises from the effect of interaction between different layers of stanene on the topological phase of the nanostructure. Based on the aforementioned motivation, we investigate the effect of fluorine absorption on the electronic and topological phase of bilayer stanene. The fluorinated bilayer stanene (FBSn) is studied by using density function theory (DFT) with considering the van der Waals (vdw) interaction between different layers of structure. Owing to the heavy atomic mass of tin atoms, fully relativistic pseudopotential together with spin-orbit coupling is considered for DFT calculations.. Model and Method For our calculations we have used the first principles density functional calculations as implemented in the Quantum Espresso code [9]. The electronic calculations have been performed within the Generalized Gradient Approximation (GGA) and the Perdew-Burke-Ernzerhof (PBE) [3] exchange correlation is adopted. The semi-empirical dispersion interactions (DFT-D) model [3, 3] has been employed in order to consider the van der Waals interaction between adjacent layers of bilayer stanene. To avoid interaction between stanene layers in neighbor unit-cells, a vacuum layer 5Å is set perpendicular to the stanene plane. A Monkhorst and Pack [33] k-point mesh is adopted to sample the Brillouin zone and the cut-off energy for plane wave expansion is set to 7Ry in all calculations. Throughout geometry optimization, cell parameters as well as atomic positions are fully relaxed until the force on atoms is less than.ev/å. We consider the hexagonal unit cell for bilayer fluorinated stanene as shown in Figure. The side view represents the buckling structure of tin atoms. The relaxation process does not break the space inversion symmetry of bilayer fluorinated stanene. The topological invariant determines the topological/trivial phase of structure. According to the Fu-Kane work [34] in a crystal with space inversion symmetry one can obtain the topological invariant by using the parity of occupied electronic bands of Time Reversal Invariant Momenta (TRIM) points [35, 36]. In the hexagonal structure, there are 4 TRIM points (Γand 3 M points) in the first Brillouin zone [7, 9]. 3. Results We consider the AB stacked of bilayer stanene. Figure represents the top and side view of buckled structure of FBSn. Both figures are plotted using XCrySDen program [37]. The lattice parameter of hexagonal unit cell of FBSn is 4.75Å which is smaller than iodinated stanene [7]. In the FBSn the buckling parameter of layers is.8å which is

3 5 frequency [/cm] Phonon mode number Figure : Phonon frequency of bilayer stanene at the Γ point. around pure stanene and the distance between adjacent layers is about 3.6Å which is in quite agreement with recent epitaxial growth of stanene [7]. The positive phonon spectrum confirms the thermal stability of material. In Figure the phonon frequencies of bilayer stanene at the Γ point are plotted for all acoustic and optical branches. The maximum phonon frequency at the Γ point is 76 cm which is almost unchanged with respect to monolayer stanene [38]. For the monolayer stanene there is a phonon frequency gap between 5 cm and cm [38]. In the case of bilayer stanene the frequency gap is filled with one phonon state at 6 cm for theγpoint. Figure 3 shows the band structure and partial density of state (PDOS) for pure bilayer stanene and FBSn in the absence of spin-orbit coupling (SOC). The band structure is plotted along the high symmetry points of the first Brillouin zone. In the pure bilayer structure, the bands at the K point cross the Fermi level which indicate the semimetal behaviour. In both valence and conduction bands the main contribution comes from atomic p orbitals of tin atoms. In the FBSn structure, unlike graphene the energy valence and conduction bands at the K point are separated by a relatively large energy gap. The minimum of conduction band at theγpoint cross the Fermi level and electronic gap is zero. The PDOS indicates that the eigenstates related to the p atomic orbital of tin atoms fill the band gap at theγpoint and FBSn becomes a semimetal. Around the Fermi level, electronic states in the valence band mostly arise from the p atomic orbitals but for conduction band the s atomic orbitals contribute more effectively. The close energies as well as difference in parities for s and p atomic orbitals allows the possibility of band inversion and the topological phase in tin D structures. The presence of s atomic orbitals in the conduction band shows the failure of any single orbital model for describing FBSn. The s and p atomic orbitals of fluorine atoms are far from the Fermi level and do not contribute in the electronic structure of bilayer stanene effectively. In both structures, the contribution of d atomic orbitals of tin atoms is negligible and ignored in figures. As follows, we consider the spin-orbit coupling and the van der Waals interaction between different layers of atoms in the DFT-D model. The atomic positions of structure are fully relaxed for a certain value of strain and the variation of total energy with lattice parameter is plotted in Figure 4 for two different models.the polynomial curve is well fitted to the DFT and DFT-D results are shown in Figure 4. For a D structure one cannot define the Young s modulus due to the very large value of c axis in the hexagonal super cell. We calculate the in-plane stiffness Y s = (/A )( E/ ε ) where A is the equilibrium area of the unit cell. We obtain the Y s = 6N/m which is much smaller than graphene s 3

4 (a) Bilayer stanene (a') Total Sn(p) Sn(s) (b) Fluorinated bilayer stanene (b') Total Sn(p) Sn(s) F (p) M K M 3 Density of States [ev ] Figure 3: Band structure and PDOS of (a, a )bilayer stanene and (b, b ) FBSn in the absence of spin-orbit interaction. in-plane stiffness [7, 39, 4] and is about twice of pure single layer stanene [8] The binding energy between two fluorinated layers which is expressed as in equation, E Binding = E Bilayer f luorinated E Monolayer f luorinated () The binding energy between two layers is about.85 ev/unitcell which is increased % and reaches to.5 ev/unitcell by considering the van der Waals interaction between two tin layers. The electronic structure of relaxed FBSn in the DFT-D model with spin-orbit interaction is shown in Figure 5. In the presence of spin-orbit interaction some degeneracies in the band structure are lifted. By considering the intrinsic spin-orbit interaction, the required gap for the topological insulators does not appear in the electronic structure and FBSn still remains semimetal. For the relaxed structure the electronic band structure of system is gapless and therefore the topological invariant is ill-defined due to the mixture of valence and conduction bands at the Γ point. The small value of in-plane stiffness (Y s ) and flexibility due to the buckling[] suggest the tuning of electronic structure by applying external strain. For a structure with lattice parameter a, the applied strain is defined asε=(a a )/a. We apply tensile strain to the hexagonal structure in order to open an energy band gap at theγpoint. Figure 6 shows the band structure of bilayer FBSn in the DFT-D with spin-orbit interaction under strain up to.8. By applying external strain the atomic bonds grow weaker which leads to reduction in electronic band widths. The applied strain reduces the band curvature and electronic velocities around the Fermi level. The electronic gaps at the M and K points are decreased but the main effect of tensile strain is for electronic bands around the Γ point. The overlap between valence band maximum (VBM) and conduction band minimum (CBM) at the Γ point is decreased and finally a very small electronic indirect electronic gap is opened at the Γ point. To investigate the topological phase of FBSn we calculate the Z topological invariant by crossing the parities of occupied bands [4, 5, 35, 36]. Table indicates the occupied bands paritiy at four TRIM points for FBSn unders percent of biaxial tensile strain. For higher values of applied biaxial strain thez topological invariant is same. The parity analysis reveals that even for relatively high values of strain, the band inversion does not occur in the FBSn structure and topological invariant is zero which indicates the trivial phase in strained structures. Although the monolayer fluorinated stanene has a spin-orbit gap of.3ev [4] and is an appropriate candidate for topological insulator, the bilayer stanene is semimetal and gapless. The 4

5 .8.6 DFT+SOC DFT-D+SOC polynomial fit strain( a/a ) Figure 4: Energy difference from the ground state for bilayer FBSn in the presence of spin-orbit and vdw interactions applied tensile strain opens an indirect band gap around the point but the strained structure is a trivial insulator which is a drawback for application of fluorinated stanene as a topological insulator in electronic devices. DFT-D+SOC Figure 5: Electronic band structure of FBSn with spin-orbit and van der Waals interactions 4. Conclusion In summary, we have studied the fluorinated bilayer stanene in the DFT model and consider the intrinsic spin-orbit and van der Waals interaction between different layers. Although the fluorinated monolayer stanene was predicted as 5

6 .5 Precent 6 Percent 9 Percent Percent 5 Percent 8 Percent M K M K M K Figure 6: The electronic band structure of FBSn in the DFT-D approximation with tensile strain up to a/a =.8. Table : Z Topological invariant Γ i Parity of occupied bands σ i Γ (.,.) M (.5,.) M (.,.5) M 3 (.5,.5) Strained FBiSn Z Topological invariant Z = a room temperature topological insulator[4], the bilayer stanene and fluorinated bilayer stanene are semi-metal. The the high values of applied tensile strain causes an electronic band gap and trivial insulator phase in the fluorinated bilayer stanene. 5. Acknowledgment One of the authors (Y. Mogulkoc) acknowledges The Scientific and Technological Research Council of Turkey (TUBITAK) through BIDEB9 Postdoctoral Research Fellowship. The numerical calculations reported in this paper were partially performed at TUBITAK ULAKBIM, High Performance and Grid Computing Center (TRUBA resources). References References [] C. L. Kane and E. J. Mele, Physical Review Letters 95 (5) 468. [] F. Zhang, and C. Zhang, Nanoscale Research Letters 4 () 4. [3] C. Zhang, and S Yan, The Journal of Physical Chemistry C 6 () 463. [4] R. Zhang, C. Zhang, W. Ji, S. Li, S. Yan, S. Hu, P. Li, P. Wang, and F. Li, Scientific Reports 6 (6)

7 [5] H. Zhao, C. Zhang, R. Zhang, S. Li, S. Yan, B. Zhang, P. Li, and P. Wang, Scientific Reports 6 (6) 5. [6] MZ. Hasan, CL. Kane, Review Modern Physics 8 () 345. [7] Q. L. Qi, S. C. Zhang, Review Modern Physics 83 () 57. [8] B. Yan,S. C. Zhang, Reports on Progress in Physics 75 () 965. [9] J. E. Moore, Nature Nanotechnology 8 (3) 94. [] B. Rasche, A. Isaeva, M. Ruck et al, Nature Material (3) 4. [] S. Nigam, S. Gupta, D. Banyai, R. Pandey, and C. Majumder, Physical Chemistry Chemical Physics 7 (5) 675. [] M. Knig, S. Wiedmann, C. Brne, A. Roth, H. Buhmann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 38 (7) 766. [3] D. Hsieh, Y. Xia, D. Qian, L. Wray, J. Dil, F. Meier, J. Osterwalder, L. Patthey, J. Checkelsky, and N. Ong, Nature 46 (9). [4] Y. Tanaka, Z. Ren, T. Sato, K. Nakayama, S. Souma, T. Takahashi, K. Segawa, and Y. Ando, Nature Physics 8 () 8. [5] Y. Chen, J. Analytis, J.-H. Chu, Z. Liu, S.-K. Mo, X.-L. Qi, H. Zhang, D. Lu, X. Dai, and Z. Fang, Science 35 (9) 78. [6] K. Kuroda, M. Ye, A. Kimura, S. Eremeev, E. Krasovskii, E. Chulkov, Y. Ueda, K. Miyamoto, T. Okuda, and K. Shimada, Physical Review Letters 5 () 468. [7] B.-H. Chou, Z.-Q. Huang, C.-H. Hsu, F.-C. Chuang, Y.-T. Liu, H. Lin, and A. Bansil, New Journal of Physics 6 (4) 58. [8] G. Cao, Y. Zhang, and J. Cao, Physics Letters A 379 (5) 475. [9] R.-W. Zhang, C.-W. Zhang, W.-X. Ji, S.-S. Li, S.-J. Hu, S.-S. Yan, P. Li, P.-J. Wang, and F. Li, New Journal of Physics 7 (5) [] Y. Ren, Z. Qiao, and Q. Niu, arxiv preprint arxiv:59.96 (5). [] C.-C. Liu, W. Feng, and Y. Yao, Physical Review Letters 7 () 768. [] H. Zhang, Y. Ma, and Z. Chen, Nanoscale (5). [3] X.-F. Wang, Y. Hu, and H. Guo, Physical Review B 85 () 44. [4] Y. Xu, B. Yan, H.-J. Zhang, J. Wang, G. Xu, P. Tang, W. Duan, and S.-C. Zhang, Physical Review Letters (3) [5] M. Modarresi, W. B. Kuang, T. P. Kaloni, M. R. Roknabadi, and G. Schreckenbach, AIP Advances 6 (6) 959. [6] C.-C. Liu, H. Jiang, Y. Yao, Physical Review B 84 () [7] F.-f. Zhu, W.-j. Chen, Y. Xu, C.-l. Gao, D.-d. Guan, C.-h. Liu, D. Qian, S.-C. Zhang, and J.-f. Jia, Nature Materials 4 (5). [8] M. Modarresi, A. Kakoee, Y. Mogulkoc, and M. Roknabadi, Computational Materials Science (5) 64. [9] P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, and I. Dabo, Journal of Physics: Condensed Matter (9) [3] J. P. Perdew, K. Burke, and M. Ernzerhof, Physical review letters 77 (996) S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, The Journal of Chemical Physics 3 () 544. [3] S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, The Journal of Chemical Physics 3 () 544. [3] S. Grimme, Journal of Computational Chemistry 7 (6) 787. [33] H. J. Monkhorst and J. D. Pack, Physical Review B 3 (976) 588. [34] L. Fu and C. L. Kane, Physical Review B 76 (7) 453. [35] R. Zhang, C. Zhang, W. Ji, P. Li, P. Wang, S. Li, and S. Yan, Applied Physics Letters 9 (6) 89. [36] Y. Wang, W. Ji, C. Zhang, P. Li, F. Li, P. Wang, S. Li, and S. Yan, Applied Physics Letters 8 (6) 734, [37] A. Kokalj, Journal of Molecular Graphics and Modelling 7 (999) 76. [38] S. J. Zaveh and M. R. Roknabadi and T. Morshedloo and M. Modarresi, Superlattices and Microstructures 9 (6) 383. [39] M. Mirnezhad, M. Modarresi, R. Ansari, and M. Roknabadi, Journal of Thermal Stresses 35 () 93. [4] Q. Peng, C. Liang, W. Ji, and S. De, Physical Chemistry Chemical Physics 5 (3) 3. 7