Berry Phases and Curvatures in Electronic-Structure Theory David - - PowerPoint PPT Presentation

berry phases and curvatures in electronic structure theory
SMART_READER_LITE
LIVE PREVIEW

Berry Phases and Curvatures in Electronic-Structure Theory David - - PowerPoint PPT Presentation

Berry Phases and Curvatures in Electronic-Structure Theory David Vanderbilt Rutgers University March APS Meeting, Baltimore, March 13 2006 Rahman prize for: Theory of polarization (King-Smith & Vanderbilt) Ultrasoft


slide-1
SLIDE 1

March APS Meeting, Baltimore, March 13 2006

Berry Phases and Curvatures in Electronic-Structure Theory

David Vanderbilt Rutgers University

slide-2
SLIDE 2

March APS Meeting, Baltimore, March 13 2006

Rahman prize for: – Theory of polarization (King-Smith & Vanderbilt) – Ultrasoft pseudopotentials Three quick preliminaries:

  • Who was Aneesur Rahman?
  • Who is Dominic King-Smith?
  • A parable about referee reports…
slide-3
SLIDE 3

March APS Meeting, Baltimore, March 13 2006

Who was Aneesur Rahman?

Photo courtesy Sam Bader via Marie-Louise Saboungi

  • Born Hyderbad, India
  • Educ. Cambridge, Louvain
  • Argonne Natl. Labs 1960-85
  • U. Minnesota 1985-87
  • Died 1987
  • Rahman Prize established in

1992 with funds from IBM

“Father of Molecular Dynamics”

slide-4
SLIDE 4

March APS Meeting, Baltimore, March 13 2006

Who is Dominic King-Smith?

“Father of Bettina”

Accelrys Job title: “Product Manager, Quantum Mechanics”

  • PhD, Cambridge, UK
  • Postdoc at Rutgers `91-`93
  • Biosym/MSI/Accelrys `93-`01
  • Presently at:
slide-5
SLIDE 5

March APS Meeting, Baltimore, March 13 2006

Ultrasoft Pseudopotentials

slide-6
SLIDE 6

March APS Meeting, Baltimore, March 13 2006

Berry Phases and Curvatures in Electronic-Structure Theory

David Vanderbilt Rutgers University

slide-7
SLIDE 7

March APS Meeting, Baltimore, March 13 2006

Introduction

  • By mid-1990s, density-functional perturbation theory

allowed calculation of linear response to E-field

  • However, it was not known how to:

– Calculate polarization itself – Treat finite E-fields

  • Analogous problem of calculating orbital magnetization

also unsolved

slide-8
SLIDE 8

March APS Meeting, Baltimore, March 13 2006

Introduction

  • Solutions of these problems are now in hand

– Modern theory of polarization (1993) – Treatment of finite E-fields (2002) – Orbital magnetization (2005)

  • Solutions rely heavily on two crucial ingredients:

– Wannier functions – Berry phases and related quantities This talk: Brief survey of methods! Almost nothing on applications

slide-9
SLIDE 9

March APS Meeting, Baltimore, March 13 2006

Outline of Talk

  • Introduction
  • Berry phases, potentials, and curvatures
  • Realizations:

– Electric polarization – Wannier functions – Electric fields – Anomalous Hall conductivity – Orbital magnetization

  • Summary and prospects
slide-10
SLIDE 10

March APS Meeting, Baltimore, March 13 2006

Berry phases

u3Ò u2Ò unÒ =u1Ò u4Ò … un-1Ò

Now take limit that density of points Æ∞

slide-11
SLIDE 11

March APS Meeting, Baltimore, March 13 2006

Berry phases

ulÒ

l=0 l=1 Continuum limit

slide-12
SLIDE 12

March APS Meeting, Baltimore, March 13 2006

(Context: Molecular coordinates)

ulÒ

l=0 l=1

(z1, z2)

z1 z2

Na3

slide-13
SLIDE 13

March APS Meeting, Baltimore, March 13 2006

Context: k-space in Brillouin zone

ukÒ

l=0 l=1 kx ky

2p/a

Bloch function

slide-14
SLIDE 14

March APS Meeting, Baltimore, March 13 2006

Stokes theorem: Berry curvature

