Open-shell and Magnetic Systems with CRYSTAL Tools

Size: px
Start display at page:

Download "Open-shell and Magnetic Systems with CRYSTAL Tools"

Transcription

1 Open-shell and Magnetic Systems with CRYSTAL Tools Klaus Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D Stuttgart, Germany MSSC 2018, Turin, September 2018

2 Contents Introduction Hints for the input Examples

3 Magnetism Considered here: Ferro-, Antiferro-, Ferrimagnetism typically in periodic systems with: d- or f-electrons but also with p-electrons defects Molecules: molecules with unpaired electrons Questions: Strength : exchange interaction J magnetic moments, spin and charge densities, electrical field gradient dependence on pressure

4 Types of exchange interaction, I Direct (between neighbouring atoms) Indirect (superexchange): via bridging ion here: antiferromagnetic

5 Types of exchange interaction, II Double exchange: ferromagnetic exchange interaction which occurs because the magnetic ion can show mixed valency Indirect exchange: around defect, mediated by conduction electrons

6 Types of exchange interaction, III Itinerant magnetism: 3d metals (Fe, Co, Ni) with partially filled shells different description with models which take into account kinetic energy

7 Computational methods for periodic systems and analogously KS, UKS (Kohn-Sham)

8 An atom: oxygen Atomic configuration: 1s2 2s2 2p4 Ground state: 3P2 (Hund s rules) 1) maximise spin: 2 unpaired electrons => spin 1 (2*ħ/2), 2S+1 =3 2) maximise L: (3 up-electrons: Lup=0, 1 down electron, Ldown=1) 1 => P J = L+ S => 2 3) more than half full p-shell: Eigenvalues: S2 : S(S+1) =2 L2 : L(L+1) = 2

