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 $(\Upsilon_{3.6})$, the bar pattern speed ($\Omega_p$), and the dark matter mass in the composite bulge ($M^B_{DM}$) within 3.2kpc. Considering Einasto dark matter profiles, we find the best models for $\Upsilon_{3.6}=0.72\pm0.02\,M_\odot/L_\odot$, $M^B_{DM}=1.2^{+0.2}_{-0.4}\times10^{10}M_\odot$ and $\Omega_p=40\pm5\,km/s/kpc$. These models have a dynamical bulge mass of $M_{dyn}^B=4.25^{+0.10}_{-0.29}\times10^{10}M_{\odot}$ including a stellar mass of $M^B=3.09^{+0.10}_{-0.12}\times10^{10}M_\odot$(73%), of which the CB has $M^{CB}=1.18^{+0.06}_{-0.07}\times10^{10}M_\odot$(28%) and the BPB $M^{BPB}=1.91\pm0.06\times10^{10}M_\odot$(45%). 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^B_{DM}$ 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.

We continue the analysis of the dataset of our spectroscopic observation campaign of M31, by deriving simple stellar population properties (age metallicity and alpha-elements overabundance) from the measurement of Lick/IDS absorption line indices. We describe their two-dimensional maps taking into account the dust distribution in M31. 80\% of the values of our age measurements are larger than 10 Gyr. The central 100 arcsec of M31 are dominated by the stars of the classical bulge of M31. They are old (11-13 Gyr), metal-rich (as high as [Z/H]~0.35 dex) at the center with a negative gradient outwards and enhanced in alpha-elements ([alpha/Fe]~ 0.28+- 0.01 dex). The bar stands out in the metallicity map, where an almost solar value of [Z/H] (~0.02+-0.01 dex) with no gradient is observed along the bar position angle (55.7 deg) out to 600 arcsec from the center. In contrast, no signature of the bar is seen in the age and [alpha/Fe] maps, that are approximately axisymmetric, delivering a mean age and overabundance for the bar and the boxy-peanut bulge of 10-13 Gyr and 0.25-0.27 dex, respectively. The boxy/peanut-bulge has almost solar metallicity (-0.04+- 0.01 dex). The mass-to-light ratio of the three components is approximately constant at M/LV ~ 4.4-4.7 Msol/Lsol. The disk component at larger distances is made of a mixture of stars, as young as 3-4 Gyr, with solar metallicity and smaller M/LV (~3+-0.1 Msol/Lsol). We propose a two-phase formation scenario for the inner region of M31, where most of the stars of the classical bulge come into place together with a proto-disk, where a bar develops and quickly transforms it into a boxy-peanut bulge. Star formation continues in the bulge region, producing stars younger than 10 Gyr, in particular along the bar, enhancing its metallicity. The disk component appears to build up on longer time-scales.

As the nearest large spiral galaxy, M31 provides a unique opportunity to learn about the structure and evolutionary history of this galaxy type in great detail. Among the many observing programs aimed at M31 are microlensing studies, which require good three-dimensional models of the stellar mass distribution. Possible non-axisymmetric structures like a bar need to be taken into account. Due to M31's high inclination, the bar is difficult to detect in photometry alone. Therefore, detailed kinematic measurements are needed to constrain the possible existence and position of a bar in M31. We obtained $\approx$ 220 separate fields with the optical IFU spectrograph VIRUS-W, covering the whole bulge region of M31 and parts of the disk. We derive stellar line-of-sight velocity distributions from the stellar absorption lines, as well as velocity distributions and line fluxes of the emission lines H$\beta$, [OIII] and [NI]. Our data supersede any previous study in terms of spacial coverage and spectral resolution. We find several features that are indicative of a bar in the kinematics of the stars, we see intermediate plateaus in the velocity and the velocity dispersion, and correlation between the higher moment $h3$ and the velocity. The gas kinematics is highly irregular, but is consistent with non-triaxial streaming motions caused by a bar. The morphology of the gas shows a spiral pattern, with seemingly lower inclination than the stellar disk. We also look at the ionization mechanisms of the gas, which happens mostly through shocks and not through starbursts.

Andromeda is our nearest neighbouring disk galaxy and a prime target for detailed modelling of the evolutionary processes that shape galaxies. We analyse the nature of M31's triaxial bulge with an extensive set of N-body models, which include Box/Peanut (B/P) bulges as well as initial classical bulges (ICBs). Comparing with IRAC 3.6$\mu m$ data, only one model matches simultaneously all the morphological properties of M31's bulge, and requires an ICB and a B/P bulge with 1/3 and 2/3 of the total bulge mass respectively. We find that our pure B/P bulge models do not show concentrations high enough to match the S\'ersic index ($n$) and the effective radius of M31's bulge. Instead, the best model requires an ICB component with mass $M^{\rm ICB}=1.1\times10^{10}{\rm M_{\odot}}$ and three-dimensional half-mass radius $r_{\rm half}^{\rm ICB}$=0.53 kpc (140 arcsec). The B/P bulge component has a mass of $M^{\rm B/P}=2.2\times10^{10}{\rm M_{\odot}}$ and a half-mass radius of $r_{\rm half}^{\rm B/P}$=1.3 kpc (340 arcsec). The model's B/P bulge extends to $r^{\rm B/P}$=3.2 kpc (840 arcsec) in the plane of the disk, as does M31's bulge. In this composite bulge model, the ICB component explains the velocity dispersion drop observed in the centre within $R<$190 pc (50 arcsec), while the B/P bulge component reproduces the observed rapid rotation and the kinematic twist of the observed zero velocity line. This model's pattern speed is $\Omega_p$=38 km/s/kpc, placing corotation at $r_{\rm cor}$=5.8 kpc (1500 arcsec). The outer Lindblad resonance (OLR) is then at $r_{\rm OLR}$=10.4kpc, near the 10kpc-ring of M31, suggesting that this structure may be related to the bar's OLR. By comparison with an earlier snapshot, we estimate that M31's thin bar extends to $r_{\rm bar}^{\rm thin}\sim$4.0 kpc (1000 arcsec) in the disk plane, and in projection extends to $R_{\rm bar}^{\rm thin}\sim$2.3 kpc (600 arcsec).

