ASSESSMENT OF DFT METHODS FOR SOLIDS

Size: px
Start display at page:

Download "ASSESSMENT OF DFT METHODS FOR SOLIDS"

Transcription

1 MSSC Ab Initio Modeling in Solid State Chemistry ASSESSMENT OF DFT METHODS FOR SOLIDS Raffaella Demichelis Università di Torino Dipartimento di Chimica IFM 1 MSSC September, 10 th 2009

2 Table of Contents DFT: a quick overview PBE s type functionals Hybrid GGA functionals DFT & CRYSTAL CRYSTAL input for DFT calculation DFT functionals available in CRYSTAL Integration grid and numerical accuracy control Performance of DFT functionals for describing periodic systems quasi-1d systems: SWCNTs and Chrysotile Aluminosilicates Relative stability of polymorphs 2

3 DFT: a quick overview E=E[ρ(r)]=E nuclear +E kinetic +E nucl-electron +E coulomb +E XC DFT energy Hohenberg, Kohn; Phys. Rev. 1964, 136, 864 Kohn, Sham, Phys. Rev. 1965, 140, 1133 The Jacob s Ladder (Perdew, Schmidt; AIP Conf. Proc. 2001, 577, 1-20) DFT HEAVEN Exact exchange, exact partial correlation Exact exchange, compatible correlation meta-gga GGA LDA HARTREE WORLD 3

4 DFT: a quick overview PBE s type exchange functionals E X = % $! " # X," [! " ] F(s " [! ",&! " ])dr " Reduced spin density Spin density Exchange energy density Different forms of F(s σ ) were proposed: PBE (1996) WC* (2006) revpbe (1998) PBEsol* (2008) RPBE (1999) SOGGA* (2008) [1/2 PBE+ 1/2 RPBE] mpbe (2002) *CRYSTAL09 (B. Civalleri ) 4

5 DFT: a quick overview HF/DFT hybrid functionals E XC =a*[e X (HF)-E X (DFT)] + E XC (DFT) (e.g. PBE0) E X =(1-a)*[E X (LDA)+b*E X (XXX)]+a*E X (HF) E C =(1-c)*E C (VWN)+c*E C (YYY) (XXX, YYY: non-local exchange and correlation) a=fraction of exact exchange (e.g. B3LYP: XXX=B, YYY=LYP B3PW: XXX=B, YYY=PW91) 5

6 DFT: a quick overview HF/DFT hybrid functionals Use of hybrid functionals in solid state simulation *a matter of discussion in the scientific community claimed lacks in accuracy use of empirical parameters bad description of narrow-gap and conducting systems *largely diffuse (e.g. B3LYP, B3PW, PBE0) and often proved to provide the best results for structural, electronic, thermodynamical and dynamical properties of nonconducting systems (in particular H-containing systems) F. Corà, M. Alfredsson, G. Mallia, D.S. Middlemiss, W.C. Mackrodt, R. Dovesi, R. Orlando, Volume 113, Springer-Verlag, Heidelberg (2004). 6

7 DFT: a quick overview HF/DFT hybrid functionals The example of Brucite, Mg(OH) 2 Hybrids Differences between calculated and experimental results: a,c: lattice parameters, Å ΔV%: volumes, relative difference (calc-exp)/exp*100 OH: bond length of OH group, Å v: anharmonic stretching of the OH oscillator, cm -1 7

8 DFT: a quick overview HF/DFT hybrid functionals J/K SVWN PBE B3LYP PBE0 HF Exp KMnF 3 exchange coupling constant (J in Kelvin) Calculated by A. Meyer, CRYSTAL06, 2009 pyrope grossular andradite spessartine uvarovite SVWN PBE B3LYP PBE0 HF Exp Closed-shell Open-shell Calculated vs experimental dielectric constants of garnets. No experimental data for almandine A. Meyer, M. Ferrero, L. Valenzano, C. M. Zicovich-Wilson, R. Orlando, R. Dovesi In preparation 8

9 CRYSTAL input for DFT calculation title 1 - GEOMETRY block 2 - BASIS SET block 3 - HAMILTONIAN and SCF block Structure of CRYSTAL input In block 3: DFT options END DFT input block Options: 1 - Choice of the exchange-correlation functionals 2 - Integration grid and numerical accuracy control 3 - Atomic parameters (and add SPIN for open-shell systems) For points 2-3 default values are supplied Point 1 must be set 9

