Scaling and full counting statistics of interference between independent fluctuating condensates

1. Scaling and full counting statistics of interference between independent fluctuating condensates Anatoli Polkovnikov, Boston University Collaboration: Ehud Altman - Weizmann Eugene Demler - Harvard Vladimir Gritsev - Harvard

2. TOF d Free expansion: x Andrews et. al. 1997 Interference between two condensates.

3. TOF • Correlated phases ( = 0) •  b) Uncorrelated, but well defined phases  int(x)=0 x Hanbury Brown-Twiss Effect Work with original bosonic fields: What do we observe? c) Initial number state. No phases?

4. Interference amplitude squared. Observable! Easy to check that at large N: The interference amplitude does not fluctuate! First theoretical explanation: I. Casten and J. Dalibard (1997): showed that the measurement induces random phases in a thought experiment. Experimental observation of interference between ~ 30 condensates in a strong 1D optical lattice: Hadzibabic et.al. (2004).

5. Z. Hadzibabic et. al., Phys. Rev. Lett. 93, 180401 (2004). Polar plots of the fringe amplitudes and phases for 200 images obtained for the interference of about 30 condensates. (a) Phase-uncorrelated condensates. (b) Phase correlated condensates. Insets: Axial density profiles averaged over the 200 images.

6. Imaging beam  L What if the condensates are fluctuating? This talk: • Access to correlation functions. • Scaling of  AQ2 with L and : power-law exponents. Luttinger liquid physics in 1D, Kosterlitz-Thouless phase transition in 2D. • Probability distribution W(AQ2): all order correlation functions. • Direct simulator (solver) for interacting problems. Quantum impurity in a 1D system of interacting fermions (an example). • Potential applications to many other systems.

7. z1 z2 z AQ Identical homogeneous condensates: x Interference amplitude contains information about fluctuations within each condensate.

8. Dephased condensates: Ideal condensates: L L x x z z Interference contrast does not depend on L. Contrast scales as L-1/2. Scaling with L: two limiting cases

9. L Ideal condensate: Thermal gas: L Formal derivation:

10. L Repulsive bosons with short range interactions: Finite temperature: Intermediate case (quasi long-range order). 1D condensates (Luttinger liquids): z

11. (for the imaging beam orthogonal to the page,  is the angle of the integration axis with respect to z.) x(z2) x(z1) z   q is equivalent to the relative momentum of the two condensates (always present e.g. if there are dipolar oscillations). Angular Dependence.

12. has a cusp singularity for K<1, relevant for fermions. Angular (momentum) Dependence.

13. z CCD camera z Time of x flight y imaging laser x Elongated condensates: Lx>>Ly . Two-dimensional condensates at finite temperature (picture by Z. Hadzibabic)

14. p 0 The phase distribution of an elongated 2D Bose gas. (courtesy of Zoran Hadzibabic) Matter wave interferometry very low temperature: straight fringes which reveal a uniform phase in each plane from time to time: dislocation which reveals the presence of a free vortex higher temperature: bended fringes S. Stock, Z. Hadzibabic, B. Battelier, M. Cheneau, and J. Dalibard: Phys. Rev. Lett. 95, 190403 (2005) “atom lasers”

15. Observing the Kosterlitz-Thouless transition Ly Lx Below KT transition Universal jump of  at TKT Always algebraic scaling, easy to detect. Above KT transition LxLy

16. Interference contrast: x Contrast after integration integration over x axis z 0.4 low T z middle T 0.2 integration over x axis high T z 0 integration distance X (pixels) 0 X 10 20 30 Zoran Hadzibabic, Peter Kruger, Marc Cheneau, Baptiste Battelier, Sabine Stock,and Jean Dalibard (2006).

17. Exponent a high T low T Vortex proliferation central contrast Fraction of images showing at least one dislocation: 30% 0.5 1.0 20% 0.4 T (K) 10% 0 1.1 1.0 1.2 low T high T 0.3 0 0.2 0.3 0.4 0 0.1 central contrast 0 0.1 0.2 0.3 Z. Hadzibabic et. al. “universal jump in the superfluid density” c.f. Bishop and Reppy

18. Identical condensates. Mean: Similarly higher moments Probe of the higher order correlation functions. Distribution function (= full counting statistics): Non-interacting non-condensed regime (Wick’s theorem): Nontrivial statistics if the Wick’s theorem is not fulfilled! Higher Moments. is an observable quantum operator

