Source author record

Rodrigo A. Ibata

Rodrigo A. Ibata appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

57 published item(s)

preprint2022arXiv

A stellar stream remnant of a globular cluster below the metallicity floor

Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least ~0.2 percent of the solar metallicity ([Fe/H] > -2.7). This metallicity floor appears universal and it has been proposed that proto-galaxies that merge into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day. Here, we report the discovery of a stellar stream, C-19, whose metallicity is less than 0.05 per cent the solar metallicity ([Fe/H]=-3.38 +/- 0.06 (stat.) +/- 0.20 (syst.)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo.

preprint2022arXiv

New constraints on the dark matter density profiles of dwarf galaxies from proper motions of globular cluster streams

The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams ($σ_{v_\mathrm{Tan}}$) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that were accreted from: cuspy subhalos produce streams with larger $σ_{v_\mathrm{Tan}}$ than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their $σ_{v_\mathrm{Tan}}$ values with simulations. The measured $σ_{v_\mathrm{Tan}}$ values are consistent with both an ``in situ'' origin and with accretion inside cored subhalos of $M\sim 10^{8-9}M_{\odot}$ (or very low-mass cuspy subhalos of mass $\sim 10^8M_{\odot}$). Despite the large current uncertainties in $σ_{v_\mathrm{Tan}}$, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of $M_{\rm subhalo}\buildrel > \over \sim$ $10^9 M_{\odot}$. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles.

preprint2022arXiv

The Complexity of the Cetus Stream Unveiled from the Fusion of STREAMFINDER and StarGO

We combine the power of two stream-searching tools, STREAMFINDER and StarGO applied to the Gaia EDR3 data, to detect stellar debris belonging to the Cetus stream system that forms a complex, nearly polar structure around the Milky Way. In this work, we find the southern extensions of the northern Cetus stream as the Palca stream and a new southern stream, which overlap on the sky but have different distances. These two stream wraps extend over more than $\sim100°$ on the sky ($-60°<δ<+40°$). The current N-body model of the system reproduces both wraps in the trailing arm. We also show that the Cetus system is confidently associated with the Triangulum/Pisces, Willka Yaku, and the recently discovered C-20 streams. The association with the ATLAS-Aliqa Uma stream is much weaker. All of these stellar debris are very metal-poor, comparable to the average metallicity of the southern Cetus stream with [Fe/H] = $-2.17\pm0.2$. The estimated stellar mass of the Cetus progenitor is at least $10^{5.6}\rm\,M_\odot$, compatible with Ursa Minor or Draco dwarf galaxies. The associated globular cluster with similar stellar mass, NGC 5824 very possibly was accreted in the same group infall. The multi-wrap Cetus stream is a perfect example of a dwarf galaxy that has undergone several periods of stripping, leaving behind debris at multiple locations in the halo. The full characterization of such systems is crucial to unravel the assembly history of the Milky Way and, as importantly, to provide nearby fossils to study ancient low-mass dwarf galaxies.

preprint2022arXiv

The Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I'itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- "C-19" ([Fe/H]$=-3.4$ dex), "Sylgr" ([Fe/H]$=-2.9$ dex) and "Phoenix" ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

preprint2022arXiv

The Hough Stream Spotter: A New Method for Detecting Linear Structure in Resolved Stars and Application to the Stellar Halo of M31

Stellar streams from globular clusters (GCs) offer constraints on the nature of dark matter and have been used to explore the dark matter halo structure and substructure of our Galaxy. Detection of GC streams in other galaxies would broaden this endeavor to a cosmological context, yet no such streams have been detected to date. To enable such exploration, we develop the Hough Stream Spotter (HSS), and apply it to the Pan-Andromeda Archaeological Survey (PAndAS) photometric data of resolved stars in M31's stellar halo. We first demonstrate that our code can re-discover known dwarf streams in M31. We then use the HSS to blindly identify 27 linear GC stream-like structures in the PAndAS data. For each HSS GC stream candidate, we investigate the morphologies of the streams and the colors and magnitudes of all stars in the candidate streams. We find that the five most significant detections show a stronger signal along the red giant branch in color-magnitude diagrams (CMDs) than spurious non-stream detections. Lastly, we demonstrate that the HSS will easily detect globular cluster streams in future Nancy Grace Roman Space Telescope data of nearby galaxies. This has the potential to open up a new discovery space for GC stream studies, GC stream gap searches, and for GC stream-based constraints on the nature of dark matter.

preprint2022arXiv

The PAndAS View of the Andromeda Satellite System. III. Dwarf galaxy detection limits

We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the surveyed area because of varying levels of contamination from both the Milky Way foreground stars and the stellar halo of Andromeda. We therefore provide recovery fractions that are a function of size, luminosity, and location within the survey on a scale of 1 square degree. Overall, the effective surface brightness for a 50-percent detection rate range between 28 and 30 mag per square arcsecond. This is in line with expectations for a search that relies on photometric data that are as deep as the PAndAS survey. The derived detection limits are an essential ingredient on the path to constraining the global properties of Andromeda's system of satellite dwarf galaxies and, more broadly, to provide constraints on dwarf galaxy formation and evolution in a cosmological context.

preprint2022arXiv

The Pristine survey -- XVII. The C-19 stream is dynamically hot and more extended than previously thought

The C-19 stream is the most metal poor stellar system ever discovered, with a mean metallicity $[Fe/H] = -3.38\pm0.06$. Its low metallicity dispersion ($σ_{\rm [Fe/H]}$ $<$ 0.18 at the 95\% confidence level) as well as variations in sodium abundances strongly suggest a globular cluster origin. In this work, we use VLT/UVES spectra of seven C-19 stars to derive more precise velocity measurements for member stars, and to identify two new members with radial velocities and metallicities consistent with the stream's properties. One of these new member stars is located 30 $°$ away from the previously identified body of C-19, implying that the stream is significantly more extended than previously known and that more members likely await discovery. In the main part of C-19, we measure a radial velocity dispersion $σ_v$ = 6.2$^{+2.0}_{-1.4}$ km s$^{-1}$ from nine members, and a stream width of 0.56$°\pm0.08°$, equivalent to $\sim$158 pc at a heliocentric distance of 18 kpc. These confirm that C-19 is comparatively hotter, dynamically, than other known globular cluster streams and shares the properties of faint dwarf galaxy streams. On the other hand, the variations in the Na abundances of the three newly observed bright member stars, the variations in Mg and Al for two of them, and the normal Ba abundance of the one star where it can be measured provide further evidence for a globular cluster origin. The tension between the dynamical and chemical properties of C-19 suggests that its progenitor experienced a complex birth environment or disruption history.

preprint2019arXiv

A detailed study of Andromeda XIX, an extreme local analogue of ultra diffuse galaxies

With a central surface brightness of $μ_0=29.3$ mag. per sq. arcsec, and half-light radius of $r_{\rm half}=3.1^{+0.9}_{-1.1}$~kpc, Andromeda XIX (And XIX) is an extremely diffuse satellite of Andromeda. We present spectra for $\sim100$ red giant branch stars in this galaxy, plus 16 stars in a nearby stellar stream. With this exquisite dataset, we re-derive the properties of And XIX, measuring a systemic velocity of $v_r=-109.0\pm1.6$ km/s and a velocity dispersion of $σ_v = 7.8^{+1.7}_{-1.5}$ km/s (higher than derived in our previous work). We marginally detect a velocity gradient along the major axis of ${\rm d}v/{\rm d}r = -2.1\pm1.8$ km/s kpc$^{-1}$. We find its mass-to-light ratio is higher than galaxies of comparable stellar mass ($[M/L]_{\rm half} = 278^{+146}_{-198}M_\odot/L_\odot$), but its dynamics place it in a halo with a similar total mass to these galaxies. This could suggest that And XIX is a "puffed up" dwarf galaxy, whose properties have been altered by tidal processes, similar to its Milky Way counterpart, Antlia II. For the nearby stream, we measure $v_r=-279.2\pm3.7$ km/s, and $σ_v=13.8^{+3.5}_{-2.6}$ km/s. We measure its metallicity, and find it to be more metal rich than And XIX, implying that the two features are unrelated. Finally, And XIX's dynamical and structural properties imply it is a local analogue to ultra diffuse galaxies (UDGs). Its complex dynamics suggest that the masses of distant UDGs measured from velocity dispersions alone should be carefully interpreted

preprint2019arXiv

A SkyMapper view of the Large Magellanic Cloud: The dynamics of stellar populations

We present the first SkyMapper stellar population analysis of the Large Magellanic Cloud (hereafter LMC),including the identification of 3578 candidate Carbon Stars through their extremely red $g-r$ colours. Coupled with Gaia astrometry, we analyse the distribution and kinematics of this Carbon Star population, finding the LMC to be centred at $(R.A., Dec.) = (80.90^{\circ}\pm{0.29}, -68.74^{\circ}\pm{0.12})$, with a bulk proper motion of $(μ_α,μ_δ) = (1.878\pm0.007,0.293\pm0.018) \mathrm{mas\ yr^{-1}}$ and a disk inclination of $i = 25.6^{\circ}\pm1.1$ at position angle $θ= 135.6^{\circ}\pm 3.3^{\circ}$. We complement this study with the identification and analysis of additional stellar populations, finding that the dynamical centre for Red Giant Branch (RGB) stars is similar to that seen for the Carbon Stars, whereas for young stars the dynamical centre is significantly offset from the older populations. This potentially indicates that the young stars were formed as a consequence of a strong tidal interaction, probably with the Small Magellanic Cloud (SMC). In terms of internal dynamics, the tangential velocity profile increases linearly within $\sim3\ \mathrm{kpc}$, after which it maintains an approximately constant value of $V_{rot} = 83.6\pm 1.7 \mathrm{km\ s^{-1}}$ until $\sim7 \mathrm{kpc}$. With an asymmetric drift correction, we estimate the mass within $7 \mathrm{kpc}$ to be $M_{\rm LMC}(<7\mathrm{kpc}) = (2.5\pm0.1)\times10^{10}\ \mathrm{M}_{\odot}$ and within the tidal radius ($\sim 30\ \mathrm{kpc}$) to be $M_{\rm LMC}(<30\mathrm{kpc}) = (1.06 \pm 0.32)\times10^{11}\ \mathrm{M}_{\odot}$, consistent with other recent measurements.

preprint2019arXiv

Phase-space correlation in stellar streams of the Milky Way halo: The clash of Kshir and GD-1

We report the discovery of a 70-degree long stellar stream in the Milky Way halo, which criss-crosses the well known "GD-1" stream. We show that this new stellar structure ("Kshir") and GD-1 lie at similar distance, and are remarkably correlated in kinematics. We propose several explanations for the nature of this new structure and its possible association with GD-1. However, a scenario in which these two streams were accreted onto the Milky Way within the same dark matter sub-halo seems to provide a natural explanation for their phase-space entanglement, and other complexities of this coupled-system.

preprint2019arXiv

The Canada-France Imaging Survey: Reconstructing the Milky Way Star Formation History from its White Dwarf Population

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

preprint2019arXiv

The Pristine survey X: a large population of low-metallicity stars permeates the Galactic disk

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0σ$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

preprint2016arXiv

Is Draco II one of the faintest dwarf galaxies? First study from Keck/DEIMOS spectroscopy

We present the first spectroscopic analysis of the faint and compact stellar system Draco II (Dra II, M_V=-2.9+/-0.8, r_h=19^{+8}_{-6} pc), recently discovered in the Pan-STARRS1 3π survey. The observations, conducted with DEIMOS on the Keck II telescope, establish some of its basic characteristics: the velocity data reveal a narrow peak with 9 member stars at a systemic heliocentric velocity <v_r>=-347.6^{+1.7}_{-1.8} km/s, thereby confirming Dra II is a satellite of the Milky Way; we infer a velocity dispersion with σ_{vr}=2.9+/-2.1 km/s (<8.4 km/s at the 95% confidence level), which implies log_{10}(M_{1/2})=5.5^{+0.4}_{-0.6} and log_{10}((M/L)_{1/2})=2.7^{+0.5}_{-0.8}, in Solar units; furthermore, very weak Calcium triplet lines in the spectra of the high signal-to-noise member stars imply [Fe/H]<-2.1, whilst variations in the line strengths of two stars with similar colours and magnitudes suggest a metallicity spread in Dra II. These new data cannot clearly discriminate whether Draco II is a star cluster or amongst the faintest, most compact, and closest dwarf galaxies. However, the sum of the three --- individually inconclusive --- pieces of evidence presented here, seems to favour the dwarf galaxy interpretation.

preprint2016arXiv

Major Substructure in the M31 Outer Halo: Distances and Metallicities along the Giant Stellar Stream

We present a renewed look at M31's Giant Stellar Stream along with the nearby structures Stream C and Stream D, exploiting a new algorithm capable of fitting to the red giant branch (RGB) of a structure in both colour and magnitude space. Using this algorithm, we are able to generate probability distributions in distance, metallicity and RGB width for a series of subfields spanning these structures. Specifically, we confirm a distance gradient of approximately 20 kpc per degree along a 6 degree extension of the Giant Stellar Stream, with the farthest subfields from M31 lying ~ 120 kpc more distant than the inner-most subfields. Further, we find a metallicity that steadily increases from -0.7^{+0.1}_{-0.1} dex to -0.2^{+0.2}_{-0.1} dex along the inner half of the stream before steadily dropping to a value of -1.0^{+0.2}_{-0.2} dex at the farthest reaches of our coverage. The RGB width is found to increase rapidly from 0.4^{+0.1}_{-0.1} dex to 1.1^{+0.2}_{-0.1} dex in the inner portion of the stream before plateauing and decreasing marginally in the outer subfields of the stream. In addition, we estimate Stream C to lie at a distance between 794 and 862 kpc and Stream D between 758 kpc and 868 kpc. We estimate the median metallicity of Stream C to lie in the range -0.7 to -1.6 dex and a metallicity of -1.1^{+0.3}_{-0.2} dex for Stream D. RGB widths for the two structures are estimated to lie in the range 0.4 to 1.2 dex and 0.3 to 0.7 dex respectively. In total, measurements are obtained for 19 subfields along the Giant Stellar Stream, 4 along Stream C, 5 along Stream D and 3 general M31 spheroid fields for comparison. We thus provide a higher resolution coverage of the structures in these parameters than has previously been available in the literature.

preprint2016arXiv

On the Stability of Satellite Planes I: Effects of Mass, Velocity, Halo Shape and Alignment

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

preprint2016arXiv

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

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.

preprint2016arXiv

Triangulum II: a very metal-poor and dynamically hot stellar system

We present a study of the recently discovered compact stellar system Triangulum II. From observations conducted with the DEIMOS spectrograph on Keck II, we obtained spectra for 13 member stars that follow the CMD features of this very faint stellar system and include two bright red giant branch stars. Tri II has a very negative radial velocity (<v_r>=-383.7^{+3.0}_{-3.3} km/s) that translates to <v_{r,gsr}> ~ -264 km/s and confirms it is a Milky Way satellite. We show that, despite the small data set, there is evidence that Tri II has complex internal kinematics. Its radial velocity dispersion increases from 4.4^{+2.8}_{-2.0} km/s in the central 2' to 14.1^{+5.8}_{-4.2} km/s outwards. The velocity dispersion of the full sample is inferred to be σ_{vr}=9.9^{+3.2}_{-2.2} km/s. From the two bright RGB member stars we measure an average metallicity <[Fe/H]>=-2.6+/-0.2, placing Tri II among the most metal-poor Milky Way dwarf galaxies. In addition, the spectra of the fainter member stars exhibit differences in their line widths that could be the indication of a metallicity dispersion in the system. All these properties paint a complex picture for Tri II, whose nature and current state are largely speculative. The inferred metallicity properties of the system however lead us to favor a scenario in which Tri II is a dwarf galaxy that is either disrupting or embedded in a stellar stream.

preprint2015arXiv

A New Faint Milky Way Satellite Discovered in the Pan-STARRS1 3 pi Survey

We present the discovery of a faint Milky Way satellite, Laevens 2/Triangulum II, found in the Panoramic Survey Telescope And Rapid Response System (Pan-STARRS 1) 3 pi imaging data and confirmed with follow-up wide-field photometry from the Large Binocular Cameras. The stellar system, with an absolute magnitude of M_V=-1.8 +/-0.5, a heliocentric distance of 30 +2/-2 kpc, and a half-mass radius of 34 +9/-8 pc, shows remarkable similarity to faint, nearby, small satellites such as Willman 1, Segue 1, Segue 2, and Boötes II. The discovery of Laevens 2/Triangulum II further populates the region of parameter space for which the boundary between dwarf galaxies and globular clusters becomes tenuous. Follow-up spectroscopy will ultimately determine the nature of this new satellite, whose spatial location hints at a possible connection with the complex Triangulum-Andromeda stellar structures.

preprint2015arXiv

Eppur si muove: Positional and kinematic correlations of satellite pairs in the low Z universe

We have recently shown (Ibata et al. 2014) that pairs of satellite galaxies located diametrically opposite each other around their host possess predominantly anti-correlated velocities. This is consistent with a scenario in which $\sim 50$% of satellite galaxies belong to kinematically-coherent rotating planar structures, similar to those detected around the giant galaxies of the Local Group. Here we extend this analysis, examining the incidence of satellites of giant galaxies drawn from an SDSS photometric redshift catalog. We find that there is a $\sim 17$% overabundance ($> 3 σ$ significance) of candidate satellites at positions diametrically opposite a spectroscopically confirmed satellite. We show that cosmological simulations do not possess this property when the contamination is included, and that there are in fact, after subtracting contamination, 2 to 3 times more satellites diametrically opposed to a spectroscopically confirmed satellite than at $90°$ from it. We also examine the correlation between the satellite pair positions and the orientation of the host galaxy major axis. We find that those satellite pairs with anti-correlated velocities have a strong preference ($\sim 3:1$) to align with the major axis of the host whereas those with correlated velocities display the opposite behavior. This correlation of the satellite alignments appears to be stronger than the well-documented preference of satellites to be located close to the major axis of their host. We finally show that repeating a similar analysis to Ibata et al. (2014) with same-side satellites is generally hard to interpret, but is not inconsistent with our previous results when strong quality-cuts are applied on the sample. All these unexpected correlations strongly suggest that a substantial fraction of satellite galaxies are causally-linked in their formation and evolution.

preprint2015arXiv

Feeling the pull, a study of natural Galactic accelerometers. I: photometry of the delicate stellar stream of the Palomar 5 globular cluster

We present an analysis of wide-field photometric surveys of the Palomar 5 globular cluster and its stellar stream, based on g- and r-band measures together with narrow-band DDO51 photometry. In this first study, we use the deep (g,r) data to measure the incidence of gaps and peaks along the stream. Examining the star-counts profile of the stream plus contaminating populations, we find no evidence for significant under-densities, and find only a single significant over-density. This is at odds with earlier studies based on matched-filter maps derived from shallower SDSS data if the contaminating population possesses plausible spatial properties. The lack of substantial sub-structure along the stream may be used in future dynamical simulations to examine the incidence of dark matter sub-halos in the Galactic halo. We also present a measurement of the relative distances along the stream which we use to create the deepest wide-field map of this system to date.

preprint2015arXiv

NGC 147, NGC 185 and CassII: a genetic approach to orbital properties, star formation and tidal debris

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

preprint2015arXiv

The Spatially-Resolved Star Formation History of the M31 Outer Disc

We present deep Hubble Space Telescope Advanced Camera for Surveys observations of the stellar populations in two fields lying at 20 and 23 kpc from the centre of M31 along the south-west semi-major axis. These data enable the construction of colour-magnitude diagrams reaching the oldest main-sequence turn-offs (~13 Gyr) which, when combined with another field at 25 kpc from our previous work, we use to derive the first precision constraints on the spatially-resolved star formation history of the M31 disc. The star formation rates exhibit temporal as well as field-to-field variations, but are generally always within a factor of two of their time average. There is no evidence of inside-out growth over the radial range probed. We find a median age of ~7.5 Gyr, indicating that roughly half of the stellar mass in the M31 outer disc was formed before z ~ 1. We also find that the age-metallicity relations (AMRs) are smoothly increasing from [Fe/H]~-0.4 to solar metallicity between 10 and 3 Gyr ago, contrary to the flat AMR of the Milky Way disc at a similar number of scale lengths. Our findings provide insight on the roles of stellar feedback and radial migration in the formation and evolution of large disc galaxies.

preprint2014arXiv

A thousand shadows of Andromeda: rotating planes of satellites in the Millennium-II cosmological simulation

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

preprint2014arXiv

Comparing the observable properties of dwarf galaxies on and off the Andromeda plane

The thin, extended planes of satellite galaxies detected around both the Milky Way and Andromeda are not a natural prediction of the LCDM paradigm. Galaxies in these distinct planes may have formed and evolved in a different way (e.g., tidally) to their off-plane neighbours. If this were the case, one would expect the on- and off-plane dwarf galaxies in Andromeda to have experienced different evolutionary histories, which should be reflected by the chemistries, dynamics, and star formation histories of the two populations. In this work, we present new, robust kinematic observations for 2 on-plane M31 dSphs (And XVI and XVII) and compile and compare all available observational metrics for the on- and off-plane dwarfs to search for a signal that would corroborate such a hypothesis. We find that, barring their spatial alignment, the on- and off-plane Andromeda dwarf galaxies are indistinguishable from one another, arguing against vastly different formative and evolutionary histories for these two populations.

preprint2014arXiv

Resolving the mass--anisotropy degeneracy of the spherically symmetric Jeans equation I: theoretical foundation

A widely employed method for estimating the mass of stellar systems with apparent spherical symmetry is dynamical modelling using the spherically symmetric Jeans equation. Unfortunately this approach suffers from a degeneracy between the assumed mass density and the second order velocity moments. This degeneracy can lead to significantly different predictions for the mass content of the system under investigation, and thus poses a barrier for accurate estimates of the dark matter content of astrophysical systems. In a series of papers we describe an algorithm that removes this degeneracy and allows for unbiased mass estimates of systems of constant or variable mass-to-light ratio. The present contribution sets the theoretical foundation of the method that reconstructs a unique kinematic profile for some assumed free functional form of the mass density. The essence of our method lies in using flexible B-spline functions for the representation of the radial velocity dispersion in the spherically symmetric Jeans equation. We demonstrate our algorithm through an application to synthetic data for the case of an isotropic King model with fixed mass-to-light ratio, recovering excellent fits of theoretical functions to observables and a unique solution. The mass-anisotropy degeneracy is removed to the extent that, for an assumed functional form of the potential and mass density pair $(Φ,ρ)$, and a given set of line-of-sight velocity dispersion $σ_{los}^2$ observables, we recover a unique profile for $σ_{rr}^2$ and $σ_{tt}^2$. Our algorithm is simple, easy to apply and provides an efficient means to reconstruct the kinematic profile.

preprint2014arXiv

Resolving the mass--anisotropy degeneracy of the spherically symmetric Jeans equation II: optimum smoothing and model validation

The spherical Jeans equation is widely used to estimate the mass content of a stellar systems with apparent spherical symmetry. However, this method suffers from a degeneracy between the assumed mass density and the kinematic anisotropy profile, $β(r)$. In a previous work, we laid the theoretical foundations for an algorithm that combines smoothing B-splines with equations from dynamics to remove this degeneracy. Specifically, our method reconstructs a unique kinematic profile of $σ_{rr}^2$ and $σ_{tt}^2$ for an assumed free functional form of the potential and mass density $(Φ,ρ)$ and given a set of observed line-of-sight velocity dispersion measurements, $σ_{los}^2$. In Paper I (submitted to MNRAS: MN-14-0101-MJ) we demonstrated the efficiency of our algorithm with a very simple example and we commented on the need for optimum smoothing of the B-spline representation; this is in order to avoid unphysical variational behaviour when we have large uncertainty in our data. In the current contribution we present a process of finding the optimum smoothing for a given data set by using information of the behaviour from known ideal theoretical models. Markov Chain Monte Carlo methods are used to explore the degeneracy in the dynamical modelling process. We validate our model through applications to synthetic data for systems with constant or variable mass-to-light ratio $Υ$. In all cases we recover excellent fits of theoretical functions to observables and unique solutions. Our algorithm is a robust method for the removal of the mass-anisotropy degeneracy of the spherically symmetric Jeans equation for an assumed functional form of the mass density.

preprint2014arXiv

Spectroscopy of the three distant Andromedan satellites Cassiopeia III, Lacerta I, and Perseus I

We present Keck II/DEIMOS spectroscopy of the three distant dwarf galaxies of M31 Lacerta I, Cassiopeia III, and Perseus I, recently discovered within the Pan-STARRS1 3πimaging survey. The systemic velocities of the three systems (v_{r,helio} = -198.4 +/- 1.1 km/s, -371.6 +/- 0.7 km/s, and -326 +/- 3 km/s, respectively) confirm that they are satellites of M31. In the case of Lacerta I and Cassiopeia III, the high quality of the data obtained for 126 and 212 member stars, respectively, yields reliable constraints on their global velocity dispersions (σ_{vr} = 10.3 +/- 0.9 km/s and 8.4 +/- 0.6 km/s, respectively), leading to dynamical-mass estimates for both of ~4x10^7 Msun within their half-light radius. These translate to V-band mass-to-light ratios of 15^{+12}_{-9} and 8^{+9}_{-5} in solar units. We also use our spectroscopic data to determine the average metallicity of the 3 dwarf galaxies ([Fe/H] = -2.0 +/- 0.1, -1.7 +/- 0.1, and -2.0 +/- 0.2, respectively). All these properties are typical of dwarf galaxy satellites of Andromeda with their luminosity and size.

preprint2014arXiv

The Nature and Origin of Substructure in the Outskirts of M31 -- II. Detailed Star Formation Histories

While wide-field surveys of M31 have revealed much substructure at large radii, understanding the nature and origin of this material is not straightforward from morphology alone. Using deep HST/ACS data, we have derived further constraints in the form of quantitative star formation histories (SFHs) for 14 inner halo fields which sample diverse substructures. In agreement with our previous analysis of colour-magnitude diagram morphologies, we find the resultant behaviours can be broadly separated into two categories. The SFHs of 'disc-like' fields indicate that most of their mass has formed since z~1, with one quarter of the mass formed in the last 5 Gyr. We find 'stream-like' fields to be on average 1.5 Gyr older, with <10 percent of their stellar mass formed within the last 5 Gyr. These fields are also characterised by an age--metallicity relation showing rapid chemical enrichment to solar metallicity by z=1, suggestive of an early-type progenitor. We confirm a significant burst of star formation 2 Gyr ago, discovered in our previous work, in all the fields studied here. The presence of these young stars in our most remote fields suggests that they have not formed in situ but have been kicked-out from through disc heating in the recent past.

preprint2014arXiv

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

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.

preprint2014arXiv

Velocity anti-correlation of diametrically opposed galaxy satellites in the low redshift universe

Recent work has shown that both the Milky Way and the Andromeda galaxies possess the unexpected property that their dwarf satellite galaxies are aligned in thin and kinematically coherent planar structures. It is now important to evaluate the incidence of such planar structures in the larger galactic population, since the Local Group may not be a sufficiently representative environment. Here we report that the measurement of the velocity of pairs of diametrically opposed galaxy satellites provides a means to determine statistically the prevalence of kinematically coherent planar alignments. In the local universe (redshift $z<0.05$), we find that such satellite pairs out to a galactocentric distance of $150$ kpc are preferentially anti-correlated in their velocities (99.994% confidence level), and that the distribution of galaxies in the larger scale environment (beyond 150 kpc and up to $\approx 2$ Mpc) is strongly elongated along the axis joining the inner satellite pair ($>7σ$ confidence). Our finding may indicate that co-rotating planes of satellites, similar to that seen around the Andromeda galaxy, are ubiquitous in nature, while their coherent motion also suggests that they are a significant repository of angular momentum on $\sim 100$ kpc scales.

preprint2013arXiv

A kinematic study of the Andromeda dwarf spheroidal system

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 Andromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS spectrographs. Based on their g-i colors (taken with the CFHT MegaCam imager), physical positions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed star. Using this information, the velocity dispersions, central masses and central densities of the dark matter halos are calculated for these objects, and compared with the properties of the Milky Way dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, all of which have large half-light radii (>700 pc) and low velocity dispersions (sigma_v<5 km/s). In addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]_{half}=10.3^{+7.0}_{-6.7}, making it consistent with a simple stellar system with no appreciable dark matter component within its 1 sigma uncertainties. We suggest that the structure of the dark matter halos of these outliers have been significantly altered by tides.

