4th CP2K tutorial
August 31st -September 4th, Zurich
GPW (GAPW) electronic structure calculations
Marcella Iannuzzi
Department of Chemistry, University of Zurich
http://www.cp2k.org Outline
Density Functional Theory and the KS formalism
Gaussian and Plane Wave method (GPW)
Basis sets and pseudo potentials
Gaussian Augmented Plane Wave method (GAPW)
2 DFT
Why DFT?
Explicit inclusion of electronic structure
Predicable accuracy (unlike empirical approaches, parameter free)
Knowledge of electronic structure gives access to evaluation of many observables
Better scaling compared to many quantum chemistry approaches
Achievable improvements: development of algorithms and functionals
large systems, condensed matter, environment effects, first principle MD
3 Hohenberg-Kohn theorems
Motivation Density functional theory History Theorem I Bloch theorem / supercells Kohn-Sham method History of DFT — II Given a potential, one obtains the wave functions via Schrödinger equation
Vext(r, R) H(r, R)=T (r)+Vext(r, R)+Vee(r) ) Walter Kohn H(r, R) (r, R)=E(R) (r, R) Walter Kohn
DFT essentials CECAM NMR/EPR tutorial 2013 5 / 37
The density is the probability distribution of the wave functions
n(r) V (r, R) , ext
the potential and hence also the total energy are unique functional of the electronic density n(r)
4 HK Total energy
Theorem II: The total energy is variational
E[n] E[n ] GS
Etot[n]=Ekin[n]+Eext[n]+EH[n]+Exc[n]
Ekin QM kinetic energy of electron (TF)
Eext energy due to external potential
EH classical Hartree repulsion
Exc non classical Coulomb energy: el. correlation
5 Kohn-Sham: non-interacting electrons Electronic density 2 n(r)= fi i(r) | | no repulsion i X Kinetic energy of non interacting electrons 1 T [n]= f (r) 2 (r) s i i | 2r | i i X ⌧ Electronic interaction with the external potential
ZI Eext[n]= n(r)Vext(r)dr Vext(r)= r r RI Z XI | | 1 Exact solution s = det [ 1 2 3... N ] pN! 6 KS energy functional
Classical e-e repulsion
1 n(r)n(r0) 1 J[n]= drdr0 = n(r)V (r)dr 2 r r 2 H Zr Zr0 | 0| Zr
Kohn-Sham functional
EKS[n]=Ts[n]+Eext[n]+J[n]+EXC[n]
E [n]=E [n] T [n]+E [n] J[n] XC kin s ee non-classical︸ part 7 KS Equations
Orthonormality constraint
⌦ [ ]=E [n] ✏ ⇤(r) (r)dr KS i KS ij i j ij X Z Lagrange︸ multipliers ⌦ [ ] Variational search in the space of orbitals KS i =0 i⇤
1 H = 2 + V = ✏ KS i 2r KS i ij j ij X
VKS(r)=Vext(r)+VH(r)+VXC(r)
8 KS Equations
!ij diagonal 1 2 + V (r) (r)=✏ (r) 2r KS i i i
KS equations looking like Schrödinger equations
coupled and highly non linear
Self consistent solution required
! and ψ are help variables
KS scheme in principle exact (Exc?)
9 Self-consistency SCF Method Generate a starting density ⇒ ninit Input 3D Coordinates init Generate the KS potential ⇒ VKS of atomic nuclei
Solve the KS equations ⇒ , ψ ! Initial Guess Molecular Orbitals Fock Matrix (1-electron vectors) Calculation 1 Calculate the new density ⇒ n
1 New KS potential ⇒ VKS Fock Matrix Diagonalization
New orbitals and energies ⇒ !1 , ψ
Calculate 2 Properties Yes SCF No New density ⇒ n Converged? End …..
until self-consistency10 to required precision Basis Set Representation
KS matrix formulation when the wavefunction is expanded into a basis
System size {Nel, M}, P [MxM], C [MxN]
⇥i(r)= C i (r) Variational n(r)= fiC iC⇥i (r) ⇥(r)= P ⇥ (r) ⇥(r) principle i ⇥ ⇥ Constrained minimisation P = PSP problem KS total energy
E[ ]=T [ ]+Eext[n]+EH[n]+EXC[n]+EII { i} { i} Matrix formulation of the KS equations
H xc K(C)C = T(C)+Vext(C)+E (C)+E (C)=SC
11 Critical Tasks Construction of the Kohn-Sham matrix
Hartree potential Introduction Energy minimization and sparseness function XC potential Time reversible BOMD Summary HF/exact exchange Why are (N) methods so important? O Fast and robust minimisation of theWith energy conventional SCF methods, hardware improvements bring functional only small gains in capability due to the steep scaling of computational time with system size, N. Efficient calculation of the density matrix and construction of the MOs (C)
O(N) scaling in basis set size
Big systems: biomolecules, interfaces, material science 1000+ atoms
Long time scale: 1 ps = 1000 MD steps, processes several ps a day
12
Valéry Weber Classes of Basis Sets
Extended basis sets, PW : condensed matter
Localised basis sets centred at atomic positions, GTO
Idea of GPW: auxiliary basis set to represent the density
Mixed (GTO+PW) to take best of two worlds, GPW: over-completeness
Augmented basis set, GAPW: separated hard and soft density domains
13 GPW Ingredients linear scaling KS matrix computation for GTO
Gaussian basis sets (many terms analytic)
2 mx my mz ↵mr ⇥i(r)= C i (r) ↵(r)= dm↵gm(r) gm(r)=x y z e m X
Pseudo potentials
Plane waves auxiliary basis for Coulomb integrals
Regular grids and FFT for the density
Sparse matrices (KS and P)
Efficient screening
G. Lippert et al, Molecular Physics, 92, 477, 1997 J. VandeVondele et al, Comp. Phys. Comm.,167 (2), 103, 2005 14 Gaussian Basis Set
Localised, atom-position dependent GTO basis