Showing 1–31 of 31 results for author: Patel, E

Deep Hubble Space Telescope Photometry of LMC and Milky Way Ultra-Faint Dwarfs: A careful look into the magnitude-size relation

Authors: Hannah Richstein, Nitya Kallivayalil, Joshua D. Simon, Christopher T. Garling, Andrew Wetzel, Jack T. Warfield, Roeland P. van der Marel, Myoungwon Jeon, Jonah C. Rose, Paul Torrey, Anna Claire Engelhardt, Gurtina Besla, Yumi Choi, Marla Geha, Puragra Guhathakurta, Evan N. Kirby, Ekta Patel, Elena Sacchi, Sangmo Tony Sohn

Abstract: We present deep Hubble Space Telescope (HST) photometry of ten targets from Treasury Program GO-14734, including six confirmed ultra-faint dwarf galaxies (UFDs), three UFD candidates, and one likely globular cluster. Six of these targets are satellites of, or have interacted with, the Large Magellanic Cloud (LMC). We determine their structural parameters using a maximum-likelihood technique. Using… ▽ More We present deep Hubble Space Telescope (HST) photometry of ten targets from Treasury Program GO-14734, including six confirmed ultra-faint dwarf galaxies (UFDs), three UFD candidates, and one likely globular cluster. Six of these targets are satellites of, or have interacted with, the Large Magellanic Cloud (LMC). We determine their structural parameters using a maximum-likelihood technique. Using our newly derived half-light radius ($r_h$) and $V$-band magnitude ($M_V$) values in addition to literature values for other UFDs, we find that UFDs associated with the LMC do not show any systematic differences from Milky Way UFDs in the magnitude-size plane. Additionally, we convert simulated UFD properties from the literature into the $M_V-r_h$ observational space to examine the abilities of current dark matter (DM) and baryonic simulations to reproduce observed UFDs. Some of these simulations adopt alternative DM models, thus allowing us to also explore whether the $M_V-r_h$ plane could be used to constrain the nature of DM. We find no differences in the magnitude-size plane between UFDs simulated with cold, warm, and self-interacting dark matter, but note that the sample of UFDs simulated with alternative DM models is quite limited at present. As more deep, wide-field survey data become available, we will have further opportunities to discover and characterize these ultra-faint stellar systems and the greater low surface-brightness universe. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 42 pages, 13 figures, 13 tables, 2 appendices; Submitted to AAS Journals

Proper Motions and Orbits of Distant Local Group Dwarf Galaxies from a combination of Gaia and Hubble Data

Authors: Paul Bennet, Ekta Patel, Sangmo Tony Sohn, Andres del Pino, Roeland P. van der Marel, Mattia Libralato, Laura L. Watkins, Antonio Aparicio, Gurtina Besla, Carme Gallart, Mark A. Fardal, Matteo Monelli, Elena Sacchi, Erik Tollerud, Daniel R. Weisz

Abstract: We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used HST as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of… ▽ More We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used HST as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of the galaxies in our sample are on first infall with $>$90\% certainty. The observed star formation histories (SFHs) of these first-infall dwarfs are generally consistent with infalling dwarfs in simulations. The remaining four galaxies have crossed the virial radius of either the MW or M31. When we compare their star formation (SF) and orbital histories we find tentative agreement between the inferred pattern of SF with the timing of dynamical events in the orbital histories. For Leo~I, SF activity rises as the dwarf crosses the MW's virial radius, culminating in a burst of SF shortly before pericenter ($\approx1.7$~Gyr ago). The SF then declines after pericenter, but with some smaller bursts before its recent quenching ($\approx0.3$~Gyr ago). This shows that even small dwarfs like Leo~I can hold on to gas reservoirs and avoid quenching for several Gyrs after falling into their host, which is longer than generally found in simulations. Leo~II, NGC~6822, and IC~10 are also qualitatively consistent with this SF pattern in relation to their orbit, but more tentatively due to larger uncertainties. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 27 pages, 9 Figures, 8 Tables, Submitted to ApJ

On the co-rotation of Milky Way satellites: LMC-mass satellites induce apparent motions in outer halo tracers

Authors: Nicolas Garavito-Camargo, Adrian M. Price-Whelan, Jenna Samuel, Emily C. Cunningham, Ekta Patel, Andrew Wetzel, Kathryn V. Johnston, Arpit Arora, Robyn E. Sanderson, Lehman Garrison, Danny Horta

Abstract: Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing th… ▽ More Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing the inner regions of the halo with respect to the outer halo. Using the Latte suite of FIRE-2 cosmological simulations of MW-mass galaxies, we confirm that the apparent motion of the outer halo induced by the infall of a massive satellite changes the observed distribution of orbital poles of outer-halo tracers, including satellites. We quantify the changes in the distribution of orbital poles using the two-point angular correlation function and find that all satellites induce changes. However, the most massive satellites with pericentric passages between 30-100kpc induce the largest changes. The best LMC-like satellite analog shows the largest change in orbital pole distribution. The dispersion of orbital poles decreases by 20° during the first two pericentric passages. Even when excluding the satellites brought in with the LMC-like satellite, there is clustering of orbital poles. These results suggest that in the MW, the recent pericentric passage of the LMC should have changed the observed distribution of orbital poles of all other satellites. Therefore, studies of kinematically-coherent planes of satellites that seek to place the MW in a cosmological context should account for the existence of a massive satellite like the LMC. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 20 pages, 10 figures. ApJ submitted, Comments are welcome

A physically motivated framework to compare pair fractions of isolated low and high mass galaxies across cosmic time

Authors: Katie Chamberlain, Gurtina Besla, Ekta Patel, Vicente Rodriguez-Gomez, Paul Torrey, Garreth Martin, Kelsey Johnson, Nitya Kallivayalil, David Patton, Sarah Pearson, George Privon, Sabrina Stierwalt

