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 (M31) galaxy subtends nearly 100 sq. deg. on the sky, with severe contamination from the Milky Way halo stars whose surface density displays a steep gradient across the entire M31 field-of-view. Planetary Nebulae (PNe) are a population of stars firmly associated with M31, that are excellent tracers of light, chemistry and motion in galaxies. We present a 16 sq. deg. survey of the disk and inner halo of M31 with MegaCam@CFHT to identify PNe, characterize their luminosity-specific PN number and luminosity function (PNLF). PNe were identified based on their bright OIII 5007 $\unicode{x212B}$ emission and absence of a continuum. Subsamples of the faint PNe were independently confirmed by matching with resolved Hubble Space Telescope sources from the PHAT survey and spectroscopic follow-up observations with HectoSpec@MMT. The current survey reaches 2 mag fainter than the previous most-sensitive survey. We identify 4289 PNe, of which only 1099 were previously known. By comparing the PN number density with the surface brightness profile of M31 out to ~30 kpc along the minor-axis, we find that the stellar population in the inner halo has a 7 times larger luminosity-specific PN number value than that of the disk. It indicates that the stellar population at deprojected minor-axis radii larger than ~10 kpc is different from that in the M31 disk. We measure the PNLF and find a bright cut-off and a slope consistent with the previous determination by Ciardullo et al. (1989). Interestingly, it shows a significant rise at the faint end, present in all radial bins covered by the survey, much steeper than that observed for the Magellanic clouds and Milky Way bulge. M31 shows two major episodes of star formation and the rise in the faint end of the PNLF is possibly associated with the older stellar population. It may also be a result of varying opacity of the PNe.

M49 (NGC 4472) is the dominant galaxy in subcluster B of the Virgo Cluster, and a benchmark for studying the build-up of the extended halos of brightest group galaxies in the outskirts of galaxy clusters. We investigate the kinematics in the outer halo of M49, look for substructures, and describe the transition to the surrounding intra-group light. As kinematic tracers we use planetary nebulae (PNe), combining kinematics from the extended Planetary Nebula Spectrograph (PN.S) early-type galaxy survey with our recent deep photometric sample. We study the position-velocity-plane for bright and faint PN populations out to 95 kpc radius, and employ a multi-Gaussian model for the velocity distribution to identify stellar populations with distinct kinematics and histories. We report the detection of stellar-kinematic substructure associated with the interaction of M49 with the dwarf irregular galaxy VCC 1249. We find two kinematically distinct PN populations associated with the main M49 halo and the extended intra-group light (IGL). These have velocity dispersions $\sigma_{halo}\simeq 170$ km/s and $\sigma_{IGL} \simeq400$ km/s at 10-80 kpc radii. The overall luminosity profile and velocity dispersion at $\sim80$ kpc are consistent with a flat circular velocity curve extrapolated from X-ray observations. The dispersion of the PNe associated with the IGL joins onto that of the satellite galaxies in subcluster B at $\sim100$ kpc radius. This is the first time that the transition from halo to IGL is observed based on the velocities of individual stars. Therefore the halo of M49, consisting of at least three distinct components, has undergone an extended accretion history within its parent group potential. The blue colours of the IGL component are consistent with a population of stars formed in low-mass galaxies at redshift $\sim0.5$ that has since evolved passively, as suggested by other data.

The galaxy M49 (NGC 4472) is the brightest early-type galaxy in the Virgo Cluster. It is located in Subcluster B and has an unusually blue, metal-poor outer halo. Planetary nebulae (PNe) are excellent tracers of diffuse galaxy and intragroup light. We present a photometric survey of PNe in the galaxy's extended halo to characterise its PN population, as well as the surrounding intragroup light (IGL) of the Subcluster B. PNe were identified based on their bright [OIII]5007 \AA\ emission and absence of a broad-band continuum. We identify 738 PNe out to a radius of 155 kpc from M49's centre from which we define a complete sample of 624 PNe within a limiting magnitude of m_5007=28.8. Comparing the PN number density to the broad-band stellar surface brightness profile, we find a variation of the PN-specific frequency (alpha-parameter) with radius. The outer halo beyond 60 kpc has a 3.2 times higher alpha-parameter compared to the main galaxy halo, which is likely due to contribution from the surrounding blue IGL. We use the Planetary Nebulae Luminosity Function (PNLF) as an indicator of distance and stellar population. Its slope, which correlates empirically with galaxy type, varies within the inner halo. In the eastern quadrant of M49, the PNLF slope is shallower, indicating an additional localised, bright PN population following an accretion event, likely that of the dwarf irregular galaxy VCC1249. We also determined a distance modulus of mu = 31.29+/-0.08 for M49, corresponding to a physical distance of 18.1+/-0.6 Mpc, which agrees with a recent surface-brightness fluctuations distance. The PN populations in the outer halo of M49 are consistent with the presence of a main Sersic galaxy halo with a slight (B-V) colour gradient of 10${}^{-4}$ mag/arcsec surrounded by intragroup light with a very blue colour of (B-V)=0.25 and a constant surface brightness mu_V=28.0 mag/arcsec${}^2$.

Stellar streams result from the tidal disruption of satellites and star clusters as they orbit a host galaxy, and can be very sensitive probes of the gravitational potential of the host system. We select and study narrow stellar streams formed in a Milky-Way-like dark matter halo of the Aquarius suite of cosmological simulations, to determine if these streams can be used to constrain the present day characteristic parameters of the halo's gravitational potential. We find that orbits integrated in static spherical and triaxial NFW potentials both reproduce the locations and kinematics of the various streams reasonably well. To quantify this further, we determine the best-fit potential parameters by maximizing the amount of clustering of the stream stars in the space of their actions. We show that using our set of Aquarius streams, we recover a mass profile that is consistent with the spherically-averaged dark matter profile of the host halo, although we ignored both triaxiality and time evolution in the fit. This gives us confidence that such methods can be applied to the many streams that will be discovered by the Gaia mission to determine the gravitational potential of our Galaxy.