Measuring correlation functions in interacting systems of cold atoms

Anatoli Polkovnikov Ehud Altman Vladimir Gritsev Mikhail Lukin Eugene Demler Harvard/Boston University Harvard/Weizmann Harvard Harvard Harvard Thanks to: M. Greiner , Z. Hadzibabic, M. Oberthaler, J. Schmiedmayer, V. Vuletic

Correlation functions in condensed matter physics Most experiments in condensed matter physics measure correlation functions Example: neutron scattering measures spin and density correlation functions Neutron diffraction patterns for MnO Shull et al., Phys. Rev. 83:333 (1951)

Outline

Lecture I : Measuring correlation functions in intereference experiments Lecture II : Quantum noise interferometry in time of flight experiments Emphasis of these lectures: detection and characterization of many-body quantum states

Lecture I

Measuring correlation functions in intereference experiments 1. Interference of independent condensates 2. Interference of interacting 1D systems 3. Interference of 2D systems 4. Full distribution function of the fringe amplitudes in intereference experiments. 5. Studying coherent dynamics of strongly interacting systems in interference experiments

Lecture II

Quantum noise interferometry in time of flight experiments 1. Detection of spin order in Mott states of atomic mixtures 2. Detection of fermion pairing

Measuring correlation functions in intereference experiments

Analysis of high order correlation functions in low dimensional systems

Polkovnikov, Altman, Demler, PNAS (2006)

Interference of two independent condensates Andrews et al., Science 275:637 (1997)

r r’ Interference of two independent condensates 1 2 d r+d Clouds 1 and 2 do not have a well defined phase difference.

However each individual measurement shows an interference pattern

x 1 Interference of one dimensional condensates d Experiments: Schmiedmayer et al., Nature Physics (2005) Amplitude of interference fringes , , contains information about phase fluctuations within individual condensates x 2 x y

Interference amplitude and correlations

L For identical condensates Instantaneous correlation function

Interacting bosons in 1d at T=0 Low energy excitations and long distance correlation functions can be described by the Luttinger Hamiltonian.

K

– Luttinger parameter Connection to original bosonic particles Small K corresponds to strong quantum fluctuations

Luttinger liquids in 1d For non-interacting bosons For impenetrable bosons Correlation function decays rapidly for small K. This decay comes from strong quantum fluctuations

L Interference between 1d interacting bosons Luttinger liquid at T=0

K

– Luttinger parameter For non-interacting bosons and For impenetrable bosons and Luttinger liquid at finite temperature Analysis of can be used for thermometry

Rotated probe beam experiment

q For large imaging angle, , Luttinger parameter

K

may be extracted from the angular dependence of

Interference between two-dimensional BECs at finite temperature.

Kosteritz-Thouless transition

Interference of two dimensional condensates Experiments: Stock, Hadzibabic, Dalibard, et al., cond-mat/0506559 Gati, Oberthaler, et al., cond-mat/0601392 L x L y L x Probe beam parallel to the plane of the condensates

Interference of two dimensional condensates.

Quasi long range order and the KT transition L y L x Above Kosterlitz-Thouless transition: Vortices proliferate. Short range order Below Kosterlitz-Thouless transition: Vortices confined. Quasi long range order Above KT transition Below KT transition

Experiments with 2D Bose gas

Hadzibabic et al., Nature (2006) Time of flight

z

Typical interference patterns low temperature higher temperature

x

z D x

Experiments with 2D Bose gas

Hadzibabic et al., Nature (2006)

x

integration Contrast after integration over

x

axis

z

0.4

low

T

integration 0.2

middle

T

over

x

axis

z

integration over

x

axis

z

high

T

0 0 10 20 30 integration distance

D x

(pixels)

Experiments with 2D Bose gas

0.4

0.2

0 0 Hadzibabic et al., Nature (2006) middle

T

low

T

fit by:

C

2 ~ 1

D x



D x



g

1 ( 0 ,

x

)  2

dx

~   1