Abstract: Low mass galaxy pair fractions are understudied, and it is unclear whether low mass pair fractions evolve in the same way as more massive systems over cosmic time. In the era of JWST, Roman, and Rubin, selecting galaxy pairs in a self-consistent way will be critical to connect observed pair fractions to cosmological merger rates across all mass scales and redshifts. Utilizing the Illustris TNG100… ▽ More Low mass galaxy pair fractions are understudied, and it is unclear whether low mass pair fractions evolve in the same way as more massive systems over cosmic time. In the era of JWST, Roman, and Rubin, selecting galaxy pairs in a self-consistent way will be critical to connect observed pair fractions to cosmological merger rates across all mass scales and redshifts. Utilizing the Illustris TNG100 simulation, we create a sample of physically associated low mass ($\rm 10^8<M_*<5\times10^9\,M_\odot$) and high mass ($\rm 5\times10^9<M_*<10^{11}\,M_\odot$) pairs between $z=0$ and $4.2$. The low mass pair fraction increases from $z=0$ to $2.5$, while the high mass pair fraction peaks at $z=0$ and is constant or slightly decreasing at $z>1$. At $z=0$, the low mass major (1:4 mass ratio) pair fraction is 4$\times$ lower than high mass pairs, consistent with findings for cosmological merger rates. We show that separation limits that vary with the mass and redshift of the system, such as scaling by the virial radius of the host halo ($r_{\mathrm{sep}}< 1 R_{\rm vir}$), are critical for recovering pair fraction differences between low mass and high mass systems. Alternatively, static physical separation limits applied equivalently to all galaxy pairs do not recover the differences between low and high mass pair fractions, even up to separations of $300$ kpc. Finally, we place isolated mass-analogs of Local Group galaxy pairs, i.e., Milky Way (MW)--M31, MW--LMC, LMC--SMC, in a cosmological context, showing that isolated analogs of LMC--SMC-mass pairs and low-separation ($<50$ kpc) MW--LMC-mass pairs are $2-3\times$ more common at $z\gtrsim2-3$. △ Less

Submitted 6 February, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 18 pages, 5 figures, accepted for publication in ApJ, updated Feb 2024 to reflect accepted version

Machine Learning the Dark Matter Halo Mass of Milky Way-Like Systems

Authors: Elaheh Hayati, Peter Behroozi, Ekta Patel

Abstract: Despite the Milky Way's proximity to us, our knowledge of its dark matter halo is fairly limited, and there is still considerable uncertainty in its halo mass. Many past techniques have been limited by assumptions such as the Galaxy being in dynamical equilibrium as well as nearby galaxies being true satellites of the Galaxy, and/or the need to find large samples of Milky Way analogs in simulation… ▽ More Despite the Milky Way's proximity to us, our knowledge of its dark matter halo is fairly limited, and there is still considerable uncertainty in its halo mass. Many past techniques have been limited by assumptions such as the Galaxy being in dynamical equilibrium as well as nearby galaxies being true satellites of the Galaxy, and/or the need to find large samples of Milky Way analogs in simulations.Here, we propose a new technique based on neural networks that obtains high precision ($<0.14$ dex mass uncertainty) without assuming halo dynamical equilibrium or that neighboring galaxies are all satellites, and which can use information from a wide variety of simulated halos (even those dissimilar to the Milky Way) to improve its performance. This method uses only observable information including satellite orbits, distances to nearby larger halos, and the maximum circular velocity of the largest satellite galaxy. In this paper, we demonstrate a proof-of-concept method on simulated dark matter halos; in future papers in this series, we will apply neural networks to estimate the masses of the Milky Way's and M31's dark matter halos, and we will train variations of these networks to estimate other halo properties including concentration, assembly history, and spin axis. △ Less

Submitted 29 March, 2024; v1 submitted 12 September, 2023; originally announced September 2023.

Modeling the orbital histories of satellites of Milky Way-mass galaxies: testing static host potentials against cosmological simulations

Authors: Isaiah B. Santistevan, Andrew Wetzel, Erik Tollerud, Robyn E Sanderson, Jorge Moreno, Ekta Patel

Abstract: Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modeling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmolog… ▽ More Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modeling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass halos. When a typical surviving satellite fell in ($3.4-9.7$ Gyr ago), the host halo mass and radius were typically $26-86$ per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $r\lesssim50$ kpc, opposite to dark-matter-only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at $z=0$ and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $1.6-4.7$ Gyr ago. Most orbital properties are minimally biased, $\lesssim10$ per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are $\approx20$ per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are $\approx40-80$ per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at $z=0$. △ Less

Submitted 8 December, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 24 pages, 12 figures, 3 appendices, 2 appendix figures. Accepted for publication in MNRAS

The Large Magellanic Cloud's $\sim30$ Kiloparsec Bow Shock and its Impact on the Circumgalactic Medium

Authors: David J. Setton, Gurtina Besla, Ekta Patel, Cameron Hummels, Yong Zheng, Evan Schneider

Abstract: The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the Circumgalactic Medium (CGM) is expected to result in a bow shock that leads the LMC's gaseous disk. In this letter, we use hydrodynamic simulations of the LMC's recent infall to predict the extent of this shock and its effect on the Milky Way's (MW) CGM. The simulations clearly predict the existence of an asy… ▽ More The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the Circumgalactic Medium (CGM) is expected to result in a bow shock that leads the LMC's gaseous disk. In this letter, we use hydrodynamic simulations of the LMC's recent infall to predict the extent of this shock and its effect on the Milky Way's (MW) CGM. The simulations clearly predict the existence of an asymmetric shock with a present day stand-off radius of $\sim6.7$ kpc and a transverse diameter of $\sim30$ kpc. Over the past 500 Myr, $\sim8\%$ of the MW's CGM in the southern hemisphere should have interacted with the shock front. This interaction may have had the effect of smoothing over inhomogeneities and increasing mixing in the MW CGM. We find observational evidence of the existence of the bow shock in recent $Hα$ maps of the LMC, providing a potential explanation for the envelope of ionized gas surrounding the LMC. Furthermore, the interaction of the bow shock with the MW CGM may also explain observations of ionized gas surrounding the Magellanic Stream. Using recent orbital histories of MW satellites, we find that many satellites have likely interacted with the LMC shock. Additionally, the dwarf galaxy Ret2 is currently sitting inside the shock, which may impact the interpretation of reported gamma ray excess in Ret2. This work highlights bow shocks associated with infalling satellites are an under-explored, yet potentially very important dynamical mixing process in the circumgalactic and intracluster media. △ Less

Submitted 17 November, 2023; v1 submitted 21 August, 2023; originally announced August 2023.

Comments: Re-uploaded after acceptance to ApJ Letters. 5 figures and 1 table. Comments welcome!

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σ$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

Metallicity Distribution Functions of 13 Ultra-Faint Dwarf Galaxy Candidates from Hubble Space Telescope Narrowband Imaging

Authors: Sal Wanying Fu, Daniel R. Weisz, Else Starkenburg, Nicolas Martin, Alessandro Savino, Michael Boylan-Kolchin, Patrick Cote, Andrew E. Dolphin, Alexander P. Ji, Nicolas Longeard, Mario L. Mateo, Ekta Patel, Nathan R. Sandford

