* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download Quantum transport and spin effects in lateral semiconductor nanostructures and graphene Martin Evaldsson
Superconductivity wikipedia , lookup
Spin (physics) wikipedia , lookup
Probability amplitude wikipedia , lookup
Bell's theorem wikipedia , lookup
Electromagnetism wikipedia , lookup
Time in physics wikipedia , lookup
Renormalization wikipedia , lookup
EPR paradox wikipedia , lookup
Quantum potential wikipedia , lookup
Aharonov–Bohm effect wikipedia , lookup
Introduction to gauge theory wikipedia , lookup
Dirac equation wikipedia , lookup
Quantum vacuum thruster wikipedia , lookup
Quantum tunnelling wikipedia , lookup
Old quantum theory wikipedia , lookup
Quantum electrodynamics wikipedia , lookup
History of quantum field theory wikipedia , lookup
Photon polarization wikipedia , lookup
Hydrogen atom wikipedia , lookup
Theoretical and experimental justification for the Schrödinger equation wikipedia , lookup
Relativistic quantum mechanics wikipedia , lookup
Introduction to quantum mechanics wikipedia , lookup
Linköping studies in science and technology. Dissertations, No 1202 Quantum transport and spin effects in lateral semiconductor nanostructures and graphene Martin Evaldsson Department of Science and Technology Linköping University, SE-601 74 Norrköping, Sweden Norrköping, 2008 Quantum transport and spin effects in lateral semiconductor nanostructures and graphene c 2008 Martin Evaldsson Department of Science and Technology Campus Norrköping, Linköping University SE-601 74 Norrköping, Sweden ISBN 978-91-7393-835-8 ISSN 0345-7524 Printed in Sweden by LiU-Tryck, Linköping, 2008 Abstract This thesis studies electron spin phenomena in lateral semi-conductor quantum dots/anti-dots and electron conductance in graphene nanoribbons by numerical modelling. In paper I we have investigated spin-dependent transport through open quantum dots, i.e., dots strongly coupled to their leads, within the Hubbard model. Results in this model were found consistent with experimental data and suggest that spin-degeneracy is lifted inside the dot – even at zero magnetic field. Similar systems were also studied with electron-electron effects incorporated via Density Functional Theory (DFT) in the Local Spin Density Approximation (LSDA) in paper II and III. In paper II we found a significant spin-polarisation in the dot at low electron densities. As the electron density increases the spin polarisation in the dot gradually diminishes. These findings are consistent with available experimental observations. Notably, the polarisation is qualitatively different from the one found in the Hubbard model. Paper III investigates spin polarisation in a quantum wire with a realistic external potential due to split gates and a random distribution of charged donors. At low electron densities we recover spin polarisation and a metalinsulator transition when electrons are localised to electron lakes due to ragged potential profile from the donors. In paper IV we propose a spin-filter device based on resonant backscattering of edge states against a quantum anti-dot embedded in a quantum wire. A magnetic field is applied and the spin up/spin down states are separated through Zeeman splitting. Their respective resonant states may be tuned so that the device can be used to filter either spin in a controlled way. Paper V analyses the details of low energy electron transport through a magnetic barrier in a quantum wire. At sufficiently large magnetisation of the barrier the conductance is pinched off completely. Furthermore, if the barrier is sharp we find a resonant reflection close to the pinch off point. This feature is due to interference between a propagating edge state and quasibond state inside the magnetic barrier. Paper VI adapts an efficient numerical method for computing the surface Green’s function in photonic crystals to graphene nanoribbons (GNR). The method is used to investigate magnetic barriers in GNR. In contrast to quantum wires, magnetic barriers in GNRs cannot pinch-off the lowest propagating state. The method is further applied to study edge dislocation defects for realistically sized GNRs in paper VII. In this study we conclude that even modest edge dislocations are sufficient to explain both the energy gap in narrow GNRs, and the lack of dependance on the edge structure for electronic properties in the GNRs. iii iv Preface This thesis summarises my years as a Ph.D. student at the Department of Science and Technology (ITN). It consists of two parts, the first serving as a short introduction both to mesoscopic transport in general and to the papers included in the second part. There are several persons to whom I would like to express my gratitude for help and support during these years: First of all my supervisor Igor Zozoulenko for patiently introducing me to the field of mesoscopic physics and research in general, I have learnt a lot from you. The other members in the Mesoscopic Physics and Photonics group, Aliaksandr Rachachou and Siarhei Ihnatsenka – it is great to have people around who actually understand what I’m doing. Torbjörn Blomquist for providing an excellent C++ matrix library which has significantly facilitated my work. A lot of people at ITN for various reasons, including Mika Gustafsson (a lot of things), Michael Hörnquist and Olof Svensson (lunch company and company), Margarita Gonzáles (company and lussebak), Frédéric Cortat (company and for motivating me to run 5km/year), Steffen Uhlig (company and lussebak), Sixten Nilsson (help with teaching), Sophie Lindesvik and Lise-Lotte Lönnedahl Ragnar (help with administrative issues). Also a general thanks to all the past and present members of the ‘Fantastic Five’, and everyone who has made my coffee breaks more interesting. I would also like to thank my parents and sister for being there, and of course my family, Chamilly, Minna and Morris, for keeping my focus where it matters. Finally, financial support from The Swedish Research Council (VR) and the National Graduate School for Scientific Computations (NGSSC) is acknowledged. v vi List of publications Paper I: M. Evaldsson, I. V. Zozoulenko, M. Ciorga, P. Zawadzki and A. S. Sachrajda. Spin splitting in open quantum dots. Europhysics Letters 68, 261 (2004) Author’s contribution: All calculations for the theoretical part. Plotting of all figures. Initial draft of the of the paper (except experimental part). Paper II: M. Evaldsson and I. V. Zozoulenko. Spin polarization in open quantum dots. Physical Review B 73, 035319 (2006) Author’s contribution: Computations for all figures. Plotting of all figures. Drafting of paper and contributed to the discussion in the process of writing the paper. Paper III: M. Evaldsson, S. Ihnatsenka and I. V. Zozoulenko. Spin polarization in modulation-doped GaAs quantum wires. Physical Review B 77, 165306 (2008) Author’s contribution: Computations for all figures but figure 4. Plotting of all figures but figure 4. Drafting of paper and contributed to the discussion in the process of writing the paper. Paper IV: I. V. Zozoulenko and M. Evaldsson. Quantum antidot as a controllable spin injector and spin filter. Applied Physics Letters 85, 3136 (2004) Author’s contribution: Calculations and plotting of figure 2. Paper V: Hengyi Xu, T. Heinzel, M. Evaldsson and I. V. Zozoulenko. Resonant reflection at magnetic barriers in quantum wires. Physical Review B 75, 205301 (2007) Author’s contribution: Providing source code and technical support to the calculations. Contributed to the discussion in the process of writing the paper. Paper VI: Hengyi Xu, T. Heinzel, M. Evaldsson and I. V. Zozoulenko. Magnetic barriers in graphene nanoribbons: Theoretical study of transport properties Physical Review B 77, 245401 (2008) Author’s contribution: Part in developing the method and writing the source code. Contributed to the discussion in the process of writing the paper. Paper VII: M. Evaldsson, I. V. Zozoulenko, Hengyi Xu and T. Heinzel. Edge disorder induced Anderson localization and conduction gap in graphene nanoribbons. Submitted. Author’s contribution: Calculations for all figures but figure 2. Plotting of all figures. Drafting of paper and contributed to the discussion in the process of writing the paper. vii viii Contents Abstract iii Preface v List of publications vii 1 Introduction 1 2 Mesoscopic physics 2.1 Heterostructures . . . . . . . . . . . . . . . . . . 2.1.1 Two-dimensional heterostructures . . . . 2.1.2 (Quasi) one-dimensional heterostructures 2.2 Graphene . . . . . . . . . . . . . . . . . . . . . . 2.2.1 Basic properties: experiment . . . . . . . 2.2.2 Basic properties: theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 3 4 7 9 10 11 3 Transport in mesoscopic systems 3.1 Landauer formula . . . . . . . . . 3.1.1 Propagating modes . . . . 3.2 Büttiker formalism . . . . . . . . 3.3 Matching wave functions . . . . . 3.3.1 S-matrix formalism . . . . 3.4 Magnetic fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 15 15 17 19 20 20 . . . . . . . . . 23 23 24 26 27 28 28 29 30 31 4 Electron-electron interactions 4.0.1 What’s the problem? . . 4.1 The Hubbard model . . . . . . 4.2 The variational principle . . . . 4.3 Thomas-Fermi model . . . . . . 4.4 Hohenberg-Kohn theorems . . . 4.4.1 The first HK-theorem . 4.4.2 The second HK-theorem 4.5 The Kohn-Sham equations . . . 4.6 Local Density Approximation . ix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4.7 4.8 Local Spin Density Approximation . . . . . . . . . . . . . . . . Brief outlook for DFT . . . . . . . . . . . . . . . . . . . . . . . 32 33 5 Modelling 5.1 Tight-binding Hamiltonian . . . . . 5.1.1 Mixed representation . . . . . 5.1.2 Energy dispersion relation . . 5.2 Green’s function . . . . . . . . . . . 5.2.1 Definition of Green’s function 5.2.2 Dyson equation . . . . . . . . 5.2.3 Surface Green’s function . . . 5.2.4 Computational procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 35 36 37 39 39 40 41 43 6 Comments on 6.1 Paper I . 6.2 Paper II . 6.3 Paper III 6.4 Paper IV 6.5 Paper V . 6.6 Paper VI 6.7 Paper VII . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 47 48 50 51 52 54 55 papers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . x . . . . . . . . . . . . . . . . . . . . . . . . . . . . Chapter 1 Introduction During the second half of the 20th century, the introduction of semiconductor materials came to revolutionise modern electronics. The invention of the transistor, followed by the integrated circuit (IC) allowed an increasing number of components to be put onto a single silicon chip. The efficiency of these ICs has since then increased several times, partly by straightforward miniaturisation of components. This process was summarised by Gordon E. Moore in the now famous “Moore’s law”, which states that the number of transistors on a chip doubles every second year1 . However, as the size of devices continue to shrink, technology will eventually reach a point when quantum mechanical effects become a disturbing factor in conventional device design. From a scientific point of view this miniaturisation is not troubling but, rather, increasingly interesting. Researchers can manufacture semiconductor systems, e.g., quantum dots or wires, which are small enough to exhibit pronounced quantum mechanical behaviour and/or mimic some of the physics seen in atoms. In contrast to working with real atoms or molecules, experimenters can now exercise precise control over external parameters, such as confining potential, the number of electrons, etc.. In parallel to this novel research field, a new applied technology is emerging – “spintronics” (spin electronics). The basic idea of spintronics is to utilise the electron spin as an additional degree of freedom in order to improve existing devices or innovate entirely new ones. Existing spintronic devices are built using ferromagnetic components – the most successful example to date is probably the read head in modern hard disk drives (see e.g., [94]) based on the Giant Magneto Resistance effect[10, 15]. Because of the vast knowledge accumulated in semiconductor technology there is an interest to integrate future spintronic devices into current semiconductor ones. This necessitates a multitude of questions to be answered, e.g.: 1 Moore’s original prediction made in 1965 essentially states (here slightly reformulated), that “the number of components on chips with the smallest manufacturing costs per component doubles roughly every 12 months”[74]. It has since then been revised and taken on several different meanings. 2 Introduction can spin-polarised currents be generated and maintained in semiconductor materials, does spin-polarisation appear spontaneously in some semiconductor systems? Experiments, however, are not the only way to approach these new and interesting phenomena. The continuous improvement of computational power together with the development of electron many-body theories such as the Density Functional Theory provide a basis for investigating these questions from a theoretical/computational point of view. Theoretical work may explain phenomena not obvious from experimental results and guide further experimental work. Modelling of electron transport in lateral semiconductor nanostructures and graphene is the main topic of this thesis. Chapter 2 Mesoscopic physics The terms macroscopic and microscopic traditionally signify the part of the world that is directly accessible to the naked eye (e.g., a flat wall), and the part of the world which is to small to see unaided (e.g., the rough and weird surface of the flat wall in a scanning electron microscope). As the electronic industry has progressed from the macroscopic world of vacuum tubes towards the microscopic world of molecular electronics, the need to name an intermediate region has come about. This region is now labelled mesoscopic, where the prefix derives from the Greek word “mesos”, which means ‘in between’. Mesoscopic systems are small enough to require a quantum mechanical description but at the same time too big to be described in terms of individual atoms or molecules, thus ‘in between’ the macroscopic and the microscopic world. The mesoscopic length scale is typically of the order of: • The mean free path of the electrons. • The phase-relaxation length, the distance after which the original phase of the electron is lost. Depending on the material used, the temperature, etc., these lengths and the actual size of a mesoscopic system could vary from a few nanometres to several hundred micrometres [22]. This chapter introduces manufacturing techniques, classification and general concepts of low-dimensional and semiconductor systems. The first sections introduce laterally defined systems in heterostructures; relevant for papers I-V. In the last section we consider graphene, relevant for papers VI-VII. 2.1 Heterostructures A heterostructure is a semiconductor composed of more than one material. By mixing layers of materials with different band gaps, i.e. band-engineering, it is possible to restrict electron movement to the interface (the heterojunction) 4 Mesoscopic physics between the materials. This is typically the first step in the fabrication of low-dimensional devices. Since any defects at the interface will impair electron mobility through surface-roughness scattering, successful heterostructure fabrication techniques must yield a very fine and smooth interface. Two of the most common growth methods are molecular beam epitaxy (MBE) and metalorganic chemical vapour deposition (MOCVD). In MBE, a beam of molecules is directed towards the substrate in an ultra high vacuum chamber, while in MOCVD a gas mixture of the desired molecules are kept at specific temperatures and pressures in order to promote growth on a substrate. Both these techniques allow good control of layer thickness and keep impurities low at the interface. In addition to the problem with surface-roughness, a mechanical stress due to the lattice constant mismatch between the heterostructure materials causes dislocations at the interface. This restricts the number of useful semiconductors to those with close/similar lattice constants. A common and suitable choice, because of good lattice constant match and band gap alignment (see figure 2.1), is to grow Alx Ga1−x As (henceforth abbreviated AlGaAs, with the mixing factor x kept implicit) on top of a GaAs substrate. The position and relative size of the band gap in a GaAs-AlGaAs heterostructure is schematically shown in figure 2.1. At room temperature the band gap for GaAs is 1.424eV[82], while the band gap in AlGaAs depends on the mixing factor x and can be approximated by the formula[82] Eg (x) = 1.424 + 1.429x − 0.14x2 [eV] 0 < x < 0.441 , (2.1) i.e., it varies between 1.424-2.026eV. 2.1.1 Two-dimensional heterostructures In GaAs-AlGaAs heterostructures, a two-dimensional electron system is typically created by n-doping the AlGaAs (figure 2.2 or the left panel of figure 2.5). Some of the donor electrons will eventually migrate into the GaAs. These electrons will still be attracted by the positive donors in the AlGaAs, but be unable to go back across the heterojunction because of the conduction band discontinuity. Trapped in a narrow potential well (see figure 2.2), their energy component in this direction will be quantised. Because the potential well is very narrow (typically 10nm[24]), the available energy states will be sparsely spaced, and at sufficiently low temperatures all electrons will be in the same (the lowest) energy state with respect to motion perpendicular to the interface. I.e., electrons are free to move in the plane parallel to the heterojunction, but restricted to the same (lowest) energy state in the third dimension. In this sense we talk about a two-dimensional electron gas (2DEG). 1 For larger x the smallest band gap in AlGaAs is indirect. 2.1. Heterostructures 5 AlGaAs GaAs Conduction band Band gap Valence band Figure 2.1: Band gap for AlGaAs (left) and GaAs (right) schematically. For these materials their corresponding band gap results in a straddling alignment, i.e., the smaller band gap in GaAs is entirely enclosed by the larger band gap in AlGaAs. Energy Original n-AlGaAs band 2DEG Space Conduction band donors n-AlGaAs Original GaAs band GaAs AlGaAs (spacer-layer) Figure 2.2: Two-dimensional electron gas seen along the confining dimension. Remote or modular doping, i.e., to place the donors only in the AlGaAs layer, prevent the electrons at the heterojunction to scatter against the positive donors. Usually an additional layer of undoped AlGaAs is grown at the interface as a spacer layer. This will, at the expense of high electron density, further shield the electrons from scattering. Densities in 2DEGs typically varies between 1 − 5 × 1015 m−1 though values as low as 5 × 1013 m−1 has been reported[41]. Density of States A simple but yet powerful characterisation of a system is given by its density of states (DOS), N (E), where N (E)dE is defined as the number of states in the energy interval E → E + dE. For free electrons it is possible to determine the 2-dimensional DOS, N2D (E), exactly. This is done by considering electrons 6 Mesoscopic physics ky 1 0 0 1 1 0 0 1 00 11 0 1 0 1 00 11 0 k 1 0 1 1 0 0 1 0 01 11 k+dk 1 01 000 1 2π 0 1 kx L 1 0 111111 000000 00 11 0 100 11 n2D (E) = m ~2 π n(E) 2π L E (a) 2D k-space (b) 2D density of states Figure 2.3: (a) Occupied and unoccupied (filled/empty circles) states in 2D k-space. (b) The two-dimensional density of states, n2D (E). situated in a square of area L2 and letting L → ∞. With periodic boundary conditions the solutions are travelling waves, φ(r) = eik·r = ei(kx x+ky y) , (2.2) where k = (kx , ky ) = 2πm 2πn , Lx Ly = 2πm 2πn , L L m, n = 0, ±1, ±2 . . . . (2.3) Plotting these states in the 2-dimensional k-space, figure 2.3a, we recognise 2 , hence the density of states is that a unit cell has the area 2π L N2D (k) = 2 (2π)2 L2 −1 = L2 2π 2 (2.4) where a factor two is added to include spin. Defining the density of states per unit area, N2D (k) 1 n2D (k) = = 2, (2.5) L2 2π we get a quantity that is defined as L → ∞. In order to change variable from n2D (k) to n2D (E), we look at the annular area described by k and k+dk in figure 2.3a. This area is approximated by 2πkdk and, thus, contains n2D (k)2πkdk states. The number of states must of course be the same whether we express it in terms of wave vector k or energy E, i.e., dE 1 dk = 2πkn2D (k)dk = 2πk 2 dk dk 2π k dE −1 ⇔ n2D (E) = . π dk n2D (E)dE = n2D (E) (2.6) (2.7) 2.1. Heterostructures 7 For free electrons, with E(k) = ~2 k2 , 2m where k = |k| (2.8) n2D (E) = m , ~2 π (2.9) we finally arrive at shown in figure 2.3b. A similar derivation for a one dimensional system yields the density of states per unit length, r 1 2m . (2.10) n1D (E) = ~π E 2.1.2 (Quasi) one-dimensional heterostructures There are several techniques to further restrict the 2DEG. Two common approaches are chemical etching, or to put metallic gates on top of the sample. Figures 2.4(a)and (b) show a quantum wire and a quantum dot, respectively, defined by a potential applied to top gates. Side gate Side gate Top gate 2DEG (a) (b) Figure 2.4: Schematic figure of (a) quantum wire, (b) quantum dot, defined by potentials applied to metalic gates. The layers in the heterostructure are, from bottom to top, substrate, spacer, donor and cap layer. Etching By etching away part of the top dopant layer, the electron gas is located to the area beneath the remaining dopant as schematically illustrated in figure 2.5. Metallic gates A second alternative is to deploy metallic gates on top of the surface as shown in figure 2.4(a)-(b). By applying a negative voltage to the gates, the 2DEG will be depleted beneath them. This technique allows experimenters to control the approximate size of the system by changing the applied potential during the experiment. 8 Mesoscopic physics n-AlGaAs AlGaAs n-AlGaAs AlGaAs 11111111111 00000000000 GaAs 2DEG 2DEG GaAs Figure 2.5: Left: 2DEG in AlGaAs-GaAs heterostructure. Right: Etching away the AlGaAs everywhere but in a narrow stripe results in a one-dimensional quantum wire beneath the remaining AlGaAs. Subbands in quasi one-dimensional systems The transversal confinement in a quasi one-dimensional system, such as the quantum wire in figure 2.4(a), results in a quantisation of energies in this dimension. Denoting the confining potential U (y) (i.e., electrons are free along the x-axis), the Schrödinger equation 2 ∂ ∂2 ~2 + + U (y) Ψ(x, y) = EΨ(x, y), (2.11) − 2m ∂x2 ∂y 2 can, by introducing Ψ(x, y) = ψ(x)φ(y), be separated into − ~2 2m − d2 dy 2 ~2 d2 ψ(x) = Ex ψ(x) 2 2m dx + U (y) φ(y) = Ey φ(y). In the simple case of a hard wall potential U (y), ( 0, 0<y<w U (y) = ∞, y < 0 ∪ w < y, (2.12) (2.13) (2.14) eq. (2.12) and (2.13) yields the solutions Ψn (x, y) = ψ(x)φn (y) = eikx x sin πny w (2.15) and eigenenergies E = Ex + Ey,n ~2 = 2m kx2 + πn 2 w = ~2 2 2 kk + π 2 k⊥ . 2m (2.16) The total energy E is the sum of a continuous part Ex , and a discreet part Ey,n with corresponding continuous (kk ) and and discrete (k⊥ ) wave-vectors. 2 , consists Hence, the energy dispersion relation, E(k) = E(kk + k⊥ ) ∼ kk2 + k⊥ of subbands, as depicted in figure 2.6. 2.2. Graphene 9 E3 φ1 φ2 φ3 Energy E2 E1 EF kk (a) (b) Figure 2.6: (a) The three lowest propagating nodes, φi (y), in a hard wall potential quantum wire. (b) Corresponding subbands to φ1 (y)–φ3 (y) seen in (a). The Fermi energy is indicated by EF . 2.2 Graphene The discovery of the carbon fullerenes in 1985[58] paved the way for extensive research into the allotropes of carbon. Today, a wide range of structures such as buckyballs, carbon nanotubes, or a combination thereof (carbon nanobuds) are known. Basically, these structures consist of a single or several layers of graphite wrapped up in some configuration. Interestingly, the existence of a single layer of graphite, not wrapped up, was long thought to be thermodynamically unstable[61, 62, 88, 89]. The initial reports of such structures in 2004[79] was therefore somewhat unexpected. Single-layered graphite is now known as graphene and has been the focus for intense research during the past few years. The reason graphene does not become unstable, as theory predicts for 2D structures, seems to be the result of corrugations which stabilizes the sheet[72]. It is thus more correct to view graphene as a two-dimensional structure in a three-dimensional space rather than a strictly two-dimensional structure. Although graphene hasn’t been experimentally studied for more than a few years it has been an issue of theoretical interest for a long time. This is partly because the properties of various carbon based materials such as graphite, carbon nanotubes, etc., derive from the properties of graphene and partly because some of the properties of graphene are rather extraordinary, making graphene an interesting “toy-model” for theoreticians. The honeycomb structure of graphene (figure 2.7) makes the charge carriers in the lattice mimic relativistic particles which can be described by the Dirac equation[29, 103, 105]. This causes a number of effects not expected in non-relativistic systems to be present; an anomalous quantum hall effect[78, 124], the Klein paradox2 [54], the presence of a minimum conductivity as charge carrier concentration is 2 Where charge carriers can tunnel through an arbitrary high energy barrier with unity probability 10 Mesoscopic physics depleted[78, 113]. A-lattice B-lattice acc { Figure 2.7: The graphene lattice can be described as two sub-lattices, A and B. Primitive vectors for sub-lattice A are indicated by arrows. The interatomic distance acc is roughly 1.42Å. 2.2.1 Basic properties: experiment The discovery of graphene had an illusive air of simplicity – repeated peeling of a graphite sample with adhesive tape or rubbing it against a surface (i.e. “drawing”), resulted in flakes a single atom layer thick[78, 79]. Yet the approach is nothing but simple; monolayers are in minority of the flakes found and conventional techniques for identifying such 2D structures are either too slow for a random search (atomic force microscopy), lack clear signatures (transmission electron microscopy) or are invisible under most circumstances (optical microscopy)[78]. Only by preparing the flakes on a proper substrate, such as a 300nm thick SiO2 , they become visible in an optical microscope due to interference[40]. Mechanical extraction of graphene from graphite can produce micrometer sized samples which are sufficient for scientific purposes. For any commercial use other techniques, such as epitaxial growth[31], need to be further developed. Mobility measurements of graphene show the extraordinary qualities of its crystal structure. Room temperature mobilities µ of the order 10,000 cm2 V−1 s−1 [79] has been reported and experiments on suspended graphene at low temperatures reported a peak mobility at 230,000cm2 V−1 s−1 [17]. The mobilities remain relatively high even when the graphene sheet is doped[104]. Another intriguing property of graphene is the presence of a minimum conductivity σmin at the Dirac point. Instead of a metal to insulator transition as the charge density decreases the conductivity stabilizes at σmin ∼ 4e2 /h[40]. Although the effect isn’t unexpected for systems governed by the Dirac equation[37, 54], the theoretically expected value is smaller by a factor of π, σmin,theor = 4e2 /πh. Occasional measurements approaches the theoretical limit but it is still unclear what physics determines σmin . It has been suggested that measured σmin depends on the width/length ratio of the measured samples[20, 113] and that wide and short ribbons would give σmin closer to the theoretical value. Another suggestion is that at low electron density, 2.2. Graphene 11 ky M K b2 b1 Γ kx E[eV] 5 Dirac point 0 -5 M (a) Brillouin zone Γ Wave vector K M (b) Tight-binding dispersion relation Figure 2.8: (a) Shaded area shows the Brillouin zone in the reciprocal lattice. b1 , b2 are the reciprocal lattice vectors. (b), Tight-binding dispersion relation for ǫ2p = γ0 =0, γ1 =-2.7eV along contour in (a). Fermi level is at 0 eV charges in graphene moves in a poodle like electron-hole landscape where the underlying physics becomes more complicated than expected[121]. 2.2.2 Basic properties: theory Graphene is a single layer of carbon atoms in a honeycomb-lattice, figure 2.7. Most of the “exotic” properties of graphene derive from its hexagonal lattice which results in a linear dispersion relation close to the fermi level instead of the parabolic dispersion typical for solid state systems. The lattice can be described as two triangular sub-lattices, figure 2.7, where the interatomic distance acc roughly equals 1.42Å. Using the tight-binding approximation the dispersion relation can be computed from the secular equation[48, 98, 102, 117] HAA − ESAA HAB − ESAB (2.17) HBA − ESBA HBB − ESBB = 0 Here, HAA = hΦA |H|ΦA i where Φj (k, r) is the Bloch function for site j = A, B, X 1 Φj √ = eik·R φj (r − R) N R (2.18) with φj being the wave function localized at site j and R the position of the lattice points. Similarly, Sjj ′ = hΦj |Φj ′ i is the overlap matrix and E the eigenenergies. The solution of eq. (2.17) is of the form[48, 98, 102] E(k) = ǫ2p ± γ1 |g(k)| 1 ± γ0 |g(k)| (2.19) 12 Mesoscopic physics Figure 2.9: (a) Example of graphene nano ribbon with zigzag edge (Z-GNR). The GNR extends towards infinity in the x-direction. The width of the Z-GNR in the y-direction is characterized by the number of zigzag chains, Nz , indicated as solid dots across the GNR. Here Nz =6. (b)-(c) Tight-binding dispersion relation and transmission through a Z-GNR with (b) Nz =7 and (c) Nz =8 transversal sites. Within the tight-binding approximation Z-GNR:s are metallic for all Nz . (d) Example of graphene nano ribbon with armchair edge (A-GNR). The width of A-GNR:s is characterized by the number of dimer lines, Na , indicated as solid dots across the GNR. Here Na =9. (e)-(f ) Tight-binding dispersion relation and transmission through an A-GNR with (e) Na =7 and (f) Na =8 transversal sites. AGNR:s are metallic only if Na =3p + 1, p being an integer. where[1] v u u |g(k)| = t1 + 4 cos ! √ 3ky a kx a kx a cos + 4 cos2 2 2 2 (2.20) √ and a = 3acc ∼ 2.46nm being the length of the reciprocal lattice vector. Fitting of the parameters ǫ2p , γ0 and γ1 around the Dirac point, i.e., the 2.2. Graphene 13 K-point in figure 2.8(a),(b)) gives[98] e2p = 0 −2.5eV < γ1 < −3.0eV (2.21) γ0 < 0.1eV. γ0 is related to the the overlap between nearest neighbours, hφA |φB i, and is sometimes set to zero to simplify the model further. Figure 2.8(b) shows the tight-binding dispersion relation along the contour in the Brillouin zone in fig. 2.8(a). The tight-binding binding approximation generally deviates significantly from ab initio calculations away from the K-points[98]. Close to the K-point the dispersion relation is linear as for relativistic particles with an effective “light” velocity of vF ∼ c/300[86]. Graphene nano ribbons For application purposes of graphene both theoretical and experimental focus has been on graphene nano ribbons (GNR:s), i.e., stripes of graphene cut out of larger sheets (figure 2.9(a),(d)). Hopefully GNR:s will allow the extraordinary properties of graphene, such as high mobility at ambient temperature and high degree of doping, to be used in conjunction with existing semiconductor components. Graphene nano ribbons has been lithographically patterned down to widths of ∼20nm[20, 47] and chemically grown with a width less than 10nm[66]. Modelling of GNR:s indicate that their electronic properties are determined by the width of the wire and the edge structure[39, 98]. Typically, two kinds of edge structures are considered; zigzag edge graphene nano ribbons (Z-GNR), as in figure 2.9(a), and armchair edge graphene nano ribbons (A-GNR), as in figure 2.9(d). The width of the ribbon is characterised by a number Nz (Na ) for Z-GNR (A-GNR) which is defined by the number of sites across the ribbon along to the paths given by solid dots figure 2.9(a),(d). Within the tight-binding approximation Z-GNR:s are metallic for all widths, e.g., figure 2.9(e),(f). For A-GNR only certain widths, Na =3p+1, p being an integer are metallic (c.f., figure 2.9(e) and (f)). In contrast, Ab initio DFT calculations of GNR:s yields a band gap for all Z-GNR and A-GNR[109]. The difference between tight-binding and DFT is due to effects at the edge of the ribbon, and the discrepancy decreases with increasing ribbon width. The overall trend of the band gaps is a 1/W dependence, figure 2.10. Because graphene itself is a gapless semiconductor, very narrow GNR:s has been proposed as a way to engineer a band gap in graphene components. Conductance measurements of GNR:s show a 1/W scaling of the band gap but surprisingly no dependence on the crystal direction of the GNR[47]. Although there is no consensus on the lack of crystal dependence in conductance measurement yet, possible explanations includes rough edges[20, 47, 66, 96], atomic scale impurities[20], coulomb blockaded transport[108]. 14 Mesoscopic physics 50 Energy gap (units of “t”) 0.5 Width [Na ] 100 150 200 0.4 0.3 0.2 0.1 0 0 5 10 15 Width [nm] 20 25 Figure 2.10: Energy gap of armchair GNR versus ribbon width within the tightbinding approximation. Metallic A-GNR:s, with Na =3p+1 where p ∈ N, are excluded. Chapter 3 Transport in mesoscopic systems Understanding transport is central to the study of mesoscopic systems. Often transport characteristics work as a flexible tool to probe the electron states inside a system – well-known examples are the Kondo effect[57] and the 0.7-anomaly[112]. This chapter introduces some basics of ballistic transport in mesoscopic systems, that is, systems where electrons pass through the conductor without scattering and remain phase coherent. In brief, we will start by deriving the Landauer formula which describes transport for a system with only two leads connected. This is then generalised within the Büttiker formalism to handle an arbitrary number of leads. A key characteristic for the quantum conductance in these expressions is the transmission probability Tn,β←m,α , i.e., the probability for an electron in mode α and lead m to end up in another mode β and lead n. Finally we will look at effects on transport when a perpendicular magnetic field is applied. 3.1 Landauer formula 3.1.1 Propagating modes For simplicity we start with the case of a single propagating mode (see figure 2.6a), i.e., where only the lowest subband in both the left and right lead is occupied. Figure 3.1 shows electrons surrounding a barrier in a 1-dimensional system with some bias eV applied. Though in principle, the Fermi energy is only defined at equilibrium, we will assume this bias to be small enough to yield a near-equilibrium electron distribution that is characterised by a quasiFermi level, EF l and EF r respectively. To find an explicit expression for the net current across the barrier due to the bias we first consider the electrons approaching the barrier from the left. In an infinitesimal momentum interval 16 Transport in mesoscopic systems EF l eV Ul EF r Ur Figure 3.1: Schematic view of a 1D barrier surrounded by electrons. A small bias eV shifts the (quasi-)Fermi levels in left (EF l ) and right (EF r ) lead. dk around k, the transmitted current is Il (k)dk = 2en1D (k)v(k)Tr←l (k)f E(k), EF l dk (3.1) where a factor 2 is added to include spin, e the electron charge, n1D (k) is the one-dimensional density of states, v(k) the velocity of the electrons, Tr←l (k) the probability that an electron passes the barrier and f E(k), EF l the FermiDirac distribution. Using eq. (2.10) for the 1D-DOS and integrating both sides yields Il = Z ∞ 2e 0 / / dk 1 v(k)Tr←l (k)f E(k), EF l dk = dk = dE 2π dE Z ∞ dk 2e v(E)Tr←l (E)f E, EF l dE. (3.2) = 2π 0 dE Recognising the group velocity as v = bottom of the left lead brings Il = 2e 2π Z ∞ v(E)Tr←l (E)f E, EF l Ul dω dk = 1 dE ~ dk and integrating from the 1 dE ~v(E) Z 2e ∞ = Tr←l (E)f E, EF l dE. (3.3) h Ul There is, except for the different Fermi level and opposite direction of flow, a similar current Ir from the right to left lead. If the bias is small, the reciprocity relation Tl←r (E)=Tr←l (E) holds (see e.g. [23]), and we can skip the indices on the transmission coefficient. The net current becomes Z h i 2e ∞ I = Il + Ir = T (E) f E, EF l − f E, EF r dE. (3.4) h Ul In the case of very low bias eV , i.e., in the linear response regime, the FermiDirac functions in eq. (3.4) can be Taylor expanded around EF = 21 (EF l + EF r ) according to 1 1 f (E, EF l ) − f (E, EF l ) = f (E, EF + eV ) − f (E, EF − eV ) 2 2 3.2. Büttiker formalism ∂f (E, EF ) ∂f (E, EF ) + (eV )2 ≈ −eV , (3.5) ∂EF ∂E = eV resulting in I= ⇔ 17 2e h Z ∞ Ul ∂f (E, EF ) T (E) −eV dE ∂E Z I 2e2 G= = V h ∞ Ul ∂f (E, EF ) T (E) − ∂E (3.6) dE. (3.7) ∂f is replaced by At very low temperatures in the linear response regime, − ∂E the Dirac delta function δ(E − EF ) and evaluating the integral in eq. (3.7) gives the famous Landauer formula for quantum conductance G= 2e2 I = T (EF ). V h (3.8) The result in eq. (3.8) above is readily extended to the case of several propagating modes in the leads. An electron incoming towards a scatterer in a specific mode α might be transmitted to some mode β in the opposite lead or reflected to some mode β in the same lead. By summing over all possible modes the total conductance is found as G= 2e2 X Tβ←α . h (3.9) α,β 3.2 Büttiker formalism Büttiker formalism describes transport in systems with more than two leads. A typical example is the three lead system where a net current flows between two leads (current probes), and the third lead is used to measure the potential (voltage probe). We will assume fairly low temperature and bias. In this limit the Fermi-Dirac functions in eq. (3.4) may be replaced by step functions, 2e I= h Z ∞ Ul h i T (E) Θ EF l − E − Θ EF r − E dE = 2e h Z EF l T (E)dE. (3.10) EF r With low bias we assume T (E) to be constant in the interval EF r < E < EF l , thus the integral in eq. (3.10) evaluates to I= 2e T × (EF l − EF r ). h (3.11) 18 Transport in mesoscopic systems We now consider a system connected to an arbitrary number of leads (indexed by m and n) and propagating modes (α and β). Im denotes the total net current in lead m. It is the sum of the net incident current P Imm (Iincoming − Iref lected ) and the current injected from all other leads, n6=m Im←n . With the lowest quasi-Fermi level in all leads denoted E0 we can write 2e X X Im←n = − Tm,α←n,β × (En − E0 ). (3.12) h n6=m α,β The minus sign indicates a current away from the system and Tm,α←n,β is the transmission probability from lead n–mode β to lead m–mode α. We will abbreviate this by the notation X (3.13) T̄m←n = Tm,α←n,β , α,β and in a similar way R̄m = X (3.14) Rm,α←m,β α,β for the reflection Rm,α←m,β from from mode β to α in lead m. Lead m carries Nm propagating modes hence the net incident current in lead m is 2e Imm = Nm − R̄m (Em − E0 ). (3.15) h Furthermore, conservation of flux requires that X X Nm = R̄m + T̄n←m = R̄m + T̄m←n (3.16) m6=n n6=m where the last equality follows from the sum rule derived in eq. (3.30)–(3.35). The total current in lead m can be written as 2e X 2e Im = Imm + Im←n = Nm − R̄m (Em − E0 ) − T̄m←n (En − E0 ) h h n6=m (3.17) X X 2e 2e = T̄m←n (Em − E0 ) − T̄m←n (En − E0 ) h h n6=m n6=m (3.18) 2e X = T̄m←n Em − T̄m←n En , h (3.19) n6=m where eq. (3.16) has been used in the second step. Noting that Em = eVm we finally get 2e2 X T̄m←n (Vm − Vn ) . (3.20) Im = h n6=m The currents in a multi-lead system are then found from solving the system of equations defined by equation (3.20). 3.3. Matching wave functions 19 Ce−ik2 x Aeik1 x Deik2 x Be−ik1 x V0 x=0 Figure 3.2: Incoming and outgoing electrons at a potential step. 3.3 Matching wave functions We now turn to a simple example for computing the transmission probability T across a potential step, see figure 3.2. k1 and k2 is the wave vector for incident and transmitted electrons, respectively. The solution to the Schrödinger eq. for electrons with energy E > V0 is ( Aeik1 x + Be−ik1 x x < 0 Ψ(x) = (3.21) Ce−ik2 x + Deik2 x x > 0. The requirement that the wave function and its derivative is continuous at x = 0 yields the system of equations A+B = C +D (3.22) k1 (A − B) = k2 (D − C). (3.23) Solving for the outgoing amplitudes, B and D, results in 1 B k1 − k2 2k1 A = . D 2k1 k2 − k1 C k1 + k2 (3.24) E.g., for a free electron incoming from the left with unit amplitude (A = 1, C = 0), the transmitted amplitude (D = t2←1 ) becomes t2←1 = 2k1 , k1 + k2 (3.25) and the reflected amplitude (B = r1←1 ) r1←1 = k1 − k2 . k1 + k2 (3.26) Because the velocities in the left and right lead may differ, the flux transmission coefficient becomes v2 (k2 ) |t2←1 |2 , (3.27) T2←1 (k1 , k2 ) = v1 (k1 ) where v2 and v1 are the group velocities for the electron. 20 3.3.1 Transport in mesoscopic systems S-matrix formalism By generalizing the notation in the preceding section we can now prove the conservation of flux used in eq. (3.16). The matrix equation (3.24) above relates the amplitude of the outgoing electrons to the amplitude of the incoming electrons. The current into the system for mode i is given by Ii = −vi (k)|Ai (k)|2 where vi (k) is the group velocity and Ai (k) the amplitude. It is convenient to define the current amplitude, by p ai (k) =Ai (k) −v(k) (current amplitude for incoming mode i) p bi (k) =Bi (k) v(k) (current amplitude of outgoing mode i), such that the current carried by mode i simply is Ii = |ai (k)|2 . We now generalise eq. (3.24) to handle an arbitrary system by defining the S-matrix (scattering matrix)[22], b̄ = Sā, (3.28) which, for each energy E, relates the incident current amplitudes a = (a1 , a2 , . . . , an ) to the outgoing current amplitudes b = (b1 , b2 , . . . , bn ). The transmission probability now equals Tm←n (E) = |smn |2 , (3.29) where smn is a matrix element in the S-matrix. Now, current conservation requires that X X Iin = e |ai |2 = e |bi |2 = Iout (3.30) i i or in matrix notation † ā† ā = b̄ b̄ ⇔ (3.31) † † ā ā = (Sā) Sā (3.32) ⇔ ā† ā = ā† S† Sā (3.33) ⇒ S† S = I = SS† , (3.34) eq. (3.28) i.e., S is unitary and we arrive at X X X X Tm←n = |smn |2 = 1 = |snm |2 = Tn←m . m 3.4 m m (3.35) m Magnetic fields Magnetic field is usually incorporated into the Hamiltonian through a vector potential A defined by the equation B=∇×A (3.36) 3.4. Magnetic fields 21 n(E) (b) (a) (c) 1 2 ~ωc 3 2 ~ωc Energy 5 2 ~ωc Figure 3.3: (a)–dashed line: 2D density of states with B=0T. (b)–solid line: Ideal 2D density of states in a magnetic field. (c)–dotted line: Broadened 2D density of states in a magnetic field. as 1 (3.37) (−i~∇ + eA)2 + V. 2m Since ∇ × ∇ξ = 0 for any function ξ, there is a possibility to choose different gauges, A→A+∇ξ, without changing the physical properties of the system. Some common choices for a magnetic field Bẑ is Landau gauge A=(−By,0,0) or A=(0,Bx,0), and the circular symmetric gauge A= 12 (−By, Bx,0). We will now consider A=(0,Bx,0), and the Schrödinger equation for a 2DEG becomes H= HΨ = ~ 2 ie~Bx ∂ (eBx)2 − ∇x,y − + Ψ(x, y) = EΨ(x, y). 2m m ∂y 2m (3.38) Trying a solution in the form Ψ(x, y) = u(x)eiky we end up with an equation only in x, 2 ! ~2 d2 mωc2 ~k − + +x u(x) = ǫu(x) (3.39) 2m dx2 2 eB with ωc2 = eB m being the cyclotron frequency. This is the Schrödinger equation for a harmonic oscillator so we can write down the eigenenergies as 1 ǫn,k = (n − )~ωc , 2 n = 1, 2, 3, . . . (3.40) The energies in eq. (3.40) are independent of k, hence the density of states – which for zero magnetic field is constant (eq. (2.9) or fig. 3.3) – now collapses to very narrow stripes, equidistantly located as shown in figure 3.3. These collapsed states are called Landau levels and contain all the degenerate states with the same k. In an actual system they are broadened due to scattering. Next we consider a hard wall potential electron wave-guide of the width L, oriented along the y-axis. Trying the same solution as above, Ψ(x, y) = 22 Transport in mesoscopic systems V (x) 10 (a) 20 (b) (c) 15 5 10 5 0 -50 0 nm 0 50 -50 0 nm 50 Figure 3.4: (a): Potentials (black) and wave functions (gray) for B=1T (solid lines) and B = 3T (dashed lines) in a quantum wire, k=0. (b): Potentials (black) and wave functions (gray) for k=0.05nm−1 (solid lines) and k=0.25nm−1 (dashed lines) in a quantum wire, B = 3T . (c): Classical skipping orbits in a straight quantum wire. u(x)eiky , we arrive at mωc2 ~2 d2 + − 2m dx2 2 ~k +x eB 2 ! + V (x) u(x) = ǫu(x) (3.41) where the only difference from eq. (3.39) is the confining potential V (x), ( 0 |x| ≤ L2 V (x) = . (3.42) ∞ |x| > L2 This equation do not have an analytical solution but we may still draw some qualitative conclusions from eq. (3.41). The electrons are confined by the hard wall potential V (x) and the parabolic potential which depends on k and B. With k=0, the parabolic potential is determined only by the magnetic field B (through wc ), and an increasing magnetic field will steepen the potential thereby confining the electrons to the middle of the wire as shown in the leftmost panel of figure 3.4. For |k| >0 the vertex of the parabola is displaced, and for large |k|’s the wave function is squeezed against the hard wall confinement (middle panel of figure 3.4). Hence, currents in this system are carried along the edges of the wire by these “edge states”. The classic analogue to these states is the so called skipping orbits shown in the rightmost panel of the same figure. Due to the Lorentz force, an electron in a perpendicular magnetic field moves in a circle and do not contribute to any net current. However, along the edges of the wire electrons will bounce off the walls, resulting in a current in both directions. Chapter 4 Electron-electron interactions Some years after the Schrödinger equation and the basics of quantum mechanics had been formulated, Paul Dirac commented that[26] The underlying physical laws necessary for the mathematical theory of a large part of physics and the whole of chemistry are thus completely known, and the difficulty is only that the exact application of these laws leads to equations much too complicated to be soluble. The results from quantum mechanics were impressive already by 1929, but there was also an awareness of the huge computational effort needed to solve some (most) problems. The difficulties comes from electron-electron interactions (e-e interactions) which poses a intractable many-body problem under most circumstances. Today, a wide range of approximations exists to handle e-e interactions. This chapter introduces and discusses the methods dealing with e-e interactions relevant for the thesis. First we take a look at the Hubbard model, where e-e interactions are considered local, i.e., restricted to the interaction between electrons on the same atom or grid point. Next we introduce the Thomas-Fermi (TF) model as a precursor to Density Functional Theory (DFT). Both TF and DFT replace the many-body wave function Ψ(r1 , r2 , . . .) from the Schrödinger equation with the far simpler electron density, n(r). DFT, or any of its extensions, is today one of the most effective theory for e-e interactions around. It has both a solid theoretical foundation and can in many cases give excellent predictions on electronic properties. 4.0.1 What’s the problem? The properties of a stationary N -electron system can be found by solving the time-independent Schrödinger equation 2 X ~2 X 1 e ĤΨ = − ∇2 + + vext (r) Ψ = EΨ. (4.1) 2m i 2 ′ |ri − ri′ | i i6=i 24 Electron-electron interactions Ψ=Ψ(r1 , r2 , . . . , rN ) is the N -electron wave function in three dimensions, the terms in the brackets of eq. (4.1) are the kinetic energy for the i:th electron, the e-e interaction and finally the external potential vext (r). Because of the e-e interaction the coordinates r1 , r2 , . . . , rN are coupled and a direct solution for increasing N is a very difficult many-body problem. 4.1 The Hubbard model The approach taken by John Hubbard in 1963 was to ignore inter-atomic e-e interactions completely. In the model, electrons are only interacting with electrons situated on the same atom (or site) in the system. This approximation is best suited for systems where the electron density is concentrated near the atom nuclei and sparse between the atoms; for example the d-band of transition metals initially considered by Hubbard. For a brief derivation we will follow the same procedure as Hubbard[50], and start with Bloch functions Ψk and energies ǫk calculated in an appropriate spin independent Hartree-Fock potential. The Hamiltonian can then be approximated by H= X ǫk c†kσ ckσ kσ + 1 2 X X k1 k2 k′1 k′2 σ1 σ2 − hk1 k2 |1/r|k′1 k′2 ic†k1 σ1 c†k2 σ2 ck′ σ2 ck′ σ1 2 1 XX 2hkk′ |1/r|kk′ i − hkk′ |1/r|k′ ki νk′ c†kσ ckσ , (4.2) kk′ σ where the momentum sum is over the first Brillouin zone, c†kσ and ckσ are the creation and annihilation operators for the Bloch state k with spin σ (=±1), the e-e interaction terms are Z ψ ∗ (x)ψ ′ (x)ψ ∗ (x′ )ψ ′ (x′ ) k1 k2 k1 k2 hk1 k2 |1/r|k′1 k′2 i = e2 dxdx′ (4.3) ′ |x − x | and νk are the occupation numbers. The first term in eq. (4.2) represents the band energies of the electrons, the second term their interaction and the third term counters double counting. Next the Wannier functions are introduced as a new basis set, 1 X φ(x) = √ ψk (x), (4.4) N k where N is the number of electrons. The Bloch wave function and creation/annihilation operators can be expressed with the Wannier functions as, 1 X ikRj e φ(x − Rj ) ψk (x) = √ N j (4.5) 4.1. The Hubbard model and 25 1 X ikRj ckσ = √ e cjσ , N j 1 X ikRj † c†kσ = √ e ciσ , N j (4.6) Ri being coordinates for lattice site i. The Hamiltonian (eq. 4.2) becomes H= XX i,j Tij c†iσ cjσ σ + 1 XX hij|1/r|klic†iσ c†jσ′ clσ′ ckσ 2 ijkl σσ′ XX − 2hij|1/r|kli − hij|1/r|lki νjl c†iσ c†kσ , (4.7) σ ijkl where Tij = N −1 X ǫk eik(Ri −Rj ) , (4.8) k 2 hij|1/r|kli = e Z φ∗ (x − Ri )φ (x − Rk )φ∗ (x′ − Rj )φ (x′ − Rl ) dxdx′ , |x − x′ | (4.9) and νjl = N −1 X νk eik(Rj −Rl ) . (4.10) k Obviously the main contribution to the e-e interaction comes from the term hii|1/r|iii = U and Hubbard subsequently suggested to ignore the contribution from all other terms. The simplified Hubbard Hamiltonian is then, H= XX i,j σ X 1 X Tij c†iσ cjσ + U niσ ni,−σ − U νii niσ 2 i,σ (4.11) i,σ where niσ = c†iσ ciσ . A few more simplifications are close at hand, the last term of eq. (4.11) can be written as −U X i,σ νii niσ = −U X X N −1 νk niσ i,σ k = −U X i,σ 1 1 N −1 n niσ = − U N n2 (4.12) 2 2 which is constant and may be ignored. Furthermore, the second term in eq. (4.11) is X 1 X 1 X U niσ ni,−σ = U (ni↑ ni↓ + ni↓ ni↑ ) = U ni↑ ni↓ . 2 2 i,σ i i (4.13) 26 Electron-electron interactions Finally, by considering only nearest neighbour hopping, i.e. hopping from site i to site i + ∆, and denoting Tii = E0 , Ti,i+∆ = t we arrive at a common formulation of the Hubbard Hamiltonian, X X † X H = E0 niσ + t ciσ cjσ + U ni↑ n↓ . (4.14) iσ (i6=j),σ i The Hubbard model has been described as a ‘highly oversimplified model’ [5]. The assumption that e-e interactions are local is clearly not satisfactory for the general case. By Hubbard’s own estimate on 3d transitions metals the e-e interaction in eq. (4.11) are of the order 10-20 eV while the biggest neglected term is 2-3 eV. Nevertheless, the model is by no means simple; exact solutions are only known in one dimension[67] or for two extreme cases • The ‘Band limit’, U =0 in eq. (4.11) • The ‘Atomic limit’, t=0 in eq. (4.11) The model is also sophisticated enough to reproduce a rich variety of phenomena seen in solid state physics such as, metal-insulator transition, antiferromagnetism, ferromagnetism, Tomonaga-Luttinger liquid and superconductivity [111]. 4.2 The variational principle If we are only interested in the ground state of the system, an alternative approach to eq. (4.1) is available through the Rayleigh-Ritz variational method[18]. A trial wave function φ is introduced and an upper bound to the ground state energy E0 is given by E = E[φ] = hφ|Ĥ|φi ≥ E0 . hφ|φi (4.15) Equality holds when φ equals the true ground state wave function (Ψ0 ). In practice, a number of parameters pi may be introduced in the trial wave function and the ground state energy and wave functions can be found by minimising E = E(p1 , p2 , . . . , pm ) over the parameters pi . The accuracy of the method depends on how closely the trial wave function φ resembles the actual wave function – the more parameters introduced, the better the result. However, as pointed out in for example[56], a rough estimate of the number of parameters needed for anything but very small systems is discouraging. If we introduce three parameters per spatial variable and consider a system of N =100 electrons, the total number of parameters (Mp ) over which to perform the minimisation becomes Mp = 33N = 3300 ∼ 10143 ! (4.16) 4.3. Thomas-Fermi model 27 Clearly this is not a feasible problem and is sometimes referred to as the exponential wall, since the size of the problem grows exponentially with the number of electrons. 4.3 Thomas-Fermi model Equation (4.1) may be reformulated using the variational principle[85] δE[Ψ] = 0, (4.17) i.e., solutions Ψ to the Schrödinger equation occurs at the extremum for the functional E[Ψ]. In order to have the wave function normalised, hΨ|Ψi = N , a Lagrange multiplier E is introduced and we arrive at δ hΨ|Ĥ|Ψi − E (hΨ|Ψi − N ) = 0. (4.18) This is, however, only a rephrasing of the original Schrödinger equation (with the constraint hΨ|Ψi = N ) and by no means easier to solve. The idea in the Thomas-Fermi (TF) model is to assume that the ground state density, Z (4.19) n0 (r) = Ψ0 (r, r2 , . . . , rN )Ψ∗0 (r, r2 , . . . , rN )dr2 dr3 . . . drN minimises the energy functional E[n] = T [n] + Uint [n] + Vext [n], (4.20) and thereby replace the wave function Ψ(r1 , r2 , . . . , rN ) with the significantly simpler density n(r). The first term T [n] is the kinetic energy functional, in the original TF-model this was approximated by the expression for a uniform gas of non-interacting electrons. For the 2DEG we may find a similar expression from the 2D-density of states, n2D (ǫ)=m/(~2 π), given in eq. (2.9). The energy over some area L2 is Z ǫF m ǫ2 ∆E = ǫn2D (ǫ)L2 dǫ = 2 F L2 , (4.21) ~ π 2 0 where ǫF is the Fermi energy. At the same time the number of electrons in the area L2 is Z ǫF m N= n2D (ǫ)L2 dǫ = 2 ǫF L2 . (4.22) ~ 0 Thus, 2 1 ~2 π mL2 ~2 π 2 m ǫ2 N = n , (4.23) ∆E = 2 F L2 = 2 2 ~ π 2 2 mL ~ π 2m n being the electron density. The kinetic energy functional for a non-uniform 2DEG is now approximated by Z ~2 π T [n] ≈ TT F [n] = n(r)2 dr (4.24) 2m 28 Electron-electron interactions where the integration is carried out over two dimensions. The second term in eq. (4.20) is the e-e interaction functional and was originally approximated by the classical Hartree energy, Z Z n(r1 )n(r2 ) e2 Uint [n] ≈ UH [n] = dr1 dr2 . (4.25) 2 |r1 − r2 | Finally, the last term is the energy due to interaction with some external potential vext (r) Z Vext [n] = e n(r)vext (r)dr. (4.26) With the constraint Z n(r)dr = N (4.27) included through a Lagrange multiplier µ, which may be identified as the chemical potential, we arrive at Z δ E[n] − µ e n(r)dr − N =0 (4.28) and consequently get the Euler-Lagrange equation Z ~2 π n(r1 ) δE[n] = n(r) + e2 dr1 + evext (r), µ= δn m |r − r1 | (4.29) which is the working equation in the Thomas-Fermi model. 4.4 Hohenberg-Kohn theorems In 1964 Walter Kohn and Pierre Hohenberg proved two theorems[49], essential to any electronic state theory based on the electron density. The first theorem justifies the use of the electron density n(r) as a basic variable as it uniquely defines an external potential and a wave function. The second theorem confirms the use of the energy variational principle for these densities, i.e., for a trial density ñ(r) > 0, with the condition (4.27) fulfilled and E[n] defined in eq. (4.20), it is true that E0 ≤ E[ñ]. (4.30) The derivations below assume non-degeneracy but can be generalised to include degenerate cases[85]. 4.4.1 The first HK-theorem For an N electron system the Hamiltonian (and thereby the wave function Ψ) in eq. (4.1) is completely determined by the external potential vext (r). N is directly obtained from the density through Z N = n(r)dr, (4.31) 4.4. Hohenberg-Kohn theorems 29 while the unique mapping between densities and external potentials are shown through a proof by contradiction. First assume there exists two different ex′ (r)1 , and thereby two different wave functions ternal potentials vext (r) and vext ′ Ψ and Ψ , which yield the same electron density n(r). Using eq. (4.15) we can write hΨ|Ĥ|Ψi = E0 < hΨ′ |Ĥ|Ψ′ i = hΨ′ |Ĥ ′ |Ψ′ i + hΨ′ |Ĥ − Ĥ ′ |Ψ′ i Z ′ (r) dr, (4.32) = E0′ − n(r) vext (r) − vext and at the same time hΨ′ |Ĥ ′ |Ψ′ i = E0′ < hΨ|Ĥ ′ |Ψi = hΨ|Ĥ|Ψi + hΨ|Ĥ ′ − Ĥ|Ψi Z ′ = E0 − n(r) vext (r) − vext (r) dr, (4.33) Adding equation (4.32) and (4.33) gives E0′ + E0 < E0′ + E0 , (4.34) and our assumption that different vext (r) could yield the same density n(r) is incorrect. 4.4.2 The second HK-theorem From the first theorem we have that a density n(r) uniquely determines a wave function Ψ, hence we can we define the functional F [n(r)] = hΨ|T̂ + Ûint |Ψi, (4.35) where T̂ and Ûint signify the kinetic energy and e-e interaction operator. According to the variational principle, the energy functional Z E[Ψ] = hΨ| (T + Uint ) |Ψi + Ψ∗ vext (r)Ψdr (4.36) has a minimum for the ground state Ψ=Ψ0 . For any other Ψ=Ψ̃ Z Z E[Ψ̃] = F [ñ] + vext (r)ñ(r)dr ≥ F [n0 ] + vext (r)n0 (r)dr = E[Ψ0 ] (4.37) | {z } | {z } =E[ñ] =E[n0 ] and we arrive at (4.30). The Hohenberg-Kohn theorems do not help us solve any specific manybody electron problem, however, they do show that there is no principal error 1 ′ differing more than vext (r) − vext (r)=constant 30 Electron-electron interactions in the Thomas-Fermi approach – only an error due to the approximations done for T [n] and Uint [n]. With an exact expression for the functional F [n] = T [n] + Uint [n] we could solve our problem exactly. Furthermore, since there is no reference to the external potential in F [n], knowing this functional would allow us to solve any system. For this reason F [n] is referred to as a universal functional. 4.5 The Kohn-Sham equations Slightly over a year after the HK-theorems were published Walter Kohn and Lu Jeu Sham derived a set of equations, the Kohn-Sham equations, that made density functional calculations feasable[55]. They started by considering a system of non-interacting electrons moving in some effective potential vef f (r), which will be defined later. Because the system is non-interacting the ground state density n(r) can be found by solving the single particle equation ~2 2 ∇ + vef f (r) ϕi (r) = εi ϕi (r), (4.38) − 2m and summing n(r) = X i |ϕi (r)|2 . (4.39) The trick applied by Kohn and Sham was to compare the Euler-Lagrange equation for this non-interacting system with the one in an interacting system. Using the index s on the kinetic energy functional Ts [n] to remind us that it refers to the non-interacting (single-electron) system, the energy functional E[n] equals Z E[n] = Ts [n] + Vef f [n] = Ts [n] + n(r)vef f (r)dr, (4.40) which, with the constraint from eq. (4.27) included by a Lagrange multiplier ε, yields the Euler-Lagrange equation δE[n] δTs [n] δVef f [n] = + −ε δn δn δn δTs [n] = + vef f (r) − ε = 0. δn (4.41) Meanwhile, the energy functional for the interacting system can, with some deliberate rearrangements, be written as E[n] = Ts [n] + UH [n] + Exc [n] + Eext [n], (4.42) where UH [n] is defined in eq. (4.25) and Exc [n] are the corrections needed to make E[n] exact, Exc [n] = T [n] − Ts [n] + Uint [n] − UH [n]. (4.43) 4.6. Local Density Approximation 31 Once again applying the variational principle with the constraint (4.27), gives δE[n] δTs [n] δUH [n] δExc δVext[n] = + + + −ε δn δn δn δn δn Z δTs [n] n(r1 ) + e2 dr1 + vxc (r) + vext (r) − ε = 0. = δn |r − r1 | (4.44) Comparing eq. (4.41) and eq. (4.44), we realise they are identical if we choose Z n(r1 ) vef f (r) = (4.45) dr1 + vxc (r) + vext (r), |r − r1 | which also was the purpose with the reshuffling made in eq. (4.42). The KohnSham procedure can now be summarised in four steps, I) II) initially guess a density n0 (r) compute vef f (r) through eq. (4.45) III) solve eq. (4.38) and compute a new density n0 (r) through (4.39) IV) repeat from step II until n0 (r) is converged However, the problem to find an explicit expression for the term Exc [n] in eq. (4.43) remains. This term should include all the corrections needed to make the energy functional E[n] exact. Part of this correction is due to the Pauli principle2 – the classical e-e energy UH [n] obviously do not take this into account. Part of the correction stems from the fact that the actual wave function can not, in general, be written as some combination of the functions φi in equation (4.38). These two contributions are often written separate, namely as an exchange (Pauli principle) and correlation (exact wave function) contribution, Exc [n] = Ex [n] + Ec [n]. (4.46) One of the more successful ways to approximate them is through the Local Density Approximation (LDA). 4.6 Local Density Approximation In the case of a uniform electron gas the exchange and correlation energy per particle, εx (n) and εc (n), can be determined. In LDA it is then assumed that the total exchange-correlation energy for any system is the sum of these energies weighted with the local density, i.e., Z (4.47) Ex [n] = εx [n]n(r)dr, 2 No two electrons in a given system can be in states characterised by the same set of quantum numbers. 32 Electron-electron interactions Ec [n] = Z εc [n]n(r)dr. (4.48) An expression for εx [n] was first proposed by Dirac as an improvement to the Thomas-Fermi model[27]. For a 2DEG[110], √ 2 2e p n(r), (4.49) εx [n] = − 3 3π 2 ε0 where e is the electron charge and ε0 the permittivity. A closed expression for the correlation part is a little bit more troublesome. Using Monte Carlo methods it may be computed exactly for different densities whereupon these values are interpolated to fit some analytical expression[92, 110, 115]. E.g., with the density parameter rs defined as rs = 1 √ a0 πn and a0 being the Bohr radius, Tanatar and Cerperly proposed[110] √ 1 + C 1 rs C0 e2 εc [n] = − √ 2 ∗ 4πε0 a0 1 + C1 rs + C2 rs + C3 (rs )3/2 (4.50) (4.51) where the coefficients Ci were determined from least square fits with their Monte Carlo simulations. 4.7 Local Spin Density Approximation Density functional theory within LDA can be extended to include electron-spin effects using the local spin density approximation (LSDA). Equation (4.38) is now written as two equations, ~2 2 σ − ∇ + vef (r) ϕσi (r) = εi ϕσi (r), (4.52) f 2m where we differentiate between the two spin species σ =↑, ↓. The spin-up/down densities are given by X ↑ n↑ (r) = |ϕi (r)|2 (4.53) i and n↓ (r) = X i and the total density |ϕ↓i (r)|2 n(r) = n↑ (r) + n↓ (r). The effective potential σ (r) vef f σ 2 vef f (r) = e Z (4.54) (4.55) is n(r1 ) σ dr1 + vxσ (r) + vcσ (r) + vext (r). |r − r1 | (4.56) 4.8. Brief outlook for DFT 33 The exchange energy per particle now depends on both spin up and spin down, εx [n↑ , n↓ ]. This is usually rewritten using the polarisation parameter ζ= n↑ − n↓ n↑ + n↓ (4.57) as εx [n, ζ], which may be expressed in terms of the unpolarised εx [n]. For a 2DEG[110] 1 εx [n, ζ] = εx [n] (1 + ζ)3/2 + (1 − ζ)3/2 . (4.58) 2 The exchange potential vxσ (r) is then given by the functional derivative vxσ (r) = where Ex [n, ζ] = vcσ (r) is similarly obtained from Z vcσ (r) = δEx [n, ζ] , δnσ (4.59) n(r)εx [n, ζ]dr. (4.60) δEc [n, ζ] , δnσ (4.61) where a parametrisation for εc [n, ζ] can be found in[7]. 4.8 Brief outlook for DFT Although the local density approximation has performed remarkably well – it gives systematic errors which often are small (7-10% error in energy)[19, 56] – it is not sufficient for general quantum chemical computations and occasionally fail also in solid state calculations. An example of the latter is the (in)famous case of the ground state in iron not being magnetic within LSDA. Efforts to improve the accuracy of DFT computations focus on finding better exchangecorrelation energy functionals. This work requires ample of imagination and mathematical skill as the only guidance are known asymptotic behaviours and scaling properties of the energy functionals[19, 85]. A natural extension to LSDA is the generalized gradient approximation (GGA)[12, 63, 90, 93] which makes DFT accurate enough for quantum chemical computations[56]. GGA consider the exchange-correlation energy from GGA is of the form both the spin densities and their local gradients, i.e., Exc (c.f. with LSDA in eq. (4.59)-(4.61)) Z GGA (4.62) Exc [n↑ , n↓ ] = d3 rn(r)ǫGGA xc (n↑ , n↓ , ∆n↑ , ∆n↓ ). Some additional extensions to “traditional” DFT are 34 Electron-electron interactions • Time-dependent DFT (TD-DFT) “Traditional” DFT solves the time-independent Schrödinger equation. TD-DFT tries to extend DFT-techniques to systems where the evolution over time is important. • Hybrid functionals A hybrid functional is made up part of the exact exchange from Hartree-Fock theory, part of exchange-correlation from some other theory (e.g., LDA)[13]. Hybrid functionals has been shown to improve the predictive power of DFT for a number of molecular properties such as atomisation energies, bond lengths and vibration frequencies[91] • LDA+U The LDA+U uses the on-site Coulomb interaction, similar to the Hubbard term in eq. (4.14), instead of the averaged Coulomb energy in the energy functional[3, 4, 68]. LDA+U can be more accurate for localised electron systems, strongly correlated materials, insulators, etc. . . Chapter 5 Modelling To model electron transport in mesoscopic systems, a number of issues need to be addressed. E.g., how to discretise analytical equations such that they remain numerically manageable, how, or if, to include electron-electron interactions and how to define boundary conditions without introducing artifacts. There are numerous ways to handle these issues; this chapter only cover the techniques relevant for this thesis. 5.1 Tight-binding Hamiltonian Tight-binding model has been one of the workhorses for problems within semiconductor physics for a long time. To start, we define the notation |m, ni as the direct product, |mi ⊗ |ni=|mi|ni=|mni, representing a state centred at site m, n in the discretisation of the 2D semi-infinite waveguide shown in figure 5.1. The |m, ni fulfills the relationship for completeness and orthonormality X |m, nihm, n| = 1 and hm, n|m′ , n′ i = δmn,m′ n′ . (5.1) m,n An arbitrary state |Ψi can be expanded in this basis as X cmn |m, ni, |Ψi = (5.2) m,n where |cmn |2 is the probability to find the electron at site m, n. The tightbinding Hamiltonian for an electron moving in a perpendicular magnetic field may now be written as Xh Ĥ = |m, ni (ε0 + Vmn ) hm, n| m,n +|m, nithm, n + 1| + |m, nithm, n − 1| i +|m, nite−iqn hm + 1, n| + |m, niteiqn hm − 1, n| . (5.3) 36 Modelling n=N y x n=2 n=1 m=1 m=2 m=3 m=4 Figure 5.1: Discretisation of a two-dimensional semi-infinite quantum waveguide with transversal sites indexed from n=1 . . . N and longitudinal sites m=1 . . . ∞. Vmn is the on-site potential, ε0 the on-site energy, t the hopping integral between sites and the phase factor e±iqn comes from inclusion of a perpendicular magnetic field (Landau gauge) via Peierl’s substitution[33, 34]. With the choice of 1 eBa2 2~2 t = − ε0 , q= , (5.4) ε0 = ∗ 2 , m a 4 ~ m∗ being the effective electron mass, e the electron charge, B the magnetic field and a the lattice discretisation constant, eq. (5.3) converges to its continuous counterpart as a → 0. 5.1.1 Mixed representation Working with quantum wires, it is sometimes convenient to pass from the real space representation with the states |m, ni to a mixed state representation using the transversal lead eigenfunctions (c.f. with the continuous case in eq. 2.15), r 2 πnα ϕn = sin α ∈ N. (5.5) N +1 N +1 Denoting these states |m, αi, where m signifies the longitudinal position and α the transverse mode we can write the transformation between |m, ni and |m, αi as |m, αi = X m′ ,n′ = |m′ , n′ ihm′ , n′ | |m, αi X m′ ,n′ X |m′ , n′ i hm′ , n′ |m, αi = hn′ |αi|m, n′ i n′ = X n′ r 2 sin N +1 πn′ α N +1 |m, n′ i. (5.6) 5.1. Tight-binding Hamiltonian 37 Carrying out the transformation gives (see Ref. [127] for further details on the Hamiltonian) H= + X αβ XX m α |α, miǫα hα, m| αβ |α, mitβα r hβ, m + 1| + |α, mitl hβ, m − 1| + |α, mi(Vαβ + ǫ0 )hβ, m| (5.7) where[16] απ wZ απy 2t w βπy iqy = dye sin sin w 0 w w iqw α+β 1 − e (−1) 1 − eiqw (−1)α−β − = −itqw w2 q 2 − π 2 (α + β)2 w2 q 2 − π 2 (α − β)2 kα = tαβ r αβ ∗ tαβ l = (tr ) Z απy βπy 2 w dy sin sin Vm (y) Vαβ = w 0 w w ǫ0 = −4t απ h2 απ 2 ǫα = 2t cos ≈ + 2t. w 2m∗ w w is the width of the wire and 5.1.2 ∗ (5.8) (5.9) (5.10) (5.11) (5.12) (5.13) (5.14) denotes complex conjugate. Energy dispersion relation We will now consider the energy dispersion relation for the tight-binding Hamiltonian in zero magnetic field. Expanding the solution |Ψi to eq. (5.3) as X |Ψi = ψmn |m, ni (5.15) m,α and substituting this back into eq. (5.3) yields ε0 ψmn + tψm+1,n + tψm−1,n + tψm,n+1 + tψm,n−1 = E. (5.16) With the discussion of the quantum wire in section 2.1.2 in mind we assume a separable solution on the form ψmn = φm ϕn , (5.17) with corresponding longitudinal and transverse energy components, E = Ek + E⊥ = Ekk + Eα . (5.18) 38 Modelling 1 Energy ~2 2m∗ Continuous lattice Discrete lattice 0.5 0 −1 -0.5 0 kk a 0.5 1 Figure 5.2: Continuous (eq. (2.16)) and discrete (eq. (5.23)) energy dispersion relation for the lowest subband. After some rearrangements and writing out ε0 and t from eq. (5.4) we get φm+1 − 2φm + φm−1 ~2 = Ekk φm , (5.19) − ∗ 2m a2 ~2 ϕn+1 − 2ϕn + ϕn−1 − ∗ = Eα ϕn . (5.20) 2m a2 These are discrete derivatives, and we will get travelling wave solutions in the parallel direction and quantised solutions in the transverse direction, φm = eikk ma r 2 πnα sin ϕn = N +1 N +1 kk ∈ R (5.21) α ∈ N. (5.22) Substituting φm back into eq. (5.19) and ϕn back into eq. (5.22) we get the dispersion relations, ~2 2 − 2 cos(kk a) , ∗ 2 2m a ~2 π 2 , α2 Eα = 2m∗ w2 transverse : Ekk = longitudinal : (5.23) (5.24) where w = a(N + 1) is the width of the channel. A comparison of eq. (5.24) with eq. (2.16) identifies απ w as k⊥ . If kk a ≪ 1, eq. (5.23) approximates to the continuous dispersion relation (see figure 5.2), !! ~2 kk2 (kk a)2 ~2 2−2 1− = . (5.25) Ekk ≈ 2m∗ a2 2 2m∗ 5.2. Green’s function 39 Knowing the energy dispersion relation allows us to compute the group velocity, vg , as 1 dE 1 d dω vg = = = dkk ~ dkk ~ dkk 5.2 ~2 2 − 2 cos(kk a) 2m∗ a2 = ~ sin(kk a). m∗ a (5.26) Green’s function The S-matrix technique briefly mentioned in section 3.3.1 gives the response of a system at the boundary due to some excitation at the boundary. This section introduces the more powerful Green’s function technique which gives the response of the system at any point due to an excitation at any point. This will allow us not only to compute the transmission coefficients through the dot but also the electron density inside the system. 5.2.1 Definition of Green’s function Given a differential operator Dop which relates an excitation S in some system to the system response R by Dop R = S, (5.27) the Green’s function is defined as −1 G = Dop . (5.28) Describing an electron travelling in a lead towards a quantum system as a unit excitation I at the system boundary, the Green’s function is defined by HΨ + I =EΨ (E − H)Ψ =I ⇒ (5.29) (5.30) −1 I (5.31) −1 . (5.32) Ψ =(E − H) G = Ψ =(E − H) This definition actually allows two possible Green’s functions, the retarded and the advanced Green’s function. The former describes the wave function caused by the excitation while the latter describes the wave function that generates the excitation. Evidently our interest lies in the retarded Green’s functions and this may technically be specified by adding a small imaginary part to the energy[22], G = (E − H + iδ)−1 . (5.33) 40 Modelling Perturbation | {z H01 | } {z H02 } Figure 5.3: Two systems with respective Hamiltonian H01 and H02 separated by a perturbation. 5.2.2 Dyson equation The most straightforward way to numerically calculate the Green’s function is to perform the inverse in equation (5.32). However, matrix inversion is computationally a very expensive operation and a more manageable approach is to compute the Green’s function recursively through the Dyson equations. Considering two adjacent systems H01 and H02 , figure 5.3, and treating hopping between them as a small perturbation V , we can write the Hamiltonian for the combined system as H = H01 + H02 +V = H0 + V. | {z } (5.34) =H0 From equation (5.32) we have G −1 = E − H = E − H0 − V = G 0 −1 − V, (5.35) where G 0 is the Green’s function of H0 . Multiplying equation (5.35) with G 0 from the right and G from the left (or vice versa) yields the Dyson equation(s) G = G 0 + G 0V G 0 0 (G = G + GV G ). (5.36) (5.37) These equations implicitly define G in terms of the known, unperturbed, G 0 . For explicit expressions for the Green’s function between two specific points r,r′ , we introduce the notation G(r, r′ ) = hr|G|r′ i = hm, n|G|m′ , n′ i = Gmn,m′ n′ . (5.38) Similarly Gmm′ is interpreted as hm|G|m′ i. The matrix Gmm′ can then be written as X 0 0 Gmi Vij Gjm′ . (5.39) Gmm′ = Gmm ′ + i,j For two systems labelled 1 and 2 we have H = H01 + H02 + V = H0 + V12 + V21 (5.40) 5.2. Green’s function 41 e− 0 m (b) (a) Figure 5.4: (a): Typical system under consideration, two leads are connecting the quantum dot to an electron reservoirs at infinity. An incoming electron e− from the right may either be reflected back to the left lead or transmitted to the right. (b): A straight, homogenous channel with zero magnetic field. where V12 = |1it12 h2| = |1iteiqn h2| −iqn V21 = |2it21 h1| = |2ite h1|. (5.41) (5.42) The Green’s functions for the combined function may now be found through eq. (5.39). For example, 0 0 + G11 V12 G21 G11 = G11 G21 = 0 G21 |{z} 0 +G22 V21 G11 , (5.43) (5.44) =0 substituting back and forth yields, after some rearrangements 0 0 0 G11 = (I − G11 V12 G22 V21 )−1 G11 G21 = 0 G22 V21 (I − 0 0 0 G11 V12 G22 V21 )−1 G11 . (5.45) (5.46) In a similar way expressions to add entire sections of recursively computed Green’s functions, such as in figure 5.5 and 5.8, can be derived. 5.2.3 Surface Green’s function Using eq. (5.39) above, we can find the full Green’s function for a finite system by recursively adding parts of the system together. However, usually the system we will consider is infinite, as in figure 5.4a, and the Green’s function in the leads has to be computed in some other way. For a straight, hard wall potential channel with zero magnetic field the solutions to the Shrödinger equation is given by eq. (5.21) and (5.22). An excitation ψ0α |m, αi at slice 0 in figure 5.4b gives the response |Ψi = Gψ0α |0, αi. (5.47) 42 Modelling Because there can be no mixing between different modes α in a homogenous wire we temporarily suppress index α and write |m, αi → |mi. |Ψi is now expanded as X |Ψi = ψm |mi = Gψ0 |0i. (5.48) m Multiplying from the right with hm| gives the response at site m ψm = hm|Gψ0 |0i = Gm0 ψ0 . (5.49) From the Dyson equation (5.37) we find Gm+1,0 as 0 0 Gm+1,0 = Gm+1,0 +Gm+1,m+1 Vm+1,m Gm0 | {z } (5.50) =0 and thereby 0 ψm+1 = Gm+1,0 ψ0 = Gm+1,m+1 Vm+1,m Gm0 ψ0 = eq. (5.49) 0 Gm+1,m+1 Vm+1,m ψm . (5.51) The potential in the wire is periodic and we apply the Bloch theorem, ikkα ma ψm = e um , (5.52) where um is periodic, um = um+1 , and kkα is the wave vector from eq. (2.16) associated with the transverse mode α. From eq. (5.51) we get ikkα (m+1)a e ikkα a e ikkα ma 0 Vm+1,m e = Gm+1,m+1 0 Vm+1,m . = Gm+1,m+1 (5.53) (5.54) Defining the surface Green’s function as Γ̃ = G0m+1,m+1 , where tilde reminds us that this is in mixed representation and noting that for zero magnetic field Vm+1,m = t, we obtain, ik α a e k = tΓ̃. (5.55) With all N modes α1 , α2 , . . . , αN , the surface Green’s function is a diagonal matrix α1 ik a e k 0 ··· α ik 2 a e k 0 1 0 Γ̃ = . (5.56) . . t .. . 0 α ik N a e k which may be transformed to real space representation by Γ(mn, m′ n′ ) = hmn|Γ|m′ n′ i = hn|Γmm′ |n′ i 5.2. Green’s function 43 ′ Gm1 ′ Gm+1,M ′ G1m Γ ′ GM,m+1 ′ G11 ′ GM M Γ 1 m ′ Gmm m+1 M ′ Gm+1,m+1 Gmm Figure 5.5: Recursively adding slices by eq. (5.39) from right and left lead we may compute the Green’s function at slice m, Gmm , for the full system. The diagonal elements of ℑ Gmm is proportional to the density of states (eq. (5.58)). = X α,α′ hn||αihα|Γmm′ |α′ ihα′ ||n′ i = X α,α′ hα|ni∗ Γ̃(mα, m′ α′ )hα′ |n′ i, (5.57) where ’∗’ denotes the complex conjugate. 5.2.4 Computational procedure Figure 5.5 schematically shows the recursive Green’s function technique applied to a two-lead quantum dot. The dot is divided into M slices; the Green’s function for each slice is described by the matrix Gkk while the Green’s function between slice k and l by the matrices Gkl and Glk . Starting from the surface Green’s function Γ, we recursively add slice by slice until we have computed ′ ′ the Green’s functions Gmm and Gm+1,m+1 indicated in figure 5.5. Using the Dyson equation we may now compute the total Green’s function Gmm for slice m. The imaginary part of the diagonal elements in this matrix (ℑ[G(r, r)]) are proportional to the density of states[33], i 1 h n(E, r) = − ℑ G(r, r, E) . π (5.58) Integrating the density of states up to the Fermi-energy we find the electron density inside the dot, Z E f 1 G(r, r, E)dE . (5.59) n(r) = − ℑ π −∞ Using DFT in conjunction with the recursive Green’s function technique we can now write down a procedure to find a self-consistent density to an arbitrary two-lead quantum dot. 44 Modelling x 10 20 ℜ(G) ℑ(G) G(E, r, r) 0 -5 -10 -15 ℑ ℑ E0 -20 EF 0.8 ℜ 0.9 E/EF ℜ E0 1 EF 0.8 0.9 ℜ(E)/EF 1 Figure 5.6: Real (solid line) and imaginary (dashed line) part of the Green’s function at a single point m, n for a typical system. Left panel shows the Green’s function along the real axis (gray arrow in inset, E0 is the lowest potential in the dot, EF the Fermi level). Because of the very sharp peaks numerical integration is difficult along this path. Right panel shows the Green’s function along the path indicated by the gray arrow in the right inset. The Green’s function is analytical in the upper half of the complex plane, hence integration along this path yields the same result as along the real axis but is numerically easier. I) guess an an initial density n(r) (e.g., by Thomas-Fermi) II) compute an effective potential vef f (r) by equation (4.45) III) update the Hamiltonian (eq. (5.3)) with the new vef f (r) IV) update the Green’s function with the new Hamiltonian (eq. (5.32)) V) VI) update the density n(r) through eq. (5.59) repeat from step II until n(r) converges There are some problems with convergence in this procedure. If the density between subsequent iterations varies too much, this may cause the solution to oscillate – this problem is typically addressed by mixing the new density with the old density by some parameter ǫ, nnew (r) = ǫnnew (r) + (1 − ǫ)nold (r) 0 < ǫ ≤ 1, thereby dampening any sudden change in the density. Green’s function has some very sharp features along the left panel of figure 5.6), making numerical integration very poles to the Green’s function are below the real axis this (5.60) Furthermore, the real axis (see the difficult. Since all may be remedied x 10 45 20 2 (a) G × fF D 0 ℜ(G) ℑ(G) -2 fF D (E) ℑ(E)/kB T 5.2. Green’s function 15 10 5 (b) 1 -4 (c) 0 -6 -1 -20 -10 0 (ℜ(x) − EF ) /kB T -20 -10 0 (ℜ(x) − EF ) /kB T Figure 5.7: (a) Real and imaginary parts of Green’s function at point m, n multiplied with the Fermi-Dirac distribution. (b) Integration path in complex plane. (c) Real and imaginary part of Fermi-Dirac distribution (fF D ) along the path in (b). by carrying out the integration in the complex plane (see the right panel of figure 5.6), where the Green’s function is far smoother. This may however still not be enough for convergence. With a high density of states close to the Fermi level the Green’s function may change rapidly as we return towards the real axis, and the last part of the integration requires a very fine stepsize ∆E to capture the integral accurately. A simpler approach is to slightly increase the temperature in the system. With T > 0K some states above the Fermi level EF will be occupied and the density is given by Z ∞ (5.61) n(E, r) fF D (E, EF , T ) dE, ne (r, EF , T ) = {z } −∞ | {z } | DOS F−Ddist. where fF D is the Fermi-Dirac distribution. We may now continue the integration in the complex plane until the Green’s function is suppressed by the exponentially decaying Fermi-Dirac distribution, see figure 5.7. With the self-consistent density nef f (r) we compute the Green’s function relating the first and last slice as shown in figure 5.8, GM 1 . In mixed representation this matrix gives the transmission amplitudes between mode α in the left lead to mode β in the opposite lead. At zero magnetic field[33], r tβ←α = i2|t| i(kkα m−kkβ n) sin kkα a sin kkβ a e G̃M β,1α (E), (5.62) where G̃M β,1α (E) = hM β|G̃(E)|1αi, and tilde indicates that it is the Green’s function in mixed representation. At zero magnetic field the amplitude for the reflected electron is similarly[33] r . ia(k β +k α ) rβ←α = i sin kkβ a sin kkα a e k k 46 Modelling ′ ′ G1,M −1 GM −1,1 Γ ′ GM −1,M −1 ′ G11 ′ GM M 1 G11 M −1 Γ M GM 1 Figure 5.8: With the self-consistent density n(r), we compute the Green’s functions GM1 (G11 ) for the full system, which in mixed representation gives the transmission amplitudes by eq. (5.62) (reflection amplitudes by eq. (5.63)). i h × 2|t| sin kkβ a G̃1β,1α (E) + iδβα , (5.63) where δαβ is the Dirac delta function. Chapter 6 Comments on papers 6.1 Paper I Magnetoconductance through small open quantum dots, i.e., dots with ∼100 electrons or less which are strongly connected to leads, shows a strong correlation between the conductance through the dot and the eigenspectrum of the corresponding closed dot[126]. A computational example is shown in figure 6.1(a) where conductance is plotted against magnetic field and the number of propagating modes in the leads. Superimposed on the conductance in the same figure are the eigenvalues of the corresponding closed dot whose features roughly match the dips and peaks of the conductance. The correlation is not perfect, depending on the lifetime broadening some energy levels are missing or merged in the conductance plot. In paper I we demonstrate that the magnetoconductance of small lateral quantum dots in the strongly coupled regime (i.e. when the leads can support one or more propagating modes) shows a pronounced splitting of the conductance peaks and dips which persists over a wide range of magnetic fields (from zero field to the edge-state regime) and is virtually independent of the magnetic field strength. Some of these measurements are shown in figure 6.2 together with their respective top gate layout as insets. Our numerical analysis of the conductance based on the Hubbard Hamiltonian demonstrates that the splitting is essentially a many-body/spin effect that can be traced to a splitting of degenerate levels in the corresponding closed dot – c.f. figure 6.1(b) for the closed dot eigenspectrum with the conductance plot in figure 6.1(c)-(d). The effect in open dots can be regarded as a counterpart of the Coulomb-blockade effect in weakly coupled dots, with the difference, however, that the splitting of the peaks originates from interactions between electrons of opposite spin. 48 Comments on papers U=0 U=3t 1 1.45 (a) Closed dot 1.5 (b) 0.8 360nm k (units π/w) k (units π/w) 1.4 0.6 1.3 1.25 0.4 1.2 210nm 0.2 kF 1.05 1.1 0 0.2 B (Tesla) 0.4 0 2e2 G( h ) (c) k (units π/w) k (units π/w) Open dot 0.4 80nm U=3t U=3t 1.56 0.2 B (Tesla) 1.52 1.2 (d) 360nm 1.15 1.48 0 0.05 B (Tesla) 0.1 1.1 0.1 0.2 0.3 B (Tesla) 210nm Figure 6.1: Paper I (a) The conductance G = G(kF , B) for U = 0. Solid lines depict the eigenspectrum of the corresponding closed dot. (b) The eigenspectrum of the corresponding closed dot for the same strength of U = 3t. Dotted line indicates kF .(c),(d) The conductance G = G(B, kF ) for U = 3t. Note that because of a large computation time the conductance is shown only for two representative regions. -0.4 (a) (b) -0.62 4 3 VG VG -0.5 2 -0.6 -0.66 1 0.2 0.3 0.4 B (Tesla) 0.5 0.6 0 0.2 0.4 0.6 2e B (Tesla) G( h ) 2 Figure 6.2: Paper I Experimental conductance as a function of magnetic field B and gate voltage VG obtained from different dots. The corresponding device layouts are shown in the insets. The lithographic size of the dots is ∼ 450nm. 6.2 Paper II The existence of a spin polarized state for the two dimensional electron gas has been the extensively probed for various systems, both experimentally, e.g., Refs[8, 41, 44, 112], and theoretically, e.g., Refs[114, 118]. It has been suggested as the explanation for the “0.7”-ananomaly[118] in quantum point contacts and proposed as a working principle for spin filtering through weakly 6.2. Paper II 49 Figure 6.3: Paper II The density for the spin-up and spin-down electrons n↑ (r), n↓ (r) (left and right panels respectively) in a (a) low, (b) intermediate and (c) high density parabolic dot calculated within LSDA. These densities corresponds to an average electron density of ne =∼ 2 × 1014 m−2 , 8 × 1014 m−2 and 2 × 1015 m−2 respectively. T=2K. coupled quantum dots[69, 81]. In open quantum dots, spin polarization has been studied with conflicting results, c.f. Folk et al.[35, 36] and Evaldsson et al.[32] (Paper I). However, the dots studied by Folk et al. were significantly larger, with ∼ 1000 − 10000 electrons and an average electron density ne =∼ 1 × 1015 m−2 , compared with ne =∼ 1 × 1014 m−2 for the dots in Evaldsson et al.. The breaking of spin degeneracy at sufficiently low electron densities is a well studied phenomena[2, 6, 8, 41, 44, 51, 112, 114, 118], and it is possible that the discrepancy between Folk et al. and Evaldsson et al. could be explained by the difference in average electron density in the respective systems. The focus of paper II is coherent transport through open lateral quantum dots using recursive Greens function technique, incorporating exchangecorrelation effects within the Density Functional Theory (DFT) in the local spin-density approximation (LSDA). At low electron densities the current is spin-polarized and the electron density in the dot shows a strong spin polarization. As the electron density increases the spin polarization in the dot gradually diminishes, see figure 6.3 where the left column show the spin-up electron density, n↑ (r), and the right column the spin-down density, n↓ (r). The average density in the plots are n̄(r)/1014 m−2 =, (a) 2, (b) 8, (c) 20. 50 Comments on papers These findings are consistent with available experimental observations. 6.3 Paper III As mentioned in the earlier comments to paper II, a spin-polarised state in low-density 2DEG:s has been studied both experimentally, e.g., Refs. [8, 32, 41, 44, 95, 106, 107, 112] and theoretically, e.g., Refs. [2, 51, 114, 118]. In experimental studies a difficulty arises because spin polarisation is expected at very low electron densities, roughly 1 − 3 × 1014 m−2 [42, 107], At such low Figure 6.4: Paper III Schematic view of the system studied. The heterostructure consists of (from bottom to top), a GaAs substrate, a 3060 nm AlGaAs spacer, a 26 nm donor layer, and a 14 nm cap layer. The side gates define a quantum wire and the top gate controls the electron density in the section with randomly distributed dopants ρd (r). electron densities the electrostatic potential due to impurities in the donor layer can significantly affect the electronic and transport properties of the 2DEG. For example, Nixon et al.[25, 77] showed that a monomode quantum wire could undergo a metallic-insulator transition MIT[9, 25, 30, 52, 77] due to the a random potential from unscreened donors. In paper III we study spin polarization in a split-gate quantum wire focusing on the effect of a realistic smooth potential due to remote donors. A schematic figure is given in figure 6.4. Electron interaction and spin effects are included within the density functional theory in the local spin density approximation. We find that depending on the electron density, the spin polarization exhibits qualitatively different features. For the case of relatively high electron density, when the Fermi energy EF exceeds a characteristic strength of a long-range impurity potential Vdonors , the density spin polarization inside 6.4. Paper IV 51 eV 700 nm EF A) 600 500 −0.02 400 −0.04 300 −0.06 200 −0.08 100 −0.1 0 0 200 400 600 800 1000 1200 1400 1600 1800 2000 400 600 800 1000 nm 1200 1400 1600 1800 2000 700 B) 600 nm 500 400 300 200 100 0 0 200 Figure 6.5: Paper III Potential due to split gates and impurities with top gate potential A) Vg =0, B) Vg =-0.04, for the 30nm thick spacer sample. The first and last segment of ∼400nm show the homogeneous part of the wire. the wire is practically negligible and the wire conductance is spin-degenerate. An example of the potential profile of the wire in this regime is given in figure 6.5(A). When the density is decreased such that EF approaches Vdonors , the electron density and conductance quickly become spin polarized. With further decrease of the density the electrons are trapped inside the lakes (droplets) formed by the impurity potential and the wire conductance approaches the pinch-off regime. An example of the potential profile in this regime is given in figure 6.5(B). We discuss the limitations of the density functional theory in the local spin density approximation in this regime and compare the obtained results with available experimental data. 6.4 Paper IV The word “Spintronics” has been invented to denote the idea to use the spin property of the electron as a building block for future devices[120]. With operational sizes of solid state components approaching the meso- and nanometre sized regimes this idea will soon reach its full potential. A remaining 52 Comments on papers problem within semiconductor spintronics is how to generate and/or inject spin-currents in semi-conductors. Numerous proposals has been put forward; transport through spin blockaded lateral quantum dots[69, 80, 81], the use of spin-orbit coupling effects[100], applying an inhomogeneous in-plane magnetic field[38], magnetic barriers[46, 83, 84], quantum point contacts[99, 101] and spin injection[60, 70, 125]. ↓ ↑, ↓ ↑ Figure 6.6: Paper IV Working principle for antidot spin filtering. A perpendicular magnetic field creates spin-up/down edge states travelling in the direction indicated by arrows. Edge state may tunnel to the quasi-bound state a the antidot and then to the edge state at the opposite edge of the quantum wire. Because of the Zeeman term, ± 21 gµE σB, the maximum tunnelling rate occurs at different magnetic field strengths for spin-up/down electrons. Paper IV proposes a device based on an antidot embedded in a narrow quantum wire in the edge-state regime, figure 6.6, that can be used to inject and/or control spin-polarized current. The operational principle of the device is based on the effect of resonant backscattering from one edge state into another through localized quasibound states, combined with the effect of Zeeman splitting of the quasibound states in a sufficiently high magnetic field. We outline the device geometry, present detailed quantum-mechanical transport calculations, and suggest a possible scheme to test the device performance and functionality. 6.5 Paper V Properties of lateral structures in 2DEGs are often analysed through conductance measurements using a specific device. Additional insights to these measurement are gained by manipulating external parameter such as the bias voltage, applying a gate voltage to a top gate or a magnetic field. Paper V analyse the details of low energy electron transport through a magnetic barrier (MB), i.e., a static inhomogeneous magnetic field applied locally to the 2DEG. Typically such a device is prepared by placing a ferromagnetic film on top of the heterostrucutre[53, 59, 73], figure 6.7. The magnetisation of the film, µ0 M , induces a perpendicular magnetic field Bz (x) inhomogeneous in the xdirection and homogeneous in the y-direction, see top inset in figure 6.7. The maximum field strength of Bz (x) and its full width at half maximum (FWHM) 6.5. Paper V 53 Bz (x) x side view z x µ0 M top view y x µ0 M spacer magn. film spacer 2DEG substrate magn. film µ0 M x=0 Figure 6.7: Paper V Schematic view of a magnetic barrier (MB) structure. A magnetic film is placed on top of a 2DEG heterostructure. Depending on the magnetization µ0 M , the film thickness and the spacer thickness the shape of the magnetic field Bz (x) at the 2DEG can be manipulated. Figure 6.8: Paper V Calculated conductance versus Fermi energy EF for three samples with 15nm, 35nm and 350nm spacer thickness and µ0 M = 1.2T. The lower inset shows the shapes of the MBs present at the distances considered. Only the most localised MBs, with 15nm and 35nm spacer, displays the resonance dip. can be adjusted by changing the magnetisation strength µ0 M , the width of the spacer and the width of the magnetic film, see inset to figure 6.8. With increasing magnetic field more propagating modes in the wire are reflected 54 Comments on papers and quantized conductance steps similar to those appearing at quantum point contacts (QPC) appear[45, 122, 123]. However, in contrast to QPCs there is an additional transmission dip for some MB just before the magnetic pinch off, see conductance curves for 15nm and 35nm samples in figure 6.8. Paper V explains this resonant feature as an interference effect between an open edge state in the MB-region and a localised vortex state in the middle of the MB. 6.6 Paper VI Paper VI adapts a numerical method for efficiently computing the surface Green’s function in photonic crystals[97] to graphene nanoribbons (GNR). Present implementations of the Green’s function technique in GNRs, e.g., Refs [43, 75], use a self-consistent, typically time-consuming, procedure to compute the surface Green’s function. The method introduced in paper VI computes the surface Green’s function directly from the Bloch states by solving a 2N ×2N eigenequation where N is the number of sites across the nanoribbon. The paper further reviews steps to compute wave functions, group velocities and conductance in two-terminal GNRs within the Green’s function formalism. Next, the method is applied to study the magnetoconductance through a magnetic barrier (MB), i.e., a static inhomogeneous magnetic field applied locally to a GNR. MBs are interesting as means to control electron conductance in 2D-systems[71, 87], and for their predicted spin polarising properties[46, 83, 84]. An example of a MB-setup, for a quantum wire, is given in figure 6.7. In our study we note that MBs in GNRs locally modify the energy spectrum and thereby the conductance steps. MBs with steep gradient show pronounced Fabri-Pérot oscillations whereas smoother MB hardly affect the magnetoconductance at all. For disorderfree zigzag-GNRs, where the lowest propagating mode is a dispersionless edge state, we find that a MB could not be used as a conductance “on-off switch” as in the case with MB in quantum wires (Paper V). Analysing the wave function for this mode show that the wave function remains at the edge of the ribbon as the strength of the MB is varied. In contrast to simulated (disorderfree) GNRs, experimental measurements on ribbons of similar widths as studied here always posses an energy gap independent of the crystal direction of the GNR1 [47]. Theoretical work has proposed the the presence of edge disorder as an explanation[96] (see also paper VII) to this phenomena. Because of this we also investigate the effect of MB in GNRs with edge disorder. For this case we find it is possible that a MB could delocalise states localised by edge disorder. 1 I.e., their edge structure. 6.7. Paper VII 6.7 55 Paper VII The high mobility of graphene, even at ambient temperature and high degree of doping, has mad graphene nano ribbons (GNR) a strong candidate for building blocks in future electronic devices. There are, however, some issues that need to be addressed regarding graphene electronics. First, graphene is gapless, making it necessary to engineer a band gap if it is to replace existing semiconductor devices. Second, modelling indicates that the electronic properties of GNR are expected to depend strongly on the edge structure of the ribbon[39, 76, 116]. The first issue can be addressed by making sufficiently narrow GNRs as this will induce a band gap. The second issue might be a result of modelling of ideal GNRs. Measurements by Han et al.[47] found no directional dependence of the conductance for GNRs of widths 20-100nm. Since these measurements, a number of explanations have been put forward to address this discrepancy between experiments and modelling; rough edges[20, 66, 96], atomic scale impurities[20], coulomb blockaded transport[108] and impurity scattering[64]. In paper VII we study the effect of the edge disorder on the conductance of the graphene nanoribbons (GNRs) within the tight-binding model. The method we use, developed in paper VI, allows us to probe GNRs of similar size as in the recent experiments by Han et al.. We find that only very modest edge disorder is sufficient to induce the conduction energy gap in the otherwise metallic GNRs and to lift any difference in the conductance between nanoribbons of different edge geometry. We relate the formation of the conduction gap to the pronounced edge disorder induced Anderson-type localization which leads to the strongly enhanced density of states at the edges, formation of surface-like states and to blocking of conductive paths through the ribbons. 56 Comments on papers Bibliography [1] Tsuneya Ando. Theory of electronic states and transport in carbon nanotubes. Journal of the Physical Society of Japan, 74:777, 2005. [2] A. V. Andreev and A. Kamenev. Itinerant ferromagnetism in disordered metals: A mean-field theory. Phys. Rev. Lett., 81:3199, 1998. [3] Vladimir I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein. Firstprinciples calculations of the electronic structure and spectra of strongly correlated systems: the LDA+U method. Journal of Physics: Condensed Matter, 9:767, 1997. [4] Vladimir I. Anisimov, Jan Zaanen, and Ole K. Andersen. Band theory and Mott insulators: Hubbard U instead of Stoner I. Phys. Rev. B, 44:943, 1991. [5] Neil W. Ashcroft and N. David Mermin. Solid State Physics. Brooks Cole, first edition, 1976. [6] A. Ashok, R. Akis, D. Valsileska, and D. K. Ferry. Spin polarization in GaAs/Al0.24 Ga0.76 As heterostructures. Molecular Simulation, page 797, 2005. [7] Claudio Attaccalite, Saverio Moroni, Paola Gori-Giorgi, and Giovanni B. Bachelet. Correlation energy and spin polarization in the 2D electron gas. Phys. Rev. Lett., 88:256601, 2002. [8] O. M. Auslaender, H. Steinberg, A. Yacoby, Y. Tserkovnyak, B. I. Halperin, K. W. Baldwin, L. N. Pfeiffer, and K. W. West. Spin-charge separation and localization in one dimension. Science, 308:88, 2005. [9] M. Baenninger, A. Ghosh, M. Pepper, H. E. Beere, I. Farrer, P. Atkinson, and D. A. Ritchie. Local transport in a disorder-stabilized correlated insulating phase. Phys. Rev. B, 72:241311(R), 2005. [10] M. N. Baibich, J. M. Broto, A. Fert, F. Nguyen Van Dau, F. Petroff, P. Eitenne, G. Creuzet, A. Friederich, and J. Chazelas. Giant magnetoresistance of (001)Fe/(001)Cr magnetic superlattices. Phys. Rev. Lett., 61:2472, 1988. 58 Bibliography [11] Verónica Barone, Oded Hod, and Gustavo E. Scuseria. Electronic structure and stability of semiconducting graphene nanoribbons. Nano Letters, 6:2748, 2006. [12] A. D. Becke. Density-functional exchange-energy approximation with correct asymptotic behavior. Physical Review A, 38:3098, 1988. [13] Axel D. Becke. A new mixing of hartreefock and local density-functional theories. Journal of Chemical Physics, 98:1372, 1993. [14] R. S. Becker, J. A. Golovchenko, and B. S. Swartzentruber. Atomic-scale surface modifications using a tunnelling microscope. Nature, 325:419, 1987. [15] G. Binasch, P. Grünberg, F. Saurenbach, and W. Zinn. Enhanced magnetoresistance in layered magnetic structures with antiferromagnetic interlayer exchange. Phys. Rev. B, 39:4828, 1989. [16] Torbjörn Blomquist. Waves and Rays in Open Billiards. Ph.d. thesis, Linköping University, 2002. [17] K. I. Bolotin, K. J. Sikes, Z. Jiang, M. Klima, G. Fudenberg, J. Hone, P. Kim, and H. L. Stormer. Ultrahigh electron mobility in suspended graphene. Solid State Communications, 146:351, 2008. [18] B.H. Bransden and C.J. Joachain. Quantum Mechanics. Prentice Hall, second edition, 2000. [19] Kieron Burke. The abc of http://chem.ps.uci.edu/ kieron/dft/book/gamma/g1.pdf, 2007. dft. [20] Zhihong Chen, Yu-Ming Lina, Michael J. Rooks, and Phaedon Avouris. Graphene nano-ribbon electronics. Physica E, 40:228, 2007. [21] M. Ciorga, A. S. Sachrajda, P. Hawrylak, C. Gould, P. Zawadzki, S. Jullian, Y. Feng, and Z. Wasilewski. Addition spectrum of a lateral dot from coulomb and spin-blockade spectroscopy. Phys. Rev. B, 61:R16315, 2000. [22] Supriyo Datta. Electronic Transport in Mesoscopic Systems. The press syndicate of the university of Cambridge, 1997. [23] Supriyo Datta and Biswajit A. Das. Electronic analog of the electro-optic modulator. Appl. Phys. Lett., 56:665, 1990. [24] John H. Davies. The physics of low-dimensional semiconductors. The press syndicate of the university of Cambridge, 1998. [25] John H. Davies and John A. Nixon. Fluctuations in submicrometer semiconducting devices caused by the random positions of dopants. Phys. Rev. B, 39:3423, 1989. Bibliography 59 [26] P.A.M. Dirac. Quantum mechanics of many-electron systems. Proc. Royal Soc. (A), 123:714, 1929. [27] P.A.M. Dirac. Note on exchange phenomena in the thomas atom. Proc. Cambridge Phil. Soc., 26:376, 1930. [28] A. Dorn, A. Fuhrer, T. Ihn, T. Heinzel, K. Ensslin, W. Wegscheider, and M. Bichler. Electronic properties of antidot lattices fabricated by atomic force lithography. Appl. Phys. Lett., 80:252, 2002. [29] V. A. Miransky E. V. Gorbar, V. P. Gusynin and I. A. Shovkovy. Magnetic field driven metal-insulator phase transition in planar systems. Phys. Rev. B, 66:045108, 2002. [30] A. L. Efros, F. G. Pikus, and V. G. Burnett. Density of states of a two-dimensional electron gas in a long-range random potential. Phys. Rev. B, 47:2233, 1993. [31] Claire Berger et al. Electronic confinement and coherence in patterned epitaxial graphene. Science, 312:5777, 2006. [32] M. Evaldsson, I. V. Zozoulenko, M. Ciorga, P. Zawadzki, and A. S. Sachrajda. Spin splitting in open quantum dots. Europhysics Letters, 68:261, 2004. [33] David K. Ferry and Stephen M. Goodnick. Transport in nanostructures. The press syndicate of the university of Cambridge, 1997. [34] Richard P. Feynman, Robert B. Leighton, and Matthew Sands. The Feynman lectures on physics. Addison-Wesly Publishing Company, third edition, 1965. [35] J. A. Folk, S. R. Patel, K. M. Birnbaum, C. M. Marcus, C. I. Duruöz, and J. S. Harris, Jr. Spin degeneracy and conductance fluctuations in open quantum dots. Phys. Rev. Lett., 86:2102, 2001. [36] J. A. Folk, R. M. Potok, C. M. Marcus, and V. Umansky. A gatecontrolled bidirectional spin filter using quantum coherence. Science, 299:679, 2003. [37] E. Fradkin. Critical behavior of disordered degenerate semiconductors. Phys. Rev. B, 33:3263, 1986. [38] D. Frustaglia, M. Hentschel, and K. Richter. Quantum transport in nonuniform magnetic fields: Aharonov-bohm ring as a spin switch. Phys. Rev. Lett., 87:256602, 2001. 60 Bibliography [39] Mitsutaka Fujita, Katsunori Wakabayashi, Kyoko Nakada, and Koichi Kusakabe. Peculiar localized state at zigzag grapite edge. Journal of the Physical Society of Japan, 65:1920, 1996. [40] A. K. Geim and K. S. Novoselov. The rise of graphene. Nature materials, 6:183, 2007. [41] A. Ghosh, C. J. B. Ford, M. Pepper, H. E. Beere, and D. A. Ritchie. Possible evidence of a spontaneous spin polarization in mesoscopic twodimensional electron systems. Phys. Rev. Lett., 92:116601, 2004. [42] P. Giudici, A. R. Goñi, P. G. Bolcatto, C. R. Proetto, C. Thomsen, K. Eberl, and M. Hauser. Evidence of breakdown of the spin symmetry in diluted 2d electron gases. Europhysics Letters, 77:37003, 2007. [43] Roksana Golizadeh-Mojarad, A. N. M. Zainuddin, Gerhard Klimeck, and Supriyo Datta. Atomistic non-equilibrium Green s function simulations of graphene nano-ribbons in the quantum Hall regime. arXiv:condmat/0801.1159, 2008. [44] A. R. Goñi, P. Giudici, F. A. Reboredo, C. R. Proetto, C. Thomsen, K. Eberl, and M. Hauser. Evidence of spontaneous spin polarization in the two-dimensional electron gas. Phys. Rev. B, 70:195331, 2004. [45] Michele Governale and Daniel Boese. Magnetic barrier in confined twodimensional electron gases: Nanomagnetometers and magnetic switches. Appl. Phys. Lett., 77:3215, 2000. [46] Yong Guo, Bing-Lin Gu, Jing-Zhi Yu Zhong Zeng, and Yoshiyuki Kawazoe. Electron-spin polarization in magnetically modulated quantum structures. Phys. Rev. B, 62:2635, 2000. [47] Melinda Y. Han, Barbaros Özyilmaz, Yuanbo Zhang, and Philip Kim. Energy band-gap engineering of graphene nanoribbons. Phys. Rev. Lett., 98:206805, 2007. [48] Jonas Rahlf Hauptmann. Spin-transport in carbon nanotubes. Master’s thesis, University of Copenhagen, 2003. [49] P. Hohenberg and W. Kohn. Inhomogeneous electron gas. Physical Review, 136:B864, 1964. [50] J. Hubbard. Electron correlations in narrow energy bands. Proc. Roy. Soc., 276:238, 1963. [51] S. Ihnatsenka and I. V. Zozoulenko. Conductance of a quantum point contact based on spin-density-functional theory. Phys. Rev. B, 76:045338, 2007. Bibliography 61 [52] S. Ilani, A. Yacoby, D. Mahalu, and Hadas Shtrikman. Unexpected behavior of the local compressibility near the B = 0 metal-insulator transition. Phys. Rev. Lett., 84:3133, 2000. [53] Mark Johnson, Brian R. Bennett, M. J. Yang, M. M. Miller, and B. V. Shanabrook. Hybrid Hall effect device. Appl. Phys. Lett., 71:974, 1997. [54] M. I. Katsnelson, K. S. Novoselov, and A. K. Geim. Chiral tunnelling and the Klein paradox in graphene. Nature Physics, 2:620, 2006. [55] W. Kohn and L.J. Sham. Self-consistent equations including exchange and correlation effects. Physical Review, 140:A1133, 1965. [56] Walter Kohn. Nobel lecture: Electronic structure of matterwave functions and density functionals. Review of Modern Physics, 71, 1999. [57] Leo Kouwenhoven and Leonid Glazman. Revival of the Kondo effect. Physics World, page 33, 2001. [58] H. W. Kroto, J. R. Heath, S. C. O’Brien, R. F. Curl, and R. E. Smalley. C60: Buckminsterfullerene. Nature, 318:162, 1985. [59] V. Kubrak, F. Rahman, B. L. Gallagher, P. C. Main, , M. Henini, C. H. Marrows, and M. A. Howson. Magnetoresistance of a two-dimensional electron gas due to a single magnetic barrier and its use for nanomagnetometry. Appl. Phys. Lett., 74:2507, 1999. [60] V. P. LaBella, D. W. Bullock, Z. Ding, C. Emery, A. Venkatesan, W. F. Oliver, G. J. Salamo, P. M. Thibado, and M. Mortazavi. Spatially resolved spin-injection probability for gallium arsenide. Science, 292:1518, 2001. [61] L. D. Landau. Theory of phase transformations. II. Phys. Z. Sowjetunion, 11:26, 1937. [62] L. D. Landau and E. M. Lifshitz. Statistical Physics, Part I. Pergamon, Oxford, 1980. [63] David C. Langreth and M. J. Mehl. Beyond the local-density approximation in calculations of ground-state electronic properties. Phys. Rev. B, 28:1809, 1983. [64] Aurélien Lherbier, Blanca Biel, Yann-Michel Niquet, and Stephan Roche. Transport length scales in disordered graphene-based materials: Strong localization regimes and dimensionality effects. Phys. Rev. Lett., 100:036803, 2008. [65] T. C. Li and Shao-Ping Lu. Quantum conductance of graphene nanoribbons with edge defects. Phys. Rev. B, 77:085408, 2008. 62 Bibliography [66] Xiaolin Li, Xinran Wang, Li Zhang, Sangwon Lee, and Hongjie Dai. Chemically derived, ultrasmooth graphene nanoribbon semiconductors. Science, 319:1229, 2008. [67] Elliott H. Lieb and F. Y. Wu. Absence of Mott transition in an exact solution of the short-range, one-band model in one dimension. Phys. Rev. Lett., 20:1445, 1968. [68] A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen. Density-functional theory and strong interactions: Orbital ordering in Mott-Hubbard insulators. Phys. Rev. B, 52:R5467, 1995. [69] Ciorga M, Wensauer A, Pioro-Ladriere M, Korkusinski M, Kyriakidis J, Sachrajda AS, and Hawrylak P. Collapse of the spin-singlet phase in quantum dots. Phys. Rev. Lett., 88:256804, 2002. [70] I. Malajovich, J. J. Berry, N. Samarth, and D. D. Awschalom. Persistent sourcing of coherent spins for multifunctional semiconductor spintronics. Nature, 411:770, 2001. [71] A. Matulis, F. M. Peeters, and P. Vasilopoulos. Persistent sourcing of coherent spins for multifunctional semiconductor spintronics. Phys. Rev. Lett., 72:1518, 1994. [72] Jannik C. Meyer, A. K. Geim, M. I. Katsnelson, K. S. Novoselov, T. J. Booth, and S. Roth. The structure of suspended graphene sheets. Nature letters, 446:60, 2007. [73] F. G. Monzon, Mark Johnson, and M. L. Roukes. Strong Hall voltage modulation in hybrid ferromagnet/semiconductor microstructures. Appl. Phys. Lett., 71:3081, 1997. [74] G.E. Moore. Cramming more components onto integrated circuits. Electronics, 38, 1965. Also available at <http://www.intel.com/research/silicon/moorespaper.pdf>. [75] F. Muñoz-Rojas, D. Jacob, J. Fernández-Rossier, and J. J. Palacios. Coherent transport in graphene nanoconstrictions. Phys. Rev. B, 74:195417, 2006. [76] Kyoko Nakada, Mitsutaka Fujita, Gene Dresselhaus, and Mildred S. Dresselhaus. Edge state in graphene ribbons: Nanometer size effect and edge shape dependence. Phys. Rev. B, 54:17954, 1996. [77] John A. Nixon and John H. Davies. Potential fluctuations in heterostructure devices. Phys. Rev. B, 41:7929, 1990. Bibliography 63 [78] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov. Two-dimensional gas of massless dirac fermions in graphene. Nature, 438:197, 2005. [79] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov. Electric field effect in atomically thin carbon films. Science, 306:666, 2004. [80] Hawrylak P, Gould C, Sachrajda A, Feng Y, and Wasilewski Z. Collapse of the zeeman gap in quantum dots due to electronic correlations. Phys. Rev. B, 59:2801, 1999. [81] Recher P, Sukhorukov EV, and Loss D. Quantum dot as spin filter and spin memory. Phys. Rev. Lett., 85:1962, 2000. [82] D.W. Palmer. The semiconductors-information web-site. www.semiconductors.co.uk, 2003. [83] G. Papp and F. M. Peeters. Erratum: ”spin filtering in a magneticelectric barrier structure”. Appl. Phys. Lett., 79:3198, 2001. [84] G. Papp and F. M. Peeters. Spin filtering in a magneticelectric barrier structure. Appl. Phys. Lett., 78:2184, 2001. [85] Robert G. Parr and Weitao Yang. Density-Functional Theory of Atoms and Molecules. Oxford University Press, 1989. [86] B. Partoens and F. M. Peeters. From graphene to graphite: Electronic structure around the k point. Phys. Rev. B, 74:075404, 2006. [87] F. M. Peeters and A. Matulis. Quantum structures created by nonhomogeneous magnetic fields. Phys. Rev. B, 48:15166, 193. [88] R. E. Peierls. Bemerkungen über Umwandliungstemperaturen. Helv. Phys. Acta, 7:81, 1934. [89] R. E. Peierls. Quelques proprietes typiques des corpses solides. Ann. Inst. Henri Poincare, 5:177, 1935. [90] John P. Perdew, Kieron Burke, and Matthias Ernzerhof. Generalized gradient approximation made simple. Phys. Rev. Lett., 77:3865, 1996. [91] John P. Perdew and Matthias Ernzerhof. Rationale for mixing exact exchange with density functional approximations. Journal of Chemical Physics, 105:9982, 1996. [92] John P. Perdew and Yue Wang. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B, 45:13244, 1992. 64 Bibliography [93] John P. Perdew and Wang Yue. Accurate and simple density functional for the electronic exchange energy: Generalized gradient approximation. Phys. Rev. B, 33:8800, 1986. [94] Gary A. Prinz. Magetoelectronics. Science, 282:1660, 1998. [95] O. Prus, Y. Yaish, M. Reznikov, U. Sivan, and V. Pudalov. Thermodynamic spin magnetization of strongly correlated two-dimensional electrons in a silicon inversion layer. Phys. Rev. B, 67:205407, 2003. [96] D. Querlioz, Y. Apertet, A. Valentin, K. Huet, A. Bournel, S. GaldinRetailleau, and P. Dollfus. Suppression of the orientation effects on bandgap in graphene nanoribbons in the presence of edge disorder. Applied Phyics Letters, 92:042108, 2008. [97] A. I. Rahachou and I. V. Zozoulenko. Light propagation in finite and infinite photonic crystals: The recursive Green’s function technique. Phys. Rev. B, 72:155117, 2005. [98] S. Reich, J. Maultzsch, C. Thomsen, and P. Ordejón. Tight-binding description of graphene. Phys. Rev. B, 66:035412, 2002. [99] Potok RM, Folk JA, Marcus CM, and Umansky V. Detecting spin-polarized currents in ballistic nanostructures. Phys. Rev. Lett., 89:266602, 2002. [100] Manuel Valı́n Rodrı́guez, Antonio Puente, Llorenç Serra, and Enrico Lipparini. Spin switching in semiconductor quantum dots through spinorbit coupling. Phys. Rev. B, 66:165302, 2002. [101] L. P. Rokhinson, L. N. Pfeiffer, and K.W. West. Spontaneous spin polarization in quantum point contacts. Phys. Rev. Lett., 96:156602, 2006. [102] R. Saito, G. Dresselhaus, and M. S. Dresselhaus. Physical Properties of Carbon Nanotubes. Imperial, London, 1998. [103] Adriaan M. J. Schakel. Relativistic quantum Hall effect. Phys. Rev. D, 43:1428, 1991. [104] F. Schedin, A. K. Geim, S. V. Morozov, E. W. Hill, P. Blake, M. I. Katsnelson, and K. S. Novoselov. Detection of individual gas molecules adsorbed on graphene. Nature materials, 6:652, 2007. [105] Gordon W. Semenoff. Condensed-matter simulation of a three-dimensional anomaly. Phys. Rev. Lett., 54:2449, 1984. Bibliography 65 [106] A. A. Shashkin, S. Anissimova, M. R. Sakr, S. V. Kravchenko, V. T. Dolgopolov, and T. M. Klapwijk. Pauli spin susceptibility of a strongly correlated two-dimensional electron liquid. Phys. Rev. Lett., 96:036403, 2006. [107] Christoph Siegert, Arindam Ghosh, Michael Pepper, Ian Farrer, and David A. Ritchie. The possibility of an intrinsic spin lattice in highmobility semiconductor heterostructures. Nature Physics, 3:315, 2007. [108] F. Sols, F. Guinea, and A. H. Castro Neto. Coulomb blockade in graphene nanoribbons. Phys. Rev. Lett., 99:166803, 2007. [109] Young-Woo Son, Marvin L. Cohen, and Steven G. Louie. Energy gaps in graphene nanoribbons. Phys. Rev. Lett., 97:216803, 2006. [110] B. Tanatar and D. M. Ceperley. Ground state of the two-dimensional electron gas. Phys. Rev. B, 39:5005, 1989. [111] Hal Tasaki. The hubbard model - an introduction and selected rigorous results. J. Phys.: Condens. Matter, 10:4353, 1998. [112] K. J. Thomas, J. T. Nicholls, M. Y. Simmons, M. Pepper, D. R. Mace, and D. A. Ritchie. Transmission of interacting electrons throuh a onedimensional periodic potential. Phys. Rev. Lett., 77:135, 1996. [113] J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C. W. J. Beenakker. Sub-poissonian shot noise in graphene. Phys. Rev. Lett., 96:246802, 2006. [114] D. Varsano, S. Moroni, and G. Senatore. Spontaneous magnetization of the 2D electron gas. Europhys. Lett., 53:348, 2001. [115] U. von Barth and L. Hedin. A local exchange correlation potential for the spin poarized case:i. J. Phys. C: Solid State Phys., 5:1629, 1972. [116] Katsunori Wakabayashi, Mitsutaka Fujita, Hiroshi Ajiki, and Manfred Sigrist. Electronic and magnetic properties of nanographite ribbons. Phys. Rev. B, 59:8271, 1999. [117] P. R. Wallace. The band theory of graphite. Physical Review, 71:622, 1947. [118] Chuan-Kui Wang and K.-F. Berggren. Spin slitting of subbands in quasione-dimensional electron quantum channels. Phys. Rev. B, 54:14257, 1996. [119] Dawen Wang, Liming Tsau, , and K. L. Wang. Nanofabrication of thin chromium film deposited on si(100) surfaces by tip induced anodization in atomic force microscopy. Appl. Phys. Lett., 67:1295, 1995. 66 Bibliography [120] S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von Molnár, M. L. Roukes, A. Y. Chtchelkanova, and D. M. Treger. Spintronics: A spin-based electronics vision for the future. Science, 294:1488, 2001. [121] Y. Zhang Y.-W. Tan, K. Bolotin, Y. Zhao, S. Adam, E. H. Hwang, S. Das Sarma, H. L. Stormer, and P. Kim. Measurement of scattering rate and minimum conductivity in graphene. Phys. Rev. Lett., 99:246803, 2007. [122] Feng Zhai and H. Q. Xu. Generation of spin polarization in two-terminal electron waveguides by spin-orbit interaction and magnetic field modulations. Phys. Rev. B, 72:085314, 2005. [123] Feng Zhai and H. Q. Xu. Spin filtering in single magnetic barrier structures revisited. Appl. Phys. Lett., 88:032502, 2006. [124] Yuanbo Zhang, Yan-Wen Tan, Horst L. Stormer, and Philip Kim. Experimental observation of the quantum Hall effect and Berry’s phase in graphene. Nature, 438:201, 2005. [125] H. J. Zhu, M. Ramsteiner, H. Kostial, M. Wassermeier, H.-P. Schönherr, and K. H. Ploog. Room-temperature spin injection from Fe into GaAs. Phys. Rev. Lett., 87:016601, 2001. [126] I. V. Zozoulenko, A. S. Sachrajda, C. Gould, K.-F. Berggren, P. Zawadzki, Y. Feng, and Z. Wasilewski. Few-electron open dots: Single level transport. Phys. Rev. Lett., page 1838, 1999. [127] I.V. Zozoulenko, Frand A. Maaø, and E.H. Hauge. Coherent magnetotransport in confined arrays of antidots. I. disperion relations and current densities. Phys. Rev. B, 53:7975, 1996.