D x

  2 a Exponent a high

T

10 20 30 0.5

integration distance

D x

0.4

if

g

1 (

r

) with decays exponentially : if

g

1 (

r

) decays algebraically or exponentially with a large : 0.3

0 high

T

0.1

0.2

central contrast 0.3

low

T

“Sudden” jump!?

Experiments with 2D Bose gas

Hadzibabic et al., Nature (2006) Exponent a

c.f. Bishop and Reppy

1.0

0 1.0

1.1

1.2

T

(K) He experiments: universal jump in the superfluid density 0.5

0.4

0.3

high

T

low

T

0 0.1

0.2

central contrast 0.3

Ultracold atoms experiments: jump in the correlation function.

KT theory predicts a =1/4 just below the transition

Experiments with 2D Bose gas. Proliferation of thermal vortices Haddzibabic et al., Nature (2006) 30% Fraction of images showing at least one dislocation 20% 10% 0 0 high

T

0.1

0.2

0.3

central contrast low

T

0.4

Rapidly rotating two dimensional condensates Time of flight experiments with rotating condensates correspond to density measurements Interference experiments measure single particle correlation functions in the rotating frame

Interference between two interacting one dimensional Bose liquids

Full distribution function of the amplitude of interference fringes

Gritsev, Altman, Demler, Polkovnikov, cond-mat/0602475

L Higher moments of interference amplitude is a quantum operator. The measured value of will fluctuate from shot to shot.

Can we predict the distribution function of ?

Higher moments Changing to periodic boundary conditions (long condensates) Explicit expressions for are available but cumbersome Fendley, Lesage, Saleur, J. Stat. Phys. 79:799 (1995)

Impurity in a Luttinger liquid Expansion of the partition function in powers of

g

Partition function of the impurity contains correlation functions taken at the same point and at different times. Moments of interference experiments come from correlations functions taken at the same time but in different points. Euclidean invariance ensures that the two are the same

Relation between quantum impurity problem and interference of fluctuating condensates Normalized amplitude of interference fringes Distribution function of fringe amplitudes Relation to the impurity partition function Distribution function can be reconstructed from using completeness relations for the Bessel functions

Bethe ansatz solution for a quantum impurity can be obtained from the Bethe ansatz following Zamolodchikov, Phys. Lett. B 253:391 (91); Fendley, et al., J. Stat. Phys. 79:799 (95) Making analytic continuation is possible but cumbersome Interference amplitude and spectral determinant is related to the single particle Schroedinger equation Dorey, Tateo, J.Phys. A. Math. Gen. 32:L419 (1999) Bazhanov, Lukyanov, Zamolodchikov, J. Stat. Phys. 102:567 (2001) Spectral determinant

K=1 K=1.5

K=3 K=5 Narrow distribution for .

Approaches Gumble distribution. Width 0 1 x 2 3 4 Wide Poissonian distribution for

From interference amplitudes to conformal field theories correspond to vacuum eigenvalues of

Q

operators of CFT Bazhanov, Lukyanov, Zamolodchikov, Comm. Math. Phys.1996, 1997, 1999 When K>1, is related to Q operators of CFT with c<0. This includes 2D quantum gravity, non intersecting loop model on 2D lattice, growth of random fractal stochastic interface, high energy limit of multicolor QCD, … 2D quantum gravity, non-intersecting loops on 2D lattice Yang-Lee singularity

Studying coherent dynamics of strongly interacting systems in interference experiments

J

Coupled 1d systems

Motivated by experiments of Schmiedmayer et al.

Interactions lead to phase fluctuations within individual condensates Tunneling favors aligning of the two phases Interference experiments measure only the relative phase

J

Coupled 1d systems

Conjugate variables Relative phase Particle number imbalance Small K corresponds to strong quantum flcutuations

Hamiltonian Quantum Sine-Gordon model Imaginary time action Quantum Sine-Gordon model is exactly integrable Excitations of the quantum Sine-Gordon model soliton antisoliton breather

