CP2K: the gaussian plane wave (GPW) method
|
|
- Andrea Willis
- 5 years ago
- Views:
Transcription
1 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
2 CPMD and CP2K CPMD CP2K Licensed by IBM and MPI Stuttgart Fortran 77 with extension Plane wave basis set heavy use of FFT possibility of k-point sampling GNU license high level fortran 90 Mixed gaussians and plane waves order N construction of the KS Hamiltonian
3 The problem We want to solve numerically the self-consistent Kohn-Sham equations: v KS ( r ) = v( r ) + d 3 r n 0( r ) r r + v xc[n 0 ]( r ) with N n 0 ( r ) = ϕ 0 i ( r ) 2 i= ϕ 0 i ( r ) + v KS ( r )ϕ 0 i ( r ) = ɛ i ϕ 0 i ( r ) (1) We need 1. a representation of φ i ( r ) and n( r ), 2. to be able to compute the total Kohn-Sham energy and its derivatives 3. a method to solve the non-linear the Kohn-Sham equation.
4 Plane wave vs gaussian basis sets: plane waves pros and cons Advantages independent of the nuclei position (good for forces) no BSSE one parameter controls the basis set size orthogonal numerical efficiency through use of FFT Disadvantages large number of basis set elements needed Necessary use of pseudo-potentials loss of chemical insight Border line fill the whole simulation box naturally periodic
5 Plane wave vs gaussian basis sets: gaussians pros and cons Advantages Good already for small basis set sizes correspond to chemical insight Computationally efficient (multi-centre integrals) Possibility to perform all-electrons calculations Disadvantages non-orthogonal atomic position dependent (Pulay forces) Basis set superposition error (BSSE) Systematic improvment not straightforward linear dependencies, over-completeness wrong asymptotic behaviour Border line no implicit periodicity can be tuned for each application
6 Plane wave basis set: Definition Basis set {ϕ α } of plane waves: (Ω is the volume of the box) ϕ α ( r ) = 1 Ω e i G α r The wavefunctions are expanded as φ i ( r ) = 1 c i ( G )e i G r Ω G (2) (3) For an orthorhombic box with lengths L x, L y and L z, the wavevectors G are G = i 2π L x x + j 2π L y y + k 2π L z z ; with i, j, k Z real wavefunctions only half space is needed
7 Representation of the electronic density n( r ) n( r ) = φ i ( r ) 2 = 1 Ω i n( r ) = 1 ñ( G Ω )e i G r G <2G max i G G c i ( G )c i ( G )e i( G G ) r (4) E n cut = 4E φ cut
8 Auxiliary real space grid Sampling interval = L N Nyquist critical frequency f c = 2π 2 G max related to E cut X i = (i 1) on the real space grid We can go from the real-space grid to the reciprocal-space grid and back using Fourier transform techniques If N is a product of small prime numbers one can use Fast Fourier Transform (FFT) techniques which are very efficient (N log N instead of N 2 )
9 Density from wavefunction coefficients Real grid spacing adapted to Ecut n ñ( G ) = Ω n( R N k )e i G Rk = Ω φ i ( R N k ) 2 e i G R k (5) R k i To calculate ñ( G ) we heavily use the FFT: R k c i ( G) (S)INV FFT φ i ( R k ) n( R k ) = i φ i ( R k ) 2 FFT ñ( G )
10 Density derivatives An important quantity will be n( R k ) c i ( G) This is n( R k ) c i ( G) = Rk n( R k ) φ i ( R k ) φ i ( R k ) ci ( G) = 1 Ω φ i ( R k )e i G R k
11 Gaussian basis set Primitive cartesian gaussian functions: g( r, η, l, A) = N c (x A x ) lx (y A y ) ly (z A z ) lz e η( r A) 2 (6) Problems: r: Electron coordinate A: Atomic coordinate l: Angular momentum No cusp at the origin wrong asymptotic behaviour Contracted cartesian gaussian functions: ϕ µ ( r) = i d µi g µi ( r) (7)
12 CP2K: the GPW method Joost VandeVondele et al. Quickstep: Fast and accurate density functional calculations using a mixed gaussian and plane waves approach. Computer Physics Communications, 167: , Two representations of the density: local basis set: n( r) = µν P µν ϕ µ ( r)ϕ ν ( r), (8) P µν : density matrix and plane waves: n( r) = 1 ñ( G) Ω exp(ig r) (9) G <G max
13 P µν is the central object to be determined E KS P µν = K µν is the Kohn-Sham matrix We need to solve the generalized eigenvalue problem K µν ci ν = ɛ i S µν ci ν ν where S µν = ϕ µ ϕ ν is the overlap matrix, with the constraint n occ and P µν = c µ i ci ν i=1 c µ j S µνci ν µ,ν ν = δ ij
14 Going from local basis set to plane waves: auxiliary grid Knowledge of ñ( G) is determined by the knowledge of n( R k ) on a regular mesh R k : n( R k ) = µν P µν ϕ µ ( R k )ϕ ν ( R k ), (10) Here too, heavy use of FFT
15 Order N construction Product of two gaussians is a gaussian centered at: P = η a A + η b B η a + η b (11) with prefactor exp η a η b η a + η b A B 2 mesh points far from P are put to zero (µ, ν) pairs far apart are not considered order N achieved using linked-cell method
16 density derivatives We will need derivatives fo the density: n( R k ) P µ,ν n( R k ) P µ,ν = ϕ µ ( R k )ϕ ν ( R k ) (12) at the points R µ,ν where ϕ µ ( R k )ϕ ν ( R k ) is calculated
17 Multi-grids For some products, a lower cutoff the G m ax can be used. In principle every pair could have a different cutoff; in practice 3 to 5. n( r) = 1 ñ 1 ( G) Ω exp(i G r) + ñ 2 ( G) exp(ig r) G <G 1 G <G 2 Products with a lower cutoff are colocated on a coarser grid adapted to the lower cutoff. The coarser grid is map on to the finer grid by Fourier transforms (FFT s) (13)
18 derivatives: chain rule n( R k ) P µ,ν = n( R k ) n( R k) n( R R k) P µ,ν k = 1 exp ig N ( R k R k) ϕ µ ( R k)ϕ ν ( R k) R k G <G 2 = 1 exp ig Ω R k Ω exp ig N R k ϕ µ ( R k)ϕ ν ( R k) G <G 2 R k Use of FFT
19 Calculation of the Kohn-Sham energy and its derivatives E[n] = T s [n] + Kinetic energy d 3 r n( r )v( r ) + J[n] + E xc [n] I J Z I Z J R I R J ϕ µ ϕ ν is analytic
20 Goedecker-Tetter-Hutter (GTH) pseudopotentials Local part v PP loc LR SR (r) =vloc (r) + vloc (r) = Z r erf(αr) + 4 i=1 C PP i ( 2αr) 2i 2 exp [ (αr) 2] Non-local part with the projectors v SR nl = lm ij r p lm i r p i lm = Ni l Y lm (ˆr)r l+2i 2 exp hij p l j lm r [ 1 2 ( ) ] r 2 r l
21 all short range terms are analytic in the local basis set representation
22 Electrostatic contribution E ES = 1 2 Ω d 3 rd 3 r n( r )n( r ) r r + I Ω d 3 r vi c ( r )n( r ) I J Z I Z J R I R J Using the Fourier transformed representation and with PBC: just like previously 1 2 Ω d 3 rd 3 r n tot( r )n tot ( r ) r r = 2πΩ G 0 ñ tot ( G ) 2 G 2 ñ tot ( G ) = ñ( G ) + I ñ c I ( G ) exp( i G R I )
23 Exchange-correlation energy E xc [n] = is, just like in CPMD, approximated by E xc [n] = Ω N Differences with cpmd Ω d 3 r ɛ xc (n( r ))n( r ) R k ɛ xc (n( R k ))n( R k ) Derivatives of the energy by Fourier transformed or smoothed finite difference Density is first smoothed using nearest-neighbours (reducing the ripples )
24 Smoothing operator q 3 (S q f ) i,j,k = q 3 + 6q q q l m n f i+l,j+m,k+n q = 10 to 50 l= 1 m= 1 n= 1 Chain rule used again for derivatives smoothing of the xc potential
25 Derivatives of the electrostatic and xc contributions (E ES + E xc ) n( R k ) = Ω N (V H( R k ) + V xc ( R k )) V H ( R k ) obtained by FFT from ṼH( G) From chain rule we obtain (E ES + E xc ) P µν = Ω N R k ( V H ( R k ) + V xc ( R ) k ) ϕ µ ( R k )ϕ ν ( R k )
26 Optimization Two ways of solving the Kohn-Sham equation: 1. Reach for self-consistency (traditional scheme) Obtain lowest occupied orbitals from diagonalization of the KS matrix for a test density Density from occupied orbitals Update test density until self-consistency 2. Minimization of the energy functional Use a non-linear minimization scheme to minimize the KS energy
27 CP2K (1): traditional scheme Diagonalization Given the Kohn-Sham Hamiltonian K µ,ν constructed previously, we need to solve Kc = Scɛ (14) where, c are the wavefunction coefficient in the local basis set and S the overlap matrix. ɛ is the diagonal matrix of eigenvalues. S in the local basis set can be calculated analytically.
28 Direct diagonalization Since the number of local basis set elements is not prohibitively large (up to a few thousands) the diagonalization can be done directly using linear algebra packages (lapack or scalapack for parallelized versions) The generalized eigenvalue problem is first transformed to a normal eigenvalue problem introducing U T U = S: K c = c ɛ (15) with ( ) K = U T 1 K (U) 1 c = Uc (16) Cholesky decomposition or symmetric decomposition U = S 1 2
29 Update of the density matrix P A new density matrix is obtained: (P out ) µ,ν = N occ n=1 c n,µc n,ν The resuting density is used to update the test density matrix. Mixing scheme P i+1 in = αp i out + (1 α)p i in (17)
30 Orbital transformation: direct minimization in CP2K Use of an orbital transformation: c(x) = c 0 cos U + xu 1 sin U (18) with c O Sc 0 = I and the linear constraint c O SX = 0 and U = x T Sx
31 Forces E[n] = T s [n] + d 3 r n( r )v( r ) + J[n] + E xc [n] I J Z I Z J R I R J We need to calculate I E[n] at the optimum density n. Three types of terms: 1. v depends explicitly on R I : term d 3 r n( r ) I v( r ) 2. the basis set elements ϕ µ depend on R I : term 2 µ,ν Pµ,ν I ϕ µ E n ϕ ν 3. P µ,ν itself depends on R I : term µ,ν K µ,ν I P µ,ν If the overlap matrix is independent of R I (e.g. plane waves), the last term vanishes. Otherwise, extra force: Pulay forces
Large 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 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 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 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 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 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 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 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 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 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 informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials
More informationBasis sets for SIESTA. Emilio Artacho. Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge
Basis sets for SIESTA Emilio Artacho Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge Solving: Basis set Expand in terms of a finite set of basis functions
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 informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
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 informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
More informationElectronic Structure Calculations and Density Functional Theory
Electronic Structure Calculations and Density Functional Theory Rodolphe Vuilleumier Pôle de chimie théorique Département de chimie de l ENS CNRS Ecole normale supérieure UPMC Formation ModPhyChem Lyon,
More informationDFT with Hybrid Functionals
DFT with Hybrid Functionals Sanliang Ling University College London 4th CP2K Tutorial, 31st August 4th September 2015, Zurich What are hybrid functionals? Hybrid functionals: mixing non-local Hartree-Fock
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 informationPost Hartree-Fock: MP2 and RPA in CP2K
Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,
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 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 informationThe Plane-wave Pseudopotential Method
The Plane-wave Pseudopotential Method k(r) = X G c k,g e i(g+k) r Chris J Pickard Electrons in a Solid Nearly Free Electrons Nearly Free Electrons Nearly Free Electrons Electronic Structures Methods Empirical
More 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 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 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 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 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 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 informationProjector-Augmented Wave Method:
Projector-Augmented Wave Method: An introduction Peter E. Blöchl Clausthal University of Technology Germany http://www.pt.tu-clausthal.de/atp/ 23. Juli 2003 Why PAW all-electron wave functions (EFG s,
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 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 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 informationI. CSFs Are Used to Express the Full N-Electron Wavefunction
Chapter 11 One Must be Able to Evaluate the Matrix Elements Among Properly Symmetry Adapted N- Electron Configuration Functions for Any Operator, the Electronic Hamiltonian in Particular. The Slater-Condon
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 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 informationCP2K. New Frontiers. ab initio Molecular Dynamics
CP2K New Frontiers in ab initio Molecular Dynamics Jürg Hutter, Joost VandeVondele, Valery Weber Physical-Chemistry Institute, University of Zurich Ab Initio Molecular Dynamics Molecular Dynamics Sampling
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationDFT / 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 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 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 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 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 informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More 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 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: 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 informationLecture 6. Tight-binding model
Lecture 6 Tight-binding model In Lecture 3 we discussed the Krönig-Penny model and we have seen that, depending on the strength of the potential barrier inside the unit cell, the electrons can behave like
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
More 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 informationKey concepts in Density Functional Theory (II)
Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université
More informationNorm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria
Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized
More 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 informationarxiv: v1 [physics.comp-ph] 22 Oct 2015
arxiv:1510.06489v1 [physics.comp-ph] 22 Oct 2015 Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework II: Force, vibration, and molecular dynamics calculations
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 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 information6: Plane waves, unit cells, k- points and all that
The Nuts and Bolts of First-Principles Simulation 6: Plane waves, unit cells, k- points and all that Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview
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 informationWavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts.
Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Silvia Casassa Università degli Studi di Torino July 12, 2017 Minnesota Workshop on ab initio MSC Symmetry
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 informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
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 informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
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 informationEfficient use of CP2K
Efficient use of CP2K Fawzi Mohamed 1 1 Department of Chemisty Humboldt University of Berlin June 10, 2009 cp2k Fortran 95 opensource (GPL) general framework for different methods DFT no official release
More informationFirst-Principles Wannier Functions of Silicon and Gallium. Arsenide arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 22 Nov 1996.
First-Principles Wannier Functions of Silicon and Gallium Arsenide arxiv:cond-mat/9611176v1 [cond-mat.mtrl-sci] 22 Nov 1996 Pablo Fernández 1, Andrea Dal Corso 1, Francesco Mauri 2, and Alfonso Baldereschi
More informationPseudopotentials and Basis Sets. How to generate and test them
Pseudopotentials and Basis Sets How to generate and test them Pseudopotential idea Atomic Si Core electrons highly localized very depth energy are chemically inert 1s 2 2s 2 2p 6 3s 2 3p 2 Valence wave
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More 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 informationZurich Open Repository and Archive. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases
University of Zurich Zurich Open Repository and Archive Winterthurerstr. 190 CH-8057 Zurich http://www.zora.uzh.ch Year: 2007 Gaussian basis sets for accurate calculations on molecular systems in gas and
More informationarxiv:cond-mat/ v1 17 May 1995
Projection of plane-wave calculations into atomic orbitals Daniel Sanchez-Portal, Emilio Artacho, and Jose M. Soler Instituto de Ciencia de Materiales Nicolás Cabrera and Departamento de Física de la Materia
More informationPractical calculations using first-principles QM Convergence, convergence, convergence
Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................
More informationAccurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids
Accurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids Published as Computer Physics Communications 140, 315 (2001) Chris-Kriton Skylaris,
More informationAb initio molecular dynamics: Basic Theory and Advanced Met
Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site
More informationDGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations The recently developed discontinuous Galerkin density functional theory (DGDFT)[21] aims at reducing the number
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 informationBasics of Density Functional Theory (DFT)
Basics of Density Functional Theory (DFT) Ari Paavo SEITSONEN Ari.P.Seitsonen@iki.fi Département de Chimie École Normale Supérieure, Paris École de Sidi-BelAbbès de Nanomateriaux // Octobre 8-12, 2016
More 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 informationLecture 3: methods and terminology, part I
Eritis sicut Deus, scientes bonum et malum. Lecture 3: methods and terminology, part I Molecular dynamics, internal coordinates, Z-matrices, force fields, periodic boundary conditions, Hartree-Fock theory,
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
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 informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
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 informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
More 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 informationDensity functionals from deep learning
Density functionals from deep learning Jeffrey M. McMahon Department of Physics & Astronomy March 15, 2016 Jeffrey M. McMahon (WSU) March 15, 2016 1 / 18 Kohn Sham Density-functional Theory (KS-DFT) The
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 informationJack Simons Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
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 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 informationQUALIFYING EXAMINATION, Part 2. Solutions. Problem 1: Quantum Mechanics I
QUALIFYING EXAMINATION, Part Solutions Problem 1: Quantum Mechanics I (a) We may decompose the Hamiltonian into two parts: H = H 1 + H, ( ) where H j = 1 m p j + 1 mω x j = ω a j a j + 1/ with eigenenergies
More informationFast and accurate Coulomb calculation with Gaussian functions
Fast and accurate Coulomb calculation with Gaussian functions László Füsti-Molnár and Jing Kong Q-CHEM Inc., Pittsburgh, Pennysylvania 15213 THE JOURNAL OF CHEMICAL PHYSICS 122, 074108 2005 Received 8
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
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 information1 Construction of norm-conserving semi-local pseudopotentials for Si
1 Construction of norm-conserving semi-local pseudopotentials for Si As discussed in class, it is desirable to replace the effective interaction of the valence electrons with the ionic core, i.e. nucleus
More informationData structure of OpenMX and ADPACK
Data structure of OpenMX and ADPACK Structure of ADPACK Solving the 1D Dirac equation Pseudopotentials Localized basis functions Structure of OpenMX Input_std() Total energy Data structure for parallelization
More informationSolving Many-Body Schrödinger Equation Using Density Functional Theory and Finite Elements
Solving Many-Body Schrödinger Equation Using Density Functional Theory and Finite Elements Institute of Physics, Academy of Sciences of the Czech Republic June 21, 2008 Introduction Contens Density Functional
More informationKey concepts in Density Functional Theory (II) Silvana Botti
Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:
More information