19. 1D condensates at zero temperature: Low energy action: Then Similarly Easy to generalize to all orders.

20. Changing open boundary conditions to periodic find These integrals can be evaluated using Jack polynomials (Fendley, Lesage, Saleur, J. Stat. Phys. 79:799 (1995)) Explicit expressions are cumbersome (slowly converging series of products).

21. z1 z2 z AQ x Two simple limits: Strongly interacting Tonks-Girardeau regime (also in thermal case) Weakly interacting BEC like regime.

22. Connection to the impurity in a Luttinger liquid problem. Boundary Sine-Gordon theory: P. Fendley, F. Lesage, H. Saleur (1995). Same integrals as in the expressions for (we rely on Euclidean invariance).

23. Distribution of interference phases (and amplitudes) from two 1D condensates. T. Schumm, et. al., Nature Phys. 1, 57 (2005). • Evaluate the integral. Experimental simulation of the quantum impurity problem • Do a series of experiments and determine the distribution function. • Read the result.

24. Use a different approach based on spectral determinant: Dorey, Tateo, J.Phys. A. Math. Gen. 32:L419 (1999); Bazhanov, Lukyanov, Zamolodchikov, J. Stat. Phys. 102:567 (2001) can be found using Bethe ansatz methods for half integer K. In principle we can find W: Difficulties: need to do analytic continuation. The problem becomes increasingly harder as K increases.

25. Evolution of the distribution function.

26. Universal asymmetric distribution at large K (-1)/

27. Further extensions: is the Baxter Q-operator, related to the transfer matrix of conformal field theories with negative charge: 2D quantum gravity, non-intersecting loops on 2D lattice Yang-Lee singularity

28. Note that K+K-1  2, so and the distribution function is always Poissonian. There is a similar cusp at 2kf Spinless Fermions. However for K+K-1  3 there is a universal cusp at nonzero momentum as well as at 2kf: Higher dimensions: nesting of Fermi surfaces, CDW, … Not a low energy probe!

29. 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 Fermions in optical lattices. Possible efficient probes of superconductivity (in particular, d-wave vs. s-wave). Not yet, but coming!

30. Conclusions. • Analysis of interference between independent condensates reveals a wealth of information about their internal structure. • Scaling of interference amplitudes with L or  :correlation function exponents. Working example: detecting KT phase transition. • Probability distribution of amplitudes (= full counting statistics of atoms): information about higher order correlation functions. • Interference of two Luttinger liquids: partition function of 1D quantum impurity problem (also related to variety of other problems like 2D quantum gravity). • Vast potential applications to many other systems, e.g.: • Fermionic systems: superconductivity, CDW orders, etc.. • Rotating condensates: instantaneous measurement of the correlation functions in the rotating frame. • Correlation functions near continuous phase transitions. • Systems away from equilibrium.

31. gap  tuning parameter  Universal adiabatic dynamics across a quantum critical point Consider slow tuning of a system through a critical point. t,   0 Gap vanishes at the transition. No true adiabatic limit! How does the number of excitations scale with  ? This question is valid for isolated systems with stable excitations: conserved quantities, topological excitations, integrable models.

32. Substitute into Schrödinger equation. Use a general many-body perturbation theory. Expand the wave-function in many-body basis.

33. Use scaling relations: Find: Uniform system: can characterize excitations by momentum:

34. Caveats: • Need to check convergence of integrals (no cutoff dependence) Scaling fails in high dimensions. • Implicit assumption in derivation: small density of excitations does not change much the matrix element to create other excitations. • The probabilities of isolated excitations: should be smaller than one. Otherwise need to solve Landau-Zeener problem. The scaling argument gives that they are of the order of one. Thus the scaling is not affected.

35. Breakdown of adiabaticity: From t and we get  Simple derivation of scaling (similar to Kibble-Zurek mechanism): In a non-uniform system we find in a similar manner:

36. Example: transverse field Ising model. There is a phase transition at g=1. This problem can be exactly solved using Jordan-Wigner transformation:

37. Spectrum: Critical exponents: z==1  d/(z +1)=1/2. Correct result (J. Dziarmaga 2005): Other possible applications: quantum phase transitions in cold atoms, adiabatic quantum computations, etc.