# Adiabatic Splitting, Transport, and Self-Trapping

of a Bose-Einstein Condensate in a Double-Well Potential

###### Abstract

We show that the adiabatic dynamics of a Bose-Einstein condensate (BEC) in a double well potential can be described in terms of a dark variable resulting from the combination of the population imbalance and the spatial atomic coherence between the two wells. By means of this dark variable, we extend, to the non-linear matter wave case, the recent proposal by Vitanov and Shore [Phys. Rev. A 73, 053402 (2006)] on adiabatic passage techniques to coherently control the population of two internal levels of an atom/molecule. We investigate the conditions to adiabatically split or transport a BEC as well as to prepare an adiabatic self trapping state by the optimal delayed temporal variation of the tunneling rate via either the energy bias between the two wells or the BEC non-linearity. The emergence of non-linear eigenstates and unstable stationary solutions of the system as well as their role in the breaking down of the adiabatic dynamics is investigated in detail.

###### pacs:

03.75.Kk, 03.75.Lm## I Introduction

Bose Einstein condensates (BEC) in double well potentials have drawn a lot of attention both theoretically and experimentally for the possibilities they offer to study fundamental quantum mechanical effects at the macroscopic level as well as for potential applications like interferometry, high precision measurements or thermometry split (). Most of the experimental realizations lead to the splitting of the condensate into two independent parts. Nevertheless, recently, also weakly linked parts of a BEC in a double well potential forming a single Josephson junction josephson () have been achieved BEC_single_josephson (). In contrast to Josephson junctions realized in superconductors and superfluids classical_josephson (), in BEC the non-linear interatomic interactions play a crucial role. In the presence of the non-linearity two dynamical regimes have been predicted: (i) anharmonic Josephson oscillations anharmonic (), if the initial population imbalance of the two wells is below a critical value and (ii) macroscopic quantum self-trapping self-trapping () i.e, inhibition of large amplitude Josephson oscillations above the threshold for the population imbalance. This threshold corresponds to the population imbalance for which the difference between the two on-site interaction energies becomes larger than the tunneling energy splitting. Both dynamical regimes have been explored experimentally in a single Josephson junction BEC_single_josephson () and in arrays josephson_arrays ().

Recently, a lot of attention has been devoted to explore techniques to coherently control the non-linear dynamics of a BEC in double well potentials by modulating in time either the potential time_potential () or the non-linearity time_nonlinearity (). The two main studied regimes are the non-linear Landau-Zener WuNiu00 (); Liu03 (); The06 (); Gra06 (); Wan06 (); Iti07 (); Zha08 () and the Rosen-Zener DYe08 () regimes. In the former the tunneling rate is fixed while the energy bias is linearly varied in time, in the later the energy bias is varied in time while the tunneling rate is switching on/off following e.g., a temporal Gaussian profile. Both these regimes have been deeply investigated in double WuNiu00 (); Liu03 (); The06 (); Iti07 (); Zha08 (); DYe08 () and triple well Gra06 (); Wan06 () potentials yielding a wide variety of dynamical scenarios ranging from robust population transfer to quantum blocking, non-linear oscillations and, even, breaking down of the adiabatic dynamics.

Following a different perspective, there have been several recent proposals to coherently manipulate single atoms NosPRA04 (); NosOC06 (); Das08 () and BECs Gra06 (); Rab08 (); Nesterenko08 () in triple-well potentials by adiabatically following a particular energy eigenstate of the system, the so-called spatial dark state. This spatial dark state only involves the two ground states of the extreme traps in a close analogy to the well known quantum optical Stimulated Raman Adiabatic Passage (STIRAP) technique Bergmann1998RMP70 (). Accordingly, these tools were named Three-Level Atom Optics (TLAO) techniques NosPRA04 (). In this paper, following the work by Vitanov and Shore Vitanov2006PRA73 () on the adiabatic passage on two-level atoms, we extend the TLAO techniques to the two-level matter wave case showing that the adiabatic dynamics of a BEC, in a double well potential, can be described in terms of a dark variable resulting from the combination of the population imbalance and the spatial atomic coherence. This dark variable can be tailored varying in time the tunneling rate, the energy bias, and the non-linearity with the goal of (i) adiabatically splitting of the BEC, (ii) achieving complete BEC transfer from one well to the other; and (iii) preparing of an adiabatic self-trapping state.

