Halo Of M31 Research Articles

The Isaac Newton Telescope Monitoring Survey of Local Group Dwarf Galaxies. IV. The Star Formation History of Andromeda VII Derived from Long-period Variable Stars

Abstract We have examined the star formation history (SFH) of Andromeda VII (And VII), the brightest and most massive dwarf spheroidal (dSph) satellite of the Andromeda galaxy (M31). Although M31 is surrounded by several dSph companions with old stellar populations and low metallicity, it has a metal-rich stellar halo with an age of 6–8 Gyr. This indicates that any evolutionary association between the stellar halo of M31 and its dSph system is frail. Therefore, the question is whether And VII (a high-metallicity dSph located ∼220 kpc from M31) can be associated with M31's young, metal-rich halo. Here we perform the first reconstruction of the SFH of And VII employing long-period variable (LPV) stars. As the most evolved asymptotic giant branch and red supergiant stars, the birth mass of LPVs can be determined by connecting their near-infrared photometry to theoretical evolutionary tracks. We found 55 LPV candidates within two half-light radii, using multiepoch imaging with the Isaac Newton Telescope in the i and V bands. Based on their birth mass function, the star formation rate (SFR) of And VII was obtained as a function of cosmic time. The main epoch of star formation occurred ≃ 6.2 Gyr ago with an SFR of 0.006 ± 0.002 M ⊙ yr−1. Over the past 6 Gyr, we find slow star formation, which continued until 500 Myr ago with an SFR ∼ 0.0005 ± 0.0002 M ⊙ yr−1. We determined And VII’s stellar mass M = (13.3 ± 5.3) × 106 M ⊙ within a half-light radius and metallicity Z = 0.0007, and we also derived its distance modulus of μ = 24.38 mag.

Observing the Stellar Halo of Andromeda in Cosmological Simulations: The AURIGA2PANDAS Pipeline

Abstract We present a direct comparison of the Pan-Andromeda Archaeological Survey (PAndAS) observations of the stellar halo of M31 with the stellar halos of six galaxies from the Auriga simulations. We process the simulated halos through the Auriga2PAndAS pipeline and create PAndAS-like mocks that fold in all observational limitations of the survey data (foreground contamination from the Milky Way stars, incompleteness of the stellar catalogs, photometric uncertainties, etc.). This allows us to study the survey data and the mocks in the same way and generate directly comparable density maps and radial density profiles. We show that the simulations are overall compatible with the observations. Nevertheless, some systematic differences exist, such as a preponderance for metal-rich stars in the mocks. While these differences could suggest that M31 had a different accretion history or has a different mass compared with the simulated systems, it is more likely a consequence of an underquenching of the star formation history of galaxies, related to the resolution of the Auriga simulations. The direct comparison enabled by our approach offers avenues to improve our understanding of galaxy formation as they can help pinpoint the observable differences between observations and simulations. Ideally, this approach will be further developed through an application to other stellar halo simulations. To facilitate this step, we release the pipeline to generate the mocks, along with the six mocks presented and used in this contribution.

Dark matter explanations of the gamma-ray excesses from the Galactic Center and M31

The presence of an excess $\ensuremath{\gamma}$-ray signal toward the Galactic Center (GC) has now been well established, and is known as the GC excess. Leading explanations for the signal include mismodeling of the Galactic diffuse emission along the line of sight, an unresolved population of millisecond pulsars, and/or the annihilation of dark matter (DM). Recently, evidence for another excess $\ensuremath{\gamma}$-ray signal has been reported toward the outer halo of M31. In this work we interpret the excess signals from both the GC and outer halo of M31 in the framework of DM annihilation, and show that the two spectra are consistent with a DM origin once $J$-factors are taken into account. We further compare the excesses to models of DM annihilation, and determine the corresponding best-fit parameters. We find good fits to the spectrum both in two body and four body annihilation modes.

Dark matter interpretation of the Fermi-LAT observations toward the outer halo of M31.

