BigDFT tutorial. BigDFT tutorial FIRST YARMOUK SCHOOL. Thierry Deutsch. November, L_Sim - CEA Grenoble. Introduction DFT.
|
|
- Charleen Wilson
- 6 years ago
- Views:
Transcription
1 Big tutorial FIRST YARMOUK SCHOOL Big tutorial Thierry Deutsch L_Sim - CEA Grenoble November,
2 Outline 1 2 Density Functional Theory (quick view) 3 4
3 A basis for nanosciences: the Big project STREP European project: Big( ) Four partners, 15 contributors: CEA-INAC Grenoble, U. Basel, U. Louvain-la-Neuve, U. Kiel Aim: To develop an ab-initio code based on Daubechies Wavelets, to be integrated in ABINIT, distributed freely (GNU-GPL license) L. Genovese, A. Neelov, S. Goedecker, T. Deutsch, et al., Daubechies wavelets as a basis set for density functional pseudopotential calculations, J. Chem. Phys. 129, (2008)
4 Big version 1.3: capabilities Free, surface and periodic boundary conditions Local geometry optimizations (with constraints) Global geometry optimization Saddle point searches Vibrations External electric fields Born-Oppenheimer MD Unoccupied KS orbitals Collinear and Non-collinear magnetism All XC functionals of the ABINIT package Empirical van der Waals interactions Also available from within the ABINIT package
5 Big version 1.4 (June 2010) Better performance for periodic systems (convolutions) Exact exchange Hybrid functionals (libxc, G2-1 tests) Molecular Dynamics from ABINIT K points (periodic) Symmetries from ABINIT Preliminary version for XANES version with S_ Big/ABINIT training workshop: June 30 Juny 1
6 Spirit of Big Simple for developers: no parser, use input files with mandatory parameters Few input files Optional input files as input.kpt (for periodic systems) or input.perf Few input parameters ex. not parameters for parallel calculations
7 Big: Executables bigdft Single point calculation, geometry optimisation, molecular dynamics Execute simply: bigdft memguess number_of_cores Check input files and determine the memory amount frequencies Vibration properties abscalc XANES calculations bart Big + ART (global minima search) NEB Determination of saddle point
8 Big: input files posinp.xyz Atomic positions + boundary conditions input.dft Everything about calculations psppar.n file with the name of atom (ABINIT format) Optional files input.geopt Geometry optimisation or molecular dynamics parameters input.freq Frequencies calculations input.perf Options for performances input.kpt K-points mesh for periodic systems
9 Atomic positions: posinp.xyz or posinp.ascii units reduced, angstroem, bohr or atomic Boundary conditions periodic, surface 8 reduced periodic Si Si Si Si Si Si Si Si
10 Atomic positions: posinp.xyz or posinp.ascii units reduced, angstroem, bohr or atomic Boundary conditions periodic, surface 3 angstroem surface O H H
11 Atomic positions: posinp.xyz or posinp.ascii units reduced, angstroem, bohr or atomic Boundary conditions periodic, surface 5 atomic Si H H H H
12 Kohn-Sham formalism H-K theorem: E is an unknown functional of the density E = E[ρ] Density Functional Theory Kohn-Sham 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/2 i=1 Ψ i (r) 2 Ψ i (r)dr + 1 ρ(r)v H (r)dr 2 + E xc [ρ(r)] + V ext (r)ρ(r)dr N/2 ρ(r) = 2 i=1 Ψ i (r)ψ i (r) 2 V H (r) = 4πρ(r)
13 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, (Kohn-Sham) Ingredients An XC potential, functional of the density several approximations exists (LDA,GGA,... ) 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...
14 Exchange-Correlation Energy (input.dft) E xc [ρ] = K exact [ρ] K KS [ρ] + V e e exact [ρ] V Hartree [ρ] K KS [ρ]: Kinetic energy for a non-interaction electron gas Kohn-Sham formalism is exact (mapping of an interacting electron system into a non-interacting system) Many types from ABINIT: Use same codes (ixc: LDA=1, PBE=11, see manual) libxc: code for exchange code for correlation ( see manual) Hybrid functionals from libxc
15 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, (Kohn-Sham) Ingredients An XC potential, functional of the density several approximations exists (LDA,GGA,... ) 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...
16 Gaussian type separable s (HGH) 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 ) = 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 ) 6 } p l 1(r) = 2 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 r l+ 3 2 l r l+ 7 2 l
17 The nonlocal part of the pseudopotential H sep Ψ = P >< P Ψ > If the expansion coefficients c J of P in a Daubechies basis are known then ( )( ) < P Ψ >= c I φ(x I) Ψ J φ(x J) dx = c J Ψ J I J J Because a Gaussian is separable and a 3-dim scaling basis is a product of 1-dim scaling functions/wavelets, c i,j,k can easily be calculated with machine precision. c i,j,k = φ k i (x)φ k j (y)φ k k(z)exp( (x 2 + y 2 + z 2 )) dx dy dz = φ k i (x)exp( x 2 ) dx φ k i (y)exp( y 2 ) dy φ k i (z)exp( z 2 ) dz
18 file: psppar.n In the same directory as the calculation. Large library: Accurate and transferable Goedecker pseudopotential for N zatom,zion,pspdat pspcod,pspxc,lmax,lloc,mmax,r2well rloc nloc c1 c2 2 nnonloc rs ns hs rp np hp kp11 Save memory: projectors calculated on-the-fly (cost < 1%, experimental option)
19 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, (Kohn-Sham) Ingredients An XC potential, functional of the density several approximations exists (LDA,GGA,... ) 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...
20 s in real space: Wavelet basis sets Properties of wavelet basis sets: localized both in real and in Fourier space allow for adaptivity (for core electrons) are a systematic basis set φ(x) ψ(x) x A basis set both adaptive and systematic, real space based
21 A brief description of wavelet theory Two kind of basis functions A Multi-Resolution real space basis The functions can be classified following the resolution level they span. Scaling Functions The functions of low resolution level are a linear combination of high-resolution functions = + m φ(x) = h j φ(2x j) j= m Centered on a resolution-dependent grid: φ j = φ 0 (x j).
22 A brief description of wavelet theory Wavelets They contain the DoF needed to complete the information which is lacking due to the coarseness of the resolution. = φ(2x) = m j= m h j φ(x j) + m j= m g j ψ(x j) Increase the resolution without modifying grid space SF + W = same DoF of SF of higher resolution ψ(x) = m j= m g j φ(2x j) All functions have compact support, centered on grid points.
23 Adaptivity of the mesh Adaptivity Resolution can be refined following the grid point. The grid is divided in: Data can thus be stored in compressed form. Low resolution pts (SF, 1 DoF) High resolution pts (SF + W, 8 DoF) Points of different resolution belong to the same grid. Localization property Empty regions must not be filled with basis functions. Optimal for big inhomogeneous systems, O(N) approach.
24 Adaptivity of the mesh Atomic positions (H 2 O), hgrid
25 Adaptivity of the mesh Fine grid (high resolution, frmult)
26 Adaptivity of the mesh Coarse grid (low resolution, crmult)
27 input.dft: hgrid, crmult, frmult hgrid Step grid: twice as bigger than the lowest gaussian coefficient of pseudopotential Roughly 0.45 (only orthorhombic mesh) crmult Extension of the coarse resolution (factor 6.): depends on the extension of the HOMO for the atom frmult Extension of the fine resolution (factor 8.): depends on the pseudopotential coefficient In practice, we use only hgrid from 0.3 to 0.55.
28 No integration error Orthogonality, scaling relation dx φ k (x)φ j (x) = δ kj φ(x) = 1 2 m j= m h j φ(2x j) The hamiltonian-related quantities can be calculated up to machine precision in the given basis. The accuracy is only limited by the basis set (O ( h 14 grid) ) Exact evaluation of kinetic energy Obtained by convolution with filters: f(x) = c l φ l (x), 2 f(x) = c l φ l (x), l l c l = c j a l j, a l φ 0 (x) 2 xφ l (x), j
29 Separability in 3D The 3-dim scaling basis is a tensor product decomposition of 1-dim Scaling Functions/ Wavelets. φ e x,e y,e z j x,j y,j z (x,y,z) = φ e x j x (x)φ e y j y (y)φ e z j z (z) With (j x,j y,j z ) the node coordinates, φ (0) j and φ (1) j the SF and the W respectively. Gaussians and wavelets The separability of the basis allows us to save computational time when performing scalar products with separable functions (e.g. gaussians): Initial wavefunctions (input guess) Poisson solver Non-local pseudopotentials
30 Wavelet families used in Big code Daubechies f(x) = l c l φ l (x) Orthogonal c l = dx φ l (x)f(x) LEAST ASYMMETRIC DAUBECHIES-16 scaling function wavelet Interpolating f(x) = j f j ϕ j (x) Dual to dirac deltas f j = f(j) 1 INTERPOLATING (DESLAURIERS-DUBUC)-8 scaling function wavelet Used for wavefunctions, scalar products Used for charge density, function products Magic Filter method (A.Neelov, S. Goedecker) The passage between the two basis sets can be performed without losing accuracy
31 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, (Kohn-Sham) Ingredients An XC potential, functional of the density several approximations exists (LDA,GGA,... ) 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...
32 Kohn-Sham 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) = 2 f i ψ i (r) 2 i
33 Kohn-Sham 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
34 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) = 2 i f i ψ i (r) 2
35 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) = 2 i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 )
36 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) = 2 i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 )
37 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) = 2 i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 )
38 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) = 2 i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 )
39 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) = 2 i f i ψ i (r) 2 Poisson Equation: V Hartree = ρ (Laplacian: = 2 x y z 2 )
40 Direct Minimisation: Flowchart { ψ i = a c i aφ a } Basis: adaptive mesh (0,...,N φ ) Orthonormalized
41 Direct Minimisation: Flowchart { ψ i = a c i aφ a } inv FWT + Magic Filter Basis: adaptive mesh (0,...,N φ ) Orthonormalized occ ρ(r) = i ψ i (r) 2 Fine non-adaptive mesh: < 8N φ
42 Direct Minimisation: Flowchart { ψ i = a c i aφ a } inv FWT + Magic Filter Basis: adaptive mesh (0,...,N φ ) Orthonormalized occ ρ(r) = i ψ i (r) 2 Fine non-adaptive mesh: < 8N φ Poisson solver V H (r) = V xc [ρ(r)] G(.)ρ V NL ({ψ i }) V effective
43 Direct Minimisation: Flowchart { ψ i = a c i aφ a } inv FWT + Magic Filter Basis: adaptive mesh (0,...,N φ ) Orthonormalized occ ρ(r) = i ψ i (r) 2 Fine non-adaptive mesh: < 8N φ Poisson solver V H (r) = V xc [ρ(r)] G(.)ρ V NL ({ψ i }) V effective Kinetic Term
44 Direct Minimisation: Flowchart { ψ i = a c i aφ a } inv FWT + Magic Filter Basis: adaptive mesh (0,...,N φ ) Orthonormalized FWT occ ρ(r) = i ψ i (r) 2 Fine non-adaptive mesh: < 8N φ Poisson solver V H (r) = V xc [ρ(r)] G(.)ρ V NL ({ψ i }) FWT δca i = E total c i (a) + Λ ij ca j j V effective Kinetic Term Λ ij =< ψ i H ψ j >
45 Direct Minimisation: Flowchart { ψ i = a c i aφ a } inv FWT + Magic Filter Basis: adaptive mesh (0,...,N φ ) Orthonormalized FWT occ ρ(r) = i ψ i (r) 2 Fine non-adaptive mesh: < 8N φ Poisson solver V H (r) = V xc [ρ(r)] G(.)ρ V NL ({ψ i }) δca i = E total c i (a) + Λ ij ca j j c new,i a FWT preconditioning = c i a + h step δc i a V effective Kinetic Term Λ ij =< ψ i H ψ j > Steepest Descent, DIIS
46 Direct Minimisation: Flowchart { ψ i = a c i aφ a } inv FWT + Magic Filter Basis: adaptive mesh (0,...,N φ ) Orthonormalized FWT occ ρ(r) = i ψ i (r) 2 Fine non-adaptive mesh: < 8N φ Poisson solver V H (r) = V xc [ρ(r)] G(.)ρ V NL ({ψ i }) Stop when δc i a small δca i = E total c i (a) + Λ ij ca j j c new,i a FWT preconditioning = c i a + h step δc i a V effective Kinetic Term Λ ij =< ψ i H ψ j > Steepest Descent, DIIS
47 input.dft: idsx DIIS history Use DIIS algorithm Direct Invesion of Iterative Subspace Use idsx=6 previous iterations Memory-consuming (decrease if too big=) Steepest descent if idsx=0 itermax=50 Specify maximum number of iterations Typically 15 iterations to converge gnrm_cv=1.e-4 convergence criterion Residue δψ 2 < gnrm_cv
48 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
49 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
50 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
51 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
52 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
53 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
54 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
55 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
56 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
57 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
58 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
59 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
60 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
61 Description of the file input.dft hx, hy, hz crmult, frmult, coarse and fine radius 1 ixc: XC parameter (LDA=1,PBE=11) ncharge: charge of the system, Electric field 1 0 nspin=1 non-spin polarization, total magnetic moment 1.E-04 gnrm_cv: convergence criterion gradient itermax, nrepmax 7 6 # CG (preconditioning), length of diis history 0 dispersion correction functional (0, 1, 2, 3) write "CUDA" to use (.config file needed) 0 F 0 InputPsiId, output_wf, output_grid length of the tail, # tail CG iterations davidson treatment, # virtual orbitals, # plotted orbitals 2 verbosity of the output 0=low, 2=high T disable the symmetry detection
62 s: wavelets self-consistent () Harris-Foulkes functional relativistic pseudopotential norm-conserving all electrons PAW Hybrid functionals beyond LDA GW non-relativistic [ N 3 scaling O(N) methods periodic non-periodic non-collinear collinear +V(r) non-spin polarized spin polarized LDA,GGA ] +µ xc (r) ψ k i = ε k i ψ k i atomic orbitals plane waves augmented real space LDA+U Gaussians Slater numerical finite difference Wavelet
63 Optimal for isolated systems 10 1 Test case: cinchonidine molecule (44 atoms) Plane waves Absolute energy precision (Ha) E c = 40 Ha Wavelets E c = 90 Ha h = 0.4bohr E c = 125 Ha h = 0.3bohr Number of degrees of freedom Allows a systematic approach for molecules Considerably faster than Plane Waves codes. 10 (5) times faster than ABINIT (CPMD) Charged systems can be treated explicitly with the same time
64 Systematic basis set Two parameters for tuning the basis The grid spacing hgrid The extension of the low resolution points crmult
65 Poisson Solver using interpolating scaling function The Hartree potential is calculated in the interpolating scaling function basis. Poisson solver with interpolating scaling functions From the density ρ( j) on an uniform grid it calculates: V H ( j) = d x ρ( x) x j Very fast and accurate, optimal parallelization Can be used independently from the code Integrated quantities (energies) are easy to extract Non-adaptive, needs data uncompression Explicitly free boundary conditions No need to subtract supercell interactions charged systems.
66 Poisson Solver for surface boundary conditions A Poisson solver for surface problems Based on a mixed reciprocal-direct space representation V px,p y (z) = 4π dz G( p ;z z )ρ px,p y (z), Can be applied both in real or reciprocal space codes No supercell or screening functions More precise than other existing approaches Example of the plane capacitor: Periodic Hockney Our approach V V V y y y x x x
67 input.geopt: Relaxation or molecular dynamics DIIS 500 ncount_cluster: max steps during geometry relaxation 1.d0 frac_fluct: geometry optimization stops if force norm less than frac_fluct of noise 0.0d0 randdis: random amplitude for atoms 4.d0 5 betax: steepest descent step size, diis history Different algorithms: DIIS: Direct Inversion in the iterative subspace VSSD: Variable Size Steepest Descent LBFGS: Limited BFGD AB6MD: Use molecular dynamics routines from ABINIT 6
68 input.geopt: Relaxation or molecular dynamics DIIS 500 ncount_cluster: max steps during geometry relaxation 1.d0 frac_fluct: geometry optimization stops if force norm less than frac_fluct of noise 0.0d0 randdis: random amplitude for atoms 4.d0 5 betax: steepest descent step size, diis history Different algorithms: DIIS: Direct Inversion in the iterative subspace VSSD: Variable Size Steepest Descent LBFGS: Limited BFGD AB6MD: Use molecular dynamics routines from ABINIT 6
69 input.geopt: Relaxation or molecular dynamics DIIS 500 ncount_cluster: max steps during geometry relaxation 1.d0 frac_fluct: geometry optimization stops if force norm less than frac_fluct of noise 0.0d0 randdis: random amplitude for atoms 4.d0 5 betax: steepest descent step size, diis history Different algorithms: DIIS: Direct Inversion in the iterative subspace VSSD: Variable Size Steepest Descent LBFGS: Limited BFGD AB6MD: Use molecular dynamics routines from ABINIT 6
70 input.geopt: Relaxation or molecular dynamics DIIS 500 ncount_cluster: max steps during geometry relaxation 1.d0 frac_fluct: geometry optimization stops if force norm less than frac_fluct of noise 0.0d0 randdis: random amplitude for atoms 4.d0 5 betax: steepest descent step size, diis history Different algorithms: DIIS: Direct Inversion in the iterative subspace VSSD: Variable Size Steepest Descent LBFGS: Limited BFGD AB6MD: Use molecular dynamics routines from ABINIT 6
71 input.geopt: Relaxation or molecular dynamics DIIS 500 ncount_cluster: max steps during geometry relaxation 1.d0 frac_fluct: geometry optimization stops if force norm less than frac_fluct of noise 0.0d0 randdis: random amplitude for atoms 4.d0 5 betax: steepest descent step size, diis history Different algorithms: DIIS: Direct Inversion in the iterative subspace VSSD: Variable Size Steepest Descent LBFGS: Limited BFGD AB6MD: Use molecular dynamics routines from ABINIT 6
72 Operations performed Different numerical operations performed: For each wavefunction (Hamiltonian application) Between wavefunctions (Linear algebra) Comput. operations Convolutions with short filters BLAS routines FFT (Poisson Solver) Interpolating Daubechies
73 Massively parallel Two kinds of parallelisation By orbitals (Hamiltonian application, preconditioning) By components (overlap matrices, orthogonalisation) A few (but big) packets of data More demanding in bandwidth than in latency No need of fast network Optimal speedup (eff. 85%), also for big systems Cubic scaling code For systems bigger than 500 atomes (1500 orbitals) : orthonormalisation operation is predominant (N 3 )
74 High Performance Computing Localisation & Orthogonality Data locality Principal code operations can be intensively optimised Optimal for application on supercomputers Little communication, big packets of data No need of fast network Optimal speedup Efficiency of the order of 90%, up to thousands of processors Efficiency (%) at 65 at 173 at 257 at 1025 at (with orbitals/proc) Number of processors Data repartition optimal for material accelerators () Graphic Processing Units can be used to speed up the computation
75 Performance in parallel (Big) Distribution per orbitals (wavelets) Nothing to specify in input files
76 -ported operations in Big (double precision) 20 Convolutions (CUDA rewritten) speedups between and 20 can be 6 obtained for different 4 sections 2 70 speedup (Double prec.) locden locham precond Wavefunction size (MB) speedup (Double prec.) DGEMM DSYRK Number of orbitals Linear algebra (CUBLAS library) The interfacing with CUBLAS is immediate, with considerable speedups
77 Data repartition on a node The S_ approach S_ library manages resource within the node S_ MPI 0 MPI 1 MPI 2 MPI 3 DATA DATA DATA DATA
78 Big code on Hybrid architectures Big code can run on hybrid CPU/ supercomputers In multi- environments, double precision calculations No Hot-spot operations Different code sections can be ported on up to 20x speedup for some operations, 7x for the full parallel code (under improvements) Percent Seconds (log. scale) CPU code Hybrid code (rel.) Number of cores 1 Time (sec) Comm Other PureCPU precond locham locden LinAlg
79 CUDA + S_ Put "CUDA" in input.dft Use a.config file (affinity between CPU and ) USE_SHARED=0 Set if the repartition is shared (1) or static (0) MPI_TASKS_PER_NODE=8 Define the number of mpi tasks per node NUM_=2 Define the number of device per node _CPUS_AFF_0=0,1,2,3 Define the affinity of the processes among the s. _CPUS_AFF_1=4,5,6,7 USE BLAS=1 Use CUDABLAS
80 OpenCL (without S_) Put "OCL" in input.dft No.config needed presently (should be removed also for CUDA) Fermi is better for concurrent access
81 Functionality: Mixed between physics and chemistry approach Order N (real space basis set as wavelets) QM/MM (Quantum Mechanics/Molecular Modeling) Multi-scale approach (more than 1000 atoms feasible for a better parametrization) Numerical experience: One-day simulation Better exploration of atomic configurations Molecular Dynamics (4s per step for 32 water molecules)
82 THE END
83 Future developments, Discussions
84 O(N) approach Versatile boundary conditions (surface, periodic, isolated) Localization regions Wannier functions Easily embedded T CPU Crossover point N
85 Embedding a system (QM/MM) QM/MM Charge defects in bulk materials Semi-infinite system as surface
86 Developments Embedding atoms, O(N) method, region of localisation Many-body: TD-, GW, Bethe-Salpeter PAW pseudopotentials (better accuracy, larger step grid) Mixed basis sets
Introduction 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 informationPoisson 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 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 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 informationDaubechies wavelets as a basis set for density functional pseudopotential calculations
Daubechies wavelets as a basis set for density functional pseudopotential calculations Luigi Genovese, Alexey Neelov, 2 Stefan Goedecker, 2 Thierry Deutsch, Seyed Alireza Ghasemi, 2 Alexander Willand,
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 informationELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS
FROM RESEARCH TO INDUSTRY 32 ème forum ORAP 10 octobre 2013 Maison de la Simulation, Saclay, France ELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS APPLICATION ON HPC, BLOCKING POINTS, Marc
More informationAdaptive and Localized Basis Functions for Linear Scaling, Large Systems and Complex QM Simulations
Max Conference on the Materials Design Ecosystem at the Exascale: High-Performance and High-Throughput Computing TRIESTE Adaptive and Localized Basis Functions for Linear Scaling, Large Systems and Complex
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 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 informationReinhold Schneider, MATHEON TU Berlin. Direct Minimization for effective single particle models (DFT)
Direct Minimization for effective single particle models (DFT) Partners joint work with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
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 informationab initio Electronic Structure Calculations
ab initio Electronic Structure Calculations New scalability frontiers using the BG/L Supercomputer C. Bekas, A. Curioni and W. Andreoni IBM, Zurich Research Laboratory Rueschlikon 8803, Switzerland ab
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 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 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 informationBlock Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
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 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 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 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 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 informationElectronic structure calculations with GPAW. Jussi Enkovaara CSC IT Center for Science, Finland
Electronic structure calculations with GPAW Jussi Enkovaara CSC IT Center for Science, Finland Basics of density-functional theory Density-functional theory Many-body Schrödinger equation Can be solved
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 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 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 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 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 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 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 informationOpportunities from Accurate and Efficient Density Functional Theory Calculations for Large Systems
Seminar CENTRE FOR PREDICTIVE MODELLING, WARWICK Opportunities from Accurate and Efficient Density Functional Theory Calculations for Large Systems Luigi Genovese L_Sim CEA Grenoble October 30, 2017 http://bigdft.org
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 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 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 informationDirect Minimization in Density Functional Theory
Direct Minimization in Density Functional Theory FOCM Hongkong 2008 Partners joint paper with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
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 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 informationCharge Analysis: Atoms in Molecules
Daubechies Wavelets in Electronic Structure Calculation: BigDFT Code Tutorial CECAM - GRENOBLE : Atoms in Molecules Ali Sadeghi Basel University 21 November 2011 An output of electronic structure calculations
More informationarxiv: v1 [physics.comp-ph] 22 Nov 2012
A Customized 3D GPU Poisson Solver for Free BCs Nazim Dugan a, Luigi Genovese b, Stefan Goedecker a, a Department of Physics, University of Basel, Klingelbergstr. 82, 4056 Basel, Switzerland b Laboratoire
More informationMinimization of the Kohn-Sham Energy with a Localized, Projected Search Direction
Minimization of the Kohn-Sham Energy with a Localized, Projected Search Direction Courant Institute of Mathematical Sciences, NYU Lawrence Berkeley National Lab 29 October 2009 Joint work with Michael
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 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 informationAll-electron density functional theory on Intel MIC: Elk
All-electron density functional theory on Intel MIC: Elk W. Scott Thornton, R.J. Harrison Abstract We present the results of the porting of the full potential linear augmented plane-wave solver, Elk [1],
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 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 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 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 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 informationIntroduction to Parallelism in CASTEP
to ism in CASTEP Stewart Clark Band University of Durham 21 September 2012 Solve for all the bands/electrons (Band-) Band CASTEP solves the Kohn-Sham equations for electrons in a periodic array of nuclei:
More informationMassively parallel electronic structure calculations with Python software. Jussi Enkovaara Software Engineering CSC the finnish IT center for science
Massively parallel electronic structure calculations with Python software Jussi Enkovaara Software Engineering CSC the finnish IT center for science GPAW Software package for electronic structure calculations
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 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 informationquantum ESPRESSO stands for Quantum open-source Package for Research in Electronic Structure, Simulation, and Optimization
The quantum ESPRESSO distribution The IOM-DEMOCRITOS center of Italian CNR is dedicated to atomistic simulations of materials, with a strong emphasis on the development of high-quality scientific software
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 informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part II Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationWeile Jia 1, Long Wang 1, Zongyan Cao 1, Jiyun Fu 1, Xuebin Chi 1, Weiguo Gao 2, Lin-Wang Wang 3
A plane wave pseudopotential density functional theory molecular dynamics code on multi-gpu machine - GPU Technology Conference, San Jose, May 17th, 2012 Weile Jia 1, Long Wang 1, Zongyan Cao 1, Jiyun
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 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 informationLocalized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy
Localized Optimization: Exploiting non-orthogonality to efficiently minimize the Kohn-Sham Energy Courant Institute of Mathematical Sciences, NYU Lawrence Berkeley National Lab 16 December 2009 Joint work
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 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 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 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 informationPreconditioned Eigenvalue Solvers for electronic structure calculations. Andrew V. Knyazev. Householder Symposium XVI May 26, 2005
1 Preconditioned Eigenvalue Solvers for electronic structure calculations Andrew V. Knyazev Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver Householder
More informationHYBRID FUNCTIONALS IN ABINIT: STATE OF THE ART AND PERSPECTIVES. F. Jollet 1, M. Torrent 1, B. Amadon 1, F. Arnardi 1, X.
HYBRID FUNCTIONALS IN ABINIT: STATE OF THE ART AND PERSPECTIVES F. Jollet 1, M. Torrent 1, B. Amadon 1, F. Arnardi 1, X. Gonze 1,2 1 CEA, DAM-DIF, Bruyères le Châtel, 91297 Arpajon Cedex 2 UCL,Louvain
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 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 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 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 informationSpeed-up of ATK compared to
What s new @ Speed-up of ATK 2008.10 compared to 2008.02 System Speed-up Memory reduction Azafulleroid (molecule, 97 atoms) 1.1 15% 6x6x6 MgO (bulk, 432 atoms, Gamma point) 3.5 38% 6x6x6 MgO (k-point sampling
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 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 informationThe Quantum ESPRESSO Software Distribution
The Quantum ESPRESSO Software Distribution The DEMOCRITOS center of Italian INFM is dedicated to atomistic simulations of materials, with a strong emphasis on the development of high-quality scientific
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationOPENATOM for GW calculations
OPENATOM for GW calculations by OPENATOM developers 1 Introduction The GW method is one of the most accurate ab initio methods for the prediction of electronic band structures. Despite its power, the GW
More informationEfficient implementation of the overlap operator on multi-gpus
Efficient implementation of the overlap operator on multi-gpus Andrei Alexandru Mike Lujan, Craig Pelissier, Ben Gamari, Frank Lee SAAHPC 2011 - University of Tennessee Outline Motivation Overlap operator
More informationCP2K: Past, Present, Future. Jürg Hutter Department of Chemistry, University of Zurich
CP2K: Past, Present, Future Jürg Hutter Department of Chemistry, University of Zurich Outline Past History of CP2K Development of features Present Quickstep DFT code Post-HF methods (RPA, MP2) Libraries
More informationBreaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points
Breaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points Michael Griebel Christian Rieger Peter Zaspel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität
More informationCP2K: Selected Developments
CP2K: Selected Developments Jürg Hutter Department of Chemistry University of Zurich Outline Introduction History and Performance Current and Future Developments Post-Hartree-Fock Methods GW Methods RPA
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 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 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 information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationOptimized energy landscape exploration for nanosciences using ab-initio based methods
Optimized energy landscape exploration for nanosciences using ab-initio based methods Eduardo Machado-Charry Nanosciences Foundation & Laboratoire de simulation atomistique (L Sim), SP2M, UMR-E CEA-Grenoble
More informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
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 informationThe Abinit project. Coding is based on modern software engineering principles
The Abinit project Abinit is a robust, full-featured electronic-structure code based on density functional theory, plane waves, and pseudopotentials. Abinit is copyrighted and distributed under the GNU
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 informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
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 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 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 information- Part 4 - Multicore and Manycore Technology: Chances and Challenges. Vincent Heuveline
- Part 4 - Multicore and Manycore Technology: Chances and Challenges Vincent Heuveline 1 Numerical Simulation of Tropical Cyclones Goal oriented adaptivity for tropical cyclones ~10⁴km ~1500km ~100km 2
More informationMaking Electronic Structure Theory work
Making Electronic Structure Theory work Ralf Gehrke Fritz-Haber-Institut der Max-Planck-Gesellschaft Berlin, 23 th June 2009 Hands-on Tutorial on Ab Initio Molecular Simulations: Toward a First-Principles
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 informationModified Becke-Johnson (mbj) exchange potential
Modified Becke-Johnson (mbj) exchange potential Hideyuki Jippo Fujitsu Laboratories LTD. 2015.12.21-22 OpenMX developer s meeting @ Kobe Overview: mbj potential The semilocal exchange potential adding
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 information