The paper is organized as follows. In Section II we present the physical system consisting of a BEC in a double-well potential. We assume the two-level approximation and describe the BEC dynamics in terms of a dark variable. The conditions for the adiabatic control of tunneling by means of this dark variable are discussed in Section III from a non-linear dynamics perspective. In Section IV we present detailed numerical simulations on the adiabatic splitting, transport, and trapping of a BEC. Section V summarizes the main conclusions of the paper and briefly discusses the validity of the two-mode approximation. It also presents some possible extensions of the present work such as its formulation in second quantization.

## Ii Model

We consider a BEC trapped in a double well potential (Fig. 1(a)), whose dynamics at zero temperature is described by the time-dependent Gross-Pitaevskii equation (GPE) ():

(1) |

where is the external trapping potential and the non-linearity is given by , with the total atom number, the -wave scattering length, and the atomic mass. Assuming the two-level approximation self-trapping (); two_lev_approx (), the wave function or classical order parameter of the BEC under study can be written as with satisfying . accounts for the ground state of the corresponding isolated trap. The amplitudes and obey the non-linear two-mode dynamical equations given by:

(2) |

with the Hamiltonian:

(3) |

where are the atomic self-interaction energies, is the tunneling rate between the two wells, and are the on-site energies. In terms of the overlaps, the former parameters read:

(4) | |||||

(5) | |||||

(6) |

with .

Notice that the two-mode approximation assumes that the parameters of the problem (tunneling rate, non-linear interaction, and energy bias) can be varied independently. In general, however, the temporal modification of any of these parameters will result in the modification of the spatial mode functions affecting the rest of the parameters. From the experimental point of view, a double-well optical potential can be created with well separations on the range of few microns by means of two focused laser beams dw1 (), or by the superposition of a 3D crossed beam dipole trap with a 1D optical lattice BEC_single_josephson (). In the later, the intensity of the standing wave field that creates the optical lattice and its displacement with respect to the center of the dipole trap could be used to manipulate both the tunneling and the energy bias. In addition, the scattering length could be controlled playing around of a Feshbach resonance Bauer09 ().

Within the density matrix formalism, assuming , and taking as the energy bias between the two wells, the coherent dynamics of the BEC in the two-well potential can be written as follows LeeHai ():

(7) |

where we have introduced the real-valued variables , , being and with the corresponding trap populations and spatial coherence (see Fig. 1(b)). Note that the conservation of the norm implies .

Inspired by the work of Vitanov and Shore Vitanov2006PRA73 () on STIRAP in a system of two internal atomic levels, we extend their proposal to the external degrees of freedom of matter waves. It is easy to verify that Eqs. (7) are equivalent to the discrete non linear Schrödinger equation (DNLSE) for a BEC in a triple well potential (see, for instance, Ref. Gra06 ()) under the following identifications: , , , and , where with are the probability amplitudes to find the BEC in the left, middle, or right trap, respectively. and denote, respectively, the tunneling interaction between the left-middle and middle-right traps. Note however that for the two level system the variables and are real quantities, while for a three level system the variables are, in general, complex numbers.

In the TLAO techniques within the three-level approximation being the left and right levels resonant, the transfer process is based on adiabatically following one of the three energy eigenstates of the system , with the mixing angle defined as , and and being the ground states of the left and right wells. In the two-level system under investigation, the analogue state will be given by a combination of the population difference and the coherence , with the mixing angle given by . This analogy opens the possibility to extend the TLAO techniques to two-level systems, by appropriately engineering the time dependence of the tunneling rate , the energy bias , and/or the non-linear interaction . Note that the non-linear interaction parameter can be modified in time by the temporal variation of the scattering length using either magnetic magnetic_FR () or optical optical_FR () Feshbach resonances or varying the trap frequency, leading to a modification of the BEC spatial profile according to Eq. (5). For a review on the manipulation of Feshbach resonances see ReviewFeshbach ().

## Iii Adiabatic Control of Tunneling

In this section we study different scenarios for the coherent control of the external degrees of freedom of a BEC by adiabatically following the dark variable . In particular, we show how to adiabatically split, transport, and inhibit tunneling of a BEC via the matter wave analogues of both STIRAP and double-STIRAP techniques in two-level systems.

Let us assume that the BEC is initially prepared in the left trap with . If so, note then that meaning and will be the initial state. Then, to coherently split the condensate, should vary adiabatically from to corresponding to and (see Fig. 1(b)) by appropriately changing , , and in time. To transfer the BEC from the left to the right trap, the mixing angle should be slowly increased up to . For the adiabatic self-trapping state case, the evolution will consist in varying from to e.g., and then back to . In all cases, in order to guarantee that the BEC follows the dark variable during the whole process two conditions must be fulfilled:

Condition 1: Adiabaticity criteria. , , and should be smoothly varied in time to adiabatically follow the dark variable which, in turn, means that the mixing angle must be slowly changed from to its expected final value. Gaussian profiles for the temporal variations of the control parameters are assumed:

(8) | |||||

(9) | |||||

(10) | |||||

with either the energy bias , or the BEC non-linearity , preceding the tunneling interaction , i.e., the counter intuitive sequences or will be assumed. is a switch that takes the value for the two-level matter wave analogue of the STIRAP sequence and () for the symmetric (antisymmetric) double-STIRAP sequences. Adiabaticity means that, at any time, should be much smaller than the energy separation between the selected eigenstate and the one energetically closest. For a weak non-linear interaction, , this adiabaticity condition reads . It is worth to remark that for , is not a parameter of the system but a dynamical variable since it contains the population difference in its definition.

Condition 2: Avoiding adverse bifurcation points. For large enough values of the non-linearity, the interaction between the atoms of the BEC results in a non-linear temporal coupling producing additional non-linear stationary states yielding loop structures and a rich variety of level crossing scenarios Gra06 (); WuNiu00 (); TL2 (); TL3 (). In the matter wave STIRAP case for a triple-well potential Gra06 (), it has been shown that even in the adiabatic limit given by , the appearance of non-linear stationary states breaks down, in some cases, the adiabatic evolution. To address this issue in the double-well potential, we start first looking for those critical values giving rise to non-linear stationary states by solving the eigenvalue equation:

(11) |

with given in Eq. (3) and being the chemical potential. After some algebra one obtains the following fourth-order eigenvalue equation for :

(12) |

For and () this quartic equation gives two (four) real roots WuNiu00 ().

Additional information on the BEC dynamics can be obtained looking for the stationary solutions of the density matrix equations (7) and analyzing their stability (see Appendix A). The stationary solutions can be written as with and . The energy of these stationary solutions, up to four, is given by Eq. (12). Note that, as the dynamics is conservative, the eigenvalues of the linear stability matrix must satisfy for each stationary solution. In fact, a Linear Stability Analysis (LSA) around any of these solutions yields one eigenvalue equal to zero together with either two pure imaginary eigenvalues corresponding to a stable fixed point (or center) or two real eigenvalues accounting for an unstable saddle point. In the limit of a large non-linear interaction, i.e., for , the system presents four stationary solutions, one of them being an unstable saddle point while the other three are stable elliptical ones.

If the number of stationary solutions remains constant during the whole adiabatic dynamics, the system will successfully follow the selected energy eigenstate. However, depending on the interplay between the non-linear interaction and the tunneling rate, during the dynamics the number of stationary solutions will change through a bifurcation point from four to two or vice versa. If the selected stationary solution is not involved in the bifurcation, the system will again adiabatically follow the corresponding eigenstate. On the opposite case, whether the adiabatic dynamics will be affected will depend on the particular bifurcation scenario:

Case 1. From 4 to 2 stationary solutions through a backward pitchfork bifurcation point in which the stationary solution to be followed merges two more solutions (one unstable) and yields one stable solution. In this case, the system will evolve adiabatically. In what follows, this case will be termed PB42 case, even if the pitchfork bifurcation is perturbed and becomes imperfect.

Case 2. From 2 to 4 solutions through a forward pitchfork bifurcation point in which the stationary solution to be followed becomes unstable and yields two stable solutions corresponding, in the limit of vanishing tunneling, to . In this case, the system after reaching the bifurcation point will split into a combination of the two stationary solutions. In this case, named PB24 in what follows, the adiabatic dynamics will, in general, break down.

Case 3. From 4 to 2 solutions through a saddle-node bifurcation in which the stationary solution to be followed is annihilated by an unstable solution. The adiabatic dynamics breaks down even in the adiabatic limit. After reaching the bifurcation the system will split into a combination of two non-degenerate energy eigenstates and, therefore, will oscillate at the corresponding energy splitting. Case SNB42 in what follows.

## Iv Numerical Simulations

