Showing 1–50 of 74 results for author: Fardal, M

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

BP3M: Bayesian Positions, Parallaxes, and Proper Motions derived from the Hubble Space Telescope and Gaia data

Authors: Kevin A. McKinnon, Andrés del Pino, Constance M. Rockosi, Miranda Apfel, Puragra Guhathakurta, Roeland P. van der Marel, Paul Bennet, Mark A. Fardal, Mattia Libralato, Sangmo Tony Sohn, Eduardo Vitral, Laura L. Watkins

Abstract: We present a hierarchical Bayesian pipeline, BP3M, that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble~Space~Telescope (HST) images and Gaia -- even for sparse fields ($N_*<10$ per image) -- expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is there… ▽ More We present a hierarchical Bayesian pipeline, BP3M, that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble~Space~Telescope (HST) images and Gaia -- even for sparse fields ($N_*<10$ per image) -- expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is therefore able to put the HST images onto a global reference frame without the use of background galaxies/QSOs. Testing our publicly-available code in the Fornax and Draco dSphs, we measure accurate PMs that are a median of 8-13 times more precise than Gaia DR3 alone for $20.5<G<21~\mathrm{mag}$. We are able to explore the effect of observation strategies on BP3M astrometry using synthetic data, finding an optimal strategy to improve parallax and position precision at no cost to the PM uncertainty. Using 1619 HST images in the sparse COSMOS field (median 9 Gaia sources per HST image), we measure BP3M PMs for 2640 unique sources in the $16<G<21.5~\mathrm{mag}$ range, 25% of which have no Gaia PMs; the median BP3M PM uncertainty for $20.25<G<20.75~\mathrm{mag}$ sources is $0.44~$mas/yr compared to $1.03~$mas/yr from Gaia, while the median BP3M PM uncertainty for sources without Gaia-measured PMs ($20.75<G<21.5~\mathrm{mag}$) is $1.16~$mas/yr. The statistics that underpin the BP3M pipeline are a generalized way of combining position measurements from different images, epochs, and telescopes, which allows information to be shared between surveys and archives to achieve higher astrometric precision than that from each catalog alone. △ Less

Submitted 4 November, 2023; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 33 pages, 25 figures, 3 tables

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

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

GaiaHub: A method for combining data from the Gaia and Hubble space telescopes to derive improved proper motions for faint stars

Authors: Andrés del Pino, Mattia Libralato, Roeland P. van der Marel, Paul Bennet, Mark A. Fardal, Jay Anderson, Andrea Bellini, Sangmo Tony Sohn, Laura L. Watkins

Abstract: We present GaiaHub, a publicly available tool that combines $Gaia$ measurements with $Hubble$ $Space$ $Telescope$ ($HST$) archival images to derive proper motions (PMs). It increases the scientific impact of both observatories beyond their individual capabilities. $Gaia$ provides PMs across the whole sky, but the limited mirror size and time baseline restrict the best PM performance to relatively… ▽ More We present GaiaHub, a publicly available tool that combines $Gaia$ measurements with $Hubble$ $Space$ $Telescope$ ($HST$) archival images to derive proper motions (PMs). It increases the scientific impact of both observatories beyond their individual capabilities. $Gaia$ provides PMs across the whole sky, but the limited mirror size and time baseline restrict the best PM performance to relatively bright stars. $HST$ can measure accurate PMs for much fainter stars over a small field, but this requires two epochs of observation which are not always available. GaiaHub yields considerably improved PM accuracy compared to $Gaia$-only measurements, especially for faint sources $(G \gtrsim 18)$, requiring only a single epoch of $HST$ data observed more than $\sim 7$ years ago (before 2012). This provides considerable scientific value especially for dynamical studies of stellar systems or structures in and beyond the Milky Way (MW) halo, for which the member stars are generally faint. To illustrate the capabilities and demonstrate the accuracy of GaiaHub, we apply it to samples of MW globular clusters (GCs) and classical dwarf spheroidal (dSph) satellite galaxies. This allows us, e.g., to measure the velocity dispersions in the plane of the sky for objects out to and beyond $\sim 100$ kpc. We find, on average, mild radial velocity anisotropy in GCs, consistent with existing results for more nearby samples. We observe a correlation between the internal kinematics of the clusters and their ellipticity, with more isotropic clusters being, on average, more round. Our results also support previous findings that Draco and Sculptor dSph galaxies appear to be radially anisotropic systems. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 24 pages, 18 figures. Accepted for publication in ApJ

Kinematics and Metallicity of Red Giant Branch Stars in the Northeast Shelf of M31

Authors: Ivanna Escala, Karoline M. Gilbert, Mark Fardal, Puragra Guhathakurta, Robyn E. Sanderson, Jason S. Kalirai, Bahram Mobasher

Abstract: We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge… ▽ More We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in 4 spectroscopic fields spanning $13-31$ projected kpc along the Northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge pattern agrees with expectations of a tidal shell formed in a radial merger and provides strong evidence in favor of predictions of Giant Stellar Stream (GSS) formation models in which the NE shelf originates from the second orbital wrap of the tidal debris. The observed concentric wedge patterns of the NE, West (W), and Southeast (SE) shelves corroborate this interpretation independently of the models. We do not detect a kinematical signature in the NE shelf region corresponding to an intact progenitor core, favoring GSS formation models in which the progenitor is completely disrupted. The shelf's photometric metallicity distribution implies that it is dominated by tidal material, as opposed to the phase-mixed stellar halo or the disk. The metallicity distribution ([Fe/H]$_{\rm phot}$ = $-0.42$ $\pm$ $0.01$) also matches the GSS, and consequently the W and SE shelves, further supporting a direct physical association between the tidal features. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 21 pages main text, 15 figures, 5 tables (including Appendix). Submitted to AJ

The nature of the Eastern Extent in the outer halo of M31

Authors: Janet Preston, Michelle Collins, R. Michael Rich, Rodrigo Ibata, Nicolas F. Martin, Mark Fardal

Abstract: We present the first comprehensive spectroscopic study of the Andromeda galaxy's Eastern Extent. This ~4 degree long filamentary structure. located 70-90 kpc from the centre of M31, lies perpendicular to Andromeda's minor axis and the Giant Stellar Stream and overlaps Stream C. In this work, we explore the properties of the Eastern Extent to look for possible connections between it, the Giant Stel… ▽ More We present the first comprehensive spectroscopic study of the Andromeda galaxy's Eastern Extent. This ~4 degree long filamentary structure. located 70-90 kpc from the centre of M31, lies perpendicular to Andromeda's minor axis and the Giant Stellar Stream and overlaps Stream C. In this work, we explore the properties of the Eastern Extent to look for possible connections between it, the Giant Stellar Stream and Stream C. We present the kinematics and photometry for ~50 red giant branch stars in 7 fields along the Eastern Extent. We measure the systemic velocities for these fields and find them to be -368 km/s < v < -331 km/s with a slight velocity gradient of -0.51 +/- 0.21 km/s/kpc towards the Giant Stellar Stream. We derive the photometric metallicities for stars in the Eastern Extent finding them to be metal-poor with values of -1.0 < [Fe/H]phot < -0.7 with a <[Fe/H]phot> ~-0.9. We find consistent properties for the Eastern Extent, Stream B and one of the substructures in Stream C, Stream Cr, plausibly linking these features. Stream Cp and its associated globular cluster, EC4, have distinctly different properties indicative of a separate structure. When we compare the properties of the Eastern Extent to those of the Giant Stellar Stream, we find them to be consistent, albeit slightly more metal-poor, such that the Eastern Extent could plausibly comprise stars stripped from the progenitor of the Giant Stellar Stream. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Mapping Gaia parallax systematic errors over the sky with faint Milky Way stars

Authors: Mark A. Fardal, Roeland van der Marel, Andres del Pino, Sangmo Tony Sohn

