Energy and Forces in DFT

Size: px

Start display at page:

Download "Energy and Forces in DFT"

Philomena Moody
6 years ago
Views:

1 Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density n(r) Hellmann-Feynman forces: F i = de = dr i V (r) n(r) dr E II (2) R i R i where V {R} = external potential acting on electrons

2 DFT and Molecular Dynamics Let us assume i) classical nuclei, ii) electrons in the ground state. We can introduce a classical Lagrangian: L = 1 2 i M i Ṙ 2 i E tot ({R}) (3) describing the motion of nuclei. Equations of motion: d L dt Ṙi L R i = 0, P i = L Ṙi (4) are nothing but usual Newton s equations: P i M i V i, M i V i = F i (5)

3 Discretization of the equation of motion Verlet algorithm: Velocity Verlet: R i (t + δt) = 2R i (t) R i (t δt) + δt2 M i F i (t) + O(δt 4 ) (6) V i (t) = 1 2δt [R i(t + δt) R i (t δt)] + O(δt 3 ). (7) V i (t + δt) = V i (t) + δt 2M i [F i (t) + F i (t + δt)] (8) R i (t + δt) = R i (t) + δtv i (t) + δt2 2M i F i (t). (9) Both sample the microcanonical ensemble, or NVE: the energy is conserved.

4 Thermodynamical averages Averages (ρ = probability of a microscopic state): A = ρ(r, P)A(R, P)dRdP (10) are usually well approximated by averages over time: and the latter by discrete average over a trajectory: lim A T A (11) T A T = 1 T T 0 A(R(t), P(t))dt 1 M M n=1 A(t n ), t n = nδt, t M = Mδt = T. (12)

5 Costant-Temperature and constant-pressure dynamics Molecular Dynamics with Nosé Thermostats samples the canonical ensemble (NVT): the average temperature is fixed N i=1 P 2 i 2M i NV T = 3 2 Nk BT (13) Molecular Dynamics with variable cell samples the NPT ensemble: the average pressure is fixed In both cases, additional (fictitious) degrees of freedom are considered

6 Technicalities time step as big as possible, but small enough to follow nuclear motion with little loss of accuracy. Rule of thumb: δt δt max, where δt max = 1/ω max = period of the fastest phonon mode. calculations of forces must be very accurate at each time step (good self-consistency needed) or else a systematic energy drift will appear Note that: error for DFT energy is quadratic in the SCF error of the charge density error for DFT forces is linear in the SCF error of the charge density

7 Car-Parrinello Molecular Dynamics The idea: introduce a (fictitious) electron dynamics that keeps the electrons close to the ground state the electron dynamics is faster than the nuclear dynamics and averages out the error......but not too fast so that a reasonable time step can be used Very simple and effective! All classical MD technology can be used Requires a judicious choice of simulation parameters

8 Car-Parrinello Lagrangian L = µ dr 2 ψ k (r) M i Ṙ 2 i E tot ({R}, {ψ}) 2 k i + ( ) Λ kl ψk(r)ψ l (r)dr δ kl k,l (14) generates equations of motion: µ ψ k = Hψ k l Λ kl ψ l, M i Ri = E tot R i (15) µ = fictitious electronic mass Λ kl = Lagrange multipliers enforce orthonormality constraints.

9 Car-Parrinello Lagrangian (2) electronic degrees of freedom ψ k orbitals into PW are the expansion coefficients of KS forces on electrons are determined by the KS Hamiltonian calculated from current values of ψ k and of R i forces acting on nuclei have the Hellmann-Feynman form: E tot R i = k ψ k V R i ψ k (16) but they slightly differ from true forces

10 Car-Parrinello technicalities Starting point: bring the electrons to the ground state at fixed ions The simulation is performed using classical MD technology (Verlet) on both nuclear positions and electronic degrees of freedom Orthonormality constraints are exactly imposed at each time step, using an iterative procedure The choice of parameters (fictitious electronic mass µ and time step δt) must ensure that there is no energy transfer from nuclei to electrons, i.e. the kinetic energy of the electronic degrees of freedom does not systematically increase

11 Choice of the parameters The fictitious electronic mass µ and time step δt must be chosen in such a way that: µ is be big enough to enable the use of a reasonable time step, small enough to guarantee adiabaticity, i.e. no energy transfer from nuclei to electrons, which always remain close to the ground state correctness of the nuclear trajectory Typical values: electron masses δt should be the largest value that yields a stable dynamics (no drifts, no loss of orthonormality)

12 Why Car-Parrinello works The CP dynamics is classical, both for nuclei and electrons: the energy should equipartition! Typical frequencies associated to electron dynamics: ω el (ɛ i ɛ j )/µ. if there is a gap in the electronic spectrum, ω el min ɛ gap /µ Typical frequencies associated to nuclear dynamics: phonon frequencies ω ph If ω ph max << ω el min there is negligible energy transfer from nuclei to electrons A fast electron dynamics keeps the electrons close to the ground state and averages out the error on the forces. The slow nuclear dynamics is very close to the correct one, i.e. with electrons in the ground state. (Pastore, Smargiassi, Buda, PRA 44, 6334, 1991).

Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014

Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014 Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid