Spin-crossover molecules: puzzling systems for electronic structure methods. Andrea Droghetti School of Physics and CRANN, Trinity College Dublin

Similar documents
First- principles studies of spin-crossover molecules

Molecular spintronics with spin crossover compounds

Transport properties of topological insulators. Andrea Droghetti School of Physics and CRANN, Trinity College Dublin, IRELAND

Organometallic optoelectronically active magnetic molecules for logic and memory

Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization

Variational Monte Carlo Optimization and Excited States

QMC dissociation energy of the water dimer: Time step errors and backflow calculations

Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign

Optimization of quantum Monte Carlo wave functions by energy minimization

Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables

Modified Becke-Johnson (mbj) exchange potential

Theoretical study of electronic and atomic structures of (MnO)n

QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters

Exchange-Correlation Functional

Quantum Monte Carlo backflow calculations of benzene dimers

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects

5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry.

arxiv: v1 [cond-mat.mtrl-sci] 21 Dec 2007

Quantum Monte Carlo tutorial. Lucas K. Wagner Dept. of Physics; University of Illinois at Urbana-Champaign

Introduction to Density Functional Theory

Quantum Monte Carlo methods

Electron transport through Shiba states induced by magnetic adsorbates on a superconductor

The Overhauser Instability

Multi-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005

QUANTUM CHEMISTRY FOR TRANSITION METALS

SUPPORTING INFORMATION. Table S1: Use of different functionals and variation of HF exchange on IS/HS splitting

Chapter 20 d-metal complexes: electronic structures and properties

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

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

MAGNETIC ANISOTROPY IN TIGHT-BINDING

Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn

Molecular Orbital Theory (MOT)

Metal-insulator transitions

Density Functional Theory - II part

Wave Function Optimization and VMC

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

One-Electron Hamiltonians

Quantum Mechanical Simulations

Fixed-Node quantum Monte Carlo for Chemistry

0 belonging to the unperturbed Hamiltonian H 0 are known

Fireball.

Electronic Structure Calculations and Density Functional Theory

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

Introduction to Density Functional Theory

Size-extensive wave functions for QMC A linear-scaling GVB approach

B. Electron Deficient (less than an octet) H-Be-H. Be does not need an octet Total of 4 valence electrons

Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139

Valence electronic structure of isopropyl iodide investigated by electron momentum spectroscopy. --- Influence of intramolecular interactions

Many-body correlations in a Cu-phthalocyanine STM single molecule junction

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

Module 6 1. Density functional theory

Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems

Harju, A.; Sverdlov, V.A.; Nieminen, Risto; Halonen, V. Many-body wave function for a quantum dot in a weak magnetic field

Electronic correlations in models and materials. Jan Kuneš

Cr(II) or Cr 2+ Consider the octahedral complex Cr[(en) 3 ] 2+ Octahedral complex with 4 d electrons. Octahedral complex with 4 d electrons

Is the homogeneous electron gas homogeneous?

Uncertainty in Molecular Photoionization!

CASSCF and NEVPT2 calculations: Ground and excited states of multireference systems. A case study of Ni(CO)4 and the magnetic system NArO

Algorithmic Challenges in Photodynamics Simulations

Binding energy of bilayer graphene and Electronic properties of oligoynes

Pseudopotentials for hybrid density functionals and SCAN

H!!!! = E! Lecture 7 - Atomic Structure. Chem 103, Section F0F Unit II - Quantum Theory and Atomic Structure Lecture 7. Lecture 7 - Introduction

Binding energy of 2D materials using Quantum Monte Carlo

Phonon Transport in Nanostructures

Principles of Quantum Mechanics

Molecular Magnetic Properties. The 11th Sostrup Summer School. Quantum Chemistry and Molecular Properties July 4 16, 2010

Dept of Mechanical Engineering MIT Nanoengineering group

Designing interfaces for Spin Injection into Organic Molecular Solids: A Surface Science Approach

Bayesian Error Estimation in Density Functional Theory

Density Functional Theory

Algorithmic Challenges in Photodynamics Simulations on HPC systems

Orbital dependent correlation potentials in ab initio density functional theory

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

The Self Interaction Correction revisited

1 Density functional theory (DFT)

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Exploring Ligand Field Effects on Fe(IV)-oxo Reactivity Channels for Versatility as a Potential Catalyst for Strong C-H Bond Activation

Density Functional Theory (DFT)

Bonding in Coordination Compounds. Crystal Field Theory. Bonding in Transition Metal Complexes

