Abstract

The Andromeda galaxy (M31) contains a box/peanut bulge (BPB) entangled with a classical bulge (CB) requiring a triaxial modelling to determine the dynamics, stellar, and dark matter mass. We construct made-to-measure models fitting new VIRUS-W IFU bulge stellar kinematic observations, the IRAC-3.6 mu m photometry, and the disc's HI rotation curve. We explore the parameter space for the 3.6 mu m mass-to-light ratio ((sic)(3.6)), the bar pattern speed (Omega(p)), and the dark matter mass in the composite bulge (M-DM(B)) within 3.2 kpc. Considering Einasto dark matter profiles, we find the best models for (sic)(3.6) = 0.72 +/- 0.02 M circle dot L circle dot-1, M-DM(B) = 1.2(-0.4)(+0.2) x 10(10) M-circle dot, and Omega(p) = 40 +/- 5 km s(-1) kpc(-1). These models have a dynamical bulge mass of M-dyn(B) = 4.25(-0.29)(+0.10) x 10(10) M-circle dot including a stellar mass of M-star(B) = 3.09(-0.12)(+ 0.10) x 10(10) M-circle dot (73 per cent), of which the CB has M-star(CB) = 1.18(-0.07)(+0.06)x 10(10) M-circle dot (28 per cent) and the BPB M-star(BPB) = 1.91 +/- 0.06 x 10(10) M-circle dot (45 per cent). We also explore models with NFW haloes finding that, while the Einasto models better fit the stellar kinematics, the obtained parameters agree within the errors. The M-DM(B) values agree with adiabatically contracted cosmological NFW haloes with M31's virial mass and radius. The best model has two bulge components with completely different kinematics that only together successfully reproduce the observations (mu(3.6), upsilon(los), sigma(los), h3, h4). The modelling includes dust absorption which reproduces the observed kinematic asymmetries. Our results provide new constraints for the early formation of M31 given the lower mass found for the classical bulge and the shallow dark matter profile, as well as the secular evolution of M31 implied by the bar and its resonant interactions with the classical bulge, stellar halo and disc.