A BGK model for charge transport in graphene. Armando Majorana Department of Mathematics and Computer Science, University of Catania, Italy
|
|
- Deirdre McKenzie
- 5 years ago
- Views:
Transcription
1 A BGK model for charge transport in graphene Armando Majorana Department of Mathematics and Computer Science, University of Catania, Italy December 6, 28 arxiv:82.862v math-ph 5 Dec 28 Abstract The Boltzmann euation describes the detailed microscopic behaviour of a dilute gas, and represents the basis of the kinetic theory of gases. In order to reduce the difficulties in solving the Boltzmann euation, simple expressions of a collision operator have been proposed to replace the true Boltzmann integral term. These new euations are called kinetic models. The most popular and widely used kinetic model is the Bhatnagar-Gross-Krook BGK model. In this work we propose and analyse a BGK model for charge transport in graphene. MSC-class: 82C4-82C7 Primary 82D37 Secondary Introduction Graphene is a gapless semiconductor made of a sheet composed of a single layer of carbon atoms arranged into a honeycomb hexagonal lattice. In view of application in graphene-based electron devices, it is crucial to understand the basic transport properties of this material. In a semiclassical kinetic setting, the charge transport in graphene is described by four Boltzmann euations, one for electrons in the valence π band and one for electrons in the conductions π band, that in turn can belong to the K or K valley. In this paper we study the case of a single distribution function f for electrons belonging to a conduction band. This corresponds to a physical case, where a n-type doping or euivalently a high value of the Fermi potential is considered, and the electrons, belonging to a conduction band, do not move to the valence band. Moreover K and K are considered euivalent. Under these assumptions Boltzmann euation writes f t + kε x f e E kf = df, dt coll where t is the time, x and k are the position and the wave-vector of a charge particle, respectively. We denote by x and k the gradients with respect to the position and wave-vector, respectively. With a very good approximation the energy bands is given by εk = v F k, where v F is the constant Fermi velocity, and is the Planck constant divided by 2π. The elementary positive charge is denoted by e, and E is the electric field. The collision operator of Euation describes the interaction of electrons with acoustic, optical and K phonons. If the electric field E is not constant, then it must be self-consistently evaluated by coupling the Boltzmann euation with the Poisson euation for the electrostatic potential.
2 More recently accurate numerical solutions, based on the discontinuous Galerkin method, to the Boltzmann euation have been shown 2, 3. A comparison with results, obtained by using a direct simulation Monte Carlo approach, shows an excellent agreement, which gives a further validation of the numerical scheme. In this work the collision operator of Euation is replaced by a relaxation BGK collision operator. We propose an operator, which retains the fundamental properties of the Boltzmann euation, such the mass conservation, the same euilibrium distribution functions and it properly deals with Pauli s exclusion principle in the degenerate case. 2 The BGK model The collision operator of Euation writes 4 Sk,kft,x,k ft,x,kdk Sk,k ft,x,k ft,x,k dk, 2 where the transition rate Sk,k, related to electron-phonon scatterings, is described in detail in the following. The collision operator 2 vanishes if f is the Fermi-Dirac distribution Fk,µ =, εk µ +exp where k B is the Boltzmann constant, T the constant graphene lattice temperature, and µ is the chemical potential. If the electric field E is null, then Fermi-Dirac distributions are solutions, which do not depend on time t and space coordinates x, of the Boltzmann euation. In this paper we propose the following BGK collision operator Qft,x,k = Sk,kFk,µ ft,x,kdk Sk,k ft,x,k Fk,µ dk. 3 It is derived from Euation 2, replacing the distribution ft,x,k, inside the integrals, with a Fermi-Dirac distribution. Now the chemical potential µ becomes a new unknown. Therefore we must add an euation to the kinetic model. It is Qft,x,kdk =. 4 This euation guarantees the mass conservation, exactly as the same for the Boltzmann euation. If we define Φ k,µ = Sk,kFk,µdk and Φ k,µ = Sk,k Fk,µ dk, 5 then 3 becomes Qft,x,k = Φ k,µ ft,x,k Φ k,µft,x,k. 6 2
3 It is evident that, for every k and µ, Φ k,µ Fk,µ Φ k,µfk,µ = Sk,kFk,µ Fk,µdk Sk,k Fk,µ Fk,µ dk =. Using this identity, we have Qft,x,k = Φ k,µ ft,x,k Φ k,µ Fk,µ ft,x,k Fk, µ = Φ k,µ Fk,µ {Fk,µ ft,x,k Fk,µft,x,k} = Φ k,µ Fk,µ By defining the collision freuency of the BGK model κk,µ = Φ k,µ Fk,µ = Φ k,µ +exp the kinetic model writes Fk,µ ft,x,k. 7 εk µ, 8 f t t,x,k+ kε x ft,x,k e E kft,x,k = κk,µfk,µ ft,x,k, 9 κk,µfk,µ ft,x,kdk =. A solution of system 9- consists of the two functions ft,x,k and µt,x, which must satisfy suitable regularity conditions. The common definitions of the kernel S guarantees that the function κk,µ is a no negative continuous function in the set R. 2. The Pauli exclusion principle The distribution function f must be not negative and less or eual to one, according to the Pauli s exclusion principle. We prove that, if the electric field is null, then every spatial homogeneous solution of Euations 9- satisfies these conditions. Theorem Let be ft,k a function defined in,+, differentiable with respect to t for every k, and integrable, with respect to k, over the domain for each time t. Moreover the function µt is continuous in,+. If f and µ satisfy the euations f t,k = κk,µfk,µ ft,k and t κk,µfk,µ ft,kdk =, and f,k, for all k, then ft,k for t,k,+. Proof. Let be Therefore we have expkk,µ f t Kk, µ such that K t = κk,µ t. = expkk,µκk,µfk,µ expkk,µκk,µft,k 3
4 that is and then expkk,µft,k = expkk,µκk,µfk,µ. t expkk,µft,k expkk,µft,k. t= Hence the distribution function ft,k is always not negative, because f,k for every k. Since f,k for every k, we define gt,k = ft,k, and we note that g,k. It is evident that g satisfies the euation g t = κk,µ Fk,µ g and then, as before, taking into account that Fk,µ <, from expkk,µgt,k expkk, µ gt, k, t= it follows gt,k ft,k. 2.2 The kernel of the Boltzmann collision operator The kernel Sk,k is given by the sum of terms of the kind G ν k,k 2 n ν + δ εk εk+ ω ν +n ν δ εk εk ω ν. The index ν labels the νth phonon mode, G ν k,k is the scattering rate, which describes the scattering mechanism between phonons ν and electrons. The symbol δ denotes the Dirac distribution function, ω ν the νth phonon freuency, and n ν is the Bose-Einstein distribution for the phonon of type ν n ν =. e ων / For acoustic phonons, usually one considers the elastic approximation, and 2n ac G ac k,k 2 = πdack 2 B T 2π 2 2 σ m vp 2 +cosϑk,k, where D ac is the acoustic phonon coupling constant, v p is the sound speed in graphene, σ m the graphene areal density, and ϑ k,k is the convex angle between k and k. There are three relevant optical phonon scatterings: the longitudinal optical LO, the transversal optical TO and the K K phonons. The scattering rates are G LO k,k 2 = G TO k,k 2 = G K k,k 2 = 2π 2 πd 2 O σ m ω O cosϑk,k k +ϑ k,k k πd 2 O +cosϑk,k 2π 2 σ m ω k +ϑ k,k k O 2πD 2 K cosϑk,k 2π 2, σ m ω K where D O is the optical phonon coupling constant, ω O the optical phonon freuency, D K is the K- phonon coupling constant and ω K the K-phonon freuency. The angles ϑ k,k k and ϑ k,k k denote 4
5 the convex angles between k and k k and between k and k k, respectively. We used the same physical parameters of the Table of Ref. 3. Since the phonon freuency of longitudinal and transversal optical phonons coincide, then, as we sum the corresponding terms, the function cosϑ k,k k + ϑ k,k k can be eliminated, easily. Therefore, in this case, S is a sum of terms where now G ν k,k 2 = G ν cosϑ k,k. 2.3 The collision operator of the BGK model Taking into account Euation 9 and the definition 8, the collision operator of the BGK model writes { } εk µ Φ k,µ ft,x,k +exp. 2 This expression can be simplified. To this aim, we consider the function Φ k,µ, which is given by the sum of the following integrals n ν + δ +n ν δ G ν cosϑ k,k εk εk ω ν n ν + = εk+ ω ν G ν cosϑ k,k δ µ R +exp 2 n ν + G ν cosϑ k,k δ εk ω ν µ R + 2 +exp εk εk+ ω ν εk εk ω ν dk εk εk+ ω ν dk, Fk,µdk where z + = max{z,}. Introducing polar coordinates k = rcosϑ,sinϑ and k = r cosϑ,sinϑ, we have G ν cosϑ k,k δ = = = + 2π 2π dr 2π εk εk ω ν dϑ G ν cosϑ ϑδ G ν cosϑ ϑdϑ δ v F R G ν cosϑ dϑ r + ων v F v F dk where H is the Heaviside function. Analogously G ν cosϑ k,k δ v F r v F r ω ν r r ων H εk εk+ ω ν v F r+ ων v F r r Hr dr = 2π 2π dk = G ν cosϑ dϑ ν εk+ ω G ν cosϑ dϑ v F 2, εk ω ν + v F 2. 5
6 Hence Φ k,µ = ν v F 2 2π n ν + +exp G ν cosϑ dϑ εk+ ω ν εk+ ω ν n ν + µ +exp εk ω ν εk ω ν + µ. + In order to simplify the notation, we define So Φ k,µ = ν ψk,ξ;a,b = a+εk+b aεk b εk+b εk b+ +ξexp +ξexp 2π µ v F 2 G ν cosϑ dϑ ψ k,exp ;n ν, ω ν. 4 Then the BGK collision operator 2 is a linear combination, with constant positive coefficients, of the terms { } µ ψ k,exp ;n ν, ω ν εk µ ft,x,k +exp. 5 Since n ν and ω ν, then we can consider the function ψ, defined in 3, only for a, b, and ξ >. Lemma If ft,x,k, then the collision operator 2 is a strictly increasing function of µ on R. Proof. Since Euation 2 is a linear combination, with constant positive coefficients, of the functions 5, it is sufficient to prove that every function 5 is strictly increasing on R with respect to the variable µ. Again to simplify the notation, we consider the general application λk,ξ;a,b = ψk,ξ;a,b { ϕk +ξexp εk }, 6 where ϕk for all k. Euation 6 gives the expression 5, choosing the parameters a, b, ξ appropriately, and ϕk = ft,x,k, for fixed t and x. If we define εk b wk = exp and β = exp, then ψk,ξ;a,b = a+εk+b +ξβwk + aεk b + +ξβ wk, 6
7 where we have taken into account the identity aεk b + εk b+ +ξexp aεk b +. εk b +ξexp Hence we can write a+εk+b λk,ξ;a,b = + aεk b + +ξβwk +ξβ wk β ϕk+ ϕk = a+εk+b +ξβwk βϕk+ ϕk +aεk b + +ξβ wk ϕk ξwkϕk β ϕk βϕk which is strictly decreasing with respect to the positive variable ξ, because ϕk. µ Now the result follows for ξ = exp. It is useful, for the following, to establish some ineualities. From β β ϕk+ ϕk and βϕk+ ϕk β, we obtain β λk,ξ;a,b a+εk+b +ξβwk β ϕk +aεk b + +ξβ wk βϕk β a+εk+b ξβwk β ϕk +aεk b + ξβ wk βϕk, 7 and β λk,ξ;a,b a+εk+b +ξβwk β ϕk +aεk b + +ξβ wk βϕk a+εk+b a+εk+b β +βaεk b +ξβwk β + ϕk. 8 We remark that the variable ξ has not been involved in ineualities. 3 The mass conservation Euation, that guarantees the conservation of mass, can be written as { } εk µ Φ k,µ ft,x,k +exp dk =, i.e. ν 2π G ν cosϑ dϑ µ λ k,exp ;n ν, ω ν dk =. 9 ϕk=ft,x,k, 7
8 Theorem 2 If ft,x,k and εkft,x,k is integrable with respect to k over, for all t and x, then there exists a uniue µ satisfying Euation 9. Proof. If there exists a solution µ of Euation 9, then it must be uniue due to the Lemma. To prove the theorem, we show that every integrals of Euation 9 is negative for µ, and positive for µ +. We do not consider the meaningless case ft,x,k = almost everywhere. Taking into account the ineuality 7, we have λk,ξ;a,b ϕk=ft,x,k dk { a+εk+b ξβwk β ft,x,k +aεk b + = a+εk+b + aεk b + β ξ βwk β dk wk a+εk+bβ +aεk b + β ft,x,kdk, β ξβ wk βft,x,k where the last integral is finite due to hypotheses of the theorem. If µ, then ξ +, and lim λk,ξ;a,b µ ϕk=ft,x,k dk R 2 = a+εk+bβ +aεk b + β ft,x,kdk <. This implies that the l.h.s. of Euation 9 is negative for µ. Now we consider the ineuality 8. We have λk,ξ;a,b ϕk=ft,x,k dk β a+ εk+b +ξβwk dk where both integrals exist. Since, for µ >, we have εk+b + v F r +b dk = 2π r dr +ξβwk vf r µ +β exp = 2π + s+b v F 2 sds 2π µ v F 2 Therefore 2π µ v F 2 +β exp s µ s+b +β sds = 2π v F 2 +β } dk a+εk+b +βaεk b β + ft,x,kdk, 3 µ3 + 2 bµ2 lim λk,ξ;a,b µ + ϕk=ft,x,k dk = +. s+b sds s µ +β exp. This implies that the l.h.s. of Euation 9 is positive for µ +, and it concludes the proof. This theorem establishes that the kinetic model 9- is correctly-set. 8
9 4 Conclusions In this paper we propose a BGK model for charge transport in graphene. The collision operator of this model replaces the true non linear Boltzmann integral operator. A further euation, for a new unknown, is added to the kinetic euation in order to guarantee the mass conservation. The model would furnish good results when the electric field is not strong, so that the distribution function remains near an euilibrium Fermi-Dirac distribution. Moreover the simple collision term allows analytical investigations, which may be prohibitive for the Boltzmann euation. A numerical scheme, based on a discontinuous Galerkin method, for finding approximate solutions to the model, does not seem more simple than the scheme used for solving the Boltzmann euation in Ref. 5, due to the non-linearity of the integral euation for the mass conservation. Acknowledgements The author acknowledges the financial support provided by the project Modellistica, simulazione e ottimizzazione del trasporto di cariche in strutture a bassa dimensionalit, University of Catania. References A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Reviews of Modern Physics, vol. 8, no., pp. 9 62, A. Majorana, G. Mascali, and V. Romano, Charge transport and mobility in monolayer graphene, Journal of Mathematics in Industry, vol. 7, no. 4, V. Romano, A. Majorana, and M. Coco, DSMC method consistent with the Pauli exclusion principle and comparison with deterministic solutions for charge transport in graphene, Journal of Computational Physics, vol. 32, pp , A. Tomadin, D. Brida, G. Cerullo, A. C. Ferrari, and M. Polini, Noneuilibrium dynamics of photoexcited electrons in graphene: Collinear scattering, auger processes, and the impact of screening, Physical Review B, vol. 88, no. 3543, pp. 8, M. Coco, A. Majorana, and V. Romano, Cross validation of discontinuous Galerkin method and Monte Carlo simulations of charge transport in graphene on substrate, Ricerche di Matematica, vol. 66, no., pp. 2 22, 27. 9
arxiv:cond-mat/ v1 [cond-mat.stat-mech] 28 Mar 2000
arxiv:cond-mat/0003457v1 [cond-mat.stat-mech] 28 Mar 2000 High Field Mobility and Diffusivity of Electron Gas in Silicon Devices S. F. Liotta 1 and A. Majorana 2 Dipartimento di Matematica - University
More informationSupplementary Figure 2 Photoluminescence in 1L- (black line) and 7L-MoS 2 (red line) of the Figure 1B with illuminated wavelength of 543 nm.
PL (normalized) Intensity (arb. u.) 1 1 8 7L-MoS 1L-MoS 6 4 37 38 39 4 41 4 Raman shift (cm -1 ) Supplementary Figure 1 Raman spectra of the Figure 1B at the 1L-MoS area (black line) and 7L-MoS area (red
More informationA Zero Field Monte Carlo Algorithm Accounting for the Pauli Exclusion Principle
A Zero Field Monte Carlo Algorithm Accounting for the Pauli Exclusion Principle Sergey Smirnov, Hans Kosina, Mihail Nedjalkov, and Siegfried Selberherr Institute for Microelectronics, TU-Vienna, Gusshausstrasse
More informationA Solution of the Two-dimensional Boltzmann Transport Equation
Applied Mathematical Sciences, Vol. 9, 2015, no. 147, 7347-7355 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.510665 A Solution of the Two-dimensional Boltzmann Transport Equation F.
More informationFinite Difference Scheme for Semiconductor Boltzmann Equation with Nonlinear Collision Operator
Applied Mathematical Sciences, Vol. 7, 2013, no. 15, 735-750 Finite Difference Scheme for Semiconductor Boltzmann Equation with Nonlinear Collision Operator C. L. R. Milazzo Kore University of Enna Faculty
More informationSupplementary Figures
Supplementary Figures 8 6 Energy (ev 4 2 2 4 Γ M K Γ Supplementary Figure : Energy bands of antimonene along a high-symmetry path in the Brillouin zone, including spin-orbit coupling effects. Empty circles
More informationMonte Carlo Study of Thermal Transport of Direction and Frequency Dependent Boundaries in High Kn Systems
Monte Carlo Study of Thermal Transport of Direction and Frequency Dependent Boundaries in High Kn Systems N.A. Roberts and D.G. Walker Department of Mechanical Engineering Vanderbilt University May 30,
More informationarxiv: v2 [cond-mat.mes-hall] 7 Dec 2009
First Principles Analysis of Electron-Phonon Interactions in Graphene K. M. Borysenko, 1 J. T. Mullen, 2 E. A. Barry, 1 S. Paul, 2 Y. G. arxiv:0912.0562v2 [cond-mat.mes-hall] 7 Dec 2009 Semenov, 1 J. M.
More informationUnified theory of quantum transport and quantum diffusion in semiconductors
Paul-Drude-Institute for Solid State Electronics p. 1/? Unified theory of quantum transport and quantum diffusion in semiconductors together with Prof. Dr. V.V. Bryksin (1940-2008) A.F. Ioffe Physical
More information2D semiconductor device simulations by WENO-Boltzmann schemes: efficiency, boundary conditions and comparison to Monte Carlo methods
2D semiconductor device simulations by WENO-Boltzmann schemes: efficiency, boundary conditions and comparison to Monte Carlo methods José A. Carrillo ICREA - Department de Matemàtiques, Universitat Autònoma
More informationElectrical Transport. Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8
Electrical Transport Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8 Electrical Transport The study of the transport of electrons & holes (in semiconductors) under various conditions. A broad & somewhat specialized
More informationLecture 1 - Electrons, Photons and Phonons. September 4, 2002
6.720J/3.43J - Integrated Microelectronic Devices - Fall 2002 Lecture 1-1 Lecture 1 - Electrons, Photons and Phonons Contents: September 4, 2002 1. Electronic structure of semiconductors 2. Electron statistics
More informationNon-Continuum Energy Transfer: Phonons
Non-Continuum Energy Transfer: Phonons D. B. Go Slide 1 The Crystal Lattice The crystal lattice is the organization of atoms and/or molecules in a solid simple cubic body-centered cubic hexagonal a NaCl
More informationarxiv: v1 [cond-mat.mes-hall] 26 Sep 2013
Berry phase and the unconventional quantum Hall effect in graphene Jiamin Xue Microelectronic Research Center, The University arxiv:1309.6714v1 [cond-mat.mes-hall] 26 Sep 2013 of Texas at Austin, Austin,
More informationLow Bias Transport in Graphene: An Introduction
Lecture Notes on Low Bias Transport in Graphene: An Introduction Dionisis Berdebes, Tony Low, and Mark Lundstrom Network for Computational Nanotechnology Birck Nanotechnology Center Purdue University West
More informationTwo-phonon Raman scattering in graphene for laser excitation beyond the π-plasmon energy
Journal of Physics: Conference Series PAPER OPEN ACCESS Two-phonon Raman scattering in graphene for laser excitation beyond the π-plasmon energy To cite this article: Valentin N Popov 2016 J. Phys.: Conf.
More informationSupplementary Figure 1: The effect of relaxation times on the sideband. spectrum of a THz driven bi-layer graphene. Contour plot of the sideband
Supplementary Figure 1: The effect of relaxation times on the sideband spectrum of a THz driven bi-layer graphene. Contour plot of the sideband spectrum of THz driven BLG for different relaxation times.
More informationSession 5: Solid State Physics. Charge Mobility Drift Diffusion Recombination-Generation
Session 5: Solid State Physics Charge Mobility Drift Diffusion Recombination-Generation 1 Outline A B C D E F G H I J 2 Mobile Charge Carriers in Semiconductors Three primary types of carrier action occur
More informationPhonon II Thermal Properties
Phonon II Thermal Properties Physics, UCF OUTLINES Phonon heat capacity Planck distribution Normal mode enumeration Density of states in one dimension Density of states in three dimension Debye Model for
More informationSummary lecture VI. with the reduced mass and the dielectric background constant
Summary lecture VI Excitonic binding energy reads with the reduced mass and the dielectric background constant Δ Statistical operator (density matrix) characterizes quantum systems in a mixed state and
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationOptical Properties of Lattice Vibrations
Optical Properties of Lattice Vibrations For a collection of classical charged Simple Harmonic Oscillators, the dielectric function is given by: Where N i is the number of oscillators with frequency ω
More informationNumerical calculation of the electron mobility in ZnS and ZnSe semiconductors using the iterative method
International Journal o the Physical Sciences Vol. 5(11), pp. 1752-1756, 18 September, 21 Available online at http://www.academicjournals.org/ijps ISSN 1992-195 21 Academic Journals Full Length Research
More informationNUMERICAL CALCULATION OF THE ELECTRON MOBILITY IN GaAs SEMICONDUCTOR UNDER WEAK ELECTRIC FIELD APPLICATION
International Journal of Science, Environment and Technology, Vol. 1, No 2, 80-87, 2012 NUMERICAL CALCULATION OF THE ELECTRON MOBILITY IN GaAs SEMICONDUCTOR UNDER WEAK ELECTRIC FIELD APPLICATION H. Arabshahi,
More informationELECTRONIC ENERGY DISPERSION AND STRUCTURAL PROPERTIES ON GRAPHENE AND CARBON NANOTUBES
ELECTRONIC ENERGY DISPERSION AND STRUCTURAL PROPERTIES ON GRAPHENE AND CARBON NANOTUBES D. RACOLTA, C. ANDRONACHE, D. TODORAN, R. TODORAN Technical University of Cluj Napoca, North University Center of
More informationESE 372 / Spring 2013 / Lecture 5 Metal Oxide Semiconductor Field Effect Transistor
Metal Oxide Semiconductor Field Effect Transistor V G V G 1 Metal Oxide Semiconductor Field Effect Transistor We will need to understand how this current flows through Si What is electric current? 2 Back
More information2 Symmetry. 2.1 Structure of carbon nanotubes
2 Symmetry Carbon nanotubes are hollow cylinders of graphite sheets. They can be viewed as single molecules, regarding their small size ( nm in diameter and µm length), or as quasi-one dimensional crystals
More informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationIdentical Particles. Bosons and Fermions
Identical Particles Bosons and Fermions In Quantum Mechanics there is no difference between particles and fields. The objects which we refer to as fields in classical physics (electromagnetic field, field
More informationCalculating Band Structure
Calculating Band Structure Nearly free electron Assume plane wave solution for electrons Weak potential V(x) Brillouin zone edge Tight binding method Electrons in local atomic states (bound states) Interatomic
More informationLow Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model
Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,
More informationThermodynamics, Gibbs Method and Statistical Physics of Electron Gases
Bahram M. Askerov Sophia R. Figarova Thermodynamics, Gibbs Method and Statistical Physics of Electron Gases With im Figures Springer Contents 1 Basic Concepts of Thermodynamics and Statistical Physics...
More informationDensity of states for electrons and holes. Distribution function. Conduction and valence bands
Intrinsic Semiconductors In the field of semiconductors electrons and holes are usually referred to as free carriers, or simply carriers, because it is these particles which are responsible for carrying
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More informationDerivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics
Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics Luigi Barletti (Università di Firenze) Carlo Cintolesi (Università di Trieste) 6th MMKT Porto Ercole, june 9th
More informationTransport Properties of Semiconductors
SVNY85-Sheng S. Li October 2, 25 15:4 7 Transport Properties of Semiconductors 7.1. Introduction In this chapter the carrier transport phenomena in a semiconductor under the influence of applied external
More informationChapter 6 Free Electron Fermi Gas
Chapter 6 Free Electron Fermi Gas Free electron model: The valence electrons of the constituent atoms become conduction electrons and move about freely through the volume of the metal. The simplest metals
More informationGRAPHENE AS A NOVEL SINGLE PHOTON COUNTING OPTICAL AND IR PHOTODETECTOR
GRAPHENE AS A NOVEL SINGLE PHOTON COUNTING OPTICAL AND IR PHOTODETECTOR J.O.D Williams*, J.S. Lapington*, M. Roy*, I.B. Hutchinson* *Space Research Centre, Department of Physics and Astronomy, University
More informationQuantum Kinetic Transport under High Electric Fields
VLSI DESIGN 1998, Vol. 6, Nos. (1-4), pp. 3-7 Reprints available directly from the publisher Photocopying permitted by license only 1998 OPA (Ov+rseas Publishers Association) N.V. Published by license
More informationSoft Carrier Multiplication by Hot Electrons in Graphene
Soft Carrier Multiplication by Hot Electrons in Graphene Anuj Girdhar 1,3 and J.P. Leburton 1,2,3 1) Department of Physics 2) Department of Electrical and Computer Engineering, and 3) Beckman Institute
More information7. FREE ELECTRON THEORY.
7. FREE ELECTRON THEORY. Aim: To introduce the free electron model for the physical properties of metals. It is the simplest theory for these materials, but still gives a very good description of many
More informationLectures Graphene and
Lectures 15-16 Graphene and carbon nanotubes Graphene is atomically thin crystal of carbon which is stronger than steel but flexible, is transparent for light, and conducts electricity (gapless semiconductor).
More informationPhonons II - Thermal Properties (Kittel Ch. 5)
Phonons II - Thermal Properties (Kittel Ch. 5) Heat Capacity C T 3 Approaches classical limit 3 N k B T Physics 460 F 2006 Lect 10 1 Outline What are thermal properties? Fundamental law for probabilities
More informationSupplementary Figure 3: Interaction effects in the proposed state preparation with Bloch oscillations. The numerical results are obtained by
Supplementary Figure : Bandstructure of the spin-dependent hexagonal lattice. The lattice depth used here is V 0 = E rec, E rec the single photon recoil energy. In a and b, we choose the spin dependence
More informationModeling of noise by Legendre polynomial expansion of the Boltzmann equation
Modeling of noise by Legendre polynomial expansion of the Boltzmann equation C. Jungemann Institute of Electromagnetic Theory RWTH Aachen University 1 Outline Introduction Theory Noise Legendre polynomial
More informationIntroduction to cold atoms and Bose-Einstein condensation (II)
Introduction to cold atoms and Bose-Einstein condensation (II) Wolfgang Ketterle Massachusetts Institute of Technology MIT-Harvard Center for Ultracold Atoms 7/7/04 Boulder Summer School * 1925 History
More informationOlivier Bourgeois Institut Néel
Olivier Bourgeois Institut Néel Outline Introduction: necessary concepts: phonons in low dimension, characteristic length Part 1: Transport and heat storage via phonons Specific heat and kinetic equation
More informationWe already came across a form of indistinguishably in the canonical partition function: V N Q =
Bosons en fermions Indistinguishability We already came across a form of indistinguishably in the canonical partition function: for distinguishable particles Q = Λ 3N βe p r, r 2,..., r N ))dτ dτ 2...
More informationCondensed Matter Physics 2016 Lecture 13/12: Charge and heat transport.
Condensed Matter Physics 2016 Lecture 13/12: Charge and heat transport. 1. Theoretical tool: Boltzmann equation (review). 2. Electrical and thermal conductivity in metals. 3. Ballistic transport and conductance
More informationUltrafast Electron-Electron Dynamics in Graphene Daniele Brida
Ultrafast Electron-Electron Dynamics in Graphene Ultrafast Spectroscopy Light-matter interaction on timescale of fundamental physical processes: electron scattering, phonon emission, energy transfer D.
More informationLecture 11: Coupled Current Equations: and thermoelectric devices
ECE-656: Fall 011 Lecture 11: Coupled Current Euations: and thermoelectric devices Professor Mark Lundstrom Electrical and Computer Engineering Purdue University, West Lafayette, IN USA 9/15/11 1 basic
More informationCovariant Electromagnetic Fields
Chapter 8 Covariant Electromagnetic Fields 8. Introduction The electromagnetic field was the original system that obeyed the principles of relativity. In fact, Einstein s original articulation of relativity
More informationMolecular Dynamics Study of Thermal Rectification in Graphene Nanoribbons
Molecular Dynamics Study of Thermal Rectification in Graphene Nanoribbons Jiuning Hu 1* Xiulin Ruan 2 Yong P. Chen 3# 1School of Electrical and Computer Engineering and Birck Nanotechnology Center, Purdue
More informationMinimal Update of Solid State Physics
Minimal Update of Solid State Physics It is expected that participants are acquainted with basics of solid state physics. Therefore here we will refresh only those aspects, which are absolutely necessary
More informationSergey SMIRNOV a),hanskosina, and Siegfried SELBERHERR, Nonmembers
350 PAPER Special Issue on the IEEE International Conference on SISPAD 02 Investigation of the Electron Mobility in Strained Si 1 x Ge x at High Ge Composition Sergey SMIRNOV a),hanskosina, and Siegfried
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationSteady-State Electron Transport in Silicon Dioxide Employing Different Electronic Band-Structures
Steady-State Electron Transport in Silicon Dioxide Employing Different Electronic Band-Structures M. Hackel, H. Kosina, and S. Selberherr Institute for Microelectronics Technical University of Vienna Gusshausstrasse
More informationMONTE CARLO SIMULATION FOR ELECTRON DYNAMICS IN SEMICONDUCTOR DEVICES
Mathematical and Computational Applications, Vol., No., pp. 9-26, 25. Association for Scientific Research MONTE CARLO SIMULATION FOR ELECTRON DYNAMICS IN SEMICONDUCTOR DEVICES Mustafa Akarsu, Ömer Özbaş
More informationSemiconductor Physical Electronics
Semiconductor Physical Electronics Sheng S. Li Department of Electrical Engineering University of Florida Gainesville, Florida Plenum Press New York and London Contents CHAPTER 1. Classification of Solids
More informationPhysics 408 Final Exam
Physics 408 Final Exam Name You are graded on your work, with partial credit where it is deserved. Please give clear, well-organized solutions. 1. Consider the coexistence curve separating two different
More informationGiant Intrinsic Carrier Mobilities in Graphene and Its Bilayer
Giant Intrinsic Carrier Mobilities in Graphene and Its Bilayer S.V. Morozov 1,, K.S. Novoselov 1, M.I. Katsnelson, F. Schedin 1, D.C. Elias 1, J.A. Jaszczak, A.K. Geim 1 1 Manchester Centre for Mesoscience
More informationQuantum Statistical Derivation of a Ginzburg-Landau Equation
Journal of Modern Physics, 4, 5, 56-568 Published Online October 4 in SciRes. http://www.scirp.org/journal/jmp http://dx.doi.org/.46/jmp.4.5657 Quantum Statistical Derivation of a Ginzburg-Landau Euation
More informationSCIENCE & TECHNOLOGY
Pertanika J. Sci. & Technol. 25 (S): 205-212 (2017) SCIENCE & TECHNOLOGY Journal homepage: http://www.pertanika.upm.edu.my/ Effect of Boron and Oxygen Doping to Graphene Band Structure Siti Fazlina bt
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationKey Words Student Paper, Electrical Engineering
A GUI Program for the Calculation of Mobility and Carrier Statistics in Semiconductors Daniel Barrett, University of Notre Dame Faculty Advisor: Dr. Debdeep Jena, University of Notre Dame Student Paper
More informationAppendix 1: List of symbols
Appendix 1: List of symbols Symbol Description MKS Units a Acceleration m/s 2 a 0 Bohr radius m A Area m 2 A* Richardson constant m/s A C Collector area m 2 A E Emitter area m 2 b Bimolecular recombination
More informationModel Question Paper ENGINEERING PHYSICS (14PHY12/14PHY22) Note: Answer any FIVE full questions, choosing one full question from each module.
Model Question Paper ENGINEERING PHYSICS (14PHY1/14PHY) Time: 3 hrs. Max. Marks: 100 Note: Answer any FIVE full questions, choosing one full question from each module. MODULE 1 1) a. Explain in brief Compton
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2014 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 POP QUIZ Phonon dispersion relation:
More informationTheory of thermopower in two-dimensional graphene
PHYSICAL REVIEW B 8, 235415 29 Theory of thermopower in two-dimensional graphene E. H. Hwang, E. Rossi, and S. Das Sarma Condensed Matter Theory Center, Department of Physics, University of Maryland, College
More informationExplicit expression for Lindhard dielectric function at finite temperature
arxiv:1412.5705v1 [cond-mat.mtrl-sci] 18 Dec 2014 Explicit expression for Lindhard dielectric function at finite temperature A. V. Andrade-Neto Departamento de Física, Universidade Estadual de Feira de
More informationGraphene and Planar Dirac Equation
Graphene and Planar Dirac Equation Marina de la Torre Mayado 2016 Marina de la Torre Mayado Graphene and Planar Dirac Equation June 2016 1 / 48 Outline 1 Introduction 2 The Dirac Model Tight-binding model
More informationAn Effective Potential Approach to Modelling 25nm MOSFET Devices S. Ahmed 1 C. Ringhofer 2 D. Vasileska 1
An Effective Potential Approach to Modelling 25nm MOSFET Devices S. Ahmed 1 C. Ringhofer 2 D. Vasileska 1 1 Department of Electrical Engineering, 2 Department of Mathematics, Arizona State University,
More information3. LATTICE VIBRATIONS. 3.1 Sound Waves
3. LATTIC VIBRATIONS Atoms in lattice are not stationary even at T 0K. They vibrate about particular equilibrium positions at T 0K ( zero-point energy). For T > 0K, vibration amplitude increases as atoms
More informationChapter 6 Free Electron Fermi Gas
Chapter 6 Free Electron Fermi Gas Free electron model: The valence electrons of the constituent atoms become conduction electrons and move about freely through the volume of the metal. The simplest metals
More informationHomework 2: Solutions
Homework : Solutions. Problem. We have approximated the conduction band of many semiconductors with an isotropic, parabolic expression of the form: E(k) h k m. () However,this approximation fails at some
More informationRaman Imaging and Electronic Properties of Graphene
Raman Imaging and Electronic Properties of Graphene F. Molitor, D. Graf, C. Stampfer, T. Ihn, and K. Ensslin Laboratory for Solid State Physics, ETH Zurich, 8093 Zurich, Switzerland ensslin@phys.ethz.ch
More informationGraphene films on silicon carbide (SiC) wafers supplied by Nitride Crystals, Inc.
9702 Gayton Road, Suite 320, Richmond, VA 23238, USA Phone: +1 (804) 709-6696 info@nitride-crystals.com www.nitride-crystals.com Graphene films on silicon carbide (SiC) wafers supplied by Nitride Crystals,
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationJournal of Atoms and Molecules
Research article Journal of Atoms and Molecules An International Online Journal ISSN 77 147 Hot Electron Transport in Polar Semiconductor at Low Lattice Temperature A. K. Ghorai Physics Department, Kalimpong
More informationBasic Semiconductor Physics
Chihiro Hamaguchi Basic Semiconductor Physics With 177 Figures and 25 Tables Springer 1. Energy Band Structures of Semiconductors 1 1.1 Free-Electron Model 1 1.2 Bloch Theorem 3 1.3 Nearly Free Electron
More informationThe Boltzmann Equation and Its Applications
Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC
More informationarxiv: v1 [cond-mat.mes-hall] 15 Apr 2007
Non-linear electromagnetic response of graphene S. A. Mikhailov Institute for Theoretical Physics II, University of Augsburg, D-86135 Augsburg, Germany (Dated: May 31, 218) arxiv:74.199v1 [cond-mat.mes-hall]
More informationExponential methods for kinetic equations
Exponential methods for kinetic equations Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara, Italy http://utenti.unife.it/lorenzo.pareschi/ lorenzo.pareschi@unife.it Joint research
More informationUnit III Free Electron Theory Engineering Physics
. Introduction The electron theory of metals aims to explain the structure and properties of solids through their electronic structure. The electron theory is applicable to all solids i.e., both metals
More informationELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS
ELECTRONS AND PHONONS IN SEMICONDUCTOR MULTILAYERS В. К. RIDLEY University of Essex CAMBRIDGE UNIVERSITY PRESS Contents Introduction 1 Simple Models of the Electron-Phonon Interaction 1.1 General remarks
More informationValley Zeeman effect in elementary optical excitations of monolayerwse 2
Valley Zeeman effect in elementary optical excitations of monolayerwse 2 Ajit Srivastava 1, Meinrad Sidler 1, Adrien V. Allain 2, Dominik S. Lembke 2, Andras Kis 2, and A. Imamoğlu 1 1 Institute of Quantum
More informationPhysics 607 Exam 2. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physics 607 Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all your
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2009 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 Spring 2009 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 Spring 2009 POP QUIZ
More informationPHZ7427 Spring 2014 Scattering mechanisms and electronic transport in conductors
PHZ7427 Spring 24 Scattering mechanisms and electronic transport in conductors Dmitrii L. Maslov (Dated: March, 24) CONTENTS A. Units and notations 2 I. Boltzmann equation 3 A. Example: conductivity in
More informationSheng S. Li. Semiconductor Physical Electronics. Second Edition. With 230 Figures. 4) Springer
Sheng S. Li Semiconductor Physical Electronics Second Edition With 230 Figures 4) Springer Contents Preface 1. Classification of Solids and Crystal Structure 1 1.1 Introduction 1 1.2 The Bravais Lattice
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationThermoelectrics: A theoretical approach to the search for better materials
Thermoelectrics: A theoretical approach to the search for better materials Jorge O. Sofo Department of Physics, Department of Materials Science and Engineering, and Materials Research Institute Penn State
More informationSupporting Information
Electronic Supplementary Material (ESI) for Energy & Environmental Science. This journal is The Royal Society of Chemistry 2017 Supporting Information Large Enhancement of Thermoelectric Properties in
More informationIntroduction. Chapter The Purpose of Statistical Mechanics
Chapter 1 Introduction 1.1 The Purpose of Statistical Mechanics Statistical Mechanics is the mechanics developed to treat a collection of a large number of atoms or particles. Such a collection is, for
More informationNonlinear Electrodynamics and Optics of Graphene
Nonlinear Electrodynamics and Optics of Graphene S. A. Mikhailov and N. A. Savostianova University of Augsburg, Institute of Physics, Universitätsstr. 1, 86159 Augsburg, Germany E-mail: sergey.mikhailov@physik.uni-augsburg.de
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationSiC Graphene Suitable For Quantum Hall Resistance Metrology.
SiC Graphene Suitable For Quantum Hall Resistance Metrology. Samuel Lara-Avila 1, Alexei Kalaboukhov 1, Sara Paolillo, Mikael Syväjärvi 3, Rositza Yakimova 3, Vladimir Fal'ko 4, Alexander Tzalenchuk 5,
More informationELECTRON SCATTERING BY SHORT RANGE DEFECTS AND RESISTIVITY OF GRAPHENE
ELECTRON SCATTERING BY SHORT RANGE DEECTS AND RESISTIVITY O GRAPHENE Natalie E. irsova *, Sergey A. Ktitorov * *The Ioffe Physical-Technical Institute of the Russian Academy of Sciences, 6 Politekhnicheskaya,
More informationSupplementary Information Supplementary Figures
Supplementary Information Supplementary Figures Supplementary Fig S1: Multilayer MoS 2 FETs on SiO2/Si substrates, and contact resistance effects. (Left): Transfer curves and effective mobility of multilayer
More information