The Pan-Andromeda Archaeological Survey is a survey of $>400$ square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of a $L_\star$ galaxy group to large projected galactocentric radii. Here, we collate and summarise the current status of our knowledge of the substructures in the stellar halo of M31, and discuss connections between these features. We estimate that the 13 most distinctive substructures were produced by at least 5 different accretion events, all in the last 3 or 4 Gyrs. We suggest that a few of the substructures furthest from M31 may be shells from a single accretion event. We calculate the luminosities of some prominent substructures for which previous estimates were not available, and we estimate the stellar mass budget of the outer halo of M31. We revisit the problem of quantifying the properties of a highly structured dataset; specifically, we use the OPTICS clustering algorithm to quantify the hierarchical structure of M31's stellar halo, and identify three new faint structures. M31's halo, in projection, appears to be dominated by two `mega-structures', that can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into smaller and smaller structures based on the stellar spatial clustering. We conclude that OPTICS is a powerful algorithm that could be used in any astronomical application involving the hierarchical clustering of points. The publication of this article coincides with the public release of all PAndAS data products.

We use simulations to build an analytic probability distribution for the asphericity in the satellite distribution around Local Group (LG) type galaxies in the Lambda Cold Dark Matter (LCDM) paradigm. We use this distribution to estimate the atypicality of the satellite distributions in the LG even when the underlying simulations do not have enough systems fully resembling the LG in terms of its typical masses, separation and kinematics. We demonstrate the method using three different simulations (Illustris-1, Illustris-1-Dark and ELVIS) and a number of satellites ranging from 11 to 15. Detailed results differ greatly among the simulations suggesting a strong influence of the typical DM halo mass, the number of satellites and the simulated baryonic effects. However, there are three common trends. First, at most $2\%$ of the pairs are expected to have satellite distributions with the same asphericity as the LG; second, at most $80\%$ of the pairs have a halo with a satellite distribution as aspherical as in M31; and third, at most $4\%$ of the pairs have a halo with satellite distribution as planar as in the MW. These quantitative results place the LG at the level of a $3\sigma$ outlier in the LCDM paradigm. We suggest that understanding the reasons for this atypicality requires quantifying the asphericity probability distribution as a function of halo mass and large scale environment. The approach presented here can facilitate that kind of study and other comparisons between different numerical setups and choices to study satellites around LG pairs in simulations.

The planar arrangement of nearly half the satellite galaxies of M31 has been a source of mystery and speculation since it was discovered. With a growing number of other host galaxies showing these satellite galaxy planes, their stability and longevity have become central to the debate on whether the presence of satellite planes are a natural consequence of prevailing cosmological models, or represent a challenge. Given the dependence of their stability on host halo shape, we look into how a galaxy plane's dark matter environment influences its longevity. An increased number of dark matter subhalos results in increased interactions that hasten the deterioration of an already-formed plane of satellite galaxies in spherical dark halos. The role of total dark matter mass fraction held in subhalos in dispersing a plane of galaxies present non trivial effects on plane longevity as well. But any misalignments of plane inclines to major axes of flattened dark matter halos lead to their lifetimes being reduced to < 3 Gyrs. Distributing > 40% of total dark mass in subhalos in the overall dark matter distribution results in a plane of satellite galaxies that is prone to change through the 5 Gyr integration time period.

The recently discovered vast thin plane of dwarf satellites orbiting the Andromeda Galaxy (M31) adds to the mystery of the small scale distribution of the Local Group's galaxy population. Such well defined planar structures are apparently rare occurrences in cold dark matter cosmological simulations, and we lack a coherent explanation of their formation and existence. In this paper, we explore the long-term survivability of thin planes of dwarfs in galactic halos, focusing, in particular, on systems mimicking the observed Andromeda distribution. The key results show that, in general, planes of dwarf galaxies are fragile, sensitive to the shape of the dark matter halo and other perturbing effects. In fact, long lived planes of satellites only exist in polar orbits in spherical dark matter halos, presenting a challenge to the observed Andromeda plane which is significantly tilted with respect to the optical disk. Our conclusion is that, in standard cosmological models, planes of satellites are generally short lived, and hence we must be located at a relatively special time in the evolution of the Andromeda Plane, lucky enough to see its coherent pattern.

NGC147, NGC185 and CassiopeiaII (CassII) have similar positions in the sky, distances and measured line of sight velocities. This proximity in phase space suggests that these three satellites of M31 form a subgroup within the Local Group. Nevertheless, the differences in their star formation history and interstellar medium, and the recent discovery of a stellar stream in NGC~147, combined with the lack of tidal features in the other two satellites, are all indications of complex and diverse interactions between M31 and these three satellites. We use a genetic algorithm to explore the different orbits that these satellites can have and select six sets of orbits that could best explain the observational features of the NGC147, NGC185 and CassII satellites. The parameters of these orbits are then used as a starting point for N-body simulations. We present models for which NGC147, NGC185 and CassII are a bound group for a total time of at least one Gyr but still undergo different interactions with M31 and as a result NGC147 has a clear stellar stream whereas the other two satellites have no significant tidal features. This result shows that it is possible to find solutions that reproduce the contrasting properties of the satellites and for which NGC147-NGC185-CassII have been gravitationally bound.

In a recent contribution, Bahl \& Baumgardt investigated the incidence of planar alignments of satellite galaxies in the Millennium-II simulation, and concluded that vast thin planes of dwarf galaxies, similar to that observed in the Andromeda galaxy (M31), occur frequently by chance in $\Lambda$-Cold Dark Matter cosmology. However, their analysis did not capture the essential fact that the observed alignment is simultaneously radially extended, yet thin, and kinematically unusual. With the caveat that the Millennium-II simulation may not have sufficient mass resolution to identify confidently simulacra of low-luminosity dwarf galaxies, we re-examine that simulation for planar structures, using the same method as employed by Ibata et al. (2013) on the real M31 satellites. We find that 0.04\% of host galaxies display satellite alignments that are at least as extreme as the observations, when we consider their extent, thickness and number of members rotating in the same sense. We further investigate the angular momentum properties of the co-planar satellites, and find that the median of the specific angular momentum derived from the line of sight velocities in the real M31 structure ($1.3\times10^4$ km/s kpc) is very high compared to systems drawn from the simulations. This analysis confirms that it is highly unlikely that the observed structure around the Andromeda galaxy is due to a chance occurrence. Interestingly, the few extreme systems that are similar to M31 arise from the accretion of a massive sub-halo with its own spatially-concentrated entourage of orphan satellites.

The long term time evolution of tidal dwarf satellite galaxies with two different initial densities orbiting a host galaxy that resembles the Milky Way has been studied using a large set of Newtonian N-Body simulations. From the simulations two maps of the orbital conditions that lead to quasi-equilibrium objects were constructed. It has been found that several orbits of the satellites allow for the existence, for about 1 Gyr or more, of out-of-equilibrium bodies with high apparent mass-to-light ratios. Within this framework the satellites in the quasi-stable phase reproduce the observed satellite properties for about 16% of the orbit for high density progenitors, and for about 66% for progenitors with lower densities An additional simulation for a single satellite with initial mass of 10^7 Msun and Plummer radius of 0.15 kpc leads to remnants in the quasi- equilibrium phase that simultaneously reproduce remarkably well the observational quantities of the UFDGs of the Milky Way. This satellite in the quasi-stable phase reproduces the observed satellite properties for about 42% of the orbit. The results suggest that a fraction of the observed satellites could plausibly be galaxies without dark matter that have true M/L ratios much lower than those measured. The inflated M/L ratios arise because they are observed at the right time, along the right orbit and during the quasi-equilibrium phase of their evolution. This is a viable explanation for the high M/L ratios observed in all satellites as long as the satellites are preferentially on certain orbits and are observed at certain times. This could arise within the TDG scenario if all satellites are created at the same time along a few specific orbits that are particularly susceptible to the quasi-equilibrium phase.