Coherent dynamics of quantum Sine-Gordon model Motivated by experiments of Schmiedmayer et al.

J

Prepare a system at t=0 Take to the regime of finite tunneling and let evolve for some time Measure amplitude of interference pattern

Coherent dynamics of quantum Sine-Gordon model Amplitude of interference fringes Oscillations or decay?

time

From integrability to coherent dynamics At t=0 we have a state with for all This state can be written as a “squeezed” state Matrix can be constructed using connection to boundary SG model Calabrese, Cardy (2006); Ghoshal, Zamolodchikov (1994) Time evolution can be easily written Interference amplitude can be calculated using form factor approach Smirnov (1992), Lukyanov (1997)

Coherent dynamics of quantum Sine-Gordon model

J

Prepare a system at t=0 Take to the regime of finite tunneling and let evolve for some time Measure amplitude of interference pattern

Coherent dynamics of quantum Sine-Gordon model Amplitude of interference fringes time Amplitude of interference fringes shows oscillations at frequencies that correspond to energies of breater

Conclusions for part I Interference of fluctuating condensates can be used to probe correlation functions in one and two dimensional systems. Interference experiments can also be used to study coherent dynamics of interacting systems

Measuring correlation functions in interacting systems of cold atoms

Lecture II

Quantum noise in time of flight interferometry experiments 1. Time of flight experiments.

Second order coherence in Mott states of spinless bosons 2. Detection of spin order in Mott states of atomic mixtures 3. Detection of fermion pairing Emphasis of these lectures: detection and characterization of many-body quantum states

Bose-Einstein condensation

Cornell et al., Science 269, 198 (1995) Ultralow density condensed matter system Interactions are weak and can be described theoretically from first principles

Superfluid to Insulator transition

Greiner et al., Nature 415:39 (2002) 

U

Mott insulator

n

 1 Superfluid t/U

Time of flight experiments Quantum noise interferometry of atoms in an optical lattice Second order coherence

Second order coherence in the insulating state of bosons.

Hanburry-Brown-Twiss experiment Theory: Altman et al., PRA 70:13603 (2004) Experiment: Folling et al., Nature 434:481 (2005)

Hanburry-Brown-Twiss stellar interferometer

Hanburry-Brown-Twiss interferometer

Second order coherence in the insulating state of bosons Bosons at quasimomentum expand as plane waves with wavevectors First order coherence : Oscillations in density disappear after summing over Second order coherence : Correlation function acquires oscillations at reciprocal lattice vectors

Second order coherence in the insulating state of bosons.

Hanburry-Brown-Twiss experiment Theory: Altman et al., PRA 70:13603 (2004) Experiment: Folling et al., Nature 434:481 (2005)

Effect of parabolic potential on the second order coherence Experiment: Spielman, Porto, et al., Theory: Scarola, Das Sarma, Demler, PRA (2006) Width of the correlation peak changes across the transition, reflecting the evolution of Mott domains

Width of the noise peaks

Interference of an array of independent condensates Hadzibabic et al., PRL 93:180403 (2004) 2.5

3 Smooth structure is a result of finite experimental resolution (filtering) 1.4

1.2

2 1 1.5

0.8

1 0.6

0.5

0.4

0 0.2

-0.5

0 -1 -1.5

0 200 400 600 800 1000 1200 -0.2

0 200 400 600 800 1000 1200

Applications of quantum noise interferometry in time of flight experiments Detection of spin order in Mott states of boson boson mixtures

Engineering magnetic systems using cold atoms in an optical lattice

See also lectures by A. Georges and I. Cirac in this school

Spin interactions using controlled collisions Experiment: Mandel et al., Nature 425:937(2003) Theory: Jaksch et al., PRL 82:1975 (1999)

Two component Bose mixture in optical lattice Example: . Mandel et al., Nature 425:937 (2003) t t Two component Bose Hubbard model

Quantum magnetism of bosons in optical lattices Kuklov and Svistunov, PRL (2003) Duan et al., PRL (2003) • Ferromagnetic • Antiferromagnetic