Abstract: We present uniformly measured stellar metallicities for 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; $M_V = -7.1$ to $-0.8$) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope (HST). This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the… ▽ More We present uniformly measured stellar metallicities for 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; $M_V = -7.1$ to $-0.8$) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope (HST). This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the number of metallicities in all known MW UFDs. We provide the first well-populated MDFs for all galaxies in this sample, with [Fe/H] ranging from -3.0 dex to -2.0 dex, and $σ_{\mbox{[Fe/H]}}$ ranging from 0.3 dex to 0.7 dex. We find a nearly constant [Fe/H] ~ -2.6 over 3 decades in luminosity (~$10^2 - 10^5 L_{\odot}$), suggesting that the mass-metallicity relationship does not hold for such faint systems. We find a larger fraction (24%) of EMP ([Fe/H]<-3.0) stars across our sample compared to the literature (14%), but note that uncertainties in our most metal-poor measurements make this an upper limit. We find 19% of stars in our UFD sample to be metal-rich ([Fe/H]>-2.0), consistent with the sum of literature spectroscopic studies. MW UFDs are known to be predominantly >13 Gyr old, meaning that all stars in our sample are truly ancient, unlike metal-poor stars in the MW, which have a range of possible ages. Our UFD metallicities are not well-matched to known streams in the MW, providing further evidence that known MW substructures are not related to UFDs. We include a catalog of our stars to encourage community follow-up studies, including priority targets for ELT-era observations. △ Less

Submitted 22 November, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 37 pages, 20 figures, ApJ accepted and published

Journal ref: ApJ 958 167 (2023)

The Hubble Space Telescope Survey of M31 Satellite Galaxies II. The Star Formation Histories of Ultra-Faint Dwarf Galaxies

Authors: A. Savino, D. R. Weisz, E. D. Skillman, A. Dolphin, A. A. Cole, N. Kallivayalil, A. Wetzel, J. Anderson, G. Besla, M. Boylan-Kolchin, T. M. Brown, J. S. Bullock, M. L. M. Collins, M. C. Cooper, A. J. Deason, A. L. Dotter, M. Fardal, A. M. N. Ferguson, T. K. Fritz, M. C. Geha, K. M. Gilbert, P. Guhathakurta, R. Ibata, M. J. Irwin, M. Jeon , et al. (12 additional authors not shown)

Abstract: We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find th… ▽ More We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50\% of their stellar mass by $z=5$ (12.6~Gyr ago), similar to known UFDs around the MW, but that 10-40\% of their stellar mass formed at later times. We uncover one remarkable UFD, \A{XIII}, which formed only 10\% of its stellar mass by $z=5$, and 75\% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\times10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs. △ Less

Submitted 13 September, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures, 5 appendices, accepted for publication in ApJ

Evidence for a Massive Andromeda Galaxy Using Satellite Galaxy Proper Motions

Authors: Ekta Patel, Kaisey S. Mandel

Abstract: We present new mass estimates for Andromeda (M31) using the orbital angular momenta of four satellite galaxies (M33, NGC 185, NGC 147, IC 10) derived from existing proper motions, distances, and line-of-sight velocities. We infer two masses for M31: $M_{\rm vir}= 2.85^{+1.47}_{-0.77}\times10^{12}\, M_{\odot}$ using satellite galaxy phase space information derived with HST-based M31 proper motions… ▽ More We present new mass estimates for Andromeda (M31) using the orbital angular momenta of four satellite galaxies (M33, NGC 185, NGC 147, IC 10) derived from existing proper motions, distances, and line-of-sight velocities. We infer two masses for M31: $M_{\rm vir}= 2.85^{+1.47}_{-0.77}\times10^{12}\, M_{\odot}$ using satellite galaxy phase space information derived with HST-based M31 proper motions and $M_{\rm vir}=3.02^{+1.30}_{-0.69}\times10^{12}\, M_{\odot}$ using phase space information derived with the weighted average of HST+Gaia-based M31 proper motions. The precision of our new M31 mass estimates (23-50%) improves by a factor of two compared to previous mass estimates using a similar methodology with just one satellite galaxy and places our results amongst the highest precision M31 estimates in recent literature. Furthermore, our results are consistent with recently revised estimates for the total mass of the Local Group (LG), with the stellar mass--halo mass relation, and with observed kinematic data for both M31 and its entire population of satellites. An M31 mass $> 2.5 \times 10^{12}\, M_{\odot}$ could have major implications for our understanding of LG dynamics, M31's merger and accretion history, and our understanding of LG galaxies in a cosmological context. △ Less

Submitted 20 February, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 21 pages, 6 figures. Submitted to ApJ

DESI Observations of the Andromeda Galaxy: Revealing the Immigration History of our Nearest Neighbor

Authors: Arjun Dey, Joan R. Najita, S. E. Koposov, J. Josephy-Zack, Gabriel Maxemin, Eric F. Bell, C. Poppett, E. Patel, L. Beraldo e Silva, A. Raichoor, D. Schlegel, D. Lang, A. Meisner, Adam D. Myers, J. Aguilar, S. Ahlen, C. Allende Prieto, D. Brooks, A. P. Cooper, K. S. Dawson, A. de la Macorra, P. Doel, A. Font-Ribera, Juan Garcia-Bellido, S. Gontcho A Gontcho , et al. (23 additional authors not shown)

Abstract: We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars:… ▽ More We present DESI observations of the inner halo of M31, which reveal the kinematics of a recent merger - a galactic immigration event - in exquisite detail. Of the 11,416 sources studied in 3.75 hour of on-sky exposure time, 7,438 are M31 sources with well measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars: streams, wedges, and chevrons. While hints of coherent structures have been previously detected in M31, this is the first time they have been seen with such detail and clarity in a galaxy beyond the Milky Way. We find clear kinematic evidence for shell structures in the Giant Stellar Stream, the Northeast Shelf and Western Shelf regions. The kinematics are remarkably similar to the predictions of dynamical models constructed to explain the spatial morphology of the inner halo. The results are consistent with the interpretation that much of the substructure in the inner halo of M31 is produced by a single galactic immigration event 1 - 2 Gyr ago. Significant numbers of metal-rich stars ([Fe/H]$>-0.5$) are present in all of the detected substructures, suggesting that the immigrating galaxy had an extended star formation history. We also investigate the ability of the shells and Giant Stellar Stream to constrain the gravitational potential of M31, and estimate the mass within a projected radius of 125 kpc to be ${\rm log_{10}}\, M_{\rm NFW}(<125\,{\rm kpc})/M_\odot = 11.80_{-0.10}^{+0.12}$. The results herald a new era in our ability to study stars on a galactic scale and the immigration histories of galaxies. △ Less

Submitted 20 January, 2023; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 45 pages, 22 figures, 8 tables; Astrophysical Journal in press; Data at https://zenodo.org/record/6977494

The Hubble Space Telescope Survey of M31 Satellite Galaxies I. RR Lyrae-based Distances and Refined 3D Geometric Structure

Authors: Alessandro Savino, Daniel R. Weisz, Evan D. Skillman, Andrew Dolphin, Nitya Kallivayalil, Andrew Wetzel, Jay Anderson, Gurtina Besla, Michael Boylan-Kolchin, James S. Bullock, Andrew A. Cole, Michelle L. M. Collins, M. C. Cooper, Alis J. Deason, Aaron L. Dotter, Mark Fardal, Annette M. N. Ferguson, Tobias K. Fritz, Marla C. Geha, Karoline M. Gilbert, Puragra Guhathakurta, Rodrigo Ibata, Michael J. Irwin, Myoungwon Jeon, Evan Kirby , et al. (11 additional authors not shown)