ukÒ

kx ky

2p/a

W

slide-15
SLIDE 15

March APS Meeting, Baltimore, March 13 2006

Context: k-space in Brillouin zone

ukÒ

l=0 l=1 kx ky

2p/a

Bloch function

slide-16
SLIDE 16

March APS Meeting, Baltimore, March 13 2006

Spanning the BZ

Bloch function

ukÒ

l=0 l=1 kx ky

2p/a

slide-17
SLIDE 17

March APS Meeting, Baltimore, March 13 2006

Does any of this have any connection to real physics

  • f materials?
slide-18
SLIDE 18

March APS Meeting, Baltimore, March 13 2006

Outline of Talk

  • Introduction
  • Berry phases, potentials, and curvatures
  • Realizations:

– Electric polarization – Wannier functions – Electric fields – Anomalous Hall conductivity – Orbital magnetization

  • Summary and prospects
slide-19
SLIDE 19

March APS Meeting, Baltimore, March 13 2006

P = dcell / Vcell ?

+ – + – + – + – + – + –

  • Textbook picture

(Claussius-Mossotti)

  • But does not correspond

to reality!

slide-20
SLIDE 20

Ferroelectric PbTiO3 (Courtesy N. Marzari)

slide-21
SLIDE 21

March APS Meeting, Baltimore, March 13 2006

P = dcell / Vcell ? dcell =

slide-22
SLIDE 22

March APS Meeting, Baltimore, March 13 2006

P = dcell / Vcell ? dcell =

slide-23
SLIDE 23

March APS Meeting, Baltimore, March 13 2006

Berry-phase theory of electric polarization

slide-24
SLIDE 24

March APS Meeting, Baltimore, March 13 2006

Berry-phase theory of electric polarization

Berry potential!

slide-25
SLIDE 25

March APS Meeting, Baltimore, March 13 2006

Simplify: 1 band, 1D

ukÒ

l=0 l=1 kx ky

2p/a

slide-26
SLIDE 26

March APS Meeting, Baltimore, March 13 2006

Discrete sampling of k-space

slide-27
SLIDE 27

March APS Meeting, Baltimore, March 13 2006

Discretized formula in 3D

where

slide-28
SLIDE 28

March APS Meeting, Baltimore, March 13 2006

Sample Application: Born Z*

Paraelectric Ferroelectric

+2 e ? +4 e ? – 2 e ? – 2 e ?

slide-29
SLIDE 29

March APS Meeting, Baltimore, March 13 2006

Outline of Talk

  • Introduction
  • Berry phases, potentials, and curvatures
  • Realizations:

– Electric polarization – Wannier functions – Electric fields – Anomalous Hall conductivity – Orbital magnetization

  • Summary and prospects
slide-30
SLIDE 30

March APS Meeting, Baltimore, March 13 2006

Wannier function representation

(Marzari and Vanderbilt, 1997)

“Wannier center”

slide-31
SLIDE 31

March APS Meeting, Baltimore, March 13 2006

Mapping to Wannier centers

Wannier center rn

slide-32
SLIDE 32

March APS Meeting, Baltimore, March 13 2006

Wannier dipole theorem DP = Sion (Zione) Drion + Swf (– 2e) Drwf

  • Exact!
  • Gives local description of

dielectric response!

Mapping to Wannier centers

slide-33
SLIDE 33

Ferroelectric BaTiO3 (Courtesy N. Marzari)

slide-34
SLIDE 34

Wannier functions in a-Si

Fornari et al.

Wannier functions in l-H2O

Silvestrelli et al.

slide-35
SLIDE 35

Courtesy S. Nakhmanson

Wannier analysis of PVDF polymers and copolymers

slide-36
SLIDE 36

March APS Meeting, Baltimore, March 13 2006

Note upcoming release of public max-loc Wannier code… (Organized by Nicola Marzari)

slide-37
SLIDE 37

March APS Meeting, Baltimore, March 13 2006

Outline of Talk

  • Introduction
  • Berry phases, potentials, and curvatures
  • Realizations:

– Electric polarization – Wannier functions – Electric fields – Anomalous Hall conductivity – Orbital magnetization

  • Summary and prospects
slide-38
SLIDE 38

