Electronic Structure Calculations, Density Functional Theory and its Modern Implementations
|
|
- Scot Hicks
- 5 years ago
- Views:
Transcription
1 Tutoriel Big RENOBLE Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Thierry Deutsch L_Sim - CEA renoble 19 October 2011
2 Outline 1 of Atomistic calculations 2 Density Functional Theory (quick view) formalism equations 3 basis sets 4
3 of Atomistic simulations
4 of Atomistic Simulations Theory Experiment Simulation Hardware Computers Algorithms Atomistic Simulations Force fields (interatomic potentials) Tight Binding Methods Hartree-Fock Density Functional Theory Configuration interactions Quantum Monte-Carlo
5 of Atomistic Simulations Theory Experiment Simulation Hardware Computers Algorithms Atomistic Simulations Force fields (interatomic potentials) Tight Binding Methods Hartree-Fock Density Functional Theory Configuration interactions Quantum Monte-Carlo
6 Algorithms Search of efficient algorithms: Fast Fourier Transform (FFT) : T CPU n log(n) ; Order N atoms methods Electrostatic (Hartree) potential of an electronic density ρ(r): V(r) = V Mesh: n = = 10 6 Direct: n 2 = : 100 s FFT: n log(n) 10 6 : 0.1 ms ρ(r ) r r dr
7 Algorithms Search of efficient algorithms: Fast Fourier Transform (FFT) : T CPU n log(n) ; Order N atoms methods Electrostatic (Hartree) potential of an electronic density ρ(r): V(r) = V Mesh: n = = 10 6 Direct: n 2 = : 100 s FFT: n log(n) 10 6 : 0.1 ms ρ(r ) r r dr
8 Atomistic Simulations Two intrinsic difficulties for numerical atomistic simulations, related to complexity: Interactions The way that atoms interact is known: i ħ Ψ t = H Ψ H ψ = E 0 ψ Exploration of the configuration space (ART, minima hopping, spline search)
9 Atomistic Simulations Two intrinsic difficulties for numerical atomistic simulations, related to complexity: Interactions The way that atoms interact is known: i ħ Ψ t = H Ψ H ψ = E 0 ψ Exploration of the configuration space (ART, minima hopping, spline search) E pot R1
10 Atomistic Simulations Two intrinsic difficulties for numerical atomistic simulations, related to complexity: Interactions The way that atoms interact is known: i ħ Ψ t = H Ψ H ψ = E 0 ψ Exploration of the configuration space (ART, minima hopping, spline search) E pot R1 R2
11 Atomistic Simulations Two intrinsic difficulties for numerical atomistic simulations, related to complexity: Interactions The way that atoms interact is known: i ħ Ψ t = H Ψ H ψ = E 0 ψ Exploration of the configuration space (ART, minima hopping, spline search) E pot R1 R3 R2
12 Atomistic Simulations Two intrinsic difficulties for numerical atomistic simulations, related to complexity: Interactions The way that atoms interact is known: i ħ Ψ t = H Ψ H ψ = E 0 ψ Exploration of the configuration space (ART, minima hopping, spline search) E pot R1000 R41 R3 Rn R2 R1
13 oals of Criteria eneral (enerality) Precision geometry atomic positions lengths of bonds angles of bonds energy of bonds vibration frequencies dipolar moments UV and visible spectra photo-emission magnetic moments System size Time scale oals all atoms of the periodic table all kinds of bonds (transition state) ±0.001 angström ±0.001 angström ±1 degree ev/atom (0.1 kj/mol) 1 Hz 0.01Debye (1 Debye = C.m 0.01 ev 0.01 ev 0.01µ B 10 4 atoms for ab initio methods 10 6 atoms for parametrized methods 1 s
14 Criteria for classification 3 criteria 1 enerality (elements, alloys) 2 Precision ( r, E) 3 System size (N, t) P S
15 Quantum Mechanics Schrödinger equation for N electrons in an external potential V ext (r): antisymmetric Ψ(r 1,r 2,...,r N ) = Ψ(r 2,r 1,...,r N ) H = 2 2m N i 2 r i + 1 4πε 0 i j 1 r i r j + N i V ext (r i ) H Ψ(r 1,r 2,...,r N ) = EΨ(r 1,r 2,...,r N ) Ion-electron interactions: V ext (r) = 1 Z α 4πε 0 R α r α
16 Hartree-Fock and Post Hartree-Fock We build the most simple antisymmetric wavefunction for N electrons from N wavefunctions at 1 electron (Slater determinant): φ 1 (x 1 )... φ 1 (x N ) φ 2 (x 1 )... φ 2 (x N ) Φ(x 1,...,x N ) = φ N (x 1 )... φ N (x N ) The Hamiltonian are exact. Many Slater determinants: Configuration interaction, MP2, coupled clusters,... P P S S
17 Hartree-Fock and Post Hartree-Fock We build the most simple antisymmetric wavefunction for N electrons from N wavefunctions at 1 electron (Slater determinant): φ 1 (x 1 )... φ 1 (x N ) φ 2 (x 1 )... φ 2 (x N ) Φ(x 1,...,x N ) = φ N (x 1 )... φ N (x N ) The Hamiltonian are exact. Many Slater determinants: Configuration interaction, MP2, coupled clusters,... P P S S
18 Quantum Monte Carlo Very complicated wavefunctions for N electrons (Jastrow factor): φ 1 (x 1 )... φ 1 (x N ) Ψ(x 1,...,x N ) = φ N (x 1 )... φ N (x N ) i j J (r ij ) Integration of the total energy: E = < Ψ H Ψ > < Ψ Ψ > Use of Monte-Carlo methods for the integration. P S Diffusive Monte Carlo (fixed node, sign problem)
19 Density Functional Theory () Theorem: The total energy of the ground state in an external potential V ext is a functional of the electronic density ρ(r): E [ρ(r)] = F [ρ(r)] + V ext (r)ρ(r)dr r Variational principles: E [ρ 0 (r)] E [ρ(r)] ρ(r) Equation: (approximation) Having a one-electron hamiltonian in a mean field. [ V eff (r) ] ψ i = ε i ψ i P occupied ρ(r) = i ψ i (r) 2 S
20 Density Functional Theory () Theorem: The total energy of the ground state in an external potential V ext is a functional of the electronic density ρ(r): E [ρ(r)] = F [ρ(r)] + V ext (r)ρ(r)dr r Variational principles: E [ρ 0 (r)] E [ρ(r)] ρ(r) Equation: (approximation) Having a one-electron hamiltonian in a mean field. [ V eff (r) ] ψ i = ε i ψ i P occupied ρ(r) = i ψ i (r) 2 S
21 Tight Binding Methods Semi-empirical quantum methods We solve a one-electron parametrized hamiltonian (Slater-Koster parameter): H ψ i >= ε i ψ i > occupied Band term: E band = ions Repulsive term: E rep = I,J i ε i f (R I,J ) P E total = E rep + E band S Mainly used for electronic properties (10 6 atoms) or transport
22 Choice of Chemistry Physics P P Compromise between accuracy, system size and CPU time. Force Fields Tight Binding Hartree-Fock S S P S S P Ab initio () Configuration Interactions P P Quantum Monte-Carlo S S
23 Choice of Chemistry Physics P P Compromise between accuracy, system size and CPU time. Force Fields Tight Binding Hartree-Fock S S P S S P Ab initio () Configuration Interactions P P Quantum Monte-Carlo S S
24 Choice of Chemistry Physics P P Compromise between accuracy, system size and CPU time. Force Fields Tight Binding Hartree-Fock S S P S S P Ab initio () Configuration Interactions P P Quantum Monte-Carlo S S
25 Choice of Chemistry Physics P P Compromise between accuracy, system size and CPU time. Force Fields Tight Binding Hartree-Fock S S P S S P Ab initio () Configuration Interactions P P Quantum Monte-Carlo S S
26 Choice of Chemistry Physics P P Compromise between accuracy, system size and CPU time. Force Fields Tight Binding Hartree-Fock S S P S S P Ab initio () Configuration Interactions P P Quantum Monte-Carlo S S
27 Choice of Chemistry Physics P P Compromise between accuracy, system size and CPU time. Force Fields Tight Binding Hartree-Fock S S P S S P Ab initio () Configuration Interactions P P Quantum Monte-Carlo S S
28 Density Functional Theory equations
29 Ab initio calculations with Several advantages Ab initio: No adjustable parameters : Quantum mechanical (fundamental) treatment Main limitations Approximated approach Requires high computer power, limited to few hundreds atoms in most cases Wide range of applications: nanoscience, biology, materials
30 Quantum Mechanics Schrödinger equation for N electrons in an external potential V ext (r): antisymmetric Ψ(r 1,r 2,...,r N ) = Ψ(r 2,r 1,...,r N ) H = 2 2m N i 2 r i + 1 4πε 0 i j 1 r i r j + N i V ext (r i ) H Ψ(r 1,r 2,...,r N ) = EΨ(r 1,r 2,...,r N ) Ion-electron interactions: V ext (r) = 1 Z α 4πε 0 R α r α
31 The Hohenberg-Kohn theorem Schrödinger equation H = N i=1 ( r i + V ext (r i,{r}) ) i j Very difficult to solve for more than two electrons! 1 r i r j The fundamental variable of the problem is however not the wavefunction, but the electronic density ρ(r) = N dr 2 dr N ψ (r,r 2,,r N )ψ(r,r 2,,r N ) Hohenberg-Kohn theorem (1964) The ground state density ρ(r) of a many-electron system uniquely determines (up to a constant) the external potential. The external potential is a functional of ρ: V ext = V ext [ρ]
32 The Hohenberg-Kohn theorem Schrödinger equation H = N i=1 ( r i + V ext (r i,{r}) ) i j Very difficult to solve for more than two electrons! 1 r i r j The fundamental variable of the problem is however not the wavefunction, but the electronic density ρ(r) = N dr 2 dr N ψ (r,r 2,,r N )ψ(r,r 2,,r N ) Hohenberg-Kohn theorem (1964) The ground state density ρ(r) of a many-electron system uniquely determines (up to a constant) the external potential. The external potential is a functional of ρ: V ext = V ext [ρ]
33 formalism H-K theorem: E is an unknown functional of the density E = E[ρ] Density Functional Theory approach Mapping of an interacting many-electron system into a system with independent particles moving into an effective potential. Find a set of orthonormal orbitals Ψ i (r) that minimizes: E = 1 2 N i=1 Ψ i (r) 2 Ψ i (r)dr + 1 ρ(r)v H (r)dr 2 + E xc [ρ(r)] + V ext (r)ρ(r)dr ρ(r) = N i=1 Ψ i (r)ψ i (r) 2 V H (r) = 4πρ(r)
34 formalism H-K theorem: E is an unknown functional of the density E = E[ρ] Density Functional Theory approach Mapping of an interacting many-electron system into a system with independent particles moving into an effective potential. Find a set of orthonormal orbitals Ψ i (r) that minimizes: E = 1 2 N i=1 Ψ i (r) 2 Ψ i (r)dr + 1 ρ(r)v H (r)dr 2 + E xc [ρ(r)] + V ext (r)ρ(r)dr ρ(r) = N i=1 Ψ i (r)ψ i (r) 2 V H (r) = 4πρ(r)
35 formalism H-K theorem: E is an unknown functional of the density E = E[ρ] Density Functional Theory approach Mapping of an interacting many-electron system into a system with independent particles moving into an effective potential. Find a set of orthonormal orbitals Ψ i (r) that minimizes: E = 1 2 N i=1 Ψ i (r) 2 Ψ i (r)dr + 1 ρ(r)v H (r)dr 2 + E xc [ρ(r)] + V ext (r)ρ(r)dr ρ(r) = N i=1 Ψ i (r)ψ i (r) 2 V H (r) = 4πρ(r)
36 Performing a calculation A self-consistent equation ρ(r) = Ψ i (r)ψ i (r), where ψ i satisfies i ( 12 ) 2 + V H [ρ] + V xc [ρ] + V ext + V pseudo ψ i = E i ψ i, () Ingredients An XC potential, functional of the density several approximations exists (LDA,A,... ) A choice of the pseudopotential (if not all-electrons) (norm conserving, ultrasoft, PAW,... ) A basis set for expressing the ψ i An (iterative) algorithm for finding the wavefunctions ψ i A (good) computer...
37 Exchange-Correlation Energy E xc [ρ] = K exact [ρ] K KS [ρ] + V e e exact [ρ] V Hartree [ρ] K KS [ρ]: Kinetic energy for a non-interaction electron gas is an exact reformulation of the many-body problem formalism is exact (mapping of an interacting electron system into a non-interacting system Many types: Local Density Approximation (LDA) eneralized gradient approximation (A), meta-a, Hybrid functional (a part of exchange energy) Random Phase Approximation (RPA),...
38 Exchange-Correlation Energy E xc [ρ] = K exact [ρ] K KS [ρ] + V e e exact [ρ] V Hartree [ρ] K KS [ρ]: Kinetic energy for a non-interaction electron gas is an exact reformulation of the many-body problem formalism is exact (mapping of an interacting electron system into a non-interacting system Many types: Local Density Approximation (LDA) eneralized gradient approximation (A), meta-a, Hybrid functional (a part of exchange energy) Random Phase Approximation (RPA),...
39 Local Density Approximation (LDA) Assumption: For each point r in space with a density ρ(r), the exchange and correlation energy at that point is the same as for an uniform electron gas with that density: Exc LDA = ρ(r)ε xc(ρ)dr V V LDA xc (r) = E LDA xc ρ(r) = ε xc (ρ) + ρ(r) ε xc(ρ) ρ(r) ε xc (ρ): Exchange-correlation energy for a uniform electron gas of ( density (analytic): Exc LDA = 3 3π ) 1/3 4 ρ(r) 4/3 dr Correlation energy (interpolated analytically): Based on Quantum Monte Carlo calculations (Ceperley and Alder (1980)) for the exchange and correlation of the electron gas. Based on uniform gas of electron: ood cancellation of error
40 eneralized radient Approximations (A) Use the gradient of the density to have a more accurate exchange energy Many generalised gradient approximations: Perdew (1985), Lee, Yang and Parr (1988), Perdew and Wang (1991), PBE (1996),HCTH (1998),... Perdew and Yue (1986) functional (which is relatively simple): s = E A xc [ρ] = 3 4 ρ(r) 2k F (ρ(r)) ( 3 π ) 1/3 ρ(r) 4/3 F(s)dr k F = ((3π 2 ρ(r)) 1/3 F(s) = ( s s s 6 ) 1/15 LDA: Overbind (energies too high) underestimate bond lengths A: In general better energies specially PBE
41 eneralized radient Approximations (A) Use the gradient of the density to have a more accurate exchange energy Many generalised gradient approximations: Perdew (1985), Lee, Yang and Parr (1988), Perdew and Wang (1991), PBE (1996),HCTH (1998),... Perdew and Yue (1986) functional (which is relatively simple): s = E A xc [ρ] = 3 4 ρ(r) 2k F (ρ(r)) ( 3 π ) 1/3 ρ(r) 4/3 F(s)dr k F = ((3π 2 ρ(r)) 1/3 F(s) = ( s s s 6 ) 1/15 LDA: Overbind (energies too high) underestimate bond lengths A: In general better energies specially PBE
42 Hybrid functionals Add a small part of exchange interaction (25% to 50%) E Exchange = i<j V Many hybrid functionals: B3LYP, HSE,... φ i (r)φ j(r)φ i (r )φ j (r ) r r drdr The gap is (really) better and also cohesion energies Very time-consuming (10) Not the right physics and sometimes not so accurate...
43 Hybrid functionals Add a small part of exchange interaction (25% to 50%) E Exchange = i<j V Many hybrid functionals: B3LYP, HSE,... φ i (r)φ j(r)φ i (r )φ j (r ) r r drdr The gap is (really) better and also cohesion energies Very time-consuming (10) Not the right physics and sometimes not so accurate...
44 Performing a calculation A self-consistent equation ρ(r) = Ψ i (r)ψ i (r), where ψ i satisfies i ( 12 ) 2 + V H [ρ] + V xc [ρ] + V ext + V pseudo ψ i = E i ψ i, () Ingredients An XC potential, functional of the density several approximations exists (LDA,A,... ) A choice of the pseudopotential (if not all-electrons) (norm conserving, ultrasoft, PAW,... ) A basis set for expressing the ψ i An (iterative) algorithm for finding the wavefunctions ψ i A (good) computer...
45 s smoothening of wavefunctions For chemical properties only the valence electrons are relevant: Eliminate the chemically inactive core electrons Reduce the number of electron orbitals The pseudo-wavefunctions of the valence electrons are smooth Eliminate the rapid variations of the valence wavefunction in the core region A reasonable approximation V ext (r) = 1 4πε 0 R α r The pseudopotential approximation is less severe than the approximate nature of the exchange correlation functional α Z α
46 Non-Local s V PP (r,r ) = Y lm (r)v l (r)δ r,r Y lm (r ) l,m where Y l,m are spherical harmonics. V NL = φ I lm> V l <φ I lm I,l,m where I is the index of atom, l and m and quantum numbers. Use ab initio atomic wavefunctions to obtain pseudopotentials. The (potential) energy contribution due to the valence charge density must be substracted away to give an ionic pseudopotential: V ion l,ps(r) = V l,ps (r) + (V Hartree [ρ(r)] + V xc [ρ(r)]) free atom
47 aussian type separable s (HH) Local part V loc (r) = Z [ ] [ ion r erf + exp 1 ( ) ] 2 r r 2rloc 2 r loc { ( ) 2 ( ) 4 ( r r r C 1 + C 2 + C 3 + C 4 r loc r loc Nonlocal (separable) part H( r, r ) H sep ( r, r ) = p l 1(r) = 2 r l+ 3 2 l 2 i=1 r loc Y s,m (ˆr) p s i (r) h s i p s i (r ) Ys,m(ˆr ) m + Y p,m (ˆr) p p 1 (r) hp 1 pp 1 (r ) Yp,m(ˆr ) m r l+ 7 2 l ) 6 } r l e 1 2 ( r ) 2 r l p Γ(l l r l+2 e 1 2 ( r ) 2 r l 2(r) = ) Γ(l + 7) 2
48 Performing a calculation A self-consistent equation ρ(r) = Ψ i (r)ψ i (r), where ψ i satisfies i ( 12 ) 2 + V H [ρ] + V xc [ρ] + V ext + V pseudo ψ i = E i ψ i, () Ingredients An XC potential, functional of the density several approximations exists (LDA,A,... ) A choice of the pseudopotential (if not all-electrons) (norm conserving, ultrasoft, PAW,... ) A basis set for expressing the ψ i An (iterative) algorithm for finding the wavefunctions ψ i A (good) computer...
49 for electronic structure calculation How can we express the wavefunctions? Plane Waves Localization in Fourier space, efficient preconditioning Systematic convergence properties No localization in real space. Empty regions must be filled with PW. Non adaptive, Slater type Orbitals localized, well suited for molecules and other open structures Small number of basis functions for moderate accuracy Many different recipes for generating basis sets Over-completeness before convergence. Non systematic basis set.
50 for electronic structure? Wavelets A basis set both adaptive and systematic, real space based Wavelet basis sets Localized both in real and in Fourier space Allow for adaptivity (for internal electrons) Systematic basis set
51 Performing a calculation A self-consistent equation ρ(r) = Ψ i (r)ψ i (r), where ψ i satisfies i ( 12 ) 2 + V H [ρ] + V xc [ρ] + V ext + V pseudo ψ i = E i ψ i, () Ingredients An XC potential, functional of the density several approximations exists (LDA,A,... ) A choice of the pseudopotential (if not all-electrons) (norm conserving, ultrasoft, PAW,... ) A basis set for expressing the ψ i An (iterative) algorithm for finding the wavefunctions ψ i A (good) computer...
52 Equations: Operators Apply different operators Having a one-electron hamiltonian in a mean field. [ 1 ] V eff (r)(r) ψ i = ε i ψ i V eff (r) = V ext (r) + V ρ(r ) r r dr + µ xc (r) E[ρ] can be expressed by the orthonormalized states of one particule: ψ i (r) with the fractional occupancy number f i (0 f i 1) : ρ(r) = f i ψ i (r) 2 i
53 Equations: Computing Energies Calculate different integrals E[ρ] = K [ρ] + U[ρ] U[ρ] = K [ρ] = m e i V dr ψ i 2 ψ i dr V ext(r)ρ(r)+ 1 ρ(r)ρ(r ) dr dr V 2 V r r }{{} Hartree We minimise with the variables ψ i (r) and f i with the constraint dr ρ(r) = N el. V + E xc [ρ] }{{} exchange correlation
54 KS Equations: Self-Consistent Field Set of self-consistent equations: { 1 } V eff ψ i = ε i ψ i 2 m e with an effective potential: V eff (r) = V ext (r)+ ) dr ρ(r V r r }{{} Hartree δe xc + δρ(r) }{{} exchange correlation and: ρ(r) = i f i ψ i (r) 2
55 KS Equations: Self-Consistent Field Set of self-consistent equations: { 1 } V eff ψ i = ε i ψ i 2 m e with an effective potential: V eff (r) = V ext (r)+ ) dr ρ(r V r r }{{} Hartree δe xc + δρ(r) }{{} exchange correlation and: ρ(r) = i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 ) Real Mesh (100 3 = 10 6 ): = evaluations!
56 KS Equations: Self-Consistent Field Set of self-consistent equations: { 1 } V eff ψ i = ε i ψ i 2 m e with an effective potential: V eff (r) = V ext (r)+ ) dr ρ(r V r r }{{} Hartree δe xc + δρ(r) }{{} exchange correlation and: ρ(r) = i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 ) Real Mesh (100 3 = 10 6 ): = evaluations!
57 KS Equations: Self-Consistent Field Set of self-consistent equations: { 1 } V eff ψ i = ε i ψ i 2 m e with an effective potential: V eff (r) = V ext (r)+ ) dr ρ(r V r r }{{} Hartree δe xc + δρ(r) }{{} exchange correlation and: ρ(r) = i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 ) Real Mesh (100 3 = 10 6 ): = evaluations!
58 KS Equations: Self-Consistent Field Set of self-consistent equations: { 1 } V eff ψ i = ε i ψ i 2 m e with an effective potential: V eff (r) = V ext (r)+ ) dr ρ(r V r r }{{} Hartree δe xc + δρ(r) }{{} exchange correlation and: ρ(r) = i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 ) Real Mesh (100 3 = 10 6 ): = evaluations!
59 KS Equations: Self-Consistent Field Set of self-consistent equations: { 1 } V eff ψ i = ε i ψ i 2 m e with an effective potential: V eff (r) = V ext (r)+ ) dr ρ(r V r r }{{} Hartree δe xc + δρ(r) }{{} exchange correlation and: ρ(r) = i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 ) Real Mesh (100 3 = 10 6 ): = evaluations!
60 Ab Initio Methods: [ ] V(r) +µ xc (r) ψ k i = εk i ψ k i
61 Ab Initio Methods: [ ] V(r) +µ xc (r) ψ k i = ε k i ψ k i atomic orbitals plane waves augmented Slater numerical real space finite difference Wavelet
62 Ab Initio Methods: self-consistent (SCF) Harris-Foulkes functional pseudopotential all electrons norm-conserving PAW [ V(r) ] +µ xc (r) ψ k i = ε k i ψ k i atomic orbitals plane waves augmented Slater numerical real space finite difference Wavelet
63 Ab Initio Methods: self-consistent (SCF) Harris-Foulkes functional pseudopotential all electrons norm-conserving PAW [ V(r) periodic non-periodic ] +µ xc (r) ψ k i atomic orbitals plane waves augmented real space = ε k i ψ k i Slater numerical finite difference Wavelet
64 Ab Initio Methods: self-consistent (SCF) Harris-Foulkes functional pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W method [ V(r) periodic non-periodic LDA,A ] +µ xc (r) ψ k i atomic orbitals plane waves augmented real space LDA+U = ε k i ψ k i Slater numerical finite difference Wavelet
65 Ab Initio Methods: self-consistent (SCF) Harris-Foulkes functional pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W [ V(r) periodic non-periodic non-collinear collinear non-spin polarized spin polarized LDA,A ] +µ xc (r) ψ k i atomic orbitals plane waves augmented real space LDA+U = ε k i ψ k i Slater numerical finite difference Wavelet
66 Ab Initio Methods: self-consistent (SCF) Harris-Foulkes functional relativistic pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W non-relativistic [ periodic non-periodic non-collinear collinear +V(r) non-spin polarized spin polarized LDA,A ] +µ xc (r) ψ k i = ε k i ψ k i atomic orbitals plane waves augmented real space LDA+U Slater numerical finite difference Wavelet
67 Ab Initio Methods: self-consistent (SCF) Harris-Foulkes functional relativistic pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W non-relativistic [ N 3 scaling O(N) methods periodic non-periodic non-collinear collinear +V(r) non-spin polarized spin polarized LDA,A ] +µ xc (r) ψ k i = ε k i ψ k i atomic orbitals plane waves augmented real space LDA+U Slater numerical finite difference Wavelet
68 List of ab initio Codes Plane Waves ABINIT Louvain-la-Neuve CPMD Zurich, Lugano PWSCF Italy VASP Vienna aussian aussian DeMon CP2K Siesta Madrid (numerical basis sets) Wien Vienna (FPLAPW, all electrons) basis set ONETEP Big
69 Minimisation of the electronic density
70 Plane Wave basis sets Natural basis set for electronic structure calculations in periodic solids. 1 V exp(i r) where V is the volume of the periodicity volume and is a multiple of the reciprocal lattice vectors, i.e. there are integers l 1,l 2,l 3 such that = l 1 b 1 + l 2 b 2 + l 3 b 3.
71 Plane Wave basis sets Matrix elements of the kinetic energy part of the Hamiltonian calculated analytically exp( i r) ( 12 ) 2 exp(i r) = 2 δ, The potential energy part obtained numerically by calculating the wave-functions in real space Ψ(r) = exp(i r) and then integrating numerically in real space. The transformation on a grid in real space is done with the help of the Fast Fourier transformation at a cost of N log 2 (N) operations for a basis set of N plane waves.
72 : Planes Waves Operator approach: no need to set up Hamiltonian matrix FFT allows for quasi linear M log(m) scaling where M is number of plane waves Localization in Fourier space allows for efficient preconditioning techniques. Hence number of iterations independent of M Systematic convergence propertie No localization in real space. Empty regions have to be filled with plane waves. Only supercell approach. O(N) scaling hard to achieve Resolution cannot adaptively be refined around the nucleus. All electron calculations impossible, hard pseudopotentials difficult.
73 : Planes Waves Operator approach: no need to set up Hamiltonian matrix FFT allows for quasi linear M log(m) scaling where M is number of plane waves Localization in Fourier space allows for efficient preconditioning techniques. Hence number of iterations independent of M Systematic convergence propertie No localization in real space. Empty regions have to be filled with plane waves. Only supercell approach. O(N) scaling hard to achieve Resolution cannot adaptively be refined around the nucleus. All electron calculations impossible, hard pseudopotentials difficult.
74 : Planes Waves self-consistent (SCF) Harris-Foulkes functional relativistic pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W non-relativistic [ N 3 scaling O(N) methods periodic non-periodic non-collinear collinear +V(r) non-spin polarized spin polarized LDA,A ] +µ xc (r) ψ k i = ε k i ψ k i atomic orbitals plane waves augmented real space LDA+U Slater numerical finite difference Wavelet
75 : List of Codes ABINIT Louvain-la-Neuve PAW, linear response, W, many features, good // CPMD Zurich, Lugano USP, good //, QM/MM with romacs PWSCF Italy USP, linear response, // VASP Vienna robust, fast, hybrid functional, W, //
76 : Flowchart { } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized
77 : Flowchart { inv FFT } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2
78 : Flowchart { inv FFT } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2 FFT V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective
79 : Flowchart { inv FFT } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2 FFT V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j Kinetic Term
80 : Flowchart { inv FFT } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2 FFT V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j FFT Kinetic Term δc j = E total c j () + l Λ jl c l Λ jl =< ψ j H ψ l >
81 : Flowchart { inv FFT } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2 FFT V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j FFT Kinetic Term δc j = E total c new,j c j () + l Λ jl c l = c j + h stepδc j Λ jl =< ψ j H ψ l > Steepest Descent, Conjugate radient, Direct Inversion of the Iterative Subspace
82 : Flowchart { inv FFT } ψ j = c j ei.r Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2 FFT V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j FFT Kinetic Term δc j = E total Stop when δc j small c new,j c j () + l Λ jl c l = c j + h stepδc j Λ jl =< ψ j H ψ l > Steepest Descent, Conjugate radient, Direct Inversion of the Iterative Subspace
83 Direct Minimization: List of Codes Plane Waves aussian CPMD Zurich, Lugano CP2K basis set ONETEP PARATEC Big
84 Scheme: Flowchart { ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized
85 Scheme: Flowchart { inv FFT ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j Basis (2 3 N ): 0,...,2
86 Scheme: Flowchart { inv FFT ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j FFT Basis (2 3 N ): 0,...,2 V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective
87 Scheme: Flowchart { inv FFT ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j FFT Basis (2 3 N ): 0,...,2 V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j Kinetic Term
88 Scheme: Flowchart { inv FFT ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j FFT Basis (2 3 N ): 0,...,2 V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j FFT Kinetic Term [ V eff (r) ] ψ j = ε j ψ j
89 Scheme: Flowchart { inv FFT ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j FFT Basis (2 3 N ): 0,...,2 V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j FFT Kinetic Term [ V eff (r) ] ψ j = ε j ψ j ρ new (r) = αρ out + (1 α)ρ in ρ out (r) = f j ψ j (r) 2 j Mixing Density: Anderson, Broyden,DIIS
90 Scheme: Flowchart { inv FFT ψ j = c j }, ei.r {f j } Basis (N ): 0,..., Orthonormalized ρ(r) = f j ψ j (r) 2 j FFT Basis (2 3 N ): 0,...,2 V 2 xc [ρ(r)] V H () = ρ() V NL ({ψ j }) V effective c j FFT Kinetic Term [ V eff (r) ] ψ j = ε j ψ j Stop when ρ new ρ in small ρ new (r) = αρ out + (1 α)ρ in ρ out (r) = f j ψ j (r) 2 j Mixing Density: Anderson, Broyden,DIIS
91 Iterative diagonalization Standard diagonalization routines, such as found in software packages like LAPACK, are designed to find all the eigenvalues and eigenvectors of a matrix. Scale like n 3, where n is the dimension of the matrix We need to calculate only a small number of eigenvalues and eigenvectors of matrices of a very large dimension, another type of diagonalization, namely iterative diagonalization is used: Lanczos, Davidson, RM-DIIS, Conjugate radient, LOBPC,...
92 Diagonalization scheme: List of Codes ABINIT Louvain-la-Neuve Diagonalization: Conjugate-radient, LOBPC Mixing: Broyden, Pulay, Anderson, others based on dielectric constant CPMD Zurich, Lugano Diagonalization: Lanczos, Davidson Mixing: Broyden, Pulay, Anderson PWSCF Italy VASP Vienna Diagonalization: RM-DIIS Big
93 Atom centered are the most popular basis functions for electronic structure calculations in the quantum chemistry community. (x X i ) l 1 (y Y i ) l 2 (z Z i ) l 3 exp ( α((x X i ) 2 + (y Y i ) 2 + (z Z i ) 2 ) ) The different functions associated to R i = (X i,y i,z i ) differ by their extent, specified by α and by the exponents l 1, l 2, l 3. Because the have qualitatively the shape of atomic orbitals one needs for moderate accuracy only a small number of basis function per atom (typically 10 to 20 per atom).
94 : localization makes them well suited for molecules and other open structures. Kinetic and overlap matrix elements can be calculated analytically. ood description of core electrons. Small number of basis functions necessary for moderate accuracy No systematic convergence. Over-completeness before convergence. This leads also to problems with O(N). Many different recipes for generating basis sets. Pulay forces have been implemented.
95 : localization makes them well suited for molecules and other open structures. Kinetic and overlap matrix elements can be calculated analytically. ood description of core electrons. Small number of basis functions necessary for moderate accuracy No systematic convergence. Over-completeness before convergence. This leads also to problems with O(N). Many different recipes for generating basis sets. Pulay forces have been implemented.
96 : self-consistent (SCF) Harris-Foulkes functional relativistic pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W non-relativistic [ N 3 scaling O(N) methods periodic non-periodic non-collinear collinear +V(r) non-spin polarized spin polarized LDA,A LDA+U ] +µ xc (r) ψ= k i ε k i S ψ k i atomic orbitals plane waves augmented real space Slater numerical finite difference Wavelet
97 basis set: wavelets (Big) self-consistent (SCF) Harris-Foulkes functional relativistic pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA W method non-relativistic [ N 3 scaling O(N) methods periodic non-periodic non-collinear collinear +V(r) +V xc (r) non-spin polarized spin polarized LDA,A ] ψ k i atomic orbitals plane waves augmented real space LDA+U = ε k i ψ k i Slater numerical finite difference Wavelet
98 : List of Codes ONETEP sinc function, O(N) method, inside Material Studio, good // Big wavelets, good // PARATEC finite difference, good // OCTOPUS finite difference, TD- PAW rid-based projector-augmented wave method, finite difference
99 : List of Codes ONETEP sinc function, O(N) method, inside Material Studio, good // Big wavelets, good // PARATEC finite difference, good // OCTOPUS finite difference, TD- PAW rid-based projector-augmented wave method, finite difference
100 : List of Codes ONETEP sinc function, O(N) method, inside Material Studio, good // Big wavelets, good // PARATEC finite difference, good // OCTOPUS finite difference, TD- PAW rid-based projector-augmented wave method, finite difference
101 : List of Codes ONETEP sinc function, O(N) method, inside Material Studio, good // Big wavelets, good // PARATEC finite difference, good // OCTOPUS finite difference, TD- PAW rid-based projector-augmented wave method, finite difference
102 : List of Codes ONETEP sinc function, O(N) method, inside Material Studio, good // Big wavelets, good // PARATEC finite difference, good // OCTOPUS finite difference, TD- PAW rid-based projector-augmented wave method, finite difference
103
104 Functionals: Mixed between physics and chemistry approach QM/MM (Quantum Mechanics/Molecular Modeling) Multi-scale approach (more than 1000 atoms feasible for a better parametrization) Order N (real space basis set as wavelets) Numerical experience (high performance computing): One-day simulation Better exploration of atomic configurations Molecular Dynamics (4s per step for 32 water molecules)
Poisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
More informationBigDFT tutorial. BigDFT tutorial FIRST YARMOUK SCHOOL. Thierry Deutsch. November, L_Sim - CEA Grenoble. Introduction DFT.
Big tutorial FIRST YARMOUK SCHOOL Big tutorial Thierry Deutsch L_Sim - CEA Grenoble November, 2. 2010 Outline 1 2 Density Functional Theory (quick view) 3 4 A basis for nanosciences: the Big project STREP
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 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 informationIntroduction to Electronic Structure Calculations with BigDFT
Multiscale Modelling Methods for Applications in Materials Science CECAM JÜLICH, GERMANY run to Electronic Structure Calculations with Thierry Deutsch, Damien Caliste, Luigi Genovese L_Sim - CEA Grenoble
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 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 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 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 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 informationMolecular Mechanics: The Ab Initio Foundation
Molecular Mechanics: The Ab Initio Foundation Ju Li GEM4 Summer School 2006 Cell and Molecular Mechanics in BioMedicine August 7 18, 2006, MIT, Cambridge, MA, USA 2 Outline Why are electrons quantum? Born-Oppenheimer
More informationGEM4 Summer School OpenCourseWare
GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Molecular Mechanics by Ju Li. Given August 9, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use
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 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 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 informationVASP: running on HPC resources. University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
VASP: running on HPC resources University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria The Many-Body Schrödinger equation 0 @ 1 2 X i i + X i Ĥ (r 1,...,r
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 informationKey concepts in Density Functional Theory
From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem
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 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 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 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 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 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 informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationWavelets for density functional calculations: Four families and three. applications
Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Interpolating
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 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 informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
More informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
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 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 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 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 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 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 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 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 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 informationWhy 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 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 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 informationPerformance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26
Performance of hybrid density functional methods, screened exchange and EXX-OEP methods in the PAW approach Georg Kresse, J Paier, R Hirschl, M Marsmann Institut für Materialphysik and Centre for Computational
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 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 informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
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 informationIntroduction to Density Functional Theory
1 Introduction to Density Functional Theory 21 February 2011; V172 P.Ravindran, FME-course on Ab initio Modelling of solar cell Materials 21 February 2011 Introduction to DFT 2 3 4 Ab initio Computational
More informationSelf Consistent Cycle
Self Consistent Cycle Step 0 : defining your system namelist SYSTEM How to specify the System All periodic systems can be specified by a Bravais Lattice and and atomic basis How to specify the Bravais
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 informationSet the initial conditions r i. Update neighborlist. Get new forces F i
v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationThe Plane-wave Pseudopotential Method
The Plane-wave Pseudopotential Method k(r) = X G c k,g e i(g+k) r Chris J Pickard Electrons in a Solid Nearly Free Electrons Nearly Free Electrons Nearly Free Electrons Electronic Structures Methods Empirical
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 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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationOrbital Density Dependent Functionals
Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)
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 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 informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
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 informationKey concepts in Density Functional Theory (II)
Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université
More informationAccuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation
Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Multi-wavelets: High accuracy
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 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 informationIntroduction to First-Principles Method
Joint ICTP/CAS/IAEA School & Workshop on Plasma-Materials Interaction in Fusion Devices, July 18-22, 2016, Hefei Introduction to First-Principles Method by Guang-Hong LU ( 吕广宏 ) Beihang University Computer
More informationMODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
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 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 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 informationThe Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus
Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,
More information2.1 Introduction: The many-body problem
Chapter 2 Smeagol: Density Functional Theory and NEGF s 2.1 Introduction: The many-body problem In solid state physics one is interested in systems comprising many atoms, and consequently many electrons.
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 informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationBefore we start: Important setup of your Computer
Before we start: Important setup of your Computer change directory: cd /afs/ictp/public/shared/smr2475./setup-config.sh logout login again 1 st Tutorial: The Basics of DFT Lydia Nemec and Oliver T. Hofmann
More information1. Thomas-Fermi method
1. Thomas-Fermi method We consider a system of N electrons in a stationary state, that would obey the stationary Schrödinger equation: h i m + 1 v(r i,r j ) Ψ(r 1,...,r N ) = E i Ψ(r 1,...,r N ). (1.1)
More informationFirst-principles Molecular Dynamics Simulations
First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationComputational Nanoscience
Computational Nanoscience Applications for Molecules, Clusters, and Solids KALMÄN VARGA AND JOSEPH A. DRISCOLL Vanderbilt University, Tennessee Щ CAMBRIDGE HP UNIVERSITY PRESS Preface Part I One-dimensional
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 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 informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
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 informationOptimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics
Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Sangeeta Sharma 1,2, J. K. Dewhurst 3, C. Ambrosch-Draxl 4, S. Pittalis 2, S. Kurth 2, N. Helbig 2, S. Shallcross
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationOn the adaptive finite element analysis of the Kohn-Sham equations
On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015 Denis Davydov, LTM, Erlangen, Germany College
More informationAtomic Structure and Atomic Spectra
Atomic Structure and Atomic Spectra Atomic Structure: Hydrogenic Atom Reading: Atkins, Ch. 10 (7 판 Ch. 13) The principles of quantum mechanics internal structure of atoms 1. Hydrogenic atom: one electron
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 informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
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 informationAb-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials
ISSN 0974-9373 Vol. 15 No.3 (2011) Journal of International Academy of Physical Sciences pp. 337-344 Ab-initio Electronic Structure Calculations of α, β and γ KNO 3 Energetic Materials Pradeep Jain and
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More information2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density
1. SIESTA DFT In a Nutshell Introduction to SIESTA Atomic Orbitals Capabilities Resources 2. TranSIESTA Transport in the Nanoscale - motivation Boundary Conditions: Open systems Greens functions and charge
More informationThe Nuclear Many-Body Problem
The Nuclear Many-Body Problem relativistic heavy ions vacuum electron scattering quarks gluons radioactive beams heavy few nuclei body quark-gluon soup QCD nucleon QCD few body systems many body systems
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 informationPlane waves and pseudopotentials (Part II)
Bristol 2007 Plane waves and pseudopotentials (Part II) X. Gonze Université Catholique de Louvain, Louvain-la-neuve, Belgium Plane Waves - Pseudopotentials, Bristol March 2007 1 In view of practical calculations,
More informationSummary of the new Modelling Vocabulary
Summary of the new Modelling Vocabulary These two pages attempts to summarise in a concise manner the Modelling Vocabulary. What are Models? What are Simulations? Materials Models consist of Physics or
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More information