An excess γ-ray signal toward the outer halo of M31 has recently been reported. Although other explanations are plausible, the possibility that it arises from dark matter (DM) is valid. In this work we interpret the excess in the framework of DM annihilation, using as our representative case WIMP DM annihilating to bottom quarks, and we perform a detailed study of the systematic uncertainty in the J-factor for the M31 field. We find that the signal favors a DM particle with a mass of ~45-72 GeV. While the mass is well constrained, the systematic uncertainty in the cross section spans 3 orders of magnitude, ranging from ~5 × 10-27-5 × 10-24 cm3 s-1. This high uncertainty is due to two main factors, namely, an uncertainty in the substructure nature and geometry of the DM halos for both M31 and the Milky Way (MW), and correspondingly, an uncertainty in the contribution to the signal from the MW's DM halo along the line of sight. However, under the conditions that the minimum subhalo mass is ≲10-6 M ⊙ and the actual contribution from the MW's DM halo along the line of sight is at least ~30% of its total value, we show that there is a large overlap with the DM interpretations of both the Galactic center (GC) excess and the antiproton excess, while also being compatible with the limits for the MW dwarf spheroidals. More generally, we summarize the results from numerous complementary DM searches in the energy range 10 GeV-300 GeV corresponding to the GC excess and identify a region in parameter space that still remains viable for discovery of the DM particle.

Detections of Dust in the Outskirts of M31 and M33

Abstract M31 and M33 serve as ideal places to study distributions of dust in the outskirts of spiral galaxies. In this Letter, using about 0.2 million stars selected from the LAMOST data and combining precise photometry and parallaxes from the Gaia DR2, we have constructed a two-dimensional foreground dust reddening map toward the M31 and M33 region (1112 ≤ gl ≤ 1362, −365 ≤ gb ≤ −165). The map has a typical spatial resolution of about 12′ and precision of 0.01 mag. The complex structure of dust clouds toward M31 is revealed. By carefully removing the foreground extinction from the dust reddening map of Schlegel et al., we thus have obtained a residual map to study dust distributions in the outskirts of M31 and M33. A large amount of dust is detected in the M31 halo out to a distance of over 100 kpc. Dust in the M31 disk is found to extend out to about 2.5 times its optical radius, with a distribution that is consistent with either an exponential disk with a scale length of 7.2 kpc or two disks with a scale length of 11.1 kpc within its optical radius and 18.3 kpc beyond its optical radius. Dust in the disk of M33 is also found to extend out to about 2.5 times its optical radius, its distribution beyond one optical radius is consistent with an exponential disk with a scale length of 5.6 kpc. Our results provide new clues to the distributions and cycling of dust in galaxies.

Elemental Abundances in M31: Properties of the Inner Stellar Halo*

Abstract We present measurements of [Fe/H] and [α/Fe] for 128 individual red giant branch stars (RGB) in the stellar halo of M31, including its Giant Stellar Stream (GSS), obtained using spectral synthesis of low- and medium-resolution Keck/DEIMOS spectroscopy ( and 6000, respectively). We observed four fields in M31's stellar halo (at projected radii of 9, 18, 23, and 31 kpc), as well as two fields in the GSS (at 33 kpc). In combination with existing literature measurements, we have increased the sample size of [Fe/H] and [α/Fe] measurements from 101 to a total of 229 individual M31 RGB stars. From this sample, we investigate the chemical abundance properties of M31's inner halo, finding and . Between 8 and 34 kpc, the inner halo has a steep [Fe/H] gradient (−0.025 ± 0.002 dex kpc−1) and negligible [α/Fe] gradient, where substructure in the inner halo is systematically more metal-rich than the smooth component of the halo at a given projected distance. Although the chemical abundances of the inner stellar halo are largely inconsistent with that of present-day dwarf spheroidal (dSph) satellite galaxies of M31, we identified 22 RGB stars kinematically associated with the smooth component of the stellar halo that have chemical abundance patterns similar to M31 dSphs. We discuss formation scenarios for M31's halo, concluding that these dSph-like stars may have been accreted from galaxies of similar stellar mass and star formation history, or of higher stellar mass and similar star formation efficiency.

A bright, high rotation-measure FRB that skewers the M33 halo

