Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

Similar documents
Probabilistic method to determine electron correlation energy

Molecular Dynamics and Density Functional Theory

5.03, Inorganic Chemistry Prof. Daniel G. Nocera Lecture 2 May 11: Ligand Field Theory

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

PHYS 705: Classical Mechanics. Calculus of Variations II

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

Lecture 14: Forces and Stresses

Mechanics Physics 151

12. The Hamilton-Jacobi Equation Michael Fowler

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.

International Journal of Pure and Applied Sciences and Technology

CI/CEPA. Introduction CI Size Consistency Derivation CEPA EPV Results Remarks

Electronic Quantum Monte Carlo Calculations of Energies and Atomic Forces for Diatomic and Polyatomic Molecules

Key Concepts, Methods and Machinery - lecture 2 -

The GW approximation in 90 minutes or so. F. Bruneval Service de Recherches de Métallurgie Physique CEA, DEN

Generalized Linear Methods

5. Response properties in ab initio schemes

Advanced Quantum Mechanics

Numerical Solution of Ordinary Differential Equations

Density Functional Theory I

Lecture 12: Discrete Laplacian

Multiscale Modeling UCLA Prof. N. Ghoniem

Continuous Time Markov Chain

On a direct solver for linear least squares problems

Treatments of the exchange energy in density-functional theory

ESI-3D: Electron Sharing Indexes Program for 3D Molecular Space Partition

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

ψ ij has the eigenvalue

Mean Field / Variational Approximations

The Feynman path integral

Army Ants Tunneling for Classical Simulations

Classical Field Theory

A new Approach for Solving Linear Ordinary Differential Equations

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Carlo Adamo. Density Functional Theory

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Quantum Chemical Studies of Protein-Bound Chromophores, UV-Light Induced DNA Damages, and Lignin Formation

1 GSW Iterative Techniques for y = Ax

Lecture 10: Euler s Equations for Multivariable

Solutions to Problems Fundamentals of Condensed Matter Physics

Why BP Works STAT 232B

Complex Atoms; The Exclusion Principle and the Periodic System

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

Paolo Giannozzi Scuola Normale Superiore, Piazza dei Cavalieri 7 I Pisa, Italy

A how to guide to second quantization method.

Feature Selection: Part 1

Electronic Structure for Excited States (multiconfigurational methods) Spiridoula Matsika

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Grid Generation around a Cylinder by Complex Potential Functions

A Hybrid Variational Iteration Method for Blasius Equation

Inductance Calculation for Conductors of Arbitrary Shape

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Three views of mechanics

Chapter 9: Statistical Inference and the Relationship between Two Variables

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

Poisson brackets and canonical transformations

High resolution entropy stable scheme for shallow water equations

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

CHAPTER 14 GENERAL PERTURBATION THEORY

Beyond Zudilin s Conjectured q-analog of Schmidt s problem

ME 501A Seminar in Engineering Analysis Page 1

The non-negativity of probabilities and the collapse of state

Supplemental document

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

Chapter 1. Probability

Lagrangian Field Theory

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

The Symmetries of Kibble s Gauge Theory of Gravitational Field, Conservation Laws of Energy-Momentum Tensor Density and the

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Chapter 8. Potential Energy and Conservation of Energy

Supporting information.

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

1 Matrix representations of canonical matrices

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

Lecture 20: Noether s Theorem

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Numerical Heat and Mass Transfer

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Multi-electron atoms (11) 2010 update Extend the H-atom picture to more than 1 electron: H-atom sol'n use for N-elect., assume product wavefct.

Sampling Theory MODULE V LECTURE - 17 RATIO AND PRODUCT METHODS OF ESTIMATION

Affine and Riemannian Connections

24. Atomic Spectra, Term Symbols and Hund s Rules

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques

Section 8.3 Polar Form of Complex Numbers

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Multi-Time Euler-Lagrange Dynamics

Uncertainty in measurements of power and energy on power networks

Mathematics Intersection of Lines

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

Transcription:

Introducton to Densty Functonal Theory Jereme Zaffran nd year-msc. (anochemstry)

A- Hartree appromatons Born- Oppenhemer appromaton H H H e The goal of computatonal chemstry H e??? Let s remnd H e T e et ee Key-problem T e et Z 0 R r j H at and ee j r r j H e H at ee

Frst Hartree appromaton: Instead of consderng the operator ee let s consder each electron n the mean feld eperenced from all the other electrons. j ee r where r r r r d 3 r Effectve potental: eff et ee

Second Hartree appromaton: Due to the Paul ecluson prncple the system wave functon s supposed to be a sngle Slater determnant: S! where s r From Hartree-Fock model to DFT

