When attempting to split coherent cold atom clouds or a Bose-Einstein condensate (BEC) by bifurcation of the trap into a double well, slow adiabatic following is unstable with respect to any slight asymmetry, and the wave "collapses" to the lower well, whereas a generic fast chopping splits the wave but it also excites it. Shortcuts to adiabaticity engineered to speed up the adiabatic process through non-adiabatic transients, provide instead quiet and robust fast splitting. The non-linearity of the BEC makes the proposed shortcut even more stable.
PACS numbers:Introduction.-The splitting of a wavefunction is an important operation for matter wave interferometry [1][2][3][4]. It is a peculiar one though, as adiabatic following, rather than being robust, is intrinsically unstable with respect to a small external potential asymmetry [5]. The ground-state wavefunction "collapses" into the slightly lower well so that a very slow trap potential bifurcation in fact fails to split the wave except for perfectly symmetrical potentials. An arbitrarily fast bifurcation may remedy this but at the price of a strong excitation which is also undesired. We propose here a way out to these problems by using shortcuts to adiabaticity that speed up the adiabatic process along a non-adiabatic route. The wave splitting via shortcuts avoids the final excitation and turns out to be signifficantly more stable than the adiabatic following with respect to the asymmetric perturbation. Specifically we shall use a simple inversion method: a streamlined version [6] of the fast-forward technique of Masuda and Nakamura [7] applied to GrossPitaievski (GP) or Schrödinger equations. We have previously found some obstacles to apply the invariantsbased method (at least using quadratic-in momentum invariants [6]) and the transitionless-driving algorythm [8] (because of difficulties to implement in practice the counter-diabatic terms).Fast-forward approach.-The fast-forward method [6,7,9] may be used to generate external potentials to drive the matter wave from the initial single well to a final symmetric double well. The starting point of the streamlined version in [6] is the 3D time-dependent GP equation