Noncollinear spins in QMC: spiral Spin Density Waves in the HEG

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

Electronic structure / bonding in d-block complexes

Quantum Monte Carlo wave functions and their optimization for quantum chemistry

Manganese Chemistry. Mn : 1s 2 2s 2 2p 6 3s 2 3p 6 4s 2 3d 5

Quantum Chemical and Dynamical Tools for Solving Photochemical Problems

2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density

Bonding in transition metal complexes

Dept of Mechanical Engineering MIT Nanoengineering group

The Basics of Theoretical and Computational Chemistry

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

Density Functional Theory. Martin Lüders Daresbury Laboratory

The Schrödinger equation for many-electron systems

All-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe

X-ray absorption spectroscopy.

Density functional theory in the solid state

Supplementary Information for "Atomistic. Simulations of Highly Conductive Molecular. Transport Junctions Under Realistic Conditions"

Towards quantifying the role of exact exchange in. predictions of transition metal complex properties

Lecture 15 February 15, 2013 Transition metals

Spin crossover phenomena in transition metal oxides under high magnetic field

Atomic basis sets for first-principle studies of Si nanowires. Electronics Theory group Tyndall National Institute, University College Cork

Transcription:

Spin-crossover molecules: puzzling systems for electronic structure methods Andrea Droghetti School of Physics and CRANN, Trinity College Dublin

Introduction

Introduction Can we read a molecule spin? Can we access electronic (and vibranic) states? Can we write and read information in the spin state of a molecule?

Introduction Electronic structure of the molecule Electronic structure of the device (equlibrium properties) Transport properties of the device (out-of-equilibrium properties)

Introduction Electronic structure of the molecule DFT DMC CASSCF Electronic structure of the device (equlibrium properties) DFT beyond-dft Transport properties of the device (out-of-equilibrium properties) DFT Beyond-DFT (?) (TDDFT, MDPT?)

Introduction Electronic structure of the molecule DFT DMC CASSCF 1- Accurate results (?) Electronic structure of the device (equlibrium properties) DFT beyond-dft Transport properties of the device (out-of-equilibrium properties) DFT Beyond-DFT (?) (TDDFT, MDPT?)

Introduction Electronic structure of the molecule DFT DMC CASSCF Electronic structure of the device (equlibrium properties) DFT beyond-dft Transport properties of the device (out-of-equilibrium properties) DFT Beyond-DFT (?) (TDDFT, MDPT?)

DFT for infinite systems HM H LM H MR. HL HR

DFT for infinite systems HM ΣL+ +ΣR. www.smeagol.tcd.ie

DFT for infinite systems HM +ΣR ΣL+. www.smeagol.tcd.ie

The Theoretical Method HM +ΣR ΣL+.

The Theoretical Method HM +ΣR ΣL+.

The Theoretical Method H ' M =H 0 +V ee ΣL+ +ΣR.

Smeagol Unitary trasf. for removing orbital overlap Continuous-time quantum Monte Carlo Double-counting correction (interaction/hybridization expansion) Dyson eq. if scf.true. Else (analytic continuation) A. Droghetti & I. Rungger Spectral function

The Theoretical Method H ' M =H 0 +V ee. DMC two-particle density matrix = model two-particle density matrix (Lucas Wagner)

Electronic structure of the molecule DFT DMC CASSCF 1- Accurate results (?) - Allow to extract first-principles model Hamiltonians Electronic structure of the device (equlibrium properties) DFT beyond-dft Transport properties of the device (out-of-equilibrium properties) DFT Beyond-DFT (?) (TDDFT, MDPT?)

Spin-Crossover molecules

Spin-crossover molecules eg dx -y dz tg dxy dxz dyz

Spin-crossover molecules eg dx -y dz tg dxy dxz dyz

Spin-crossover molecules eg dx -y dz tg dxy dxz dyz

Spin-crossover molecules dx -y eg eg* dz tg dxy dxz dyz tg Fe eg eg Ligands

Spin-crossover molecules dx -y dz eg eg* tg dxy dxz dyz tg Fe eg eg Ligands

Magnetic molecules G=G HS G LS = E T S E 0 T C= S 0 E S

Transport through SCMs F. Prinz et al., Adv. Mat., 3, 1545 (011) N. Baadji & S. Sanvito, Phys. Rev. Lett. 108, 1701 (01) M.S. Alam et al., Angew. Chem. Int. Ed., 49, 1159 (010)