Abstract: We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days… ▽ More We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system. △ Less

Submitted 12 September, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 32 pages, 15 figures, accepted for publication on ApJ

Structural parameters and possible association of the Ultra-Faint Dwarfs Pegasus III and Pisces II from deep Hubble Space Telescope photometry

Authors: Hannah Richstein, Ekta Patel, Nitya Kallivayalil, Joshua D. Simon, Paul Zivick, Erik Tollerud, Tobias Fritz, Jack T. Warfield, Gurtina Besla, Roeland P. van der Marel, Andrew Wetzel, Yumi Choi, Alis Deason, Marla Geha, Puragra Guhathakurta, Myoungwon Jeon, Evan N. Kirby, Mattia Libralato, Elena Sacchi, Sangmo Tony Sohn

Abstract: We present deep Hubble Space Telescope (HST) photometry of the ultra-faint dwarf (UFD) galaxies Pegasus III (Peg III) and Pisces II (Psc II), two of the most distant satellites in the halo of the Milky Way (MW). We measure the structure of both galaxies, derive mass-to-light ratios with newly determined absolute magnitudes, and compare our findings to expectations from UFD-mass simulations. For Pe… ▽ More We present deep Hubble Space Telescope (HST) photometry of the ultra-faint dwarf (UFD) galaxies Pegasus III (Peg III) and Pisces II (Psc II), two of the most distant satellites in the halo of the Milky Way (MW). We measure the structure of both galaxies, derive mass-to-light ratios with newly determined absolute magnitudes, and compare our findings to expectations from UFD-mass simulations. For Peg III, we find an elliptical half-light radius of $a_h=1.88^{+0.42}_{-0.33}$ arcminutes ($118^{+31}_{-30}$ pc) and $M_V{=}{-4.17}^{+0.19}_{-0.22}$; for Psc II, we measure $a_h{=}1.31^{+0.10}_{-0.09}$ arcminutes ($69\pm8$ pc) and $M_V{=}{-4.28}^{+0.19}_{-0.16}$. We do not find any morphological features that indicate a significant interaction between the two has occurred, despite their close separation of only $\sim$40 kpc. Using proper motions (PMs) from Gaia early Data Release 3, we investigate the possibility of any past association by integrating orbits for the two UFDs in a MW-only and a combined MW and Large Magellanic Cloud (LMC) potential. We find that including the gravitational influence of the LMC is crucial, even for these outer-halo satellites, and that a possible orbital history exists where Peg III and Psc II experienced a close ($\sim$10-20 kpc) passage about each other just over $\sim$1 Gyr ago, followed by a collective passage around the LMC ($\sim$30-60 kpc) just under $\sim$1 Gyr ago. Considering the large uncertainties on the PMs and the restrictive priors imposed to derive them, improved PM measurements for Peg III and Psc II will be necessary to clarify their relationship. This would add to the rare findings of confirmed pairs of satellites within the Local Group. △ Less

Submitted 13 February, 2024; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 24 pages, 13 figures; 1 appendix; published in ApJ

Journal ref: The Astrophysical Journal, Volume 933, Issue 2, id.217, July 2022

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Authors: Monica Valluri, Solene Chabanier, Vid Irsic, Eric Armengaud, Michael Walther, Connie Rockosi, Miguel A. Sanchez-Conde, Leandro Beraldo e Silva, Andrew P. Cooper, Elise Darragh-Ford, Kyle Dawson, Alis J. Deason, Simone Ferraro, Jaime E. Forero-Romero, Antonella Garzilli, Ting Li, Zarija Lukic, Christopher J. Manser, Nathalie Palanque-Delabrouille, Corentin Ravoux, Ting Tan, Wenting Wang, Risa Wechsler, Andreia Carrillo, Arjun Dey , et al. (7 additional authors not shown)

Abstract: Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on thei… ▽ More Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science. △ Less

Submitted 1 July, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contributed white paper to Snowmass 2021, CF03; minor revisions

R2-D2: Roman and Rubin -- From Data to Discovery

Authors: Suvi Gezari, Misty Bentz, Kishalay De, K. Decker French, Aaron Meisner, Michelle Ntampaka, Robert Jedicke, Ekta Patel, Daniel Perley, Robyn Sanderson, Christian Aganze, Igor Andreoni, Eric F. Bell, Edo Berger, Ian Dell'Antonio, Ryan Foley, Henry Hsieh, Mansi Kasliwal, Joel Kastner, Charles D. Kilpatrick, J. Davy Kirkpatrick, Casey Lam, Karen Meech, Dante Minniti, Ethan O. Nadler , et al. (6 additional authors not shown)

Abstract: The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to… ▽ More The NASA Nancy Grace Roman Space Telescope (Roman) and the Vera C. Rubin Observatory Legacy Survey of Space and Time (Rubin), will transform our view of the wide-field sky, with similar sensitivities, but complementary in wavelength, spatial resolution, and time domain coverage. Here we present findings from the AURA Roman+Rubin Synergy Working group, charged by the STScI and NOIRLab Directors to identify frontier science questions in General Astrophysics, beyond the well-covered areas of Dark Energy and Cosmology, that can be uniquely addressed with Roman and Rubin synergies in observing strategy, data products and archiving, joint analysis, and community engagement. This analysis was conducted with input from the community in the form of brief (1-2 paragraph) "science pitches" (see Appendix), and testimony from "outside experts" (included as co-authors). We identify a rich and broad landscape of potential discoveries catalyzed by the combination of exceptional quality and quantity of Roman and Rubin data, and summarize implementation requirements that would facilitate this bounty of additional science with coordination of survey fields, joint coverage of the Galactic plane, bulge, and ecliptic, expansion of General Investigator and Target of Opportunity observing modes, co-location of Roman and Rubin data, and timely distribution of data, transient alerts, catalogs, value-added joint analysis products, and simulations to the broad astronomical community. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 29 pages, 12 figures, Table of Implementation Recommendations, Appendix of Community Science Pitches, AURA-commissioned whitepaper submitted to the Director of STScI (Ken Sembach) and the Director of NOIRLab (Pat McCarthy)

The Triangulum Extended (TREX) Survey: The Stellar Disk Dynamics of M33 as a Function of Stellar Age

Authors: A. C. N. Quirk, P. Guhathakurta, K. Gilbert, L. Chemin, J. Dalcanton, B. Williams, A. Seth, E. Patel, J. Fung, P. Tangirala, I. Yusufali

