220 likes | 389 Views
Dive into the world of Density Functional Theory (DFT) with this comprehensive guide. Learn the theory behind DFT, practical applications, and the basics of solid-state physics. Discover the complexities of DFT energy functionals, exchange-correlation potentials, and the role of electronic density in determining ground state energy. Explore crystal structures, band structures, reciprocal space, and the use of plane waves in solid-state calculations. With this crash course, you'll gain valuable insights into computational chemistry and its practical implementations.
E N D
Useful Material • Books • A chemist’s guide to density-functional theoryWolfram Koch and Max C. Holthausen (second edition, Wiley) • The theory of the cohesive energies of solidsG. P. Srivastava and D. WeaireAdvances in Physics 36 (1987) 463-517 • Gulliver among the atomsMike Gillan, New Scientist 138 (1993) 34 • Web • www.nobel.se/chemistry/laureates/1998/ • www.abinit.org Version 4.2.3 compiled for windows, install and good tutorial
Outline: Part 1, The Framework of DFT • DFT: the theory • Schroedinger’s equation • Hohenberg-Kohn Theorem • Kohn-Sham Theorem • Simplifying Schroedinger’s • LDA, GGA • Elements of Solid State Physics • Reciprocal space • Band structure • Plane waves • And then ? • Forces (Hellmann-Feynman theorem) • E.O., M.D., M.C. …
Outline: Part2Using DFT • Practical Issues • Input File(s) • Output files • Configuration • K-points mesh • Pseudopotentials • Control Parameters • LDA/GGA • ‘Diagonalisation’ • Applications • Isolated molecule • Bulk • Surface
The Basic Problem Dangerously classical representation Cores Electrons
Schroedinger’s Equation Wave function Potential Energy Kinetic Energy Coulombic interaction External Fields Energy levels Hamiltonian operator Very Complex many body Problem !! (Because everything interacts)
First approximations • Adiabatic (or Born-Openheimer) • Electrons are much lighter, and faster • Decoupling in the wave function • Nuclei are treated classically • They go in the external potential
H.K. Theorem The ground state is unequivocally defined by the electronic density Universalfunctional • Functional ?? Function of a function • No more wave functions here • But still too complex
K.S. Formulation Use an auxiliary system • Non interacting electrons • Same Density • => Back to wave functions, but simpler this time(a lot more though) N K.S. equations (ONE particle in a box really) (KS1) Exchange correlation potential (KS2) (KS3)
Self consistent loop Initial density From density, work out Effective potential Solve the independents K.S. =>wave functions Deduce new density from w.f. New density ‘=‘ input density ?? NO YES Finita la musica
DFT energy functional • Exchange correlation funtional • Contains: • Exchange • Correlation • Interacting part of K.E. Electrons are fermions (antisymmetric wave function)
Exchange correlation functional At this stage, the only thing we need is: Still a functional (way too many variables) • #1 approximation, Local Density Approximation: • Homogeneous electron gas • Functional becomes function !! (see KS3) • Very good parameterisation for LDA Generalised Gradient Approximation: GGA
DFT: Summary • The ground state energy depends only on the electronic density (H.K.) • One can formally replace the SE for the system by a set of SE for non-interacting electrons (K.S.) • Everything hard is dumped into Exc • Simplistic approximations of Exc work ! LDA or GGA
And now, for something completely different:A little bit of Solid State Physics Crystal structure Periodicity
Reciprocal space (Inverting effect) sin(k.r) Reciprocal Space bi Real Space ai Brillouin Zone k-vector (or k-point) See X-Ray diffraction for instance Also, Fourier transform and Bloch theorem
Band structure E Energy levels (eigenvalues of SE) Crystal Molecule
The k-point mesh Corresponds to a supercell 36 time bigger than the primitive cell Brillouin Zone Question: Which require a finer mesh, Metals or Insulators ?? (6x6) mesh
Plane waves • Project the wave functions on a basis set • Tricky integrals become linear algebra • Plane Wave for Solid State • Could be localised (ex: Gaussians) + + = Sum of plane waves of increasing frequency (or energy) One has to stop: Ecut
Solid State: Summary • Quantities can be calculated in the direct or reciprocal space • k-point Mesh • Plane wave basis set, Ecut
Now what ? We have access to the energy of a system, without any empirical input With little efforts, the forces can be computed, Hellman-Feynman theorem • Then, the methodologies discussed for atomistic potential can be used • Energy Optimisation • Monte Carlo • Molecular dynamics