Two-step melting of three-sublattice order Kedar Damle, TIFR, Mumbai - - PowerPoint PPT Presentation

two step melting of three sublattice order
SMART_READER_LITE
LIVE PREVIEW

Two-step melting of three-sublattice order Kedar Damle, TIFR, Mumbai - - PowerPoint PPT Presentation

Two-step melting of three-sublattice order Kedar Damle, TIFR, Mumbai JNU-ICTP Workshop Feb 10 2015 Symmetry breaking transitions: Generalities Symmetry-breaking state characterized by long-range correlations of order-parameter O


slide-1
SLIDE 1

Two-step melting of three-sublattice order

Kedar Damle, TIFR, Mumbai JNU-ICTP Workshop Feb 10 2015

slide-2
SLIDE 2

Symmetry breaking transitions: Generalities

◮ Symmetry-breaking state characterized by long-range

correlations of “order-parameter” ˆ O

◮ phenomenological Landau free energy density F[ˆ

O] Expanding F in powers of ˆ O (symmetry allowed terms)

◮ Neglecting derivatives (fluctuations):

phase transition → change in minimum of F

slide-3
SLIDE 3

Fluctuation effects at continuous transitions:

◮ More complete description of long-wavelength physics:

Include (symmetry allowed) gradient terms in F

◮ In most cases: Corrections to mean-field exponents ◮ In rare cases: Fluctuation-induced first-order behaviour

slide-4
SLIDE 4

Symmetries are (usually) decisive:

◮ Transformation properties of ˆ

O determine nature of continuous transition

slide-5
SLIDE 5

In this talk...

◮ Two well-known scenarios for continuous melting of

three-sublattice order in frustrated triangular and Kagome-lattice easy-axis antiferromagnets: Two-step melting with intermediate power-law ordered phase with power-law exponent η(T) ∈ ( 1

9, 1 4)

OR Three-state Potts transition

◮ Main message of talk—

Thermodynamic signature of two-step melting process distinguishes between the two kinds of continuous transitions

slide-6
SLIDE 6

Frustrated easy-axis antiferromagnets

◮ Easy-axis n and triangular motifs...

? +n −n

slide-7
SLIDE 7

Wannier’s triangular lattice Ising antiferromagnet

◮ HIsing = J

ij σz i σz j on the triangular lattice

◮ T → 0 limit characterized by power-law correlations:

σz

rσz 0 ∼ cos(Q·r) r1/2

Incipient order at three-sublattice wavevector Q = (2π/3, 2π/3) Stephenson (1964) Power-law spin-liquid in the T → 0 limit

slide-8
SLIDE 8

Lattice-gas models for monolayers on graphite

◮ Three-sublattice long-range order of noble-gas monolayers on

graphite HJ1J2 = J

ij σz i σz j − J1

  • ij σz

i σz j − J2 · · · − B i σz i

Long-range three-sublattice ordering (wavevector Q) at low temperature

  • D. P

. Landau (1983)

slide-9
SLIDE 9

Prototypical example of order-by-quantum fluctuations

◮ HTFIM = J

ij σz i σz j − Γ i σx i on the triangular lattice

Long-range order at three-sublattice wavevector Q

◮ Equivalent: Plaquette-ordered valence-bond-solid state of

honeycomb lattice quantum dimer model Moessner, Sondhi, Chandra (2001), Isakov & Moessner (2003)

slide-10
SLIDE 10

Ferri vs antiferro three-sublattice order

x

e

y

e R a b c b a b a c c c c c c

π

6

θ +S +S −S −S +S

ψ = |ψ|eiθ = −

  • R eiQ·

RSz

  • R

Ferri vs antiferro order distinguished by the choice of phase θ Ferri: θ = 2πm/6, Antiferro: θ = (2m + 1)π/6 (m = 0, 1, 2 . . . 5)

slide-11
SLIDE 11

S = 1 antiferromagnets with single-ion anisotropy

◮ HAF = J

ij

Si · Sj − D

i(Sz i)2 on triangular lattice

◮ Low-energy physics for D ≫ J:

