Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin...

32
Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre et Marie Curie”, LCT-CNRS, Paris

Transcript of Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin...

Page 1: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Vallico, july 2013

Correlated maximum probability domains in crystalsMauro Causa’ Andreas SavinUniversita’ di Napoli “Federico II”Universitè “Pierre et Marie Curie”, LCT-CNRS, Paris

Page 2: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

C

O

Chemical “observable” Observable ρ (r)

In so far quantum mechanics is correct, chemical questions are

problems in “applied” mathematics (Eyring, Walter, Kimball, 1944)

Page 3: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Chemists like a local topology

Typical questions:

Atomic charges ?

Oxidation numbers ?

Bonds (order)?

Is more covalent or ionic?

Page 4: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

ELF

Electron

Localization

Function

(Becke)

AIM

Atoms in

Molecule

(Bader)

MPD

Maximum

Probability

Domains

(Savin)

Page 5: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Atoms in Molecules R.Bader, 1960->

ρ(r) , number electronic density

Ω

Ω

N=drrρ

Bader, R. F. W. “Atoms in molecules: a quantum theory” Oxford

University Press: Oxford,1990

Page 6: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Analysis of the ρ(r) scalar field

Page 7: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 8: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Electron Localization Function (ELF)

1. (Local (Ek) cost for respecting Pauli principle)normalized

Page 9: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Elf region definition

Page 10: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Topological analysis and bifurcation

Page 11: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

TiO2: ionic(0.6) and covalent(0.3) elf

Page 12: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Probability functions

p2(bond)=0.3-0.5

Pn(core)=0.7-1.0

Page 13: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 14: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Probability and number density

Ω

Ω

N=drrρ

Page 15: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Maximum Probability Domains (MPD)unique definition of a domain

Page 16: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

molecules

Page 17: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

MPD for solidsCrystalline solids only shows extreme and

peculiar chemical bonds

(ionic, metallic, covalent)

Based ol Localized Orbitals (for insulators)

For delocalized system->QMC

Page 18: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 19: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 20: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 21: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 22: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Ionic crystals

MC, AS JPC 115 2011 13139

Page 23: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Anions

Page 24: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 25: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 26: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.
Page 27: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Covalent crystals

Page 28: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Voronoi ELF MPD ELF-MPD

MC, AS ZAAC 637 2011 882

Page 29: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Semi-ionic crystals: α-SiO2

P(cov)=0.38;

p(O-)=0.37,p(O2-)=0.58

Page 30: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

MPD-QMC(CASINO)

25 crystals in test

•Covalent

•Ionic

•Semi-ionic

•Metals

•Rare gases

•Easy work on configuration file for

VMC and DMC

•Configurations are distributed like

•|Ψvmc|2 |Ψdmc Ψvmc |

•Try to reducing numerical noise of

Shape Derivative

Page 31: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

MPD-QMC preliminary results (VMC)

•For covalent and ionic crystal the shape of the domain

• doesn’t change (accuracy 0.1%)

•In covalent systems the core and bond probability increases:

•P2(bond,Si): 0.42 ---> 0.44

•P10(core,Si): 0.69 ---> 0.73

•In ionic crystals the Ionic probability seems unaffected by

electronic correlation

Page 32: Vallico, july 2013 Correlated maximum probability domains in crystals Mauro Causa’ Andreas Savin Universita’ di Napoli “Federico II” Universitè “Pierre.

Next steps

Complete the data collection on all the 25 test

systems

Metals !!!!!! (QMC only)

DMC

Better accuracy (number of configurations)

Multiple domain probabilities (QMC only)