9 O with Molpro ( Input file, part 1: gprint,basis,orbital basis=vdz geometry={o} {rhf,throrb=0; occ,2,1,1,0,1;closed,2,0,0,0,1 wf,8,4,2 } i=1 meth(i) = 'RHF' e(i)= energy {multi;noextra occ,2,1,1,0,1;closed,2,0,0,0,1;frozen,2,0, 0,0,1 wf,8,4,2 expec2,lxx,lyy,lzz} i=i+1 meth(i) = 'MULTI - 1 state' e(i)= energy part 2 (copy+paste below part 1): {multi occ,2,1,1,0,1;closed,2,0,0,0,0;frozen,0,0, 0,0,0 wf,8,4,2; wf,8,6,2; wf,8,7,2; expec2,lxx,lyy,lzz} i=i+1 meth(i) = 'MULTI-stateaveraged' e(i)= energy {uhf,throrb=0; occ,2,1,1,0,1;closed,2,0,0,0,1 wf,8,4,2 } i=i+1 meth(i) = 'UHF' e(i)= energy table,meth,e;

10 O with Molpro Output: Energies: METH RHF MULTI - 1 state MULTI-stateaveraged UHF E RHF: symmetry-broken: 2 pz electrons, 1 px, 1 py electron Not spherically symmetric, no L2 eigenstate! <1.4 L**2 1.4> MCSCF state-averaged: 3 states, one with 2 px electrons, one with 2 py, one with 2pz!MCSCF expec <1.4 L**2 1.4> !MCSCF expec <1.6 L**2 1.6> !MCSCF expec <1.7 L**2 1.7> => spherically symmetric! UHF: symmetry-broken, not an eigenstate of S2 EXPECTATION VALUE OF S**2:

11 O with CRYSTAL Input file, part 1: Oxygen Atom MOLECULE END part 2 (copy+paste below part 1): END UHF ATOMSPIN 1 11 SPINLOCK 23 MAXCYCLE 60 FMIXING 70 PPAN END spin-polarised atom with spin initial spin state (2 more up electrons than down electrons, fix for 3 cycles)

12 O with CRYSTAL Output: Energies: ATOMIC WAVEFUNCTION(S) NUCLEAR CHARGE 8.0 SYMMETRY SPECIES S P N. ELECTRONS 8.0 NUMBER OF PRIMITIVE GTOS 19 5 NUMBER OF CONTRACTED GTOS 3 2 NUMBER OF CLOSED SHELLS 2 0 OPEN SHELL OCCUPATION 0 4 ZNUC SCFIT TOTAL HF ENERGY KINETIC ENERGY VIRIAL THEOREM ACCURACY E E E E-06 Identical to Molpro, state-averaged ( ) == SCF ENDED - CONVERGENCE ON ENERGY Identical to Molpro, UHF ( ) E(AU) E+01 CYCLES 7

13 Exchange interaction: mathematically Triplet: Singlet: [ ](Φa (1) Φb (2) Φb (1)Φ a (2)) [ ](Φa (1) Φb (2)+Φb (1) Φa (2)) compute energy difference Esinglet-Etriplet fit models: 1 S 2) Heisenberg: H= J ( S Ising: H= J (S1 S 2)

14 Calculations for solids: only Sz Eigenstate feasible, fit Ising model ferromagnet versus antiferromagnet: energy difference: njsz2 n: neighbors energy scale: ~J (~TN) order of mev total energy: H-atom 13.6 ev ; Ni atom ~ ev => code must be accurate (little numerical noise)

15 Keywords general: look at test cases, tutorials for examples Always: UHF or DFT SPIN Geometry: often a SUPERCEL is necessary Usually: ATOMSPIN Often: SPINLOCK (keep number of cycles smaller than number necessary to converge otherwise artificial solution) Sometimes: EIGSHIFT LEVSHIFT MODISYMM

16 NiO: band gap with various functionals LDA: K. Terakura, T. Oguchi, A. R. Williams, J. Kübler, Phys. Rev. B 30, 4734 (1984) B3LYP: I. Moreira et al, Phys. Rev. B 65, (2002) Hartree-Fock: M. D. Towler et al, Phys. Rev. B 50, 5041 (1994)

17 Band gaps: difficult with density functional theory GaAs NiO ZnO TiO2 LDA B3LYP HF Expt. (ev) Hartree-Fock: correct self-interaction, but Coulomb-repulsion unscreened LDA: self-interaction error (e.g. hydrogen: 1 electron interacts with itself) B3LYP: mixed Hartree-Fock and LDA, is in between Practical approach: use several methods, compare, this gives approximate result Luckily not all properties are so difficult to compute! Geometries, frequencies, charges usually match well occupied bands, unoccupied bands match well, but not the gap

18 Superexchange calculation for NiO FM: ev AF2: ev energy difference: ev = 3JS2 => J = 5.5 mev exp.: J=20 mev M. D. Towler, N.L. Allan, N.M.Harrison, V.R. Saunders, W. C. Mackrodt, E. Aprà, Phys. Rev. B 50, 5041 (1994) today: ~ 1 minute CPU time for NiO

19 Superexchange calculation for KMnF3 J=0.1 mev J ~ d-14 exp.:j=0.3 mev d -12 +/- 2 J. M. Ricart, R. Dovesi, C. Roetti, V.R. Saunders, Phys. Rev. B 52, 2381 (1995) R. Dovesi, F. Freyria-Fava, C. Roetti, V. R. Saunders, Faraday Discuss. 106, 173 (1997)

20 NiO: spin density plot Hartree-Fock magn. moment J Ni O mev B3LYP mev LDA mev exp.: -20 mev

21 Some more examples Hartree-Fock B3LYP LDA exp. NiO -5.5 mev -29 mev -94 mev -20 mev KNiF3-2.6 mev -15 mev La2CuO4-36 mev -130 mev MnF2 0.2 mev 0.05 mev mev -0.3 mev J1 J2-9 mev systematic deviation; LDA too delocalized; HF too localized prediction possible for 180º angle

22 Fe(pyrimidine)2Cl2:a molecule based magnet canted =>weakly ferromagnetic Orientation of the magnetic moment from Mössbauer spectroscopy and calculation (electrical field gradient) R. Feyerherm, A. Loose, T. Ishida, T. Nogami, J. Kreitlow, D. Baabe, F. J. Litterst, S. Süllow, H.-H. Klauss, K. Doll, Phys. Rev. B 69, (2004)

23 Exchange interaction FePM2Cl2 NiPM2Cl2 J exp. th. (B3LYP) mev mev mev -0.6 mev Increase of J when compressing by 5%: ~15% (Experiment and theory)

24 Bulk modulus exp.: very soft: 15 GPa calculated: 18 GPa What happens under pressure? a;c Fe-N Fe-Cl N-C C-H ; ; C-C (Å) A. U. B. Wolter, H.-H. Klauss, F. J. Litterst, T. Burghardt, A. Eichler, R. Feyerherm, S. Süllow, Polyhedron 22, 2139 (2003) J. Kreitlow, D. Menzel, A. U. B. Wolter, J. Schoenes, S. Süllow, K. Doll, Phys. Rev. B 72, (2005)

25 Rb4O6: magnetism in p-shell most favorable solution: Rb4(O2-)2 (O22-) O22- (peroxide) nonmagnetic, O2- (superoxide) magnetic moment ~1 μb 1 peroxide, bond length 1.56 Å 2 superoxides, bond length 1.36 Å symmetry becomes I-42d (in case of ferromagnetic order) J. Winterlik, G. H. Fecher, C. A. Jenkins, C. Felser, C. Mühle, K. Doll, M. Jansen, L. M. Sandratskii, J. Kübler, Phys. Rev. Lett. 102, (2009)

26 Raman spectra Mn3O4: Raman Tetragonal at room temperature Becomes ferrimagnetic below 42 K red: Oxygen green: Mn, octahedral site (6-fold coordinated) sky blue: Mn, tetrahedral site (4-fold coordinated) simulated as ferromagnetic (all Mn up spin) T. Larbi, K. Doll, T. Manoubi, J. Alloys Compd. 688, 692 (2016)

27 Raman spectra Mn3O4: ferrimagnetic red: Oxygen green: Mn, octahedral site, up spin dark blue: Mn, octahedral site, down spin sky blue: Mn, tetrahedral site, up spin symmetry reduced in simulation, becomes orthorhombic magnetic states discussed in A. Chartier, Ph. D Arco, R. Dovesi, V. R. Saunders, Phys. Rev. B 60, (1999)

28 Symmetry change of modes Mn3O4: Raman Space group, high symmetry: 141 Point group: Oh (4/mmm) Point group, low symmetry: mmm

29 Example for f-elements: GdN LDA: metallic, very good agreement with previous calculations, e.g. P. Larson and W. R. L. Lambrecht, Phys. Rev. B 74, (2006) B3LYP: only points of band structure touch Fermi surface K. Doll, J. Phys.: Cond. Matt. 20, (2008)

30 Conclusion Methods: UHF, spin-polarised DFT magnetism: systematic deviation for exchange interaction J, hybrid functionals give relatively good agreement with experiment vibrational spectra well reproduced unusual magnetism (p-electrons involved)

Gaussian Basis Sets for Solid-State Calculations

Gaussian Basis Sets for Solid-State Calculations Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction

More information

Luigi Paolasini

Luigi Paolasini Luigi Paolasini paolasini@esrf.fr LECTURE 4: MAGNETIC INTERACTIONS - Dipole vs exchange magnetic interactions. - Direct and indirect exchange interactions. - Anisotropic exchange interactions. - Interplay

More information

Principles of Quantum Mechanics

Principles of Quantum Mechanics Principles of Quantum Mechanics - indistinguishability of particles: bosons & fermions bosons: total wavefunction is symmetric upon interchange of particle coordinates (space,spin) fermions: total wavefuncftion

More information

Introduction to Heisenberg model. Javier Junquera

Introduction to Heisenberg model. Javier Junquera Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics

More information

LUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G

LUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G Q1: Display all the MOs for N2 in your report and classify each one of them as bonding, antibonding or non-bonding, and say whether the symmetry of the orbital is σ or π. Sketch a molecular orbital diagram

More information

J 12 J 23 J 34. Driving forces in the nano-magnetism world. Intra-atomic exchange, electron correlation effects: Inter-atomic exchange: MAGNETIC ORDER

J 12 J 23 J 34. Driving forces in the nano-magnetism world. Intra-atomic exchange, electron correlation effects: Inter-atomic exchange: MAGNETIC ORDER Driving forces in the nano-magnetism world Intra-atomic exchange, electron correlation effects: LOCAL (ATOMIC) MAGNETIC MOMENTS m d or f electrons Inter-atomic exchange: MAGNETIC ORDER H exc J S S i j

More information

Exchange Mechanisms. Erik Koch Institute for Advanced Simulation, Forschungszentrum Jülich. lecture notes:

Exchange Mechanisms. Erik Koch Institute for Advanced Simulation, Forschungszentrum Jülich. lecture notes: Exchange Mechanisms Erik Koch Institute for Advanced Simulation, Forschungszentrum Jülich lecture notes: www.cond-mat.de/events/correl Magnetism is Quantum Mechanical QUANTUM MECHANICS THE KEY TO UNDERSTANDING

More information

Metal-insulator transitions

Metal-insulator transitions Metal-insulator transitions What can we learn from electronic structure calculations? Mike Towler mdt26@phy.cam.ac.uk www.tcm.phy.cam.ac.uk/ mdt26 Theory of Condensed Matter Group Cavendish Laboratory

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

QUANTUM CHEMISTRY FOR TRANSITION METALS

QUANTUM CHEMISTRY FOR TRANSITION METALS QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline

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

2 B B D (E) Paramagnetic Susceptibility. m s probability. A) Bound Electrons in Atoms