This section is devoted to the numerical simulations of the BEC adiabatic splitting, transport, and self-trapping by means of the temporal variation of the energy bias, the non-linearity, and the tunneling rate. We will show, for the most relevant cases, the eigenvalues and stationary states predicted by Eqs. (12) and (7) as well as their linear stability. The main goal of these numerical simulations will consist in illustrating the different dynamical scenarios to adiabatically control the BEC while looking for those parameter values that prevent the system to reach the two unwanted bifurcation cases PB24 and SNB42 previously described.

### iv.1 Adiabatic Splitting of a BEC

We start first by fixing the non-linear interaction control parameter to a constant value while, counter intuitively, time-varying the energy bias and the tunneling rate (see the inset of Fig. 2(b))) with the goal of achieving equal population in the two wells, i.e., starting with the BEC being in the left trap, i.e., , or, equivalently, from . Fig. 2(a) shows , after the numerical integration of Eqs. (7) with and the rest of parameters values given in the figure caption. Clearly, there is a large set of parameters where the adiabatic splitting process takes place with a high fidelity (see the plateau in Fig. 2(a)). Fig. 2(b) shows the plateau for and . To give more insight into the dynamics of the system, we have divided the non-linear interaction range studied in five regions, from to , as shown at the top of Fig. 2(b). In the following we will discuss in detail the dynamics for each of these regions. Note that in all cases, the dynamics starts and ends with which, for , implies that for the system presents four fixed points or stationary solutions. For intermediate times and large enough tunneling rates, there are only two stationary solutions.

In Fig. 3 we show the temporal dynamics corresponding to the parameter values of case (i) in region (c) of Fig. 2(b) with (see Fig. 3(a)). The adiabatic energy eigenvalues and stationary solutions are shown in Figs. 3(b) and (c), respectively, with the dashed curve accounting for the unstable ones. Short arrows in Figs. 3(b) and 3(c) indicate the stationary solution to be adiabatically followed. The dynamical variable , after integration of Eqs. (7), is plotted in Fig. 3(d). For this set of parameter values the adiabaticity condition is fulfilled and the selected eigenstate is involved in a single bifurcation, at , corresponding to the PB42 case previously described. After the bifurcation point, the system follows the eigenstate that, at the end of the process, yields and . This non-linear dynamical scenario holds for the whole plateau, region in Fig. 2(b), where the adiabatic splitting succeeds.

For the parameter range corresponding to region in Fig. 2(b), the final state of the BEC strongly depends on the parameter values. For a slight modification of , the final state alternates between and accounting for the BEC being located at the right or left trap, respectively. To understand the origin of this behavior, we plot in Fig. 4 the detailed dynamics for corresponding to case (ii) in Fig. 2(b). At , the solution that the system is adiabatically following is annihilated with an unstable one through a saddle-node bifurcation, i.e., case SNB42 previously described. At this point, the system splits into two non-degenerate components (being the largest one the energetically closest stable solution) and starts to oscillate at the frequency corresponding to the energy separation of the two solutions. At , the system reaches a bifurcation point of the PB24 type. This bifurcation yields two energetically degenerated stable solutions corresponding to the two final states . Thus, at , the largest component of the system chooses one out of the two new stable solutions. We have numerically verified that the selected solution strongly depends on the previous oscillatory dynamics.

For corresponding to regions and of Fig. 2(b), there are four stationary solutions during the whole dynamics that do not present any bifurcation. In this case, the non-linear coupling forces the adiabatic evolution from to instead of evolving to the desired . Then, the initial and the final state of the system coincide.

Finally, for , region in Fig. 2, the energy eigenstate to be followed reaches, during the ramping down of the tunneling interaction , a pitchfork bifurcation of the PB24 type. After this bifurcation point, the two stable stationary solutions are energetically degenerated and, therefore, the system splits into a non-oscillatory combination of them that for () yields ().

Up to now, we have discussed the possibility to adiabatically split a BEC by fixing and time varying and . Note, however, that it is also possible to adiabatically split the BEC by appropriately time varying the non-linear interaction and the tunneling rate with a fixed energy bias. Thus, for , , and the rest of the parameters given in the figure caption, Fig. 5 shows as a function of the time delay between the modulation of the tunneling rate and the non-linear interaction. We again obtain a plateau of parameter values, , where the adiabatic splitting, corresponding to , becomes successful. A detailed analysis of the temporal dynamics reveals that, as in Fig. 2, the plateau corresponds to the case where the selected solution reaches during its dynamics a single bifurcation point of the PB42 type.

### iv.2 Adiabatic Transport