preprint2013arXiv

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

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.

preprint2013arXiv

Dynamical Modeling of NGC 6809: Selecting the best model using Bayesian Inference

The precise cosmological origin of globular clusters remains uncertain, a situation hampered by the struggle of observational approaches in conclusively identifying the presence, or not, of dark matter in these systems. In this paper, we address this question through an analysis of the particular case of NGC 6809. While previous studies have performed dynamical modeling of this globular cluster using a small number of available kinematic data, they did not perform appropriate statistical inference tests for the choice of best model description; such statistical inference for model selection is important since, in general, different models can result in significantly different inferred quantities. With the latest kinematic data, we use Bayesian inference tests for model selection and thus obtain the best fitting models, as well as mass and dynamic mass-to-light ratio estimates. For this, we introduce a new likelihood function that provides more constrained distributions for the defining parameters of dynamical models. Initially we consider models with a known distribution function, and then model the cluster using solutions of the spherically symmetric Jeans equation; this latter approach depends upon the mass density profile and anisotropy $β$ parameter. In order to find the best description for the cluster we compare these models by calculating their Bayesian evidence. We find smaller mass and dynamic mass-to-light ratio values than previous studies, with the best fitting Michie model for a constant mass-to-light ratio of $Υ= 0.90^{+0.14}_{-0.14}$ and $M_{\text{dyn}}=6.10^{+0.51}_{-0.88} \times 10^4 M_{\odot}$. We exclude the significant presence of dark matter throughout the cluster, showing that no physically motivated distribution of dark matter can be present away from the cluster core.