2 B B D (E) Paramagnetic Susceptibility. m s probability. A) Bound Electrons in Atoms Paramagnetic Susceptibility A) Bound Electrons in Atoms m s probability B +½ p ½e x Curie Law: 1/T s=½ + B ½ p + ½e +x With increasing temperature T the alignment of the magnetic moments in a B field is

More information

Electronic structure calculations results from LDA+U method

Electronic structure calculations results from LDA+U method Electronic structure calculations results from LDA+U method Vladimir I. Anisimov Institute of Metal Physics Ekaterinburg, Russia LDA+U method applications Mott insulators Polarons and stripes in cuprates

More information

CHEM6085: Density Functional Theory Lecture 10

CHEM6085: Density Functional Theory Lecture 10 CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired

More information

Calculation of exchange integrals and electronic structure for manganese ferrite

Calculation of exchange integrals and electronic structure for manganese ferrite Calculation of exchange integrals and electronic structure for manganese ferrite Xu Zuo and Carmine Vittoria Department of Electrical and Computer Engineering, Northeastern University, Boston, Massachusetts

More information

Electronic correlation and Hubbard approaches

Electronic correlation and Hubbard approaches Electronic correlation and Hubbard approaches Matteo Cococcioni Department of Chemical Engineering and Materials Science University of Minnesota Notable failures of LDA/GGA: transition-metal oxides Introduction

