We present a general method for obtaining the exact static solutions and collective excitation frequencies of a trapped Bose-Einstein condensate (BEC) with dipolar atomic interactions in the Thomas-Fermi regime. The method incorporates analytic expressions for the dipolar potential of an arbitrary polynomial density profile, thereby reducing the problem of handling nonlocal dipolar interactions to the solution of algebraic equations. We comprehensively map out the static solutions and excitation modes, including non-cylindrically-symmetric traps, and also the case of negative scattering length where dipolar interactions stabilize an otherwise unstable condensate. The dynamical stability of the excitation modes gives insight into the onset of collapse of a dipolar BEC. We find that global collapse is consistently mediated by an anisotropic quadrupolar collective mode, although there are two trapping regimes in which the BEC is stable against quadrupole fluctuations even as the ratio of the dipolar to $s$-wave interactions becomes infinite. Motivated by the possibility of a fragmented condensate in a dipolar Bose gas due to the partially attractive interactions, we pay special attention to the scissors modes, which can provide a signature of superfluidity, and identify a long-range restoring force which is peculiar to dipolar systems. As part of the supporting material for this paper we provide the computer program used to make the calculations, including a graphical user interface.

• Received 28 February 2010

DOI:https://doi.org/10.1103/PhysRevA.82.033612