preprint2013arXiv

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

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.

preprint2013arXiv

The Comparative Chemical Evolution of an Isolated Dwarf Galaxy: A VLT and Keck Spectroscopic Survey of WLM

Building on our previous spectroscopic and photometric analysis of the isolated Local Group dwarf irregular (dIrr) galaxy WLM, we present a comparison of the metallicities of its RGB stars with respect to the well studied Local Group dwarf spheroidal galaxies (dSphs) and Magellanic Clouds. We calculate a mean metallicity of [Fe/H]$ = -1.28 \pm 0.02$, and intrinsic spread in metallicity of $σ= 0.38 \pm 0.04$ dex, similar to the mean and spread observed in the massive dSph Fornax and the Small Magellanic Cloud. Thus, despite its isolated environment the global metallicity still follows expectations for WLM's mass and its global chemical evolution is similar to other nearby luminous dwarf galaxies (gas-rich or gas-poor). The data also show a radial gradient in [Fe/H] of $\rm{d[Fe/H]/dr_{c}} = -0.04 \pm 0.04$ dex $\rm{r_{c}^{-1}}$, which is flatter than that seen in the unbiased and spatially extended surveys of dSphs. Comparison of the spatial distribution of [Fe/H] in WLM, the Magellanic Clouds, and a sample of Local Group dSphs, shows an apparent dichotomy in the sense that the dIrrs have statistically flatter radial [Fe/H] gradients than the low angular momentum dSphs. The correlation between angular momentum and radial metallicity gradient is further supported when considering the Local Group dEs. This chemodynamic relationship offers a new and useful constraint for environment driven dwarf galaxy evolution models in the Local Group.

