Abstract

Attractive Bose-Einstein condensates are investigated with numerical continuation methods capturing stationary solutions of the Gross-Pitaevskii equation. The branches of stable (elliptic) and unstable (hyperbolic) solutions are found to meet at a critical particle number through a generic Hamiltonian saddle node bifurcation. The condensate decay rates corresponding to macroscopic quantum tunneling, two and three body inelastic collisions, and thermally induced collapse are computed from the exact numerical solutions. These rates show experimentally significant differences with previously published rates. Universal scaling laws stemming from the bifurcation are derived.