Self-consistent models of triaxial galaxies in MOND gravity

The Bekenstein- Milgrom gravity theory with a modified Poisson equation is tested here for the existence of triaxial equilibrium solutions. Using the nonnegative least- square method, we show that self- consistent triaxial galaxies exist for baryonic models with a mild density cusp, rho similar to Sigma/r. Self- consistency is achieved for a wide range of central concentrations, Sigma similar to 10-1000 M-circle dot pc(-2), representing low to high surface brightness galaxies. Our results demonstrate for the first time that the orbit superposition technique is fruitful for constructing galaxy models beyond Newtonian gravity, and triaxial cuspy galaxies might exist without the help of cold dark matter.