Abstract: Parallaxes measured by the Gaia mission have huge significance for astronomy, but parallaxes in Gaia DR2 are known to have systematic errors that depend on the source position and other quantities. We use the abundant information in faint Milky Way stars, along with the GOG simulation of the Gaia catalog, to probe the spatial dependence of Gaia DR2 parallax systematic errors in an empirical way. T… ▽ More Parallaxes measured by the Gaia mission have huge significance for astronomy, but parallaxes in Gaia DR2 are known to have systematic errors that depend on the source position and other quantities. We use the abundant information in faint Milky Way stars, along with the GOG simulation of the Gaia catalog, to probe the spatial dependence of Gaia DR2 parallax systematic errors in an empirical way. The parallax signal, concentrated in thick disk turnoff stars with magnitude G ~ 17, is sufficient to construct maps of the parallax systematic error over the majority of the sky. These maps show a locally regular "waffle pattern" on ~1 degree scales following Gaia scan directions, stronger linear "scar" features, and coherent variations on larger scales. The parallax bias maps also retain traces of astrophysical effects such as dust clouds. The waffle pattern, known from earlier maps of the Magellanic Clouds, extends over the entire sky; its local rms amplitude averages 15 microarcsec and varies by about a factor of two. The strength of this pattern increases by a factor ~6 from magnitude G = 13 to G = 20. Correlations with parallaxes of quasars and of stars with independent distance estimates support our bias estimates. Using similar methods, we map systematic errors in the proper motion and examine the relationship with the parallax systematics. We provide a code package to access and query our bias maps. Similar tests on the general stellar population should be useful in quantifying systematic errors in future Gaia releases. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: Accepted in AJ. 17 pages, 11 figures. Higher resolution figures available in journal. Eagerly anticipating Gaia EDR3

Revealing the Structure and Internal Rotation of the Sagittarius Dwarf Spheroidal Galaxy with Gaia and Machine Learning

Authors: Andrés del Pino, Mark A. Fardal, Roeland P. van der Marel, Ewa L. Łokas, Cecilia Mateu, Sangmo Tony Sohn

Abstract: We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than… ▽ More We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than $1.2\times10^5$ member stars in the core of the galaxy. Our results show that Sgr has a bar structure $\sim 2.5$ kpc long, and that tidal tails emerge from its tips to form what it is known as the Sgr stream. The main body of the galaxy, strongly sheared by tidal forces, is a triaxial (almost prolate) ellipsoid with its longest principal axis of inertia inclined $43^\circ\pm6^\circ$ with respect to the plane of the sky and axis ratios of 1:0.67:0.60. Its external regions are expanding mainly along its longest principal axis, yet the galaxy conserves an inner core of about $500\times330\times300$ pc that shows no net expansion and is rotating at $v_{\rm rot} = 4.13 \pm 0.16$ ${\rm{ km \ s^{-1}}}$. The internal angular momentum of Sgr forms an angle $θ= 18^\circ\pm6^\circ$ with respect to its orbital angular momentum, meaning that the galaxy is in an inclined prograde orbit around the Milky Way. We compared our results with predictions from $N$-body models with spherical, pressure-supported progenitors and a model whose progenitor is a flattened rotating disk. Only the rotating model, based on preexisting simulations aimed at reproducing the line-of-sight velocity gradients observed in Sgr, was able to reproduce the observed properties in the core of the galaxy. △ Less

Submitted 15 February, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 37 pages, 26 figures, 2 tables. Accepted for publication in the Astrophysical Journal

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 absolute proper motions of the Arches and Quintuplet clusters

Authors: M. Libralato, M. Fardal, D. Lennon, R. P. van der Marel, A. Bellini

Abstract: Arches and Quintuplet are two young, massive clusters projected near the Galactic Center. To date, studies focused on understanding their origin have been based on proper motions (PMs) derived in the clusters' reference frames and required some assumptions about their 3D motion. In this paper, we combine public PM catalogs of these clusters with the Gaia DR2 catalog and, for the first time, transf… ▽ More Arches and Quintuplet are two young, massive clusters projected near the Galactic Center. To date, studies focused on understanding their origin have been based on proper motions (PMs) derived in the clusters' reference frames and required some assumptions about their 3D motion. In this paper, we combine public PM catalogs of these clusters with the Gaia DR2 catalog and, for the first time, transform the relative PMs of the Arches and Quintuplet clusters onto an absolute reference system. We find that the absolute PM of the Arches is $(μ_α\cosδ,μ_δ)$ $=$ $(-1.45\ \pm\ 0.23,-2.68\ \pm\ 0.14)$ mas yr$^{-1}$, and that of the Quintuplet is $(μ_α\cosδ,μ_δ)$ $=$ $(-1.19\ \pm\ 0.09,-2.66\ \pm \ 0.18)$ mas yr$^{-1}$. These values suggest that these systems are moving almost parallel to the Galactic plane. A measurement of the clusters' distances is still required to meaningfully constrain the clusters' orbits and shed light on the origin of the Arches and Quintuplet. △ Less

Submitted 19 August, 2020; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: 9 pages, 8 figures, accepted for publication in MNRAS. v2: Fixed Bibliography

The Impact of Wind Scalings on Stellar Growth and the Baryon Cycle in Cosmological Simulations

Authors: Shuiyao Huang, Neal Katz, Romeel Davé, Benjamin D. Oppenheimer, David H. Weinberg, Mark Fardal, Juna A. Kollmeier, Molly S. Peeples

Abstract: Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stella… ▽ More Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $η$ with velocity dispersion $σ$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $η\sim σ^5$ in low mass haloes with $σ< 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique. △ Less

Submitted 9 January, 2020; v1 submitted 17 October, 2019; originally announced October 2019.

Comments: 31 pages, 19 figures, accepted by MNRAS

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 course of the Orphan Stream in the Northern Galactic Hemisphere traced with Gaia DR2

Authors: Mark A. Fardal, Roeland P. van der Marel, Sangmo Tony Sohn, Andres del Pino Molina

Abstract: The Orphan Stream is one of the most prominent tidal streams in the Galactic halo. Using data on red giants, RR Lyrae, and horizontal branch stars from Gaia and other surveys, we determine the proper motion of the Orphan Stream over a path of more than 90 degrees on the sky. We also provide updated tracks for the sky position, distance, and radial velocity of the stream. Our tracks in these latter… ▽ More The Orphan Stream is one of the most prominent tidal streams in the Galactic halo. Using data on red giants, RR Lyrae, and horizontal branch stars from Gaia and other surveys, we determine the proper motion of the Orphan Stream over a path of more than 90 degrees on the sky. We also provide updated tracks for the sky position, distance, and radial velocity of the stream. Our tracks in these latter dimensions mostly agree with previous results. However, there are significant corrections to the earlier distance and latitude tracks as the stream approaches the Galactic disk. Stream stars selected with three-dimensional kinematics display a very tight red giant sequence. Concordantly, we find that applying a proper motion cut removes the most metal-rich stars from earlier spectroscopic samples of stream stars, though a significant dispersion remains indicating a dwarf galaxy origin. The deceleration of the stream towards its leading end suggests a circular velocity of ~200 km/s at a galactocentric radius ~30 kpc, consistent with other independent evidence. However, the track of the stream departs significantly from an orbit; the spatial track does not point along the same direction as the velocity vector, and it exhibits a lateral wiggle that is unlikely to match any reasonable orbit. The low metallicity and small dispersion of the stream in the various coordinates point to a progenitor with a relatively low dynamical mass ~$10^8$ Msun. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 14 pages, 10 figures, 1 table; submitted to MNRAS

MSC Class: The course of the Orphan Stream in the Northern Galactic Hemisphere traced with Gaia DR2

The Robustness of Cosmological Hydrodynamic Simulation Predictions to Changes in Numerics and Cooling Physics

Authors: Shuiyao Huang, Neal Katz, Romeel Davé, Mark Fardal, Juna Kollmeier, Benjamin D. Oppenheimer, Molly S. Peeples, Shawn Roberts, David H. Weinberg, Philip F. Hopkins, Robert Thompson

Abstract: We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The P… ▽ More We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$α$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII. △ Less

Submitted 9 January, 2020; v1 submitted 30 October, 2018; originally announced October 2018.

Comments: 29 pages, 25 figures, accepted by MNRAS

The internal rotation of globular clusters revealed by Gaia DR2

Authors: P. Bianchini, R. P. van der Marel, A. del Pino, L. L. Watkins, A. Bellini, M. A. Fardal, M. Libralato, A. Sills

