arxiv: v1 [cond-mat.other] 5 Dec 2008

Transcription

1 Phase-diffsion dynamics in weakly copled Bose-Einstein condensates Erez Bokobza a, Maya Chchem b, Doron Cohen b, and Amichay Vardi a Departments of Chemistry a and Physics b, Ben-Grion University of the Negev, P.O.B. 653, Beer-Sheva 845, Israel We stdy the phase-sensitivity of collisional phase-diffsion between weakly copled Bose-Einstein condensates, sing a semiclassical pictre of the two-mode Bose-Hbbard model. When weakcopling is allowed, zero relative phase locking is attained in the Josephson-Fock transition regime, whereas a π relative phase is only locked in Rabi-Josephson point. Or analytic semiclassical estimates agree well with the nmerical reslts. arxiv:8.4v [cond-mat.other] 5 Dec 8 Bose-Einstein condensates (BECs) of dilte, weaklyinteracting gases are crrently sed as a testing grond for the qantitative stdy of varios condensed-matter models. In addition to offering a remarkable degree of controllability, they open the way for exploring noneqilibrim dynamics, far beyond small pertrbations of the grond state. Whereas this regime is inaccessible in the eqivalent condensed matter realizations, de to the high Fermi energy, highly excited states are natrally prodced in BEC experiments and their dynamics can be traced with great precision. One sch example is the phase-diffsion [,, 3, 4, 5, 6] between two BECs prepared in a coherent state and conseqently separated. The physics of this process is captred by the two site Bose-Hbbard Hamiltonian (BHH). Defining â i, â i as bosonic annihilation and creation operators for particles in condensate i =,, the corresponding particle nmber operators are ˆn i = â iâi, and the BHH takes the form [7, 8], Ĥ = εˆl z J ˆL x + U ˆL z, () where ˆL x = (â â + â â)/, ˆL y = (â â â â)/(i), and ˆL z = (ˆn ˆn )/ are the generators of an SU() Lie algebra. Bias potential, copling, and collisional interaction energies are denoted as ε, and J, and U, respectively. We have eliminated c-nmber terms that depend on the conserved total particle nmber N = ˆn + ˆn. Below we se for representation the Fock space basis states l, m z, which are the joint eigenstates of L and L z, with l = N/. Experimental procedres allows to prepare the system in any desired SU() coherent state θ, ϕ = e iϕˆl z e iθˆl y l, l z state. The θ rotation is realized by strong initial copling J NU, while the phase ϕ is tned by switching the bias ε, ths indcing coherent phase oscillations [4, 9]. In particlar it is possible to prepare all the eqal-poplation coherent states of the form π, ϕ l l m= l ( e iϕ ) ( ) / l+m l l, m l + m z (). In phase-diffsion experiments, the preparation stage is followed by a sdden separation of the condensates, so as to obtain two eqally poplated modes of a symmetric (ε=) doble well with a definite relative phase ϕ. De to the interaction terms, different L z Fock states oscillate with different freqencies, and the evoltion nder the Hamiltonian () with = NU/J leads to the loss of single-particle coherence, qantified by the mean fringe-visibility of ensemble-averaged interference patterns g () = L x /l. For flly separated condensates, fringe visibility decays as g () (t) = e (t/t d) with t d = (U l) [, 3, 6]. Long time revivals are obtained at t r = π/u, as confirmed experimentally [3]. This behavior is insensitive to the initial phase ϕ, and merely reflects the Binomial Gassian-like distribtion of the occpation which is implied by Eq. (). In reality we still may have some inter-mode copling J, which is small compared with the collisional energy NU [4, 5, 6, 9]. Recent experiments have explored the loss of single-particle coherence nder sch circmstances [5], reslting in the decay to a non-vanishing eqilibrim vale of g (). Other experiments have fond ϕ- sensitivity in the fll merging of two separated condensates [9], where a relative ϕ=π phase led to significant heating losses compared to ϕ= phase. The common expectation in these experiments is that small copling between the condensates will lead to their phase-locking and sppress the loss of single-particle coherence. Motivated by this experimental interest we here stdy the ϕ dependence of the fringe visibility evoltion when the two condensates remain weakly copled dring the hold time. In Fig. we plot the nmerically calclated for both the π/, and the π/, π preparations. For the ϕ= preparation, the expected phase-locking is obtained even by a very weak copling, leading to a nonvanishing asymptotic vale of the fringe visibility, with some small oscillations. However, if the system is prepared in the ϕ=π state, phase locking is not attained if is moderately large, and the fringe visibility exhibits apparently complex dynamics with a somewhat different characteristic timescale. Or prpose in the present work is to provide both qalitative and qantitative analysis of this phase-sensitivity of the phase-diffsion process. g () Otline In order to have this presentation self contained we start with a terse review of the spectrm sing a phase space langage. Then we extend the analysis