ABSTRACT We report the detection of a bright fast radio burst, FRB 191108, with Apertif on the Westerbork Synthesis Radio Telescope. The interferometer allows us to localize the FRB to a narrow 5 arcsec × 7 arcmin ellipse by employing both multibeam information within the Apertif phased-array feed beam pattern, and across different tied-array beams. The resulting sightline passes close to Local Group galaxy M33, with an impact parameter of only 18 kpc with respect to the core. It also traverses the much larger circumgalactic medium (CGM) of M31, the Andromeda Galaxy. We find that the shared plasma of the Local Group galaxies could contribute ∼10 per cent of its dispersion measure of 588 pc cm−3. FRB 191108 has a Faraday rotation measure (RM) of +474 $\pm \, 3$ rad m−2, which is too large to be explained by either the Milky Way or the intergalactic medium. Based on the more moderate RMs of other extragalactic sources that traverse the halo of M33, we conclude that the dense magnetized plasma resides in the host galaxy. The FRB exhibits frequency structure on two scales, one that is consistent with quenched Galactic scintillation and broader spectral structure with Δν ≈ 40 MHz. If the latter is due to scattering in the shared M33/M31 CGM, our results constrain the Local Group plasma environment. We found no accompanying persistent radio sources in the Apertif imaging survey data.

The Assembly History of M87 through Radial Variations in Chemical Abundances of Its Field Star and Globular Cluster Populations

Abstract We present an extensive study of spectroscopically derived chemical abundances for M87 and its globular cluster (GC) system. Using observations from the Mitchell spectrograph at McDonald, LRIS at Keck, and Hectospec on the MMT, we derive new metallicity gradients from ∼2 to 140 kpc. We use a novel hierarchical statistical framework to simultaneously separate the GC system into subpopulations while measuring the metallicity gradients of those subpopulations. We create physically motivated spectral stacks of the GC subpopulations by leveraging the output of this statistical framework to perform the first application of abundance tagging in a massive early-type galaxy to better constrain the origins of the GC subpopulations and thus the assembly history of M87. We find a metal-poor, α-enhanced population of GCs in both the inner and outer halos unanticipated by current cosmological simulations of galaxy evolution. We use the remarkably flat metallicity gradients we find for both the metal-rich and metal-poor GC subpopulations in the inner halo as tentative evidence that some amount of the metal-poor GCs formed directly in the halo of M87 at high redshift.

Keck Cosmic Web Imager (KCWI) spectra of globular clusters and ultracompact dwarfs in the halo of M87

ABSTRACT Using the Keck Cosmic Web Imager, we obtain spectra of several globular clusters (GCs), ultracompact dwarfs (UCDs), and the inner halo starlight of M87, at a similar projected galactocentric radius of ∼5 kpc. This enables us, for the first time, to apply the same stellar population analysis to the GCs, UCDs, and starlight consistently to derive ages, metallicities, and alpha-element abundances in M87. We find evidence for a dual stellar population in the M87 halo light, i.e. an ∼80 per cent component by mass that is old and metal-rich and a ∼20 per cent component that is old but metal-poor. Two red GCs share similar stellar populations to the halo light suggesting they may have formed contemporaneously with the dominant halo component. Three UCDs, and one blue GC, have similar stellar populations, with younger mean ages, lower metallicities, and near solar alpha-element abundances. Combined with literature data, our findings are consistent with the scenario that UCDs are the remnant nucleus of a stripped galaxy. We further investigate the discrepancy in the literature for M87’s kinematics at large radii, favouring a declining velocity dispersion profile. This work has highlighted the need for more self-consistent studies of galaxy haloes.

Elemental Abundances in M31: A Comparative Analysis of Alpha and Iron Element Abundances in the the Outer Disk, Giant Stellar Stream, and Inner Halo of M31

Abstract We measured [Fe/H] and [α/Fe] using spectral synthesis of low-resolution stellar spectroscopy for 70 individual red-giant-branch stars across four fields spanning the outer disk, Giant Stellar Stream (GSS), and inner halo of M31. Fields at M31-centric projected distances of 23 kpc in the halo, 12 kpc in the halo, 22 kpc in the GSS, and 26 kpc in the outer disk are α-enhanced, with [α/Fe] = 0.43, 0.50, 0.41, and 0.58, respectively. The 23 and 12 kpc halo fields are relatively metal-poor, with [Fe/H] = −1.54 and −1.30, whereas the 22 kpc GSS and 26 kpc outer disk fields are relatively metal-rich with [Fe/H] = −0.84 and −0.92, respectively. For fields with substructure, we separated the stellar populations into kinematically hot stellar halo components and kinematically cold components. We did not find any evidence of a radial [α/Fe] gradient along the high surface brightness core of the GSS between ∼17 and 22 kpc. However, we found tentative suggestions of a negative radial [α/Fe] gradient in the stellar halo, which may indicate that different progenitor(s) or formation mechanisms contributed to the build up of the inner versus outer halo. Additionally, the [α/Fe] distribution of the metal-rich ([Fe/H] > −1.5), smooth inner stellar halo (rproj ≲ 26 kpc) is inconsistent with having formed from the disruption of a progenitor(s) similar to present-day M31 satellite galaxies. The 26 kpc outer disk is most likely associated with the extended disk of M31, where its high α-enhancement provides support for an episode of rapid star formation in M31's disk possibly induced by a major merger.