March APS Meeting, Baltimore, March 13 2006

Electric Fields: The Problem

Easy to do in practice: For small E-fields, tZener >> tUniverse ; is it OK? But ill-defined in principle:

Zener tunneling

slide-39
SLIDE 39

March APS Meeting, Baltimore, March 13 2006

Electric Fields: The Problem

  • is not periodic
  • Bloch’s theorem does not apply

y(x) is very

messy

slide-40
SLIDE 40

March APS Meeting, Baltimore, March 13 2006

Electric Fields: The Solution

  • Seek long-lived resonance
  • Described by Bloch functions
  • Minimizing the electric enthalpy functional

(Nunes and Gonze, 2001)

Usual EKS Berry phase polarization

Souza, Iniguez, and Vanderbilt, PRL 89, 117602 (2002);

  • P. Umari and A. Pasquarello, PRL 89, 157602 (2002).
slide-41
SLIDE 41

March APS Meeting, Baltimore, March 13 2006

Electric Fields: Implementation

As long as Dk is not too small:

  • Can use standard methods to find minimum
  • The e · P term introduces coupling between k-points

p/a

–p/a k

slide-42
SLIDE 42

March APS Meeting, Baltimore, March 13 2006

Sample Application: Born Z*

Can check that previous results for BaTiO3 are reproduced

slide-43
SLIDE 43

March APS Meeting, Baltimore, March 13 2006

Outline of Talk

  • Introduction
  • Berry phases, potentials, and curvatures
  • Realizations:

– Electric polarization – Wannier functions – Electric fields – Anomalous Hall conductivity – Orbital magnetization

  • Summary and prospects
slide-44
SLIDE 44

March APS Meeting, Baltimore, March 13 2006

Anomalous Hall effect

Ferromagnetic Material

slide-45
SLIDE 45

March APS Meeting, Baltimore, March 13 2006

Anomalous Hall effect

  • Karplus-Luttinger theory (1954)

– Scattering-free, intrinsic

  • Skew-scattering mechanism (1955)

– Impurity scattering

  • Side-jump mechanism (1970)

– Impurity or phonon scattering

  • Berry-phase theory (1999)

– Restatement of Karplus-Luttinger

Semiclassical equations

  • f motion:

Sundaram and Niu, PRB 59, 14925 (1999).

slide-46
SLIDE 46

March APS Meeting, Baltimore, March 13 2006

Stokes theorem: Berry curvature

ukÒ

kx ky

2p/a

W

slide-47
SLIDE 47

March APS Meeting, Baltimore, March 13 2006

  • Z. Fang et al, Science 302,

92 (2003).

Wz for kz=0

Anomalous Hall conductivity of SrRuO3

slide-48
SLIDE 48

March APS Meeting, Baltimore, March 13 2006

  • X. Wang, J. Yates, I. Souza, and D. Vanderbilt, G23.00001

(Tuesday 8am).

See also Y.G. Yao et al., PRL 92, 037204 (2004).

Wz(kx,kz) in bcc Fe

slide-49
SLIDE 49

March APS Meeting, Baltimore, March 13 2006

Outline of Talk

  • Introduction
  • Berry phases, potentials, and curvatures
  • Realizations:

– Electric polarization – Wannier functions – Electric fields – Anomalous Hall conductivity – Orbital magnetization

  • Summary and prospects
slide-50
SLIDE 50

March APS Meeting, Baltimore, March 13 2006

Orbital Magnetization

M is a bulk property?

fl fl K = M x n

K is only apparently a surface property?

K +s

  • s

P is a bulk property

fl fl s = P ⋅ n

s is only apparently a surface property

slide-51
SLIDE 51

March APS Meeting, Baltimore, March 13 2006

Theory of orbital magnetization

  • T. Thonhauser, D. Ceresoli, D. Vanderbilt, and R. Resta,
  • Phys. Rev. Lett. 95, 137205 (2005).
  • Context:

– Ferromagnetic insulators – Single-particle approximation – Vanishing magnetic field

  • Used Wannier representation to derive a

formula for the orbital magnetization

slide-52
SLIDE 52

March APS Meeting, Baltimore, March 13 2006

Orbital currents in Wannier representation r r

= +