2 (e) (f) () g (g).95.9 (h) ν J t ν J t FIG. : (Color online) Evoltion of fringe visibility with l = 5 (N = particles), starting from the spin coherent states π/, (a-d) and π/, π (e-h). The copling parameter is = 4 (a,e), 3, (b,f) (c,g), and (d,h). The evoltion of the π/, state is characterized by small oscillations at the Josephson freqency ν J = NUJ/π, and a transition to coherent evoltion taking place between = N and = N. By contrast for π/, π cohererent dynamics is only approached for < and the characteristic freqency has a log correction. to inclde the dynamics in the separatrix region, which is relevant for the nderstanding of the evoltion of the ϕ=π preparation. This is contrasted with the case of ϕ= preparation. In particlar we discss the characteristics of the flctations which are displayed in Fig.. Interaction Regimes It is instrctive to se a Wigner fnction ρ(θ, φ) for the representation of the qantm states []. The Wigner fnction that represents a θ, ϕ coherent state resembles a minimal Gassian centered at the corresponding angle on the BHH s spherical phase space. Using semiclassical reasoning one can arge that the Wigner fnctions of the BHH eigenstates (with ε = ) are concentrated along the contor lines of the Gross-Pitaevskii classical energy fnctional, E(θ, ϕ) = NJ [ ] (cosθ) sin θ cosϕ (3) Frthermore, within the framework of the WKB approximation the qantization of the energy is implied by the condition A(E n ) = (4π/N)n where A(E) is the phase-space area enclosed by a fixed energy E contor, and 4π/N is the Planck cell. [For simplicity we disregard Maslov indexes, and below we measre A(E) in Planck cell nits]. The phase space pictre implies that one has to distingish between three regimes according to the vale of the dimensionless interaction parameter N U/J [7, 8]: The linear, weak-interaction Rabi regime < ; the intermediate Josephson regime < < N ; and the extreme strong-interaction Fock regime > N. FIG. : (Color online) Schematic drawing of the classical phase-space trajectories and the corresponding semiclassical spectrm in the Josephson regime ( = ). Dotted ble lines mark the low-energy sea trajectories, leading to eqally-spaced levels, solid green lines mark the trajectories within the nonlinear islands, reslting in a Um high-energy spectrm of doblets. The separatrix trajectory and corresponding energies are marked by thick solid red lines. In the Josephson regime the spherical phase-space is split by a figre-eight separatrix trajectory (see Fig ), to a sea of Rabi-like trajectories (ble dotted lines), and two interaction-dominated, nonlinear islands (solid green lines). In the Fock regime the phase-space area of the sea becomes smaller than Planck cell, and therefore effectively disappears. Spectrm The energy landscape for, as implied by Eq.(3), consists of sea levels that extend from the bottom energy E = lj, and of island levels that extend p to the pper energy E + l U. The border between the sea and the islands is the separatrix energy E x = lj. The oscillation freqency ω(e) = [A (E)] arond the minimm-energy (stable) fixed point (π/, ) is the plasma / Josephson freqency ω J ω(e ) = (J + NU)J NUJ (4) If we had < then the opposite phase-space fixed point (π/, π) wold be also stable with ω + = (J NU)J. Bt for > it bifrcates, and replaced by an nstable fixed point at the separatrix energy, accompanied by the two twin stable fixed points (arcsin(/), π), located within the islands. For these stationary points approach the poles (arcsin(/) =, π) and the oscillation freqency becomes ω + NU. The significance of the varios freqencies is implied by the WKB qantization condition. At low-energy we have a non-degenerate set of Josephson levels with spacing ω J. De to the nonlinearity the spacing becomes smaller as one approaches the separatrix energy (see frther analysis below), while the high-energy levels are dobly-degenerate with spacing ω + (the high energy part of the spectrm consists of Um energy levels whose Um spacing approaches ω + as

3 3 m l). Nonlinearity The motion in the vicinity of the separatrix is highly non-linear. This is implied by the non-linear variation of the phase space area near the separatrix: A(E) A(E x ) = E E x ω J log NJ E E x (5) We stress that Planck cell nits are sed for A(E). On the basis of this expression, the WKB qantization condition implies that the level spacing at the separatrix is finite and given by the expression ω x = [ ( )] N log ω J (6) It is important to notice that in the strict classical limit this freqency becomes zero, while for finite N it differs from the Josephson freqency only by a logarithmic factor. The non-linearity can be characterized by a parameter α that reflects the E dependence of the oscillation freqency ω(e), and hence pon WKB qantization it reflects the dependence of the level spacing E n+ E n ω(e n ) on the rnning index n. The implied definition of this parameter and its vale in the vicinity of the separatrix are: α(e) [ ( )] dω N ω dn = ω(e) A (E) log (7) Note that this parameter is meaningfl only deep in the Josephson regime, where it is less than nity. As shown below, the parameter α controls the amplitde of the flctations in the case of a ϕ=π preparation. Strangely enogh, in the classical limit (N with fixed ) the non-linear effect becomes more prononced as far as Eq.(6) is concerned, bt weaker with regard to Eq.(7). Eigenstates In the Rabi regime, the eigenstates of Hamiltonian () approach the ˆL x eigenstates l, m x, with eigenvales Jm proportional to the relative nmber difference between the even and odd sperpositions of the modes. In particlar, the lowest eigenstate l, l x (Fig. 3a) and the highest eigenstate l, l x (Fig. 3b), are the coherent states π/, and π/, π, respectively, with binomial m distribtions in the l, m z basis, approaching the normal (Gassian/Poissonian) distribtion for large l. As increases beyond nity, a transition is made into the Josephson regime. The relative-nmber variance L z of the grond state decreases continosly in this regime from its coherent N/ vale (Fig. 3d). However, the coherence g () of the grond state remains close to nity (Fig. 3c) and the relative-phase is welldefined throghot the Josephson regime, jstifying the se of mean-field theory to depict the Josephson dynamics of grond state pertrbations. The coherence of the grond state is only lost in the Fock regime, where the P m () g m Rabi Josephson Fock P m log ( L z ) 4 N/4 4 5 m N /4 Rabi Josephson Fock FIG. 3: (Color online) Grond and most excited eigenstates of the qantized Josephson model with l = 5, as a fnction of the interaction parameter : nmber distribtion of the grond state, nmber distribtion of most excited state, fringe-visibility for grond state (dashed) and most excited state (dash-dot), relative-nmber variance for grond (dashed) and most excited (dash-dot) states. tre grond state approaches the relative site-nmber state l, z. ϕ= preparation The above discssion implies that in the case of a coherent preparation π/,, any weak copling beyond J > U/N wold sffice to lock the relative phase and to prevent the loss of single-particle coherence [5]. Using a phase space langage this means that the ϕ= preparation is located at the minimm of the sea region, arond a stable fixed point (π/, ), and resembles the grond state of the Hamiltonian. The expected small-oscillation freqency for the evoltion of this coherent state is ω J. ϕ=π preparation For the coherent preparation π/, π the pictre is qite different. In the Rabi regime, this state coincides with the l, l x most excited eigenstate, and hence it is possible to have phase locking arond the stable fixed point (π/, π), where the characteristic freqency is ω +. However, as soon as >, this fixed point bifrcates, and the Wigner fnction of the eigenstates has to eqally poplate the two twin islands in phase space, hence forming a cat state []. The corresponding m distribtion becomes bimodal, as shown in Fig. 3b. This crossover is reflected by the sharp increase in the relative-nmber variance to the cat-state vale L z = N/, and by the eqally sharp decrease of its single-particle coherence. The evoltion of the coherent preparation π/, π in the vicinity of the nstable point is frther analyzed below. Dynamics The phase-sensitivity of phase-diffsion is clearly explained by the semiclassical approach. The coherent state π/, is a sperposition of a few lowenergy states. Since the level-spacing is fixed on the Josephson freqency ω J, the fringe visibility carries ot clean harmonic oscillations. By contrast, the state 5