10 DFT functionals available in CRYSTAL EXCHANGE CORRELATION LDA LDA (=S) PWLSD VWN PZ VBH VBH GGA BECKE LYP PWGGA (=PW91) PWGGA (=PW91) PBE PBE PBESOL* PBESOL* WCGGA* (=WC) P86 SOGGA* WL *CRYSTAL09 (References in the CRYSTAL users manual) 10

11 CRYSTAL input for DFT calculation 1 - Choice of the functional DFT EXCHANGE.. CORRELAT.. If not set, E X (HF) used If not set, no E C HYBRID NONLOCAL Percentage of E X (HF), if not set: pure DFT Non-local weighting parameters for E C 11

12 CRYSTAL input for DFT calculation 1 - Choice of the functional GLOBAL KEYWORDS EXCHANGE CORRELAT HYBRID NONLOCAL SOGGAXC* SOGGA PBE - - PBE1PBE* PBE PBE 25 - (PBE0) B3PW BECKE PWGGA B3LYP BECKE LYP *CRYSTAL09 (References in the CRYSTAL users manual) 12

13 Integration grid and numerical accuracy control E XC and its gradient are evaluated by numerical integration over the cell volume: an integration grid is required - Grid based on atomic partition method (Becke), V cell partitioned into atomic volumes: V A centered on atom A - Radial and angular components (Gauss-Legendre radial quadrature, Lebedev 2-dimensional angular point distribution; points associated to a weight, w) Uniform (or unpruned) grid 13

14 Integration grid and numerical accuracy control For reducing the grid size and maintaining its effectiveness: pruned grid (not radially uniform): the number of angular points is maximum in the valence region and gradually decreases in both directions 14

15 Integration grid and numerical accuracy control O O The example of MgO (R. Orlando ) Mg Grid of Mg Grid of O Radial points: Gauss formula (n r =number of radial points) O 001 plane of a unit cell of MgO O Angular points: Lebedev distribution (L=Lebedev accuracy parameter)

16 Integration grid and numerical accuracy control unpruned grid pruned grid n r =55 L=13 n r =75 L=13 n r =75 L= (31878) (42956) (96038) (11552) (20542) (41530) (R. Orlando )

17 Integration grid and numerical accuracy control: use of symmetry 55, 13 75, 13 75, 16 size of grid t elapsed (s) (R. Orlando )

18 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control RADIAL ANGULAR Number of intervals, integration interval limits, number of radial points in the i th interval, n r Number of intervals, integration interval limits, accuracy level (L) of each integral 18

19 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control Default grid: (55,434)p angular points in the valence region RADIAL ANGULAR ( ) Angular points generated n r L 19

20 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control GLOBAL KEYWORDS: LGRID --> n r =75, L=13; (75,434)p XLGRID --> n r =75, L=16; (75,974)p XXLGRID --> n r =99, L=18; (99,1454)p (see CRYSTAL users manual for details) Numerically integrated density (NID) Diaspore, 4AlOOH in the unit cell, 120 electrons Grid NID Δ DEFAULT E-03 XLGRID E-05 XXLGRID E-05 (n r =99,L=22) E-06 20

21 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control Calcite: effect of the grid on geometry and vibrational modes default LGRID XLGRID 75,16 unpr 75,18 unpr 0 Reference (IR+Raman+silent) Average of Abs diff. Maximum difference Minimum difference Grid points µ-electron (TOT: 100 e) µ-hartree Å Å M. Prencipe, F. Pascale, C. M. Zicovich-Wilson, V. R. Saunders, R. Orlando, R. Dovesi; Phys. Chem. Miner., 2004, 31,

22 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control OH vibrational data of brucite: Effect of pruned and unpruned grid size n r,l N. points O-H ν(oh)h Pruned Np O-H ν(oh)h 55, , (LGRID) , , , (XLGRID) , , (uniform) S. Tosoni, F. Pascale, P. Ugliengo, R. Orlando, V. R. Saunders, R. Dovesi; Mol. Phys., 2005, 103,

