Zurich, July 12th 2017
CP2K : GPW and GAPW
Marcella Iannuzzi
Department of Chemistry, University of Zurich
http://www.cp2k.org 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
2 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-consistency 3to required precision 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
Augmented basis set, GAPW: separated hard and soft density domains
4 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 5 Gaussian Basis Set
Localised, atom-position dependent GTO basis