Abstract: Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal ro… ▽ More Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal rotation on the plane-of-the-sky for a large sample of globular clusters using proper motions from Gaia DR2. We detect signatures of rotation in the tangential component of proper motions for 11 out of 51 clusters at a $>$3-sigma confidence level, confirming the detection reported in Gaia collaboration et al. (2018) for 8 clusters, and additionally identify 11 GCs with a 2-sigma rotation detection. For the clusters with a detected global rotation, we construct the two-dimensional rotation maps and proper motion rotation curves, and we assess the relevance of rotation with respect to random motions ($V/σ\sim0.08-0.51$). We find evidence of a correlation between the degree of internal rotation and relaxation time, highlighting the importance of long-term dynamical evolution in shaping the clusters current properties. This is a strong indication that angular momentum must have played a fundamental role in the earliest phases of cluster formation. Finally, exploiting the spatial information of the rotation maps and a comparison with line-of-sight data, we provide an estimate of the inclination of the rotation axis for a subset of 8 clusters. Our work demonstrates the potential of Gaia data for internal kinematic studies of globular clusters and provides the first step to reconstruct their intrinsic three-dimensional structure. △ Less

Submitted 30 August, 2018; v1 submitted 7 June, 2018; originally announced June 2018.

Comments: Accepted for publication in 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

Connecting the Milky Way potential profile to the orbital timescales and spatial structure of the Sagittarius Stream

Authors: Mark A. Fardal, Roeland P. van der Marel, David R. Law, Sangmo Tony Sohn, Branimir Sesar, Nina Hernitschek, Hans-Walter Rix

Abstract: Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic p… ▽ More Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic potentials. We find that debris at the apocenters must be dynamically young, in the sense of being stripped off in the last two pericentric passages, while the Sagittarius dwarf is currently experiencing a third passage. The ratio of apocenters is sensitive to both dynamical friction and the outer slope of the Galactic rotation curve. These dependences can be understood with simple regularities connecting the apocentric radii, circular velocities, and orbital period of the progenitor. The effect of dynamical friction on the stream can be constrained using substructure within the leading apocenter. Our ensembles of models are not intended as statistically proper fits to the stream. Nevertheless, out of the range of models we consider, we consistently find the mass within 100 kpc to be $\sim 7 \times 10^{11} \, M_{\odot}$, with a nearly flat rotation curve between 50 and 100 kpc. This points to a more extended Galactic halo than assumed in some current models. As in previous work, we find prolate or triaxial halos ease agreement with the track of the leading stream. We display the behavior of our models in various observational spaces and characterize the substructure expected within the stream. In particular, the young trailing stream visible near trailing apocenter should exhibit a tight trend of velocity with distance separate from the older debris, and we suggest that this will serve as an especially useful probe of the outer Galactic potential. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Submitted to MNRAS

Absolute HST Proper Motion (HSTPROMO) of Distant Milky Way Globular Clusters: Galactocentric Space Velocities and the Milky Way Mass

Authors: Sangmo Tony Sohn, Laura L. Watkins, Mark A. Fardal, Roeland P. van der Marel, Alis J. Deason, Gurtina Besla, Andrea Bellini

