Transcript Lecture 1
Quantum Theory of Solids Mervyn Roy (S6) www2.le.ac.uk/departments/physics/people/mervynroy PA4311 Quantum Theory of Solids Course Outline 1. Introduction and background 2. The many-electron wavefunction - Introduction to quantum chemistry (Hartree, HF, and CI methods) 3. Introduction to density functional theory (DFT) - Periodic solids, plane waves and pseudopotentials 4. Linear combination of atomic orbitals 5. Effective mass theory 6. ABINIT computer workshop (LDA DFT for periodic solids) Assessment: 70% final exam 30% coursework – mini ‘project’ report for ABINIT calculation (Set problems are purely formative) PA4311 Quantum Theory of Solids Resources Lecture notes www2.le.ac.uk/departments/physics/people/mervynroy Abinit website - including comprehensive help and tutorials • www.abinit.org Books • Electronic Structure, RM Martin • Solid State Physics, Ashcroft and Mermin • Solid State Physics, Hook and Hall Plus many other relevant text books and online references – see the library Prerequisites • 3210 Quantum Mechanics • 2230 Condensed Matter Physics • Also – 214 Fourier Series, 372 Fourier Transforms, etc. PA4311 Quantum Theory of Solids The problem electron-ion interaction electron KE 𝐻Ψ = 𝐸Ψ 𝐻 −ℏ2 = 2𝑚 ℏ2 − ion 2𝑚 KE 𝐼 𝛻𝑖2 − 𝑖 𝐼 𝑖 1 2 𝛻𝐼 + 2 𝐼≠𝐽 𝐼 𝑍𝐼 𝑒 2 1 + 4𝜋𝜖0 𝒓𝑖 − 𝑹𝐼 2 𝑍𝐼 𝑍𝐽 𝑒 2 4𝜋𝜖0 |𝑹𝐼 − 𝑹𝐽 | electronelectron interaction 𝑒2 𝑖≠𝑗 4𝜋𝜖0 𝒓𝑖 − 𝒓𝑗 ion-ion interaction 𝑁 > 1023 electrons and ions Ψ is a function of 𝑁 electron co-ordinates, 𝒓𝑖 , (and spins), and 𝑁 ion co-ordinates, 𝑹𝐼 (and spins) But, mI ≫ 𝑚 so ion KE term is small PA4311 Quantum Theory of Solids Timescales From CA Ullrich, TimeDependent DensityFunction Theory, Oxford University Press (2012) PA4311 Quantum Theory of Solids The problem Born-Oppenheimer approximation - electrons react instantaneously to changes in nuclear positions −ℏ2 𝐻= 2𝑚 𝛻𝑖2 − 𝑖 𝑖 𝐼 𝑍𝐼 𝑒 2 1 + 4𝜋𝜖0 𝒓𝑖 − 𝑹𝐼 2 𝑒2 𝑖≠𝑗 4𝜋𝜖0 𝒓𝑖 − 𝒓𝑗 Or, in atomic units, 1 𝐻=− 2 𝛻𝑖2 𝑖 − 𝑖 1 𝑣 𝒓𝑖 + 2 1 𝑖≠𝑗 𝒓𝑖 − 𝒓𝑗 + 𝐸𝐼𝐼 But, still have N > 1023 electrons Ψ is a function of 𝑁 electron co-ordinates, 𝒓𝑖 , (and spins). Need to develop some approximations! PA4311 Quantum Theory of Solids + 𝐸𝐼𝐼 Constant depends on ion positions Why bother? The modern world is build upon our understanding of the electronic properties of solids Solid state (nano) physics, materials physics, space technology etc. etc. Spectroscopy e.g. astrophysics, Earth observation science – ExoMol line lists (TDDFT) Plasma physics… PA4311 Quantum Theory of Solids Highest cited papers in Physical Review suite of journals (2014) Citations Generalized Gradient Approximation Made Simple, JP Perdew, K Burke, and M Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996) 25 083 Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density, C Lee, W Yang, and RG Parr, Phys. Rev. B 37, 785 (1988) 24 292 Self-Consistent Equations Including Exchange and Correlation Effects, W Kohn and LJ Sham, Phys. Rev. 140, A1133 (1965) 18 399 Inhomogeneous Electron Gas, P Hohenberg and W Kohn, Phys. Rev. 136, B864 (1964) 15 629 Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, G Kresse and J Furthmüller, Phys. Rev. B 54, 11169 (1996) 14 495 Density-functional exchange-energy approximation with correct asymptotic behavior, AD Becke, Phys. Rev. A 38, 3098 (1988) 14 191 Special points for Brillouin-zone integrations, HJ Monkhorst, JD Pack, Phys. Rev. B 13, 5188 (1976) 12 938 From ultrasoft pseudopotentials to the projector augmented-wave method, G Kresse and D Joubert, Phys. Rev. B 59, 1758 (1999) 10 351 Helical microtubules of graphitic carbon, S. Ijima, Nature 354 , 56 (1991) (24 225) Electric field effect in atomically thin carbon films, K.S. Novoselov, A.K. Geim, et al. Science 306 (2004) (15 139) Highest cited paper on ADS astrophysics (2014) Maps of Dust Infrared Emission for Use in Estimation of Reddening and Cosmic Microwave Background Radiation Foregrounds, D.J. Schlegel et al, APJ 500, 525 (1998) 8920 Comparison with other areas - data from ADS Black hole 6765 Particle physics 4164 Metamaterials 2289 Aurora 1785 Gamma Ray Burst 1513 Density-Functional Thermochemistry 3. The Role of Exact Exchange, AD Becke, J. Chemical Physics 98, 5648 (1993) 46 280 Chemistry (data from WOK) PA4311 Quantum Theory of Solids 3210 revision (see Rae) Schrödinger equation - 𝐻𝜓(𝒓) = 𝐸𝜓(𝒓) Probability density - |𝜓(𝒓)|2 𝜓 ∗ 𝑂𝜓𝑑𝒓 Expectation values - 𝑂 = 𝑉 𝜙𝐻𝜙 Variational principle - 𝐸 = ≥ 𝐸0 𝜙𝜙 If 𝜓 = 𝑐𝑛 𝜙𝑛 , then 𝑛 PA4311 Quantum Theory of Solids 𝜕𝐸 ∗ =0= 𝜕𝑐𝑚 𝐻𝑚𝑛 − 𝐸𝛿𝑚𝑛 𝑐𝑛 𝑛 Question 1.1 If 𝜓 = 𝑛 𝑐𝑛 𝜙𝑛 is normalised and the 𝜙𝑛 are orthonormal, show that 𝑛 𝑐𝑛 2 = 1. ∗ 𝑐 𝜙 𝐻 𝜙 subject If we wish to minimise 𝐸 = 𝑚 𝑛 𝑐𝑚 𝑛 𝑚 𝑛 to the constraint that the 𝜙𝑛 are normalised, show that the appropriate Lagrange multiplier, 𝜆 = 𝐸. PA4311 Quantum Theory of Solids Full variation and functionals 𝐸 is a functional of 𝜓 ∗ and 𝜓 𝐸 𝜓∗, 𝜓 = 𝜓 𝐻 𝜓 − 𝐸 𝜓 𝜓 − 1 𝛿𝐸 𝜓 ∗ , 𝜓 = 𝛿𝜓 (𝐻 − 𝐸) 𝜓 = ∫ 𝛿𝜓 ∗ 𝐻 − 𝐸 𝜓𝑑𝒓 = 0 𝐻−𝐸 𝜓 =0 A functional maps a function onto a value, 𝐹 𝑓 𝑥 𝐹𝑓 𝑥 = 𝑥2 = 𝑎. For example, 𝑤 𝑥 𝑓 𝑥 𝑑𝑥 𝑥1 𝛿𝐹 𝑓 = 𝐹 𝑓 + 𝛿𝑓 − 𝐹 𝑓 = 𝑥2 𝑤 𝑥 𝛿𝑓 𝑥 𝑑𝑥 𝑥1 See RM Martin App. A, or e.g. GC Evans, Functionals and their applications, Dover, New York, 1964 PA4311 Quantum Theory of Solids The N-electron wavefunction The 𝑁-electron wavefunction depends on N spatial coordinates (and spins) Ψ(𝒓1 , 𝒓2 , … , 𝒓N ) Electrons are indistinguishable: Ψ 𝒓1 , 𝒓2 , … , 𝒓N 2 = Ψ 𝒓2 , 𝒓1 , … , 𝒓N 2 Fermions are anti-symmetric: Ψ 𝒓1 , 𝒓2 , … , 𝒓N = −Ψ 𝒓2 , 𝒓1 , … , 𝒓N - they obey the Pauli exclusion principle See Tipler (4th Ed Sec. 36.6 on ‘The Schrödinger equation for 2 identical particles’) Expectation values 𝑂 = Ψ𝑂Ψ = 𝑉 Ψ ∗ 𝒓1 , 𝒓2 , … , 𝒓𝑁 𝑂Ψ(𝒓1 , 𝒓2 , … , 𝒓𝑁 ) 𝑑𝒓𝟏 𝑑𝒓2 … 𝑑𝒓𝑁 PA4311 Quantum Theory of Solids