Bose-Einstein condensates and the thin-shell limit in anisotropic bubble traps 2
for confining atomic gases into shells in order to study its properties, but some new alternatives such
as gravity compensation mechanisms are arising [23,24]. Also, an interesting substitute to the usual
procedure of radio-frequency dressing was proposed for dual-species atomic mixtures [25], which led to
the creation of a BEC on Earth’s gravity [26].
Confinements in three-dimensional shell shaped condensates inspired some theoretical works worth
mentioning here. For instance, in [27], the authors apply analytical methods to investigate the ground
state wave function of a BEC and its collective modes. In [28], both analytical and numerical techniques
are used in order to obtain expansion properties. The interesting paper [29] employs thermodynamic
arguments in order to survey the formation of clusters. The thermodynamics of a BEC on a spherical
shell is analyzed, including the critical temperature, in Refs. [30–32]. The topological hollowing
transition from a full sphere to a thin-shell was studied in [33] and [34], thus finding some universal
properties. The ground state and collective excitations of a dipolar BEC was considered in [35]. The
general physical relevance of cold atoms on curved manifolds is also addressed in [36]. The contribution
of [37] plays with the idea that the external potential is equivalent to the harmonic trap for a large
radius thin-shell. Universal scaling relations are found for topological superfluid transitions in bubble
traps in [38]. And non-Hermitian phase transitions are meticulously described in [39]. Although it
is not the focus of this work, it is also worth citing vortex studies on spherical-like surfaces with no
holes [40–47].
In this paper, we work out an explicit relation between the confinement of the particles in the
thin-shell limit and the geometrical distortion of a bubble trap for a family of confinement potentials.
They are chosen in such a way that the current experiments are included as special cases. All potentials
in this family turn out to exhibit the same angular dependency for the confinement strength. Section 2
describes the mathematical background by defining the concepts in which our theory is developed. In a
second step we calculate in section 3 the harmonic radial frequency in the Gaussian Normal Coordinate
System (GNCS) and expose its resulting angular dependency. Furthermore, the definition of the thin-
shell limit is discussed in a more rigorous way by elucidating how it depends on the geometrical
distortion of the bubble trap. The general Gross-Pitaevskii Hamiltonian of the system is deduced
using a perturbative approach near the thin-shell limit in section 4. In section 5, special topics of
the spherical shell and the Thomas-Fermi approximation are considered. Finally, in section 6, the
corresponding effective Gross-Pitaevskii equation is perturbatively solved for small distortions from a
spherical bubble trap in order to determine both the ground-state wave function and the oscillation
frequencies.
2. Gaussian Normal Coordinate System
In this section, some preliminary concepts are defined and explained in order to establish the
mathematical background in which our theory is developed. One of such main concepts concerns the
manifold [48, 49] considered here. In this work, we study 2D surfaces embedded into a 3D Euclidean
space. More specifically, ellipsoidal surfaces [50–52] are considered since they correspond to the bubble
trap potentials in BEC experiments. Therefore, our manifolds are compact, smooth and differentiable
everywhere. In order to describe the 3D region around the 2D manifold we choose as a suitable
coordinate system, the so-called Gaussian Normal Coordinate System (GNCS) [53,54]. Further features
and particularities on its application can be found in [27].
It is always possible to describe the region around smooth manifolds with the aid of a GNCS.
The main idea is to consider two coordinates x1and x2over the 2D manifold M, also called tangent
coordinates, describing arbitrary points in the manifold. Thus, any point pof this manifold Mis
portrayed by the position vector −→
p(x1, x2). Any point qin such a vicinity of Mcan be represented
by a coordinate x0referred to as the orthogonal coordinate, and a normal unit vector ˆnat the point p
through the following equation
−→
q(x0, x1, x2) = −→
p(x1, x2) + x0ˆn(x1, x2).(1)