23 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control Details on the grid size and its performance: M. Prencipe, F. Pascale, C. M. Zicovich-Wilson, V. R. Saunders, R. Orlando, R. Dovesi; Phys. Chem. Miner., 2004, 31, F. Pascale, C. M. Zicovich-Wilson, R. Orlando, C. Roetti, P. Ugliengo, R. Dovesi; J. Phys. Chem. B, 2005, 109, S. Tosoni, F. Pascale, P. Ugliengo, R. Orlando, V. R. Saunders, R. Dovesi; Mol. Phys., 2005, 103,

24 CRYSTAL input for DFT calculation 2 - Integration grid and numerical accuracy control - unpruned grid can be set - thresholds on grid weigths and electronic density can be increased or decreased with respect to default - different methods for the point weights - keywords for memory menagement 3 - Atomic parameters - set atomic radius and formal charge different from the computed one (see CRYSTAL users manual for details) 24

25 Performance of DFT functionals for describing periodic systems quasi-1d TEST SYSTEMS: + SW-carbon nanotubes + SW-chrysotile POSTER N. 4 R. Demichelis, B. Civalleri,Y. Noel and R. Dovesi 3D + aluminium hydroxides semi-ionic, H-containing, HBs, layered and non-layered + brucite semi-ionic, H-containing, no HBs, layered + orthosilicates semi-ionic, H-free + garnets + olivine + Si and Al oxides (+ diamond, silicon pure-covalent, ) work in progress (+ MgO pure-ionic, ) work in progress TEST FUNCTIONALS: SWVN, SPWLSD (LDA) PBE, PW91, PBEsol,WC-PBE,SOGGA-PBE (GGA) B1WC, WC1LYP, B3PW, B3LYP, PBE0 (hybrids) 25

26 Performance of DFT functionals for describing periodic systems Quasi-1D systems SWCNTs: all the tested DFT functionals (LDA, GGA, hybrids) provide reasonable results for structure, energetics and vibrational frequencies. The choice of Hamiltonian is crucial for the description of the band gap. R. Demichelis, Y. Noel, P. D'Arco, R. Dovesi (2009) in preparation. Chrysotile: comparison PBE0 vs B3LYP. Both provide good structures and relative stability. PBE0 structural data are however more accurate. P. D'Arco, Y. Noel, R. Demichelis, R. Dovesi (2009) J. Chem. Phys., accepted 26

27 Performance of DFT functionals for describing periodic systems Aluminosilicates 27

28 Performance of DFT functionals for describing periodic systems ΔV% Aluminosilicates 1- R. Demichelis, B. Civalleri, M Ferrabone and R. Dovesi; Int. J. Quantum Chem. (2009) DOI: /qua R. Demichelis, B. Civalleri, P. DʼArco, R. Dovesi; (2009) in preparation 28

29 29

30 Performance of DFT functionals for describing periodic systems The relative stability of the polymorphs 30

31 Performance of DFT functionals for describing periodic systems In general: The relative stability of the polymorphs 1- R. Demichelis, B. Civalleri, P. DʼArco, R. Dovesi; (2009) in preparation + 2- L. Valenzano, M. Ferrabone, B. Civalleri, R. Dovesi; (2009) in preparation 31

The Schrödinger equation for many-electron systems

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

More information

Electronic Supplementary Information

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

More information

Vibrational frequencies in solids: tools and tricks

Vibrational frequencies in solids: tools and tricks Vibrational frequencies in solids: tools and tricks Roberto Dovesi Gruppo di Chimica Teorica Università di Torino Torino, 4-9 September 2016 This morning 3 lectures: R. Dovesi Generalities on vibrations

More information

AB INITIO MODELING OF ALKALI METAL CHALCOGENIDES USING SOGGA THEORY

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

More information

Applications: Molecular crystals Graphite MgO(001)/CO MIL-53(Al) 2

Applications: Molecular crystals Graphite MgO(001)/CO MIL-53(Al) 2 Bartolomeo Civalleri Voice: Loredana Valenzano B3LYP augmented with an empirical dispersion term (B3LYP-D*) as applied to solids Università di Torino Dipartimento di Chimica IFM & NIS Torino - MSSC2009-10/09/09

More information