Consequences o aratonnal prncple o Orthonormalsaton of the spn orbtals bass set E e H e H at j jj j j jj r j Jj Coulomb ntegral j j r j j Kj Echange ntegral Hartree-Fock potental: HF J K

Lmtatons of the Hartree-Fock Model. Hartree-Fock model deals wth a non-nteractng reference system and thus the correlaton energy s not taken n account.. The wave functon s a functons reles on 3 varables where s the number of electrons n the system. ery tme-consumng only small systems.b: The correlaton energy could be reached usng post Hartree-Fock methods epandng the wave functon on a bass of several Slater determnants (Confguraton Interacton-CI ) or perturbaton method (MP ) 3. The wave functon has no physcal sense only ts square has one!

B- Densty Functonal Theory (DFT) Let s set bass H e T e et ee Functonal: Mathematcal applcaton gong from the functons space to the scalars space. otaton: F[f]= whch means f F Any chemcal system s utterly defned provded one knows ts electrons number and ts eternal potental et and thus ground state energy could be reached such that E 0 E et

Hohenberg-Kohn Theorems (964) ) The eternal potental ets a unque functonal of the ground state densty 0 Consequence: Provded one knows the ground state densty one gets n turn the eternal potental and thus the hamltonan resultng n the ground state wave functon and energy (and all the system propertes). et H E 0 0 0 (and all other propertes)

Hohenberg-Kohn Theorems (964) E ) The ground state energy wll be reached f and only f one use the ground state densty 0 0 n the energy functonal. n other words the well known varatonal prncple! E E E whch means tral tral 0 0 or E 0 mn E

DFT key ponts From Hartree-Fock model to DFT The electronc densty becomes the fundamental varable! Interest: s only a functon of 4 varables (yzs) and no more of 3 varables as wth. s an observable. Any DFT algorthm should am to reach only the ground state and no ected state! The energy mnmsaton algorthms have to take care about two man constrants lyng on the densty: must be -representable whch means assocated to an acceptable wave functon : square ntegrable functons The Slater determnant s only an eample of such a set! must be et -representable whch means gvng rse to a fnte eternal potental. ote that to ths date we don t know what makes a densty et -representable on the mathematcal pont of vew. Levy constraned search scheme

Epresson of the energy functonal and lmtatons of the Hohenberg-Kohn theorems E T et E ee F HK Feature of the system Unversal functonal et e dr F HK???? ( T s not a functonal of the densty and E ee s not completely known)

Kohn-Sham approach F T HK E ee J E c Coloumbc repulson (known)? Idea The major part of the Hohenberg-Kohn functonal s the knetc energy the remander could be just appromated. So let s fnd a way to epress T

Owng to the Hartree-Fock theory T s eactly known for a non-nteractng reference system S! : Kohn-Sham (KS) orbtals And thus (teratve resoluton) Kohn-Sham equatons KS f Where the KS operator s j j S j KS r f Effectve or Sham potental ee et j S r

Hghlghts of KS approach KS-orbtals have no physcal meanng! The target s only to reach the densty. The KS-orbtals could be epressed as atomc orbtals or as Bloch waves accordng to the calculaton code. An ntal electronc densty nput s necessary Self-Consstent-Feld (SCF)

SCF scheme Intal nput densty Sham potental calculaton Sham equatons resoluton Densty output etracton Self-Consstent- Feld? O YES OUT

Focus on correlaton energy Correlaton: Mathematcal defnton: electron at r and electron at r are correlated f the followng relaton s OT verfed r r r r ~ Physcal meanng: Classcal and non classcal effects due to the many-body nteractng system. To ths date no useful epresson of the correlaton s known! If an eact functonal of the correlaton energy was known the Schrodnger equaton could have been solved EXACTLY-wthout any appromaton

But unfortunately the only mathematcal (and not so useful) formalsm we have s E c T TS Eee J T E C ncl wth T TS T C E ncl : Knetc energy of the real system : Knetc energy of the reference system : Resdual knetc energy : on-classcal energy The Kohn-Sham approach s eact only the echange-correlaton functonal to be appromated! E c has

C- The echange-correlaton problem How to appromate ths functonal? Local Densty Appromaton (LDA): Based on the homogeneous electrons gas model. Echange-correlaton densty functonal s eactly known owng to the Thomas-Ferm model. Gradent Generalzed Appromaton (GGA): PBE Applcaton of the gradent operator on the prevous model. Meta-GGA: BB95 Applcaton also of the laplacan operator. Hybrd functonal: HSE06 B3LYP Introducton of an eact Hartree-Fock part n the Echange functonal. E c % E ( %) E HF X GGA X E GGA C

Jacob s ladder (Pedrew metaphor) Heaven: Eact soluton? Hybrd MGGA GGA LDA Earth: HF model