More information

University of Bristol. 1 Naval Research Laboratory 2 II. Physikalisches Institut, Universität zu Köln

University of Bristol. 1 Naval Research Laboratory 2 II. Physikalisches Institut, Universität zu Köln Charge ordering as alternative to Jahn-Teller distortion In collaboration with Michelle Johannes 1, Daniel Khomskii 2 (theory) and Mohsen Abd-Elmeguid et al 2, Radu Coldea et al 3 (experiment) 1 Naval

More information

Metal-insulator and magnetic transition of NiO at high pressures

Metal-insulator and magnetic transition of NiO at high pressures Metal-insulator and magnetic transition of NiO at high pressures Xiao-Bing Feng Department of Physics, Dalian Railway Institute, Dalian 116028, Peoples Republic of China N. M. Harrison Department of Chemistry,

More information

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National

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

Lecture B6 Molecular Orbital Theory. Sometimes it's good to be alone.

Lecture B6 Molecular Orbital Theory. Sometimes it's good to be alone. Lecture B6 Molecular Orbital Theory Sometimes it's good to be alone. Covalent Bond Theories 1. VSEPR (valence shell electron pair repulsion model). A set of empirical rules for predicting a molecular geometry

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

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

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

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

The Gutzwiller Density Functional Theory

The Gutzwiller Density Functional Theory The Gutzwiller Density Functional Theory Jörg Bünemann, BTU Cottbus I) Introduction 1. Model for an H 2 -molecule 2. Transition metals and their compounds II) Gutzwiller variational theory 1. Gutzwiller

