Brillouin Zone Viewer Tool

Total Page:16

File Type:pdf, Size:1020Kb

Brillouin Zone Viewer Tool Network for Computational Nanotechnology (NCN) Purdue, Norfolk State, Northwestern, MIT, Molecular Foundry, UC Berkeley, Univ. of Illinois, UTEP Brillouin Zone Viewer Tool Arun Goud Sebastian Steiger Michael Povolotskyi Tillmann Kubis Hong-Hyun Park Gerhard Klimeck Electrical and Computer Engineering Purdue University, West Lafayette IN, USA Brillouin Zone Tool Purpose of the tool Users can visualize 1st Brillouin zones of the following common crystal lattice systems: - Cubic systems (Simple, BCC & FCC) - Hexagonal (Wurtzite) - Rhombohedral (Bi2TeO3) - Graphene (Honeycomb) Tab1 -> Choose lattice type Tab2 -> Define lattice vector coordinates & number of translations along each coordinate needed to create the unit cell Theory Crystal lattice Bravais lattice atomic basis Real-space lattice: As Ga 2 atoms zincblende fcc A,B,C = Bravais lattice vectors Reciprocal lattice: Atomic basis has no influence on rec. lattice 1st Brillouin zone = set of points in reciprocal space which are closer to 0 than to any other reciprocal lattice vector Features Both primitive & conventional (cubic) lattice translation vectors can be specified for the cubic systems Only primitive translation vectors are currently supported for - Graphene - Hexagonal - Rhombohedral Output of the tool: • Brillouin zone • VTK file • Real space structure & unit cell • Input deck for NEMO5 • Output log Inputs - Lattice Vectors & Translations FCC example: Cubic indices Primitive indices . Translation determines multiples of the A, B & C vectors . Shape of unit cell & hence the Brillouin zone will be altered by changing the translation field entry . Input lattice parameters (depends on lattice type) History of the Tool The Brillouin Zone viewer tool is currently in its 1.0 version. Rappture GUI for the tool by Arun Goud. NEMO5 simulation engine by Sebastian Steiger, Michael Povolotskyi, Tillmann Kubis and Hong-Hyun Park. General supervision: Dr. Gerhard Klimeck. Last update of this document: April 2011.
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]