Self-consistent implementation of meta-gga exchange-correlation functionals within the ONETEP linear-scaling DFT code

Size: px
Start display at page:

Download "Self-consistent implementation of meta-gga exchange-correlation functionals within the ONETEP linear-scaling DFT code"

Transcription

1 Self-consistent implementation of meta-gga exchange-correlation functionals within the ONETEP linear-scaling DFT code James C. Womack, Chris-Kriton Skylaris Chemistry, Faculty of Natural & Environmental Sciences, University of Southampton, Highfield, Southampton, SO17 1BJ, UK RSC Theoretical Chemistry Group Conference, Nottingham, 2016

2 Overview Background & Motivation Theory Linear-scaling DFT ONETEP formalism Density functional approximations Meta-GGA functionals Approaches to self-consistency Why? Implementation General meta-gga framework in ONETEP Results & Conclusions Numerical validation Demonstration of linear-scaling Specific functionals implemented 2

3 Linear-scaling density functional theory Orthogonalize KS orbitals extending over entire system Conventional Kohn-Sham DFT Cubically scaling computational cost Linear-scaling density matrix DFT The density matrix is nearsighted For insulators:* Exploit this to obtain linear scaling computational cost *see e.g. S. Ismail-Beigi and T.A. Arias, Phys. Rev. Lett. 82, 2127 (1999) 3

4 Approximate exchange-correlation functionals No practical exact form is known Ideally, these should be Accurate Computationally efficient Hierarchical series New ingredients on each rung double-hybrid non-local Jacob's ladder of approximations Heaven of chemical accuracy Approximations are required hybrid meta-gga (semi-)local GGA Higher rungs are more flexible Which is the best tool for the job? LDA Hartree world J.P. Perdew and K. Schmidt, AIP Conf. Proc. 577, 1 (2001) 4

5 Motivation: The meta-gga Goldilocks zone Non-local ingredients (e.g. exact exchange) pose a challenge to linear scaling* Increasing accuracy Meta-GGAs sit in a Goldilocks zone Well-balanced cost vs. accuracy Simple semi-local form Add kinetic energy density, τ More flexible functional forms Greater constraint satisfaction** Modern meta-ggas offer impressive accuracy e.g. B97M-V, M15-L, SCAN... Goldilocks zone * Linear scaling exact exchange has been achieved in ONETEP but with a relatively large prefactor: J. Dziedzic, Q. Hill, and C.-K. Skylaris, J. Chem. Phys. 139, (2013) Decreasing computational cost ** For more on constraint satifaction: J.P. Perdew, A. Ruzsinszky, J. Tao, V.N. Staroverov, G.E. Scuseria, and G.I. Csonka, J. Chem. Phys. 123, (2005) 5

6 Order-N Electronic Total Energy Package NGWFs Non-orthogonal generalized Wannier functions Strictly localized orbitals Self-consistently optimized Periodic cardinal sinc psinc basis Density kernel Sparse matrix Truncated above cut-off length, i.e. C.-K. Skylaris, P.D. Haynes, A.A. Mostofi, and M.C. Payne, J. Chem. Phys. 122, (2005) 6

7 Order-N Electronic Total Energy Package Total energy minimization Via two nested loops, i.e. Initial guess Optimize K Converged? Yes No Optimize NGWFs Psinc basis No Converged? Orthogonal Bandwidth-limited delta functions on a regular real-space grid Finite sum over plane waves, i.e. Yes Energy minimized 7

8 Self-consistent meta-gga evaluation Unknown dependence of τ on the density prevents direct functional derivative evaluation There are three main approaches for overcoming this: 1. Give up! 2. Optimized effective potential Perform meta-gga calculations nonself-consistently Obtain a local, multiplicative potential which satisfies Use orbitals optimized using a non-τdependent functional (LDA, GGA) Yields surprisingly good results Example: J.P. Perdew, S. Kurth, A. Zupan, and P. Blaha, Phys. Rev. Lett. 82, 2544 (1999). by solving the OEP integral equation Formidable mathematical machinery (even with approximations) Example: A.V. Arbuznikov and M. Kaupp, Chem. Phys. Lett. 381, 495 (2003). 8

