Density functional theory
|
|
- Gavin Thornton
- 6 years ago
- Views:
Transcription
1 Density functional theory Seminar Ib - 1. letnik, II. stopnja Author: Matic Pečovnik Mentor: Doc. Dr. Tomaž Rejec Faculty of Mathematics and Physics, University in Ljubljana Abstract Many-body problems in physics have been the focal point of quantum mechanics in the previous decades. Such problems are found in solid state physics, complex molecules and so on. Most of the theoretical work must be done numerically. Here we have an issue with the computational power of computers as the solving of the Schroedinger many-body equation requires computational power that scales exponential with the number of particles in the system. In this seminar I present to the reader the need to transition from wave-centered to density-centered approach when performing theoretical analysis on complex many-body systems. Then I present the basic concepts of the density functional theory (DFT) first introduced by Kohn and Sham that offers a solution to problems presented above. The second half of my seminar will be focused on simulating the behaviour of a two dimensional electron gas with the use of density functional theory. 1
2 Contents 1 Introduction to density functional theory Wave-centered approach Density functional theory - A workaround the an leck catastrophe The self-consistent Kohn-Sham equations Uniqueness of the external potential - Hohenberg-Kohn Theorem Obtaining the ground state Local density approximation Improvements of the local density approximation Numerical packages for DFT calculations An application of DFT - 2D electron gas 7 4 Conclusion Introduction to density functional theory In this first chapter I will illuminate the problem we encounter when dealing with the many-body wave function. Then I will introduce a possible solution in the form of density functional theory (DFT) which plays a major role in numerically simulating the behaviour of many-body systems in modern day physics Wave-centered approach Let D be the dimensionality of our system. In wave-centered computation the wave function derived from the Schroedinger s equation is the main body of focus. Ψ = Ψ(r 1, r 2,..., r N ) (1) It contains N coordinates r i = (x i, y i, z i ; σ i ) so as a whole it contains (D + 1)N parameters with which we describe the system. It has to satisfy the normalization restraint: dr 1 dr 2... Ψ(r 1, r 2,..., r N ) 2 dr N = 1 (2) It has to also satisfy the antisymmetricity of a fermion function(pauli principle): P i,j Ψ = Ψ (3) Here P i,j is the parity operator that switches the i-th and j-th fermion. In systems with only a couple of particles the wave-centered approach performs really well at predicting the results of experiments and has therefore become the pillar of quantum mechanics but it faces huge problems with many-body systems such as we find in nature. Molecules, crystals, DNA molecules and other complex matter are all described with the Schroedinger equation but as we can not numerically solve it on a realistic time scale it becomes useless. Let us demonstrate this on an example. 2
3 In this example we will consider a molecule with N interacting electrons. We will disregard all spin degrees of freedom and ignore symmetries. In numerical calculations we normally minimize a functional of energy in accordance to the Rayleigh-Ritz minimal principle so we obtain the ground state of a system. E = min Ψ Ψ H Ψ (4) As we ignore the spin variable of the electrons the wave function consists of 3N continuous variables which we have to represent discretely by M parameters p. In this way the energy is a function of the parameters E = E(p 1, p 2,..., p M ). In practical calculations it turns out that we require 3 n 10 parameters per variable for the desired accuracy od δe We therefore need: M = n 3N 3 n 10 (5) parameters of the wave function. If M = 10 9 is the maximum computable number of parameters and n = 3 then this means the maximum number of computable electrons is: N = 1 3 ln M n. = 6 (6) With the same thinking we can compute that for a system with a 100 electrons we would need M = = which is more than the predicted baryons in the whole universe [1]. This is what we call the an leck catastrophe Density functional theory - A workaround the an leck catastrophe The first density functional theory was published by Thomas [2] and Fermi [3] in the 1920s. In the centre of this approach was the density of electrons n(r) that is also easily measurable by various techniques. It was not accurate enough as it neglects the quantum corrections to the Couloumb potential which I will later introduce as exchange and correlation. It proposed an expected value of the kinetic term in the H operator that depends solely on the density. This approach is acceptable only in high electron density metals where the quantum corrections to the Couloumb potential can be neglected. Only in the 1960s a new complete formalism was proposed by Kohn and Sham [4] that sparked a new era in physics. It provides many results that are in accordance to results from experiments in various fields such as crystals, biophysics, molecular physics and many more. These systems would be otherwise inaccessible by numerical computations with wave functions as these systems contain N 1 electrons. 2. The self-consistent Kohn-Sham equations In the centre of our approach will be the density of electrons n(r) as the ground state n(r) uniquely (to a trivial constant) defines the external potential v(r). This means that if we know the external potential we in theory know n(r) and vice-versa. In the case of systems with degenerate ground states any degenerate ground state n i (r) uniquely describes the same v(r) - This is widely known as the Hohenberg-Kohn Theorem [5]. 3
4 2.1. Uniqueness of the external potential - Hohenberg-Kohn Theorem Let n(r) be the ground-state of a system described by the hamiltonian H 1. It corresponds to a ground-state wave function Ψ 1 with energy E 1, which we assume is non-degenerate. E 1 = Ψ 1 H 1 Ψ 1 = v 1 (r)n(r) d + Ψ 1 T + U Ψ 1 (7) Here U is the interaction between electrons and T is the kinetic term. Now lets assume that there exists a second potential v 2 (r) = v 1 (r) + v 0 with its own ground-state wave function Ψ 2 = exp iθ Ψ 1 that gives the same density n(r). We can write the energy as: E 2 = Ψ 2 H 2 Ψ 2 = v 2 (r)n(r) d + Ψ 2 T + U Ψ 2 (8) Since the system is non-degenerate the Rayleigh-Ritz minimal principle applies. E 1 < Ψ 2 H 1 Ψ 2 = With the same reasoning we can also conclude: E 2 < Ψ 1 H 2 Ψ 1 = If we add these equations we get a contradiction. Therefore it must apply that v 1 (r) = v 2 (r). [1] v 1 (r)n(r) d + Ψ 2 T + U Ψ 2 = E 2 + (v 1 (r) v 2 (r))n(r) d (9) v 2 (r)n(r) d + Ψ 1 T + U Ψ 1 = E 1 + (v 2 (r) v 1 (r))n(r) d (10) E 1 + E 2 < E 1 + E 2 (11) 2.2. Obtaining the ground state In this subsection we will be considering an inhomogeneous gas of electrons without spin. We can write the energy functional E[n(r)] for this system as: E[n(r)] = v(r)n(r) n(r)n(r ) r r dr dr + F[n(r)] (12) The equation is written in standard atomic units, meaning the all distances are measured in Bohrs radius (a 0 ) and energies in Hartrees (Ha). Ha = 2Ry = m ee 4 (2ɛ 0 h) 2 = 27.2e a 0 = ɛ 0h 2 = 53pm (13) πm e e2 The first part of the functional presents the interactions of electrons and the external potential and the second part is the classical Coulomb interaction between the electrons also called the Hartree interaction. F[n(r)] = T s [n(r)] + E xc [n(r)] is a functional uniquely described by the electron density and is independent of the system that is defined by the external potential v(r). T s [n(r)] is the functional of kinetic energy for a system of non-interacting electrons with the same density and E xc [n(r)] = E x [n(r)] + E c [n(r)] is by our definition the functional that describes exchange E x [n(r)] and correlations E c [n(r)] between electrons and contains everything about not already explicitly written out. 4
5 We are already familiar with the exchange part from our classes in Condensed matter on the matter of the Hartree-Fock approximation and we know the exchange comes from using only one Slater determinant to describe our wave-function. The correlation part comes from using more Slater determinants to describe the interaction between electrons. Let f xc be the functional derivative of E xc. So f xc (r) = δe xc[n(r)] δn(r). The functional E xc [n(r)] is not analytically known and neither is f xc (r). Because we have written the functional in a way where we separate the non-interacting kinetic term T s [n(r)] from the interacting kinetic term which is hidden in E xc [n(r)] we know that the minimum of the functional will be the solution of the corresponding Schroedinger equation which is in this context knows as the Kohn-Sham equation. ( 1 ) n(r ) v(r) + r r dr + f xc(r) ɛ i ϕ i (r) = 0 (14) ɛ i is a Lagrange multiplication factor ensuring the conservation of particles. Because the effective potential, in which the interacting kinetic energy is hidden, is dependent on n(r) we must solve this equation self-consistently meaning that: n(r) = N ϕ i (r) 2 n(r) dr = N (15) i=1 Here the sum goes over all electrons in the system and ϕ i are filled in the ascending order of their respective energies ɛ i. As of yet there is no clear indication of the meaning of the Kohn-Sham orbitals ϕ i and their energies ɛ i Local density approximation To use this formalism effectively we must come by some approximation for the functional E xc [n(r)]. The first approximation is the so-called local density approximation or LDA made by Kohn and Sham in their original paper [4]. It assumes that the density n(r) is slowly varying on the scale of the Fermi wavelength that we derive from k F which corresponds to the radius of the Fermi sphere. It also works very well if the density is very high because the Fermi vector is small in this regime and the variations of density on this scale are therefore small. We can express the LDA in the following way: Exc LDA [n(r)] = n(r)(e x [n(r)] + e c [n(r)]) dr (16) e xc [n(r)] = e x [n(r)] + e c [n(r)] is the exchange-correlation energy per particle of a uniform electron gas of density n. This functional is well known. The exchange (e x ) can be analytically calculated [6] and we can estimate the correlation (e c ) with an analytical fit obtained with interpolation. A rough estimate is presented here: e x = r s +... e c = 0.43 r s (17) Here r s is defined as the radius of a sphere that contains one electrons r s = ( 4πn 3 ) 3 1 and is a dimensionless quantity. In recent decades a better approximation for the correlation per electron was obtained by using an analytical fit for the data gathered by Monte-Carlo simulations that has the required behaviour when r s 1 and r s 1 which are analytically known. 5
6 2.4. Improvements of the local density approximation As we have seen in the previous section we obtain adequate results with the LDA approximation if the system has a very high density or if the local fluctuations of density are very small. But in realistic example this is rarely the case. So we need a new better approximation to describe these systems. At the centre of our approach in LDA was that we can write the energy functional simply as a functional of density as we can see in equation 16 but we must realize this is only the first order of expansion of the general functional E xc that describes effects of correlation and exchange between electrons. We can get a better approximation with the next orders in the expansion series that also include the change of density with n(r). We can write the expansion to the second order as: Exc GGA [n(r)] = f (n(r), n(r)) dr (18) This second order is called the Generalized Gradient Approximation (GGA for short). For it to be of any use we must find a analytical fit to the f (n(r), n(r)) function. The best choice is still a matter of debate. We can see the improvement that we obtain with the GGA approximation in Fig. 1. Figure 1: Errors of bond lengths computation for a GGA functional and LDA functional. We can see that GGA is better or comparable than LDA in these calculations. Taken from [7]. Also we can generalize the DFT by including a finite temperature. You can see the whole approach in [4], while here I will only cite the results. The Kohn-Sham equation is the same only how we self-consistently produce the density changes: n(r) = N i=1 ϕ i (r) 2 1 exp ɛ i µ kt +1 (19) k is the Boltzmann constant, T is the temperature and µ is the chemical potential derived from the number of particles by integrating the equation 19. We can see that the only change is the Fermi function. These approaches can also be further generalized to study excitation energies, photo-absorption spectra and frequency responses using time-dependent external potentials like oscillating electric fields and magnetic fields in a theory we call time-dependent density functional theory - TD-DFT for short. From experience it is known that the accuracy of modelling the energy of bonds in molecules is about 10-20% and the sizes of the bonds is about 1%[1] as we can see in Fig. 1. We get the 6
7 bond lengths by introducing these lengths as parameters into the energy functional E[n(r)] and minimizing the corresponding energy function of parameters by using the Rayleigh-Ritz principle. We can not make a general statement for the quality of the DFT approach for all systems. For instance in semiconductor systems the DFT gives quantitatively wrong results for the energy gap while accurately predicting the energy band structure. In some systems like two dimensional electron gas in specific cases the LDA is good enough to get quantitative results. Also systems with very strong electron-electron interactions are better described by dynamical mean-field theory [8] Numerical packages for DFT calculations There are many numerical packages that one can acquire for DFT calculations. To know which package suits one best we should note the important differences between these packages. Electrons We can divide the programs by how they treat electrons far from the Fermi surface. Some consider all electrons as equally important for understanding the behaviour of the system. Others use pseudo-potentials. When using pseudo-potentials we consider only the electrons near the Fermi surface accurately while the electrons far bellow the Fermi level are considered only as a background that effects the electrons near the Fermi surface. Basis set The numerical calculations are done in some basis set of functions. There are two widely used variants. The Gaussian basis set (SIESTA) is local and the functions differ if we are in solids or molecules. The plane-wave basis set (Quantum-Espresso, ASP, ABINIT) is non-local and the functions are universal in all systems. 3. An application of DFT - 2D electron gas We can fabricate a 2D electron gas at the junction of two semiconductors. One such example is at the heterojunction of AlGaAs and GaAs. The semiconductors have a different energy gap between the valence and conducting bands and also their chemical potentials are not equal. When we put the semiconductors in contact a current starts to flow, because of different chemical potentials. It flows until both chemical potentials are equal. The result is that the energy bands bend to accommodate this requirement. In this process of bending a potential well about 10 nm wide forms at the junction and some electrons get trapped as their movement in the direction perpendicular to the junction gets quantized. The chemical potential is usually lower then the energy of the first excited state so electrons get trapped in the ground state of the well. Therefore their movement is reduced to two dimensions in which they are almost completely free. We can see the junction in Fig. 2. This system of electrons is of particular interest to scientist because of its high mobility in the two dimensions in which the electrons are free as it is comparable or even higher than in most metals. We can study this system using a quantum point contact (QPC) which is a narrow constriction made by the potential of two metallic gate electrodes that separate two electron reservoirs of 2D electron gas as see in Fig. 3a. The gate electrodes lie about 100 nm above the electron gas and are separated from it by an insulator so the entire effect on the electron gas is the potential the electrodes form. With such a device we can measure the conductance G in the QPC, by measuring 7
8 Figure 2: A junction of two semiconductors. On the top we can see the energy band structure before we create the junction. The red lines represent their respective chemical potentials. In the lower part of the figure we can see the energy band structure after we create the junction with the chemical potential in equilibrium. We can see the trapped electrons depicted in blue. the current through such a device, that is proportional to the conductance, while changing the voltage: I = G (a) Figure 3: (a) An image of a QPC through a microscope. (b) An ELS (highlighted with the white circle) that formed in a QPC calculated by DFT using LSDA approximation. Taken from [9]. In the constriction a saddle in the potential forms. When the Fermi energy of the electrons reaches the value of the potential saddle the density of electrons at the saddle is low and high in the reservoirs. Now lets look at a 2D jellium system, where the density of electrons is constant, in the 2D Hartree-Fock approximation. Also let the electrons have spin. We can define the polarization of such a system as ζ = n n n +n where n, are the densities of electrons with spin up or down. We can write the energy of such a system with the parameters r s and ζ [10] as: E(r s, ζ) = 1 + ζ2 r 2 s 4 2 3π (b) (1 ζ) (1 + ζ) 3 2 r s (20) 8
9 Lets look at both limits for polarization ζ = 0 and ζ = ±1. Their respective energies are: E(r s, ζ = 0) = 1 r 2 s 8 2 3π 1 r s E(r s, ζ = ±1) = 2 r 2 s π 1 (21) r s We can see that for low densities (r s 1) the system will prefer to be polarized because the energy of the polarized state is lower. Now lets return to the QPC. Even though this system is not consistent with the jellium model as the density is not constant we expect that we should get qualitatively similar behaviour. We said that when the Fermi energy approaches the potential of the saddle the density there will be low so we expect that we will get a spin-polarized local state at the saddle. We can study systems where spin is relevant with DFT using a generalization of LDA named local spin density approximation or LSDA where we can write the exchange and correlation energies as: Exc LDA [n (r), n (r)] = n(r)e xc (n (r), n (r)) dr (22) This DFT approximation also adds correlations previously mentioned in Sec. 2.3 to the Hartree- Fock approximation we studied earlier. An article [9] has reported that by using LSDA the authors have found a single electron that gets caught in the constriction and so forms a localized spin-polarized state often named emergent localized state - ELS. We can see such an ELS in Fig. 3b. Now we can generalize our findings from this relatively simple system to a more complex system. To get high values of conductivity in this system one of the mentioned semiconductors is doped with a n-type material like Si. These dopants contribute even more electrons into the 2D electron gas but it comes at a cost. The dopants form a random potential seen in Fig. 4. Such a disordered potential also contains saddles on which ELS could form. We can see these ELS in Fig. 4. Figure 4: Upper left: A random disordered potential. Other pictures: Polarization calculated with the Octopus code using the LSDA approximation as a function of the space coordinates for various electron densities represented by the number of electrons N on top of every picture. For different densities different ELS are prevalent. This is because an ELS forms when the Fermi energy is the same as the energy of the saddle in the potential. Taken from [11]. 9
10 Such systems are still under study as we want to understand how the disordered potential influences the conductivity of the 2D electron gas. 4. Conclusion The transition from the wave-centered approach to the density approach is paramount to being able to quantify complex systems such as crystals and molecules as it allows us to predict behaviours of these complex systems. With this new-found understanding we can put these systems to use from nano physics such as the conductance of a 2D electron gas on the border of two semiconductors to various other examples not mentioned in this seminar. It is also very useful in chemistry and biology. Many improvements can still be made to the DFT formalism as the exact functional that describes the kinetic energy of interacting particles and the exchange and correlation of these particles is not known yet. It may never be known exactly but more elaborate methods of approximating this functional are to be expected in the future as more precision is demanded of the DFT formalism. References [1] W. Kohn, Nobel lecture: Electronic structure of matter wave functions and density functionals, Reviews of Modern Physics, vol. 71, pp , October [2] L. H. Thomas, The calculation of atomic fields, Proceedings of the Cambridge Philosophical Society, vol. 23, no. 5, p , [3] E. Fermi, Un metodo statistico per la determinazione di alcune prioprietà dell atomo, Rend. Accad. Naz. Lincei., vol. 6, pp , [4] W. Kohn and L. J. Sham, Self-consistent equations including exchange and correlation effects, Physical Review Letters, vol. 140, pp , November [5] P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Physical Review Letters, vol. 136, November [6] G. D. Mahan, Many particle physics. Springer, [7] M. Fuchs, Comparison of exchange-correlation functionals: from lda to gga and beyond, in A Hands-On Computer Course, [8] D. ollhardt, K. Byczuk, and M. Kollar, Dynamical Mean-Field Theory. Springer, [9] T. Rejec and Y. Meir, Magnetic impurity formation in quantum point contacts, Nature, vol. 442, pp , August [10] B. Tanatar and D. M. Ceperley, Ground state of the two-dimensional electron gas, Physical Review B, [11] R. Levy, Spontaneous Formation of Magnetic Moments and Dephasing in Two-Dimensional Disordered Systems. PhD thesis, Ben-Gurion University of the Negev,
1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationAll electron optimized effective potential method for solids
All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22
More informationThe electronic structure of materials 2 - DFT
Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationIntroduction to Density Functional Theory
1 Introduction to Density Functional Theory 21 February 2011; V172 P.Ravindran, FME-course on Ab initio Modelling of solar cell Materials 21 February 2011 Introduction to DFT 2 3 4 Ab initio Computational
More information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationOptimized Effective Potential method for non-collinear Spin-DFT: view to spindynamics
Optimized Effective Potential method for non-collinear Spin-DFT: view to spindynamics S. Sharma and E. K. U. Gross nstitut für Theoretische Physik, FU Berlin Fritz-Haber nstitut, Berlin 21 Sept 2006 Synopsis
More informationDensity Functional Theory (DFT)
Density Functional Theory (DFT) An Introduction by A.I. Al-Sharif Irbid, Aug, 2 nd, 2009 Density Functional Theory Revolutionized our approach to the electronic structure of atoms, molecules and solid
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationCLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122
CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationKey concepts in Density Functional Theory
From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationKohn Sham density functional theory [1 3] is. Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW.
Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW Department of Physics and Quantum Theory Group, Tulane University, New Orleans, Louisiana 70118 Received 11 March 1999;
More informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationIs the homogeneous electron gas homogeneous?
Is the homogeneous electron gas homogeneous? Electron gas (jellium): simplest way to view a metal homogeneous and normal Hartree-Fock: simplest method for many-electron systems a single Slater determinant
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationJournal of Theoretical Physics
1 Journal of Theoretical Physics Founded and Edited by M. Apostol 53 (2000) ISSN 1453-4428 Ionization potential for metallic clusters L. C. Cune and M. Apostol Department of Theoretical Physics, Institute
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationIndependent electrons in an effective potential
ABC of DFT Adiabatic approximation Independent electrons in an effective potential Hartree Fock Density Functional Theory MBPT - GW Density Functional Theory in a nutshell Every observable quantity of
More informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_
More informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
More informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationNotes on Density Functional Theory
Notes on Density Functional Theory Rocco Martinazzo E-mail: rocco.martinazzo@unimi.it Contents 1 Introduction 1 Density Functional Theory 7 3 The Kohn-Sham approach 11 1 Introduction We consider here a
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More informationAdvanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University
Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation
More informationGEM4 Summer School OpenCourseWare
GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Molecular Mechanics by Ju Li. Given August 9, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationADIABATIC CONNECTION AND UNIFORM DENSITY SCALING: CORNERSTONES OF A DENSITY FUNCTIONAL THEORY
ADIABATIC CONNECTION AND UNIFORM DENSITY SCALING: CORNERSTONES OF A DENSITY FUNCTIONAL THEORY BY RUDOLPH J. MAGYAR A dissertation submitted to the Graduate School Piscataway Rutgers, The State University
More informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
More informationA. F. J. Levi 1 EE539: Engineering Quantum Mechanics. Fall 2017.
A. F. J. Levi 1 Engineering Quantum Mechanics. Fall 2017. TTh 9.00 a.m. 10.50 a.m., VHE 210. Web site: http://alevi.usc.edu Web site: http://classes.usc.edu/term-20173/classes/ee EE539: Abstract and Prerequisites
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » Recap of HK theorems and KS equations» The physical meaning of the XC energy» Solution of a one-particle Schroedinger equation» Pseudo Potentials»
More informationBohr s Model, Energy Bands, Electrons and Holes
Dual Character of Material Particles Experimental physics before 1900 demonstrated that most of the physical phenomena can be explained by Newton's equation of motion of material particles or bodies and
More informationElectronic Structure: Density Functional Theory
Electronic Structure: Density Functional Theory S. Kurth, M. A. L. Marques, and E. K. U. Gross Institut für Theoretische Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany (Dated:
More informationMesoscopic Nano-Electro-Mechanics of Shuttle Systems
* Mesoscopic Nano-Electro-Mechanics of Shuttle Systems Robert Shekhter University of Gothenburg, Sweden Lecture1: Mechanically assisted single-electronics Lecture2: Quantum coherent nano-electro-mechanics
More informationChapter 3 Properties of Nanostructures
Chapter 3 Properties of Nanostructures In Chapter 2, the reduction of the extent of a solid in one or more dimensions was shown to lead to a dramatic alteration of the overall behavior of the solids. Generally,
More informationKey concepts in Density Functional Theory (II) Silvana Botti
Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:
More informationFundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering
Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)
More informationElectronic structure of correlated electron systems. Lecture 2
Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No
More informationKey concepts in Density Functional Theory (II)
Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université
More informationCondensed matter physics FKA091
Condensed matter physics FKA091 Ermin Malic Department of Physics Chalmers University of Technology Henrik Johannesson Department of Physics University of Gothenburg Teaching assistants: Roland Jago &
More informationSolid State Physics. Lecture 10 Band Theory. Professor Stephen Sweeney
Solid State Physics Lecture 10 Band Theory Professor Stephen Sweeney Advanced Technology Institute and Department of Physics University of Surrey, Guildford, GU2 7XH, UK s.sweeney@surrey.ac.uk Recap from
More informationBasics of density-functional theory and fast guide to actual calculations Matthias Scheffler
Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler http://www.fhi-berlin.mpg.de/th/th.html I. From the many-particle problem to the Kohn-Sham functional II. From
More informationOptical Properties of Semiconductors. Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India
Optical Properties of Semiconductors 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India http://folk.uio.no/ravi/semi2013 Light Matter Interaction Response to external electric
More informationFirst-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design
First-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design Peter Kratzer 5/2/2018 Peter Kratzer Abeokuta School 5/2/2018 1 / 34 Outline
More informationStudy of Carbon Nanomaterials Based on Density Functional Theory. Mohammad Shafiul Alam
Study of Carbon Nanomaterials Based on Density Functional Theory Mohammad Shafiul Alam July 2013 Dissertation Study of Carbon Nanomaterials Based on Density Functional Theory Graduate School of Natural
More informationFYS-6306 QUANTUM THEORY OF MOLECULES AND NANOSTRUCTURES
i FYS-6306 QUANTUM THEORY OF MOLECULES AND NANOSTRUCTURES Credit units: 6 ECTS Lectures: 48 h Tapio Rantala, prof. Tue 10 12 SC203 SG219 8 10 SG312 FirstName.LastName@tut.fi http://www.tut.fi/~trantala/opetus/
More informationBasics of DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Basics of DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Basics of DFT Lausanne12 1
More informationA FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY
A FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY JOHN P. PERDEW PHYSICS & CHEMISTRY, TEMPLE UNIVERSITY PHILADELPHIA, PENNSYLVANIA, USA SUPPORTED BY THE U.S. DEPARTMENT OF ENERGY, EFRC
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationDensity functional theory in the solid state
Density functional theory in the solid state Ari P Seitsonen IMPMC, CNRS & Universités 6 et 7 Paris, IPGP Department of Applied Physics, Helsinki University of Technology Physikalisch-Chemisches Institut
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationLecture 7: Extrinsic semiconductors - Fermi level
Lecture 7: Extrinsic semiconductors - Fermi level Contents 1 Dopant materials 1 2 E F in extrinsic semiconductors 5 3 Temperature dependence of carrier concentration 6 3.1 Low temperature regime (T < T
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationIntroduction to spin and spin-dependent phenomenon
Introduction to spin and spin-dependent phenomenon Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. May 16th, 2007
More informationThe calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationELECTRONIC STRUCTURE IN THE ORGANIC-BASED MAGNET LINEAR COMBINATION OF ATOMIC ORBITALS METHOD WITH HUBBARD U. A THESIS IN Physics
ELECTRONIC STRUCTURE IN THE ORGANIC-BASED MAGNET [Fe II (TCNE)(NCMe)2][Fe III Cl4] USING THE ORTHOGONALIZED LINEAR COMBINATION OF ATOMIC ORBITALS METHOD WITH HUBBARD U A THESIS IN Physics Presented to
More informationarxiv:cond-mat/ v1 [cond-mat.other] 24 Mar 2005
arxiv:cond-mat/0503601v1 [cond-mat.other] 4 Mar 005 Abstract Vortices as Yukawa bosons in the two-dimensional Harmonic confinement K. K. Rajagopal A system of N vortices in a two-dimensional harmonic trap
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 informationIntroduction to Sources: Radiative Processes and Population Inversion in Atoms, Molecules, and Semiconductors Atoms and Molecules
OPTI 500 DEF, Spring 2012, Lecture 2 Introduction to Sources: Radiative Processes and Population Inversion in Atoms, Molecules, and Semiconductors Atoms and Molecules Energy Levels Every atom or molecule
More informationLecture 1. OUTLINE Basic Semiconductor Physics. Reading: Chapter 2.1. Semiconductors Intrinsic (undoped) silicon Doping Carrier concentrations
Lecture 1 OUTLINE Basic Semiconductor Physics Semiconductors Intrinsic (undoped) silicon Doping Carrier concentrations Reading: Chapter 2.1 EE105 Fall 2007 Lecture 1, Slide 1 What is a Semiconductor? Low
More informationThe potential of Potential Functional Theory
The potential of Potential Functional Theory IPAM DFT School 2016 Attila Cangi August, 22 2016 Max Planck Institute of Microstructure Physics, Germany P. Elliott, AC, S. Pittalis, E.K.U. Gross, K. Burke,
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
More informationSpins and spin-orbit coupling in semiconductors, metals, and nanostructures
B. Halperin Spin lecture 1 Spins and spin-orbit coupling in semiconductors, metals, and nanostructures Behavior of non-equilibrium spin populations. Spin relaxation and spin transport. How does one produce
More informationAdvanced Solid State Theory SS Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt
Advanced Solid State Theory SS 2010 Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt i 0. Literatur R. M. Martin, Electronic Structure: Basic Theory and
More informationThe Semiconductor in Equilibrium
Lecture 6 Semiconductor physics IV The Semiconductor in Equilibrium Equilibrium, or thermal equilibrium No external forces such as voltages, electric fields. Magnetic fields, or temperature gradients are
More informationOptical Properties of Solid from DFT
Optical Properties of Solid from DFT 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology, University of Oslo, Norway http://folk.uio.no/ravi/cmt15
More informationSemiconductor Physics. Lecture 3
Semiconductor Physics Lecture 3 Intrinsic carrier density Intrinsic carrier density Law of mass action Valid also if we add an impurity which either donates extra electrons or holes the number of carriers
More informationABC of ground-state DFT
ABC of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) ABC of ground-state
More informationSemiconductor Module
Semiconductor Module Optics Seminar July 18, 2018 Yosuke Mizuyama, Ph.D. COMSOL, Inc. The COMSOL Product Suite Governing Equations Semiconductor Schrödinger Equation Semiconductor Optoelectronics, FD Semiconductor
More informationCHEM6085: Density Functional Theory
Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron
More informationSpin Superfluidity and Graphene in a Strong Magnetic Field
Spin Superfluidity and Graphene in a Strong Magnetic Field by B. I. Halperin Nano-QT 2016 Kyiv October 11, 2016 Based on work with So Takei (CUNY), Yaroslav Tserkovnyak (UCLA), and Amir Yacoby (Harvard)
More informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More informationCorrelated 2D Electron Aspects of the Quantum Hall Effect
Correlated 2D Electron Aspects of the Quantum Hall Effect Magnetic field spectrum of the correlated 2D electron system: Electron interactions lead to a range of manifestations 10? = 4? = 2 Resistance (arb.
More informationψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.
1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product
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 informationECE 340 Lecture 39 : MOS Capacitor II
ECE 340 Lecture 39 : MOS Capacitor II Class Outline: Effects of Real Surfaces Threshold Voltage MOS Capacitance-Voltage Analysis Things you should know when you leave Key Questions What are the effects
More informationLecture 18: Semiconductors - continued (Kittel Ch. 8)
Lecture 18: Semiconductors - continued (Kittel Ch. 8) + a - Donors and acceptors J U,e e J q,e Transport of charge and energy h E J q,e J U,h Physics 460 F 2006 Lect 18 1 Outline More on concentrations
More informationSurfaces, Interfaces, and Layered Devices
Surfaces, Interfaces, and Layered Devices Building blocks for nanodevices! W. Pauli: God made solids, but surfaces were the work of Devil. Surfaces and Interfaces 1 Interface between a crystal and vacuum
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics University of York Lecturer: Rex Godby Notes by Victor Naden Robinson Lecture 1: TDSE Lecture 2: TDSE Lecture 3: FMG Lecture 4: FMG Lecture 5: Ehrenfest s Theorem and the Classical
More informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
More informationPrinciples of Quantum Mechanics
Principles of Quantum Mechanics - indistinguishability of particles: bosons & fermions bosons: total wavefunction is symmetric upon interchange of particle coordinates (space,spin) fermions: total wavefuncftion
More informationFREQUENTLY ASKED QUESTIONS February 21, 2017
FREQUENTLY ASKED QUESTIONS February 21, 2017 Content Questions How do you place a single arsenic atom with the ratio 1 in 100 million? Sounds difficult to get evenly spread throughout. Yes, techniques
More informationLecture 20: Semiconductor Structures Kittel Ch 17, p , extra material in the class notes
Lecture 20: Semiconductor Structures Kittel Ch 17, p 494-503, 507-511 + extra material in the class notes MOS Structure Layer Structure metal Oxide insulator Semiconductor Semiconductor Large-gap Semiconductor
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More information