Strain-Induced Large Band-Gap Topological Insulator in a New Stable Silicon Allotrope: Dumbbell Silicene

Total Page:16

File Type:pdf, Size:1020Kb

Strain-Induced Large Band-Gap Topological Insulator in a New Stable Silicon Allotrope: Dumbbell Silicene Strain-induced large band-gap topological insulator in a new stable silicon allotrope: dumbbell silicene Tian Zhang1, Yan Cheng1,∗ Xiang-Rong Chen1,† and Ling-Cang Cai2 1Institute of Atomic and Molecular Physics, College of Physical Science and Technology, Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Sichuan University, Chengdu 610065, China 2National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621900, China (Dated: Oct. 30, 2015) By the generalized gradient approximation in framewok of density functional theory, we investigate a 2D topological insulator of new silicon allotrope (call dumbbell silicene synthesized recently by Cahangirov et al) through tuning external compression strain, and find a topological quantum phase transition from normal to topological insulator, i.e., the dumbbell silicene can turn a two-dimensional topological insulator with an inverted band gap. The obtained maximum topological nontrivial band gap about 12 meV under isotropic strain is much larger than that for previous silicene, and can be further improved to 36 meV by tuning anisotropic strain, which is sufficiently large to realize quantum spin Hall effect even at room-temperature, and thus is beneficial to the fabrication of high- speed spintronics devices. Furthermore, we confirm that the boron nitride sheet is an ideal substrate for the experimental realization of the dumbbell silicene under external strain, maintaining its nontrivial topology. These properties make the two-dimensional dumbbell silicene a good platform to study novel quantum states of matter, showing great potential for future applications in modern silicon-based microelectronics industry. PACS numbers: 73.22.-f, 73.43.-f, 71.70.Ej, 85.75.-d Two-dimensional (2D) materials have been a focus of mobility24,25. Unfortunately, its tiny band gap (about intense research in recent years1–7. As opposed to three- 8×10−4 meV26) opened by spin-orbit coupling (SOC) ef- dimensional (3D) one, the optical, electronic, mechanical fect seriously limits its device applications. Subsequently, and thermal properties of the 2D materials are easily ad- silicene with a low buckled honeycomb lattice was syn- justed by external strains, defects, electric field, or stack- thesized and predicted quickly to be a new 2D TI with a ing orders8–11, and thus its realistic performance can be relatively large spin-orbit gap of 1.55 meV27,28. Almost readily improved through current microfabrication tech- every striking property of graphene can be transferred to nology. More fascinatingly, 2D modes12,13 was first pre- this innovative material. Indeed, these features together dicted with quantum spin Hall (QSH) effect, and recently with the natural compatibility with current silicon-based more and more 2D materials have been confirmed as 2D microelectronics industry make silicene a promising can- topological insulators (TIs)14–18, also known as QSH in- didate for future nanoelectronics application. However, sulators. 2D TIs are novel materials characterized by its topological nontrivial band gap is still small and limits a bulk energy gap and gapless spin-filtered edge states. its room-temperature application in spintronics. Hence, Different from surface states of 3D TIs, which is only free there is great interest in searching for room temperature from exact 1800-backscattering and suffers from scatter- 2D TI (topological band gap about 26 meV), especially ing of other angles, the special edges of 2D TIs are topo- for silicon-based 2D TI. logically protected by the time reversal symmetry and In this work, we will investigate a 2D topological insu- can immune to nonmagnetic scattering and geometry lator of new silicon allotrope (call dumbbell silicene syn- perturbations, thus open new ways for backscattering- thesized recently by Cahangirov et al29) through tuning free transport. Clearly, 2D TIs is better than 3D TIs external strain, then find a topological quantum phase for coherent nondissipative spin transport related appli- transition from normal to topological insulator, accompa- arXiv:1510.08338v1 [cond-mat.mtrl-sci] 28 Oct 2015 cations. nied by a band inversion at Γ point that changes Z2 num- The 2D TIs have stimulated enormous research ac- ber from 0 to 1. Isotropic strain can induce the quantum tivities in condensed matter physics due to their novel topological phase transition in the dumbbell silicene and QSH effect and the potential application in quan- modulate its topological nontrivial band gap regularly. tum computation and spintronics19,20. Well-known The maximum topological nontrivial band gap about 12 2D TI graphene, a monolayer of carbon atoms form- meV under isotropic strain can be further improved to ing a similar honeycomb lattice, hosts a miraculous 35 meV by tuning anisotropic strain, which is sufficiently electronic system, and thus becomes perfect breeding large to realize QSH effect of the dumbbell silicene even ground for a variety of exotic quantum phenomena, such at room temperature. The pristine dumbbell silicene also as quantum anomalous Hall effect (QAHE), Majorana owns a Dirac cone at K point originated from pz orbital, fermions and superconductor21–23. Furthermore, mass- like in previous silicene, but the Dirac zone can be easily less Dirac fermions endow graphene with superior carrier destroyed as a result of spatial symmetry breaking under 2 (a)a1 (b) dumbbell silicene silicene (a) (b) 16 16 a1 14 14 12 12 2 10 10 a2 a 8 8 6 6 s = 4.29 Å s = 2.58Å Frequency(THz) Frequency(THz) 4 4 Si(α) h= 0.44Å 2 2 h= 2.71Å Si(β) Si 0 0 dumbbell silicene silicene Γ M K Γ Γ M K Γ (c) 3 (c) (d) 4.46 5.94 7.43 8.92 10.40 (d) 3 b1 -2.5 2 2 dumbbell silicene M Γ -3.0 silicene 1 1 K -3.5 0 0 -4.0 Energy(eV) HB -1 Energy(eV) -1 E nergy (eV/atom) E nergy -4.5 L B -5.0 b2 -2 -2 2.30 3.07 3.84 4.61 5.38 -3 -3 L attice constant (angstrom) M Γ K M M Γ K M FIG. 1: (Color online) (a) Top and side view of the dumbbell FIG. 2: (Color online) Vibrational band structure ω(~k) of (a) silicene structure. (b) Top and side view of the silicene struc- the dumbbell silicene and (b) the silicene. The slope of the red ture. Pink dashed lines mark the hexagonal lattice; ~a1 and dashed lines along the longitudinal acoustic branches near Γ ~a2 are the lattice unit vectors. (c) The first Brillouin zone. corresponds to the speed of sound and the in-plane stiffness. ~b1 and ~b2 are the reciprocal lattice unit vectors. Γ, M and K Electronic band structure of (c) the dumbbell silicene and refer to special points in the first Brillouin zone. (d) Total (d) the silicene without SOC. Dirac cone is marked by black energies of the dumbbell silicene and silicene with different rectangular box. The green dots corresponding to Si-pz. The lattice constants. Fermi level is set to zero. anisotropic strain. It implies that the hexagonal symme- ditional, we calculated the phonon vibration spectra by try is the main precondition for the present of Dirac cone using the frozen phonon method as implemented in the in the dumbbell silicene. In additional, we further show PHONOPY code37. that the boron nitride (BN) substrate is suitable to sup- Crystal structure of the dumbbell silicene is presented port the dumbbell silicen under external strain, and the in Fig. 1(a). Different from silicene [see Fig. 1(b)], in the topologically nontrivial properties of the strained dumb- dumbbell silicene there are two types of Si atoms: one is bell silicene can be retained well in the realistic growth fourfold coordinated [denoted by Si(α)], and the other on the boron nitride (BN) substrate. one is only threefold coordinated with a dangling bond The first-principle calculations based on the density [denoted by Si(β)]. It is easy to find that all the Si(α) functional theory were performed using the Vienna ab atoms keep on the same plane, like in graphene, which simulation package (VASP)30,31. We use the general- is surrounded by the Si(β) atoms to form the dumbbell ized gradient approximations (GGA) of Perdew-Burke- structure. The optimized dumbbell silicene has a hexag- 32 Ernzerhof (PBE) for electron-electron interactions onal structure with space group D6h (P 6/mmm) and ten and the projector-augmented-wave (PAW)33 pseudo- Si atoms in one unit cell, which is different from the op- potentials in the plane-wave basis with an energy cutoff timized silicene with space group D3d (P 3m1) and only of 700 eV. The Brillouin zone (BZ) was sampled using two Si atoms in one unit cell. The buckling height of the an 13×13×1 gamma-centered Monkhorst-Pack grid34. A dumbbell silicene [h in Fig. 1(a) and (b)] about 2.71 A˚ is slab model, together with a vacuum layer larger than 30 much larger than that in the silicene about 0.44 A,˚ which A, was employed to avoid spurious interactions due to may result in the larger SOC strength for the dumbbell the nonlocal nature of the correlation energy35. The lat- silicene. The obtained lattice constant and the distance tice constants are relaxed until the residual force on each between the neighboring dumbbell units (s in Fig. 1(a)) atom is less than 0.01 eV/A˚ and the total energy is less for the dumbbell silicene are 7.43 A˚ and 4.29 A,˚ respec- than 1×10−6 eV. Spin-orbit coupling (SOC)36 is included tively. These results including those for the silicene are in in the calculations after the structural relaxations.
Recommended publications
  • Lecture Notes
    Solid State Physics PHYS 40352 by Mike Godfrey Spring 2012 Last changed on May 22, 2017 ii Contents Preface v 1 Crystal structure 1 1.1 Lattice and basis . .1 1.1.1 Unit cells . .2 1.1.2 Crystal symmetry . .3 1.1.3 Two-dimensional lattices . .4 1.1.4 Three-dimensional lattices . .7 1.1.5 Some cubic crystal structures ................................ 10 1.2 X-ray crystallography . 11 1.2.1 Diffraction by a crystal . 11 1.2.2 The reciprocal lattice . 12 1.2.3 Reciprocal lattice vectors and lattice planes . 13 1.2.4 The Bragg construction . 14 1.2.5 Structure factor . 15 1.2.6 Further geometry of diffraction . 17 2 Electrons in crystals 19 2.1 Summary of free-electron theory, etc. 19 2.2 Electrons in a periodic potential . 19 2.2.1 Bloch’s theorem . 19 2.2.2 Brillouin zones . 21 2.2.3 Schrodinger’s¨ equation in k-space . 22 2.2.4 Weak periodic potential: Nearly-free electrons . 23 2.2.5 Metals and insulators . 25 2.2.6 Band overlap in a nearly-free-electron divalent metal . 26 2.2.7 Tight-binding method . 29 2.3 Semiclassical dynamics of Bloch electrons . 32 2.3.1 Electron velocities . 33 2.3.2 Motion in an applied field . 33 2.3.3 Effective mass of an electron . 34 2.4 Free-electron bands and crystal structure . 35 2.4.1 Construction of the reciprocal lattice for FCC . 35 2.4.2 Group IV elements: Jones theory . 36 2.4.3 Binding energy of metals .
    [Show full text]
  • Introduction to the Physical Properties of Graphene
    Introduction to the Physical Properties of Graphene Jean-No¨el FUCHS Mark Oliver GOERBIG Lecture Notes 2008 ii Contents 1 Introduction to Carbon Materials 1 1.1 TheCarbonAtomanditsHybridisations . 3 1.1.1 sp1 hybridisation ..................... 4 1.1.2 sp2 hybridisation – graphitic allotopes . 6 1.1.3 sp3 hybridisation – diamonds . 9 1.2 Crystal StructureofGrapheneand Graphite . 10 1.2.1 Graphene’s honeycomb lattice . 10 1.2.2 Graphene stacking – the different forms of graphite . 13 1.3 FabricationofGraphene . 16 1.3.1 Exfoliatedgraphene. 16 1.3.2 Epitaxialgraphene . 18 2 Electronic Band Structure of Graphene 21 2.1 Tight-Binding Model for Electrons on the Honeycomb Lattice 22 2.1.1 Bloch’stheorem. 23 2.1.2 Lattice with several atoms per unit cell . 24 2.1.3 Solution for graphene with nearest-neighbour and next- nearest-neighour hopping . 27 2.2 ContinuumLimit ......................... 33 2.3 Experimental Characterisation of the Electronic Band Structure 41 3 The Dirac Equation for Relativistic Fermions 45 3.1 RelativisticWaveEquations . 46 3.1.1 Relativistic Schr¨odinger/Klein-Gordon equation . ... 47 3.1.2 Diracequation ...................... 49 3.2 The2DDiracEquation. 53 3.2.1 Eigenstates of the 2D Dirac Hamiltonian . 54 3.2.2 Symmetries and Lorentz transformations . 55 iii iv 3.3 Physical Consequences of the Dirac Equation . 61 3.3.1 Minimal length for the localisation of a relativistic par- ticle ............................ 61 3.3.2 Velocity operator and “Zitterbewegung” . 61 3.3.3 Klein tunneling and the absence of backscattering . 61 Chapter 1 Introduction to Carbon Materials The experimental and theoretical study of graphene, two-dimensional (2D) graphite, is an extremely rapidly growing field of today’s condensed matter research.
    [Show full text]
  • Phys 446: Solid State Physics / Optical Properties
    Phys 446: Solid State Physics / Optical Properties Fall 2015 Lecture 2 Andrei Sirenko, NJIT 1 Solid State Physics Lecture 2 (Ch. 2.1-2.3, 2.6-2.7) Last week: • Crystals, • Crystal Lattice, • Reciprocal Lattice Today: • Types of bonds in crystals Diffraction from crystals • Importance of the reciprocal lattice concept Lecture 2 Andrei Sirenko, NJIT 2 1 (3) The Hexagonal Closed-packed (HCP) structure Be, Sc, Te, Co, Zn, Y, Zr, Tc, Ru, Gd,Tb, Py, Ho, Er, Tm, Lu, Hf, Re, Os, Tl • The HCP structure is made up of stacking spheres in a ABABAB… configuration • The HCP structure has the primitive cell of the hexagonal lattice, with a basis of two identical atoms • Atom positions: 000, 2/3 1/3 ½ (remember, the unit axes are not all perpendicular) • The number of nearest-neighbours is 12 • The ideal ratio of c/a for Rotated this packing is (8/3)1/2 = 1.633 three times . Lecture 2 Andrei Sirenko, NJITConventional HCP unit 3 cell Crystal Lattice http://www.matter.org.uk/diffraction/geometry/reciprocal_lattice_exercises.htm Lecture 2 Andrei Sirenko, NJIT 4 2 Reciprocal Lattice Lecture 2 Andrei Sirenko, NJIT 5 Some examples of reciprocal lattices 1. Reciprocal lattice to simple cubic lattice 3 a1 = ax, a2 = ay, a3 = az V = a1·(a2a3) = a b1 = (2/a)x, b2 = (2/a)y, b3 = (2/a)z reciprocal lattice is also cubic with lattice constant 2/a 2. Reciprocal lattice to bcc lattice 1 1 a a x y z a2 ax y z 1 2 2 1 1 a ax y z V a a a a3 3 2 1 2 3 2 2 2 2 b y z b x z b x y 1 a 2 a 3 a Lecture 2 Andrei Sirenko, NJIT 6 3 2 got b y z 1 a 2 b x z 2 a 2 b x y 3 a but these are primitive vectors of fcc lattice So, the reciprocal lattice to bcc is fcc.
    [Show full text]
  • Part 1: Supplementary Material Lectures 1-5
    Crystal Structure and Dynamics Paolo G. Radaelli, Michaelmas Term 2014 Part 1: Supplementary Material Lectures 1-5 Web Site: http://www2.physics.ox.ac.uk/students/course-materials/c3-condensed-matter-major-option Contents 1 Frieze patterns and frieze groups 3 2 Symbols for frieze groups 4 2.1 A few new concepts from frieze groups . .5 2.2 Frieze groups in the ITC . .7 3 Wallpaper groups 8 3.1 A few new concepts for Wallpaper Groups . .8 3.2 Lattices and the “translation set” . .9 3.3 Bravais lattices in 2D . .9 3.3.1 Oblique system . 10 3.3.2 Rectangular system . 10 3.3.3 Square system . 10 3.3.4 Hexagonal system . 11 3.4 Primitive, asymmetric and conventional Unit cells in 2D . 11 3.5 The 17 wallpaper groups . 12 3.6 Analyzing wallpaper and other 2D art using wallpaper groups . 12 1 4 Fourier transform of lattice functions 13 5 Point groups in 3D 14 5.1 The new generalized (proper & improper) rotations in 3D . 14 5.2 The 3D point groups with a 2D projection . 15 5.3 The other 3D point groups: the 5 cubic groups . 15 6 The 14 Bravais lattices in 3D 19 7 Notation for 3D point groups 21 7.0.1 Notation for ”projective” 3D point groups . 21 8 Glide planes in 3D 22 9 “Real” crystal structures 23 9.1 Cohesive forces in crystals — atomic radii . 23 9.2 Close-packed structures . 24 9.3 Packing spheres of different radii . 25 9.4 Framework structures . 26 9.5 Layered structures .
    [Show full text]
  • Valleytronics: Opportunities, Challenges, and Paths Forward
    REVIEW Valleytronics www.small-journal.com Valleytronics: Opportunities, Challenges, and Paths Forward Steven A. Vitale,* Daniel Nezich, Joseph O. Varghese, Philip Kim, Nuh Gedik, Pablo Jarillo-Herrero, Di Xiao, and Mordechai Rothschild The workshop gathered the leading A lack of inversion symmetry coupled with the presence of time-reversal researchers in the field to present their symmetry endows 2D transition metal dichalcogenides with individually latest work and to participate in honest and addressable valleys in momentum space at the K and K′ points in the first open discussion about the opportunities Brillouin zone. This valley addressability opens up the possibility of using the and challenges of developing applications of valleytronic technology. Three interactive momentum state of electrons, holes, or excitons as a completely new para- working sessions were held, which tackled digm in information processing. The opportunities and challenges associated difficult topics ranging from potential with manipulation of the valley degree of freedom for practical quantum and applications in information processing and classical information processing applications were analyzed during the 2017 optoelectronic devices to identifying the Workshop on Valleytronic Materials, Architectures, and Devices; this Review most important unresolved physics ques- presents the major findings of the workshop. tions. The primary product of the work- shop is this article that aims to inform the reader on potential benefits of valleytronic 1. Background devices, on the state-of-the-art in valleytronics research, and on the challenges to be overcome. We are hopeful this document The Valleytronics Materials, Architectures, and Devices Work- will also serve to focus future government-sponsored research shop, sponsored by the MIT Lincoln Laboratory Technology programs in fruitful directions.
    [Show full text]
  • Electronic Structure Theory for Periodic Systems: the Concepts
    Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation • There are 1020 atoms in 1 mm3 of a solid. • It is impossible to track all of them. • But for most solids, atoms are arranged periodically Christian Ratsch, UCLA IPAM, DFT Hands-On Workshop, July 22, 2014 Outline Part 1 Part 2 • Periodicity in real space • The supercell approach • Primitive cell and Wigner Seitz cell • Application: Adatom on a surface (to calculate for example a diffusion barrier) • 14 Bravais lattices • Surfaces and Surface Energies • Periodicity in reciprocal space • Surface energies for systems with multiple • Brillouin zone, and irreducible Brillouin zone species, and variations of stoichiometry • Bandstructures • Ab-Initio thermodynamics • Bloch theorem • K-points Christian Ratsch, UCLA IPAM, DFT Hands-On Workshop, July 22, 2014 Primitive Cell and Wigner Seitz Cell a2 a1 a2 a1 primitive cell Non-primitive cell Wigner-Seitz cell = + 푹 푁1풂1 푁2풂2 Christian Ratsch, UCLA IPAM, DFT Hands-On Workshop, July 22, 2014 Primitive Cell and Wigner Seitz Cell Face-center cubic (fcc) in 3D: Primitive cell Wigner-Seitz cell = + + 1 1 2 2 3 3 Christian푹 Ratsch, UCLA푁 풂 푁 풂 푁 풂 IPAM, DFT Hands-On Workshop, July 22, 2014 How Do Atoms arrange in Solid Atoms with no valence electrons: Atoms with s valence electrons: Atoms with s and p valence electrons: (noble elements) (alkalis) (group IV elements) • all electronic shells are filled • s orbitals are spherically • Linear combination between s and p • weak interactions symmetric orbitals form directional orbitals, • arrangement is closest packing • no particular preference for which lead to directional, covalent of hard shells, i.e., fcc or hcp orientation bonds.
    [Show full text]
  • Electronic and Structural Properties of Silicene and Graphene Layered Structures
    Wright State University CORE Scholar Browse all Theses and Dissertations Theses and Dissertations 2012 Electronic and Structural Properties of Silicene and Graphene Layered Structures Patrick B. Benasutti Wright State University Follow this and additional works at: https://corescholar.libraries.wright.edu/etd_all Part of the Physics Commons Repository Citation Benasutti, Patrick B., "Electronic and Structural Properties of Silicene and Graphene Layered Structures" (2012). Browse all Theses and Dissertations. 1339. https://corescholar.libraries.wright.edu/etd_all/1339 This Thesis is brought to you for free and open access by the Theses and Dissertations at CORE Scholar. It has been accepted for inclusion in Browse all Theses and Dissertations by an authorized administrator of CORE Scholar. For more information, please contact [email protected]. Electronic and Structural Properties of Silicene and Graphene Layered Structures A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science by Patrick B. Benasutti B.S.C.E. and B.S., Case Western Reserve and Wheeling Jesuit, 2009 2012 Wright State University Wright State University SCHOOL OF GRADUATE STUDIES August 31, 2012 I HEREBY RECOMMEND THAT THE THESIS PREPARED UNDER MY SUPERVI- SION BY Patrick B. Benasutti ENTITLED Electronic and Structural Properties of Silicene and Graphene Layered Structures BE ACCEPTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master of Science in Physics. Dr. Lok C. Lew Yan Voon Thesis Director Dr. Douglas Petkie Department Chair Committee on Final Examination Dr. Lok C. Lew Yan Voon Dr. Brent Foy Dr. Gregory Kozlowski Dr. Andrew Hsu, Ph.D., MD Dean, Graduate School ABSTRACT Benasutti, Patrick.
    [Show full text]
  • Symmetry in the Solid State - Part IV: Brillouin Zones and the Symmetry of the Band Structure
    Lecture 4 — Symmetry in the solid state - Part IV: Brillouin zones and the symmetry of the band structure. 1 Symmetry in Reciprocal Space — the Wigner-Seitz construc- tion and the Brillouin zones Non-periodic phenomena in the crystal (elastic or inelastic) are described in terms of generic (non-RL) reciprocal-space vectors and give rise to scattering outside the RL nodes. In describing these phenomena, however, one encounters a problem: as one moves away from the RL origin, symmetry-related “portions” of reciprocal space will become very distant from each other. In order to take full advantage of the reciprocal-space symmetry, it is therefore advanta- geous to bring symmetry-related parts of the reciprocal space together in a compact form. This is exactly what the Wigner-Seitz (W-S) constructions accomplish very cleverly. ”Brillouin zones” is the name that is given to the portions of the extended W-S contructions that are “brought back together”. A very good description of the Wigner-Seitz and Brillouin constructions can be found in [3]. 1.1 The Wigner-Seitz construction The W-S construction is essentially a method to construct, for every Bravais lattice, a fully- symmetric unit cell that has the same volume of a primitive cell. As such, it can be applied to both real and reciprocal spaces, but it is essentially employed only for the latter. For a given lattice node τ , the W-S unit cell containing τ is the set of points that are closer to τ than to any other lattice node. It is quite apparent that: • Each W-S unit cell contains one and only one lattice node.
    [Show full text]
  • Reciprocal Space and Brillouin Zones in Two and Three Dimensions As
    Reciprocal Space and Brillouin Zones in Two and Three Dimensions As briefly stated at the end of the first section, Bloch’s theorem has the following form in two and three dimensions: 2ikR k (r + R) = e k(r) . In this expression, R is a lattice vector between a pair of unit cells: R = ua + vb + wc ; u, v, and w are integers and the dot product k R = kau + kbv + kcw . (In two dimensions, R = ua + vb and k R = kau + kbv .) At this point, we need to clarify the meaning of the vector k and find a way to define the two- and three-dimensional Brillouin zones. To do this, let’s consider the general, oblique, two dimensional lattice below. For this lattice, the basis vectors, a and b, are not orthogonal and there is no symmetry relationship between their lengths. An operation that translates by the lattice vector R (which we may call tR ) involves independent and commutative translation operations in the a and b directions. By an extension of the reasoning used for the one-dimensional translation group, we can deduce that if we operate on a basis function that belongs to the irreducible representation labeled by k, we should obtain the following: R vb ua ua vb 2 i k u 2 i k v t (k) = t t (k) = t t (k) = e a e b (k) where R = ua + vb. R 2 ik R The form t (k) = e (k) makes sense if we define k-space (often called reciprocal space) basis vectors, a* and b*, such that k R = (kaa + kbb ) (ua + vb) = kau + kbv for all u, v, ka, kb .
    [Show full text]
  • PHYSICS 231 Homework 4, Question 4, Graphene
    PHYSICS 231 Homework 4, Question 4, Graphene It has recently become possible to produce samples of graphene one atom thick. A single layer of graphene consists carbon atoms in the form of a honeycomb lattice. There are four valence electrons (two 2s and two 2p electrons). Three of those participate in the chemical bonding and so are in bands well below the Fermi energy. We therefore just consider the bands formed by the one remaining electron. We assume a tight-binding model in which the electron hops between neighboring atoms. The lattice structure is as shown in Fig. 1. We denote the spacing between neighboring atoms by a. B A Figure 1: The solid lines indicate the crystal structure of graphene. The basis has two atoms, labeled A and B. The Bravais lattice (consider, e.g, the lattice formed by the A atoms shown by dashed lines) is triangular with a Bravais lattice spacing 2 sin 60◦ a = √3a, where a is the spacing between × × neighboring atoms. a. There are two atoms per unit cell so 1 band will be filled. b. The Bravais lattice is the same as the lattice formed by all the A atoms, say. As shown in Fig. 1 this is a triangular lattice with lattice spacing √3a. c. From Fig. 1, we see that two basis vectors of the Bravais lattice are a1 = √3 a (1, 0), a2 = a (√3/2, 3/2). (1) The Bravais lattices, b1, b2, are defined such that b a = 2πδ . (2) · j ij To determine them, either consider the formulae discussed in class for three dimensions, with the third basis vector a3 =z ˆ, a unit vector in the z direction, or just figure it out.
    [Show full text]
  • The Two-Dimensional Brillouin Zone of Uniaxially Strained Graphene
    The two-dimensional Brillouin zone of uniaxially strained graphene Marcel Mohr,1, ∗ Konstantinos Papagelis,2 Janina Maultzsch,1 and Christian Thomsen1 1Institut f¨ur Festk¨orperphysik, Technische Universit¨at Berlin, Hardenbergstr. 36, 10623 Berlin, Germany 2 Materials Science Department, 26504 Rio, Patras (Dated: October 22, 2018) We present an in-depth analysis of the electronic and vibrational band structure of uniaxially strained graphene by ab-initio calculations. Depending on the direction and amount of strain, the Fermi crossing moves away from the K-point. However, graphene remains semimetallic under small strains. The deformation of the Dirac cone near the K-point gives rise to a broadening of the 2D Raman mode. In spite of specific changes in the electronic and vibrational band structure the strain-induced frequency shifts of the Raman active E2g and 2D modes are independent of the direction of strain. Thus, the amount of strain can be directly determined from a single Raman measurement. PACS numbers: 71.15.-m, 63.22.-m, 31.15.E- The discovery of graphene in 2004 has led to strong re- specific changes in the electronic and phononic band structure search activities in the last years[1]. Graphenehas been shown for different strain directions are found, the frequencies of the to possess unique material properties. In graphene the quan- Raman active E2g and 2D modes are independent of the di- tum Hall effect could be observed at room temperature [2, 3]. rection of strain. Thus the amount of strain can be directly Due to the specific band structure at the Fermi level, the elec- determined from a single Raman measurement.
    [Show full text]
  • Graphene Bilayers with a Twist
    FOCUS | REVIEW ARTICLE FOCUS | https://doi.org/10.1038/s41563-020-00840-0REVIEW ARTICLE Graphene bilayers with a twist Eva Y. Andrei! !1 and Allan H. MacDonald2 Near a magic twist angle, bilayer graphene transforms from a weakly correlated Fermi liquid to a strongly correlated two-dimensional electron system with properties that are extraordinarily sensitive to carrier density and to controllable envi- ronmental factors such as the proximity of nearby gates and twist-angle variation. Among other phenomena, magic-angle twisted bilayer graphene hosts superconductivity, interaction-induced insulating states, magnetism, electronic nematic- ity, linear-in-temperature low-temperature resistivity and quantized anomalous Hall states. We highlight some key research results in this field, point to important questions that remain open and comment on the place of magic-angle twisted bilayer graphene in the strongly correlated quantum matter world. he two-dimensional (2D) materials field can be traced back to zone. Since σ·pK is the chirality operator, and σy pK σ pK À Á 0 ¼TðÞÁ the successful isolation of single-layer graphene sheets exfoli- is its time-reversal partner ( is the time-reversalI operator), states T ated from bulk graphite1. Galvanized by this breakthrough, in the K and Kʹ Dirac cones haveI opposite chirality. K and Kʹ define T 2 researchers have by now isolated and studied dozens of 2D crystals a valley degree of freedom, which together with spin produces a and predicted the existence of many more3. Because all the atoms in four-fold band degeneracy. The Dirac points separate the conduc- any 2D crystal are exposed to our 3D world, this advance has made tion and valence bands so that at charge neutrality the Fermi sur- it possible to tune the electronic properties of a material without face (pink plane in Fig.
    [Show full text]