9 3. FDO approximation Functional Derivatives (of τ-dependent functionals) with respect to the Orbitals* Use the functional derivative chain rule to obtain The functional derivative allows us to derive the following relationship * AKA NNH approach after the creators: R. Neumann, R.H. Nobes, and N.C. Handy, Mol. Phys. 87, 1 (1996) This functional derivative is approximate, since it neglects the dependence of the KS orbitals on each other via the density, i.e. For a detailed account, see: F. Zahariev, S.S. Leang, and M.S. Gordon, J. Chem. Phys. 138, (2013) 9

10 3. FDO approximation We know the functional dependence of τ on the KS orbitals FDO Non-multiplicative, orbital-specific potential and can therefore evaluate and obtain our FDO XC potential using This approach has the advantages of Theoretical simplicity Computational simplicity Wide use in other code (i.e. well-tested) A.V. Arbuznikov and M. Kaupp, Chem. Phys. Lett. 381, 495 (2003). 10

11 Implementation overview Kinetic energy density Exchange-correlation potential Self-consistent meta-gga functionals in ONETEP Integral over τ-dependent potential Energy gradients Functional forms 11

12 Kinetic energy density We require τ in terms of the NGWFs and density kernel, i.e. In practice, we evaluate τ via a multi-step process: 1 Evaluate intermediate quantity, summing over NGWFs which overlap with α 2 Form part of τ local to each α NGWF 3 Deposit each α-component in simulation cell Based on the approach used to evaluate the charge density in ONETEP: C.-K. Skylaris, P.D. Haynes, A.A. Mostofi, and M.C. Payne, Phys. Stat. Sol. (B) 243, 973 (2006) 12

13 Exchange-correlation potential We require an XC potential in terms of NGWFs, not KS orbitals: FDO Expand KS orbital in terms of NGWFs Per-NGWF exchange-correlation potential 13

14 Integral over τ-dependent potential Hamiltonian matrix elements contain integrals over the τ-dependent potential: Gradient of this function avoided using integration-by-parts Form originally used by NNH and in later implementations No problem for ONETEP! Apply divergence to vector field in reciprocal space Form used in ONETEP Integrand is zero outside of localization region of α 14

15 Energy gradients Gradient expressions are required for ONETEP's energy minimization scheme. Only minor changes and additions to the existing expressions are needed. New integral type Q matrix depends on New τ-dependent potential term 15

16 Functionals PKZB An early meta-gga (1999) Modern meta-gga (2015) 1 empirical parameter Combinatorially designed form Evolved into TPSS (2003) Empirically fitted to very large data set Simple to implement Based on PBE, but with τ Incorporates non-local VV10 functional Widely available B97M-V Stable implementations to compare and test against Excellent accuracy Comparable to popular hybrid functionals Broadly applicable Inconsistent accuracy Superseded by modern meta-ggas J.P. Perdew, S. Kurth, A. Zupan, and P. Blaha, Phys. Rev. Lett. 82, 2544 (1999) Available in ONETEP Code provided by N. Mardirossian N. Mardirossian and M. Head-Gordon, J. Chem. Phys. 142, (2015) 16

17 Aside: How good is B97M-V? Thermochemical data Non-covalent data 2460 data points RMSD Hybrid M06-2X 3.21 ωb97x-v 3.60 ωb97x-d 3.61 B97M-V 3.93 B97-D M11 RMSD Hybrid B97M-V 0.22 ωb97x-v 0.32 M06-L 0.42 B97-D ωb97x-d M M M B3LYP-D M06-L TPSS-D VV PBE-D PBE-D VV Zero Zero Performs well against hybrid functionals! M06-2X NLC is not sufficient for good performance B3LYP-D3 in either set TPSS-D aqz basis set atz basis set RMSDs in kcal/mol for thermochemical and non-covalent datasets (combined training and test data). Based on Table V of N. Mardirossian and M. Head-Gordon, J. Chem. Phys. 142, (2015). 17