To adiabatically transfer the BEC from the left to the right trap, the mixing angle must be slowly varied from up to . Taking and the rest of parameters as in Fig. 2, Fig. 6 shows the population difference at the end of the antisymmetric double STIRAP process (shown in the inset) as a function of the non-linear interaction . Clearly, there is a wide parameter setting, the plateau , where the BEC transfer process takes places with a high fidelity. In the parameter region , the system starts at following an energy eigenstate of the system that does not involve any bifurcation point. For the system crosses first a PB42 bifurcation to later on reaching a PB24 bifurcation point that, for these parameters, yields eventually . In regions and the non-linearity yields forcing the initial and final state to be the same.

### iv.3 Adiabatic Self Trapping

It is also possible to completely inhibit the BEC transport between the two wells by means of the matter wave analogue to the double STIRAP process (see inset of Fig. 7). Starting with the BEC located at the left trap, , we show in Fig. 7 the population difference at the time when the tunneling rate is maximum, , and at the end of the double-STIRAP sequence, , as a function of the energy bias . The adiabatic self-trapping process succeeds for a wide domain of parameter values, except for the region . In this last region, the solution to be followed reaches during its dynamics a saddle-point bifurcation while, for the rest of values of the energy bias, the selected solution does not involve any bifurcation point.

## V Conclusions

By means of a dark variable that results from the combination of the population imbalance and the spatial atomic coherence, we have investigated in detail the adiabatic dynamics of a BEC in a double well potential in the framework of the two-level approximation. We have shown that it is possible to robustly split, transport or trap a BEC by appropriate temporal variation of either the energy bias or the non-linear interaction together with the tunneling rate. All these proposals have been studied from a non-linear dynamics perspective deriving the stationary solutions of the system, evaluating their linear stability, and discussing the bifurcation scenarios. For the eventual implementation of the techniques here discussed, we want to highlight the recent work by M. Bauer et al. Bauer09 () reporting an accurate control on magnetic Feshbach resonances by means of laser light.

In order to check the validity of the two-mode approximation, we have numerically integrated the one-dimensional Gross-Pitaevskii (GP) equation for the non-linear adiabatic splitting of a BEC in a double-well potential, obtaining the characteristic ’plateau’ of the STIRAP protocol shown in Fig. 2(b). However, a detailed numerical investigation of the GP equation still remains to be performed to validate that the matter wave non-linear STIRAP techniques here derived in the two-mode approximation could be used for matter wave interferometry or coherent transport of a BEC. In performing this analysis, the global landscape provided by the results of the two-mode approximation should be a useful roadmap.

Although it is out of the scope of the present paper, it would be very interesting also to extend the present non-linear matter wave STIRAP techniques to the second-quantization formalism secondq (). Within this formalism, we could elucidate whether the adiabatic splitting of the BEC results in a macroscopic superposition where the entire condensate localizes in one of the wells, i.e. in a NOON state NOON (), or in a perfect 50% spatial splitting of all the atoms. Note that in both cases the final population difference reads =0.

###### Acknowledgements.

We acknowledge support from the Spanish Ministry of Education and Science under contracts FIS2008-02425 and CSD2006-00019, and for the ”Juan de la Cierva” postdoctoral fellowship (C. O.). Support from the Catalan Government under contract SGR2009-00347 is also acknowledged. We thank A. Benseny, G. Birkl, Y. Loika, and G. Orriols for valuable and clarifying discussions.## Appendix A Linear Stability Analysis

The stationary solutions of the density matrix equations (7) take the form with and . Let us consider now a perturbation around any of these solutions in the form of:

(13) | |||||

(14) | |||||

(15) |

where the real part of determines the stability of the solution. The stability is governed by the equation where and is the linear stability matrix defined as:

(16) |

Solving the secular equation , one obtains the following three eigenvalues:

(17) | |||||

(18) |

As expected from a conservative system, for every stationary solution. One eigenvalue is equal to zero while the other two are either pure imaginary corresponding to a stable fixed point (or center), or pure real accounting for an unstable one (or saddle point).

## References