More information

An ab initio study of Fe-doped Nickel Oxide. Kamaludin Dingle

An ab initio study of Fe-doped Nickel Oxide. Kamaludin Dingle An ab initio study of Fe-doped Nickel Oxide Kamaludin Dingle Summer 2006 Abstract In the following pages we investigate the magnetic properties of the dilute magnetic semiconductor Fe-doped NiO from an

More information

Transition Elements. pranjoto utomo

Transition Elements. pranjoto utomo Transition Elements pranjoto utomo Definition What is transition metal? One of which forms one or more stable ions which have incompletely filled d orbitals. 30Zn? Definition Zink is not transition elements

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

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

2.1 Experimental and theoretical studies

2.1 Experimental and theoretical studies Chapter 2 NiO As stated before, the first-row transition-metal oxides are among the most interesting series of materials, exhibiting wide variations in physical properties related to electronic structure.

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

Magnetic Oxides. Gerald F. Dionne. Department of Materials Science and Engineering Massachusetts Institute of Technology

Magnetic Oxides. Gerald F. Dionne. Department of Materials Science and Engineering Massachusetts Institute of Technology Magnetic Oxides Gerald F. Dionne Department of Materials Science and Engineering Massachusetts Institute of Technology Spins in Solids Summer School University of Virginia Charlottesville, VA 21 June 2006

More information

MOLECULAR MAGNETISM. Leigh Jones Room 133 School of Chemistry NUI Galway. Introduction to Molecular Magnetism

MOLECULAR MAGNETISM. Leigh Jones Room 133 School of Chemistry NUI Galway. Introduction to Molecular Magnetism 4 th year undergraduate course ecture 5 MOECUAR MAGNETISM eigh Jones Room 133 School of Chemistry NUI Galway ecture 6: 5: Outcomes Introduction to Molecular Magnetism To understand the difference between

More information

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 3502/4502. Final Exam Part I. May 14, 2005 Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle

More information

Electron Correlation

Electron Correlation Series in Modern Condensed Matter Physics Vol. 5 Lecture Notes an Electron Correlation and Magnetism Patrik Fazekas Research Institute for Solid State Physics & Optics, Budapest lb World Scientific h Singapore

More information

First principle calculations of plutonium and plutonium compounds: part 1

First principle calculations of plutonium and plutonium compounds: part 1 First principle calculations of plutonium and plutonium compounds: part 1 A. B. Shick Institute of Physics ASCR, Prague, CZ Outline: u Lecture 1: Methods of Correlated band theory DFT and DFT+U u Lecture

More information

Mott insulators. Mott-Hubbard type vs charge-transfer type

Mott insulators. Mott-Hubbard type vs charge-transfer type Mott insulators Mott-Hubbard type vs charge-transfer type Cluster-model description Chemical trend Band theory Self-energy correction Electron-phonon interaction Mott insulators Mott-Hubbard type vs charge-transfer

More information

First-Principles Calculation of Exchange Interactions

First-Principles Calculation of Exchange Interactions Chapter 2 First-Principles Calculation of Exchange Interactions Before introducing the first-principles methods for the calculation of exchange interactions in magnetic systems we will briefly review two

More information

7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig

7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space

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

Electronic structure of correlated electron systems. G.A.Sawatzky UBC Lecture

Electronic structure of correlated electron systems. G.A.Sawatzky UBC Lecture Electronic structure of correlated electron systems G.A.Sawatzky UBC Lecture 6 011 Influence of polarizability on the crystal structure Ionic compounds are often cubic to maximize the Madelung energy i.e.

More information

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012 2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.

More information

Supplementary Information