18 Results I: Numerical testing Water dimer Å Comparison of relative energies for B97M-V QChem (Gaussian basis, BSSE corrected) ONETEP (NGWFs/psincs) B97M-rV* Complicated by several possible sources of error Pseudopotential vs. all-electron Periodic vs. open boundary conditions Binding energy Program Basis PBE B97M-rV E(ONETEP) E(QChem) PBE B97M-rV - - QChem aug-pc ONETEP 900eV / 8 bohr ONETEP 900eV / 12 bohr ONETEP 900eV / 15 bohr ONETEP 1200eV / 8 bohr ONETEP 1200eV / 12 bohr ONETEP 1200eV / 15 bohr Water dimer binding energies in kcal/mol. QChem results provided by N. Mardirossian (UC Berkeley). * rv refers to revised VV10: R. Sabatini, T. Gorni, and S. de Gironcoli, Phys. Rev. B 87, (2013) 18

19 Results I: Numerical testing Water dimer The 8 bohr results seem too good to be true! Further ONETEP calculation details: Why are the results worse for larger NGWF radii? Use ONETEP's extended NGWF functionality. Possibly a favourable cancellation of errors (unfortunately). Basis 24Å cubic simulation cell Norm-conserving pseudopotentials* E(ONETEP) E(QChem) E(ONETEP,local) E(ONETEP,extended) PBE PBE B97M-rV B97M-rV 900eV / 8.0 bohr eV / 12.0 bohr eV / 15.0 bohr eV / Extended eV / 8.0 bohr eV / 12.0 bohr eV / 15.0 bohr eV / Extended Energy differences in kcal/mol. QChem results provide by N. Mardirossian (UC Berkeley). *Rappe-Bennett pseudopotential library: 19

20 Results II: Computational efficiency atom structure Iridis 4 supercomputer 128 MPI processes 4 OpenMP threads Structure provided by authors of J.T. Berryman, S.E. Radford, and S.A. Harris, Biophys. J. 97, 1 (2009) 20

21 Conclusions Self-consistent meta-gga functional evaluation in ONETEP Theoretical and computational framework implemented and tested Further functional forms can easily be added Meta-GGAs can be applied to very large systems Provisional numerical agreement with other implementations Further testing and validation is ongoing Linear-scaling computational cost Self-consistent evaluation of meta-gga functionals Demonstrated for systems with 1000s of atoms 21

22 The future Extensions to the framework Support for projector augmented waves (PAW) Support for hybrid meta-gga functionals using existing exact exchange in ONETEP* Additional meta-gga functional forms TPSS, SCAN, M15-L are of particular interest Public availability Self-consistent meta-gga support will be made available in a future version of ONETEP *J. Dziedzic, Q. Hill, and C.-K. Skylaris, J. Chem. Phys. 139, (2013) 22

23 Acknowledgements People Martin Head-Gordon (Berkeley) Narbe Mardirossian (Berkeley) Jacek Dziedzic (Southampton) Compute resources Funding University of Southampton IRIDIS HPC facility Engineering and Physical Sciences Research Council University of Southampton 23

24 Original artwork by Lauren Womack

Self-consistent implementation of meta-gga exchange-correlation functionals within the ONETEP linear-scaling DFT code

Self-consistent implementation of meta-gga exchange-correlation functionals within the ONETEP linear-scaling DFT code Self-consistent implementation of meta-gga exchange-correlation functionals within the linear-scaling DFT code, August 2016 James C. Womack,1 Narbe Mardirossian,2 Martin Head-Gordon,2,3 Chris-Kriton Skylaris1

More information

Self-consistent implementation of meta-gga functionals for the ONETEP linear-scaling electronic structure package

Self-consistent implementation of meta-gga functionals for the ONETEP linear-scaling electronic structure package Self-consistent implementation of meta-gga functionals for the ONETEP linear-scaling electronic structure package James C. Womack, 1 Narbe Mardirossian, 2 Martin Head-Gordon, 2, 3 and Chris-Kriton Skylaris

More information

Electronic energy optimisation in ONETEP

