Abstract:We study the maximally supersymmetric plane wave matrix model (the BMN model) at finite temperature, T , and locate the high temperature phase boundary in the (µ, T ) plane, where µ is the mass parameter. We find the first transition, as the system is cooled from high temperatures, is from an approximately SO(9) symmetric phase to one where three matrices expand to form fuzzy spheres. For µ > 3.0 there is a second distinct transition at a lower temperature. The two transitions approach one another at smaller µ and merge in the vicinity of µ = 3.0. The resulting single transition curve then approaches the gauge/gravity prediction as µ is further decreased. We find a rough estimate of the transition, for all µ, is given by a Padé resummation of the large-µ, three-loop perturbative predictions. We find evidence that the transition at small µ is to an M5-brane phase of the theory.