Supplementary Information Supplementary Information Supplementary Figure S1: Structure and composition of Teflon tape. (a) XRD spectra of original Teflon tape and Teflon tape subjected to annealing at 150 o C under Ar atmosphere.

More information

Electronic structure of correlated electron systems. Lecture 2

Electronic structure of correlated electron systems. Lecture 2 Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No

More information

v(r i r j ) = h(r i )+ 1 N

v(r i r j ) = h(r i )+ 1 N Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j

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

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

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

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

Mott insulators. Introduction Cluster-model description Chemical trend Band description Self-energy correction

Mott insulators. Introduction Cluster-model description Chemical trend Band description Self-energy correction Mott insulators Introduction Cluster-model description Chemical trend Band description Self-energy correction Introduction Mott insulators Lattice models for transition-metal compounds Hubbard model Anderson-lattice

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

Molecules in strong magnetic fields

Molecules in strong magnetic fields Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway

More information

Electronic communication through molecular bridges Supporting Information

Electronic communication through molecular bridges Supporting Information Electronic communication through molecular bridges Supporting Information Carmen Herrmann and Jan Elmisz Institute of Inorganic and Applied Chemistry, University of Hamburg, Martin-Luther-King-Platz 6,

More information

Lecture contents. Magnetic properties Diamagnetism Band paramagnetism Atomic paramagnetism Ferromagnetism. Molecular field theory Exchange interaction

Lecture contents. Magnetic properties Diamagnetism Band paramagnetism Atomic paramagnetism Ferromagnetism. Molecular field theory Exchange interaction 1 Lecture contents Magnetic properties Diamagnetism and paramagnetism Atomic paramagnetism Ferromagnetism Molecular field theory Exchange interaction NNSE 58 EM Lecture #1 [SI] M magnetization or magnetic

More information

Correlation effects in MgO and CaO: Cohesive energies and lattice constants

Correlation effects in MgO and CaO: Cohesive energies and lattice constants PHYSICAL REVIEW B VOLUME 54, NUMBER 19 15 NOVEMBER 1996-I Correlation effects in MgO and CaO: Cohesive energies and lattice constants Klaus Doll and Michael Dolg Max-Planck-Institut für Physik Komplexer

More information

X-Ray Magnetic Dichroism. S. Turchini ISM-CNR

X-Ray Magnetic Dichroism. S. Turchini ISM-CNR X-Ray Magnetic Dichroism S. Turchini SM-CNR stefano.turchini@ism.cnr.it stefano.turchini@elettra.trieste.it Magnetism spin magnetic moment direct exchange: ferro antiferro superexchange 3d Ligand 2p 3d

More information

Hybrid density functional theory applied to magnetite: Crystal structure, charge order, and phonons

Hybrid density functional theory applied to magnetite: Crystal structure, charge order, and phonons Hybrid density functional theory applied to magnetite: Crystal structure, charge order, and phonons Andrew D. Rowan and Charles H. Patterson School of Physics, Trinity College Dublin, Dublin 2, Ireland

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

CHAPTER 2 MAGNETISM. 2.1 Magnetic materials

CHAPTER 2 MAGNETISM. 2.1 Magnetic materials CHAPTER 2 MAGNETISM Magnetism plays a crucial role in the development of memories for mass storage, and in sensors to name a few. Spintronics is an integration of the magnetic material with semiconductor

More information

Competing Ferroic Orders The magnetoelectric effect

Competing Ferroic Orders The magnetoelectric effect Competing Ferroic Orders The magnetoelectric effect Cornell University I would found an institution where any person can find instruction in any study. Ezra Cornell, 1868 Craig J. Fennie School of Applied

More information

Quantum chemical modelling of molecular properties - parameters of EPR spectra

Quantum chemical modelling of molecular properties - parameters of EPR spectra Quantum chemical modelling of molecular properties - parameters of EPR spectra EPR ( electric paramagnetic resonance) spectra can be obtained only for open-shell systems, since they rely on transitions

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

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

IFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations

IFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Introduction: The objective of the second computer laboratory exercise is to get acquainted with a program for performing quantum chemical