4 4 P n E / J n (e) (f) (g) (h) E / J n FIG. 4: (Color online) Nmerical LDOS expansion of the states π/, (a-d) and π/, π (e-f) in the eigenstates basis, for l = 5, = (a,e), (b,d), (c,g), and (d,h). Bar labels denote the rnning index n = l,..., l. The π/, state is a good approximation to an eigenstate p to N whereas the π/, π state corresponds to narrow distribtion of eigenstates, starting from. Solid lines correspond to the semiclassical estimates of Eq. (8) in -, and Eq. (9) in (e)-(h). π/, π consists of levels arond the separatrix energy, which are not eqally spaced, reslting in the qasiperiodic oscillations of Fig. (e)-(h). The observed timescale of these oscillations agrees well with ω x. In order to have a qantitative estimate for the amplitde of the oscillation we have to expand the initial state ψ in the E n basis. A reasonable estimate for the envelope fnction P(E n ) = E n ψ = trace(ρ (n) ρ (ψ) ) can be obtained sing a semiclassical approximation. The Wigner fnction of the nth eigenstate is approximated by a microcanonical distribtion ρ (n) (ϕ, θ) = ω(e n )δ(e(θ, ϕ) E n ) and the coherent state is approximated by a minimal Gassian. With appropriate approximations for the phase space integration gives for the ϕ= preparation [( P(Ẽ) = I ) ] Ẽ e (+/)Ẽ, (8) with Ẽ (E E )/J, and for the ϕ=π preparation P(Ẽ) = ( [( ω(e) )K + ) ] π ω J Ẽ e ( /)Ẽ, (9) with Ẽ (E E x)/j. In the above, I and K are the modified Bessel fnctions. The comparison of this analytical estimate with the exact nmerical diagonalization is illstrated in Fig. 4. De to the even parity of both preparations ( π/, is always even and π/, π has ( ) N parity), the occpation of odd n states vanishes and the occpation of the even n states is twice the semiclassical estimate. Using these expressions we can estimate the participation ratio of the preparation, which comes ot M NU/ω J. Flctations Regarding ψ as a sperposition of M energy states we can estimate the flctations sing the following reasoning: In the linear case, if the energy levels are eqally spaced, there is only one M = basic freqency and the motion is strictly periodic. If the energies were qasi-random, then there are M = M independent freqencies, and the motion is qasiperiodic with relative variance /M. Bt what we have for the π preparation are energy levels that are characterized by a nonlinearity parameter α. Then, if the nonlinearity parameter is not very small, the effective nmber of independent freqencies is M = αm, which interpolates between the linear and the qasi-random estimates. This estimate agrees with the nmerical simlation in Fig.. Smmary The vast majority of recent BEC interference experiments related to phase-diffsion, are carried ot in the Josephson regime. Within it, the dynamics of single-particle coherence is expected to strongly depend on the initial relative-phase between the partially separated condensates. Zero relative phase leads to phase-locking even for as large as N, whereas a π relative phase preparation reslts in the loss of fringevisibility all the way down to. Using a semiclassical phase-space pictre we related this behavior to the classical phase-space strctre and the ensing WKB spectrm, obtaining analytical expressions for the envelope fnctions of these two coherent preparations in the eigenstate basis, for in the Josephson regime. These expansions agree well with nmerical calclations and enable the accrate prediction of the amplitde of fringevisibility oscillations. This work was spported by the Israel Science Fondation (Grant 58/7), and by a grant from the USA-Israel Binational Science Fondation (BSF). [] Y. Castin and J. Dalibard, Phys. Rev. A 55, 433 (997); M. Lewenstein and L. Yo, Phys. Rev. Lett. 77, 3489 (996); E. M. Wright, D. F. Walls and J. C. Garrison Phys. Rev. Lett. 77, 58 (996); M. Javanainen and M. Wilkens, Phys. Rev. Lett. 78, 4675 (997). [] A. Vardi and J. R. Anglin, Phys. Rev. Lett. 86, 568 (); J. R. Anglin and A. Vardi, Phys. Rev. A 64, 365 (); Y. Khodorkovsky, G. Krizki, and A. Vardi, Phys. Rev. Lett., 43 (8) [3] M. Greiner, M. O. Mandel, T. Hänsch, and I. Bloch Natre 49, 5 (). [4] T. Schmm et al., Nat. Phys., 57 (5). [5] S. Hofferberth et al., Natre (London) 449, 34 (7). [6] G.-B. Jo et al., Phys. Rev. Lett. 98, 347 (7). [7] Y. Makhlin, G. Schön, and A. Shnirman, Rev. Mod. Phys. 73, 357 (); R. Gati and M. K. Oberthaler, J. Phys. B 4, R6 (7). [8] Gh-S. Paraoan et al., At. Mol. Opt. Phys. 34, 4689 (); A. J. Leggett, Rev. Mod. Phys. 73, 37 (). [9] G.-B. Jo et al., Phys. Rev. Lett. 98, 84 (7). [] G. S. Agarwal, Phys. Rev. A. 4, 889 (98). J. P. Dowling, G. S. Agarwal,W. P Schleich, Phys. Rev. A 49, 4 (994).

5 [] P. J. Y. Lois, P. M. R. Brydon, and C. M. Savage, Phys. Rev. A 64, 5363 (); A. Micheli et al., Phys. Rev. A 67, 367 (3); Y. P. Hang and M. G. Moore, Phys. Rev. A 73, 366 (6). 5