VASP: running on HPC resources. University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
|
|
- Randolph Ellis
- 5 years ago
- Views:
Transcription
1 VASP: running on HPC resources University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
2 The Many-Body Schrödinger equation 1 2 X i i + X i Ĥ (r 1,...,r N )=E (r 1,...,r N ) 1 V (r i )+ X 1 A (r 1,...,r N )=E (r 1,...,r N ) r i r j i6=j For instance, many-body WF storage demands are prohibitive: (r 1,...,r N ) (#grid points) N 5 electrons on a grid ~ 10 PetaBytes! A solution: map onto one-electron theory: (r 1,...,r N )! { 1 (r), 2 (r),..., N (r)}
3 Hohenberg-Kohn-Sham DFT Map onto one-electron theory: (r 1,...,r N )! { 1 (r), 2 (r),..., N (r)} (r 1,...,r N )= NY i i(r i ) Total energy is a functional of the density: E[ ] =T s [{ i [ ]}]+E H [ ]+E xc [ ] + E Z [ ]+U[Z] The density is computed using the one-electron orbitals: 1 2 (r) = NX i (r) 2 The one-electron orbitals are the solutions of the Kohn-Sham equation: + V Z (r)+v H [ ](r)+v xc [ ](r) i i(r) = i i (r) BUT: E xc [ ] =??? V xc [ ](r) =???
4 Exchange-Correlation E xc [ ] =??? V xc [ ](r) =??? Exchange-Correlation functionals are modeled on the uniform-electron-gas (UEG): The correlation energy (and potential) has been calculated by means of Monte- Carlo methods for a wide range of densities, and has been parametrized to yield a density functional. LDA: we simply pretend that an inhomogeneous electronic density locally behaves like a homogeneous electron gas. Many, many, many different functionals available: LDA, GGA, meta-gga, van-der-waals functionals, etc etc
5 An N-electron system: N = O(10 23 ) Hohenberg-Kohn-Sham DFT takes us a long way: (r 1,...,r N )! { 1 (r), 2 (r),..., N (r)} (#grid points) N N (#grid points) Nice for atoms and molecules, but in a realistic piece of solid state material N= O(10 23 )!
6 Translational invariance: Periodic Boundary Conditions Translational invariance implies: and nk(r + R) = nk (r)e ikr nk(r) =u nk (r)e ikr u nk (r + R) =u nk (r) All states are labeled by Bloch vector k and the band index n: The Bloch vector k is usually constrained to lie within the first Brillouin zone of the reciprocal space lattice. The band index n is of the order if the number of electrons per unit cell.
7 Reciprocal space & the first Brillouin zone L z 4pi/L z z b 2 b 1 b 2 b 1 x a 2 a 1 a 3 y x y x y b 3 b 3 A B C b 1 = 2 a 2 a 3 b 2 = 2 a 3 a 1 b 3 = 2 a 1 a 2 = a 1 a 2 a 3 a i b j =2 ij
8 Sampling the 1 st BZ The evaluation of many key quantities involves an integral over the 1 st BZ. For instance the charge density: (r) = 1 BZ X n Z BZ f nk nk (r) 2 dk We exploit the fact that the orbitals at Bloch vectors k that are close together are almost identical and approximate the integral over the 1 st BZ by a weighted sum over a discrete set of k-points: (r) = X nk w k f nk nk (r) 2 dk, Fazit: the intractable task of determining (r with N= ,...,r N ), has been reduced to calculating nk(r) at a discrete set of k-points in the 1 st BZ, for a number of bands that is of the order if the number of electrons in the unit cell.
9 The total energy E[, {R,Z}] =T s [{ nk [ ]}]+E H [, {R,Z}]+E xc [ ]+U({R,Z}) The kinetic energy The Hartree energy where T s [{ nk [ ]}] = X nk E H [, {R,Z}] = 1 2 w k f nk h nk 1 2 nki ZZ ez (r) ez (r 0 ) r r 0 dr 0 dr ez (r) = (r)+ X i Z i (r R i ) (r) = X nk w k f nk nk (r) 2 dk The Kohn-Sham equations V H [ ez ](r)+v xc [ ](r) nk(r) = nk nk (r) The Hartree potential V H [ ez ](r) = Z ez (r 0 ) r r 0 dr0
10 A plane wave basis set nk(r) =u nk (r)e ikr u nk (r + R) =u nk (r) All cell-periodic functions are expanded in plane waves (Fourier analysis): u nk (r) = 1 1/2 X G C Gnk e igr nk(r) = 1 1/2 X G C Gnk e i(g+k)r (r) = X G G e igr V (r) = X G V G e igr The basis set includes all plane waves for which C rnk = X G 1 2 G + k 2 <E cuto Transformation by means of FFT between real space and reciprocal space: C Gnk e igr FFT! CGnk = 1 N FFT X r C rnk e igr
11 real space FFT reciprocal space G cut τ N 1 0 b τ 1 x = n / N 1 1 C rnk = X G τ g = n 2 π / τ b 1 N/2+1 C Gnk e igr FFT! CGnk = 1 X C rnk e igr N FFT r N/2
12 The charge density G cut FFT ψ G 2 G cut ψ r ψ r FFT ρ G ρ r
13 The Self-Consistency-Cycle (cont.) Two sub-problems: Optimization of { n } Iterative Diagonalization e.g. RMM-DIIS or Blocked Davidson in Construction of Density Mixing e.g. Broyden mixer
14 The self-consistency-cycle A naïve algorithm: express the Hamilton matrix in a plane wave basis and diagonalize it: Self-consistency-cycle: H = hg Ĥ[ ] G0 i!diagonalize H! { i, i } i =1,..,N FFT 0! H 0! 0! 1 = f( 0, 0 )! H 1!... Iterate until: = 0 BUT: we do not need N FFT eigenvectors of the Hamiltonian (at a cost of O(N FFT3 )). Actually we only the N b lowest eigenstates of H, where N b is of the order of the number of electrons per unit cell (N b << N FFT ). Solution: use iterative matrix diagonalization techniques to find the N b lowest Eigenvector of the Hamiltonian: RMM-DIIS, blocked-davidson, etc.
15 Key ingredients: Subspace diagonalization and the Residual Rayleigh-Ritz: diagonalization of the N b x N b subspace with X m H nm B mk = X m app k S nm B mk H nm = h n Ĥ mi Snm = h n Ŝ mi yields N b eigenvectors with eigenvalues ε app. ki = X m B mk m i These eigenstates are the best approximation to the exact N b lowest eigenstates of H within the subspace spanned by the current orbitals. The Residual: R( n )i =(Ĥ appŝ) ni app = h n Ĥ ni h n Ŝ ni (its norm is measure for the error in the eigenvector)
16 Blocked-Davidson Take a subset of all bands: { n n =1,..,N} ) { 1 k k =1,..,n 1 } Extend this subset by adding the (preconditioned) residual vectors to the presently considered subspace: { 1 k/g 1 k = K(H app S) 1 k k =1,..,n 1 } Rayleigh-Ritz optimization ( sub-space rotation ) in the 2n 1 dimensional subspace to determine the n 1 lowest eigenvectors: diag{ 1 k/g 1 k} Extend subspace with the residuals of { 2 k} { 2 k k =1,..,n 1 } { 1 k/g 1 k/g 2 k = K(H app S) 2 k k =1,..,n 1 } Rayleigh-Ritz optimization ) { 3 k k =1,..,n 1 } Etc The optimized set replaces the original subset: { m k k =1,..,n 1 } { n n =1,..,n 1 } Continue with next subset: { k k 1 = n 1 +1,..,n 2 }, etc, After treating all bands: Rayleigh-Ritz optimization of { n n =1,..,N}
17 The action of the Hamiltonian The action H nk i Using the convention V (r) nk(r) hr G + ki = 1 1/2 ei(g+k)r! h G + k nk i = C Gnk Kinetic energy: hg + k hg + k V nk i = 1 N FFT X 1 2 nki = 1 2 G + k 2 C Gnk Local potential: Exchange-correlation: easily obtained in real space FFT to reciprocal space {V xc,r }! {V xc,g } Hartree potential: solve Poisson eq. in reciprocal space Add all contributions V G = V H,G + V xc,g + V ext,g FFT to real space {V G }! {V r } The action r V r C rnk e igr N NPLW V = V H [ ]+V xc [ ]+V ext V xc,r = V xc[ r ] V H,G = 4 G 2 G N FFT log N FFT
18 Solving the KS equations V ext (r)+v H (r)+v xc (r) nk(r) = nk nk (r) FFTs extensively used to evaluate H n i We actually use a mixed basis set (Projector-Augmented-Waves): n i = e ni + X i ( ii e ii)hep i e ni which involves projection of the pseudo-wave functions on local projection operators (DGEMM). One needs to keep the solutions (bands) at each k-point orthonormal: essentially done by a Choleski decomposition (LU) of the overlap matrix, followed by an inversion of U (LAPACK/scaLAPACK) and a transformation between the wave function (ZGEMM). Diagonalization of the Hamiltonian in the subspace of the current wave functions (LAPACK/scaLAPACK or ELPA), followed by a unitary transformation between the wave function (ZGEMM).
19 A typical workload Action: H ni FFT( n) N ln N 8 n N 2 ln N fftw V (r) n(r) N 8 n N 2 h p i ni N 8 i, n N 3 const. N 2 (real sp.) BLAS3 (DGEMM) LREAL= A Subspace rotation: H nm = h n H mi 8 n, m N 3 BLAS3 (ZGEMM) diag(h) N 3 (sca)lapack 0 ni = X m U nm mi 8 n N 3 BLAS3 (ZGEMM)
20 Scaling with system size (N) Self-Consistency Cycle (SCC): RMM-DIIS running on 8 Intel X5550 quadcore procs. (total: 32 x 2.66 GHz cores)
21 Distribution of work and data Distribute work and data over-orbitals Default NCORE = 1 (or equivalently: NPAR = #-of-mpi-ranks) KPAR = 1 The Kohn-Sham equation: V ext (r)+v H (r)+v xc (r) nk(r) = nk nk (r) Orbital index n
22 Distribution of work and data Distribute work and data over-orbitals Default NCORE = 1 (or equivalently: NPAR = #-of-mpi-ranks) KPAR = 1 Distribution work and data over-plane-waves NCORE = #-of-mpi-ranks (or equivalently: NPAR = 1) KPAR = 1
23 Distribution of work and data Combinations of over-orbitals and over-plane-wave distributions are allowed as well
24 Distribution of work and data Additionally work may be distributed over-k-points KPAR = n (n>1) m = (#-of-mpi-ranks / n) must be an integer Work is distributed in a round-robin fashion over groups of m MPI-ranks Data is duplicated! 1 + V ext (r)+v H (r)+v xc (r) nk(r) = nk nk (r) 2 Orbital index n, k-point index k
25 All-2-All communication
26 Parallel FFT: Ball-2-Cube x z y 4G 4G G FFT G 2G FFT 4G FFT 4G (A) (B) (C) (D) In-house parallel Ball-2-Cube FFT: Less 1D-FFTs (reduction: 1.76 ) BUT: communication from (B) (C) and (C) (D) For small to medium sized FFT grids a highly optimized 3D-FFT (e.g. from Intel s mkl) is equally fast
27 Hardware considerations Typical configuration: N interconnected nodes 2 packages/node M cores/package Distribution over-plane-waves : MPI-ranks that share an orbital should reside on the same node (better even on the same package). NCORE = n 2M (2M / n) should be an integer Typically: n = M or n = M/2 Distribution over-k-points : KPAR = n #-of-k-points (NKPTS) (#-of-mpi-ranks / n) should be an integer If memory allows: KPAR = NKPTS Default placement of MPI-ranks on the nodes/packages/cores depends on the particulars of the MPI implementation and its configuration!
28 For instance: 2 interconnected nodes 2 packages/node 4 cores/package Default placement of MPI-ranks on the nodes/packages/cores depends on the particulars of the MPI implementation and its configuration!
29 For instance: 2 interconnected nodes 2 packages/node 4 cores/package Good: Place subsequent MPI-ranks close together, i.e., first on subsequent cores of the same package, then moving on the second package of the same node, before starting to fill the next node.
30 For instance: 2 interconnected nodes 2 packages/node 4 cores/package Bad: Distribute subsequent MPI-ranks in a round-robin fashion over the packages.
31 For instance: 2 interconnected nodes 2 packages/node 4 cores/package Worse: Distribute subsequent MPI-ranks in a round-robin fashion over the nodes.
32 NCORE, KPAR, and the #-of-mpi-ranks N interconnected nodes 2 packages/node M cores/package 2NM MPI-ranks Distribution over-plane-waves : MPI-ranks that share an orbital should reside on the same node (better even on the same package). NCORE = n 2M (2M / n) should be an integer Typically: n = M or n = M/2 Distribution over-k-points : KPAR = n #-of-k-points (NKPTS) (#-of-mpi-ranks / n) should be an integer If memory allows: KPAR = NKPTS Distribution over-orbitals : The number of orbitals (NBANDS) is such that NBANDS / (2NM / NCORE / KPAR) is an integer number Increasing #-of-mpi-ranks may lead to an unnecessarily large NBANDS (i.e., adding empty orbitals) Some algorithms converge faster when each MPI rank owns (part of) a few orbitals (e.g. blocked-davidson) Generally speaking: having lots of MPI ranks and very little work/data per rank is never a good idea since (all-2-all) communication becomes unreasonably expensive.
33 Strong/Weak scaling (Si N )
34 Scaling under MPI (on a Cray XC-40) 5713 e- 2 k-points RMM-DIIS PBE KPAR=2 Courtesy of P. Saxe, Materials Design Inc. (and Cray).
35 Scaling under MPI (on a Cray XC-40) 864 e- 8 k-points HSE06 Davidson KPAR=8 Courtesy of P. Saxe, Materials Design Inc. (and Cray).
36 Large systems Use real-space PAW projectors (instead of reciprocal space projectors) h p i ni N 8 i, n N 3 const. N 2 (real sp.) BLAS3 (DGEMM) LREAL= A If you can limit k-point sampling to the Gamma (Γ k = 0) point: use the gamma-only version of VASP In the gamma-only version of VASP the orbitals are stored as real quantities in real-space: real-2-complex FFTs DGEMMs instead of ZGEMMs
37 Compilation Fast FFTs: the FFTs from Intel s mkl-library seem to be unbeatably fast scalapackfor large systems A compiler that effectively generates AVX2 instructions and libraries (e.g. BLAS) that are optimized for AVX2 (up to 20% performance gain)
38 The End Thank you!
Introduction 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 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 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 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 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 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 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 informationThe electronic structure of materials 2 - DFT
Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie
More 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 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 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 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 informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More informationIterative Diagonalization in Augmented Plane Wave Based Methods in Electronic Structure Calculations
* 0. Manuscript Click here to view linked References 1 1 1 1 0 1 0 1 0 1 0 1 0 1 Iterative Diagonalization in Augmented Plane Wave Based Methods in Electronic Structure Calculations P. Blaha a,, H. Hofstätter
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 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 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 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 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 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 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 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 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 informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
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 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 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 informationwoptic: Transport Properties with Wannier Functions and Adaptive k-integration
woptic: Transport Properties with Wannier Functions and Adaptive k-integration Elias Assmann Institute of Solid State Physics, Vienna University of Technology Split, 2013-09-29 About this presentation
More information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
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 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 informationSelf-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT
Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,
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 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 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 informationParallel Eigensolver Performance on High Performance Computers
Parallel Eigensolver Performance on High Performance Computers Andrew Sunderland Advanced Research Computing Group STFC Daresbury Laboratory CUG 2008 Helsinki 1 Summary (Briefly) Introduce parallel diagonalization
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More 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 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 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 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 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 informationBand Structure Calculations; Electronic and Optical Properties
; Electronic and Optical Properties Stewart Clark University of Outline Introduction to band structures Calculating band structures using Castep Calculating optical properties Examples results Some applications
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 informationElectron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46
Electron levels in a periodic potential: general remarks Daniele Toffoli January 11, 2017 1 / 46 Outline 1 Mathematical tools 2 The periodic potential 3 Bloch s theorem and Born-von Karman boundary conditions
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 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 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 informationDomain specific libraries. Material science codes on innovative HPC architectures Anton Kozhevnikov, CSCS December 5, 2016
Domain specific libraries Material science codes on innovative HPC architectures Anton Kozhevnikov, CSCS December 5, 2016 Part 1: Introduction Kohn-Shame equations 1 2 Eigen-value problem + v eff (r) j(r)
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 informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
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 informationDegenerate Perturbation Theory. 1 General framework and strategy
Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying
More informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
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 informationCRYSTAL in parallel: replicated and distributed (MPP) data. Why parallel?
CRYSTAL in parallel: replicated and distributed (MPP) data Roberto Orlando Dipartimento di Chimica Università di Torino Via Pietro Giuria 5, 10125 Torino (Italy) roberto.orlando@unito.it 1 Why parallel?
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 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 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 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 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 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 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 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 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 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 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 informationApplications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices
Applications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices Vahid Dehdari and Clayton V. Deutsch Geostatistical modeling involves many variables and many locations.
More informationFast Eigenvalue Solutions
Fast Eigenvalue Solutions Techniques! Steepest Descent/Conjugate Gradient! Davidson/Lanczos! Carr-Parrinello PDF Files will be available Where HF/DFT calculations spend time Guess ρ Form H Diagonalize
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More 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 informationThe Quantum Heisenberg Ferromagnet
The Quantum Heisenberg Ferromagnet Soon after Schrödinger discovered the wave equation of quantum mechanics, Heisenberg and Dirac developed the first successful quantum theory of ferromagnetism W. Heisenberg,
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 informationH ψ = E ψ. Introduction to Exact Diagonalization. Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden
H ψ = E ψ Introduction to Exact Diagonalization Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml laeuchli@comp-phys.org Simulations of
More informationELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers Victor Yu and the ELSI team Department of Mechanical Engineering & Materials Science Duke University Kohn-Sham Density-Functional
More informationElectron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals
Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University
More informationVASP: Dielectric response Perturbation theory, linear response, and finite electric fields
VASP: Dielectric response Perturbation theory, linear response, and finite electric fields University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria Outline
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 informationKronig-Penney model. 2m dx. Solutions can be found in region I and region II Match boundary conditions
Kronig-Penney model I II d V( x) m dx E Solutions can be found in region I and region II Match boundary conditions Linear differential equations with periodic coefficients Have exponentially decaying solutions,
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationMagnetism in transition metal oxides by post-dft methods
Magnetism in transition metal oxides by post-dft methods Cesare Franchini Faculty of Physics & Center for Computational Materials Science University of Vienna, Austria Workshop on Magnetism in Complex
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 informationDFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics,
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 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 informationProjector augmented wave Implementation
Projector augmented wave Implementation Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 = Projector augmented wave +
More 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 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 informationSelf-consistent implementation of meta-gga exchange-correlation functionals within the ONETEP linear-scaling DFT code
Self-consistent implementation of meta-gga exchange-correlation functionals within the linear-scaling DFT code, August 2016 James C. Womack,1 Narbe Mardirossian,2 Martin Head-Gordon,2,3 Chris-Kriton Skylaris1
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 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 informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
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 informationGW quasiparticle energies
Chapter 4 GW quasiparticle energies Density functional theory provides a good description of ground state properties by mapping the problem of interacting electrons onto a KS system of independent particles
More informationPerformance optimization of WEST and Qbox on Intel Knights Landing
Performance optimization of WEST and Qbox on Intel Knights Landing Huihuo Zheng 1, Christopher Knight 1, Giulia Galli 1,2, Marco Govoni 1,2, and Francois Gygi 3 1 Argonne National Laboratory 2 University
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