Hb = − J2

D

  • ij(b†

i bj + h.c.) + J ij(ni − 1 2)(nj − 1 2) + . . . ◮ Low-temperature state for D ≫ J: “supersolid” state of hard-core

bosons at half-filling on triangular lattice with unfrustrated hopping t = J2/D and frustrating nearest-neighbour repulsion U = J

◮ Implies: Coexisting three-sublattice order in Sz and

“ferro-nematic” order in S2

⊥ (KD & Senthil 2006)

(Simple easy-axis version of Chandra-Coleman (1991) “spin-nematic” ideas)

slide-12
SLIDE 12

Is three-sublattice ordering of Sz in HAF ferri or antiferro?

◮ Natural expectation: Quantum fluctuations induce antiferro order

→ Ordering will be antiferro three-sublattice order (like transverse field Ising antiferromagnet)

  • e. g. Melko et. al. (2005)
slide-13
SLIDE 13

QMC evidence: Ferri three-sublattice order of Sz

  • 0.005

0.005

δρ=ρ−1/2

0.1

(δρ) β=20 (a) U=10, L=48 , P

Heidarian and KD (2005)

slide-14
SLIDE 14

Ising models for “Artificial Kagome-ice”

◮ HKagome = J

ij σz i σz j − J1

  • ij σz

i σz j − J2 . . .

◮ Only nearest-neighbour couplings → classical short-range spin

liquid (Kano & Naya 1950)

◮ Further neighbour couplings destabilize spin liquid →

three-sublattice order at low T (Wolff & Schotte 1988)

◮ “Artificial Kagome-ice: Moments Mi = σz

i ni

(ni at different sites non-collinear) Expt: Tanaka et. al. (2006), Qi et. al. (2008), Ladak et. al. (2010,11) Theory: Moller, Moessner (2009), Chern, Mellado, Tchernyshyov (2011)

slide-15
SLIDE 15

Three-sublattice order on the Kagome lattice

π

6

θ +S +S −S −S +S

y

e

x

e R R

x

e +

y

e R+ 2 2 a b c 2 1 1 1 2 1 b c 2 c b

ψ = |ψ|eiθ = −

  • R
  • α=0,1,2 eiQ·

R−2πi α

3 Sz

  • R,α

Ferri vs antiferro distinguished by the choice of phase θ Ferri: θ = 2πm/6, Antiferro: θ = (2m + 1)π/6 (m = 0, 1, 2 . . . 5)

slide-16
SLIDE 16

Landau-theory for melting of three-sublattice order

◮ F = K|∇ψ|2 + r|ψ|2 + u|ψ|4 + λ6(ψ6 + ψ∗6) + . . .

Connection to physics of six-state clock models Z =

{pi} exp[ ij V( 2π 6 (pi − pj))]

Each pi = 0, 1, 2, ...5 V(x) = K1 cos(x) + K2 cos(2x) + K3 cos(3x) Cardy (1980)

slide-17
SLIDE 17

Melting scenarios for three-sublattice order

◮ Analysis (Cardy 1980) of generalized six-state clock models

→ Three generic possibilities of relevance here: Two-step melting, with power-law ordered intermediate phase OR 3-state Potts transition OR First-order transition (always possible!) Both these continuous melting scenarios realized in one or more examples on triangular and kagome lattices

slide-18
SLIDE 18

Nature of melting transition in triangular lattice supersolid?

◮ Clearly: Nature of melting transition not a priori obvious ◮ Prediction of Boninsegni & Prokofiev (2005)

Three-state Potts transition Prediction based on argument about relative energies of different kinds of domain walls hard to get right at quantitative level

slide-19
SLIDE 19

Our answer from large-scale QMC simulations

KD & Heidarian (in preparation)

slide-20
SLIDE 20

Detecting power-law order?

Need extremely sensitive scattering experiment to detect power-law version of Bragg peaks Or High resolution real-space data by scanning some local probe + Lots

  • f image-processing

difficult!

slide-21
SLIDE 21

Alternate thermodynamic signature(!)