Two major accretion epochs in M31 from two distinct populations of globular clusters.

Large galaxies grow through the accumulation of dwarf galaxies1,2. In principle it is possible to trace this growth history viathe properties of a galaxy's stellar halo3-5. Previous investigations of the galaxy Messier 31 (M31, Andromeda) have shown that outside a galactocentric radius of 25kiloparsecs the population of halo globular clusters is rotating in alignment with the stellar disk6,7, as are more centrally located clusters8,9. The M31 halo also contains coherent stellar substructures, along with a smoothly distributed stellar component10-12. Many of the globular clusters outside a radius of 25 kiloparsecs are associated with the most prominent substructures, but some are part of the smooth halo13. Here we report an analysis of the kinematics of these globular clusters. We find two distinct populations rotating perpendicular to each other. The rotation axis for the population associated with the smooth halo is aligned with the rotation axis for the plane of dwarf galaxies14 that encircles M31. We interpret these separate cluster populations as arising from two major accretion epochs, probably separated by billions of years. Stellar substructures from the first epoch are gone, but those from the more recent second epoch still remain.

Elemental Abundances in M31: First Alpha and Iron Abundance Measurements in M31's Giant Stellar Stream∗

Abstract We present the first measurements of [Fe/H] and abundances, obtained using spectral synthesis modeling, for red giant branch stars in M31's giant stellar stream (GSS). The spectroscopic observations, obtained at a projected distance of 17 kpc from M31's center, yielded 61 stars with [Fe/H] measurements, including 21 stars with measurements, from 112 targets identified as M31 stars. The [Fe/H] measurements confirm the expectation from photometric metallicity estimates that stars in this region of M31's halo are relatively metal rich compared to stars in the Milky Way’s inner halo: more than half the stars in the field, including those not associated with kinematically identified substructure, have [Fe/H] abundances . The stars in this field are α-enhanced at lower metallicities, while decreases with increasing [Fe/H] above metallicities of [Fe/H] ≳ −0.9. Three kinematical components have been previously identified in this field: the GSS, a second kinematically cold feature of unknown origin, and M31's kinematically hot halo. We compare probabilistic [Fe/H] and distribution functions for each of the components. The GSS and the second kinematically cold feature have very similar abundance distributions, while the halo component is more metal poor. Although the current sample sizes are small, a comparison of the abundances of stars in the GSS field with abundances of M31 halo and dSph stars from the literature indicate that the progenitor of the stream was likely more massive, and experienced a higher efficiency of star formation, than M31's existing dSphs or the dEs NGC 147 and NGC 185.

A rogues gallery of Andromeda’s dwarf galaxies – II. Precise distances to 17 faint satellites

ABSTRACT We present new horizontal branch (HB) distance measurements to 17 of the faintest known M31 satellites (−6 ≲ MV ≲ −13) based on deep Hubble Space Telescope (HST) imaging. The colour–magnitude diagrams extend ∼1–2 mag below the HB, which provides for well-defined HBs, even for faint galaxies in which the tip of the red giant branch (TRGB) is sparsely populated. We determine distances across the sample to an average precision of 4 per cent (∼30 kpc at 800 kpc). We find that the majority of these galaxies are in good agreement, though slightly farther (0.1–0.2 mag) when compared to recent ground-based TRGB distances. Two galaxies (And IX and And XVII) have discrepant HST and ground-based distances by ∼0.3 mag (∼150 kpc), which may be due to contamination from Milky Way foreground stars and/or M31 halo stars in sparsely populated TRGB regions. We use the new distances to update the luminosities and structural parameters for these 17 M31 satellites. The new distances do not substantially change the spatial configuration of the M31 satellite system. We comment on future prospects for precise and accurate HB distances for faint galaxies in the Local Group and beyond.