Model systems [Fe(HO)6]+ [Fe(NH3)6]+ [Fe(NCH)6]+ lower crystal field [Fe(CO)6]+ higher crystal field Spectrochemical series

Model systems [Fe(HO)6]+ [Fe(NH3)6]+ [Fe(NCH)6]+ lower crystal field [Fe(CO)6]+ higher crystal field Expected ground state: S= Δ E adia <0 S= Δ E adia <0 S= Δ E adia <0 S=0 Δ E adia >0

[Fe(HO)6]+ [Fe(NH3)6]+ ΔEadia (ev) LDA -0.49 GGA -1.04 B3LYP -1.37 PBE0-1.74 HH -.6 DMC -.54(1) ΔEadia (ev) LDA 0.96 GGA 0.08 B3LYP -0.59 PBE0-0.88 HH -1.68 DMC -1.59(1) E adia 0 E HS E LS [Fe(NCH)6]+ ΔEadia (ev) [Fe(CO)6]+ ΔEadia (ev) LDA 5.18 LDA.37 GGA 3.4 GGA 1.04 B3LYP 1.7 B3LYP -0.0 PBE0 1.3 PBE0-0.44 HH 0.4 HH -1.49 DMC -1.37(3) DMC 0.37(3) Caculations with CASINO Dirac-Fock Trail-Needs PP (small core for Fe) Slater-Jastrow (LDA orbitals) A. Droghetti & S. Sanvito and in collaboration with D. Alfe' (UCL)

Spin Crossover [FeL](BF4) (L=,6-dypirazol-1-yl-4-hydroxymethylpyridine) V.A. Money et al., Dalton Trans. 10, 1516 (004) adia ΔE ΔEadia (ev) (ev) GGA 1.4 LDA -0.49 B3LYP GGA 0.01-1.04 PBE0 B3LYP -0.3-1.36 HH -1.33 PB0-1.74 DMC -1.1(4) DFT A. Droghetti, D. Alfe' and S. Sanvito, J. Chem. Phys. 137, 14303 (01) A. Droghetti, D. Alfe' and S. Sanvito, Phys. Rev. B 87, 05114 (013)

Spin Crossover adia ΔE ΔEadia (ev) (ev) GGA 0.44 LDA -0.49 B3LYP GGA -0.33-1.04 PBE0 B3LYP -0.54-1.36 DMC -0.69(8) DFT CASPT 0.10 DMC with QWALK, (new) Burkatzki-Filippi-Dong PP Slater-Jastrow Calculations performed by A. Droghetti and M. Fumanal (Barcelona) Thanks to Lucas Wagner for the code and technical support

Spin Crossover [Fe(NCH)6]+ CAS(10,1) 4d-Fe dx -y eg dz eg* tg dxy dxz dyz tg Fe eg eg Ligands

Spin Crossover [Fe(NCH)6]+ dx -y tg eg dz eg* tg dxy eg CAS(10,1) 4d-Fe eg* dxz dyz tg Fe eg

Spin Crossover [Fe(NCH)6]+ dx -y tg eg dz eg* HF tg dxy eg CAS(10,1) 4d-Fe eg* dxz dyz tg Fe eg ΔEadia -3.8 ev CASSCF -1.96 ev CASPT -0.33 ev

Spin Crossover [Fe(NCH)6]+ ΔEadia HF -3.8 ev B3LYP -0.9 ev PBE 0.76 ev CASSCF -1.96 ev CASPT -0.33 ev ΔEadia DMC (HF) -1.16(5) ev DMC -1.05(5) ev (B3LYP) DMC -0.87(4) ev (PBE) DMC (Slater-Jastrow)

Spin Crossover [Fe(NCH)6]+ Energy (a.u) CSF1-0.698(1) CSF7-0.700(1) CSF13-0.706(1) CSF1-0.700(1)

Spin Crossover Energy (a.u) CSF1-0.645(1) CSF7-0.648(1) CSF13-0.647(1) CSF3-0.649(1)

Conclusions DFT, DMC and CASPT return different results (different order of magnitude!!) CASPT chemsists say that we should trust it as reference DMC...trial wave-function seems not to matter much Pseudo-potential...hard to say...but can a bad pseudo-potential be responsible for an systematic error of about more than 0.5 ev on an energy difference? What else?

The algorithm Acknowledgements - European Union for funding (ACMOL projects) - Irish Centre For High-End Computing (ICHEC) and Trinity Centre for High Performing Computing (TCHPC) for computational resources - Alin Elena (ICHEC) for technical support Thank you!