◮ Singular thermodynamic susceptibility to uniform easy-axis field

B: χu(B) ∼

1 |B|p(T)

◮ p(T) = 4−18η(T)

4−9η(T) for η(T) ∈ ( 1 9, 2 9)

So p(T) varies from 1/3 to 0 as T increases from T1 to just below T2 (KD 2014, with referees)

slide-22
SLIDE 22

Review: picture for power-law ordered phase

◮ In state with long-range three-sublattice order, θ feels λ6 cos(6θ)

potential. Locks into values 2πm/6 (resp. (2m + 1)π/6) in ferri (resp. antiferro) three-sublattice ordered state for T < T1

◮ In power-law three-sublattice ordered state for T ∈ (T1, T2), λ6

does not pin phase θ θ spread uniformly (0, 2π) Distinction between ferri and antiferro three-sublattice order lost for T ∈ (T1, T2)

slide-23
SLIDE 23

Review: more formal RG description

◮ Fixed point free-energy density: FKT

kBT = 1 4πg(∇θ)2

with g(T) ∈ ( 1

9, 1 4) corresponding to T ∈ (T1, T2)

◮ λ6 cos(6θ) irrelevant along fixed line ◮ ψ∗(r)ψ(0) ∼

1 rη(T)

with η(T) = g(T) Jose, Kadanoff, Kirkpatrick, Nelson (1977)

slide-24
SLIDE 24

General argument for result—I

◮ Landau theory admits term λ3m(ψ3 + ψ∗3)

m is uniform magnetization mode

◮ Formally irrelevant along fixed line FKT

→ Physics of two-step melting unaffected—m “goes for a ride...” But ...

slide-25
SLIDE 25

General argument for result—II

◮ m “inherits” power-law correlations of cos(3θ):

Cm(r) = m(r)m(0) ∼

1 r9η(T)

◮ χL ∼

L d2rCm(r) in a finite-size system at B = 0

◮ χL ∼ L2−9η(T) for η(T) ∈ ( 1

9, 2 9)

Diverges with system size at B = 0

slide-26
SLIDE 26

General argument for result—III

◮ Uniform field B > 0 → additional term h3 cos(3θ) in FKT ◮ Strongly relevant along fixed line, with RG eigenvalue 2 − 9g/2 ◮ Implies finite correlation length ξ(B) ∼ |B|−

2 4−9η

◮ χu(B) ∼ |B|− 4−18η

4−9η for η(T) ∈ ( 1

9, 2 9)

slide-27
SLIDE 27

The proof of the pudding...I

0.75 1.5 3 0.5 1 2 4 50 100 200 400 Cψ L

3

  • L

10−5 10−3 10−1 50 100 200 400 Cσ L 3

  • L

R = −1; κ = −1 T = 4.3 η = 0.14 T = 4.5 η = 0.17 T = 4.7 η = 0.21 T = 4.8 η = 0.24

In power-law ordered phase of HJ1J2 (R = −(J1 + J2)/J and κ = (J2 − J1)/J) (Ghanshyam, KD (in preparation))

slide-28
SLIDE 28

The proof of the pudding...II

1 10

χ0

L (Linear dimension of system)

data a L2-9η 100 1000 10000 36 48 60 72 84

χQ

L (Linear dimension of system)

η=0.142(2)

Data a L2-η

0.1 1 10

χ0

L (Linear dimension of system)

data a L2-9η 100 1000 10000 36 48 60 72 84

χQ

L (Linear dimension of system)

η=0.171(2)

Data a L2-η

In power-law ordered phase of HTFIM (Biswas, KD (in preparation))

slide-29
SLIDE 29

The proof of the pudding...III

In power-law ordered phase of Hb (KD, Heidarian (in preparation))

slide-30
SLIDE 30

Acknowledgements

◮ Collaborators:

QMC on Hb: Dariush Heidarian (Toronto) QMC on HTFIM: Geet Ghanshyam and Sounak Biswas (TIFR) Classical MC of HJ1J2: Geet Ghanshyam (TIFR)

◮ Computational resources of Dept. of Theoretical Physics, TIFR