Self-consistent models of cuspy triaxial galaxies with dark matter halos

We have constructed realistic, self-consistent models of triaxial elliptical galaxies embedded in triaxial dark matter halos. We examined three different models for the shape of the dark matter halo: (1) with the same axis ratios as the luminous matter ( 0.7:0.86:1), (2) with a more prolate shape (0.5:0.66:1), and (3) with a more oblate shape (0.7:0.93:1). Self-consistent solutions by means of the standard orbital superposition technique introduced by Schwarzschild were found in each of the three cases. The equilibrium velocity distribution is reproduced by a Lorentzian function better than by a Gaussian. Chaotic orbits were found to be important in all of the models, and their presence was shown to imply a possible slow evolution of the shapes of the halos. Our results demonstrate for the first time that triaxial dark matter halos can coexist with triaxial galaxies.