- (1) M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. Durfee, D. M. Kurn, and W. Ketterle, Science 275, 637 (1997); T. G. Tiecke, M. Kemmann, Ch. Buggle, I. Schvarchuck, W. von Klitzing, and J. T. M. Walraven, J. Opt. B 5, S119 (2003); Y. Shin, M. Saba, A. Schirotzek, T. A. Pasquini, A. E. Leanhardt, D. E Pritchard, and W. Ketterle, Phys. Rev. Lett. 92, 150401 (2004); J. Estève, J.-B. Trebbia, T. Schumm, A. Aspect, C. I. Westbrook, and I. Bouchoule, Phys. Rev. Lett. 96, 130403 (2006); B. V. Hall, S. Whitlock, R. Anderson, P. Hannaford, and A. I. Sidorov, Phys. Rev. Lett. 98, 030402 (2007); U. Hohenester, P. K. Rekdal, A. Borzi, and J. Schmiedmayer, Phys. Rev. A 75, 023602 (2007); J. Esteve, C. Gross, A. Weller, S. Giovanazzi, and M. K. Oberthaler, Nature 455, 1216 (2008).
- (2) B.D. Josephson, Phys. Lett. 1, 251 (1962).
- (3) M. Albiez, R. Gati, J. Fölling, S. Hunsmann, M. Cristiani, and M. K. Oberthaler, Phys. Rev. Lett. 95, 010402 (2005).
- (4) P. L. Anderson, and J. W. Rowell, Phys. Rev. Lett. 10, 230 (1963); S. V. Pereverzev, A. Loshak, S. Backhaus, J. C. Davis, and R. E. Packard, Nature 388, 449 (1997); A. K. Sukhatme, Y. Mukharsky, T. Chui, and D. Pearson, Nature 411, 280 (2001).
- (5) J. Javanainen, Phys. Rev. Lett. 57, 3164 (1986); M. W. Jack, M. J. Collett, and D. F. Walls, Phys. Rev. A 54, R4625 (1996); I. Zapata, F. Sols, and A. J. Leggett, Phys. Rev. A 57, R28 (1998).
- (6) G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, Phys. Rev. A 55, 4318 (1997); A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, Phys. Rev. Lett. 79, 4950 (1997); S. Raghavan, A. Smerzi, S. Fantoni, and S. R. Shenoy, Phys. Rev. A 59, 620 (1999).
- (7) F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, Science 293, 843 (2001); F. S. Cataliotti, L. Fallani, F. Ferlaino, C. Fort, P. Maddaloni, and M. Inguscio, New Journal of Physics 5, 71 (2003); T. Anker, M. Albiez, R. Gati, S. Hunsmann, B. Eiermann, A. Trombettoni, and M. K. Oberthaler, Phys. Rev. Lett. 94, 020403 (2005).
- (8) F. Kh. Abdullaev and R. A. Kraenkel, Phys. Rev. A 62, 023613 (2000); M. Holthaus, Phys. Rev. A 64, 011601(R) (2001); C. Lee, W. Hai, L. Shi, X. Zhu, and K. Gao, Phys. Rev. A 64, 053604 (2001); G. F. Wang, L. B. Fu , and J. Liu , Phys. Rev. A 73, 013619 (2006); Q. Xie and W. Hai, Phys. Rev. A 75, 015603 (2007).
- (9) F. K. Abdullaev and R. A. Kraenkel, Phys. Lett. A 272, 395 (2000); F. K. Abdullaev, J. S. Shaari, and M. R. B. Wahiddin, Phys. Lett. A 345, 237 (2005); C. Lee, Phys. Rev. Lett. 102, 070401 (2009).
- (10) B. Wu and Q. Niu, Phys. Rev. A 61, 023402 (2000).
- (11) J. Liu, B. Wu, and Q. Niu, Phys. Rev. Lett. 90, 170404 (2003).
- (12) G. Theocharis, P. G. Kevrekidis, D. J. Frantzeskakis, and P. Schmelcher, Phys. Rev. E 74, 056608 (2006).
- (13) E.M. Graefe, H. J. Korsch, and D. Witthaut, Phys. Rev. A 73 013617 (2006).
- (14) G.-F. Wang, D.-F. Ye, L.-B. Fu, X.-Z. Chen, and J. Liu, Phys. Rev. A 74, 033414 (2006).
- (15) A. P. Itin and S. Watanabe, Phys. Rev. E 76, 026218 (2007).
- (16) Q. Zhang, P. Hänggi, and J. Gong, Phys. Rev. A 77, 053607 (2008).
- (17) D.-F. Ye, L.-B. Fu and J. Liu, Phys. Rev. A 77, 013402 (2008)
- (18) K. Eckert, M. Lewenstein, R. Corbalán, G. Birkl, W. Ertmer, J. Mompart, Phys. Rev. A 70, 023606 (2004); Th. Busch, K. Deasy, and S. Nic Chormaic, Journal of Physics: Conference Series 84 012002 (2007).
- (19) K. Eckert, J. Mompart, R. Corbalán, M. Lewenstein, G. Birkl, Opt. Commun. 264, 264 (2006).
- (20) T. Opatrný and K. K. Das, Phys. Rev. A 79, 012113 (2009).
- (21) M. Rab, J. H. Cole, N. G. Parker, A. D. Greentree, L. C. L. Hollenberg, and A. M. Martin, Phys. Rev. A 77, 061602(R) (2008).
- (22) V. O. Nesterenko, A. N. Nikonov, F. F. de Souza Cruz, and E. L. Lapolli, Laser Phys. 19, 616 (2009).
- (23) K. Bergmann, H. Theuer, and B. W. Shore, Rev. Mod. Phys. 70, 1003-1025 (1998).
- (24) N. V. Vitanov and B. W. Shore, Phys. Rev. A 73, 053402 (2006).
- (25) I. Zapata, F. Sols, and A. J. Leggett, Phys. Rev. A 57, R28 (1998); G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, Phys. Rev. A 55, 4318 (1997); R. W. Spekkens and J. E. Sipe Phys. Rev. A 59, 3868 (1999); S. Giovanazzi, A. Smerzi, and S. Fantoni, Phys. Rev. Lett. 84, 4521 (2000); D. Ananikian and T. Bergeman, Phys. Rev. A 73, 013604 (2006).
- (26) Y. Shin, M. Saba, T. A. Pasquini, W. Ketterle, D. E. Pritchard, and A. E. Leanhardt, Phys. Rev. Lett. 92, 050405 (2004).
- (27) D. M. Bauer, M. Lettner, C. Vo, G. Rempe, and S. Dürr, Nature Physics 5, 339 (2009).
- (28) C. Lee, W. Hai, L. Shi, and K. Gao, Phys. Rev. A 69, 033611 (2004).
- (29) N. V. Vitanov and S. Stenholm, Phys. Rev. A 55, 648 (1997).
- (30) S. Inouye, M. R. Andrews, J. Stenger, H.-J. Miesner, D. M. Stamper-Kurn, and W. Ketterle, Nature 392, 151 (1998); E. A. Donley, N. R. Claussen, S. L. Cornish, J. L. Roberts, E. A. Cornell, and C. E. Wieman, Nature 412, 295 (2001); H. Saito and M. Ueda, Phys. Rev. A 65, 033624 (2002).
- (31) P. O. Fedichev, Yu. Kagan, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 77, 2913 (1996); J. L. Bohn and P. S. Julienne, Phys. Rev. A 56, 1486 (1997); F. K. Fatemi, K. M. Jones, and P. D. Lett, Phys. Rev. Lett. 85, 4462 (2000); M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, G. Ruff, R. Grimm, and J. H. Denschlag, Phys. Rev. Lett. 93, 123001 (2004); R. Ciurylo, E. Tiesinga, and P. S. Julienne, Phys. Rev. A 71, 030701(R) (2005).
- (32) C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, arXiv: 0812.1496 [quant-ph].
- (33) M. Holthaus, Phys. Rev. A 64, 011601(R) (2001).
- (34) R. D’Agosta and C. Presilla, Phys. Rev. A 65, 043609 (2002).
- (35) G.-B. Jo, Y. Shin, S. Will, T. A. Pasquini, M. Saba, W. Ketterle, D. E. Pritchard, M. Vengalattore, and M. Prentiss, Phys. Rev. Lett. 98, 030407 (2007); G.-B. Jo, J.-H. Choi, C. A. Christensen, T. A. Pasquini, Y.-R. Lee, W. Ketterle, and D. E. Pritchard, Phys. Rev. Lett. 98, 180401 (2007); J. Estève, C. Gross, A. Weller, S. Giovanazzi, and M. K. Oberthaler, Nature 455, 1216â1219 (2008).
- (36) D. W. Hallwood, K. Burnett, and J. Dunningham, New Journal of Physics 8, 180 (2006); T. J. Haigh, A. J. Ferris, and M. K. Olsen, arXiv:0907.1333 [quant-phys] (2009); D. W. Hallwood, A. Stokes, J. J. Cooper, and J. Dunningham, New Journal of Physics 11, 103040 (2009).