
Fast transport of Bose-Einstein condensates in anharmonic traps
Jing Li,1, ∗Xi Chen,2, 3 and Andreas Ruschhaupt1
1Department of Physics, University College Cork, Cork, Ireland, T12 H6T1
2Department of Physical Chemistry, University of the Basque Country UPV/EHU, Apartado 644, 48080 Bilbao, Spain
3EHU Quantum Center, University of the Basque Country UPV/EHU, 48940 Leioa, Spain
(Dated: October 18, 2022)
We present a method to transport Bose-Einstein condensates (BECs) in anharmonic traps and in the presence
of atom-atom interactions in short times without residual excitation. Using a combination of a variational ap-
proach and inverse engineering methods, we derive a set of Ermakov-like equations that take into account the
coupling between the center of mass motion and the breathing mode. By an appropriate inverse engineering
strategy of those equations, we then design the trap trajectory to achieve the desired boundary conditions. Nu-
merical examples for cubic or quartic anharmonicities are provided for fast and high-fidelity transport of BECs.
Potential applications are atom interferometry and quantum information processing.
I. INTRODUCTION
The accurate manipulation of ultracold atoms is a key pre-requisite to implement quantum technologies within atomic, molec-
ular and optical science [1]. In particular, the transport of individual atoms and of thermal or Bose-condensed clouds using
moving traps has been demonstrated in many experiments [2–12] for different goals in quantum information processing and
metrology. In all quantum technologies, preserving quantum coherence and achieving high final fidelities in short times is of
crucial importance. One possibility is called shortcuts to adiabaticity (STA) [13,14] which provides a toolbox to control both
the internal and external degrees of freedom of a quantum system in faster-than adiabatic times.
Various shortcuts to adiabatic transport have been proposed: Lewis-Riesenfeld invariant-based inverse engineering [15–19],
enhanced STA scheme [20–22], the Fourier optimization [23], fast-forward scaling method [24,25], and the counter-diabatic
driving [26] have been theoretically put forward, and experimentally demonstrated for various systems [7,10,11,27]. The
possibility to operate with short times not only reduces the sensitivity to low frequency noise, but also allows for improved
measurement statistics in the total time available for the experiment.
Different approaches for transporting particles have been implemented. Neutral atoms have been transported as Bose-Einstein
condensates (BECs) [2], thermal atomic clouds [28], or individually [5], using magnetic or optical traps. The commonly used
traps for ultracold atoms based on electromagnetic fields are never perfectly harmonic. The weak cubic anharmonicity plays
a role when a BEC is transported perpendicular to the atom chip surface [29]. The quartic anharmonicity is significant when
approximating the potential of an optical tweezers for transport [7,16]. Thus cancelling the anharmonic contributions of the
trapping potential is vital for useful control schemes and is already a difficult technical challenge for a static trap [30].
Anharmonicities can have an important impact on the dynamics as observed in atom cooling [31], collective modes [32], or
wave packet dynamics [33]. In most cases, the anharmonic traps are considered as a perturbation of a harmonic one. Perturbation
theory has been used to design shortcut protocols for expansion/compression [34] and transport [35]. Of course the results are
limited by the premises of perturbation theory, i.e., by small anharmonicities. Considering a non-perturbative scenario is thus of
much interest.
In this paper, we propose to inverse engineer rapid and robust transport of an interacting Bose-Einstein condensate (BEC) in
anharmonic traps using a variational approach. The method relies on a variational formulation of the dynamics to derive a set of
coupled Ermakov-like and Newton-like equations, from which the trap trajectory is inferred interpolating between the desired
boundary conditions. In Sec. II, we explain the variational formalism. In Sec. III, we work out the explicit solutions for quartic
and cubic anharmonicities of the confining potential, and illustrate the efficiency of the method with various numerical examples.
In section IV, we will discuss the results.
II. MODEL, HAMILTONIAN AND METHOD
For a cigar-shaped trap with strong transverse confinement, e.g. ω⊥>> ω, it is appropriate to consider a 1D dimensionless
formula by freezing the transverse dynamics to the respective ground state and integrating over the transverse variables[36]. The
∗jli@ucc.ie
arXiv:2210.03788v1 [cond-mat.quant-gas] 7 Oct 2022