Stationary states, dynamical stability, and vorticity of Bose-Einstein condensates in tilted rotating harmonic traps Srivatsa B. Prasad,1, ∗ Brendan C. Mulkerin,2 and Andrew M. Martin1 1School of Physics, University of Melbourne, Melbourne, 3010, Australia 2Centre for Quantum and Optical Science, Swinburne University of Technology, Melbourne, 3122, Australia (Dated: May 13, 2020) We theoretically investigate a Bose-Einstein condensate confined by a rotating harmonic trap whose rotation axis is not aligned with any of its principal axes. The principal axes of the Thomas- Fermi density profiles of the resulting stationary solutions are found to be tilted with respect to those of the rotating trap, representing an extra degree of freedom that is associated with the existence of additional branches of stationary solutions for any given rotation axis alignment. By linearizing the time-dependent theory about the stationary states, we obtain a semi-analytical prediction of their dynamical instability at high rotation frequencies against collective modes arising from environmen- tal perturbations. Comparing the stationary states to direct simulations of the Gross-Pitaevskii equation, we predict the nucleation of quantum vortices in the dynamically unstable rotational regime. These vortex lines are aligned along the rotation axis despite the tilting of the rotating trap although the background density profile is tilted with respect to the trapping and rotation axes. I. INTRODUCTION in the condensate [11]. Alternate methods for produc- ing vortices in BECs also exist, such as stirring with a For more than 70 years the behaviour of rotating su- Gaussian laser beam [12], dragging a laser configuration perfluids have been of considerable experimental and the- through a trapped condensate (or, equivalently, a con- oretical interest to several generations of physicists. Ini- densate through a laser configuration) [13, 14], applying tially, it was recognized by Onsager [1] and Feynman [2] oscillatory perturbations to the trapping [15, 16], con- that superfluid flow is characterized by a quantized cir- densing a rotating thermal (non-condensed) atomic va- culation, which was subsequently experimentally veri- por [17], and utilizing the Kibble-Zurek mechanism by fied by Hall and Vinen [3,4]. Since then, it has been quenching a thermal vapor across the BEC critical tem- recognized that the presence of a nonzero angular mo- perature [18]. The theoretical and analytical study of the mentum in a superfluid leads to complex, nontrivial be- resulting vortices have uncovered phenomena rich in vari- haviour, spurring discoveries such as the recent realiza- ety; some examples that are relevant to scalar, nondipo- tion of negative-temperature Onsager vortex clusters in lar, single-component condensates at zero temperature two-dimensional Bose gases [5,6]. In particular, Bose- includes Kelvin waves [19{21], Abrikosov vortex lattices Einstein condensates (BECs) offer a uniquely flexible and their Tkachenko modes [22{24], quantum Hall-like platform for the study of superfluid rotation and have physics [25{28], vortex reconnections [29, 30], quantum been the focus of intense study for several years [7{9]. analogs of classical fluid instabilities [14, 31, 32], and hys- Consider a superfluid in a rotating bucket. This su- teresis [33, 34]. perfluid cannot support rigid-body rotation if the bucket By contrast, previous studies relating to tilting effects is symmetric about the rotation axis, due to the absence in rotating BECs have mainly focussed on the collec- of shear forces, and so its angular momentum manifests tive modes of vortices in response to tilting perturba- in the form of quantum vortices above a certain critical tions of the trap [35{38], while the literature concerning rotation frequency [10]. However, when the symmetry of the steady rotation of the external confinement about a the bucket about the rotation axis is broken, it is able to non-principal axis is chiefly limited to the stability of the transfer angular momentum to the superfluid even in the centre-of-mass oscillations in the rotating frame [39, 40]. absence of shear forces, and the condensate exhibits solid- Given that such tilted rotating traps may be experimen- body rotation at slow rotation frequencies and quantized tally generated in a similar manner to the excitation of vortices above a critical rotation frequency. For Bose- the tilting modes, and that roughly analogous systems Einstein condensates, in which the `bucket' is replaced such as dipolar BECs with tilted rotating dipole moments by an atomic trap generated by electromagnetic fields, have been realized experimentally [41, 42], a systematic arXiv:2004.02568v3 [cond-mat.quant-gas] 12 May 2020 angular momentum may be transferred from the trap to study of BECs confined by a tilted rotating trap is war- the condensate by modulating the applied fields such that ranted. In this paper, we analytically obtain stationary the condensate is confined by a rotating potential that is solutions of the Gross-Pitaevskii equation in the Thomas- asymmetrical about the rotation axis [7]. At sufficiently Fermi limit for a condensate subject to a range of differ- high rotation frequencies, this method induces vorticity ent tilting angles and harmonic trapping regimes. For all but the most trivial cases, the stationary solution den- sities are found to be tilted about the rotation axis by a different angle than the trap itself. One of the conse- ∗ [email protected] quences of this additional degree of freedom is the exis- 2 tence of two previously unknown branches of stationary (GPE) [7{9, 52]: solutions. These exist even when the trap is not tilted @ 1 2 2 away from the rotation axis, a result that is analogous to i VT r; t g~ iΩ r : (2) the tilted triaxial ellipsoids that are rotating-frame sta- @t 2 tionary solutions for self-gravitating irrotational classical Here we= − define∇ g ~+ 4πNa( s) l⊥+asS theS effective+ ⋅ ( strength× ∇) of a fluids [43]. Focusing on the stationary solution branch mean-field, two-body interaction, with a corresponding s- existing in the nonrotating limit, we semi-analytically lin- wave scattering length= given~ by as, and denote the time- earize the condensate's fluctuations in response to small dependent harmonic trapping potential by VT. perturbations and thus predict a dynamical instability at Previously, theoretical and experimental studies of the higher rotation frequencies, where the amplitude of one angular momentum of trapped BECs have tended to as- or more collective modes is expected to amplify expo- sume that the rotation axis of the confinement coincides nentially in time. In the regions of dynamical instabil- with one of its symmetry axes. The rotation of a har- ity we show via numerical Gross-Pitaevskii simulations monic trap about an arbitrary axis can effectively be that untilted vortices are nucleated from a tilted conden- modeled by fixing Ω Ω^z, without loss of generality, sate, despite the background condensate density still be- and assuming that the trapping potential is not symmet- ing tilted. The theoretical formalism utilized represents a ric under the transformation= z z. In the co-rotating generalization of existing theoretical methods for study- reference frame, this potential can be specified by ing BECs in asymmetric, untilted rotating traps [44{49] → − and are readily conducive to experimental investigation 1 2 2 VT r 1 " x cos θ z sin θ 1 " y along the lines of previous studies that have probed the 2 untilted regime [11, 50, 51]. 1 2 ( ) = γ(2 x−sin)(θ z cos+θ ; ) + ( + ) (3) This paper is structured as follows. Section II defines 2 the concept of a tilted rotating trap and introduces the where " + 1; 1 (and γ − R. This) external potential is relevant coordinate reference frames, while Sec. III dis- equivalent to cusses the methodology for solving for the vortex-free sta- ∈ (− ) ∈ tionary solutions in the Thomas-Fermi limit. In Sec. IV, 1 V R 1 " X2 1 " Y 2 γ2Z2 ; (4) we examine the features of these stationary solutions for T 2 two distinct trapping regimes, and in Sec.V, the time- via a rotation( ) = of the( co-rotating− ) + ( coordinates+ ) + as given by dependent theory is linearized in order to characterize the dynamical stability of the vortex-free stationary so- X cos θ 0 sin θ x lutions during a quasi-adiabatic rampup of the trapping Y 0 1 0 y : (5) rotation frequency. Finally, Sec.VI contains a discussion ⎛Z ⎞ ⎛ sin θ 0 cos θ⎞ ⎛z⎞ of the outcomes of a series of numerical simulations of ⎜ ⎟ = ⎜ ⎟ ⎜ ⎟ such rampups, where the dynamical route to vortex nu- By inspection,⎝ Eq.⎠ (3⎝)− is equivalent to⎠ ⎝ Eq.⎠ (4) when, for cleation in a vorticity-free condensate in a tilted rotating integer n, the tilting angle obeys θ nπ. A similar equiv- trap is demonstrated. alence, albeit with modified values of γ and ", holds when n takes on half-integer values. = To simulate the stationary state and dynamics of a II. THE TILTED, ROTATING HARMONIC BEC in this trap via numerical methods, it is sufficient TRAP to use Eqs. (2) and (3). However, the vorticity-free sta- tionary solutions of Eq. (3), and their linear response to environmental perturbations, are well-described in the In order to describe a dilute, scalar BEC, at zero tem- θ 0 limit by purely semi-analytical methods [44, 45]. perature in a rotating, tilted harmonic trap, we utilize To utilise these methods for an arbitrary value of θ it is the Gross-Pitaevskii equation for the condensate order necessary= to set up a hydrodynamic formalism. This in- parameter, . We assume that N condensed bosons, volves the definition of the condensate's density, n, phase, each with a mass m, are confined in the trap and that S, and superfluid velocity, v, via the relations [7]: the root mean squared harmonic trapping frequency in the x-y plane is given by !⊥.