More information

X-ray absorption spectroscopy.

X-ray absorption spectroscopy. X-ray absorption spectroscopy www.anorg.chem.uu.nl/people/staff/frankdegroot/ X-ray absorption spectroscopy www.anorg.chem.uu.nl/people/staff/frankdegroot/ Frank de Groot PhD: solid state chemistry U Nijmegen

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

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

Interaction of matter with magnetic fields

Interaction of matter with magnetic fields LN07-1 Interaction of matter with magnetic fields All substances have magnetic properties, which can be determined by examining their behaviour in the presence of an external magnetic field, H. N S When

More information

Magnetism at finite temperature: molecular field, phase transitions

Magnetism at finite temperature: molecular field, phase transitions Magnetism at finite temperature: molecular field, phase transitions -The Heisenberg model in molecular field approximation: ferro, antiferromagnetism. Ordering temperature; thermodynamics - Mean field

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

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

CHAPTER 4. ELECTRONIC AND MAGNETIC PROPERTIES OF MX 2 (M = V, Nb; X = Al, Ga, In, Cl, Br AND I) COMPOUNDS IN CdI 2 -TYPE STRUCTURE

CHAPTER 4. ELECTRONIC AND MAGNETIC PROPERTIES OF MX 2 (M = V, Nb; X = Al, Ga, In, Cl, Br AND I) COMPOUNDS IN CdI 2 -TYPE STRUCTURE 84 CHAPTER 4 ELECTRONIC AND MAGNETIC PROPERTIES OF MX 2 (M = V, Nb; X = Al, Ga, In, Cl, Br AND I) COMPOUNDS IN CdI 2 -TYPE STRUCTURE 4.1 INTRODUCTION As ideal materials for use in spintronic devices, the

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

Supporting Information for Ultra-narrow metallic armchair graphene nanoribbons

Supporting Information for Ultra-narrow metallic armchair graphene nanoribbons Supporting Information for Ultra-narrow metallic armchair graphene nanoribbons Supplementary Figure 1 Ribbon length statistics. Distribution of the ribbon lengths and the fraction of kinked ribbons for

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

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 3502/4502. Final Exam Part I. May 14, 2005 Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The

More information

College of Chemistry, Peking University, Beijing, China. Fritz-Haber-Institut der MPG, Berlin, Germany

College of Chemistry, Peking University, Beijing, China. Fritz-Haber-Institut der MPG, Berlin, Germany KITP Program Excitations in Condensed Matter Localized and Itinerant States in a Unified Picture beyond Density Functional Theory Hong Jiang 1, Patrick Rinke 2 and Matthias Scheffler 2 1 College of Chemistry,

More information

Density Functional Theory

Density Functional Theory Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function

More information

Time-Dependent Density-Functional Theory

Time-Dependent Density-Functional Theory Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique

More information

Crystalline and Magnetic Anisotropy of the 3d Transition-Metal Monoxides

Crystalline and Magnetic Anisotropy of the 3d Transition-Metal Monoxides Crystalline and of the 3d Transition-Metal Monoxides Institut für Festkörpertheorie und -optik Friedrich-Schiller-Universität Max-Wien-Platz 1 07743 Jena 2012-03-23 Introduction Crystalline Anisotropy

More information

Winter School for Quantum Magnetism EPFL and MPI Stuttgart Magnetism in Strongly Correlated Systems Vladimir Hinkov

Winter School for Quantum Magnetism EPFL and MPI Stuttgart Magnetism in Strongly Correlated Systems Vladimir Hinkov Winter School for Quantum Magnetism EPFL and MPI Stuttgart Magnetism in Strongly Correlated Systems Vladimir Hinkov 1. Introduction Excitations and broken symmetry 2. Spin waves in the Heisenberg model

More information

MODULE 2: QUANTUM MECHANICS. Principles and Theory

MODULE 2: QUANTUM MECHANICS. Principles and Theory MODULE 2: QUANTUM MECHANICS Principles and Theory You are here http://www.lbl.gov/cs/html/exascale4energy/nuclear.html 2 Short Review of Quantum Mechanics Why do we need quantum mechanics? Bonding and

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