preprint2013arXiv

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

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.

preprint2013arXiv

The masses of Local Group dwarf spheroidal galaxies: The death of the universal mass profile

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

preprint2013arXiv

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

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.

preprint2013arXiv

The Three-Dimensional Structure of the M31 Satellite System; Strong Evidence for an Inhomogeneous Distribution of Satellites

We undertake an investigation into the spatial structure of the M31 satellite system utilizing the distance distributions presented in a previous publication. These distances make use of the unique combination of depth and spatial coverage of the Pan-Andromeda Archaeological Survey (PAndAS) to provide a large, homogeneous sample consisting of 27 of M31's satellites, as well as M31 itself. We find that the satellite distribution, when viewed as a whole, is no more planar than one would expect from a random distribution of equal size. A disk consisting of 15 of the satellites is however found to be highly significant, and strikingly thin, with a root-mean-square thickness of just $12.34^{+0.75}_{-0.43}$ kpc. This disk is oriented approximately edge on with respect to the Milky Way and almost perpendicular to the Milky Way disk. It is also roughly orthogonal to the disk like structure regularly reported for the Milky Way satellite system and in close alignment with M31's Giant Stellar Stream. A similar analysis of the asymmetry of the M31 satellite distribution finds that it is also significantly larger than one would expect from a random distribution. In particular, it is remarkable that 20 of the 27 satellites most likely lie on the Milky Way side of the galaxy, with the asymmetry being most pronounced within the satellite subset forming the aforementioned disk. This lopsidedness is all the more intriguing in light of the apparent orthogonality observed between the satellite disk structures of the Milky Way and M31.