We investigate the correlations between the black hole mass $M_{BH}$, the velocity dispersion $\sigma$, the bulge mass $M_{Bu}$, the bulge average spherical density $\rho_h$ and its spherical half mass radius $r_h$, constructing a database of 97 galaxies (31 core ellipticals, 17 power-law ellipticals, 30 classical bulges, 19 pseudo bulges) by joining 72 galaxies from the literature to 25 galaxies observed during our recent SINFONI black hole survey. For the first time we discuss the full error covariance matrix. We analyse the well known $M_{BH}-\sigma$ and $M_{BH}-M_{Bu}$ relations and establish the existence of statistically significant correlations between $M_{Bu}$ and $r_h$ and anti-correlations between $M_{Bu}$ and $\rho_h$. We establish five significant bivariate correlations ($M_{BH}-\sigma-\rho_h$, $M_{BH}-\sigma-r_h$, $M_{BH}-M_{Bu}-\sigma$, $M_{BH}-M_{Bu}-\rho_h$, $M_{BH}-M_{Bu}-r_h$) that predict $M_{BH}$ of 77 core and power-law ellipticals and classical bulges with measured and intrinsic scatter as small as $\approx 0.36$ dex and $\approx 0.33$ dex respectively, or 0.26 dex when the subsample of 45 galaxies defined by Kormendy and Ho (2013) is considered. In contrast, pseudo bulges have systematically lower $M_{BH}$, but approach the predictions of all the above relations at spherical densities $\rho_h\ge 10^{10} M_\odot/kpc^3$ or scale lengths $r_h\le 1$ kpc. These findings fit in a scenario of co-evolution of BH and classical-bulge masses, where core ellipticals are the product of dry mergers of power-law bulges and power-law Es and bulges the result of (early) gas-rich mergers and of disk galaxies. In contrast, the (secular) growth of BHs is decoupled from the growth of their pseudo bulge hosts, except when (gas) densities are high enough to trigger the feedback mechanism responsible for the existence of the correlations between $M_{BH}$ and galaxy structural parameters.

Most Milky Way globular clusters (GCs) exhibit measurable flattening, even if on a very low level. Both cluster rotation and tidal fields are thought to cause this flattening. Nevertheless, rotation has only been confirmed in a handful of GCs, based mostly on individual radial velocities at large radii. We are conducting a survey of the central kinematics of Galactic GCs using the new Integral Field Unit instrument VIRUS-W. We detect rotation in all 11 GCs that we have observed so far, rendering it likely that a large majority of the Milky Way GCs rotate. We use published catalogs of the ACS survey of GCs to derive central ellipticities and position angles. We show that in all cases where the central ellipticity permits an accurate measurement of the position angle, those angles are in excellent agreement with the kinematic position angles that we derive from the VIRUS-W velocity fields. We find an unexpected tight correlation between central rotation and outer ellipticity, indicating that rotation drives flattening for the objects in our sample. We also find a tight correlation between central rotation and published values for the central velocity dispersion, most likely due to rotation impacting the old dispersion measurements.

We present the first results of an analysis of the properties of the molecular gas in the nuclear regions (r < 300 pc) of a sample of six nearby galaxies, based on new high spatial resolution observations obtained in the K-band with the near-infrared integral field spectrograph SINFONI at the Very Large Telescope. We derive two-dimensional distributions of the warm molecular and ionized gas from the H2, Br_gamma and HeI emission lines present in the spectra of the galaxies. We find a range of morphologies, including bar- and ring-like distributions and either centrally peaked or off-centre emission. The morphologies of the molecular and the ionized gas are not necessarily coincident. The observed emission-line ratios point towards thermal processes as the principal mechanism responsible for the H2 excitation in the nuclear and circumnuclear regions of the galaxies, independently of the presence of an active nucleus. We find that a rescaling of the H2 2.12 microns emission-line luminosity by a factor beta~1200 gives a good estimate (within a factor of 2) of the total (cold) molecular gas mass. The galaxies of the sample contain large quantities of molecular gas in their centres, with total masses in the ~ 105 - 108 Msol range. Never the less, these masses correspond to less than 3 per cent of the stellar masses derived for the galaxies in these regions, indicating that the presence of gas should not affect black hole mass estimates based on the dynamical modelling of the stars. The high-spatial resolution provided by the SINFONI data allowed us to resolve a circumnuclear ring (with a radius of ~270 pc) in the galaxy NGC 4536. The measured values of the Br_gamma equivalent width and the HeI/Br_gamma emission-line ratio suggests that bursts of star formation occurred throughout this ring as recently as 6.5 Myr ago.