To understand the history and formation mechanisms of galaxies it is crucial to determine their current multidimensional structure. Here we focus on stellar population properties, such as metallicity and [$\alpha$/Fe] enhancement. We devise a new technique to recover the distribution of these parameters using spatially resolved, line-of-sight averaged data. Our chemodynamical method is based on the made-to-measure (M2M) framework and results in an $N$-body model for the abundance distribution. We test our method on a mock dataset and find that the radial and azimuthal profiles are well-recovered, however only the overall shape of the vertical profile matches the true profile. We apply our procedure to spatially resolved maps of mean [Z/H] and [$\alpha$/Fe] for the Andromeda Galaxy, using an earlier barred dynamical model of M31. We found that the metallicity is enhanced along the bar, with possible maxima at the ansae. In the edge-on view the [Z/H] distribution has an X shape due to the boxy/peanut bulge; the average vertical metallicity gradient is equal to $-0.132\pm0.006$ dex/kpc. We identify a metallicity-enhanced ring around the bar which also has relatively lower [$\alpha$/Fe]. The highest [$\alpha$/Fe] is found in the centre, due to the classical bulge. Away from the centre, the $\alpha$-overabundance in the bar region increases with height, which could be an indication of a thick disc. We argue that the galaxy assembly resulted in a sharp peak of metallicity in the central few hundred parsecs and a more gentle negative gradient in the remaining disc, but no [$\alpha$/Fe] gradient. The formation of the bar lead to the re-arrangement of the [Z/H] distribution, causing a flat gradient along the bar. Subsequent star formation close to the bar ends may have produced the metallicity enhancements at the ansae and the [Z/H] enhanced lower-$\alpha$ ring.

The age-velocity dispersion relation is an important tool to understand the evolution of the disc of the Andromeda galaxy (M31) in comparison with the Milky Way. We use Planetary Nebulae (PNe) to obtain the age-velocity dispersion relation in different radial bins of the M31 disc. We separate the observed PNe sample based on their extinction values into two distinct age populations. The observed velocities of our high- and low-extinction PNe, which correspond to higher and lower mass progenitors respectively, are fitted in de-projected elliptical bins to obtain their rotational velocities, $V_{\phi}$, and corresponding dispersions, $\rm\sigma_{\phi}$. We assign ages to the two PNe populations by comparing central-star properties of an archival sub-sample of PNe, having models fitted to their observed spectral features, to stellar evolution tracks. For the high- and low-extinction PNe, we find ages of $\sim2.5$ Gyr and $\sim4.5$ Gyr respectively, with distinct kinematics beyond a deprojected radius R$\rm_{GC}= 14$ kpc. At R$\rm_{GC}$=17--20 kpc, which is the equivalent distance in disc scale lengths of the Sun in the Milky Way disc, we obtain $\rm\sigma_{\phi,~2.5~Gyr}= 61\pm 14$ km s$^{-1}$ and $\rm\sigma_{\phi,~4.5~Gyr}= 101\pm 13$ km s$^{-1}$. The age-velocity dispersion relation for the M31 disc is obtained in two radial bins, R$\rm_{GC}$=14--17 and 17--20 kpc. The high- and low-extinction PNe are associated with the young thin and old thicker disc of M31 respectively, whose velocity dispersion values increase with age. These values are almost twice and thrice that of the Milky Way disc stellar population of corresponding ages. From comparison with simulations of merging galaxies, we find that the age-velocity dispersion relation in the M31 disc measured using PNe is indicative of a single major merger that occurred 2.5 -- 4.5 Gyr ago with an estimated merger mass ratio $\approx$ 1:5.

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.

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.

We attempt to find a progenitor for the ultra-faint object Segue 1 under the assumption that it formed as a dark matter free star cluster in the past. We look for orbits, using the elongation of Segue 1 on the sky as a tracer of its path. Those orbits are followed backwards in time to find the starting points of our N-body simulations. The successful orbit, with which we can reproduce Segue 1 has a proper motion of mu_alpha = -0.19 mas/yr and mu_delta = -1.9 mas/yr, placing Segue 1 near its apo-galacticon today. Our best fitting model has an initial mass of 6224 Msun and an initial scale-length of 5.75 pc.

