We study Bose-Einstein condensates with purely dipolar interactions in oblate traps. We find that the condensate always becomes unstable to collapse when the number of particles is sufficiently large. We analyze the instability, and find that it is the trapped-gas analogue of the "roton-maxon" instability previously reported for a gas that is unconfined in 2D. In addition, we find that under certain circumstances the condensate wave function attains a biconcave shape, with its maximum density away from the center of the gas. These biconcave condensates become unstable due to azimuthal excitation--an angular roton.

We report on measurements of the excitation spectrum of a strongly interacting Bose-Einstein condensate. A magnetic-field Feshbach resonance is used to tune atom-atom interactions in the condensate and to reach a regime where quantum depletion and beyond mean-field corrections to the condensate chemical potential are significant. We use two-photon Bragg spectroscopy to probe the condensate excitation spectrum; our results demonstrate the onset of beyond mean-field effects in a gaseous Bose-Einstein condensate.

The calculation of properties of Bose-Einstein condensates with dipolar interactions has proven a computationally intensive problem due to the long range nature of the interactions, limiting the scope of applications. In particular, the lowest lying Bogoliubov excitations in three dimensional harmonic trap with cylindrical symmetry were so far computed in an indirect way, by Fourier analysis of time dependent perturbations, or by approximate variational methods. We have developed a very fast and accurate numerical algorithm based on the Hankel transform for calculating properties of dipolar Bose-Einstein condensates in cylindrically symmetric traps. As an application, we are able to compute many excitation modes by directly solving the Bogoliubov-De Gennes equations. We explore the behavior of the excited modes in different trap geometries. We use these results to calculate the quantum depletion of the condensate by a combination of a computation of the exact modes and the use of a local density approximation