Ab initio structure prediction for molecules and solids

Ab initio structure prediction for molecules and solids Ab initio structure prediction for molecules and solids Klaus Doll Max-Planck-Institute for Solid State Research Stuttgart Chemnitz, June/July 2010 Contents structure prediction: 1) global search on potential

More information

Institut Néel Institut Laue Langevin. Introduction to electronic structure calculations

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

More information

DFT: Exchange-Correlation

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,

More information

Advanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University

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

More information

Ab-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials

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

More information

Density Functional Theory - II part

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

More information

Density Func,onal Theory (Chapter 6, Jensen)

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

More information

Minnesota Functional Module Version 1.8

Minnesota Functional Module Version 1.8 1 Minnesota Functional Module Version 1.8 Subroutines for evaluating the M05, M05-2X, M06-L, M06-HF, M06, M06-2X, M08-HX, M08-SO, M11, M11-L, MN12-L, SOGGA, SOGGA11, SOGGA11-X, N12, N12-SX Functionals

More information

Supplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals

Supplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals Supplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals Yuan Ping 1, Giulia Galli 2 and William A. Goddard III 3 1 Joint Center for Artificial Photosynthesis, Lawrence

More information

One-Electron Properties of Solids

One-Electron Properties of Solids One-Electron Properties of Solids Alessandro Erba Università di Torino alessandro.erba@unito.it most slides are courtesy of R. Orlando and B. Civalleri Energy vs Wave-function Energy vs Wave-function Density

More information

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 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 information

Zeitschrift fr physikalische Chemie

Zeitschrift fr physikalische Chemie Zeitschrift fr physikalische Chemie Ab initio study of the vibrational spectrum and related properties of crystalline compounds; the case of CaCO 3 calcite Dr. Loredana Valenzano: Dipartimento di Chimica

More information

The high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations

The high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations J. Phys.: Condens. Matter 8 (1996) 3993 4000. Printed in the UK The high-pressure phase transitions of silicon and gallium nitride: a comparative study of Hartree Fock and density functional calculations

More information

Computational Methods. Chem 561

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

More information

Walter 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 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 information

Pseudopotentials for hybrid density functionals and SCAN

Pseudopotentials 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 information

DFT calculations of NMR indirect spin spin coupling constants

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

More information

One-Electron Hamiltonians

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:

More information

Geometry optimization of solids

Geometry optimization of solids The Minnesota Workshop on ab Initio Modeling in Solid State Chemistry with CRYSTAL Minneapolis, MN(U.S.A.) 9-14 July 2017 7 1 Geometry optimization of solids Bartolomeo Civalleri Dip. di Chimica IFM, Via

More information

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 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 information

Computer Laboratory DFT and Frequencies

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,

More information

International Journal of Quantum Chemistry

International Journal of Quantum Chemistry International Journal of Quantum Chemistry First-principles calculation of second-order elastic constants and equations of state for Lithium Azide, LiN, and Lead Azide, Pb(N ) Journal: International Journal

More information

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 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 information

CRYSTAL: Input / Output

CRYSTAL: Input / Output 1 ISAMS14 International School on Ab initio Modelling of Solids with CRYSTAL14 Regensburg, 20 25/07/2014 CRYSTAL: Input / Output & A. Erba Dipartimento di Chimica and Centre of Excellence NIS, Università

More information

Toward an accurate ab initio estimation of compressibility and thermal expansion of diamond in

Toward an accurate ab initio estimation of compressibility and thermal expansion of diamond in 1 2 Toward an accurate ab initio estimation of compressibility and thermal expansion of diamond in the [0, 3000K] temperature, and [0, 30GPa] pressures ranges, at the hybrid HF/DFT theoretical level 3

More information

First-Principles Vibrational spectroscopy and lattice dynamics of materials in the solid state

First-Principles Vibrational spectroscopy and lattice dynamics of materials in the solid state First-Principles Vibrational spectroscopy and lattice dynamics of materials in the solid state Keith Refson Computational Science and Engineering Department STFC Rutherford Appleton Laboratory First principles

More information

MANUAL Minnesota Functional Module