Exchange Interactions in Solids

antibonding bonding Kinetic energy dominates: antiferromagnetic state Coulomb energy dominates: ferromagnetic state

Hysteresis Two component Bose mixture in optical lattice.

Mean field theory + Quantum fluctuations Altman et al., NJP 5:113 (2003) 1 st order

Probing spin order of bosons Correlation Function Measurements

 Engineering exotic phases • Optical lattice in 2 or 3 dimensions: polarizations & frequencies of standing waves can be different for different directions YY ZZ • Example: exactly solvable model Kitaev (2002), honeycomb lattice with

H



J x



i

, 

j



x



i x



x j



J y



i

, 

j



y



i y



y j



J z



i

, 

j



z



i z



z j

• • Can be created with 3 sets of standing wave light beams !

Non-trivial topological order, “spin liquid” + non-abelian anyons …those has not been seen in controlled experiments

Applications of quantum noise interferometry in time of flight experiments Detection of fermion pairing

Fermionic atoms in optical lattices Pairing in systems with repulsive interactions. Unconventional pairing. High Tc mechanism

Fermionic atoms in a three dimensional optical lattice

Kohl et al ., PRL 94:80403 (2005)

See also lectures of T. Esslinger and W. Ketterle in this school

Fermions with repulsive interactions t t

Possible d-wave pairing of fermions

U

High temperature superconductors

Picture courtesy of UBC Superconductivity group Superconducting Tc 93 K Hubbard model – minimal model for cuprate superconductors P.W. Anderson, cond-mat/0201429 After many years of work we still do not understand the fermionic Hubbard model

Positive U Hubbard model

Possible phase diagram. Scalapino, Phys. Rep. 250:329 (1995) Antiferromagnetic insulator D-wave superconductor

Second order correlations in the BCS superfluid

n(k) n(r’) k F BCS BEC k n(r)

Expansion of atoms in TOF maps

k

into

r



n

(

r

,

r

' ) 

n

(

r

) 

n

(

r

' ) 

n

(

r

, 

r

) 

BCS

 0

Momentum correlations in paired fermions

Greiner et al., PRL 94:110401 (2005)

Fermion pairing in an optical lattice

Second Order Interference In the TOF images Normal State Superfluid State measures the Cooper pair wavefunction One can identify unconventional pairing

Simulation of condensed matter systems: Hubbard Model and high Tc superconductivity U t t Personal opinion : The fermionic Hubbard model contains 90% of the physics of cuprates. The remaining 10% may be crucial for getting high Tc superconductivity. Understanding Hubbard model means finding what these missing 10% are. Electron-phonon interaction?

Mesoscopic structures (stripes)?

Using cold atoms to go beyond “plain vanilla” Hubbard model a) Boson-Fermion mixtures: Hubbard model + phonons b) Inhomogeneous systems, role of disorder

Boson Fermion mixtures

Fermions interacting with phonons

Boson Fermion mixtures

See lectures by T. Esslinger and G. Modugno in this school BEC

Bosons provide cooling for fermions and mediate interactions. They create non-local attraction between fermions

Charge Density Wave Phase Periodic arrangement of atoms Non-local Fermion Pairing P-wave, D wave, …

Boson Fermion mixtures

“Phonons” : Bogoliubov (phase) mode Effective fermion ”phonon” interaction Fermion ”phonon” vertex Similar to electron-phonon systems

Boson Fermion mixtures in 1d optical lattices

Cazalila et al., PRL (2003); Mathey et al., PRL (2004) Spinless fermions Spin ½ fermions

Boson Fermion mixtures in 2d optical lattices

Wang et al., PRA (2005) 40K -- 87Rb 40K -- 23Na (a) (b) =1060nm =765.5nm

=1060 nm

Conclusions Interference of extended condensates is a powerful tool for analyzing correlation functions in one and two dimensional systems Noise interferometry can be used to probe quantum many-body states in optical lattices