preprint2012arXiv

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

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...

preprint2012arXiv

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

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.

preprint2012arXiv

The Resolved Structure and Dynamics of an Isolated Dwarf Galaxy: A VLT and Keck Spectroscopic Survey of WLM

We present spectroscopic data for 180 red giant branch stars in the isolated dwarf irregular galaxy WLM. Observations of the Calcium II triplet lines in spectra of RGB stars covering the entire galaxy were obtained with FORS2 at the VLT and DEIMOS on Keck II allowing us to derive velocities, metallicities, and ages for the stars. With accompanying photometric and radio data we have measured the structural parameters of the stellar and gaseous populations over the full galaxy. The stellar populations show an intrinsically thick configuration with $0.39 \leq q_{0} \leq 0.57$. The stellar rotation in WLM is measured to be $17 \pm 1$ km s$^{-1}$, however the ratio of rotation to pressure support for the stars is $V/σ\sim 1$, in contrast to the gas whose ratio is seven times larger. This, along with the structural data and alignment of the kinematic and photometric axes, suggests we are viewing WLM as a highly inclined oblate spheroid. Stellar rotation curves, corrected for asymmetric drift, are used to compute a dynamical mass of $4.3\pm 0.3\times10^{8} $M$_{\odot}$ at the half light radius ($r_{h} = 1656 \pm 49$ pc). The stellar velocity dispersion increases with stellar age in a manner consistent with giant molecular cloud and substructure interactions producing the heating in WLM. Coupled with WLM's isolation, this suggests that the extended vertical structure of its stellar and gaseous components and increase in stellar velocity dispersion with age are due to internal feedback, rather than tidally driven evolution. These represent some of the first observational results from an isolated Local Group dwarf galaxy which can offer important constraints on how strongly internal feedback and secular processes modulate SF and dynamical evolution in low mass isolated objects.