Elemental Abundances in M31: Alpha and Iron Element Abundances from Low-resolution Resolved Stellar Spectroscopy in the Stellar Halo

Abstract Measurements of [Fe/H] and [α/Fe] can probe the minor merging history of a galaxy, providing a direct way to test the hierarchical assembly paradigm. While measurements of [α/Fe] have been made in the stellar halo of the Milky Way (MW), little is known about detailed chemical abundances in the stellar halo of M31. To make progress with existing telescopes, we apply spectral synthesis to low-resolution DEIMOS spectroscopy (R ∼ 2500 at 7000 Å) across a wide spectral range (4500 Å < λ < 9100 Å). By applying our technique to low-resolution spectra of 170 giant stars in five MW globular clusters, we demonstrate that our technique reproduces previous measurements from higher resolution spectroscopy. Based on the intrinsic dispersion in [Fe/H] and [α/Fe] of individual stars in our combined cluster sample, we estimate systematic uncertainties of ∼0.11 dex and ∼0.09 dex in [Fe/H] and [α/Fe], respectively. We apply our method to deep, low-resolution spectra of 11 red giant branch stars in the smooth halo of M31, resulting in higher signal-to-noise ratios per spectral resolution element compared to DEIMOS medium-resolution spectroscopy, given the same exposure time and conditions. We find [α/Fe] = 0.49 ± 0.29 dex and [Fe/H] = −1.59 ± 0.56 dex for our sample. This implies that—much like the MW—the smooth halo field of M31 is likely composed of disrupted dwarf galaxies with truncated star formation histories that were accreted early in the halo’s formation.

Probing Galactic Halos with Fast Radio Bursts

The precise localization (<1'') of multiple fast radio bursts (FRBs) to z>0.1 galaxies has confirmed that the dispersion measures (DMs) of these enigmatic sources afford a new opportunity to probe the diffuse ionized gas around and in between galaxies. In this manuscript, we examine the signatures of gas in dark matter halos (aka halo gas) on DM observations in current and forthcoming FRB surveys. Combining constraints from observations of the high velocity clouds, OVII absorption, and the DM to the Large Magellanic Cloud with hydrostatic models of halo gas, we estimate that our Galactic halo will contribute ${\rm DM}_{\rm MW,halo} \approx 50-80 \rm pc/cm^{-3}$ from the Sun to 200 kpc independent of any contribution from the Galactic ISM. Extending analysis to the Local Group, we demonstrate that M31's halo will be easily detected by high-sample FRB surveys (e.g. CHIME) although signatures from a putative Local Group medium may compete. We then review current empirical constraints on halo gas in distant galaxies and discuss the implications for their DM contributions. We further examine the DM probability distribution function of a population of FRBs at z >> 0 using an updated halo mass function and new models for the halo density profile. Lastly, we illustrate the potential of FRB experiments for resolving the baryonic fraction of halos by analyzing simulated sightlines through the CASBaH survey. All of the code and data products of our analysis are available at https://github.com/FRBs.

The outer halo globular cluster system of M31 – III. Relationship to the stellar halo

We utilise the final catalogue from the Pan-Andromeda Archaeological Survey to investigate the links between the globular cluster system and field halo in M31 at projected radii $R_p=25-150$ kpc. In this region the cluster radial density profile exhibits a power-law decline with index $\Gamma=-2.37\pm0.17$, matching that for the stellar halo component with [Fe/H] $<-1.1$. Spatial density maps reveal a striking correspondence between the most luminous substructures in the metal-poor field halo and the positions of many globular clusters. By comparing the density of metal-poor halo stars local to each cluster with the azimuthal distribution at commensurate radius, we reject the possibility of no correlation between clusters and field overdensities with high confidence. We use our stellar density measurements and previous kinematic data to demonstrate that $\approx35-60\%$ of clusters exhibit properties consistent with having been accreted into the outskirts of M31 at late times with their parent dwarfs. Conversely, at least $\sim40\%$ of remote clusters show no evidence for a link with halo substructure. The radial density profile for this subgroup is featureless and closely mirrors that observed for the apparently smooth component of the metal-poor stellar halo. We speculate that these clusters are associated with the smooth halo; if so, their properties appear consistent with a scenario where the smooth halo was built up at early times via the destruction of primitive satellites. In this picture the entire M31 globular cluster system outside $R_p=25$ kpc comprises objects accumulated from external galaxies over a Hubble time of growth.

