Lecture 16: DFT for Metallic Systems
|
|
- Derrick Bailey
- 5 years ago
- Views:
Transcription
1 The Nuts and Bolts of First-Principles Simulation Lecture 16: DFT for Metallic Systems Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network
2 Overview of talk What is a metal? Problems with metals Finite temperature DFT Density mixing Ensemble DFT Conclusions Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 2
3 1. What is a metal? For our purposes a metal is any system with unoccupied states very close to the Fermi level This means that several bands may cross over near the Fermi surface Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 3
4 2. Problems with metals Band crossings at Fermi level Charge sloshing These manifest themselves in several different ways Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 4
5 Review of Orthogonalisation We use a Gram-Schmidt method to orthogonalise each band of the search direction to all bands of the wavefunction We do not orthogonalise with respect to bands above the valence band For most systems these bands are much higher in energy, so the energy minimisation will remove them from the top band Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 5
6 Orthogonalisation Problems For metals there are bands very close in energy to the valence band, so the energy minimisation will take a long time to remove them from our trial wavefunction By including higher, unoccupied bands we can ensure they are orthogonal to the lower bands Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 6
7 Unoccupied Bands We want to include some unoccupied bands in our calculation The cost of our calculation scales quadratically with bands, so we need to keep these to a minimum Unfortunately we cannot determine beforehand which bands are unoccupied We must run the calculation, and then check that the top band is unoccupied Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 7
8 Band Occupancies Since we re including bands above the Fermi level, we need to know which bands are important when constructing the density Introduce the concept of band occupancies, {f i }, which are 0 if the band is unoccupied and 1 if occupied n(r) = f i Ψ i 2 i Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 8
9 Assigning Occupancies Our energy minimisation algorithm gives us wavefunctions which are orthogonal mixtures of the Kohn-Sham eigenstates In order to determine which states are occupied or not, we need the true eigenstates of the Kohn-Sham Hamiltonian We diagonalise the Hamiltonian in the subspace of the bands to get the true eigenstates and eigenenergies Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 9
10 The new line search Since the density is now a function of the occupancies, as well as the wavefunctions, we need to modify our search algorithm. Ψ i + λφ i As before, we take a trial step We now recalculate the occupancies as well as the density Fit parabola and move to parabolic minimum Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 10
11 Unfortunately the occupancies are discontinuous and so even very small steps can completely change which bands are important for the density. The result is that this algorithm is often unstable. Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 11
12 3. Finite Temperature DFT We know that at nonzero temperatures the occupancies are no longer discontinuous Mermin extended the use of LDA DFT to systems with finite temperatures, in which case the bands are smeared in energy We now need to minimise the free energy of the system F = E σ n S Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 12 ε σ n E F
13 Partial Occupancies We can use the same idea for our metallic calculations to improve the conditioning We smear the bands in energy so that the occupancies become continuous There is an additional entropic contribution to the energy which must be calculated. In general we have to approximate this term Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 13
14 A Smooth Operator Because we re only using the smearing as a way to improve the conditioning of our problem, we re not restricted to physical smearing schemes We can use any smooth operator to smear our bands, and if we can accurately calculate the entropic contribution then we can always recover the zero-temperature result Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 14
15 Smearing Schemes Gaussian Fermi- Dirac Cold Smearing (Methfessel-Paxton) Methfessel and Paxton expanded the deltafunction in Hermite polynomials. The entropic contribution due to this smearing can be calculated accurately, allowing good zero-temperature energies to be obtained using large smearing widths. Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 15
16 Charge Sloshing Restrict our attention to the Hartree potential (usually the most important) Response of the system to a perturbation is given by the dielectric matrix, J: J = 1 χu Where χ is the susceptibility and U is given by: G U G = δ GG 4 π G e 2 2 Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 16
17 Thus we have: V H n(g) G 2 Sloshing instabilities arise for large simulation cells, where G is small For metals, degenerate states at the Fermi level can also lead to instabilities, since macroscopic changes in the density can occur for little change in energy Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 17
18 4. Density Mixing Wavefunctions found for fixed Hamiltonian Trial density calculated from wavefunctions The density is then mixed with previous densities, and the mixed density used to construct the new Hamiltonian n ( i+ 1) ( r ) = n ( i) ( r ) +α R[ n ( i ) ( r)] Where R is the density residual, defined as: R [ ( )] ( i ) 2 ( i ) n r = f Ψ n ( r ) n Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 18 n n
19 We can use a more general mixing scheme: n i ( i + 1) ( i+ 1) ( r ) = n ( r) + α j j = 1 R [ ] ( j ) n Density evolution is now decoupled from that of the wavefunctions The wavefunction and density searches can be preconditioned separately The scheme is no longer variational Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 19
20 Density Evolution By mixing the density with past densities, we damp out charge oscillations If the damping is small, charge sloshing may still occur If the damping is large, the system will not converge rapidly to the groundstate Various mixing schemes exist which attempt to optimally mix the density at each step Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 20
21 5. Ensemble DFT The instabilities arise because we do not properly account for the self-consistent variation of the density with our trial step We can avoid these instabilities if we take the trial step not for a fixed Hamiltonian, but for fixed occupancies Once we have found the optimum wavefunctions we minimise the energy again, this time with respect to the occupancies Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 21
22 The EDFT Algorithm We have split our minimisation problem into a search over the wavefunctions followed by a search over the occupancies We need to perform an occupancy search every time the wavefunctions are updated Each time we update the occupancies we must reapply the Hamiltonian Scheme is fully variational Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 22
23 6. Conclusions For metallic systems we must include unoccupied bands The bands must be smeared with a typical smearing width O(0.1)eV Sloshing instabilities can arise, and must be quenched using density mixing, or circumvented using ensemble DFT Metallic calculations are more expensive than for semiconductors, but good convergence can be achieved. Nuts and Bolts 2001 Lecture 16: DFT for Metallic Systems 23
3: 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 informationPractical Density Functional Theory
Practical Density Functional Theory I. Abdolhosseini Sarsari 1, 1 Department of Physics Isfahan University of Technology April 29, 2015 Collaborators at IUT : M. Alaei and S. J. Hashemifar April 29, 2015
More informationMaking Electronic Structure Theory work
Making Electronic Structure Theory work Ralf Gehrke Fritz-Haber-Institut der Max-Planck-Gesellschaft Berlin, 23 th June 2009 Hands-on Tutorial on Ab Initio Molecular Simulations: Toward a First-Principles
More informationElectronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory
Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.
More informationBehind the "exciting" curtain: The (L)APW+lo method
Behind the "exciting" curtain: The (L)APW+lo method Aug 7, 2016 Andris Gulans Humboldt-Universität zu Berlin Kohn-Sham equation Potential due to nuclei Exchange-correlation potential Potential due to electron
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
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 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 informationDFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics,
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 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 informationLS coupling. 2 2 n + H s o + H h f + H B. (1) 2m
LS coupling 1 The big picture We start from the Hamiltonian of an atomic system: H = [ ] 2 2 n Ze2 1 + 1 e 2 1 + H s o + H h f + H B. (1) 2m n e 4πɛ 0 r n 2 4πɛ 0 r nm n,m Here n runs pver the electrons,
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
More information6: Plane waves, unit cells, k- points and all that
The Nuts and Bolts of First-Principles Simulation 6: Plane waves, unit cells, k- points and all that Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview
More informationSection 9 Variational Method. Page 492
Section 9 Variational Method Page 492 Page 493 Lecture 27: The Variational Method Date Given: 2008/12/03 Date Revised: 2008/12/03 Derivation Section 9.1 Variational Method: Derivation Page 494 Motivation
More informationGeometry Optimisation
Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background
More informationBand Structure Calculations; Electronic and Optical Properties
; Electronic and Optical Properties Stewart Clark University of Outline Introduction to band structures Calculating band structures using Castep Calculating optical properties Examples results Some applications
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 information10: Testing Testing. Basic procedure to validate calculations
The Nuts and Bolts of First-Principles Simulation 10: Testing Testing. Basic procedure to validate calculations Durham, 6th-13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network
More informationPage 404. Lecture 22: Simple Harmonic Oscillator: Energy Basis Date Given: 2008/11/19 Date Revised: 2008/11/19
Page 404 Lecture : Simple Harmonic Oscillator: Energy Basis Date Given: 008/11/19 Date Revised: 008/11/19 Coordinate Basis Section 6. The One-Dimensional Simple Harmonic Oscillator: Coordinate Basis Page
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationDFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
More informationMODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
More informationIntroduction to multiconfigurational quantum chemistry. 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. Notations
More informationPerturbation Theory 1
Perturbation Theory 1 1 Expansion of Complete System Let s take a look of an expansion for the function in terms of the complete system : (1) In general, this expansion is possible for any complete set.
More informationHandbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield
Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
More informationN independent electrons in a volume V (assuming periodic boundary conditions) I] The system 3 V = ( ) k ( ) i k k k 1/2
Lecture #6. Understanding the properties of metals: the free electron model and the role of Pauli s exclusion principle.. Counting the states in the E model.. ermi energy, and momentum. 4. DOS 5. ermi-dirac
More informationReferences. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
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 information10 Time-Independent Perturbation Theory
S.K. Saiin Oct. 6, 009 Lecture 0 0 Time-Independent Perturbation Theory Content: Non-degenerate case. Degenerate case. Only a few quantum mechanical problems can be solved exactly. However, if the system
More informationPreconditioned Eigenvalue Solvers for electronic structure calculations. Andrew V. Knyazev. Householder Symposium XVI May 26, 2005
1 Preconditioned Eigenvalue Solvers for electronic structure calculations Andrew V. Knyazev Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver Householder
More informationBrief review of Quantum Mechanics (QM)
Brief review of Quantum Mechanics (QM) Note: This is a collection of several formulae and facts that we will use throughout the course. It is by no means a complete discussion of QM, nor will I attempt
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
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 informationLecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 6 Fermion Pairing WS2010/11: Introduction to Nuclear and Particle Physics Experimental indications for Cooper-Pairing Solid state physics: Pairing of electrons near the Fermi surface with antiparallel
More information11 Perturbation Theory
S.K. Saikin Oct. 8, 009 11 Perturbation Theory Content: Variational Principle. Time-Dependent Perturbation Theory. 11.1 Variational Principle Lecture 11 If we need to compute the ground state energy of
More informationElectronic energy optimisation in ONETEP
Electronic energy optimisation in ONETEP Chris-Kriton Skylaris cks@soton.ac.uk 1 Outline 1. Kohn-Sham calculations Direct energy minimisation versus density mixing 2. ONETEP scheme: optimise both the density
More informationDFT: Exchange-Correlation
DFT: Exchange-Correlation Local functionals, exact exchange and other post-dft methods Paul Tulip Centre for Materials Physics Department of Physics University of Durham Outline Introduction What is exchange
More informationIntroduction to Density Functional Theory in a Plane Wave + Pseudopotential Framework
Introduction to Density Functional Theory in a Plane Wave + Pseudopotential Framework Shobhana Narasimhan Theoretical Sciences Unit JNCASR, Bangalore shobhana@jncasr.ac.in School- From the chemical bond
More informationS. Bellucci, A. Sindona, D. Mencarelli, L. Pierantoni Electrical conductivity of graphene: a timedependent density functional theory study
S. Bellucci, A. Sindona, D. Mencarelli, L. Pierantoni Electrical conductivity of graphene: a timedependent density functional theory study INFN Laboratori Nazionali Frascati (LNF), Italy Univ. Calabria,
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
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 informationMany-Body Perturbation Theory. Lucia Reining, Fabien Bruneval
, Fabien Bruneval Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Belfast, 27.6.2007 Outline 1 Reminder 2 Perturbation Theory
More information1. For the case of the harmonic oscillator, the potential energy is quadratic and hence the total Hamiltonian looks like: d 2 H = h2
15 Harmonic Oscillator 1. For the case of the harmonic oscillator, the potential energy is quadratic and hence the total Hamiltonian looks like: d 2 H = h2 2mdx + 1 2 2 kx2 (15.1) where k is the force
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 informationCHAPTER 8. EXPLORING THE EFFECTS OF STRUCTURAL INSTABILITIES AND OF TRIPLET INSTABILITY ON THE M-I PHASE TRANSITION IN (EDO-TTF)2PF6.
CHAPTER 8. EXPLORING THE EFFECTS OF STRUCTURAL INSTABILITIES AND OF TRIPLET INSTABILITY ON THE M-I PHASE TRANSITION IN (EDO-TTF)2PF6. The thermal metal-insulator phase transition in (EDO-TTF)2PF6 is attributed
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 informationSpin-Orbit Interactions in Semiconductor Nanostructures
Spin-Orbit Interactions in Semiconductor Nanostructures Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, U.S.A. http://www.physics.udel.edu/~bnikolic Spin-Orbit Hamiltonians
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
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 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 informationAdiabatic connection for near degenerate excited states
PHYSICAL REVIEW A 69, 052510 (2004) Adiabatic connection for near degenerate excited states Fan Zhang Department of Physics and Astronomy, Rutgers University, 136 Frelinghuysen Road, Piscataway, New Jersey
More informationThe Gutzwiller Density Functional Theory
The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller
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 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 informationLecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor
Lecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor It turns out that the boundary condition of the wavefunction going to zero at infinity is sufficient to quantize the value of energy that
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 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 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 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 informationElectronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier Pôle de chimie théorique Département de chimie de l ENS CNRS Ecole normale supérieure UPMC Formation ModPhyChem Lyon,
More informationThe free electron gas: Sommerfeld model
The free electron gas: Sommerfeld model Valence electrons: delocalized gas of free electrons (selectrons of alkalies) which is not perturbed by the ion cores (ionic radius of Na + equals 0.97 Å, internuclear
More informationAdvanced Spectroscopy. Dr. P. Hunt Rm 167 (Chemistry) web-site:
Advanced Spectroscopy Dr. P. Hunt p.hunt@imperial.ac.uk Rm 167 (Chemistry) web-site: http://www.ch.ic.ac.uk/hunt Maths! Coordinate transformations rotations! example 18.1 p501 whole chapter on Matrices
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 informationMassachusetts Institute of Technology Physics Department
Massachusetts Institute of Technology Physics Department Physics 8.32 Fall 2006 Quantum Theory I October 9, 2006 Assignment 6 Due October 20, 2006 Announcements There will be a makeup lecture on Friday,
More informationHandbook of Computational Quantum Chemistry
Handbook of Computational Quantum Chemistry David B. Cook Dept. of Chemistry University of Sheffield DOVER PUBLICATIONS, INC. Mineola, New York F Contents 1 Mechanics and molecules 1 1.1 1.2 1.3 1.4 1.5
More informationLecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas
Lecture on First-principles Computation (2): The Hartree-Fock theory and the homogeneous electron gas 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2016.9.18 Recall:
More informationSIMPLE QUANTUM SYSTEMS
SIMPLE QUANTUM SYSTEMS Chapters 14, 18 "ceiiinosssttuu" (anagram in Latin which Hooke published in 1676 in his "Description of Helioscopes") and deciphered as "ut tensio sic vis" (elongation of any spring
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 informationGW quasiparticle energies
Chapter 4 GW quasiparticle energies Density functional theory provides a good description of ground state properties by mapping the problem of interacting electrons onto a KS system of independent particles
More informationMetal Conquest. Lianheng Tong. 2 nd March 2011
Metal Conquest Lianheng Tong nd March 011 Abstract This report describes the work done in the one year Distributed Computational Science and Engineering (dcse) project aimed to develop an ab initio Density
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 informationMany electron effects in semiconductor quantum dots
Bull. Mater. Sci., Vol. 26, No. 1, January 2003, pp. 63 67. Indian Academy of Sciences. Many electron effects in semiconductor quantum dots K PANDEY*, MANOJ K HABOLA, V ANJAN and VIJAY A SINGH Physics
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
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 informationQUALIFYING EXAMINATION, Part 2. Solutions. Problem 1: Quantum Mechanics I
QUALIFYING EXAMINATION, Part Solutions Problem 1: Quantum Mechanics I (a) We may decompose the Hamiltonian into two parts: H = H 1 + H, ( ) where H j = 1 m p j + 1 mω x j = ω a j a j + 1/ with eigenenergies
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 informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
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 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 informationNext topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization"
Next topic: Quantum Field Theories for Quantum Many-Particle Systems; or "Second Quantization" Outline 1 Bosons and Fermions 2 N-particle wave functions ("first quantization" 3 The method of quantized
More informationIntroduction to Parallelism in CASTEP
to ism in CASTEP Stewart Clark Band University of Durham 21 September 2012 Solve for all the bands/electrons (Band-) Band CASTEP solves the Kohn-Sham equations for electrons in a periodic array of nuclei:
More informationnanohub.org learning module: Prelab lecture on bonding and band structure in Si
nanohub.org learning module: Prelab lecture on bonding and band structure in Si Ravi Vedula, Janam Javerhi, Alejandro Strachan Center for Predictive Materials Modeling and Simulation, School of Materials
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 informationLecture. Ref. Ihn Ch. 3, Yu&Cardona Ch. 2
Lecture Review of quantum mechanics, statistical physics, and solid state Band structure of materials Semiconductor band structure Semiconductor nanostructures Ref. Ihn Ch. 3, Yu&Cardona Ch. 2 Reminder
More informationIntroduction to Quantum Mechanics PVK - Solutions. Nicolas Lanzetti
Introduction to Quantum Mechanics PVK - Solutions Nicolas Lanzetti lnicolas@student.ethz.ch 1 Contents 1 The Wave Function and the Schrödinger Equation 3 1.1 Quick Checks......................................
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 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 informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
More informationPhysics 127c: Statistical Mechanics. Weakly Interacting Fermi Gas. The Electron Gas
Physics 7c: Statistical Mechanics Wealy Interacting Fermi Gas Unlie the Boson case, there is usually no ualitative change in behavior going from the noninteracting to the wealy interacting Fermi gas for
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 informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More informationNext generation extended Lagrangian first principles molecular dynamics
Next generation extended Lagrangian first principles molecular dynamics Anders M. N. Niklasson Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (Dated: June 1, 017) LA-UR-17-3007
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 informationIn-class exercises. Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 8 Exercises due Mon March 19 Last correction at March 5, 2018, 8:48 am c 2017, James Sethna,
More informationLinear-response excitations. Silvana Botti
from finite to extended systems 1 LSI, CNRS-CEA-École Polytechnique, Palaiseau, France 2 LPMCN, CNRS-Université Lyon 1, France 3 European Theoretical Spectroscopy Facility September 3, 2008 Benasque, TDDFT
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 information