HECToR CSE technical meeting, Oxford Parallel Algorithms for the Materials Modelling code CRYSTAL
|
|
|
- Rudolf Glenn
- 5 years ago
- Views:
Transcription
1 HECToR CSE technical meeting, Oxford 2009 Parallel Algorithms for the Materials Modelling code CRYSTAL Dr Stanko Tomi Computational Science & Engineering Department, STFC Daresbury Laboratory, UK
2 Acknowledgements Dr Ian Bush (NAG) Dr Barry Searle (DL-STFC) Prof Nic Harrison (Imperial-DL-STFC) dcse-epsrc Initiative
3 Outline Motivation CRYSTAL code: Theoretical background CRYSTAL code: Numerical implementation CRYSTAL code: Parallelisation CRYSTAL code: Memory economisation Concussions
4 Motivation for using highly parallel CRYSTAL code a) Miniaturisation of devices and advance in current technology, requires material description at atomistic level and detailed quantum mechanical knowledge of its properties. b) New materials can not be described from empirical method because lack of parameters. Experimentally to expensive or parameter space is too big to be completely examined. c) Ab-initio model (parameter free) is the only option then d) Community that I support is interested in very big (computationally) systems: 1. bulks with dilute level of impurities (supercells) as a new materials, a building blocks, for the active region of novel PV or optoelectronics devices 2. Band alignments (offsets) at the surface interfaces: new materials, transport properties 3. Electronic structure of the QD (nano-clusters) form the first principles: new materials, PV, chemistry, biology e) The (1-3) are very challenging task and NEW algorithms are needed now!
5 Outline Motivation CRYSTAL code: Theoretical background CRYSTAL code: Numerical implementation CRYSTAL code: Parallelisation CRYSTAL code: Memory economisation Concussions
6 Hohenberg-Kohn Theory For any legal wave-function (anti-symmetric and normalised) * ˆ * E[ Ψ ] = Ψ HΨ dr = Ψ Hˆ Ψ The total energy is functional of E[ Ψ ] > E Search all (3N-coordinates) to minimize E in order to find ground state E 0!!! Theorem 1: External potential and the number of electrons in the system are uniquely determined by the density (r), so that the total energy is a unique functional of density: E=E[ (r)] 0 Theorem 2: The density that minimises the total energy is the ground state density and the minimum energy is the ground state energy, min{e[ (r)]}=e 0
7 Density Functional Theory-LDA Z α i + + Ψ ( r1,..., rn ) = EΨ( r1,..., rn ) 2 i i, j ri rj i, α ri Rα For the ground state energy and density there is an exact mapping between many body system and fictitious non-interacting system 1 ρ( r ) Z dr V ( r) Eψ ( r) 2 α XC ψ i = i i 2 r r α r Rα The fictitious system is subject to an unknown potential V XC derived from exchange-correlation functional In LDA the energy functional can be approximated as a local function of density: δ E ρ and 2 = [ ( )] XC V [ ] XC ρ r = ρ( r) ψ i ( r) δρ i The DFT Eq. describe non-interacting electrons under the influence of a mean field potential consisting of the classical Coulomb potential and local exchange-correlation potential
8 Hartree-Fock Theory Specific structure of HF WF: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) * 1 Z 1 ψ r ψ r ψ r ψ r 1 ψ r ψ r ψ r ψ r α EHF = ψ i ( r) + ψ i ( r) dr + dr dr dr dr 2 α r Rα 2 r r 2 r r * * * * i 1 i 1 j 2 j 2 i 1 j 1 j 2 i i, j i j i, j i j Coulomb energy Exchange energy After application of the variation theorem and under: ψ ψ = δ i j ij The ground state Hartree-Fock Equation: 1 2 ρ( r ) Z α + dr ψ i ( r) vx ( r, r ) ψ i ( r ) dr Eiψ i ( r) = r r α r Rα Exact non-local exchange potential: Ψ 1 det[ ψ ψ... ψ N ] N! HF 1 2 = ψ * j ( r) ψ j ( r ) ψ ( r ) dr i j r r The Hartree-Fock Eq. describe non-interacting electrons under the influence of a mean field potential consisting of the classical Coulomb potential and NON-LOCAL exchange potential
9 Self Interaction Energy In DFT-LDA each individual electron contributes to the total density An electron in an occupied orbital interacts with N electrons, when it should interact with N-1 this is the self interaction DFT-LDA: all occupied bands are pushed up in energy by this interaction consequence of an on-site/diagonal Coulomb/Exchange repulsion. DFT- LDA band gap are TOO SMALL Hartree-Fock theory: Coulomb and exchange interactions cancels exactly, i.e. J ii = K ii, no self interaction, + absent of correlations: HF band gaps TOO BIG Hybrid functional 20% HF and 80% DFT Energy gaps appear to be good!!!
10 Functionals: PBE0 & B3LYP Hybrid Functional: An exchange-correlation functional in DFT that incorporates portion of exact exchange from HF (E HF x exact) with exchange and correlation from other sources (LDA, GGA, etc.) PBE0 = Perdew-Burke-Ernzerhof E = E +0.25( E E ) PBE0 GGA HF GGA xc xc x x B3LYP = Becke exchange + (Lee-Yang-Parr + Vosko-Wilk-Nusair) correlation E = E + a ( E E ) + a ( E E ) + a ( E E ) B3LYP LDA HF GGA LDA GGA LDA GGA xc xc 1 x x 2 x x 3 c c E = 0.8 E +0.2E ( E E ) B3LYP LDA HF B88( GGA) LDA x x x x x E = 0.19 E +0.81E B3LYP VWN 3( LDA) LYP( GGA) c c c HF E x exact = LDA 3 Exc d r f [ ρ( r)] 3 E = d r f [ ρ( r), ρ( r)] GGA xc
11 Outline Motivation CRYSTAL code: Theoretical background CRYSTAL code: Numerical implementation CRYSTAL code: Parallelisation CRYSTAL code: Memory economisation Concussions
12 Implementation: Gaussian Basis Set Anti-symmetric many body wave function Φ = Â N ψ ( r, k) N i i i= 1 One-electron wave-function=lc of Bloch functions = ψ ( r, k) a ( k) χ ( r, k) i i µ, i µ i µ Bloch function=lc of Atomic Orbitals χ ( r, k) = ϕ ( r r R) i µ i µ i µ e kr R Atomic Orbitals = LC of n G Gaussian type functions (orbitals) (GTO) α n µ µ G i d jg α j i µ j= 1 ϕ ( r r R) = (, r r R) α 2 G (, µ ) = exp[ ( µ ) ]( µ ) l ( µ ) m ( µ ) n j ri r j ri r xi x yi y zi z G = G G G r x y z needs to be found coeficients exponents
13 Linear dependency problem in basis set Typically Gaussian Atomic Orbitals, (r i -r -R), are localized and non-orthogonal!!! Very diffused Gaussians (very small exponent j ) can cause linear dependency problem and trouble SCF cycles!!! Two ways to cure this problem: 1. To limit basis set to less diffused Gaussians: typically up to ~0.1 Bohr To project out several smallest eigenvalues of the overlap matrix S(k) that enters Fock equation: F(k)A(k) = S(k)E(k)A(k). In this way possible singularities in Fock matrix F(k) are avoided during SCF cycles
14 Implementation: Gaussian Basis Set Solution of secular equation in k- space gives: F( k) A( k) = S( k) E( k) A( k) a µ ( ), i k While we have basis functions (Gaussians) in r- space = F( k) F( g) i k e g Fock matrix sum of one- and two- electron contributions F ( g) = H ( g) + B ( g) i, j i, j i, j g One electron contributions: kinetic & nuclear H ( g) = T ( g) + Z ( g) = i Tˆ j + i Zˆ j i, j i, j i, j Two electron contribution: Coulomb and exchange 1 Bi, j ( g) = Ci, j ( g) + X i, j( g) = Pk, l[ i, j k, l i, k j, l ] 2 k, l NON-LOCAL term
15 Outline Motivation CRYSTAL code: Theoretical background CRYSTAL code: Numerical implementation CRYSTAL code: Parallelisation CRYSTAL code: Memory economisation Concussions
16 CRYSTAL scalability on HECToR Test 1: Bulk material (periodic solid) with dilute amount of impurities (1.62%): Ga 432 Sb 426 N 6 No of atoms: 864 Basis sets: m-pvdz-pp relativistic eff. core for Ga, Sb; and m-6-311g for N No of AO: Symmetry : 1 K-points: 1 X R + 1 x C B: 96 SCF (s) 1024 SHELLX (s) MPP_DIAG (s) B: MPP_DIAG : wall time MPP_DAIG-MPP_SIML : elapsed No of processors No of processors No of processors Problem identified: system run only when: #PBS -l mppnppn=1 Bad memory management => unnecessary expensive to run
17 CRYSTAL scalability on HECToR Test 1: New command is implemented in CRYSTAL:. MPP_BLOCK if omitted in INPUT file 16 block is calculated in. CRYSTAL code MPP_DIAG (s) No. procesors = 3584 Default MPP_BLOCK Elapse : MPP_BACK + MPP_SIML (s) No. procesors = MPP_BLOCK
18 CRYSTAL scalability on HECToR Test 2: Slab CdSe/ZnS No of atoms: 128 Basis sets: All electron basis sett for all: Cd, Se, Zn, and S No of AO: 4024 Symmetry : 2 K-points: 4 X R + 5 x C NON-DC DIAG routine SCF (s) No of processors SHELLX (s) No of processors 1 MPP_DIAG : wall time MPP_DAIG-MPP_SIML : elapsed MPP_SIML + MPP_BACK (s) CF 2.5 ideal (CF 2.5) CF 2 CF 3 CF No of processes CMPLXFAC was changed in order to find optimal load balance between complex and real k-points diagonalization MPP_DIAG (s) CF 2.5 ideal (CF 2.5) CF 2 CF 3 CF 4 No of processors
19 CRYSTAL scalability on HECToR Test 3: Cluster (QD) Cd 159 Se 166 No of atoms: 325 Basis sets: m-pvdz-pp relativistic eff. core for Cd and Se No of AO: 8747 Symmetry : 6 K-points: 1 X R SCF (s) B16 B32 Default = B96 SHELLX (s) MONMO3 (s) No of processors MPI_DIAG (s) B16 B32 Default = B96 MPP_DIAG : wall time MPP_DAIG-MPP_SIML : elapsed No of processors N o of processors No of processors
20 Outline Motivation CRYSTAL code: Theoretical background CRYSTAL code: Numerical implementation CRYSTAL code: Parallelisation CRYSTAL code: Memory economisation Concussions
21 CRYSTAL memory economization on HECToR Test 1: Bulk material (periodic solid) with dilute amount of impurities (1.62%): Ga 432 Sb 426 N 6 No of atoms: 864 No of AO: Symmetry : 1 Problem identified: system run ONLY when: #PBS -l mppnppn=1 Bad memory management => unnecessary expensive to run!!! Solution: Replicated data structure of the Fock and density matrices is replaced with direct space IRREDUCIBLE representation of the Fock and density matrices for all class of SCF calculations
22 CRYSTAL memory economization on HECToR grep CYC mppnppn_1/out.dat INFORMATION **** MAXCYCLE **** MAX NUMBER OF SCF CYCLES SET TO 2 MAX NUMBER OF SCF CYCLES 2 CONVERGENCE ON DELTAP 10**-16 CYC 0 ETOT(AU) E+05 DETOT -2.15E+05 tst 0.00E+00 PX 0.00E+00 CYC 1 ETOT(AU) E+05 DETOT 1.19E+02 tst 0.00E+00 PX 0.00E+00 CYC 2 ETOT(AU) E+05 DETOT -2.90E+00 tst 1.60E-04 PX 0.00E+00 == SCF ENDED - TOO MANY CYCLES E(AU) E+05 CYCLES 2 #PBS -l mppwidth=896 #PBS -l mppnppn=1 st53@nid00015:~/work/bulk_gasbn/896_sym> grep CYC mppnppn_2/out.dat INFORMATION **** MAXCYCLE **** MAX NUMBER OF SCF CYCLES SET TO 2 MAX NUMBER OF SCF CYCLES 2 CONVERGENCE ON DELTAP 10**-16 CYC 0 ETOT(AU) E+05 DETOT -2.15E+05 tst 0.00E+00 PX 0.00E+00 CYC 1 ETOT(AU) E+05 DETOT 1.19E+02 tst 0.00E+00 PX 0.00E+00 CYC 2 ETOT(AU) E+05 DETOT -2.90E+00 tst 1.60E-04 PX 0.00E+00 == SCF ENDED - TOO MANY CYCLES E(AU) E+05 CYCLES 2 st53@nid00015:~/work/bulk_gasbn/896_sym> grep CYC mppnppn_4/out.dat INFORMATION **** MAXCYCLE **** MAX NUMBER OF SCF CYCLES SET TO 2 MAX NUMBER OF SCF CYCLES 2 CONVERGENCE ON DELTAP 10**-16 CYC 0 ETOT(AU) E+05 DETOT -2.15E+05 tst 0.00E+00 PX 0.00E+00 CYC 1 ETOT(AU) E+05 DETOT 1.19E+02 tst 0.00E+00 PX 0.00E+00 CYC 2 ETOT(AU) E+05 DETOT -2.91E+00 tst 1.60E-04 PX 0.00E+00 == SCF ENDED - TOO MANY CYCLES E(AU) E+05 CYCLES 2 #PBS -l mppwidth=896 #PBS -l mppnppn=2 #PBS -l mppwidth=896 #PBS -l mppnppn=4
23 CRYSTAL memory economization on HAPU test_si]$ grep CYC std/out.dat INFORMATION **** MAXCYCLE **** MAX NUMBER OF SCF CYCLES SET TO 200 MAX NUMBER OF SCF CYCLES 200 CONVERGENCE ON DELTAP 10**-17 CYC 0 ETOT(AU) E+03 DETOT -4.57E+03 tst 0.00E+00 PX 1.00E+00 CYC 1 ETOT(AU) E+03 DETOT -2.96E+00 tst 0.00E+00 PX 1.00E+00 == SCF ENDED - CONVERGENCE ON ENERGY E(AU) E+03 CYCLES 9 [st53@hapu33 test_si]$ grep CYC sym/out.dat INFORMATION **** MAXCYCLE **** MAX NUMBER OF SCF CYCLES SET TO 200 MAX NUMBER OF SCF CYCLES 200 CONVERGENCE ON DELTAP 10**-17 CYC 0 ETOT(AU) E+03 DETOT -4.57E+03 tst 0.00E+00 PX 0.00E+00 CYC 1 ETOT(AU) E+03 DETOT -2.96E+00 tst 0.00E+00 PX 0.00E+00 == SCF ENDED - CONVERGENCE ON ENERGY E(AU) E+03 CYCLES 9 [ [st53@hapu33 test_si]$ grep CYC sym_nolm/out.dat INFORMATION **** MAXCYCLE **** MAX NUMBER OF SCF CYCLES SET TO 200 MAX NUMBER OF SCF CYCLES 200 CONVERGENCE ON DELTAP 10**-17 CYC 0 ETOT(AU) E+03 DETOT -4.57E+03 tst 0.00E+00 PX 0.00E+00 CYC 1 ETOT(AU) E+03 DETOT -2.96E+00 tst 0.00E+00 PX 0.00E+00 == SCF ENDED - CONVERGENCE ON ENERGY E(AU) E+03 CYCLES 9 Old code: replicated data in memory MPP New code: irreducible representation MPP LOWMEM New code: irreducible representation MPP NOLOWMEM
24 CRYSTAL memory economization on HECToR HECToR phase 1 HECToR phase 2 1 core/node = 6 GB/core 2 core/node = 3 GB/core 1 core/node = 8 GB/core 2 core/node = 4 GB/core 4 core/node = 2 GB/core Memory economisation achieved: from more then 3 GB/core required to less the 2 GB/core, i.e., the whole node can be used efficiently.
25 Outline Motivation CRYSTAL code: Theoretical background CRYSTAL code: Numerical implementation CRYSTAL code: Parallelisation CRYSTAL code: Memory economisation Concussions
26 Conclusions 1) Memory bottleneck in the CRYSTAL code caused by replicated data structure of the Fock and density matrix is removed by replacing those data structures with its irreducible representations [saving: problem size/(2 x symmetry of the system)]. 2) For big-ish systems the optimal load balance between complex and real k-points is achieved with CMPLXFAC=3, that provide for best scalability of the CRYSTAL code up to 4864 cores on HECToR. 3) Implementation of the new command MPP_BLOCK provides for better control of the block size of distributed data. Reducing block size from default value of 96 to 64 or 32 it can be achieved more than 10% speed-up in diagonalisation part and about 20% speed-up in back and similarity transform for the system size rank=19837 run on 3584 cores. 4) For systems when integral calculations time prevails (MONMO3, SHELLX) over diagonalisation time (MPP_DIAG) better synchronisation of integral calculations is needed: possible global or shared counter implementations investigation. 5) For systems where diagonalisation (MPP_DIAG) dominates over integral calculations (MONMO3, SHELLX) careful examination of blocking is needed
Advanced 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
Quantum 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
DFT 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
An 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
Density 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
Computational 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
Institut 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
Teoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
Ab-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
Electron 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
Exchange 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,
The 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
DENSITY 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
Computational 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
Introduction to Computational Quantum Chemistry: Theory
Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction
Self-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,
Electronic 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.
Introduction 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
Oslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
Multi-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
CRYSTAL 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?
Density 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
Same idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)
XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low
Density Func,onal Theory (Chapter 6, Jensen)
Chem 580: DFT Density Func,onal Theory (Chapter 6, Jensen) Hohenberg- Kohn Theorem (Phys. Rev., 136,B864 (1964)): For molecules with a non degenerate ground state, the ground state molecular energy and
Introduction to Computational Chemistry Computational (chemistry education) and/or. (Computational chemistry) education
Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools to help increase student understanding of material
Orbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
Finite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals. Travis Sjostrom. IPAM 2012 Workshop IV
1 of 45 Finite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals Travis Sjostrom Quantum Theory Project Depts. of Physics and Chemistry IPAM 2012 Workshop IV 2012 2 of
Density 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
OVERVIEW 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
Density 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
Modified 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
Introduction to Computational Chemistry: Theory
Introduction to Computational Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC andrew.gilbert@anu.edu.au 3023 Course Lectures Introduction Hartree Fock Theory Basis Sets Lecture 1 1 Introduction
Computational Modeling Software and their applications
Computational Modeling Software and their applications June 21, 2011 Damilola Daramola Center for Electrochemical Engineering Research ABC s of electrochemistry Introduction Computational Modeling the
Introduction 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
DFT: 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
CLIMBING 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
1 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
Basics of DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Basics of DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Basics of DFT Lausanne12 1
How to generate a pseudopotential with non-linear core corrections
How to generate a pseudopotential with non-linear core corrections 14 12 AE core charge AE valence charge PS core charge PS valence charge 10 8 6 4 2 Objectives 0 0 0.5 1 1.5 2 2.5 3 Check whether the
ABC of ground-state DFT
ABC of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) ABC of ground-state
Module 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
2.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.
Introduction 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
Fast 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
Bert de Groot, Udo Schmitt, Helmut Grubmüller
Computergestützte Biophysik I Bert de Groot, Udo Schmitt, Helmut Grubmüller Max Planck-Institut für biophysikalische Chemie Theoretische und Computergestützte Biophysik Am Fassberg 11 37077 Göttingen Tel.:
DFT: Exchange-Correlation
DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,
Study of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal
Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
Exchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
Density functional theory in the solid state
Density functional theory in the solid state Ari P Seitsonen IMPMC, CNRS & Universités 6 et 7 Paris, IPGP Department of Applied Physics, Helsinki University of Technology Physikalisch-Chemisches Institut
Density 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
2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
Supporting Information for. Predicting the Stability of Fullerene Allotropes Throughout the Periodic Table MA 02139
Supporting Information for Predicting the Stability of Fullerene Allotropes Throughout the Periodic Table Qing Zhao 1, 2, Stanley S. H. Ng 1, and Heather J. Kulik 1, * 1 Department of Chemical Engineering,
Electronic Supplementary Information
Electronic Supplementary Material (ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2014 Electronic Supplementary Information Configurational and energetical study of the (100) and
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen
AB INITIO MODELING OF ALKALI METAL CHALCOGENIDES USING SOGGA THEORY
Int. J. Chem. Sci.: 13(4), 215, 163-1638 ISSN 972-768X www.sadgurupublications.com AB INITIO MODELING OF ALALI METAL CHALCOGENIDES USING SOGGA THEORY HITESH CHANDRA SWARNAR and GUNJAN ARORA a,* Department
Spring 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
ASSESSMENT OF DFT METHODS FOR SOLIDS
MSSC2009 - Ab Initio Modeling in Solid State Chemistry ASSESSMENT OF DFT METHODS FOR SOLIDS Raffaella Demichelis Università di Torino Dipartimento di Chimica IFM 1 MSSC2009 - September, 10 th 2009 Table
Introduction 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
Short Course on Density Functional Theory and Applications VII. Hybrid, Range-Separated, and One-shot Functionals
Short Course on Density Functional Theory and Applications VII. Hybrid, Range-Separated, and One-shot Functionals Samuel B. Trickey Sept. 2008 Quantum Theory Project Dept. of Physics and Dept. of Chemistry
Session 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
This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
One-Electron Hamiltonians
One-Electron Hamiltonians Hartree-Fock and Density Func7onal Theory Christopher J. Cramer @ChemProfCramer 2017 MSSC, July 10, 2017 REVIEW A One-Slide Summary of Quantum Mechanics Fundamental Postulate:
Quantum Chemical and Dynamical Tools for Solving Photochemical Problems
2.165430 3.413060 3.889592 9 H 3.413060 2.165430 1.099610 2.165430 3.413060 10 H 3.889592 3.413060 2.165430 1.099610 2.165430 11 H 3.413060 3.889592 3.413060 2.165430 1.099610 12 H 2.165430 3.413060 3.889592
QMC dissociation energy of the water dimer: Time step errors and backflow calculations
QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.
Electronic 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,
Density Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
Calculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
NWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory
NWChem: Hartree-Fock, Density Functional Theory, Time-Depent Density Functional Theory Hartree-Fock! Functionality! Input! Wavefunctions! Initial MO vectors! Direct and semidirect algorithms! Convergence,
Electrochemistry 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
Answers 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
Lecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
Parallel 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
XYZ of ground-state DFT
XYZ of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) XYZ of ground-state
7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional
Ab initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
Electronic structure calculations: fundamentals George C. Schatz Northwestern University
Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic Structure (often called Quantum Chemistry) calculations use quantum mechanics to determine the wavefunctions
Introduction 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
Lecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
Post 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,
Basic introduction of NWChem software
Basic introduction of NWChem software Background NWChem is part of the Molecular Science Software Suite Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry
Density 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
Set the initial conditions r i. Update neighborlist. Get new forces F i
v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )
Gaussian Basis Sets for Solid-State Calculations
Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction
Structure 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
Computer Laboratory DFT and Frequencies
Computer Laboratory DFT and Frequencies 1. DFT KEYWORDS FOR DFT METHODS Names for the various pure DFT models are given by combining the names for the exchange and correlation functionals. In some cases,
Density 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
Reactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci)
Reactivity and Organocatalysis (Adalgisa Sinicropi and Massimo Olivucci) The Aldol Reaction - O R 1 O R 1 + - O O OH * * H R 2 R 1 R 2 The (1957) Zimmerman-Traxler Transition State Counterion (e.g. Li
Algorithms 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
Al 2 O 3. Al Barrier Layer. Pores
Overview of DFT Calculations for Nanoscience and Nanotechnology Oğuz Gülseren Al 2 O 3 Al Barrier Layer Pores Jules Verne A Trip to the Moon, 1865 Nanostructures: a. Contain a countable number of atoms
Translation 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
Intermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Intermediate DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Intermediate DFT Lausanne12
DFT 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,
MBPT and TDDFT Theory and Tools for Electronic-Optical Properties Calculations in Material Science
MBPT and TDDFT Theory and Tools for Electronic-Optical Properties Calculations in Material Science Dott.ssa Letizia Chiodo Nano-bio Spectroscopy Group & ETSF - European Theoretical Spectroscopy Facility,
ABC of ground-state DFT
ABC of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA July 31, 2014 Kieron (UC Irvine) ABC of ground-state DFT HoW
Introduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_
Time-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