MANUAL Minnesota Functional Module 1 MANUAL Minnesota Functional Module Version 4.0 Subroutines for evaluating the following exchange-correlation functionals: GAM, M05, M05-2X, M06, M06-2X, M06-HF, M06-L, M08-HX, M08-SO, M11, M11-L, MN12-L,

More information

Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity

Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Torino, Italy, September 4-9, 2016 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it

More information

Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations

Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations, USA, July 9-14, 2017 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it 2017 Outline -

More information

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)

More information

Teoría del Funcional de la Densidad (Density Functional Theory)

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

More information

DFT: Exchange-Correlation

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

More information

Electronic structure of solid FeO at high pressures by quantum Monte Carlo methods

Electronic structure of solid FeO at high pressures by quantum Monte Carlo methods Physics Procedia 3 (2010) 1437 1441 www.elsevier.com/locate/procedia Electronic structure of solid FeO at high pressures by quantum Monte Carlo methods Jindřich Kolorenč a and Lubos Mitas a a Department

More information

Minnesota Functional Module Version 2.0

Minnesota Functional Module Version 2.0 1 Minnesota Functional Module Version 2.0 Subroutines for evaluating the M05, M05-2X, M06-L, M06-HF, M06, M06-2X, M08-HX, M08-SO, M11, M11-L, MN12-L, GAM, MN15-L, MN15, SOGGA, SOGGA11, SOGGA11-X, N12,

More information

An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method

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

More information

Spin densities and related quantities in paramagnetic defects

Spin densities and related quantities in paramagnetic defects Spin densities and related quantities in paramagnetic defects Roberto Orlando Dipartimento di Scienze e Tecnologie Avanzate Università del Piemonte Orientale Via G. Bellini 25/G, Alessandria roberto.orlando@unipmn.it

More information

Module 6 1. Density functional theory

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

More information

Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić

Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824

More information

Oslo node. Highly accurate calculations benchmarking and extrapolations

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

More information

NWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory

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,

More information

The Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal Structure

The Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal Structure Int. J. Mol. Sci. 2004, 5, 93-100 International Journal of Molecular Sciences ISSN 1422-0067 2004 by MDPI www.mdpi.net/ijms/ The Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal

More information

Introduction to Density Functional Theory

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

More information

Spring College on Computational Nanoscience May Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor

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

More information

ELECTRONIC STRUCTURE OF MAGNESIUM OXIDE

ELECTRONIC STRUCTURE OF MAGNESIUM OXIDE Int. J. Chem. Sci.: 8(3), 2010, 1749-1756 ELECTRONIC STRUCTURE OF MAGNESIUM OXIDE P. N. PIYUSH and KANCHAN LATA * Department of Chemistry, B. N. M. V. College, Sahugarh, MADHIPUR (Bihar) INDIA ABSTRACT

More information

Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory

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

More information

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 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

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)

More information

XYZ of ground-state DFT

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

More information

NANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT

NANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT NANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT Quantum Chemistry Laboratory Dipartimento di Scienza dei Materiali Università Milano-Bicocca http://www.mater.unimib.it/utenti/pacchioni

More information

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC

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

More information

Solid State Theory: Band Structure Methods

Solid 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 information

1 Density functional theory (DFT)

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

More information

DFT augmented with an empirical dispersion term as applied to solids Bartolomeo Civalleri

DFT augmented with an empirical dispersion term as applied to solids Bartolomeo Civalleri MSSC2011 - Ab initio Modelling in Solid State Chemistry Torino, 05/09/2011 1 DFT augmented with an empirical dispersion term as applied to solids Bartolomeo Civalleri Theoretical Chemistry Group Department

More information

A theoretical study of stability, electronic, and optical properties of GeC and SnC

A theoretical study of stability, electronic, and optical properties of GeC and SnC JOURNAL OF APPLIED PHYSICS VOLUME 88, NUMBER 11 1 DECEMBER 2000 A theoretical study of stability, electronic, and optical properties of GeC and SnC Ravindra Pandey a) Department of Physics, Michigan Technological

More information

The rôle of surfaces in prebiotic chemistry. Piero Ugliengo University of Torino Dip. Chimica Via P. Giuria, 7 - Torino