preprint2011arXiv

AAOmega spectroscopy of 29 351 stars in fields centered on ten Galactic globular clusters

Galactic globular clusters have been pivotal in our understanding of many astrophysical phenomena. Here we publish the extracted stellar parameters from a recent large spectroscopic survey of ten globular clusters. A brief review of the project is also presented. Stellar parameters have been extracted from individual stellar spectra using both a modified version of the Radial Velocity Experiment (RAVE) pipeline and a pipeline based on the parameter estimation method of RAVE. We publish here all parameters extracted from both pipelines. We calibrate the metallicity and convert this to [Fe/H] for each star and, furthermore, we compare the velocities and velocity dispersions of the Galactic stars in each field to the Besançon Galaxy model. We find that the model does not correspond well with the data, indicating that the model is probably of little use for comparisons with pencil beam survey data such as this.

preprint2011arXiv

Density Variations in the NW Star Stream of M31

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.

preprint2011arXiv

PAndAS' progeny: extending the M31 dwarf galaxy cabal

We present the discovery of five new dwarf galaxies, Andromeda XXIII-XXVII, located in the outer halo of M31. These galaxies were found in the second year of data from the Pan-Andromeda Archaeological Survey (PAndAS) of the M31/M33 subgroup. This survey now provides an almost complete panoramic view of the M31 halo out to an average projected radius of ~150 kpc. The metal-poor stellar density map for this whole region serves, not only as an illustration of the discovery space for satellite galaxies, but also gives a birds-eye view of the ongoing assembly process of an L* disk galaxy. Four of the new dwarfs appear as well-defined spatial over-densities of stars lying on the expected locus of metal-poor red giant branch stars at the distance of M31. The fifth over-density, And XXVII, is embedded in an extensive stream of such stars and is possibly the remnant of a strong tidal disruption event. All five satellites have metallicities and luminosities typical of dwarf spheroidal galaxies and continue the trend whereby the brighter dwarf spheroidal satellites of M31 generally have much larger half-light radii than their Milky Way counterparts. With an extended sample of M31 satellite galaxies we also revisit the spatial distribution of this population and find that, within the current limits of the PAndAS survey, the surface density of satellites is essentially constant out to 150 kpc. This corresponds to a radial density distribution of satellites varying as 1/r, a result seemingly in conflict with the predictions of cosmological simulations. (Abridged)

preprint2011arXiv

Stellar Streams as Probes of Dark Halo Mass and Morphology: A Bayesian Reconstruction