Abstract: We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing li… ▽ More We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{\rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{\rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{\rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $β= 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (<39.5 \rm{kpc}) = 0.61^{+0.18}_{-0.12} \times 10^{12}$ M$_{\odot}$, and the virial mass $M_\rm{vir} = 2.05^{+0.97}_{-0.79} \times 10^{12}$ M$_{\odot}$, are consistent with, but on the high side among recent mass estimates in the literature. △ Less

Submitted 14 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 20 pages, 10 figures, accepted for publication in ApJ

SMHASH: Anatomy of the Orphan Stream using RR Lyrae stars

Authors: David Hendel, Victoria Scowcroft, Kathryn V. Johnston, Mark A. Fardal, Roeland P. van der Marel, Sangmo Tony Sohn, Adrian M. Price-Whelan, Rachael L. Beaton, Gurtina Besla, Giuseppe Bono, Maria-Rosa L. Cioni, Gisella Clementini, Judith G. Cohen, Michele Fabrizio, Wendy L. Freedman, Alessia Garofalo, Carl J. Grillmair, Nitya Kallivayalil, Juna A. Kollmeier, David R. Law, Barry F. Madore, Steven R. Majewski, Massimo Marengo, Andrew J. Monson, Jillian R. Neeley , et al. (8 additional authors not shown)

Abstract: Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc… ▽ More Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 μm}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Comments: 18 pages, 10 figures. Submitted to MNRAS; comments welcome

The $>100$ kpc Distant Spur of the Sagittarius Stream and the Outer Virgo Overdensity, as Seen in PS1 RR Lyrae Stars

Authors: Branimir Sesar, Nina Hernitschek, Marion I. P. Dierickx, Mark A. Fardal, Hans-Walter Rix

Abstract: We report the detection of spatially distinct stellar density features near the apocenters of the Sagittarius (Sgr) stream's main leading and trailing arm. These features are clearly visible in a high-fidelity stellar halo map that is based on RR Lyrae from Pan-STARRS1: there is a plume of stars 10 kpc beyond the apocenter of the leading arm, and there is a "spur" extending to 130 kpc, almost 30 k… ▽ More We report the detection of spatially distinct stellar density features near the apocenters of the Sagittarius (Sgr) stream's main leading and trailing arm. These features are clearly visible in a high-fidelity stellar halo map that is based on RR Lyrae from Pan-STARRS1: there is a plume of stars 10 kpc beyond the apocenter of the leading arm, and there is a "spur" extending to 130 kpc, almost 30 kpc beyond the previously detected apocenter of the trailing arm. Such apocenter substructure is qualitatively expected in any Sgr stream model, as stars stripped from the progenitor at different pericenter passages become spatially separated there. The morphology of these new Sgr stream substructures could provide much-needed new clues and constraints for modeling the Sgr system, including the level of dynamical friction that Sgr has experienced. We also report the discovery of a new, presumably unrelated halo substructure at 80 kpc from the Sun and $10^\circ$ from the Sgr orbital plane, which we dub the Outer Virgo Overdensity. △ Less

Submitted 30 June, 2017; originally announced June 2017.

Comments: Accepted to ApJ Letters, this is Branimir Sesar's last paper in astronomy and he would like to thank his friends and colleagues for supporting him during the last 15 years

The PAndAS view of the Andromeda satellite system - II. Detailed properties of 23 M31 dwarf spheroidal galaxies

Authors: Nicolas F. Martin, Rodrigo A. Ibata, Geraint F. Lewis, Alan McConnachie, Arif Babul, Nicholas F. Bate, Edouard Bernard, Scott C. Chapman, Michelle M. L. Collins, Anthony R. Conn, Denija Crnojević, Mark A. Fardal, Annette M. N. Ferguson, Michael Irwin, A. Dougal Mackey, Brendan McMonigal, Julio F. Navarro, R. Michael Rich

Abstract: We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda's known satellite dwarf galaxies and cover a wide range in luminosity ($-11.6<M_V<-5.8$ or $10^{4.2}< L <10^{6.5} L_\odot$) and surface… ▽ More We present a comprehensive analysis of the structural properties and luminosities of the 23 dwarf spheroidal galaxies that fall within the footprint of the Pan-Andromeda Archaeological Survey (PAndAS). These dwarf galaxies represent the large majority of Andromeda's known satellite dwarf galaxies and cover a wide range in luminosity ($-11.6<M_V<-5.8$ or $10^{4.2}< L <10^{6.5} L_\odot$) and surface brightness ($25.1<μ_0<29.3$ mag/arcsec$^2$). We confirm most previous measurements, but find And XIX to be significantly larger than before ($r_h=3065^{+1065}_{-935}$ pc, $M_V=-10.1^{+0.8}_{-0.4}$) and cannot derive parameters for And XXVII as it is likely not a bound stellar system. We also significantly revise downward the luminosities of And~XV and And~XVI, which are now $M_V\sim-7.5$ or $L\sim10^5 L_\odot$. Finally, we provide the first detailed analysis of Cas II/And XXX, a fairly faint system ($M_V=-8.0^{+0.4}_{-0.3}$) of typical size ($r_h=270\pm50$ pc), located in close proximity to the two bright elliptical dwarf galaxies NGC 147 & 185. Combined with the set of homogeneous distances published in an earlier contribution, our analysis dutifully tracks all relevant sources of uncertainty in the determination of the properties of the dwarf galaxies from the PAndAS photometric catalogue. We further publish the posterior probability distribution functions of all the parameters we fit for in the form of MCMC chains available online; these inputs should be used in any analysis that aims to remain truthful to the data and properly account for covariance between parameters. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 30 pages, 39 figures. Accepted for publication in ApJ

The Growth and Enrichment of the Intragroup Gas

Authors: Lichen Liang, Fabrice Durier, Arif Babul, Romeel Davé, Benjamin D. Oppenheimer, Neal Katz, Mark Fardal, Tom Quinn

Abstract: The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/super… ▽ More The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback. △ Less

Submitted 3 December, 2015; originally announced December 2015.

Comments: 26 pages, 16 figures, accepted by MNRAS

Generation of mock tidal streams

Authors: Mark A. Fardal, Shuiyao Huang, Martin D. Weinberg

Abstract: In this paper we discuss a method for the generation of mock tidal streams. Using an ensemble of simulations in an isochrone potential where the actions and frequencies are known, we derive an empirical recipe for the evolving satellite mass and the corresponding mass loss rate, and the ejection conditions of the stream material. The resulting stream can then be quickly generated either with direc… ▽ More In this paper we discuss a method for the generation of mock tidal streams. Using an ensemble of simulations in an isochrone potential where the actions and frequencies are known, we derive an empirical recipe for the evolving satellite mass and the corresponding mass loss rate, and the ejection conditions of the stream material. The resulting stream can then be quickly generated either with direct orbital integration, or by using the action-angle formalism. The model naturally produces streaky features within the stream. These are formed due to the radial oscillation of the progenitor and the bursts of stars emitted near pericenter, rather than clumping at particular oscillation phases as sometimes suggested. When detectable, these streaky features are a reliable diagnostic for the stream's direction of motion and encode other information on the progenitor and its orbit. We show several tests of the recipe in alternate potentials, including a case with a chaotic progenitor orbit which displays a marked effect on the width of the stream. Although the specific ejection recipe may need adjusting when elements such as the orbit or satellite density profile are changed significantly, our examples suggest that model tidal streams can be quickly and accurately generated by models of this general type for use in Bayesian sampling. △ Less

Submitted 7 October, 2014; originally announced October 2014.

Comments: 20 pages, 11 figures, 1 table; submitted to MNRAS

A PAndAS view of M31 dwarf elliptical satellites: NGC147 and NGC185

Authors: D. Crnojević, A. M. N. Ferguson, M. J. Irwin, A. W. McConnachie, E. J. Bernard, M. A. Fardal, R. A. Ibata, G. F. Lewis, N. F. Martin, J. F. Navarro, N. E. D. Noël, S. Pasetto

Abstract: We exploit data from the Pan-Andromeda Archaeological Survey (PAndAS) to study the extended structures of M31's dwarf elliptical companions, NGC147 and NGC185. Our wide-field, homogeneous photometry allows to construct deep colour-magnitude diagrams (CMDs) which reach down to $\sim3$ mag below the red giant branch (RGB) tip. We trace the stellar components of the galaxies to surface brightness of… ▽ More We exploit data from the Pan-Andromeda Archaeological Survey (PAndAS) to study the extended structures of M31's dwarf elliptical companions, NGC147 and NGC185. Our wide-field, homogeneous photometry allows to construct deep colour-magnitude diagrams (CMDs) which reach down to $\sim3$ mag below the red giant branch (RGB) tip. We trace the stellar components of the galaxies to surface brightness of $μ_g \sim 32$ mag arcsec$^{-2}$ and show they have much larger extents ($\sim5$ kpc radii) than previously recognised. While NGC185 retains a regular shape in its peripheral regions, NGC147 exhibits pronounced isophotal twisting due to the emergence of symmetric tidal tails. We fit single Sersic models to composite surface brightness profiles constructed from diffuse light and star counts and find that NGC147 has an effective radius almost 3 times that of NGC185. In both cases, the effective radii that we calculate are larger by a factor of $\sim2$ compared to most literature values. We also calculate revised total magnitudes of $M_g=-15.36\pm0.04$ for NGC185 and $M_g=-16.36\pm0.04$ for NGC147. Using photometric metallicities computed for RGB stars, we find NGC185 to exhibit a metallicity gradient of [Fe/H]$\sim-0.15$ dex/kpc over the radial range 0.125 to 0.5 deg. On the other hand, NGC147 exhibits almost no metallicity gradient, $\sim-0.02$ dex/kpc from 0.2 to 0.6 deg. The differences in the structure and stellar populations in the outskirts of these systems suggest that tidal influences have played an important role in governing the evolution of NGC147. △ Less

Submitted 24 September, 2014; originally announced September 2014.

Comments: 17 pages, 13 figures; MNRAS accepted

The outer halo globular cluster system of M31 - II. Kinematics

Authors: J. Veljanoski, A. D. Mackey, A. M. N. Ferguson, A. P. Huxor, P. Cote, M. J. Irwin, N. R. Tanvir, J. Penarrubia, E. J. Bernard, M. Fardal, N. F. Martin, A. McConnachie, G. F. Lewis, S. C. Chapman, R. A. Ibata, A. Babul

Abstract: We present a detailed kinematic analysis of the outer halo globular cluster (GC) system of M31. Our basis for this is a set of new spectroscopic observations for 78 clusters lying at projected distances between Rproj ~20-140 kpc from the M31 centre. These are largely drawn from the recent PAndAS globular cluster catalogue; 63 of our targets have no previous velocity data. Via a Bayesian maximum li… ▽ More We present a detailed kinematic analysis of the outer halo globular cluster (GC) system of M31. Our basis for this is a set of new spectroscopic observations for 78 clusters lying at projected distances between Rproj ~20-140 kpc from the M31 centre. These are largely drawn from the recent PAndAS globular cluster catalogue; 63 of our targets have no previous velocity data. Via a Bayesian maximum likelihood analysis we find that GCs with Rproj > 30 kpc exhibit coherent rotation around the minor optical axis of M31, in the same direction as more centrally- located GCs, but with a smaller amplitude of 86+/-17 km s-1. There is also evidence that the velocity dispersion of the outer halo GC system decreases as a function of projected distance from the M31 centre, and that this relation can be well described by a power law of index ~ -0.5. The velocity dispersion profile of the outer halo GCs is quite similar to that of the halo stars, at least out to the radius up to which there is available information on the stellar kinematics. We detect and discuss various velocity correlations amongst subgroups of GCs that lie on stellar debris streams in the M31 halo. Many of these subgroups are dynamically cold, exhibiting internal velocity dispersions consistent with zero. Simple Monte Carlo experiments imply that such configurations are unlikely to form by chance, adding weight to the notion that a significant fraction of the outer halo GCs in M31 have been accreted alongside their parent dwarf galaxies. We also estimate the M31 mass within 200 kpc via the Tracer Mass Estimator, finding (1.2 - 1.6) +/- 0.2 10^{12}M_sun. This quantity is subject to additional systematic effects due to various limitations of the data, and assumptions built in into the TME. Finally, we discuss our results in the context of formation scenarios for the M31 halo. △ Less

Submitted 1 June, 2014; originally announced June 2014.

Comments: 24 pages, 12 figures, 7 tables; Accepted for publication in MNRAS

The Photon Underproduction Crisis

Authors: Juna A. Kollmeier, David H. Weinberg, Benjamin D. Oppenheimer, Francesco Haardt, Neal Katz, Romeel A. Davé, Mark Fardal, Piero Madau, Charles Danforth, Amanda B. Ford, Molly S. Peeples, Joseph McEwen

Abstract: We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed pro… ▽ More We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision. △ Less

Submitted 10 April, 2014; originally announced April 2014.

Comments: Submitted to ApJ Letters; 6 pages including 3 figures

Parametrising Star Formation Histories

Authors: Vimal Simha, David H. Weinberg, Charlie Conroy, Romeel Dave, Mark Fardal, Neal Katz, Benjamin D. Oppenheimer

Abstract: We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/ta… ▽ More We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies. △ Less

Submitted 1 April, 2014; originally announced April 2014.

Comments: 28 pages, 18 figures

The PAndAS Field of Streams: stellar structures in the Milky Way halo toward Andromeda and Triangulum

Authors: Nicolas F. Martin, Rodrigo A. Ibata, R. Michael Rich, Michelle L. M. Collins, Mark A. Fardal, Michael J. Irwin, Geraint F. Lewis, Alan W. McConnachie, Arif Babul, Nicholas F. Bate, Scott C. Chapman, Anthony R. Conn, Denija Crnojević, Annette M. N. Ferguson, A. Dougal Mackey, Julio F. Navarro, Jorge Peñarrubia, Nial T. Tanvir, David Valls-Gabaud

Abstract: We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defin… ▽ More We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of Σ_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances. △ Less

Submitted 19 March, 2014; originally announced March 2014.

Comments: 11 pages, 8 figures, accepted for publication in the ApJ, Figure 3 is the money plot

The large-scale structure of the halo of the Andromeda Galaxy Part I: global stellar density, morphology and metallicity properties

Authors: Rodrigo A. Ibata, Geraint F. Lewis, Alan W. McConnachie, Nicolas F. Martin, Michael J. Irwin, Annette M. N. Ferguson, Arif Babul, Edouard J. Bernard, Scott C. Chapman, Michelle Collins, Mark Fardal, A. D. Mackey, Julio Navarro, Jorge Penarrubia, R. Michael Rich, Nial Tanvir, Lawrence Widro

Abstract: We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into… ▽ More We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into stream-like populations and a smooth halo component. Fitting a three-dimensional halo model reveals that the most metal-poor populations ([Fe/H]<-1.7) are distributed approximately spherically (slightly prolate with ellipticity c/a=1.09+/-0.03), with only a relatively small fraction (42%) residing in discernible stream-like structures. The sphericity of the ancient smooth component strongly hints that the dark matter halo is also approximately spherical. More metal-rich populations contain higher fractions of stars in streams (86% for [Fe/H]>-0.6). The space density of the smooth metal-poor component has a global power-law slope of -3.08+/-0.07, and a non-parametric fit shows that the slope remains nearly constant from 30kpc to 300kpc. The total stellar mass in the halo at distances beyond 2 degrees is 1.1x10^10 Solar masses, while that of the smooth component is 3x10^9 Solar masses. Extrapolating into the inner galaxy, the total stellar mass of the smooth halo is plausibly 8x10^9 Solar masses. We detect a substantial metallicity gradient, which declines from [Fe/H]=-0.7 at R=30kpc to [Fe/H]=-1.5 at R=150kpc for the full sample, with the smooth halo being 0.2dex more metal poor than the full sample at each radius. While qualitatively in-line with expectations from cosmological simulations, these observations are of great importance as they provide a prototype template that such simulations must now be able to reproduce in quantitative detail. △ Less

Submitted 22 November, 2013; originally announced November 2013.

Comments: 21 pages, 17 figures, accepted for publication in the Astrophysical Journal

Major Substructure in the M31 Outer Halo: the South-West Cloud

Authors: N. F. Bate, A. R. Conn, B. McMonigal, G. F. Lewis, N. F. Martin, A. W. McConnachie, J. Veljanoski, A. D. Mackey, A. M. N. Ferguson, R. A. Ibata, M. J. Irwin, M. Fardal, A. P. Huxor, A. Babul

Abstract: We undertake the first detailed analysis of the stellar population and spatial properties of a diffuse substructure in the outer halo of M31. The South-West Cloud lies at a projected distance of ~100 kpc from the centre of M31, and extends for at least ~50 kpc in projection. We use Pan-Andromeda Archaeological Survey photometry of red giant branch stars to determine a distance to the South-West Cl… ▽ More We undertake the first detailed analysis of the stellar population and spatial properties of a diffuse substructure in the outer halo of M31. The South-West Cloud lies at a projected distance of ~100 kpc from the centre of M31, and extends for at least ~50 kpc in projection. We use Pan-Andromeda Archaeological Survey photometry of red giant branch stars to determine a distance to the South-West Cloud of 793 +/- 45 kpc. The metallicity of the cloud is found to be [Fe/H] = -1.3 +/- 0.1. This is consistent with the coincident globular clusters PAndAS-7 and PAndAS-8, which have metallicities determined using an independent technique of [Fe/H] = -1.35 +/- 0.15. We measure a brightness for the Cloud of M_V = -12.1 mag; this is ~75 per cent of the luminosity implied by the luminosity-metallicity relation. Under the assumption that the South-West Cloud is the visible remnant of an accreted dwarf satellite, this suggests that the progenitor object was amongst M31's brightest dwarf galaxies prior to disruption. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

The PAndAS view of the Andromeda satellite system - I. A Bayesian search for dwarf galaxies using spatial and color-magnitude information

Authors: Nicolas F. Martin, Rodrigo A. Ibata, Alan W. McConnachie, A. Dougal Mackey, Annette M. N. Ferguson, Michael J. Irwin, Geraint F. Lewis, Mark A. Fardal

Abstract: We present a generic algorithm to search for dwarf galaxies in photometric catalogs and apply it to the Pan-Andromeda Archaeological Survey (PAndAS). The algorithm is developed in a Bayesian framework and, contrary to most dwarf-galaxy-search codes, makes use of both the spatial and color-magnitude information of sources in a probabilistic approach. Accounting for the significant contamination fro… ▽ More We present a generic algorithm to search for dwarf galaxies in photometric catalogs and apply it to the Pan-Andromeda Archaeological Survey (PAndAS). The algorithm is developed in a Bayesian framework and, contrary to most dwarf-galaxy-search codes, makes use of both the spatial and color-magnitude information of sources in a probabilistic approach. Accounting for the significant contamination from the Milky Way foreground and from the structured stellar halo of the Andromeda galaxy, we recover all known dwarf galaxies in the PAndAS footprint with high significance, even for the least luminous ones. Some Andromeda globular clusters are also recovered and, in one case, discovered. We publish a list of the 143 most significant detections yielded by the algorithm. The combined properties of the 39 most significant isolated detections show hints that at least some of these trace genuine dwarf galaxies, too faint to be individually detected. Follow-up observations by the community are mandatory to establish which are real members of the Andromeda satellite system. The search technique presented here will be used in an upcoming contribution to determine the PAndAS completeness limits for dwarf galaxies. Although here tuned to the search of dwarf galaxies in the PAndAS data, the algorithm can easily be adapted to the search for any localised overdensity whose properties can be modeled reliably in the parameter space of any catalog. △ Less

Submitted 30 July, 2013; v1 submitted 29 July, 2013; originally announced July 2013.

Comments: 20 pages, 16 figures, 1 table; accepted for publication in ApJ. High res pdf available at https://www.dropbox.com/s/7zk7pme2wunwkjv/PAndAS_dwarf_galaxies.pdf

Inferring the Andromeda Galaxy's mass from its giant southern stream with Bayesian simulation sampling

Authors: Mark A. Fardal, Martin D. Weinberg, Arif Babul, Mike J. Irwin, Puragra Guhathakurta, Karoline M. Gilbert, Annette M. N. Ferguson, Rodrigo A. Ibata, Geraint F. Lewis, Nial R. Tanvir, Avon P. Huxor

Abstract: M31 has a giant stream of stars extending far to the south and a great deal of other tidal debris in its halo, much of which is thought to be directly associated with the southern stream. We model this structure by means of Bayesian sampling of parameter space, where each sample uses an N-body simulation of a satellite disrupting in M31's potential. We combine constraints on stellar surface densit… ▽ More M31 has a giant stream of stars extending far to the south and a great deal of other tidal debris in its halo, much of which is thought to be directly associated with the southern stream. We model this structure by means of Bayesian sampling of parameter space, where each sample uses an N-body simulation of a satellite disrupting in M31's potential. We combine constraints on stellar surface densities from the Isaac Newton Telescope survey of M31 with kinematic data and photometric distances. This combination of data tightly constrains the model, indicating a stellar mass at last pericentric passage of log(M_s / Msun) = 9.5+-0.1, comparable to the LMC. Any existing remnant of the satellite is expected to lie in the NE Shelf region beside M31's disk, at velocities more negative than M31's disk in this region. This rules out the prominent satellites M32 or NGC 205 as the progenitor, but an overdensity recently discovered in M31's NE disk sits at the edge of the progenitor locations found in the model. M31's virial mass is constrained in this model to be log(M200) = 12.3+-0.1, alleviating the previous tension between observational virial mass estimates and expectations from the general galactic population and the timing argument. The techniques used in this paper, which should be more generally applicable, are a powerful method of extracting physical inferences from observational data on tidal debris structures. △ Less

Submitted 11 July, 2013; originally announced July 2013.

Comments: 27 pages, 10 figures. Accepted by MNRAS

Kinematics of Outer Halo Globular Clusters in M31

Authors: J. Veljanoski, A. M. N. Ferguson, A. D. Mackey, A. P. Huxor, M. J. Irwin, P. Côté, N. R. Tanvir, E. J. Bernard, S. C. Chapman, R. A. Ibata, M. Fardal, G. F. Lewis, N. F. Martin, A. McConnachie, J. Peñarrubia

Abstract: We present the first kinematic analysis of the far outer halo globular cluster (GC) population in the Local Group galaxy M31. Our sample contains 53 objects with projected radii of ~20-130 kpc, of which 44 have no previous spectroscopic information. GCs with projected radii >30 kpc are found to exhibit net rotation around the minor axis of M31, in the same sense as the inner GCs, albeit with a sma… ▽ More We present the first kinematic analysis of the far outer halo globular cluster (GC) population in the Local Group galaxy M31. Our sample contains 53 objects with projected radii of ~20-130 kpc, of which 44 have no previous spectroscopic information. GCs with projected radii >30 kpc are found to exhibit net rotation around the minor axis of M31, in the same sense as the inner GCs, albeit with a smaller amplitude of 79 +/-19 km/s. The rotation-corrected velocity dispersion of the full halo GC sample is 106 +/-12 km/s, which we observe to decrease with increasing projected radius. We find compelling evidence for kinematic-coherence amongst GCs which project on top of halo substructure, including a clear signature of infall for GCs lying along the North-West stream. Using the tracer mass estimator, we estimate the dynamical mass of M31 within 200 kpc to be M_M31 = (1.2-1.5) +/- 0.2 x 10^12 M_sun. This value is highly dependent on the chosen model and assumptions within. △ Less

Submitted 29 March, 2013; originally announced March 2013.

Comments: 5 pages, 3 figures, accepted for publication in ApJL

A Vast Thin Plane of Co-rotating Dwarf Galaxies Orbiting the Andromeda Galaxy

Authors: Rodrigo A. Ibata, Geraint F. Lewis, Anthony R. Conn, Michael J. Irwin, Alan W. McConnachie, Scott C. Chapman, Michelle L. Collins, Mark Fardal, Annette M. N. Ferguson, Neil G. Ibata, A. Dougal Mackey, Nicolas F. Martin, Julio Navarro, R. Michael Rich, David Valls-Gabaud, Lawrence M. Widrow

Abstract: Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, a… ▽ More Dwarf satellite galaxies are thought to be the remnants of the population of primordial structures that coalesced to form giant galaxies like the Milky Way. An early analysis noted that dwarf galaxies may not be isotropically distributed around our Galaxy, as several are correlated with streams of HI emission, and possibly form co-planar groups. These suspicions are supported by recent analyses, and it has been claimed that the apparently planar distribution of satellites is not predicted within standard cosmology, and cannot simply represent a memory of past coherent accretion. However, other studies dispute this conclusion. Here we report the existence (99.998% significance) of a planar sub-group of satellites in the Andromeda galaxy, comprising approximately 50% of the population. The structure is vast: at least 400 kpc in diameter, but also extremely thin, with a perpendicular scatter <14.1 kpc (99% confidence). Radial velocity measurements reveal that the satellites in this structure have the same sense of rotation about their host. This finding shows conclusively that substantial numbers of dwarf satellite galaxies share the same dynamical orbital properties and direction of angular momentum, a new insight for our understanding of the origin of these most dark matter dominated of galaxies. Intriguingly, the plane we identify is approximately aligned with the pole of the Milky Way's disk and is co-planar with the Milky Way to Andromeda position vector. The existence of such extensive coherent kinematic structures within the halos of massive galaxies is a fact that must be explained within the framework of galaxy formation and cosmology. △ Less

Submitted 3 January, 2013; originally announced January 2013.

Comments: Published in the 3rd Jan 2013 issue of Nature. 19 pages, 4 figures, 1 three-dimensional interactive figure. To view and manipulate the 3-D figure, an Adobe Reader browser plug-in is required; alternatively save to disk and view with Adobe Reader

Journal ref: Nature 493, 62-65 (2013)

PAndAS in the mist: The stellar and gaseous mass within the halos of M31 and M33

Authors: Geraint F. Lewis, Robert Braun, Alan W. McConnachie, Michael J. Irwin, Rodrigo A. Ibata, Scott C. Chapman, Annette M. N. Ferguson, Nicolas F. Martin, Mark Fardal, John Dubinski, Larry Widrow, A. Dougal Mackey, Arif Babul, Nial Tanvir, Michael Rich

Abstract: Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the… ▽ More Large scale surveys of the prominent members of the Local Group have provided compelling evidence for the hierarchical formation of massive galaxies, revealing a wealth of substructure that is thought to be the debris from ancient and on-going accretion events. In this paper, we compare two extant surveys of the M31-M33 subgroup of galaxies; the Pan-Andromeda Archaeological Survey (PAndAS) of the stellar structure, and a combination of observations of the HI gaseous content, detected at 21cm. Our key finding is a marked lack of spatial correlation between these two components on all scales, with only a few potential overlaps between stars and gas.The paucity of spatial correlation significantly restricts the analysis of kinematic correlations, although there does appear to the HI kinematically associated with the Giant Stellar Stream where it passes the disk of M31. These results demonstrate that that different processes must significantly influence the dynamical evolution of the stellar and HI components of substructures, such as ram pressure driving gas away from a purely gravitational path. Detailed modelling of the offset between the stellar and gaseous substructure will provide a determination of the properties of the gaseous halo of M31 and M33. △ Less

Submitted 16 November, 2012; originally announced November 2012.

Comments: 11 pages, 6 figures. Accepted for publication in the Astrophysical Journal. Figure quality reduced. High quality version available at http://www.physics.usyd.edu.au/~gfl/Arxiv_Papers/PAndAS_Mist/

Report number: GFL-001

Unearthing Foundations of a Cosmic Cathedral: Searching the Stars for M33's Halo

Authors: Robert Cockcroft, Alan W. McConnachie, William E. Harris, Rodrigo Ibata, Mike J. Irwin, Annette M. N. Ferguson, Mark A. Fardal, Arif Babul, Scott C. Chapman, Geraint F. Lewis, Nicolas F. Martin, Thomas H. Puzia

Abstract: We use data from the Pan-Andromeda Archaeological Survey (PAndAS) to search for evidence of an extended halo component belonging to M33 (the Triangulum Galaxy). We identify a population of red giant branch (RGB) stars at large radii from M33's disk whose connection to the recently discovered extended "disk substructure" is ambiguous, and which may represent a "bona-fide" halo component. After firs… ▽ More We use data from the Pan-Andromeda Archaeological Survey (PAndAS) to search for evidence of an extended halo component belonging to M33 (the Triangulum Galaxy). We identify a population of red giant branch (RGB) stars at large radii from M33's disk whose connection to the recently discovered extended "disk substructure" is ambiguous, and which may represent a "bona-fide" halo component. After first correcting for contamination from the Milky Way foreground population and misidentified background galaxies, we average the radial density of RGB candidate stars over circular annuli centered on the galaxy and away from the disk substructure. We find evidence of a low-luminosity, centrally concentrated component that is everywhere in our data fainter than mu_V ~ 33 mag arcsec^(-2). The scale length of this feature is not well constrained by our data, but it appears to be of order r_exp ~ 20 kpc; there is weak evidence to suggest it is not azimuthally symmetric. Inspection of the overall CMD for this region that specifically clips out the disk substructure reveals that this residual RGB population is consistent with an old population with a photometric metallicity of around [Fe/H] ~ -2 dex, but some residual contamination from the disk substructure appears to remain. We discuss the likelihood that our findings represent a bona-fide halo in M33, rather than extended emission from the disk substructure. We interpret our findings in terms of an upper limit to M33's halo that is a few percent of its total luminosity, although its actual luminosity is likely much less. △ Less

Submitted 15 October, 2012; originally announced October 2012.

Comments: 57 pages, 13 figures, 2 tables. Accepted by the Monthly Notices of the Royal Astronomical Society

A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (Part II); Distances to the Satellites of M31

Authors: Anthony R. Conn, Rodrigo A. Ibata, Geraint F. Lewis, Quentin A. Parker, Daniel B. Zucker, Nicolas F. Martin, Alan W. McConnachie, Mike J. Irwin, Nial Tanvir, Mark A. Fardal, Annette M. N. Ferguson, Scott C. Chapman, David Valls-Gabaud

Abstract: In `A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (PART I),' a new technique was introduced for obtaining distances using the TRGB standard candle. Here we describe a useful complement to the technique with the potential to further reduce the uncertainty in our distance measurements by incorporating a matched-filter weighting scheme into the model likelihood calculations. In t… ▽ More In `A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (PART I),' a new technique was introduced for obtaining distances using the TRGB standard candle. Here we describe a useful complement to the technique with the potential to further reduce the uncertainty in our distance measurements by incorporating a matched-filter weighting scheme into the model likelihood calculations. In this scheme, stars are weighted according to their probability of being true object members. We then re-test our modified algorithm using random-realization artificial data to verify the validity of the generated posterior probability distributions (PPDs) and proceed to apply the algorithm to the satellite system of M31, culminating in a 3D view of the system. Further to the distributions thus obtained, we apply a satellite-specific prior on the satellite distances to weight the resulting distance posterior distributions, based on the halo density profile. Thus in a single publication, using a single method, a comprehensive coverage of the distances to the companion galaxies of M31 is presented, encompassing the dwarf spheroidals Andromedas I - III, V, IX-XXVII and XXX along with NGC147, NGC 185, M33 and M31 itself. Of these, the distances to Andromeda XXIV - XXVII and Andromeda XXX have never before been derived using the TRGB. Object distances are determined from high-resolution tip magnitude posterior distributions generated using the Markov Chain Monte Carlo (MCMC) technique and associated sampling of these distributions to take into account uncertainties in foreground extinction and the absolute magnitude of the TRGB as well as photometric errors. The distance PPDs obtained for each object both with, and without the aforementioned prior are made available to the reader in tabular form... △ Less

Submitted 21 September, 2012; originally announced September 2012.

Comments: 18 pages, 13 Figures

Journal ref: Published in the Astrophysical Journal 20 September 2012. ApJ, 758, 11

A spectroscopic survey of Andromeda's Western Shelf

Authors: Mark A. Fardal, Puragra Guhathakurta, Karoline M. Gilbert, Erik J. Tollerud, Jason S. Kalirai, Mikito Tanaka, Rachael Beaton, Masashi Chiba, Yutaka Komiyama, Masanori Iye

Abstract: The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius,… ▽ More The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius, the sample shows the wedge pattern expected from a radial shell, which is detected clearly here for the first time. This confirms predictions from an earlier model of formation of the GSS, which proposed that the W Shelf is a shell from the third orbital wrap of the same tidal debris stream that produces the GSS, with the main body of the progenitor lying in the second wrap. We calculate the distortions in the shelf wedge pattern expected from its outward expansion and angular momentum, and show that these effects are echoed in the data. In addition, a hot, relatively smooth spheroid population is clearly present. We construct a bulge-disk-halo N-body model that agrees with surface brightness and kinematic constraints, and combine it with a simulation of the GSS. From the contrasting kinematic signatures of the hot spheroid and shelf components, we decompose the observed stellar metallicity distribution into contributions from each component using a non-parametric mixture model. The shelf component's metallicity distribution matches previous observations of the GSS superbly, further strengthening the evidence they are connected and bolstering the case for a massive progenitor of this stream. △ Less

Submitted 12 June, 2012; originally announced June 2012.

Comments: Accepted for publication in MNRAS

The M31 Velocity Vector. II. Radial Orbit Towards the Milky Way and Implied Local Group Mass

Authors: Roeland P. van der Marel, Mark Fardal, Gurtina Besla, Rachael L. Beaton, Sangmo Tony Sohn, Jay Anderson, Tom Brown, Puragra Guhathakurta

Abstract: We determine the velocity vector of M31 with respect to the Milky Way and use this to constrain the mass of the Local Group, based on HST proper-motion measurements presented in Paper I. We construct N-body models for M31 to correct the measurements for the contributions from stellar motions internal to M31. We also estimate the center-of-mass motion independently, using the kinematics of satellit… ▽ More We determine the velocity vector of M31 with respect to the Milky Way and use this to constrain the mass of the Local Group, based on HST proper-motion measurements presented in Paper I. We construct N-body models for M31 to correct the measurements for the contributions from stellar motions internal to M31. We also estimate the center-of-mass motion independently, using the kinematics of satellite galaxies of M31 and the Local Group. All estimates are mutually consistent, and imply a weighted average M31 heliocentric transverse velocity of (v_W,v_N) = (-125.2+/-30.8, -73.8+/-28.4) km/s. We correct for the reflex motion of the Sun using the most recent insights into the solar motion within the Milky Way. This implies a radial velocity of M31 with respect to the Milky Way of V_rad = -109.3+/-4.4 km/s, and a tangential velocity V_tan = 17.0 km/s (<34.3 km/s at 1-sigma confidence). Hence, the velocity vector of M31 is statistically consistent with a radial (head-on collision) orbit towards the Milky Way. We revise prior estimates for the Local Group timing mass, including corrections for cosmic bias and scatter. Bayesian combination with other mass estimates yields M_LG = M_MW(vir) + M_M31(vir) = (3.17 +/- 0.57) x 10^12 solar masses. The velocity and mass results imply at 95% confidence that M33 is bound to M31, consistent with expectation from observed tidal deformations. (Abridged) △ Less

Submitted 30 May, 2012; originally announced May 2012.

Comments: 38 pages, 4 figures, to be published in ApJ. Version with high resolution figures and N-body movies available at http://www.stsci.edu/~marel/M31 . Press materials, graphics, and visualizations available at http://hubblesite.org/newscenter/archive/releases/2012/20

The SPLASH Survey: Kinematics of Andromeda's Inner Spheroid

Authors: Claire E. Dorman, Puragra Guhathakurta, Mark A. Fardal, Dustin Lang, Marla C. Geha, Kirsten M. Howley, Jason S. Kalirai, James S. Bullock, Jean-Charles Cuillandre, Julianne J. Dalcanton, Karoline M. Gilbert, Anil C. Seth, Erik J. Tollerud, Benjamin F. Williams, Basilio Yniguez

Abstract: The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Mar… ▽ More The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples. △ Less

Submitted 19 April, 2012; originally announced April 2012.

Comments: Accepted for publication in ApJ

The Recent Stellar Archeology of M31 - The Nearest Red Disk Galaxy

Authors: T. J. Davidge, A. W. McConnachie, M. A. Fardal, J. Fliri, D. Valls-Gabaud, S. C. Chapman, G. F. Lewis, R. M. Rich

Abstract: We examine the star-forming history (SFH) of the M31 disk during the past few hundred Myr. The luminosity functions (LFs) of main sequence stars at distances R_GC > 21 kpc (i.e. > 4 disk scale lengths) are matched by models that assume a constant star formation rate (SFR). However, at smaller R_GC the LFs suggest that during the past ~10 Myr the SFR was 2 - 3 times higher than during the preceding… ▽ More We examine the star-forming history (SFH) of the M31 disk during the past few hundred Myr. The luminosity functions (LFs) of main sequence stars at distances R_GC > 21 kpc (i.e. > 4 disk scale lengths) are matched by models that assume a constant star formation rate (SFR). However, at smaller R_GC the LFs suggest that during the past ~10 Myr the SFR was 2 - 3 times higher than during the preceding ~100 Myr. The rings of cool gas that harbor a significant fraction of the current star-forming activity are traced by stars with ages ~100 Myr, indicating that (1) these structures have ages of at least 100 Myr, and (2) stars in these structures do not follow the same relation between age and random velocity as their counterparts throughout the disks of other spiral galaxies, probably due to the inherently narrow orbital angular momentum distribution of the giant molecular clouds in these structures. The distribution of evolved red stars is not azimuthally symmetric, in the sense that their projected density along the north east segment of the major axis is roughly twice that on the opposite side of the galaxy. The north east arm of the major axis thus appears to be a fossil star-forming area that dates to intermediate epochs. Such a structure may be the consequence of interactions with a companion galaxy. △ Less

Submitted 27 March, 2012; originally announced March 2012.

Comments: To appear in The Astrophysical Journal

A Bayesian Approach to Locating the Red Giant Branch Tip Magnitude (Part I)

Authors: A. R. Conn, G. F. Lewis, R. A. Ibata, Q. A. Parker, D. B. Zucker, A. W. McConnachie, N. F. Martin, M. J. Irwin, N. Tanvir, M. A. Fardal, A. M. N. Ferguson

Abstract: We present a new approach for identifying the Tip of the Red Giant Branch (TRGB) which, as we show, works robustly even on sparsely populated targets. Moreover, the approach is highly adaptable to the available data for the stellar population under study, with prior information readily incorporable into the algorithm. The uncertainty in the derived distances is also made tangible and easily calcul… ▽ More We present a new approach for identifying the Tip of the Red Giant Branch (TRGB) which, as we show, works robustly even on sparsely populated targets. Moreover, the approach is highly adaptable to the available data for the stellar population under study, with prior information readily incorporable into the algorithm. The uncertainty in the derived distances is also made tangible and easily calculable from posterior probability distributions. We provide an outline of the development of the algorithm and present the results of tests designed to characterize its capabilities and limitations. We then apply the new algorithm to three M31 satellites: Andromeda I, Andromeda II and the fainter Andromeda XXIII, using data from the Pan-Andromeda Archaeological Survey (PAndAS), and derive their distances as $731^{(+ 5) + 18}_{(- 4) - 17}$ kpc, $634^{(+ 2) + 15}_{(- 2) - 14}$ kpc and $733^{(+ 13)+ 23}_{(- 11) - 22}$ kpc respectively, where the errors appearing in parentheses are the components intrinsic to the method, while the larger values give the errors after accounting for additional sources of error. These results agree well with the best distance determinations in the literature and provide the smallest uncertainties to date. This paper is an introduction to the workings and capabilities of our new approach in its basic form, while a follow-up paper shall make full use of the method's ability to incorporate priors and use the resulting algorithm to systematically obtain distances to all of M31's satellites identifiable in the PAndAS survey area. △ Less

Submitted 16 July, 2011; originally announced July 2011.

Comments: 11 pages, 18 figures

Journal ref: Approved for publication in the Astrophysical Journal ApJ83086 2011

Density Variations in the NW Star Stream of M31

Authors: R. G. Carlberg, Harvey B. Richer, Alan W. McConnachie, Mike Irwin, Rodrigo A. Ibata, Aaron L. Dotter, Scott Chapman, Mark Fardal, A. M. N. Ferguson, G. F. Lewis, Julio F. Navarro, Thomas H. Puzia, David Valls-Gabaud

Abstract: The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 system has found a star stream which extends about 120 kpc NW from the center of M31. The great length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means that it is unusually well suited for a measurement of stream gaps and clumps along its length as a test for the pr… ▽ More The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 system has found a star stream which extends about 120 kpc NW from the center of M31. The great length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means that it is unusually well suited for a measurement of stream gaps and clumps along its length as a test for the predicted thousands of dark matter sub-halos. The main result of this paper is that the density of the stream varies between zero and about three times the mean along its length on scales of 2 to 20 kpc. The probability that the variations are random fluctuations in the star density is less than 10^-5. As a control sample we search for density variations at precisely the same location in stars with metallicity higher than the stream, [Fe/H]=[0, -0.5] and find no variations above the expected shot noise. The lumpiness of the stream is not compatible with a low mass star stream in a smooth galactic potential, nor is it readily compatible with the disturbance caused by the visible M31 satellite galaxies. The stream's density variations appear to be consistent with the effects of a large population of steep mass function dark matter sub-halos, such as found in LCDM simulations, acting on an approximately 10Gyr old star stream. The effects of a single set of halo substructure realizations are shown for illustration, reserving a statistical comparison for another study. △ Less

Submitted 17 February, 2011; originally announced February 2011.

Comments: ApJ revised version submitted

Journal ref: 2011 ApJ, 731, 124

Testing Subhalo Abundance Matching in Cosmological Smoothed Particle Hydrodynamics Simulations

Authors: Vimal Simha, David Weinberg, Romeel Dave, Mark Fardal, Neal Katz, Benjamin D. Oppenheimer

Abstract: Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearl… ▽ More Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales. △ Less

Submitted 2 May, 2012; v1 submitted 22 November, 2010; originally announced November 2010.

Comments: 21 pages, 11 figures, Accepted by MNRAS

Galactic outflows and the kinematics of damped Lyman alpha absorbers

Authors: Sungryong Hong, Neal Katz, Romeel Davé, Mark Fardal, Dušan Kereš, Benjamin D. Oppenheimer

Abstract: The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previ… ▽ More The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies. △ Less

Submitted 27 August, 2010; v1 submitted 25 August, 2010; originally announced August 2010.

Comments: submitted to MNRAS

On the algorithms of radiative cooling in semi-analytic models

Authors: Yu Lu, Dušan Kereš, Neal Katz, H. J. Mo, Mark Fardal, Martin D. Weinberg

Abstract: We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothl… ▽ More We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations. △ Less

Submitted 28 September, 2011; v1 submitted 5 August, 2010; originally announced August 2010.

Comments: 1 figure added, MORGANA model added into the comparison, published in MNRAS

Journal ref: 2011MNRAS.416..660L

Intergalactic Dust Extinction in Hydrodynamic Cosmological Simulations

Authors: Ying Zu, David H. Weinberg, Romeel Davé, Mark Fardal, Neal Katz, Dusan Keres, Benjamin D. Oppenheimer

Abstract: Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volum… ▽ More Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust. △ Less

Submitted 5 January, 2011; v1 submitted 24 May, 2010; originally announced May 2010.

Comments: Matches version accepted by MNRAS

Feedback and Recycled Wind Accretion: Assembling the z=0 Galaxy Mass Function

Authors: Benjamin D. Oppenheimer, Romeel Davé, Dušan Kereš, Mark Fardal, Neal Katz, Juna A. Kollmeier, David H. Weinberg

Abstract: We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically d… ▽ More We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF. △ Less

Submitted 13 May, 2010; v1 submitted 2 December, 2009; originally announced December 2009.

Comments: 16 pages, 7 figures, accepted by MNRAS

Asymmetric Warfare: M31 and its Satellites

Authors: Mark A. Fardal

Abstract: Photometric surveys of M31's halo vividly illustrate the wreckage caused by hierarchical galaxy formation. Several of M31's satellites are being disrupted by M31's tidal field, among them M33 and And I, while other tidal structures are the corpses of satellites already destroyed. The extent to which M31's satellites have left battle scars upon it is unknown; to answer this we need accurate orbit… ▽ More Photometric surveys of M31's halo vividly illustrate the wreckage caused by hierarchical galaxy formation. Several of M31's satellites are being disrupted by M31's tidal field, among them M33 and And I, while other tidal structures are the corpses of satellites already destroyed. The extent to which M31's satellites have left battle scars upon it is unknown; to answer this we need accurate orbits and masses of the perturbers. I focus here on M31's 150-kpc-long Giant Southern Stream (GSS) as an example of how these can be determined even in the absence of a visible progenitor. Comparing N-body models to photometric and spectroscopic data, I find this stream resulted from the disruption of a large satellite galaxy by a close passage about 750 Myr ago. The GSS is connected to several other debris structures in M31's halo. Bayesian sampling of the simulations estimates the progenitor's initial mass as log(Mstar/Msun) = 9.5 +- 0.2, showing it was one of the most massive Local Group galaxies until quite recently. The stream model constrains M31's halo mass to be (1.8+-0.5)x10^{12} Msun. While these small uncertainties neglect several important degrees of freedom, they are likely to remain good even with a more complete model. Future work on M31's satellites and streams will provide independent constraints on M31's mass, and reveal the shared history of M31 and its halo components. △ Less

Submitted 19 October, 2009; v1 submitted 15 October, 2009; originally announced October 2009.

Comments: 6 pages, 3 figures. To appear in the Proceedings of "Galaxy Wars: Stellar Populations and Star Formation in Interacting Galaxies" Conference, B. Smith, N. Bastian, J. Higdon and S. Higdon eds, in press. Minor figure and caption edits