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.

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

We predict a new kind of instability in a Bose-Einstein condensate composed of dipolar particles. Namely, a comparatively weak dipole moment can produce a large, negative two-body scattering length that can collapse the Bose-Einstein condensate. To verify this effect, we validate mean-field solutions to this problem using exact, diffusion Monte Carlo methods. We show that the diffusion Monte Carlo energies are reproduced accurately within a mean-field framework if the variation of the s-wave scattering length with the dipole strength is accounted for properly.