Tidal streams provide a powerful tool by means of which the matter distribution of the dark matter halos of their host galaxies can be studied. However, the analysis is not straightforward because streams do not delineate orbits, and for most streams, especially those in external galaxies, kinematic information is absent. We present a method wherein streams are fit with simple corrections made to possible orbits of the progenitor, using a Bayesian technique known as Parallel Tempering to efficiently explore the parameter space. We show that it is possible to constrain the shape of the host halo potential or its density distribution using only the projection of tidal streams on the sky, if the host halo is considered to be axisymmetric. By adding kinematic data or the circular velocity curve of the host to the fitting data, we are able to recover other parameters of the matter distribution such as its mass and profile. We test our method on several simulated low mass stellar streams and also explore the cases for which additional data are required.

preprint2011arXiv

The Star Formation History and Dust Content in the Far Outer Disc of M31

We present a detailed analysis of two fields located 26 kpc (~5 scalelengths) from the centre of M31. One field samples the major axis populations--the Outer Disc field--while the other is offset by ~18' and samples the Warp in the stellar disc. The CMDs based on HST/ACS imaging reach old main-sequence turn-offs (~12.5 Gyr). We apply the CMD-fitting technique to the Warp field to reconstruct the star formation history (SFH). We find that after undergoing roughly constant SF until about 4.5 Gyr ago, there was a rapid decline in activity and then a ~1.5 Gyr lull, followed by a strong burst lasting 1.5 Gyr and responsible for 25% of the total stellar mass in this field. This burst appears to be accompanied by a decline in metallicity which could be a signature of the inflow of metal-poor gas. The onset of the burst (~3 Gyr ago) corresponds to the last close passage of M31 and M33 as predicted by detailed N-body modelling, and may have been triggered by this event. We reprocess the deep M33 outer disc field data of Barker et al. (2011) in order to compare consistently-derived SFHs. This reveals a similar duration burst that is exactly coeval with that seen in the M31 Warp field, lending further support to the interaction hypothesis. The complex SFHs and the smoothly-varying age-metallicity relations suggest that the stellar populations observed in the far outer discs of both galaxies have largely formed in situ rather than migrated from smaller galactocentric radii. The strong differential reddening affecting the CMD of the Outer Disc field prevents derivation of the SFH. Instead, we quantify this reddening and find that the fine-scale distribution of dust precisely follows that of the HI gas. This indicates that the outer HI disc of M31 contains a substantial amount of dust and therefore suggests significant metal enrichment in these parts, consistent with inferences from our CMD analysis.

preprint2010arXiv

AAOmega Observations of 47 Tucanae: Evidence for a Past Merger?

The globular cluster 47 Tucanae is well studied but it has many characteristics that are unexplained, including a significant rise in the velocity dispersion profile at large radii, indicating the exciting possibility of two distinct kinematic populations. In this Letter we employ a Bayesian approach to the analysis of the largest available spectral dataset of 47 Tucanae to determine whether this apparently two-component population is real. Assuming the two models were equally likely before taking the data into account, we find that the evidence favours the two-component population model by a factor of ~3x10^7. Several possible explanations for this result are explored, namely the evaporation of low-mass stars, a hierarchical merger, extant remnants of two initially segregated populations, and multiple star formation epochs. We find the most compelling explanation for the two-component velocity distribution is that 47 Tuc formed as two separate populations arising from the same proto-cluster cloud which merged <7.3 +/- 1.5 Gyr ago. This may also explain the extreme rotation, low mass-to-light ratio and mixed stellar populations of this cluster.

preprint2010arXiv

Halo globular clusters observed with AAOmega: dark matter content, metallicity and tidal heating

Globular clusters have proven to be essential to our understanding of many important astrophysical phenomena. Here we analyse spectroscopic observations of ten Halo globular clusters to determine their dark matter content, their tidal heating by the Galactic disc and halo, describe their metallicities and the likelihood that Newtonian dynamics explain their kinematics. We analyse a large number of members in all clusters, allowing us to address all these issues together, and we have included NGC 288 and M30 to overlap with previous studies. We find that any flattening of the velocity dispersion profiles in the outer regions of our clusters can be explained by tidal heating. We also find that all our GCs have M/L_V < 5, therefore, we infer the observed dynamics do not require dark matter, or a modification of gravity. We suggest that the lack of tidal heating signatures in distant clusters indicates the Halo is not triaxial. The isothermal rotations of each cluster are measured, with M4 and NGC 288 exhibiting rotation at a level of 0.9 +/- 0.1 km/s and 0.25 +/- 0.15 km/s, respectively. We also indirectly measure the tidal radius of NGC 6752, determining a more realistic figure for this cluster than current literature values. Lastly, an unresolved and intriguing puzzle is uncovered with regard to the cooling of the outer regions of all ten clusters.

preprint2009arXiv

PAndAS' cubs: discovery of two new dwarf galaxies in the surroundings of the Andromeda and Triangulum galaxies

We present the discovery of two new dwarf galaxies, Andromeda XXI and Andromeda XXII, located in the surroundings of the Andromeda and Triangulum galaxies (M31 and M33). These discoveries stem from the first year data of the Pan-Andromeda Archaeological Survey (PAndAS), a photometric survey of the M31/M33 group conducted with the Megaprime/MegaCam wide-field camera mounted on the Canada-France-Hawaii Telescope. Both satellites appear as spatial overdensities of stars which, when plotted in a color-magnitude diagram, follow metal-poor, [Fe/H]=-1.8, red giant branches at the distance of M31/M33. Andromeda XXI is a moderately bright dwarf galaxy (M_V=-9.9+/-0.6), albeit with low surface brightness, emphasizing again that many relatively luminous M31 satellites still remain to be discovered. It is also a large satellite, with a half-light radius close to 1 kpc, making it the fourth largest Local Group dwarf spheroidal galaxy after the recently discovered Andromeda XIX, Andromeda II and Sagittarius around the Milky Way, and supports the trend that M31 satellites are larger than their Milky Way counterparts. Andromeda XXII is much fainter (M_V=-6.5+/-0.8) and lies a lot closer in projection to M33 than it does to M31 (42 vs. 224 kpc), suggesting that it could be the first Triangulum satellite to be discovered. Although this is a very exciting possibility in the context of a past interaction of M33 with M31 and the fate of its satellite system, a confirmation will have to await a good distance estimate to confirm its physical proximity to M33. Along with the dwarf galaxies found in previous surveys of the M31 surroundings, these two new satellites bring the number of dwarf spheroidal galaxies in this region to 20.