The Large-scale Structure of the Halo of the Andromeda Galaxy. II. Hierarchical Structure in the Pan-Andromeda Archaeological Survey

Abstract The Pan-Andromeda Archaeological Survey is a survey of &gt;400 square degrees centered on the Andromeda (M31) and Triangulum (M33) galaxies that has provided the most extensive panorama of an L ⋆ galaxy group to large projected galactocentric radii. Here, we collate and summarize 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 Gyr. We suggest that a few of the substructures farthest 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 data set; 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,” which can be considered as the two most significant branches of a merger tree produced by breaking M31's stellar halo into increasingly 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.

The Next Generation Virgo Cluster Survey (NGVS). XXXI. The Kinematics of Intracluster Globular Clusters in the Core of the Virgo Cluster

Abstract Intracluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy M87. This study focuses on the Virgo core, for which the combination of Next Generation Virgo Cluster Survey photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at R &gt; 320 kpc. By kinematically tagging M87 halo and intracluster GCs, we find that (1) the M87 halo has a smaller fraction (52 ± 3%) of blue clusters with respect to the IC counterpart (77 ± 10%), (2) the (g′−r′)0 versus (i′−z′)0 color–color diagrams reveal a galaxy population that is redder than the IC population, which may be due to a different composition in chemical abundance and progenitor mass, and (3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, S N,ICL = 10.2 ± 4.8, is consistent with what is observed for the population of quenched, low-mass galaxies within 1 Mpc from the cluster’s center. The IC population at Virgo’s center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of .

Tidally Induced Morphology of M33 in Hydrodynamical Simulations of Its Recent Interaction with M31

Abstract We present a hydrodynamical model of M33 and its recent interaction with M31. This scenario was previously proposed in the literature in order to explain the distorted gaseous and stellar disks of M33, as well as the increased star formation rate in both objects around 2 Gyr ago. We used an orbit integration scheme to find which estimate of the transverse velocity of M31 favors the interaction scenario more and then tried to reproduce it in our simulations. M33 was modeled as a stellar and gaseous disk embedded in a live dark matter halo, while M31 was approximated only with a live dark halo. In the simulations, the two galaxies passed each other with a pericenter distance of 37 kpc. Tides excited a two-armed spiral structure in the M33 disk, which is found to be the predominant spiral signal in the observed galaxy and has long been known as a feature easily induced by tidal interactions. We found that the gaseous warp produced by the interaction did not resemble enough the observed one, and we performed an additional simulation including the hot gas halo of M31 to show that this feature can be properly reproduced by tidal forces and ram pressure stripping acting simultaneously on the gaseous disk. In addition to the spiral arms, tidal forces produced a stellar stream similar to the one observed and triggered a star formation burst at radii similar to where it is observed.

ΛCDM predictions for the satellite population of M33

Triangulum (M33) is the most massive satellite galaxy of Andromeda (M31), with a stellar mass of about $3\times10^9\; M_{\odot}$. Based on abundance matching techniques, M33's total mass at infall is estimated to be of order $10^{11}\; M_{\odot}$. $\Lambda$CDM simulations predict that M33-mass halos host several of their own low mass satellite companions, yet only one candidate M33 satellite has been detected in deep photometric surveys to date. This `satellites of satellites' hierarchy has recently been explored in the context of the dwarf galaxies discovered around the Milky Way's Magellanic Clouds in the Dark Energy Survey. Here we discuss the number of satellite galaxies predicted to reside within the virial radius ($\sim$160 kpc) of M33 based on $\Lambda$CDM simulations. We also calculate the expected number of satellite detections in N fields of data using various ground--based imagers. Finally, we discuss how these satellite population predictions change as a function of M33's past orbital history. If M33 is on its first infall into M31's halo, its proposed satellites are expected to remain bound to M33 today. However, if M33 experienced a recent tidal interaction with M31, the number of surviving satellites depends strongly on the distance achieved at pericenter due to the effects of tidal stripping. We conclude that a survey extending to $\sim$100 kpc around M33 would be sufficient to constrain its orbital history and a majority of its satellite population. In the era of WFIRST, surveys of this size will be considered small observing programs.