We examine the effects of gas-expulsion on initially substructured distributions of stars. We perform N-body simulations of the evolution of these distributions in a static background potential to mimic the gas. We remove the static potential instantaneously to model gas-expulsion. We find that the exact dynamical state of the cluster plays a very strong role in affecting a cluster's survival, especially at early times: they may be entirely destroyed or only weakly affected. We show that knowing both detailed dynamics and relative star-gas distributions can provide a good estimate of the post-gas expulsion state of the cluster, but even knowing these is not an absolute way of determining the survival or otherwise of the cluster.

Dwarf spheroidal (dSph) galaxies are considered the basic building blocks of the galaxy formation process in the LCDM (Lambda Cold Dark Matter) hierarchical cosmological model. These galaxies are believed to be the most dark matter (DM) dominated systems known, have the lowest stellar content, and are poor in gas. Many theories attempt to explain the formation of dSph galaxies resorting to the fact that these galaxies are mainly found orbiting large galaxies or invoking other mechanisms of interactions. Here we show the full set of simulation as an extension of our fiducial model, where we study the formation of classical dSph galaxies in isolation by dissolving star clusters within the DM halo of the dwarf galaxy. In our parameter survey we adopt cored and cusped DM halo profiles and consider different numbers of dissolving star clusters. We investigate the dependency of observable quantities with different masses and scale-lengths of the DM halo and different star formation efficiencies (SFE). We find that our proposed scenario explains many features of the classical dSph galaxies of the Milky Way, like their morphology and their dynamics. We see trends how the surface brightness and the scale-length of the luminous component vary with the parameters of our simulations. We also identify how irregularities in their shape, i.e. clumpiness and ellipticity vary in our simulations. In velocity space, we identify the parameters leading to flat velocity dispersions curves. We recognize kinematically cold substructures in velocity space, named fossil remnants and stemming from our unique initial conditions, which alter the expected results. These streaming motions are considered as a key feature for future observation with high resolution to validate our scenario.

Recent deep photometry of the dwarf spheroidal Ursa Major II's morphology, and spectroscopy of individual stars, have provided a number of new constraints on its properties. With a velocity dispersion $\sim$6 km s$^{-1}$, and under the assumption that the galaxy is virialised, the mass-to-light ratio is found to be approaching $\sim$2000 - apparently heavily dark matter dominated. Using N-Body simulations, we demonstrate that the observed luminosity, ellipticity, irregular morphology, velocity gradient, and the velocity dispersion can be well reproduced through processes associated with tidal mass loss, and in the absence of dark matter. These results highlight the considerable uncertainty that exists in measurements of the dark matter content of Ursa Major II. The dynamics of the inner tidal tails, and tidal stream, causes the observed velocity dispersion of stars to be boosted to values of $>$5 km s$^{-1}$ ($>$20 km s$^{-1}$ at times). This effect is responsible for raising the velocity dispersion of our model to the observed values in UMaII. We test an iterative rejection technique for removing unbound stars from samples of UMaII stars whose positions on the sky, and line-of-sight velocities, are provided. We find this technique is very effective at providing an accurate bound mass from this information, and only fails when the galaxy has a bound mass less than 10$%$ of its initial mass. However when $<2%$ mass remains bound, mass overestimation by $>$3 orders of magnitude are seen. Additionally we find that mass measurements are sensitive to measurement uncertainty in line-of-sight velocities. Measurement uncertainties of 1-4 km s$^{-1}$ result in mass overestimates by a factor of $\sim$1.3-5.7.

The last few years have seen the discovery of many faint and ultra-faint dwarf spheroidal galaxies around the Milky Way. Among these is a pair of satellites called Leo IV and Leo V. This pair is found at large distances from the Milky Way (154 and 175 kpc respectively). The rather small difference in radial distance, and the fact that they also show a close projected distance on the sky, has led to the idea that we might be seeing a new pair of bound galaxies - like the Magellanic Clouds. In this paper we investigate this speculation by means of a simple integration code (confirming the results with full N-body simulations). As the luminous mass of both faint dwarfs is far too low to allow them to be bound, we simulate the pair assuming extended dark matter haloes. Our results show that the minimum dark matter mass required for the pair to be bound is rather high - ranging from 1.6 x 10^10 Msun to 5.4 x 10^10 Msun (within the virial radii). Computing the mass of dark matter within a commonly adopted radius of 300 pc shows that our models are well within the predicted range of dark matter content for satellites so faint. We therefore conclude that it could be possible that the two galaxies constitute a bound pair.