preprint2009arXiv

Testing Newtonian Gravity with AAOmega: Mass-To-Light Profiles and Metallicity Calibrations From 47 Tuc and M55

Globular clusters are an important test bed for Newtonian gravity in the weak-acceleration regime, which is vital to our understanding of the nature of the gravitational interaction. Recent claims have been made that the velocity dispersion profiles of globular clusters flatten out at large radii, despite an apparent paucity of dark matter in such objects, indicating the need for a modification of gravitational theories. We continue our investigation of this claim, with the largest spectral samples ever obtained of 47 Tucanae and M55. Furthermore, this large sample allows for an accurate metallicity calibration based on the equivalent widths of the calcium triplet lines and K band magnitude of the Tip of the Red Giant Branch. Assuming an isothermal distribution, the rotations of each cluster are also measured with both clusters exhibiting clear rotation signatures. The global velocity dispersions of NGC 121 and Kron 3, two globular clusters in the Small Magellanic Cloud, are also calculated. By applying a simple dynamical model to the velocity dispersion profiles of 47 Tuc and M55, we calculate their mass-to-light profiles, total masses and central velocity dispersions. We find no statistically significant flattening of the velocity dispersion at large radii for M55, and a marked increase in the profile of 47 Tuc for radii greater than approximately half the tidal radius. We interpret this increase as an evaporation signature, indicating that 47 Tuc is undergoing, or has undergone, core-collapse, but find no requirement for dark matter or a modification of gravitational theories in either cluster.

preprint2009arXiv

Testing Newtonian Gravity with AAOmega: Mass-to-Light Profiles of Four Globular Clusters

Testing Newtonian gravity in the weak-acceleration regime is vital to our understanding of the nature of the gravitational interaction. It has recently been claimed that the velocity dispersion profiles of several globular clusters flatten out at large radii, reminiscent of galaxy rotation curves, even though globular clusters are thought to contain little or no dark matter. We investigate this claim, using AAOmega observations of four globular clusters, namely M22, M30, M53 and M68. M30, one such cluster that has had this claim made for its velocity dispersion, was included for comparison with previous studies. We find no statistically significant flattening of the velocity dispersion at large radii for any of our target clusters and therefore we infer the observed dynamics do not require that globular clusters are dark matter dominated, or a modification of gravity. Furthermore, by applying a simple dynamical model we determine the radial mass-to-light profiles for each cluster. The isothermal rotations of each cluster are also measured, with M22 exhibiting clear rotation, M68 possible rotation and M30 and M53 lacking any rotation, within the uncertainties.

preprint2009arXiv

The remnants of galaxy formation from a panoramic survey of the region around M31

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

preprint2002arXiv

An investigation of gravitational lens determinations of H_o in quintessence cosmologies

There is growing evidence that the majority of the energy density of the universe is not baryonic or dark matter, rather it resides in an exotic component with negative pressure. The nature of this `quintessence' influences our view of the universe, modifying angular diameter and luminosity distances. Here, we examine the influence of a quintessence component upon gravitational lens time delays. As well as a static quintessence component, an evolving equation of state is also considered. It is found that the equation of state of the quintessence component and its evolution influence the value of the Hubble's constant derived from gravitational lenses. However, the differences between evolving and non-evolving cosmologies are relatively small. We undertake a suite of Monte Carlo simulations to examine the potential constraints that can be placed on the universal equation of state from the monitoring of gravitational lens system, and demonstrate that at least an order of magnitude more lenses than currently known will have to be discovered and analysed to accurately probe any quintessence component.

preprint2002arXiv

Cosmology on a Mesh

An adaptive multi grid approach to simulating the formation of structure from collisionless dark matter is described. MLAPM (Multi-Level Adaptive Particle Mesh) is one of the most efficient serial codes available on the cosmological 'market' today. As part of Swinburne University's role in the development of the Square Kilometer Array, we are implementing hydrodynamics, feedback, and radiative transfer within the MLAPM adaptive mesh, in order to simulate baryonic processes relevant to the interstellar and intergalactic media at high redshift. We will outline our progress to date in applying the existing MLAPM to a study of the decay of satellite galaxies within massive host potentials.

preprint2001arXiv

Was SN1997ff at z~1.7 magnified by gravitational lensing?

The quest for the cosmological parameters has come to fruition with the identification of a number of supernovae at a redshift of $z\sim1$. Analyses of the brightness of these standard candles reveal that the Universe is dominated by a large cosmological constant. The recent identification of the $z\sim1.7$ SN1997ff in the northern Hubble Deep Field has provided further evidence for this cosmology. Here we examine the case for gravitational lensing of SN1997ff due to the presence of galaxies lying along our line of sight. We find that, while the alignment of SN1997ff with foreground masses was not favorable for it to be multiply imaged and strongly magnified, two galaxies did lie close enough to result in significant magnification: $μ\sim1.4$ for the case where these elliptical galaxies have velocity dispersion $200 {\rm km/s}$. Given the small difference between supernova brightnesses in different cosmologies, detailed modeling of the gravitational lensing properties of the intervening matter is therefore required before the true cosmological significance of SN1997ff can be deduced.

preprint1998arXiv

Archer of the Galactic disk? The effect on the outer HI disk of the Milky Way of collisional encounters with the Sagittarius dwarf galaxy

Hydrodynamical calculations undertaken to simulate the collisional interaction between the Sgr dwarf galaxy and the Galactic outer HI disk are presented, constrained by recently derived orbital and mass parameters for this dwarf galaxy. It is found that a significant distortion to the structure of the Galactic HI disk will be induced by the collision if the mass of the dwarf exceeds ~10^9 M_Sun; this value is consistent with an estimate derived by requiring that the dwarf galaxy is sufficiently robust to survive tidal disruption until the present time. Though the precise details of the interaction are compromised in our simulations by the lack of a live Galactic halo, we find that for model masses > 5x10^9 M_Sun, prominent spiral arms and a substantial lopsidedness in the outer disk are produced. Furthermore, a noticeable warp-like structure is induced in the disk. Thus the Sgr dwarf may have significantly affected the star formation history and structure of the outer Galaxy. These simulations confirm the possibility of determining the current merging rate of low surface brightness, gas-poor dwarf galaxies of mass > 10^9 M_Sun onto giant spiral galaxies from careful analysis of observations of the structure of HI disks.