Abstract: Triangulum, M33, is a low mass, relatively undisturbed spiral galaxy that offers a new regime in which to test models of dynamical heating. In spite of its proximity, the dynamical heating history of M33 has not yet been well constrained. In this work, we present the TREX Survey, the largest stellar spectroscopic survey across the disk of M33. We present the stellar disk kinematics as a function o… ▽ More Triangulum, M33, is a low mass, relatively undisturbed spiral galaxy that offers a new regime in which to test models of dynamical heating. In spite of its proximity, the dynamical heating history of M33 has not yet been well constrained. In this work, we present the TREX Survey, the largest stellar spectroscopic survey across the disk of M33. We present the stellar disk kinematics as a function of age to study the past and ongoing dynamical heating of M33. We measure line of sight velocities for ~4,500 disk stars. Using a subset, we divide the stars into broad age bins using Hubble Space Telescope and Canada-France-Hawaii-Telescope photometric catalogs: massive main sequence stars and helium burning stars (~80 Myr), intermediate mass asymptotic branch stars (~1 Gyr), and low mass red giant branch stars (~4 Gyr). We compare the stellar disk dynamics to that of the gas using existing HI, CO, and Halpha kinematics. We find that the disk of M33 has relatively low velocity dispersion (~16 km/s), and unlike in the Milky Way and Andromeda galaxies, there is no strong trend in velocity dispersion as a function of stellar age. The youngest disk stars are as dynamically hot as the oldest disk stars and are dynamically hotter than predicted by most M33 like low mass simulated analogs in Illustris. The velocity dispersion of the young stars is highly structured, with the large velocity dispersion fairly localized. The cause of this high velocity dispersion is not evident from the observations and simulated analogs presented here. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 27 pages, 15 figures, 6 tables

Metallicity Distribution Function of the Eridanus~II Ultra-Faint Dwarf Galaxy from Hubble Space Telescope Narrow-band Imaging

Authors: Sal Wanying Fu, Daniel R. Weisz, Else Starkenburg, Nicolas Martin, Alexander P. Ji, Ekta Patel, Michael Boylan-Kolchin, Patrick Cote, Andrew E. Dolphin, Nicolas Longeard, Mario L. Mateo, Nathan R. Sandford

Abstract: We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more… ▽ More We use deep narrowband Ca H&K ($F395N$) imaging taken with the Hubble Space Telescope (HST) to construct the metallicity distribution function (MDF) of Local Group (LG) ultra-faint dwarf (UFD) galaxy Eridanus II (Eri II). When combined with archival $F475W$ and $F814W$ data, we measure metallicities for 60 resolved red giant branch stars as faint as $m_{F475W}\sim24$ mag, a factor of $\sim4$x more stars than current spectroscopic MDF determinations. We find that Eri II has a mean metallicity of [Fe/H]$=$-2.50$^{+0.07}_{-0.07}$ and a dispersion of $σ_{\mbox{[Fe/H]}}=0.42^{+0.06}_{-0.06}$, which are consistent with spectroscopic MDFs, though more precisely constrained owing to a larger sample. We identify a handful of extremely metal-poor star candidates (EMP; [Fe/H] $< -3$) that are marginally bright enough for spectroscopic follow up. Eri II's MDF appears well-described by a leaky box chemical evolution model. We also compute an updated orbital history for Eri II using Gaia eDR3 proper motions, and find that it is likely on first infall into the Milky Way. Our findings suggest that Eri II underwent an evolutionary history similar to that of an isolated galaxy. Compared to MDFs for select cosmological simulations of similar mass galaxies, we find that Eri II has a lower fraction of stars with [Fe/H] $< -3$, though such comparisons should currently be treated with caution due to a paucity of simulations, selection effects, and known limitations of Ca H&K for EMPs. This study demonstrates the power of deep HST CaHK imaging for measuring the MDFs of UFDs. △ Less

Submitted 29 October, 2021; originally announced November 2021.

Comments: Resubmitted to ApJ with revisions following a positive referee report of the initial draft; 29 pages, 19 figures, and 3 tables, including appendix. Comments welcome!

The Clustering of Orbital Poles Induced by the LMC: Hints for the Origin of Planes of Satellites

Authors: Nicolas Garavito-Camargo, Ekta Patel, Gurtina Besla, Adrian M. Price-Whelan, Facundo A. Gomez, Chervin F. P Laporte, Kathryn V. Johnston

Abstract: A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the… ▽ More A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter (DM) particles from the present-day random distribution. Instead, the orbital poles of particles beyond $R\approx 50$kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor $δρ_{max}$=30\%(50\%) with respect to underdense regions, and $δρ_{iso}$=15\%(30\%) relative to the isolated MW simulation (no LMC) between 50-150 kpc (150-300 kpc). The clustering appears after the LMC's pericenter ($\approx$ 50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: 1) the LMC shifts the velocity and position of the central density of the MW's halo and disk; 2) the DM dynamical friction wake and collective response induced by the LMC changes the kinematics of particles; 3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least $\approx$ 1:10 mass ratio), causing the clustering of orbital poles of halo tracers. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ. 23 pages and 16 figures. Comments are welcome

Star Formation Histories of Ultra-Faint Dwarf Galaxies: environmental differences between Magellanic and non-Magellanic satellites?

Authors: Elena Sacchi, Hannah Richstein, Nitya Kallivayalil, Roeland van der Marel, Mattia Libralato, Paul Zivick, Gurtina Besla, Thomas M. Brown, Yumi Choi, Alis Deason, Tobias Fritz, Marla Geha, Puragra Guhathakurta, Myoungwon Jeon, Evan Kirby, Steven R. Majewski, Ekta Patel, Joshua D. Simon, Sangmo Tony Sohn, Erik Tollerud, Andrew Wetzel

Abstract: We present the color-magnitude diagrams and star formation histories (SFHs) of seven ultra-faint dwarf galaxies: Horologium 1, Hydra 2, Phoenix 2, Reticulum 2, Sagittarius 2, Triangulum 2, and Tucana 2, derived from high-precision Hubble Space Telescope photometry. We find that the SFH of each galaxy is consistent with them having created at least 80% of the stellar mass by $z\sim6$. For all galax… ▽ More We present the color-magnitude diagrams and star formation histories (SFHs) of seven ultra-faint dwarf galaxies: Horologium 1, Hydra 2, Phoenix 2, Reticulum 2, Sagittarius 2, Triangulum 2, and Tucana 2, derived from high-precision Hubble Space Telescope photometry. We find that the SFH of each galaxy is consistent with them having created at least 80% of the stellar mass by $z\sim6$. For all galaxies, we find quenching times older than 11.5 Gyr ago, compatible with the scenario in which reionization suppresses the star formation of small dark matter halos. However, our analysis also reveals some differences in the SFHs of candidate Magellanic Cloud satellites, i.e., galaxies that are likely satellites of the Large Magellanic Cloud and that entered the Milky Way potential only recently. Indeed, Magellanic satellites show quenching times about 600 Myr more recent with respect to those of other Milky Way satellites, on average, even though the respective timings are still compatible within the errors. This finding is consistent with theoretical models that suggest that satellites' SFHs may depend on their host environment at early times, although we caution that within the error bars all galaxies in our sample are consistent with being quenched at a single epoch. △ Less