The rôle of surfaces in prebiotic chemistry. Piero Ugliengo University of Torino Dip. Chimica Via P. Giuria, 7 - Torino The rôle of surfaces in prebiotic chemistry Piero Ugliengo University of Torino Dip. Chimica Via P. Giuria, 7 - Torino 1 A collaborative work Dip. Chimica - NIS Centre Torino University TO BCN Dep. Quimica

More information

Introduction to Tutorial Sessions

Introduction to Tutorial Sessions Michigan Tech University - Washington State University - University of Torino Ab initio Simulation of Crystalline Systems ASCS2006 September 17-22, 2006 - Spokane, Washington (USA) Introduction to Tutorial

More information

Supporting Information. for

Supporting Information. for Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting Information for Isoreticular Zirconium-Based Metal-Organic Frameworks:

More information

How to generate a pseudopotential with non-linear core corrections

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

More information

Orbital dependent correlation potentials in ab initio density functional theory

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

More information

André Schleife Department of Materials Science and Engineering

André Schleife Department of Materials Science and Engineering André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic

More information

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 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

More information

Note that if DFTTYP=NONE, an ab initio calculation will be performed, rather than density functional theory.

Note that if DFTTYP=NONE, an ab initio calculation will be performed, rather than density functional theory. Input Description $DFT 2-52 ========================================================== $DFT group (relevant if DFTTYP is chosen) (relevant if SCFTYP=RHF,UHF,ROHF) Note that if DFTTYP=NONE, an ab initio

More information

ABC of ground-state DFT

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

More information

Hamiltonians: HF & DFT

Hamiltonians: HF & DFT ISAMS14 International School on Ab initio Modelling of Solids with CRYSTAL14 Regensburg, 0 5/07/014 1 Hamiltonians: HF & DFT and he dreamed that there was a ladder set up on the earth, the top of it reaching

More information

Exchange-Correlation Functional

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

More information

DFT EXERCISES. FELIPE CERVANTES SODI January 2006

DFT EXERCISES. FELIPE CERVANTES SODI January 2006 DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a

More information

Multi-Scale Modeling from First Principles

Multi-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 information

One-electron Hamiltonians: HF & DFT

One-electron Hamiltonians: HF & DFT MSSC06 Ab initio Modelling in Solid State Chemistry Torino, 4 9/09/06 One-electron Hamiltonians: HF & DFT Bartolomeo Civalleri Department of Chemistry IS Centre of xcellence University of Torino bartolomeo.civalleri@unito.it

More information

Intermolecular Forces in Density Functional Theory

Intermolecular Forces in Density Functional Theory Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect

More information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

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

More information

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 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 information

Orbital Density Dependent Functionals

Orbital 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 information

HECToR CSE technical meeting, Oxford Parallel Algorithms for the Materials Modelling code CRYSTAL

HECToR CSE technical meeting, Oxford Parallel Algorithms for the Materials Modelling code CRYSTAL 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 Acknowledgements

More information

The vibrational spectrum of calcite (CaCO 3 ): an ab initio quantum-mechanical calculation

The vibrational spectrum of calcite (CaCO 3 ): an ab initio quantum-mechanical calculation Phys Chem Minerals (2004) 31: 559 564 Ó Springer-Verlag 2004 DOI 10.1007/s00269-004-0418-7 ORIGINAL PAPER M. Prencipe Æ F. Pascale Æ C. M. Zicovich-Wilson V. R. Saunders Æ R. Orlando Æ R. Dovesi The vibrational

More information

Introduction to Computational Chemistry Computational (chemistry education) and/or. (Computational chemistry) education

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

More information

The electronic structure of materials 2 - DFT

The 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 information

A local MP2 periodic study of crystalline argon

A local MP2 periodic study of crystalline argon Journal of Physics: Conference Series A local MP2 periodic study of crystalline argon To cite this article: S Casassa et al 2008 J. Phys.: Conf. Ser. 117 012007 Recent citations - Laplace transformed MP2

More information

Method development at SUNCAT in general

Method development at SUNCAT in general Bayesian error estimation functionals and further method development at SUNCAT GPAW 2016 University of Jyväskylä June 8th, 2016 Johannes Voss vossj@stanford.edu Method development at SUNCAT in general

More information

FIRST-PRINCIPLES MODELING OF OXYGEN INTERACTION WITH ABO 3 -TYPE PEROVSKITE SURFACES