Electronic energy optimisation in ONETEP Electronic energy optimisation in ONETEP Chris-Kriton Skylaris cks@soton.ac.uk 1 Outline 1. Kohn-Sham calculations Direct energy minimisation versus density mixing 2. ONETEP scheme: optimise both the density

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

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

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

THE JOURNAL OF CHEMICAL PHYSICS 127,

THE JOURNAL OF CHEMICAL PHYSICS 127, THE JOUNAL OF CHEMICAL PHYSICS 127, 214103 2007 Avoiding singularity problems associated with meta-gga generalized gradient approximation exchange and correlation functionals containing the kinetic energy

More information

Density functional theory in solution: Implementing an implicit solvent model for CASTEP and ONETEP

Density functional theory in solution: Implementing an implicit solvent model for CASTEP and ONETEP Density functional theory in solution: Implementing an implicit solvent model for CASTEP and ONETEP, November 2016 James C. Womack University of Southampton, Southampton, UK Overview Background The model

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

Rapid and precise thermochemical calculations by quantum chemical methods

Rapid and precise thermochemical calculations by quantum chemical methods Rapid and precise thermochemical calculations by quantum chemical methods Ph.D. thesis By: Adrienn Ruzsinszky Supervisor: Dr. Gábor Csonka Budapest University of Technology and Economics Department of

More information

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole ONETEP PB/SA: Application to G-Quadruplex DNA Stability Danny Cole Introduction Historical overview of structure and free energy calculation of complex molecules using molecular mechanics and continuum

More information

Energies of Organic Molecules and Atoms in Density Functional Theory

Energies of Organic Molecules and Atoms in Density Functional Theory Energies of Organic Molecules and Atoms in Density Functional Theory Gábor I. Csonka,*,1 Adrienn Ruzsinszky, 1 Jianmin Tao, 2 and John P. Perdew,*,2 1 Department of Inorganic Chemistry, Budapest University

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

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

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

A FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY

A FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY A FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY JOHN P. PERDEW PHYSICS & CHEMISTRY, TEMPLE UNIVERSITY PHILADELPHIA, PENNSYLVANIA, USA SUPPORTED BY THE U.S. DEPARTMENT OF ENERGY, EFRC

More information

The Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus

The Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,

More information

Introduction to DFT and its Application to Defects in Semiconductors

Introduction to DFT and its Application to Defects in Semiconductors Introduction to DFT and its Application to Defects in Semiconductors Noa Marom Physics and Engineering Physics Tulane University New Orleans The Future: Computer-Aided Materials Design Can access the space

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

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

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

Modified Becke-Johnson (mbj) exchange potential

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

More information

Introducing ONETEP: Linear-scaling density functional simulations on parallel computers

Introducing ONETEP: Linear-scaling density functional simulations on parallel computers THE JOURNAL OF CHEMICAL PHYSICS 122, 084119 2005 Introducing ONETEP: Linear-scaling density functional simulations on parallel computers Chris-Kriton Skylaris, a) Peter D. Haynes, Arash A. Mostofi, and

More information

Role of van der Waals Interactions in Physics, Chemistry, and Biology

Role of van der Waals Interactions in Physics, Chemistry, and Biology Role of van der Waals Interactions in Physics, Chemistry, and Biology How can we describe vdw forces in materials accurately? Failure of DFT Approximations for (Long-Range) Van der Waals Interactions 1

More information

Kohn Sham density functional theory [1 3] is. Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW.

Kohn Sham density functional theory [1 3] is. Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW. Role of the Exchange Correlation Energy: Nature s Glue STEFAN KURTH, JOHN P. PERDEW Department of Physics and Quantum Theory Group, Tulane University, New Orleans, Louisiana 70118 Received 11 March 1999;

More information

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 Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,

More 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

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

More accurate generalized gradient approximation for solids

More accurate generalized gradient approximation for solids More accurate generalized gradient approximation for solids Zhigang Wu* and R. E. Cohen Carnegie Institution of Washington, Washington, DC 20015, USA Received 24 February 2006; published 20 June 2006 We

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

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

Magnetism in transition metal oxides by post-dft methods