Submitted 28 September, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 7 pages, 3 figures, 2 tables. Accepted for publication in ApJL

HST Proper Motions of NGC 147 and NGC 185: Orbital Histories and Test of Dynamically Coherent Andromeda Satellite Plane

Authors: Sangmo Tony Sohn, Ekta Patel, Mark A. Fardal, Gurtina Besla, Roeland P. van der Marel, Marla Geha, Puragra Guhathakurta

Abstract: We present the first proper motion (PM) measurements for the dwarf elliptical galaxies NGC 147 and NGC 185, two satellite galaxies of M31, using multi-epoch HST imaging data with time baselines of $\sim 8$ years. For each galaxy, we take an error-weighted average of measurements from HST ACS/WFC and WFC3/UVIS to determine the PMs. Our final results for the PMs are… ▽ More We present the first proper motion (PM) measurements for the dwarf elliptical galaxies NGC 147 and NGC 185, two satellite galaxies of M31, using multi-epoch HST imaging data with time baselines of $\sim 8$ years. For each galaxy, we take an error-weighted average of measurements from HST ACS/WFC and WFC3/UVIS to determine the PMs. Our final results for the PMs are $(μ_\mathrm{W}, μ_\mathrm{N})_\mathrm{N147} = (-0.0232, 0.0378) \pm (0.0143, 0.0146)\,\mathrm{mas}\,\mathrm{yr}^{-1}$ for NGC 147, and $(μ_\mathrm{W}, μ_\mathrm{N})_\mathrm{N185} = (-0.0242, 0.0058) \pm (0.0141, 0.0147)\,\mathrm{mas}\,\mathrm{yr}^{-1}$ for NGC 185. The 2-dimensional direction of motion for NGC 147 about M31 is found to be aligned with its tidal tails. The 3-d positions and velocities of both galaxies are transformed into a common M31-centric coordinate system to study the detailed orbital histories of the combined M31+NGC 147+NGC 185 system via numerical orbit integration. We find that NGC 147 (NGC 185) had its closest passage to M31 0.3-0.5~Gyr ($\gtrsim 1.6$~Gyr) within the past 6 Gyr at distances of $\sim 70$ kpc (70-260 kpc). The pericentric times of NGC 147/NGC 185 correlate qualitatively well with the presence/absence of tidal tails seen around the galaxies. Our PMs show that the orbital poles of NGC 147, and also NGC 185 albeit to a lesser degree, agree within the uncertainties with the normal of the Great Plane of Andromeda (GPoA). These are the first measurements of the 3-d angular momentum vector of any satellite identified as original GPoA members. Our results strengthen the hypothesis that the GPoA may be a dynamically coherent entity. We revisit previous claims that NGC 147 and NGC 185 are binary galaxies and conclude that it is very unlikely the two galaxies were ever gravitationally bound to each other. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 20 pages, 9 figures, 4 tables, accepted for publication in ApJ

The Orbital Histories of Magellanic Satellites Using Gaia DR2 Proper Motions

Authors: Ekta Patel, Nitya Kallivayalil, Nicolas Garavito-Camargo, Gurtina Besla, Daniel R. Weisz, Roeland P. van der Marel, Michael Boylan-Kolchin, Marcel S. Pawlowski, Facundo A. Gómez

Abstract: With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint… ▽ More With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: i.) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); ii.) Magellanic satellites that were recently captured by the MCs $<$ 1 Gyr ago (Reticulum 2, Phoenix 2); iii.) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and iv.) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC's gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic satellites. △ Less

Submitted 1 March, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Comments: 32 pages, 11 figures. Accepted to ApJ

Asymmetric Drift of Andromeda Analogs in the IllustrisTNG Simulation

Authors: Amanda C. N. Quirk, Ekta Patel

Abstract: We analyze the kinematics as a function of stellar age for Andromeda (M31) mass analogs from the IllustrisTNG cosmological simulation. We divide the star particles into four age groups: less than 1 Gyr, 1 to 5 Gyr, 5 to 10 Gyr, and greater 10 Gyr, and compare the kinematics of these groups to that of the neutral gas cells. We calculate rotation curves for the stellar and gaseous components of each… ▽ More We analyze the kinematics as a function of stellar age for Andromeda (M31) mass analogs from the IllustrisTNG cosmological simulation. We divide the star particles into four age groups: less than 1 Gyr, 1 to 5 Gyr, 5 to 10 Gyr, and greater 10 Gyr, and compare the kinematics of these groups to that of the neutral gas cells. We calculate rotation curves for the stellar and gaseous components of each analog from 2 kpc to 20 kpc from the center of mass. We find that the lag, or asymmetric drift (AD), between the gas rotation curve and the stellar rotation curve on average increases with stellar age. This finding is consistent with observational measurements of AD in the disk of the Andromeda galaxy. When the M31 analogs are separated into groups based on merger history, we find that there is a difference in the AD of the analogs that have had a 4:1 merger the last 4 Gyr, 8 Gyr, or 12 Gyr compared to analogs that have not experienced a 4:1 merger in the same time frame. The subset of analogs that have had a 4:1 merger within the last 4 Gyr are also similar to AD measurements of stars in the disk of M31, providing evidence that M31 may in fact have recently merged with a galaxy nearly one fourth of its mass. Further work using high resolution zoom in simulations is required to explore the contribution of internal heating to AD. △ Less

Submitted 20 July, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: 14 pages, 11 figures, accepted to MNRAS for publication July 20, 2020

Astro2020 Science White Paper: Construction of an L* Galaxy: the Transformative Power of Wide Fields for Revealing the Past, Present and Future of the Great Andromeda System

Authors: Karoline M. Gilbert, Erik J. Tollerud, Jay Anderson, Rachael L. Beaton, Eric F. Bell, Alyson Brooks, Thomas M. Brown, James Bullock, Jeffrey L. Carlin, Michelle Collins, Andrew Cooper, Denija Crnojevic, Julianne Dalcanton, Andres del Pino, Richard D'Souza, Ivanna Escala, Mark Fardal, Andreea Font, Marla Geha, Puragra Guhathakurta, Evan Kirby, Geraint F. Lewis, Jennifer L. Marshall, Nicolas F. Martin, Kristen McQuinn , et al. (12 additional authors not shown)

Abstract: The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover,… ▽ More The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today. △ Less

Submitted 1 April, 2019; originally announced April 2019.

Comments: Submitted as a science white paper to the Astro2020 Decadal Survey

The Frequency of Dwarf Galaxy Multiples at Low Redshift in SDSS vs. Cosmological Expectations

