Basics of Density Functional Theory (DFT)
|
|
- Damon Wilcox
- 6 years ago
- Views:
Transcription
1 Basics of Density Functional Theory (DFT) Ari Paavo SEITSONEN Département de Chimie École Normale Supérieure, Paris École de Sidi-BelAbbès de Nanomateriaux // Octobre 8-12, 2016 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 1 / 94
2 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 2 / 94
3 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 3 / 94
4 Summary 1 Density functional theory Motivation Kohn-Sham method Exchange-correlation functional 2 Bloch theorem / supercells 3 Localised basis sets 4 Plane wave basis set 5 Pseudo potentials apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 4 / 94
5 Motivation: Why use DFT? Explicit inclusion of electronic structure Predictable accuracy (unlike fitted/empirical approaches) Knowledge of the electron structure can be used for the analysis; many observables can be obtained directly Preferable scaling compared to many quantum chemistry methods apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 5 / 94
6 History of DFT I There were already methods in the early 20th century Thomas-Fermi-method Hartree-Fock-method apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 6 / 94
7 History of DFT II Walter Kohn apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 7 / 94
8 History of DFT III: Foundations apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 8 / 94
9 Hohenberg-Kohn theorems: Theorem I Given a potential, one obtains the wave functions via Schrödinger equation: V (r) ψ i (r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 9 / 94
10 Hohenberg-Kohn theorems: Theorem I Given a potential, one obtains the wave functions via Schrödinger equation: V (r) ψ i (r) The density is the probability distribution of the wave functions: n (r) = i ψ i (r) 2 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 9 / 94
11 Hohenberg-Kohn theorems: Theorem I Given a potential, one obtains the wave functions via Schrödinger equation: V (r) ψ i (r) The density is the probability distribution of the wave functions: n (r) = i ψ i (r) 2 Thus V (r) ψ i (r) n (r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 9 / 94
12 Hohenberg-Kohn theorems: Theorem I Theorem The potential, and hence also the total energy, is a unique functional of the electron density n(r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 10 / 94
13 Hohenberg-Kohn theorems: Theorem I Theorem The potential, and hence also the total energy, is a unique functional of the electron density n(r) Thus V (r) ψ i (r) n (r) V (r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 10 / 94
14 Hohenberg-Kohn theorems: Theorem I Theorem The potential, and hence also the total energy, is a unique functional of the electron density n(r) Thus V (r) ψ i (r) n (r) V (r) The electron density can be used to determine all properties of a system apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 10 / 94
15 Hohenberg-Kohn theorems: Theorem II Theorem The total energy is variational: In the ground state the total energy is minimised apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 11 / 94
16 Hohenberg-Kohn theorems: Theorem II Theorem The total energy is variational: In the ground state the total energy is minimised Thus E[n] E[n GS ] apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 11 / 94
17 History of DFT IV: Foundations apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 12 / 94
18 History of DFT V: The reward... in 1998: apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 13 / 94
19 Kohn-Sham method: Total energy Let us write the total energy as: E tot [n] = E kin [n]+e ext [n]+e H [n]+e xc [n] E kin [n] = QM kinetic energy of electrons E ext [n] = energy due to external potential (usually ions) E H [n] = classical Hartree repulsion (e e ) E xc [n] = exchange-correlation energy apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 14 / 94
20 Kohn-Sham method: Noninteracting electrons To solve the many-body Schrödinger equation as such is an unformidable task apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 15 / 94
21 Kohn-Sham method: Noninteracting electrons To solve the many-body Schrödinger equation as such is an unformidable task Let us write the many-body wave function as a determinant of single-particle equations Then kinetic energy of electrons becomes E kin,s = i 1 2 f i ψi (r) 2 ψ i (r) f i = occupation of orbital i (with spin-degeneracy included) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 15 / 94
22 Kohn-Sham method: External energy Energy due to external potential; usually V ext = I E ext = r n (r) V ext (r) dr n (r) = f i ψ i (r) 2 i Z I r R I apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 16 / 94
23 Kohn-Sham method: Hartree energy Classical electron-electron repulsion E H = 1 n (r) n (r ) 2 r r r r dr dr = 1 n (r) V H (r) dr 2 r n (r ) V H (r) = r r r dr apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 17 / 94
24 Kohn-Sham method: Exchange-correlation energy The remaining component: Many-body complications combined = Will be discussed later apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 18 / 94
25 Total energy expression Kohn-Sham (total 1 ) energy: E KS [n] = 1 2 f i ψi 2 ψ i + n (r) V ext (r) dr i r + 1 n (r) n (r ) 2 r r r r dr dr + E xc 1 without ion-ion interaction apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 19 / 94
26 Kohn-Sham equations Vary the Kohn-Sham energy E KS with respect to ψ j (r ): δe KS δψ j (r ) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 20 / 94
27 Kohn-Sham equations Vary the Kohn-Sham energy E KS with respect to ψ j (r ): δe KS δψ j (r ) Kohn-Sham equations { 1 } V KS (r) ψ i (r) = ε i ψ i (r) n (r) = i f i ψ i (r) 2 V KS (r) = V ext (r) + V H (r) + V xc (r) V xc (r) = δexc δn(r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 20 / 94
28 Kohn-Sham equations: Notes { 1 } V KS (r) ψ i (r) = ε i ψ i (r) ; n (r) = i f i ψ i (r) 2 Equation looking like Schrödinger equation The Kohn-Sham potential, however, depends on density The equations are coupled and highly non-linear Self-consistent solution required ε i and ψ i are in principle only help variables (only ε HOMO has a meaning) The potential V KS is local The scheme is in principle exact apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 21 / 94
29 Kohn-Sham equations: Self-consistency 1 Generate a starting density n init 2 Generate the Kohn-Sham potential V init KS 3 Solve the Kohn-Sham equations ψ init i 4 New density n 1 5 Kohn-Sham potential V 1 KS 6 Kohn-Sham orbitals ψ 1 i 7 Density n apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 22 / 94
30 Kohn-Sham equations: Self-consistency 1 Generate a starting density n init 2 Generate the Kohn-Sham potential V init KS 3 Solve the Kohn-Sham equations ψ init i 4 New density n 1 5 Kohn-Sham potential V 1 KS 6 Kohn-Sham orbitals ψ 1 i 7 Density n until self-consistency is achieved (to required precision) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 22 / 94
31 Kohn-Sham equations: Self-consistency Usually the density coming out from the wave functions is mixed with the previous ones, in order to improve convergence In metals fractional occupations numbers are necessary The required accuracy in self-consistency depends on the observable and the expected apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 23 / 94
32 Kohn-Sham energy: Alternative expression Take the Kohn-Sham equation, multiply from the left with f i ψi and integrate: 1 2 f i ψ i (r) 2 ψ i (r) dr + f i V KS (r) ψ i (r) 2 dr = f i ε i r Sum over i and substitute into the expression for Kohn-Sham energy: E KS [n] = f i ε i E H + E xc n (r) V xc dr i r r apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 24 / 94
33 Exchange-correlation functional The Kohn-Sham scheme is in principle exact however, the exchange-correlation energy functional is not known explicitly Exchange is known exactly, however its evaluation is usually very time-consuming, and often the accuracy is not improved over simpler approximations (due to error-cancellations in the latter group) Many exact properties, like high/low-density limits, scaling rules etc are known Famous approximations: Local density approximation, LDA Generalised gradient approximations, GGA Hybrid functionals... apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 25 / 94
34 Kohn-Sham method: Exchange-correlation energy The remaining component: Many-body complications combined apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 26 / 94
35 Total energy expression Kohn-Sham (total 2 ) energy: E KS [n] = 1 2 f i ψi 2 ψ i + n (r) V ext (r) dr i r + 1 n (r) n (r ) 2 r r r r dr dr + E xc 2 without ion-ion interaction apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 27 / 94
36 Kohn-Sham equations Vary the Kohn-Sham energy E KS with respect to ψ j (r ): δe KS δψ j (r ) Kohn-Sham equations { 1 } V KS (r) ψ i (r) = ε i ψ i (r) n (r) = i f i ψ i (r) 2 V KS (r) = V ext (r) + V H (r) + V xc (r) V xc (r) = δexc δn(r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 28 / 94
37 Exchange-correlation functional The Kohn-Sham scheme is in principle exact however, the exchange-correlation energy functional is not known explicitly Exchange is known exactly, however its evaluation is usually very time-consuming, and often the accuracy is not improved over simpler approximations (due to error-cancellations in the latter group) Many exact properties, like high/low-density limits, scaling rules etc are known Famous approximations: Local density approximation, LDA Generalised gradient approximations, GGA Hybrid functionals... apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 29 / 94
38 Exchange-correlation functional: LDA Local density approximation: Use the exchange-correlation energy functional for homogeneous electron gas at each point of space: E xc r n (r) eheg xc [n(r)]dr Works surprisingly well even in inhomogeneous electron systems, thanks fulfillment of certain sum rules Energy differences over-bound: Cohesion, dissociation, adsorption energies Lattice constants somewhat (1-3 %) too small, bulk modulii too large Asymptotic potential decays exponentially outside charge distribution, leads to too weak binding energies for the electrons, eg. ionisation potentials Contains self-interaction in single-particle case apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 30 / 94
39 Exchange-correlation functional: GGA Generalised gradient approximation: The gradient expansion of the exchange-correlation energy does not improve results; sometimes leads to divergencies Thus a more general approach is taken, and there is room for several forms of GGA: E xc r n (r) e xc[n(r), n(r) 2 ]dr Works reasonably well, again fulfilling certain sum rules Energy differences are improved Lattice constants somewhat, 1-3 % too large, bulk modulii too small Contains self-interaction in single-particle case Again exponential asymptotic decay in potential Negative ions normally not bound Usually the best compromise between speed and accuracy in large systems apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 31 / 94
40 Exchange-correlation functional: Hybrid functionals Hybrid functionals Include partially the exact (Hartree-Fock) exchange: E xc αe HF + (1 α)ex GGA + Ec GGA ; again many variants Works in general well Energy differences are still improved Vibrational properties slightly Improved magnetic moments in some systems Partial improvement in asymptotic form Usually the best accuracy if the computation burden can be handled Calculations for crystals appearing apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 32 / 94
41 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
42 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
43 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
44 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
45 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
46 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
47 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
48 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
49 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
50 Exchange-correlation functional: Observations The accuracy can not be systematically improved! van der Waals interactions still a problem (tailored approximations in sight) Most widely used parametrisations: GGA Perdew-Burke-Ernzerhof, PBE (1996); among phycisists Beck Lee-Yang-Parr, BLYP (early 1990 s); among chemists Hybrid functionals PBE0; among phycisists B3LYP; among chemists van der Waals included DFT-D3; empirical (pair-)potential vdw-df*; rvv10? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 33 / 94
51 Summary 1 Density functional theory 2 Bloch theorem / supercells 3 Localised basis sets 4 Plane wave basis set 5 Pseudo potentials apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 34 / 94
52 Periodic systems In realistic systems there are atoms in cubic millimetre unformidable to treat by any numerical method At this scale the systems are often repeating (crystals)... or the observable is localised and the system can be made periodic Choices: Periodic boundary conditions or isolated (saturated) cluster apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 35 / 94
53 Periodic systems apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 36 / 94
54 Periodic systems apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 36 / 94
55 Periodic systems apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 36 / 94
56 Periodic systems apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 36 / 94
57 Periodic systems Is it possible to replace the summation over translations L with a modulation? apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 37 / 94
58 Periodic systems Is it possible to replace the summation over translations L with a modulation? Bloch s theorem For a periodic potential V (r + L) = V (r) the eigenfunctions can be written in the form ψ i (r) = e ik r u ik (r), u ik (r + L) = u ik (r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 37 / 94
59 Periodic systems: Reciprocal space Reciprocal lattice vectors: b 1 = 2π a 2 a 3 a 1 a 2 a 3 b 2 = 2π a 3 a 1 a 2 a 3 a 1 b 3 = 2π a 1 a 2 a 3 a 1 a 2 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 38 / 94
60 Periodic systems: Brillouin zone First Brillouin zone: Part of space closer to the origin than to any integer multiple of the reciprocal lattice vectors, K = n 1 b 1 + n 2 b 2 + n 3 b 3 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 39 / 94
61 Integration over reciprocal space Thus the summation over infinite number of translations becomes an integral over the first Brillouin zone: L k 1.BZ In practise the integral is replaced by a weighted sum of discrete points: dk Thus eg. w k dk k k n (r) = w k f ik ψ ik (r) 2 k i apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 40 / 94
62 Periodic systems: Dispersion apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 41 / 94
63 Band structure: Example Pb/Cu(111) Photoemission vs DFT calculations for a free-standing layer Felix Baumberger, Anna Tamai, Matthias Muntwiler, Thomas Greber and Jürg Osterwalder; Surface Science (2003) doi: /s (03) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 42 / 94
64 Monkhorst-Pack algorithm Approximate the integral with an equidistance grid of k vectors with identical weight: apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 43 / 94
65 Monkhorst-Pack algorithm Approximate the integral with an equidistance grid of k vectors with identical weight: n = 2p q 1 2q, p = 1... q k ijk = n 1 b 1 + n 2 b 2 + n 3 b 3 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 43 / 94
66 Symmetry operations If the atoms are related by symmetry operation S (Sψ (r) = ψ (Sr)) the integration over the whole 1st Brillouin zone can be reduced into the irreducible Brillouin zone, IBZ apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 44 / 94
67 Symmetry operations If the atoms are related by symmetry operation S (Sψ (r) = ψ (Sr)) the integration over the whole 1st Brillouin zone can be reduced into the irreducible Brillouin zone, IBZ Sψ ik (r) = ψ ik (Sr) = e ik Sr u ik (Sr) = e ik r u ik (r), k = S 1 k apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 44 / 94
68 Symmetry operations If the atoms are related by symmetry operation S (Sψ (r) = ψ (Sr)) the integration over the whole 1st Brillouin zone can be reduced into the irreducible Brillouin zone, IBZ Sψ ik (r) = ψ ik (Sr) = e ik Sr u ik (Sr) = e ik r u ik (r), k = S 1 k dk w k = k k BZ k IBZ S w Sk apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 44 / 94
69 Irreducible Brillouin zone: Examples apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 45 / 94
70 Doubling the unit cell apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 46 / 94
71 Doubling the unit cell (super-cells) If one doubles the unit cell in one direction, it is enough to take only half of the k points in the corresponding direction in the reciprocal space apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 47 / 94
72 Doubling the unit cell (super-cells) If one doubles the unit cell in one direction, it is enough to take only half of the k points in the corresponding direction in the reciprocal space And has to be careful when comparing energies in cells with different size unless either equivalent sampling of k points is used apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 47 / 94
73 Doubling the unit cell (super-cells) If one doubles the unit cell in one direction, it is enough to take only half of the k points in the corresponding direction in the reciprocal space And has to be careful when comparing energies in cells with different size unless either equivalent sampling of k points is used or one is converged in the total energy in both cases apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 47 / 94
74 Summary 1 Density functional theory 2 Bloch theorem / supercells 3 Localised basis sets 4 Plane wave basis set 5 Pseudo potentials apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 48 / 94
75 Implementation of Kohn-Sham equations: Discretisation Kohn-Sham orbitals expanded in a basis set: ψ ik = α c ikαχ α Different choices of basis set: Localised: Gaussian (x j y k z l e (x/σ)2 ), Slater (x j y k z l e (x/σ) ), wavelet,... Extended: Plane waves Combination: FLAPW, GPW/GAPW,... apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 49 / 94
76 Summary 1 Density functional theory 2 Bloch theorem / supercells 3 Localised basis sets 4 Plane wave basis set Basics of plane wave basis set Operators Energy terms in plane wave basis set 5 Pseudo potentials apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 50 / 94
77 Kohn Sham method The ground state energy is obtained as the solution of a constrained minimisation of the Kohn-Sham energy: min E KS[{Φ i (r)}] {Φ} apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 51 / 94
78 Kohn Sham method The ground state energy is obtained as the solution of a constrained minimisation of the Kohn-Sham energy: min E KS[{Φ i (r)}] {Φ} Φ i (r)φ j (r)dr = δ ij apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 51 / 94
79 Expansion using a basis set For practical purposes it is necessary to expand the Kohn-Sham orbitals using a set of basis functions Basis set {ϕ α (r)} M α=1 Usually a linear expansion ψ i (r) = M c αi ϕ α (r) α=1 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 52 / 94
80 Plane waves Philosophy Assemblies of atoms are slight distortions to free electrons ϕ α (r) = 1 Ω e igα r (... = cos(g α r) + i sin(g α r)) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 53 / 94
81 Plane waves Philosophy Assemblies of atoms are slight distortions to free electrons ϕ α (r) = 1 Ω e igα r + orthogonal apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 53 / 94
82 Plane waves Philosophy Assemblies of atoms are slight distortions to free electrons ϕ α (r) = 1 Ω e igα r + orthogonal + independent of atomic positions apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 53 / 94
83 Plane waves Philosophy Assemblies of atoms are slight distortions to free electrons ϕ α (r) = 1 Ω e igα r + orthogonal + independent of atomic positions + no BSSE apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 53 / 94
84 Plane waves Philosophy Assemblies of atoms are slight distortions to free electrons ϕ α (r) = 1 Ω e igα r + orthogonal + independent of atomic positions + no BSSE ± naturally periodic apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 53 / 94
85 Plane waves Philosophy Assemblies of atoms are slight distortions to free electrons ϕ α (r) = 1 Ω e igα r + orthogonal + independent of atomic positions + no BSSE ± naturally periodic many functions needed apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 53 / 94
86 Computational box Box matrix : h = [a 1, a 2, a 3 ] Box volume : Ω = det h apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 54 / 94
87 Lattice vectors Direct lattice h = [a 1, a 2, a 3 ] Translations in direct lattice: L = i a 1 + j a 2 + k a 3 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 55 / 94
88 Lattice vectors Direct lattice h = [a 1, a 2, a 3 ] Translations in direct lattice: L = i a 1 + j a 2 + k a 3 Reciprocal lattice 2π(h t ) 1 = [b 1, b 2, b 3 ] b i a j = 2πδ ij Reciprocal lattice vectors : G = i b 1 + j b 2 + k b 3 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 55 / 94
89 Expansion of Kohn-Sham orbitals Plane wave expansion ψ ik (r) = G c ik (G)e i(k+g) r To be solved: Coefficients c ik (G) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 56 / 94
90 Expansion of Kohn-Sham orbitals Plane wave expansion ψ ik (r) = G c ik (G)e i(k+g) r To be solved: Coefficients c ik (G) Different routes: Direct optimisation of total energy Iterative diagonalisation/minimisation apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 56 / 94
91 Dependence on position Translation: φ(r) φ(r R I ) φ(r R I ) = G = G φ(g)e ig (r R I ) ( φ(g)e ig r) e ig R I apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 57 / 94
92 Dependence on position Translation: φ(r) φ(r R I ) φ(r R I ) = G = G φ(g)e ig (r R I ) ( φ(g)e ig r) e ig R I Structure Factor: S I (G) = e ig R I apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 57 / 94
93 Dependence on position Translation: φ(r) φ(r R I ) φ(r R I ) = G = G φ(g)e ig (r R I ) ( φ(g)e ig r) e ig R I Structure Factor: Derivatives: φ(r R I ) R I,s S I (G) = e ig R I = i G G s (φ(g)e ig r) S I (G) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 57 / 94
94 Plane waves: Kinetic energy Kinetic energy operator in the plane wave basis: ϕ G (r) = 1 2 (ig)2 1 Ω e ig r = 1 2 G 2 ϕ G (r) Thus the operator is diagonal in the plane wave basis set E kin (G) = 1 2 G 2 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 58 / 94
95 Cutoff: Finite basis set Choose all basis functions into the basis set that fulfill 1 2 G 2 E cut a cut-off sphere apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 59 / 94
96 Cutoff: Finite basis set Choose all basis functions into the basis set that fulfill 1 2 G 2 E cut a cut-off sphere apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 59 / 94
97 Cutoff: Finite basis set Choose all basis functions into the basis set that fulfill 1 2 G 2 E cut a cut-off sphere N PW 1 3/2 ΩE 2π2 cut [a.u.] Basis set size depends on volume of box and cutoff only apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 59 / 94
98 Cutoff: Finite basis set Choose all basis functions into the basis set that fulfill 1 2 G 2 E cut a cut-off sphere N PW 1 3/2 ΩE 2π2 cut [a.u.] Basis set size depends on volume of box and cutoff only and is variational! apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 59 / 94
99 Plane waves: Fast Fourier Transform The information contained in ψ(g) and ψ(r) are equivalent ψ(g) ψ(r) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 60 / 94
100 Plane waves: Fast Fourier Transform The information contained in ψ(g) and ψ(r) are equivalent ψ(g) ψ(r) Transform from ψ(g) to ψ(r) and back is done using fast Fourier transforms (FFT s) Along one direction the number of operations N log[n] apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 60 / 94
101 Plane waves: Fast Fourier Transform The information contained in ψ(g) and ψ(r) are equivalent ψ(g) ψ(r) Transform from ψ(g) to ψ(r) and back is done using fast Fourier transforms (FFT s) Along one direction the number of operations N log[n] 3D-transform = three subsequent 1D-transforms apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 60 / 94
102 Plane waves: Fast Fourier Transform The information contained in ψ(g) and ψ(r) are equivalent ψ(g) ψ(r) Transform from ψ(g) to ψ(r) and back is done using fast Fourier transforms (FFT s) Along one direction the number of operations N log[n] 3D-transform = three subsequent 1D-transforms Information can be handled always in the most appropriate space apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 60 / 94
103 Plane waves: Integrals Parseval s theorem Ω G A (G)B(G) = Ω A (r i )B(r i ) N i apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 61 / 94
104 Plane waves: Integrals Parseval s theorem Ω G A (G)B(G) = Ω A (r i )B(r i ) N i Proof. I = A (r)b(r)dr Ω = A (G)B(G) GG exp[ ig r] exp[ig r]dr = GG A (G)B(G) Ω δ GG = Ω G A (G)B(G) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 61 / 94
105 Plane waves: Electron density n(r) = ik w k f ik ψ ik (r) 2 = 1 w k f ik cik Ω (G)c ik(g )e i(g G ) r ik G,G n(r) = 2G max G= 2G max n(g)e ig r apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 62 / 94
106 Plane waves: Electron density n(r) = ik w k f ik ψ ik (r) 2 = 1 w k f ik cik Ω (G)c ik(g )e i(g G ) r ik G,G n(r) = 2G max G= 2G max n(g)e ig r The electron density can be expanded exactly in a plane wave basis with a cut-off four times the basis set cutoff. N PW (4E cut ) = 8N PW (E cut ) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 62 / 94
107 Plane waves: Operators The Kohn-Sham equations written in reciprocal space: { 1 } V KS (G, G ) ψ ik (G) = ε i ψ ik (G) However, it is better to do it like Car and Parrinello (1985) suggested: Always use the appropriate space (via FFT) There one needs to apply an operator on a wave function: O(G, G ( )ψ G ) = ( c G ) G O G G G Matrix representation of operators in: O(G, G ) = G O G Eg. Kinetic energy operator T G,G = G G = 1 2 G 2 δ G,G apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 63 / 94
108 Kohn Sham energy E KS = E kin + E ES + E pp + E xc E kin Kinetic energy E ES Electrostatic energy (sum of electron-electron interaction + nuclear core-electron interaction + ion-ion interaction) E pp Pseudo potential energy not included in E ES E xc Exchange correlation energy apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 64 / 94
109 Kinetic energy E kin = ik w k f ik ψ ik ψ ik = w k f ik cik (G)c ik(g ) k + G k + G ik GG = w k f ik cik (G)c ik(g ) Ω 1 2 k + G 2 δ G,G ik GG = Ω 1 w k f ik 2 k + G 2 c ik (G) 2 ik G apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 65 / 94
110 Periodic Systems Hartree-like terms are most efficiently evaluated in reciprocal space via the Poisson equation 2 V H (r) = 4πn tot (r) V H (G) = 4π n(g) G 2 V H (G) is a local operator with same cutoff as n tot apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 66 / 94
111 Electrostatic energy E ES = 1 2 n(r)n(r ) r r dr dr + I n(r)vcore(r)dr I + 1 Z I Z J 2 R I r J I J apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 67 / 94
112 Electrostatic energy E ES = 1 2 n(r)n(r ) r r dr dr + I n(r)vcore(r)dr I + 1 Z I Z J 2 R I r J I J The isolated terms do not converge; the sum only for neutral systems Gaussian charge distributions a la Ewald summation: [ nc(r) I = Z ( ) ] I r 2 ( ) R c 3 π 3/2 RI exp I Electrostatic potential due to nc: I n Vcore(r) I I = c (r ) r r dr = Z [ ] I r r R I erf RI R c I R c I apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 67 / 94
113 Electrostatic energy Electrostatic energy E ES = 2π Ω n tot (G) 2 G 2 G 0 + E ovrl E self E ovrl = Z I Z J R I r J L erfc R I r J L R c I 2 + R c J 2 I,J L E self = I 1 Z I 2 2π Sums expand over all atoms in the simulation cell, all direct lattice vectors L; the prime in the first sum indicates that I < J is imposed for L = 0. apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 68 / 94 R c I
114 Exchange-correlation energy E xc = n(r)ε xc (r)dr = Ω r G ε xc (G)n (G) ε xc (G) is not local in G space; the calculation in real space requires very accurate integration scheme. If the function ε xc (r) requires the gradients of the density, they are calculated using reciprocal space, otherwise the calculation is done in real space (for LDA and GGA; hybrid functionals are more intensive) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 69 / 94
115 Plane waves: Basic self-consistent cycle n(r) V xc [n](r) V KS (r) = V xc [n](r) + V ES (r) ψ ik (r) V KS (r)ψ ik (r) ψ ik (r) n (r) = ik w kf ik ψ ik (r) 2 FFT FFT N i FFT N i FFT N i FFT n(g) V ES (G) V ES (G) ψ ik (G) [V KS ψ ik ] (G) update ψ ik (G) ψ ik (G) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 70 / 94
116 Plane waves: Calculation of forces With the plane wave basis set one can apply the Hellmann-Feynman theorem (F I =) d dr I Ψ H KS Ψ = Ψ R I H KS Ψ apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 71 / 94
117 Plane waves: Calculation of forces With the plane wave basis set one can apply the Hellmann-Feynman theorem (F I =) d dr I Ψ H KS Ψ = Ψ R I H KS Ψ All the terms where R I appear explicitly are in reciprocal space, and are thus very simple to evaluate: R I e ig R I = ige ig R I apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 71 / 94
118 Plane waves: Summary Plane waves are delocalised, periodic basis functions apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
119 Plane waves: Summary Plane waves are delocalised, periodic basis functions Plenty of them are needed, however the operations are simple apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
120 Plane waves: Summary Plane waves are delocalised, periodic basis functions Plenty of them are needed, however the operations are simple The quality of basis set adjusted using a single parametre, the cut-off energy apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
121 Plane waves: Summary Plane waves are delocalised, periodic basis functions Plenty of them are needed, however the operations are simple The quality of basis set adjusted using a single parametre, the cut-off energy Fast Fourier-transform used to efficiently switch between real and reciprocal space apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
122 Plane waves: Summary Plane waves are delocalised, periodic basis functions Plenty of them are needed, however the operations are simple The quality of basis set adjusted using a single parametre, the cut-off energy Fast Fourier-transform used to efficiently switch between real and reciprocal space Forces and Hartree term/poisson equation are trivial apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
123 Plane waves: Summary Plane waves are delocalised, periodic basis functions Plenty of them are needed, however the operations are simple The quality of basis set adjusted using a single parametre, the cut-off energy Fast Fourier-transform used to efficiently switch between real and reciprocal space Forces and Hartree term/poisson equation are trivial The system has to be neutral! Usual approach for charged states: Homogeneous neutralising background apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
124 Plane waves: Summary Plane waves are delocalised, periodic basis functions Plenty of them are needed, however the operations are simple The quality of basis set adjusted using a single parametre, the cut-off energy Fast Fourier-transform used to efficiently switch between real and reciprocal space Forces and Hartree term/poisson equation are trivial The system has to be neutral! Usual approach for charged states: Homogeneous neutralising background The energies must only be compared with the same E cut apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 72 / 94
125 Summary 1 Density functional theory 2 Bloch theorem / supercells 3 Localised basis sets 4 Plane wave basis set 5 Pseudo potentials Introduction to pseudo potentials Generation of pseudo potentials Usage of pseudo potentials apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 73 / 94
126 Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10 instead of 78 Unnecessary Why bother? They are inert anyway... Inclusion of relativistic effects relativistic effects can be included partially into effective potentials apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 74 / 94
127 Why pp? Estimate for number of plane waves plane wave cutoff most localized function 1s Slater type function exp[ Zr] Z: effective nuclear charge φ 1s (G) 16πZ 5/2 G 2 + Z 2 Cutoff Plane waves H 1 1 Li 4 8 C 9 27 Si Ge apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 75 / 94
128 Pseudo potential What is it? Replacement of the all-electron, Z/r problem with a Hamiltonian containing an effective potential It should reproduce the necessary physical properties of the full problem at the reference state The potential should be transferable, ie. also be accurate in different environments The construction consists of two steps of approximations Frozen core approximation Pseudisation apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 76 / 94
129 Frozen core approximation Core electrons are chemically inert Core/valence separation is sometimes not clear Core wavefunctions are from atomic reference calculation Core electrons of different atoms do not overlap apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 77 / 94
130 Remaining problems Valence wavefunctions must be orthogonal to core states nodal structures high plane wave cutoff Thus pseudo potential should produce node-less functions and include Pauli repulsion Pseudo potential replaces Hartree and XC potential due to the core electrons apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 78 / 94
131 Remaining problems Valence wavefunctions must be orthogonal to core states nodal structures high plane wave cutoff Thus pseudo potential should produce node-less functions and include Pauli repulsion Pseudo potential replaces Hartree and XC potential due to the core electrons XC functionals are not linear: approximation E XC (n c + n v ) = E XC (n c ) + E XC (n v ) This assumes that core and valence electrons do not overlap. This restriction can be overcome with the non linear core correction (NLCC) discussed later apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 78 / 94
132 Atomic pseudo potentials: Basic idea 1 Create a pseudo potential for each atomic species separately, in a reference state (usually a neutral or slightly ionised atom) 2 Use these in the Kohn-Sham equations for the valence only ( T + V (r, r ) + V H (n v ) + V XC (n v ) ) Φ v i (r) = ɛ i Φ v i (r) Goal Pseudo potential V (r, r ) has to be chosen such that the main properties of the all-electron atom are imitated as well as possible apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 79 / 94
133 Pseudisation of valence wave functions apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 80 / 94
134 Generation of pseudo potentials: General recipe 1 Atomic all electron calculation (reference state) Φ v l (r) and ɛ l. 2 Pseudize Φ v l Φ PS l 3 Calculate potential from [T + V l (r)] Φ PS l (r) = ɛ l Φ PS (r) 4 Calculate pseudo potential by unscreening of V l (r) l V PS l V PS is dependent l (r) = on V l l! (r) Otherwise V H (n PS poor ) V transferability XC (n PS ) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 81 / 94
135 Norm-conserving pseudo potentials Hamann-Schlüter-Chiang-Recipe (HSC) conditions D R Hamann, M Schlüter & C Chiang, Phys Rev Lett 43, 1494 (1979) 1 Real and pseudo valence eigenvalues agree for a chosen prototype atomic configuration: ɛ l = ˆɛ l 2 Real and pseudo atomic wave functions agree beyond a chosen core radius r c : Ψ l (r) = Φ l (r) for r r c 3 Integral from 0 to R of the real and pseudo charge densities agree for R r c (norm conservation): Φ l Φ l R = Ψ l Ψ l R for R r c, Φ Φ R = R 0 r 2 φ(r) 2 dr 4 The logarithmic derivatives of the real and pseudo wave function and their first energy derivatives agree for r r c. Points 3) and 4) are related via 1 [ 2 (rφ) 2 d d dɛ dr lnφ] R = R 0 r 2 Φ 2 dr apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 82 / 94
136 Semi-local pseudo potential Final semi-local form L max V PS (r, r ) = Vloc PS (r) + L=0 V PS L (r) Y L Y L Local pseudo potential Vloc PS Non-local pseudo potential VL PS = VL PS Vloc PS Any L quantum number can have a non-local part This form leads, however, to large computing requirements in plane wave basis sets apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 83 / 94
137 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 84 / 94
138 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 85 / 94
139 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 86 / 94
140 Kleinman Bylander form Write the pseudo potential terms as Kleinman-Bylander operator ( r, r ) = V L φ L ω L V L φ L ˆV KB PS,nl φ L is the pseudo wave function E PS,nl = ik L ω L = φ L V L φ L w k f ik ψ ik V L φ L ω L V L φ L ψ ik L Now the expression is fully non-local {f (r ) f (r)} and thus computationally efficient also in plane wave basis set! apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 87 / 94
141 Ghost states Problem: In Kleinman Bylander form, the node-less wave function is no longer the solution with the lowest energy Solution: Carefully tune the local part of the pseudo potential until the ghost states disappear How to detect ghost states: Look for following properties Deviations of the logarithmic derivatives of the energy of the KB pseudo potential from those of the respective semi-local pseudo potential or all electron potential. Comparison of the atomic bound state spectra for the semi-local and KB pseudo potentials. Ghost states below the valence states are identified by a rigorous criteria by Gonze et al apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 88 / 94
142 Ultra soft pseudo potentials and PAW method Many elements require high cutoff for plane wave calculations First row elements: O, F Transition metals: Cu, Zn f elements: Ce relax norm-conservation condition n PS (r)dr + Q(r)dr = 1 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 89 / 94
143 Ultra soft pseudo potentials and PAW method Augmentation functions Q(r) depend on environment. No full un-screening possible, Q(r) has to be recalculated for each atom and atomic position. Complicated orthogonalisation and force calculations. Allows for larger r c, reduces cutoff for all elements to about 30 Ry apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 90 / 94
144 Non-linear core correction (NLCC) For many atoms (eg. alkali atoms, transition metals) core states overlap with valence states; thus linearisation of XC energy breaks down Possible cures Add additional states to valence adds more electrons needs higher cutoff Add core charge to valence charge in XC energy non linear core correction (NLCC) S.G. Louie et al., Phys. Rev. B, (1982) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 91 / 94
145 Testing of pseudo potentials Calculation of other atomic valence configurations (occupations), comparison with all-electron results Calculation of transferability functions, logarithmic derivatives, hardness Calculation of small molecules, comparison with all electron calculations (geometry, harmonic frequencies, dipole moments) Calculation of other test systems Check of basis set convergence (cutoff requirements) apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 92 / 94
146 Convergence test for pseudo potentials in O 2 apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 93 / 94
147 Pseudo potentials: Summary Pseudo potential are necessary when using plane wave basis sets in order to keep the number of the basis function manageable Pseudo potentials are generated at the reference state; transferability is the quantity describing the accuracy of the properties at other conditions The mostly used scheme in plane wave calculations is the Troullier-Martins pseudo potentials in the fully non-local, Kleinman-Bylander form Once created, a pseudo potential must be tested, tested, tested!!! apsi (ENS, UPMC, CNRS) DFT Sidi-BelAbbès 94 / 94
Why use pseudo potentials?
Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10
More informationDensity functional theory in the solid state
Density functional theory in the solid state Ari P Seitsonen IMPMC, CNRS & Universités 6 et 7 Paris, IPGP Department of Applied Physics, Helsinki University of Technology Physikalisch-Chemisches Institut
More information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationThe electronic structure of materials 2 - DFT
Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie
More informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
More informationDFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More informationElectronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory
Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
More informationDFT: Exchange-Correlation
DFT: Exchange-Correlation Local functionals, exact exchange and other post-dft methods Paul Tulip Centre for Materials Physics Department of Physics University of Durham Outline Introduction What is exchange
More informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationDFT: Exchange-Correlation
DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,
More informationCLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122
CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationPlane waves, pseudopotentials and PAW. X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium
Plane waves, pseudopotentials and PAW X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium 1 Basic equations in DFT Solve self-consistently the Kohn-Sham equation H ψ n = ε n ψ n!!! ρ(r
More informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » Recap of HK theorems and KS equations» The physical meaning of the XC energy» Solution of a one-particle Schroedinger equation» Pseudo Potentials»
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
More informationPractical Guide to Density Functional Theory (DFT)
Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationHow to generate a pseudopotential with non-linear core corrections
How to generate a pseudopotential with non-linear core corrections 14 12 AE core charge AE valence charge PS core charge PS valence charge 10 8 6 4 2 Objectives 0 0 0.5 1 1.5 2 2.5 3 Check whether the
More informationDFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
More informationAri P Seitsonen CNRS & Université Pierre et Marie Curie, Paris
Self-organisation on noble metal surfaces Ari P Seitsonen CNRS & Université Pierre et Marie Curie, Paris Collaborations Alexandre Dmitriev, Nian Lin, Johannes Barth, Klaus Kern,... Thomas Greber, Jürg
More informationCHAPTER 3 WIEN2k. Chapter 3 : WIEN2k 50
CHAPTER 3 WIEN2k WIEN2k is one of the fastest and reliable simulation codes among computational methods. All the computational work presented on lanthanide intermetallic compounds has been performed by
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationAB INITIO MODELLING TECHNIQUES APPLIED TO SILICON. P.R.BRIDDON Department of Physics, University of Newcastle, Newcastle upon Tyne, NE1 7RU, UK
AB INITIO MODELLING TECHNIQUES APPLIED TO SILICON P.R.BRIDDON Department of Physics, University of Newcastle, Newcastle upon Tyne, NE1 7RU, UK August 13, 1999 1 Introduction In this chapter we will consider
More informationIntegrated Computational Materials Engineering Education
Integrated Computational Materials Engineering Education Lecture on Density Functional Theory An Introduction Mark Asta Dept. of Materials Science and Engineering, University of California, Berkeley &
More informationAdvanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University
Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation
More informationAdvanced Solid State Theory SS Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt
Advanced Solid State Theory SS 2010 Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt i 0. Literatur R. M. Martin, Electronic Structure: Basic Theory and
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationJournal of Theoretical Physics
1 Journal of Theoretical Physics Founded and Edited by M. Apostol 53 (2000) ISSN 1453-4428 Ionization potential for metallic clusters L. C. Cune and M. Apostol Department of Theoretical Physics, Institute
More informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationThe Schrödinger equation for many-electron systems
The Schrödinger equation for many-electron systems Ĥ!( x,, x ) = E!( x,, x ) 1 N 1 1 Z 1 Ĥ = " $ # " $ + $ 2 r 2 A j j A, j RAj i, j < i a linear differential equation in 4N variables (atomic units) (3
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationElectron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals
Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University
More informationAll electron optimized effective potential method for solids
All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22
More informationStrategies for Solving Kohn- Sham equations
Strategies for Solving Kohn- Sham equations Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 1 Appetizer: high-k oxides
More informationGW quasiparticle energies
Chapter 4 GW quasiparticle energies Density functional theory provides a good description of ground state properties by mapping the problem of interacting electrons onto a KS system of independent particles
More informationTerminal exercise: DFT via Quantum ESPRESSO
Terminal exercise: DFT via Quantum ESPRESSO Ari Paavo SEITSONEN Ari.P.Seitsonen@iki.fi Département de Chimie École Normale Supérieure, Paris École de Sidi-BelAbbès de Nanomateriaux // Octobre 8-12, 2016
More informationDFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
More informationThe Linearized Augmented Planewave (LAPW) Method
The Linearized Augmented Planewave (LAPW) Method David J. Singh Oak Ridge National Laboratory E T [ ]=T s [ ]+E ei [ ]+E H [ ]+E xc [ ]+E ii {T s +V ks [,r]} I (r)= i i (r) Need tools that are reliable
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationThe Projector Augmented Wave method
The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An
More informationElectronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier Pôle de chimie théorique Département de chimie de l ENS CNRS Ecole normale supérieure UPMC Formation ModPhyChem Lyon,
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationCHEM6085: Density Functional Theory
Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron
More informationGeneralized generalized gradient approximation: An improved density-functional theory for accurate orbital eigenvalues
PHYSICAL REVIEW B VOLUME 55, NUMBER 24 15 JUNE 1997-II Generalized generalized gradient approximation: An improved density-functional theory for accurate orbital eigenvalues Xinlei Hua, Xiaojie Chen, and
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationReferences. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
More informationNorm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria
Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized
More informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
More informationProjector augmented wave Implementation
Projector augmented wave Implementation Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 = Projector augmented wave +
More informationElectronic Structure Calculations, Density Functional Theory and its Modern Implementations
Tutoriel Big RENOBLE Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Thierry Deutsch L_Sim - CEA renoble 19 October 2011 Outline 1 of Atomistic calculations
More informationSpring College on Computational Nanoscience May Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor
2145-25 Spring College on Computational Nanoscience 17-28 May 2010 Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor Stefano BARONI SISSA & CNR-IOM DEMOCRITOS Simulation Center
More informationPseudopotential methods for DFT calculations
Pseudopotential methods for DFT calculations Lorenzo Paulatto Scuola Internazionale Superiore di Studi Avanzati and CNR-INFM DEMOCRITOS National Simulation Center Tieste Italy July 9, 2008 Outline pseudopotential
More informationDensity Functional Theory (DFT)
Density Functional Theory (DFT) An Introduction by A.I. Al-Sharif Irbid, Aug, 2 nd, 2009 Density Functional Theory Revolutionized our approach to the electronic structure of atoms, molecules and solid
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationFundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering
Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationBasics of DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Basics of DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Basics of DFT Lausanne12 1
More informationDensity Functional Theory (DFT) modelling of C60 and
ISPUB.COM The Internet Journal of Nanotechnology Volume 3 Number 1 Density Functional Theory (DFT) modelling of C60 and N@C60 N Kuganathan Citation N Kuganathan. Density Functional Theory (DFT) modelling
More informationStudy of Carbon Nanomaterials Based on Density Functional Theory. Mohammad Shafiul Alam
Study of Carbon Nanomaterials Based on Density Functional Theory Mohammad Shafiul Alam July 2013 Dissertation Study of Carbon Nanomaterials Based on Density Functional Theory Graduate School of Natural
More informationAb initio study of Mn doped BN nanosheets Tudor Luca Mitran
Ab initio study of Mn doped BN nanosheets Tudor Luca Mitran MDEO Research Center University of Bucharest, Faculty of Physics, Bucharest-Magurele, Romania Oldenburg 20.04.2012 Table of contents 1. Density
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationIntermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Intermediate DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Intermediate DFT Lausanne12
More informationABC of ground-state DFT
ABC of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) ABC of ground-state
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More information3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional
More informationElectrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation
Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
More informationNotes on Density Functional Theory
Notes on Density Functional Theory 1 Basic Theorems The energy, E, of a system with a given Hamiltonian H is a functional of the (normalized, many-particle) wave function Ψ. We write this functional as
More information