Magnetism in transition metal oxides by post-dft methods Magnetism in transition metal oxides by post-dft methods Cesare Franchini Faculty of Physics & Center for Computational Materials Science University of Vienna, Austria Workshop on Magnetism in Complex

More 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

MODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO

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

FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY

FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY C.A. Madu and B.N Onwuagba Department of Physics, Federal University of Technology Owerri, Nigeria

More information

Tutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30

Tutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30 Tutorial I: IQ MOL and Basic DFT and MP2 Calculations Q-Chem User Workshop, Denver March 21, 2015 1 / 30 2 / 30 Introduction to IQMOL DFT and MP2 Calculations 3 / 30 IQMOL and Q-CHEM IQMOL is an open-source

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

Comparison of various abinitio codes used in periodic calculations

Comparison of various abinitio codes used in periodic calculations Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,

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

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

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

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

Key concepts in Density Functional Theory (II)

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

Two implementations of the Projector Augmented Wave (PAW) formalism

Two implementations of the Projector Augmented Wave (PAW) formalism Introduction The tools available for detailed first-principles studies of materials have benefited enormously from the development of several international collaborations engaged in developing open source

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

Density Functional Theory

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

More information

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

A self-interaction free non-local correlation functional derived from exact criteria

A self-interaction free non-local correlation functional derived from exact criteria A self-interaction free non-local correlation functional derived from exact criteria Tobias Schmidt and Stephan Kümmel Theoretical Physics IV University of Bayreuth Germany - International Workshop DFT

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

Orbital currents in the Colle-Salvetti correlation energy functional and the degeneracy problem. Abstract

Orbital currents in the Colle-Salvetti correlation energy functional and the degeneracy problem. Abstract Orbital currents in the Colle-Salvetti correlation energy functional and the degeneracy problem S. Pittalis 1, S. Kurth 1, S. Sharma 1,2 and E.K.U. Gross 1 1 Institut für Theoretische Physik, Freie Universität

More information

ASSESSMENT OF DFT METHODS FOR SOLIDS

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

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

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

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

Performance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26

Performance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26 Performance of hybrid density functional methods, screened exchange and EXX-OEP methods in the PAW approach Georg Kresse, J Paier, R Hirschl, M Marsmann Institut für Materialphysik and Centre for Computational

More information

The Linearized Augmented Planewave (LAPW) Method

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

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

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

Density functional theory in the solid state

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

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

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

Understanding electron correlation energy through density functional theory

Understanding electron correlation energy through density functional theory 1 Understanding electron correlation energy through density functional theory by Teepanis Chachiyo, 1,2,* and Hathaithip Chachiyo 2 1 Department of Physics, Faculty of Science, Naresuan University, Phitsanulok

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

Practical Guide to Density Functional Theory (DFT)

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

Finite-Temperature Hartree-Fock Exchange and Exchange- Correlation Free Energy Functionals. Travis Sjostrom. IPAM 2012 Workshop IV

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

More information

Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water

Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical

More information

Computational Modeling Software and their applications

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

More information

The Projector Augmented Wave method

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

Computational Chemistry I

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

More information

Integrated Computational Materials Engineering Education

Integrated Computational Materials Engineering Education Integrated Computational Materials Engineering Education Lecture on Density Functional Theory An Introduction Mark Asta Dept. of Materials Science and Engineering, University of California, Berkeley &

More information

Table of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures

Table of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures Table of Contents Table of Contents Spin-orbit splitting of semiconductor band structures Relavistic effects in Kohn-Sham DFT Silicon band splitting with ATK-DFT LSDA initial guess for the ground state

More information

Self-averaging stochastic Kohn-Sham density functional theory

Self-averaging stochastic Kohn-Sham density functional theory Self-averaging stochastic Kohn-Sham density functional theory Roi Baer *, Daniel Neuhauser and Eran Rabani * Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, Hebrew University, Jerusalem

More information

Density Functional Theory for Electrons in Materials

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

More information

Comparison of exchange-correlation functionals: from LDA to GGA and beyond