Authors: Gurtina Besla, David R. Patton, Sabrina Stierwalt, Vicente Rodriguez-Gomez, Ekta Patel, Nitya J. Kallivayalil, Kelsey Johnson, Sarah Pearson, George Privon, Mary E. Putman

Abstract: We quantify the frequency of companions of low redshift ($0.013 < z < 0.0252$), dwarf galaxies ($2 \times 10^8$ M$_\odot <$ M$_{*} < 5 \times 10^9$ M$_\odot$) that are isolated from more massive galaxies in SDSS and compare against cosmological expectations using mock observations of the Illustris simulation. Dwarf multiples are defined as 2 or more dwarfs that have angular separations > 55'', pro… ▽ More We quantify the frequency of companions of low redshift ($0.013 < z < 0.0252$), dwarf galaxies ($2 \times 10^8$ M$_\odot <$ M$_{*} < 5 \times 10^9$ M$_\odot$) that are isolated from more massive galaxies in SDSS and compare against cosmological expectations using mock observations of the Illustris simulation. Dwarf multiples are defined as 2 or more dwarfs that have angular separations > 55'', projected separations r$_p < 150$ kpc and relative line-of-sight velocities $ΔV_{\rm LOS} < 150$ km/s. While the mock catalogs predict a factor of 2 more isolated dwarfs than observed in SDSS, the mean number of observed companions per dwarf is $N_c \sim 0.04$, in good agreement with Illustris when accounting for SDSS sensitivity limits. Removing these limits in the mock catalogs predicts $N_c\sim 0.06$ for future surveys (LSST, DESI), which will be complete to M$_* = 2\times 10^8$ M$_\odot$. The 3D separations of mock dwarf multiples reveal a contamination fraction of $\sim$40% in observations from projection effects. Most isolated multiples are pairs; triples are rare and it is cosmologically improbable that bound groups of dwarfs with more than 3 members exist within the parameter range probed in this study. We find that $<$1% of LMC-analogs in the field have an SMC-analog companion. The fraction of dwarf "Major Pairs'' (stellar mass ratio $>$1:4) steadily increases with decreasing Primary stellar mass, whereas the cosmological "Major Merger rate'' (per Gyr) has the opposite behaviour. We conclude that cosmological simulations can be reliably used to constrain the fraction of dwarf mergers across cosmic time. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: 21 pages, 16 Figures, submitted to MNRAS (Referee stage ; Comments are welcome)

$Λ$CDM Predictions for the Satellite Population of M33

Authors: Ekta Patel, Jeffrey L. Carlin, Erik J. Tollerud, Michelle L. M. Collins, Gregory A. Dooley

Abstract: 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}$. $Λ$CDM simulations predict that M33-mass halos host several of their own low mass satellite companions, yet only one candidate M33 satellite has been de… ▽ More 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}$. $Λ$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 $Λ$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. △ Less

Submitted 17 July, 2018; v1 submitted 13 July, 2018; originally announced July 2018.

Comments: 15 pages, 3 figures, submitted to MNRAS

First Gaia Dynamics of the Andromeda System: DR2 Proper Motions, Orbits, and Rotation of M31 and M33

Authors: Roeland P. van der Marel, Mark A. Fardal, Sangmo Tony Sohn, Ekta Patel, Gurtina Besla, Andrés del Pino-Molina, Johannes Sahlmann, Laura L. Watkins

Abstract: The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is… ▽ More The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is confidently detected, consistent with the known line-of-sight rotation curves: $V_{\rm rot} = -206\pm86$ km s$^{-1}$ (counter-clockwise) for M31, and $V_{\rm rot} = 80\pm52$ km s$^{-1}$ (clockwise) for M33. We measure the center-of-mass PM of each galaxy relative to surrounding background quasars in DR2. This yields that $(μ_{α*},μ_δ)$ equals $(65 \pm 18 , -57 \pm 15)$ $μ$as yr$^{-1}$ for M31, and $(31 \pm 19 , -29 \pm 16)$ $μ$as yr$^{-1}$ for M33. In addition to the listed random errors, each component has an additional residual systematic error of 16 $μ$as yr$^{-1}$. These results are consistent at 0.8$σ$ and 1.0$σ$ with the (2 and 3 times higher-accuracy) measurements already available from Hubble Space Telescope (HST) optical imaging and VLBA water maser observations, respectively. This lends confidence that all these measurements are robust. The new results imply that the M31 orbit towards the Milky Way is somewhat less radial than previously inferred, $V_{\rm tan, DR2+HST} = 57^{+35}_{-31}$ km s$^{-1}$, and strengthen arguments that M33 may be on its first infall into M31. The results highlight the future potential of Gaia for PM studies beyond the Milky Way satellite system. △ Less

Submitted 28 January, 2019; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 15 pages, 7 figures, ApJ, in press

Estimating the Mass of the Milky Way Using the Ensemble of Classical Satellite Galaxies

Authors: Ekta Patel, Gurtina Besla, Kaisey Mandel, Sangmo Tony Sohn

Abstract: High precision proper motion (PM) measurements are available for approximately 20% of all known dwarf satellite galaxies of the Milky Way (MW). Here we extend the Bayesian framework of Patel et al. (2017b) to include all MW satellites with measured 6D phase space information and apply it with the Illustris-Dark simulation to constrain the MW's mass. Using the properties of each MW satellite indivi… ▽ More High precision proper motion (PM) measurements are available for approximately 20% of all known dwarf satellite galaxies of the Milky Way (MW). Here we extend the Bayesian framework of Patel et al. (2017b) to include all MW satellites with measured 6D phase space information and apply it with the Illustris-Dark simulation to constrain the MW's mass. Using the properties of each MW satellite individually, we find that the scatter among mass estimates is reduced when the magnitude of specific orbital angular momentum (j) is adopted rather than their combined instantaneous positions and velocities. We also find that high j satellites (i.e. Leo II) constrain the upper limits for the MW's mass and low j satellites rather than the highest speed satellites (i.e. Leo I and LMC), set the lower mass limits. When j of all classical satellites is used to simultaneously estimate the MW's mass, we conclude the halo mass is $0.85^{+0.23}_{-0.26}\times 10^{12} \rm \, M_{\odot}$ (including Sagittarius dSph) and $0.96^{+0.29}_{-0.28} \times 10^{12} \rm \, M_{\odot}$ (excluding Sagittarius dSph), cautioning that low j satellites on decaying orbits like Sagittarius dSph may bias the distribution. These estimates markedly reduce the current factor of two spread in the mass range of the MW. We also find a well-defined relationship between host halo mass and satellite j distribution, which yields the prediction that upcoming PMs for ultra-faint dwarfs should reveal j within $ 5\times10^3 - 10^4$ kpc km s$^{-1}$. This is a promising method to significantly constrain the cosmologically expected mass range for the MW and eventually M31 as more satellite PMs become available. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 20 pages, 7 figures, submitted to ApJ