FIRST-PRINCIPLES MODELING OF OXYGEN INTERACTION WITH ABO 3 -TYPE PEROVSKITE SURFACES FIRST-PRINCIPLES MODELING OF OXYGEN INTERACTION WITH ABO 3 -TYPE PEROVSKITE SURFACES S. Piskunov, E.A. Kotomin, Yu.F. Zhukovskii, and V. Alexandrov FMNT 2010, Riga, March 17 Motivation: Oxygen adsorption

More information

DFT with Hybrid Functionals

DFT 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 information

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

More information

Vibrational properties of isotopically enriched materials: the case of calcite Ben Xu, Anna Hirsch, Leeor Kronik, and Kristin M.

Vibrational properties of isotopically enriched materials: the case of calcite Ben Xu, Anna Hirsch, Leeor Kronik, and Kristin M. Electronic Supplementary Material (ESI) for RSC Advances. This journal is The Royal Society of Chemistry 8 Electronic Supporting Information Vibrational properties of isotopically enriched materials: the

More information

Density Functional Theory. Martin Lüders Daresbury Laboratory

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

More information

This is an author version of the contribution published on:

This is an author version of the contribution published on: This is an author version of the contribution published on: M. De La Pierre, C. Carteret, L. Maschio, E. Andre, R. Orlando, R. Dovesi The Raman spectrum of CaCO3 polymorphs calcite and aragonite: A combined

More information

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 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_

More information

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 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

More information

Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids

Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as Benchmarking Tools for Complex Solids Kevin P. Driver, Ph.D. Defense, February 1, 011 Outline 1) Introduction: Quantum mechanics

More information

Spring College on Computational Nanoscience

Spring College on Computational Nanoscience 2145-29 Spring College on Computational Nanoscience 17-28 May 2010 At the Fifth Rung of Jacob's Ladder: A Discussion of Exact Exchange plus Local- and Nonlocal-density Approximations to the Correlation

More information

THE PERFORMANCE OF DENSITY FUNCTIONALS WITH RESPECT TO THE CORRELATION CONSISTENT BASIS SETS. Xuelin Wang, B.S., M.S

THE PERFORMANCE OF DENSITY FUNCTIONALS WITH RESPECT TO THE CORRELATION CONSISTENT BASIS SETS. Xuelin Wang, B.S., M.S THE PERFORMANCE OF DENSITY FUNCTIONALS WITH RESPECT TO THE CORRELATION CONSISTENT BASIS SETS Xuelin Wang, B.S., M.S Dissertation Prepared for the Degree of DOCTOR OF PHILOSOPHY UNIVERSITY OF NORTH TEXAS

More information

A BADER S TOPOLOGICAL APPROACH FOR THE CHARACTERIZATION OF PRESSURE INDUCED PHASE TRANSITIONS

A BADER S TOPOLOGICAL APPROACH FOR THE CHARACTERIZATION OF PRESSURE INDUCED PHASE TRANSITIONS A BADER S TOPOLOGICAL APPROACH FOR THE CHARACTERIZATION OF PRESSURE INDUCED PHASE TRANSITIONS FILIPPO PARISI Dipartimento di Matematica e Geoscienze, Università di Trieste, Via E. Weiss 2, 34128 Trieste,

More information

Band calculations: Theory and Applications

Band calculations: Theory and Applications Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation

More information

Pseudopotentials: design, testing, typical errors

Pseudopotentials: 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 information

arxiv:physics/ v1 [physics.comp-ph] 19 Dec 2005

arxiv:physics/ v1 [physics.comp-ph] 19 Dec 2005 Analytical Hartree-Fock gradients with respect to the cell parameter: systems periodic in one and two dimensions K. Doll Institut für Mathematische Physik, TU Braunschweig, Mendelssohnstraße 3, D-38106

More information

Short Course on Density Functional Theory and Applications III. Implementations

Short Course on Density Functional Theory and Applications III. Implementations Short Course on Density Functional Theory and Applications III. Implementations Samuel B. Trickey Sept. 2008 Quantum Theory Project Dept. of Physics and Dept. of Chemistry trickey@qtp.ufl.edu KS E xc and

More information