Comparison of exchange-correlation functionals: from LDA to GGA and beyond Comparison of ehange-correlation functionals: from LDA to GGA and beyond Martin Fuchs Fritz-Haber-Institut der MPG, Berlin, Germany Density-Functional Theory Calculations for Modeling Materials and Bio-Molecular

More information

Density Functional Theory (DFT)

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

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

Introduction to Density Functional Theory

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

More information

Binding energy of bilayer graphene and Electronic properties of oligoynes

Binding energy of bilayer graphene and Electronic properties of oligoynes Binding energy of bilayer graphene and Electronic properties of oligoynes E. Mostaani and N. Drummond Thursday, 31 July 2014 Van der Waals interaction Important contributions to the description of binding

More information

The Quantum ESPRESSO Software Distribution

The Quantum ESPRESSO Software Distribution The Quantum ESPRESSO Software Distribution The DEMOCRITOS center of Italian INFM is dedicated to atomistic simulations of materials, with a strong emphasis on the development of high-quality scientific

More 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

Norm-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 Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized

More information

Accurate energies calculated by empirical corrections to the local spin density approximation

Accurate energies calculated by empirical corrections to the local spin density approximation Theor Chem Acc (2004) 112: 1 6 DOI 10.1007/s00214-003-0555-y Accurate energies calculated by empirical corrections to the local spin density approximation Rene Fournier, Luba Boroukhovskaia Department

More information

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

Explaining the apparent arbitrariness of the LDA-1/2 self-energy. correction method applied to purely covalent systems

Explaining the apparent arbitrariness of the LDA-1/2 self-energy. correction method applied to purely covalent systems Explaining the apparent arbitrariness of the LDA-1/2 self-energy correction method applied to purely covalent systems Kan-Hao Xue, 1,2 Leonardo R. C. Fonseca, 3 and Xiang-Shui Miao 1,2 1 School of Optical

More information

Adiabatic connection from accurate wave-function calculations

Adiabatic connection from accurate wave-function calculations JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 12 22 MARCH 2000 Adiabatic connection from accurate wave-function calculations Derek Frydel and William M. Terilla Department of Chemistry, Rutgers University,

More information

CHEM6085: Density Functional Theory

CHEM6085: Density Functional Theory Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron

More information

Fourteen Easy Lessons in Density Functional Theory

Fourteen Easy Lessons in Density Functional Theory Fourteen Easy Lessons in Density Functional Theory JOHN P. PERDEW, ADRIENN RUSINSKY Department of Physics and Quantum Theory Group, Tulane University, New Orleans, LA 70118 Received 1 April 2010; accepted

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

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

Wavelets for density functional calculations: Four families and three. applications

Wavelets for density functional calculations: Four families and three. applications Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Interpolating

More information

College of Science, Xi an University of Science and Technology, Xi an *Corresponding author

College of Science, Xi an University of Science and Technology, Xi an *Corresponding author 2016 International Conference on Advanced Manufacture Technology and Industrial Application (AMTIA 2016) ISBN: 978-1-60595-387-8 The Study of Coordination Adsorption Effect that CO Adsorption on 4H-SiC

More information

Understanding the structure and reactivity of NiCu nanoparticles: an atomistic modeling

Understanding the structure and reactivity of NiCu nanoparticles: an atomistic modeling Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2017 Understanding the structure and reactivity of NiCu nanoparticles: an atomistic

More information

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 17 Nov 2006

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 17 Nov 2006 Simple implementation of complex functionals: scaled selfconsistency Matheus P. Lima, Luana S. Pedroza, Antonio J. R. da Silva and A. Fazzio Instituto de Física, Universidade de São Paulo, São Paulo, Brazil

More information

Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals

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

First-principles calculations of structural, electronic and optical properties of HfZn 2

First-principles calculations of structural, electronic and optical properties of HfZn 2 ~ 1 ~ First-principles calculations of structural, electronic and optical properties of HfZn 2 Md. Atikur Rahman *1, Md. Afjalur Rahman 2, Md. Zahidur Rahaman 3 1, 2, 3 Department of Physics, Pabna University

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