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.

We investigate the evolution of mass segregation in initially sub-structured young embedded star clusters with two different background potentials mimicking the gas. Our clusters are initially in virial or sub-virial global states and have different initial distributions for the most massive stars: randomly placed, initially mass segregated or even inverse segregation. By means of N-body simulation we follow their evolution for 5 Myr. We measure the mass segregation using the minimum spanning tree method Lambda_MSR and an equivalent restricted method. Despite this variety of different initial conditions, we find that our stellar distributions almost always settle very fast into a mass segregated and more spherical configuration, suggesting that once we see a spherical or nearly spherical embedded star cluster, we can be sure it is mass segregated no matter what the real initial conditions were. We, furthermore, report under which circumstances this process can be more rapid or delayed, respectively.

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).

Hercules is a dwarf spheroidal satellite of the Milky Way, found at a distance of about 138 kpc, and showing evidence of tidal disruption. It is very elongated and exhibits a velocity gradient of 16 +/- 3 km/s/kpc. Using this data a possible orbit of Hercules has previously been deduced in the literature. In this study we make use of a novel approach to find a best fit model that follows the published orbit. Instead of using trial and error, we use a systematic approach in order to find a model that fits multiple observables simultaneously. As such, we investigate a much wider parameter range of initial conditions and ensure we have found the best match possible. Using a dark matter free progenitor that undergoes tidal disruption, our best-fit model can simultaneously match the observed luminosity, central surface brightness, effective radius, velocity dispersion, and velocity gradient of Hercules. However, we find it is impossible to reproduce the observed elongation and the position angle of Hercules at the same time in our models. This failure persists even when we vary the duration of the simulation significantly, and consider a more cuspy density distribution for the progenitor. We discuss how this suggests that the published orbit of Hercules is very likely to be incorrect.