ÔwsÒ ÔwsÒ ÔwsÒ

Local Circulation (LC) Itinerant Circulation (IC) ·vÒ

slide-53
SLIDE 53

March APS Meeting, Baltimore, March 13 2006

Berry curvature Something new

  • T. Thonhauser, H6.00001 (Tuesday 11:15am)

(invited talk)

See also D. Xiao, J. Shi and Q. Niu, PRL 95, 137204 (2005).

slide-54
SLIDE 54

March APS Meeting, Baltimore, March 13 2006

Summary and Prospects

  • Berry phases are everywhere!
  • We discussed:

– Electric polarization – Electric fields – Anomalous Hall coefficient – Orbital magnetization

  • Other “hot topics”:

– Multiferroics and magnetoelectric effects – Single graphene sheets – Spin Hall effect and spin injection

  • More Berry phases lurking around the corner?
slide-55
SLIDE 55

March APS Meeting, Baltimore, March 13 2006

Extras

slide-56
SLIDE 56

March APS Meeting, Baltimore, March 13 2006

Electric Fields: Justification

Seek long-lived metastable periodic solution

slide-57
SLIDE 57

March APS Meeting, Baltimore, March 13 2006

Electric Fields: The Hitch

  • There is a hitch!
  • For given E-field, there is a limit on k-point sampling
  • Length scale LC = 1/Dk
  • Meaning: LC = supercell dimension

Nk = 8 Lc = 8a

  • Solution: Keep Dk > 1/Lt = e/Eg
slide-58
SLIDE 58

March APS Meeting, Baltimore, March 13 2006

  • X. Wang, J. Yates, I. Souza, and D. Vanderbilt, G23.00001

(Tuesday 8am).

Anomalous Hall conductivity of bcc Fe

See also Y.G. Yao et al., PRL 92, 037204 (2004).

slide-59
SLIDE 59

March APS Meeting, Baltimore, March 13 2006

Orbital Magnetization

K = M x n Is M a bulk property? Is K only apparently a surface property? Definition: If K is predetermined at all surfaces in such a way that K = M x n for some vector M, then M is the bulk magnetization.

K

slide-60
SLIDE 60

March APS Meeting, Baltimore, March 13 2006

Orbital Magnetization

Clarification:

  • Microscopic M(r) defined via — x M(r) = J(r)
  • M(r) ill-defined: M(r) fi M(r) + M0 + —h
  • Therefore, cannot define M as cell average of M(r)

Just as: P is not, even in principle, a functional

  • f the bulk charge density distribution r(r)

Conclusion: M is not, even in principle, a functional

  • f the bulk current distribution J(r)

(Hirst, RMP, 1997)

slide-61
SLIDE 61

March APS Meeting, Baltimore, March 13 2006

Strong reasons to expect bulk M

  • Nearsightedness:

Surface current depends only

  • n local environment
  • Stationary quantum state:

dr/dt = 0

  • Conservation of charge:

—⋅J = 0

Edge of type A Edge of type B

Iy

(A)

Iy

(B)

So: Iy

(A) = Iy (B) = Mz

Mz

slide-62
SLIDE 62

March APS Meeting, Baltimore, March 13 2006

Comparison: P vs. M

Defined for insulators only r(r) insufficient in principle; need access to Berry physics r operator “Quantum of polarization” Derivable from adiabatic theory Derivable from Wannier rep. Insulators and metals with broken TR symmetry J(r) insufficient in principle; need access to Berry physics r ¥ v operator No quantum (no monopoles) No obvious adiabatic theory Derivable from Wannier rep.?

Electric Polarization Orbital Magnetization

slide-63
SLIDE 63

March APS Meeting, Baltimore, March 13 2006

Ultrasoft Pseudopotentials

Then, the good news:

Sidney Redner, Physics Today, June 2005.

(A hot paper is…) “defined as a nonreview paper with 350 or more citations, an average ratio of citation age to publication age greater than two-thirds, and a citation rate increasing with time.”

* *

slide-64
SLIDE 64

March APS Meeting, Baltimore, March 13 2006

Ultrasoft Pseudopotentials

Then, the good news:

Sidney Redner, APS talk, March, 2004; Physics Today, June 2005.