Space Motions of the Dwarf Spheroidal Galaxies Draco and Sculptor based on HST Proper Motions with ~10-year Time Baseline

Authors: Sangmo Tony Sohn, Ekta Patel, Gurtina Besla, Roeland P. van der Marel, James S. Bullock, Louis E. Strigari, Glenn van de Ven, Matthew G. Walker

Abstract: We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu… ▽ More We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu_W,mu_N)_Scl = (-0.0296+/-0.0209,-0.1358 +/-0.0214) mas/yr. The implied Galactocentric velocity vectors for Draco and Sculptor have radial and tangential components: (V_rad,V_tan)_Dra = (-88.6,161.4) +/- (4.4,5.6) km/s; and (V_rad,V_tan)_Scl = (72.6,200.2) +/- (1.3,10.8) km/s. We study the detailed orbital history of both Draco and Sculptor via numerical orbit integrations. Orbital periods of Draco and Sculptor are found to be 1-2 and 2-5 Gyrs, respectively, accounting for uncertainties in the MW mass. We also study the influence of the Large Magellanic Cloud (LMC) on the orbits of Draco and Sculptor. Overall, the inclusion of the LMC increases the scatter in the orbital results. Based on our calculations, Draco shows a rather wide range of orbital parameters depending on the MW mass and inclusion/exclusion of the LMC, but Sculptor's orbit is very well constrained with its most recent pericentric approach to the MW being 0.3-0.4 Gyr ago. Our new PMs imply that the orbital trajectories of both Draco and Sculptor are confined within the Disk of Satellites (DoS), better so than implied by earlier PM measurements, and likely rule out the possibility that these two galaxies were accreted together as part of a tightly bound group. △ Less

Submitted 3 October, 2017; v1 submitted 9 July, 2017; originally announced July 2017.

Comments: 17 pages, 8 figures, 6 tables. Accepted for publication in ApJ

Orbits of massive satellite galaxies: II. Bayesian Estimates of the Milky Way and Andromeda masses using high precision astrometry and cosmological simulations

Authors: Ekta Patel, Gurtina Besla, Kaisey Mandel

Abstract: Space observatories like the Hubble Space Telescope and Gaia are providing unprecedented 6D phase space information of satellite galaxies. Such measurements can shed light on the structure and assembly history of the Local Group, but improved statistical methods are needed to use them efficiently. Here we illustrate such a method using analogs of the Local Group's two most massive satellite galaxi… ▽ More Space observatories like the Hubble Space Telescope and Gaia are providing unprecedented 6D phase space information of satellite galaxies. Such measurements can shed light on the structure and assembly history of the Local Group, but improved statistical methods are needed to use them efficiently. Here we illustrate such a method using analogs of the Local Group's two most massive satellite galaxies, the Large Magellanic Cloud (LMC) and Triangulum (M33), from the Illustris dark-matter-only cosmological simulation. We use a Bayesian inference scheme combining measurements of positions, velocities, and specific orbital angular momenta ( j ) of the LMC/M33 with importance sampling of their simulated analogs to compute posterior estimates of the Milky Way (MW) and Andromeda's (M31) halo masses. We conclude the resulting host halo mass is more susceptible to bias when using measurements of the current position and velocity of satellites, especially when satellites are at short-lived phases of their orbits (i.e. at pericentre). Instead, the j value of a satellite is well-conserved over time and provides a more reliable constraint on host mass. The inferred virial mass of the MW (M31) using j of the LMC (M33) is $\rm M_{vir, MW}=1.02^{+0.77}_{-0.55}\times10^{12}\; M_{\odot}$ ($\rm M_{vir, M31}=1.37^{+1.39}_{-0.75}\times10^{12}\; M_{\odot}$). Choosing simulated analogs whose j values are consistent with the conventional picture of a previous (< 3 Gyr ago), close encounter (< 100 kpc) of M33 about M31 results in a very low virial mass for M31 ($\rm\sim\!10^{12}\; M_{\odot}$). This supports the new scenario put forth in Patel et al. (2017), wherein M33 is on its first passage about M31 or on a long period orbit. We conclude that this Bayesian inference scheme, utilising satellite j, is a promising method to reduce the current factor of two spread in the mass range of the MW and M31 moving forward. △ Less

Submitted 16 March, 2017; originally announced March 2017.

Comments: 24 pages, 9 figures. Submitted to MNRAS

Orbits of Massive Satellite Galaxies: I. A Close Look at the Large Magellanic Cloud and a New Orbital History for M33

Authors: Ekta Patel, Gurtina Besla, Tony Sohn

Abstract: The Milky Way (MW) and M31 both harbor massive satellite galaxies, the Large Magellanic Cloud (LMC) and M33, which may comprise up to 10 per cent of their host's total mass. Massive satellites can change the orbital barycentre of the host-satellite system by tens of kiloparsecs and are cosmologically expected to harbor dwarf satellite galaxies of their own. Assessing the impact of these effects de… ▽ More The Milky Way (MW) and M31 both harbor massive satellite galaxies, the Large Magellanic Cloud (LMC) and M33, which may comprise up to 10 per cent of their host's total mass. Massive satellites can change the orbital barycentre of the host-satellite system by tens of kiloparsecs and are cosmologically expected to harbor dwarf satellite galaxies of their own. Assessing the impact of these effects depends crucially on the orbital histories of the LMC and M33. Here, we revisit the dynamics of the MW-LMC system and present the first detailed analysis of the M31-M33 system utilizing high precision proper motions and statistics from the dark matter-only Illustris cosmological simulation. With the latest Hubble Space Telescope proper motion measurements of M31, we reliably constrain M33's interaction history with its host. In particular, like the LMC, M33 is either on its first passage (t_{inf} < 2 Gyr ago) or if M31 is massive (>=2x10^12 Msun), it is on a long period orbit of about 6 Gyr. Cosmological analogs of the LMC and M33 identified in Illustris support this picture and provide further insight about their host masses. We conclude that, cosmologically, massive satellites like the LMC and M33 are likely completing their first orbits about their hosts. We also find that the orbital energies of such analogs prefer a MW halo mass ~1.5x10^12 Msun and an M31 halo mass >=1.5x10^12 Msun. Despite conventional wisdom, we conclude it is highly improbable that M33 made a close (< 100 kpc) approach to M31 recently ( t_{peri} < 3 Gyr ago). Such orbits are rare (< 1 per cent) within the 4$σ$ error space allowed by observations. This conclusion cannot be explained by perturbative effects through four body encounters between the MW, M31, M33, and the LMC. This surprising result implies that we must search for a new explanation for M33's strongly warped gas and stellar discs. △ Less

Submitted 16 November, 2016; v1 submitted 15 September, 2016; originally announced September 2016.

Comments: 28 pages, 17 figures, 8 tables, accepted to MNRAS Oct 7 2016