Crystal structure prediction

Crystal structure prediction Crystal structure prediction Klaus Doll Institute for Mathematical Physics, TU Braunschweig Max Planck Institute for Solid State Research, Stuttgart MSSC2011, Turin, September 2011 Motivation structure

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

Computational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica

Computational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Ab Initio Implementations of Hartree-Fock Molecular Orbital Theory Fundamental assumption of HF theory: each electron

More information

c E If photon Mass particle 8-1

c E If photon Mass particle 8-1 Nuclear Force, Structure and Models Readings: Nuclear and Radiochemistry: Chapter 10 (Nuclear Models) Modern Nuclear Chemistry: Chapter 5 (Nuclear Forces) and Chapter 6 (Nuclear Structure) Characterization

More information

Lecture 19: Building Atoms and Molecules

Lecture 19: Building Atoms and Molecules Lecture 19: Building Atoms and Molecules +e r n = 3 n = 2 n = 1 +e +e r y even Lecture 19, p 1 Today Nuclear Magnetic Resonance Using RF photons to drive transitions between nuclear spin orientations in

More information

1.1 Atoms. 1.1 Atoms

1.1 Atoms. 1.1 Atoms 1. Chemical bonding and crystal structure 19 21 Hydrogen atom Scanning electron microscopy Ni surface Cleaved surface ZnO, TiO 2, NiO, NaCl, Si, Ge, GaAs, InP Crystals are build by small repeating units

More information

CHAPTER 8. EXPLORING THE EFFECTS OF STRUCTURAL INSTABILITIES AND OF TRIPLET INSTABILITY ON THE M-I PHASE TRANSITION IN (EDO-TTF)2PF6.

CHAPTER 8. EXPLORING THE EFFECTS OF STRUCTURAL INSTABILITIES AND OF TRIPLET INSTABILITY ON THE M-I PHASE TRANSITION IN (EDO-TTF)2PF6. CHAPTER 8. EXPLORING THE EFFECTS OF STRUCTURAL INSTABILITIES AND OF TRIPLET INSTABILITY ON THE M-I PHASE TRANSITION IN (EDO-TTF)2PF6. The thermal metal-insulator phase transition in (EDO-TTF)2PF6 is attributed

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

Disturbing the dimers: electron- and hole-doping in the nonmagnetic intermetallic insulator FeGa 3

Disturbing the dimers: electron- and hole-doping in the nonmagnetic intermetallic insulator FeGa 3 Disturbing the dimers: electron- and hole-doping in the nonmagnetic intermetallic insulator FeGa 3 Antia S. Botana, Yundi Quan, Warren E. Pickett University of California-Davis, Davis, California, USA

More information

Quantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno

Quantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion

More information

Magnetism (FM, AFM, FSM)

Magnetism (FM, AFM, FSM) Magnetism (FM, AFM, FSM) Karlheinz Schwarz Institute of Materials Chemistry TU Wien Localized vs. itinerant systems In localized systems (e.g. some rare earth) the magnetism is mainly governed by the atom

More information

RFSS: Lecture 8 Nuclear Force, Structure and Models Part 1 Readings: Nuclear Force Nuclear and Radiochemistry:

RFSS: Lecture 8 Nuclear Force, Structure and Models Part 1 Readings: Nuclear Force Nuclear and Radiochemistry: RFSS: Lecture 8 Nuclear Force, Structure and Models Part 1 Readings: Nuclear and Radiochemistry: Chapter 10 (Nuclear Models) Modern Nuclear Chemistry: Chapter 5 (Nuclear Forces) and Chapter 6 (Nuclear

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) Silvana Botti

Key concepts in Density Functional Theory (II) Silvana Botti Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:

More information

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

The electronic structure of materials 1

The electronic structure of materials 1 Quantum mechanics 2 - Lecture 9 December 18, 2013 1 An overview 2 Literature Contents 1 An overview 2 Literature Electronic ground state Ground state cohesive energy equilibrium crystal structure phase

More information