Search

We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins., Comment: Submitted to ApJ, 23 pages, 13 figures 4 tables

The dark matter subhalos orbiting in a galactic halo perturb the orbits of stars in thin stellar streams. Over time the random velocities in the streams develop non-Gaussian wings. The rate of velocity increase is approximately a random walk at a rate proportional to the number of subhalos, primarily those in the mass range $\approx 10^{6-7} M_\odot$. The distribution of random velocities in long, thin, streams is measured in simulated Milky Way-like halos that develop in representative WDM and CDM cosmologies. The radial velocity distributions are well modeled as the sum of a Gaussian and an exponential. The resulting MCMC fits find Gaussian cores of 1-2 km/sec and exponential wings that increase from 3 km/sec for 5.5 keV WDM, 4 km/sec for 7 keV WDM, to 6 km/sec for a CDM halo. The observational prospects to use stream measurements to constrain the nature of galactic dark matter are discussed., Comment: ApJ accepted

Context. Gaia DR3 has offered the scientific community a remarkable dataset of approximately one million spectra acquired with the Radial Velocity Spectrometer (RVS) in the Calcium II triplet region, that is well-suited to identify very metal-poor (VMP) stars. However, over 40% of these spectra have no released parameters by Gaia's GSP Spec pipeline in the domain of VMP stars, whereas VMP stars are key tracers of early Galactic evolution. Aims. We aim to provide spectroscopic metallicities for VMP stars using Gaia RVS spectra, thereby producing a catalogue of bright VMP stars distributed over the full sky that can serve as the basis to study early chemical evolution throughout the Galaxy. Methods. We select VMP stars using photometric metallicities from the literature and analyse the Gaia RVS spectra to infer spectroscopic metallicities for these stars. Results. The inferred metallicities agree very well with literature high-resolution metallicities with a median systematic offset of 0.1 dex and standard deviation of $\sim$0.15 dex. The purity of this sample in the VMP regime is $\sim$80% with outliers representing a mere $\sim$3%. Conclusions. We make available an all-sky catalogue of $\sim$1500 stars with reliable spectroscopic metallicities down to [Fe/H]$\sim$-4.0, of which $\sim$1000 are VMP stars. More than 75% of these stars have either no metallicity value in the literature to date or are flagged to be unreliable in their literature metallicity estimates. This catalogue of bright (G<13) VMP stars is three times larger than the current sample of well-studied VMP stars in the literature in this magnitude range, making it ideal for high-resolution spectroscopic follow-up and to study the properties of VMP stars in different parts of our Galaxy., Comment: Submitted to A&A. 11 pages, 4 pages appendix, 5 figures. Catalogue is made available public here: https://astroakshara.github.io/rvs-paper/Gaia-RVS-VMP-catalogue-AV23b-vSep23.csv Feel free to use it to do amazing science with the amazing data Gaia has given us

Dark matter sub-halos that pass near or through a thin tidal star stream locally increase its velocity dispersion. Subsequent orbital evolution further increases the velocity dispersion and stream width, lowering the surface density of a stream. The kinematic properties of streams are measured in cosmological Milky Way-like halo simulations. The distance along a stream is a proxy for the time a star has spent in the stream, although there are a range of ages at any distance. Power law fits to the velocity dispersion with angular distance for the average of the streams in the 10-60 kpc range finds sigma_theta=6 phi^{0.25} km/s, sigma_phi=8 phi^{0.39} km/s, and sigma_r=10 phi^{0.44} km/s for |phi|< 34 degrees, for stars within theta=+/-5 degrees of the stream equator. The errors of the coefficients are about 10% for these streams, with comparable systematic errors depending on exactly which streams are selected and the stream width and length fitted. The stream velocity dispersions close to the clusters generally increase with the sub-halo numbers., Comment: accepted for publication in the ApJ

The investigation of the metal-poor tail in the Galactic bulge provides unique information on the early Milky Way assembly and evolution. A chemo-dynamical analysis of 17 very metal-poor stars (VMP, [Fe/H] $<-2.0$) selected from the Pristine Inner Galaxy Survey was carried out based on Gemini/GRACES spectra. The chemistry suggests that the majority of our stars are very similar to metal-poor stars in the Galactic halo. Orbits calculated from {\it Gaia} EDR3 imply these stars are brought into the bulge during the earliest Galactic assembly. Most of our stars have large [Na,Ca/Mg] abundances, and thus show little evidence of enrichment by pair-instability supernovae. Two of our stars (P171457, P184700) have chemical abundances compatible with second-generation globular cluster stars, suggestive of the presence of ancient and now dissolved globular clusters in the inner Galaxy. One of them (P171457) is extremely metal-poor ([Fe/H] $<-3.0$) and well below the metallicity floor of globular clusters, which supports the growing evidence for the existence of lower-metallicity globular clusters in the early Universe. A third star (P180956, [Fe/H] $\sim-2$) has low [Na,Ca/Mg] and very low [Ba/Fe] for its metallicity, which are consistent with formation in a system polluted by only one or a few low-mass supernovae. Interestingly, its orbit is confined to the Galactic plane, like other very metal-poor stars found in the literature, which have been associated with the earliest building blocks of the Milky Way., Comment: 22 pages + appendix, 14 Figures + 3 in the appendix, Accepted for publication in MNRAS

SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is intermediate between two subtypes of normal Type Ia: Core-Normal and Broad-Line. The excess emission could have contributions from the radioactive decay of surface iron-peak elements as well as ejecta interaction with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H$\alpha$ and He~I favour a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O~I] and He~I disfavours a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN~2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe~II] and [Ni~II]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1-D models are incompatible with the infant-phase excess emission, $B_{\rm max}-V_{\rm max}$ color, and absence of nebular-phase [Ca~II]. Although the explosion processes of SN~2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear normal after $\sim$ 1 day., Comment: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

Using {\it Herschel}-SPIRE imaging and the Canada-France Imaging Survey (CFIS) Low Surface Brightness data products from the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS), we present a cross-correlation between the cosmic far-infrared background and cosmic optical background fluctuations. The cross-spectrum is measured for two cases: all galaxies are kept in the images; or all individually-detected galaxies are masked to produce `background' maps. We report the detection of the cross-correlation signal at $\gtrsim 18\,\sigma$ ($\gtrsim 14\,\sigma$ for the background map). The part of the optical brightness variations that are correlated with the submm emission translates to an rms brightness of $\simeq 32.5\,{\rm mag}\,{\rm arcsec}^{-2}$ in the $r$ band, a level normally unreachable for individual sources. A critical issue is determining what fraction of the cross-power spectrum might be caused by emission from Galactic cirrus. For one of the fields, the Galactic contamination is 10 times higher than the extragalactic signal; however, for the other fields, the contamination is around 20~per cent. An additional discriminant is that the cross-power spectrum is of the approximate form $P(k)\propto 1/k$, much shallower than that of Galactic cirrus. We interpret the results in a halo-model framework, which shows good agreement with independent measurements for the scalings of star-formation rates in galaxies. The approach presented in this study holds great promise for future surveys such as FYST/CCAT-prime combined with {\it Euclid} or the Vera Rubin Observatory (LSST), which will enable a detailed exploration of the evolution of star formation in galaxies., Comment: 40 pages, 26 figures, published in MNRAS

M54 is a prototype of a globular cluster embedded in a dark matter halo. Gaia EDR3 photometry and proper motions separate the old, metal-poor stars from the more metal rich and younger dwarf galaxy stars. The metal poor stars dominate the inner 50 pc, with a velocity dispersion profile that declines to a minimum around 30 pc then rises back to nearly the central velocity dispersion, as expected for a globular cluster at the center of a dark matter halo. The Jeans analysis of the three separate stellar populations give consistent masses that rise approximately linearly with radius to 1 kpc, implying a small core or cuspy halo. These data are compatible with an infalling CDM dark matter halo reduced to 3x10^8 M_sun at the 50 kpc apocenter 2.3 Gyr ago, with a central globular cluster surrounded by the remnant of a dwarf galaxy. Tides gradually remove material beyond 1 kpc but have little effect on the stars and dark matter within 300 pc of the center. M54 appears to be a transitional system between globular clusters with and without local dark halos, whose evolution within the galaxy depends sensitively on the time of accretion and orbital pericenter., Comment: ApJ accepted

Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions., Comment: Accepted for publication in Nature Astronomy. Main text = 8 pages, 4 figures, 1 table; Full document = 46 pages, with Methods, Supplementary Information, 7 Supplementary figures, 2 Supplementary tables and references. Nat Astron (2022)

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., Comment: Published in ApJ. Key results in Table 3. The paper contains 38 pages, 17 figures, 4 Tables

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., Comment: 14 pages, 6 figures, 5 tables. Nature, accepted

We present a new spectroscopic study of the dwarf galaxy Bootes I (Boo I) with data from the Anglo-Australian Telescope and its AAOmega spectrograph together with the Two Degree Field multi-object system. We observed 36 high-probability Boo I stars selected using Gaia Early Data Release 3 proper motions and photometric metallicities from the Pristine survey. Out of those, 27 are found to be Boo I's stars, resulting in an excellent success rate of 75% at finding new members. Our analysis uses a new pipeline developed to estimate radial velocities and equivalent widths of the calcium triplet lines from Gaussian and Voigt line profile fits. The metallicities of 16 members are derived, including 3 extremely metal-poor stars ([Fe/H] < -3.0), which translates into a success rate of 25% at finding them with the combination of Pristine and Gaia. Using the large spatial extent of our new members that spans up to 4.1 half-light radii and spectroscopy from the literature, we find a systemic velocity gradient of 0.40 +/- 0.10 km/s/arcmin and a small but resolved metallicity gradient of -0.008 +/- 0.003 dex/arcmin. Finally, we show that Boo I is more elongated than previously thought with an ellipticity of epsilon = 0.68 +/- 0.15. Its velocity and metallicity gradients as well as its elongation suggest that Boo I may have been affected by tides, a result supported by direct dynamical modelling.

We evaluate the rotational velocity of stars observed by the Pristine survey towards the Galactic anticentre, spanning a wide range of metallicities from the extremely metal-poor regime ($\mathrm{[Fe/H]}<-3$) to nearly solar metallicity. In the Galactic anticentre direction, the rotational velocity ($V_{\phi}$) is similar to the tangential velocity in the galactic longitude direction ($V_{\ell}$). This allows us to estimate $V_{\phi}$ from Gaia early data-release 3 (Gaia EDR3) proper motions for stars without radial velocity measurements. This substantially increases the sample of stars in the outer disc with estimated rotational velocities. Our stellar sample towards the anticentre is dominated by a kinematical thin disc with a mean rotation of $\sim -220$ km $\mathrm{s}^{-1}$. However, our analysis reveals the presence of more stellar substructures. The most intriguing is a well populated extension of the kinematical thin disc down to $\mathrm{[Fe/H]} \sim -2$. A scarser fast rotating population reaching the extremely metal-poor regime, down to $\mathrm{[Fe/H]} \sim -3.5$ is also detected, but without statistical significance to unambiguously state whether this is the extremely metal-poor extension of the thin disc or the high rotating tail of hotter structures (like the thick disc or the halo). In addition, a more slowly rotating kinematical thick disc component is also required to explain the observed $V_{\ell}$ distribution at $\mathrm{[Fe/H]} > -1.5$. Furthermore, we detect signatures of a "heated disc", the so-called Splash, at metallicities higher than $\sim-1.5$. Finally, at $\mathrm{[Fe/H]} < -1.5$ our anticentre sample is dominated by a kinematical halo with a net prograde motion., Comment: 18 pages, 14 figures. Accepted for publication in MNRAS

The proper motions of stars in the outskirts of globular clusters are used to estimate cluster velocity dispersion profiles as far as possible within their tidal radii. We use individual color-magnitude diagrams to select high probability cluster stars for 25 metal-poor globular clusters within 20 kpc of the sun, 19 of which have substantial numbers of stars at large radii. Of the 19, 11 clusters have a falling velocity dispersion in the 3-6 half mass radii range, 6 are flat, and 2 plausibly have a rising velocity dispersion. The profiles are all in the range expected from simulated clusters started at high redshift in a zoom-in cosmological simulation. The 11 clusters with falling velocity dispersion profiles are consistent with no dark matter above the Galactic background. The 6 clusters with approximately flat velocity dispersion profiles could have local dark matter, but are ambiguous. The 2 clusters with rising velocity dispersion profiles are consistent with a remnant local dark matter halo, but need membership confirmation and detailed orbital modeling to further test these preliminary results., Comment: AAS revised version submitted

A cosmological zoom-in simulation which develops into a Milky Way-like halo is started at redshift 7. The initial dark matter distribution is seeded with dense star clusters, median mass $5\times 10^5 M_\sun$, placed in the largest sub-halos present, which have a median peak circular velocity of 25 \kms. Three simulations are initialized using the same dark matter distribution, with the star clusters started on approximately circular orbits having initial median radii 6.8 kpc, 0.14 kpc, and, at the exact center of the sub-halos. The simulations are evolved to the current epoch at which time the median galactic orbital radii of the three sets of clusters are 30, 5 and 16 kpc, with the clusters losing about 2, 50 and 15\% of their mass, respectively. Clusters started at small orbital radii have so much tidal forcing that they are often not in equilibrium. Clusters started at larger sub-halo radii have a velocity dispersion that declines smoothly to $\simeq$20\% of the central value at $\simeq$20 half mass radii. The clusters started at the sub-halo centers can show a rise in velocity dispersion beyond 3-5 half mass radii. That is, the clusters formed without local dark matter always have stellar mass dominated kinematics at all radii, whereas about 25\% of the clusters started at sub-halo centers have remnant local dark matter., Comment: ApJ revised submission

The currently known Galactic globular cluster population extends out to a maximum galactocentric distance of $\sim$ 145 kpc, with the peculiarity that the outermost clusters predominantly have an inward velocity. Orbit averaging finds that this configuration occurs by chance about $6\%$ of the time, suggesting that several globular clusters with positive radial velocities remain undiscovered. We evaluate the expected number of undiscovered clusters at large distances under the assumption that the cluster population has a smooth radial distribution and is in equilibrium within the Milky Way's virial radius. By comparing the present day kinematic properties of outer clusters to random orbital configurations of the Galactic globular cluster system through orbit averaging, we estimate a likelihood of $73\%$ of there being at least one undiscovered globular cluster within the Milky Way. This estimate assumes the current population is complete out to 50 kpc, and increases to $91\%$ if the population is complete out to 150 kpc. The likelihood of there being two undiscovered clusters is between $60\%$ and $70\%$, with the likelihood of there being three undiscovered clusters being on the order of $50\%$. The most likely scenario is that the undiscovered clusters are moving outwards, which results in the outer cluster population being consistent with an equilibrium state. Searches for distant and possibly quite low concentration and very low metallicity globular clusters will be enabled with upcoming deep imaging surveys., Comment: 13 pages, 9 figures, 4 tables, resubmitted to MNRAS for publication

Metal-poor stars are important tools for tracing the early history of the Milky Way, and for learning about the first generations of stars. Simulations suggest that the oldest metal-poor stars are to be found in the inner Galaxy. Typical bulge surveys, however, lack low metallicity ([Fe/H] < -1.0) stars because the inner Galaxy is predominantly metal-rich. The aim of the Pristine Inner Galaxy Survey (PIGS) is to study the metal-poor and very metal-poor (VMP, [Fe/H] < -2.0) stars in this region. In PIGS, metal-poor targets for spectroscopic follow-up are selected from metallicity-sensitive CaHK photometry from the CFHT. This work presents the ~250 deg^2 photometric survey as well as intermediate-resolution spectroscopic follow-up observations for ~8000 stars using AAOmega on the AAT. The spectra are analysed using two independent tools: ULySS with an empirical spectral library, and FERRE with a library of synthetic spectra. The comparison between the two methods enables a robust determination of the stellar parameters and their uncertainties. We present a sample of 1300 VMP stars -- the largest sample of VMP stars in the inner Galaxy to date. Additionally, our spectroscopic dataset includes ~1700 horizontal branch stars, which are useful metal-poor standard candles. We furthermore show that PIGS photometry selects VMP stars with unprecedented efficiency: 86%/80% (lower/higher extinction) of the best candidates satisfy [Fe/H] < -2.0, as do 80%/63% of a larger, less strictly selected sample. We discuss future applications of this unique dataset that will further our understanding of the chemical and dynamical evolution of the innermost regions of our Galaxy., Comment: accepted for publication in MNRAS, 17 pages, 9 figures

Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for tens of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles.

We present a new spectroscopic study of the faint Milky Way satellite Sagittarius II. Using multi-object spectroscopy from the Fibre Large Array Multi Element Spectrograph, we supplement the dataset of Longeard et al. (2020) with 47 newly observed stars, 19 of which are identified as members of the satellite. These additional member stars are used to put tighter constraints on the dynamics and the metallicity properties of the system. We find a low velocity dispersion of SgrII v = 1.7 +/- 0.5 km s-1, in agreement with the dispersion of Milky Way globular clusters of similar luminosity. We confirm the very metal-poor nature of the satellite ([Fe/H]_SgrII = -2.23 +/- 0.07) and find that the metallicity dispersion of Sgr II is not resolved, reaching only 0.20 at the 95% confidence limit. No star with a metallicity below -2.5 is confidently detected. Therefore, despite the unusually large size of the system (rh = 35.5 +1.4-1.2 pc), we conclude that Sgr II is an old and metal-poor globular cluster of the Milky Way.

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events., Comment: 23 pages, 13 figures including appendices; accepted to ApJ with minor revisions

Tidally limited star clusters are started on nearly circular orbits in the dark matter sub-halos present at high redshift and evolved with an n-body code augmented with gravitational interactions in the clusters. The tidally removed stars and the remaining high redshift clusters from a redshift 8 start are more concentrated than the dark matter, as expected. However, the subset of stars from the clusters that began in the lower mass sub-halos have a distribution somewhat more extended than the dark matter halo, with a mean galactic radius of about 60 kpc inside 150 kpc. The clusters from low mass sub-halos, those with a peak circular velocity of $ 12-18$ \kms, also produce most of the population's thin stellar streams. The dependence of the stellar population distribution on sub-halo mass is not seen in simulations that start clusters at lower redshift. The half mass radii of the clusters are set by the tidal fields of the initial cluster orbits, which leads to the average half mass radius decreasing with increased redshift of formation. Starting clusters at greater than redshift 4 leads to cluster half mass radii approximately as seen in the Milky Way, with lower redshifts leading to clusters that are too large. If extremely metal poor globular clusters were preferentially formed in high redshift, low mass, dwarf galaxies, then the clusters should be spread nearly uniformly over the dark halo, along with their tidally removed stars., Comment: ApJ accepted, edited

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\sigma$ 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., Comment: Submitted to MNRAS letter, 5 pages with 1 figure. Appendices A and B with 3 figures

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.

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., Comment: 6 pages, 3 figures, 1 Table, accepted for publication in ApJ

We use the Pristine survey CaHK narrow-band photometry, combined with the SDSS ugr photometry, to provide a cleaner sample of blue horizontal branch stars in the Galactic halo out to large distances. We demonstrate a completeness of 91% and a purity of 93% with respect to available spectroscopic classifications. We subsequently use our new clean sample of these standard candles to investigate the substructure in the Galactic halo over the Pristine footprint. Among other features, this allows for a careful tracing of multiple parts of the Sagittarius stream, providing a measurement independent from other tracers used and reaching larger distances. Moreover, we demonstrate with this clean and complete sample that the halo follows a density profile with a negative power-law slope of 3.5 - 4.0. As the relatively shallow SDSS u-band is the limiting factor in this technique, we foresee large potential for combining Pristine survey photometry with the much deeper u-band photometry from the Canada-France-Imaging Survey., Comment: 14 pages, 11 figures. Accepted for publication in MNRAS

We present direct $N$-body simulations of tidally filling 30,000 ${\rm M}_\odot$ star clusters orbiting between 10 kpc and 100 kpc in galaxies with a range of dark matter substructure properties. The time-dependent tidal force is determined based on the combined tidal tensor of the galaxy's smooth and clumpy dark matter components, the latter of which causes fluctuations in the tidal field that can heat clusters. The strength and duration of these fluctuations are sensitive to the local dark matter density, substructure fraction, sub-halo mass function, and the sub-halo mass-size relation. Based on the cold dark matter framework, we initially assume sub-halos are Hernquist spheres following a power-law mass function between $10^5$ and $10^{11} {\rm M}_\odot$ and find that tidal fluctuations are too weak and too short to affect star cluster evolution. Treating sub-halos as point masses, to explore how denser sub-halos affect clusters, we find that only sub-halos with masses greater than $10^{6} {\rm M}_\odot$ will cause cluster dissolution times to decrease. These interactions can also decrease the size of a cluster while increasing the velocity dispersion and tangential anisotropy in the outer regions via tidal heating. Hence increased fluctuations in the tidal tensor, especially fluctuations that are due to low-mass halos, do not necessarily translate into mass loss. We further conclude that the tidal approximation can be used to model cluster evolution in the tidal fields of cosmological simulations with a minimum cold dark matter sub-halo mass of $10^{6} {\rm M}_\odot$, as the effect of lower-mass sub-halos on star clusters is negligible., Comment: 17 pages, 10 figures, Accepted for publication in MNRAS

Tidally disrupted globular cluster streams are usually observed, and therefore perceived, as narrow, linear and one-dimensional structures in the 6D phase-space. Here we show that the GD-1 stellar stream ($\approx$ 30 pc wide), which is the tidal debris of a disrupted globular cluster, possesses a secondary diffuse and extended stellar component ($\approx$ 100 pc wide) around it, detected at >5$\sigma$ confidence level. Similar morphological properties are seen in synthetic streams that are produced from star clusters that are formed within dark matter sub-halos and then accrete onto a massive host galaxy. This lends credence to the idea that the progenitor of the highly retrograde GD-1 stream was originally formed outside of the Milky Way in a now defunct dark satellite galaxy. We deem that in future studies, this newly found $cocoon$ component may serve as a structural hallmark to distinguish between the in-situ and ex-situ (accreted) formed globular cluster streams., Comment: 11 pages, 10 figures, accepted for publication in The Astrophysical Journal

Collisionless dark matter galactic halos are expected to exhibit damped oscillations as a result of ongoing late time accretion. An n-body model of the cosmological assembly of a Milky Way-like halo is used to quantify the time dependence of its gravitational field. The simulation contains stellar streams whose incommensurate perpendicular velocities are found to have an approximately exponential distribution with a scale of 10-20\kms, depending on how the stars are selected, comparable to those reported for the Orphan stream. The fluctuations in the quadrupole moment of the dark matter halo are sufficient to largely explain the tangential velocities. If velocity measurements of a larger sample of Milky Way streams finds (or does not find) the expected distribution of transverse velocities it will lead to limits on the cross-section of self-interacting dark matter, in which kinetic viscosity can damp the oscillations more rapidly than the mixing processes of collisionless dark matter alone., Comment: ApJ revised

Simulations are run with and without a normal cold dark matter sub-halo population below 4x10^8 M_sun to examine the role of the lower mass sub-halos in the creation of density variations, "gaps", within thin tidal star streams. Dense star clusters composed of stellar mass particles are started within a redshift 3 sub-galactic dark matter halo distribution which assembles to create a Milky Way like halo. Hierarchical assembly causes the early-time tidal star streams from the clusters to later be dispersed as thick streams around the thin streams that develop in the main late-time halo as it settles towards equilibrium. The two simulations are set up to be essentially dynamically identical and result in about 140 star clusters within a radius of 50 kpc, of which 40 some have thin tidal streams longer than 40\degr\ with a separation of more than 20\degr\ from the cluster as viewed from the halo center. Averaged over about 20 separate snapshots, the fraction of streams with significant density variations inside 50 kpc is $0.58\pm0.02$ in the simulation with low mass sub-halos and $0.52\pm0.02$ in the simulation without low mass halos. That is, only $12\pm5$\% of the stream density variations can be attributed to the sub-halos below 4x10^8 M_sun. The dominant source of density variations in streams is the varying tidal fields along the complex trajectories of the clusters and their streams as the galaxy assembles., Comment: ApJ accepted. Fig1 movies available at https://www.astro.utoronto.ca/~carlberg/movie60_461.mp4 and movie10_461.mmp4 (+/-60 and 10 kpc boxes centered on star cluster)

We use Gaia DR2 astrometric and photometric data, published radial velocities and MESA models to infer distances, orbits, surface gravities, and effective temperatures for all ultra metal-poor stars ($\FeH<-4.0$ dex) available in the literature. Assuming that these stars are old ($>11\Gyr$) and that they are expected to belong to the Milky Way halo, we find that these 42 stars (18 dwarf stars and 24 giants or sub-giants) are currently within $\sim20\kpc$ of the Sun and that they map a wide variety of orbits. A large fraction of those stars remains confined to the inner parts of the halo and was likely formed or accreted early on in the history of the Milky Way, while others have larger apocentres ($>30\kpc$), hinting at later accretion from dwarf galaxies. Of particular interest, we find evidence that a significant fraction of all known UMP stars ($\sim26$\%) are on prograde orbits confined within $3\kpc$ of the Milky Way plane ($J_z < 100 \kms \kpc$). One intriguing interpretation is that these stars belonged to the massive building block(s) of the proto-Milky Way that formed the backbone of the Milky Way disc. Alternatively, they might have formed in the early disc and have been dynamically heated, or have been brought into the Milky Way by one or more accretion events whose orbit was dragged into the plane by dynamical friction before disruption. The combination of the exquisite Gaia DR2 data and surveys of the very metal-poor sky opens an exciting era in which we can trace the very early formation of the Milky Way., Comment: V2: Accepted in MNRAS; minor changes in Scenario 1 and references. New appendix B

We present a detailed study of the faint Milky Way satellite Draco II (Dra II) from deep CFHT/MegaCam broadband $g$ and $i$ photometry and narrow-band metallicity-sensitive CaHK observations, along with follow-up Keck II/DEIMOS multi-object spectroscopy. Forward modeling of the deep photometry allows us to refine the structural and photometric properties of Dra II: the distribution of stars in colour-magnitude space implies Dra II is old (13.5 $ \pm 0.5 $ Gyr), very metal poor, very faint ($ L_V = 180 ^{+124}_{-72} L_\odot $), and at a distance $d = 21.5 \pm 0.4$ kpc. The narrow-band, metallicity-sensitive CaHK Pristine photometry confirms this very low metallicity ([Fe/H]$ = -2.7 \pm 0.1$ dex). Even though our study benefits from a doubling of the spectroscopic sample size compared to previous investigations, the velocity dispersion of the system is still only marginally resolved ($\sigma_{vr}<5.9$ km s$^{-1}$ at the 95 per cent confidence level) and confirms that Dra II is a dynamically cold stellar system with a large recessional velocity ($\langle v_{r}\rangle = -342.5^{+1.1}_{-1.2} $ km s$^{-1}$). We further show that the spectroscopically confirmed members of Dra~II have a mean proper motion of $(\mu_\alpha^*,\mu_\delta)=(1.26 \pm 0.27,0.94 \pm 0.28) $ mas yr$^{-1}$ in the Gaia DR2 data, which translates to an orbit with a pericenter and an apocenter of $21.3 ^{+0.7}_{-1.0}$ and $153.8 ^{+56.7}_{-34.7}$ kpc, respectively. Taken altogether, these properties favour the scenario of Dra~II being a potentially disrupting dwarf galaxy. The low-significance extra-tidal features we map around the satellite tentatively support this scenario., Comment: Accepted for publication in MNRAS

The early Universe presented a star formation environment that was almost devoid of heavy elements. The lowest metallicity stars thus provide a unique window into the earliest Galactic stages, but are exceedingly rare and difficult to find. Here we present the discovery of an ultra-metal-poor star, Pristine_221.8781+9.7844, using narrow-band Ca H&K photometry from the Pristine survey. Follow-up medium and high-resolution spectroscopy confirms the ultra-metal-poor nature of Pristine_221.8781+9.7844 ([Fe/H] = -4.66 +/- 0.13 in 1D LTE) with an enhancement of 0.3-0.4 dex in alpha-elements relative to Fe, and an unusually low carbon abundance. We derive an upper limit of A(C) = 5.6, well below typical A(C) values for such ultra metal-poor stars. This makes Pristine_221.8781+9.7844 one of the most metal-poor stars; in fact, it is very similar to the most metal-poor star known (SDSS J102915+172927). The existence of a class of ultra metal-poor stars with low(er) carbon abundances suggest that there must have been several formation channels in the early Universe through which long-lived, low-mass stars were formed., Comment: Re-submitted to MNRAS after minor revisions

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project., Comment: 47 pages, 18 figures; accepted for publication in the Astronomical Journal

Stellar dynamical model globular clusters are introduced into reconstituted versions of the dark matter halos of the Via-Lactea II (VL-2) simulation to follow the star cluster tidal mass loss and stellar stream formation. The clusters initially evolve within their local sub-galactic halo, later being accreted into the main halo. Stars are continually removed from the clusters, but those that emerged in the sub-galactic halos are dispersed in a wide stream when accreted into the main halo. Thin tidal streams that survive to the present can begin to form once a cluster is in the main halo. A higher redshift start places the star clusters in denser halos where they are subject to stronger tides leading to higher average mass loss rates. A z=3 start leads to a rich set of star streams with nearly all within 100 kpc having a remnant progenitor star cluster in the stream. In contrast, in a z=8 start, all star clusters that are accreted onto the main halo are completely dissolved. These results are compared to the available data on Milky-Way streams, where the majority of streams do not have clearly associated globular clusters. which, if generally true, suggests that there were at least twice as many massive globular clusters at high redshift., Comment: AAS submitted, revised

We present the Pristine survey, a new narrow-band photometric survey focused on the metallicity-sensitive Ca H & K lines and conducted in the northern hemisphere with the wide-field imager MegaCam on the Canada-France-Hawaii Telescope (CFHT). This paper reviews our overall survey strategy and discusses the data processing and metallicity calibration. Additionally we review the application of these data to the main aims of the survey, which are to gather a large sample of the most metal-poor stars in the Galaxy, to further characterise the faintest Milky Way satellites, and to map the (metal-poor) substructure in the Galactic halo. The current Pristine footprint comprises over 1,000 deg2 in the Galactic halo ranging from b~30 to 78 and covers many known stellar substructures. We demonstrate that, for SDSS stellar objects, we can calibrate the photometry at the 0.02-magnitude level. The comparison with existing spectroscopic metallicities from SDSS/SEGUE and LAMOST shows that, when combined with SDSS broad-band g and i photometry, we can use the CaHK photometry to infer photometric metallicities with an accuracy of ~0.2 dex from [Fe/H]=-0.5 down to the extremely metal-poor regime ([Fe/H]<-3.0). After the removal of various contaminants, we can efficiently select metal-poor stars and build a very complete sample with high purity. The success rate of uncovering [Fe/H]SEGUE<-3.0 stars among [Fe/H]Pristine<-3.0 selected stars is 24% and 85% of the remaining candidates are still very metal poor ([Fe/H]<-2.0). We further demonstrate that Pristine is well suited to identify the very rare and pristine Galactic stars with [Fe/H]<-4.0, which can teach us valuable lessons about the early Universe., Comment: 22 pages, 17 figures, accepted for publication in MNRAS

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

Thin halo star streams originate from the evaporation of globular clusters and therefore provide information about the early epoch globular cluster population. The observed tidal tails from halo globular clusters in the Milky Way are much more shorter than expected from a star cluster orbiting for 10~Gyr. The discrepancy is likely the result of the assumption that the clusters have been orbiting in a non-evolving galactic halo for a Hubble time. As a first step towards more realistic stream histories, a toy model that combines an idealized merger model with a simplified model of the internal collisional relaxation of individual star clusters is developed. On the average, the velocity dispersion increases with distance causing the density of the stream to decline with distance. Consequently, the streams visible in current data will normally be some fraction of the entire stream. Nevertheless, the high surface density segment of the stellar streams created from the evaporation of the more massive globular clusters should all be visible in low obscuration parts of the sky if closer than about 30~kpc. The Pan-STARRS1 halo volume is used to compare the numbers of halo streams and globular clusters., Comment: revised version

We report on the detection in the combined $Gaia$-DR1/RAVE data of a lack of disk stars in the solar neighbourhood with velocities close to zero angular momentum. We propose that this may be caused by the scattering of stars with very low angular momentum onto chaotic, halo-type orbits when they pass through the Galactic nucleus. We model the effect in a Milky-Way like potential and fit the resulting model directly to the data, finding a likelihood ($\sim2.7\sigma$) of a dip in the distribution. Using this effect, we can make a dynamical measurement of the Solar rotation velocity around the Galactic center: $v_{\odot}=239\pm9$ km s$^{-1}$. Combined with the measured proper motion of Sgr A$^*$, this measurement gives a measurement of the distance to the Galactic centre: $R_0=7.9\pm0.3$ kpc., Comment: 5 pages, 4 figures, submitted to ApJL

Tidal debris streams from galaxy satellites can provide insight into the dark matter distribution in halos. This is because we have more information about stars in a debris structure than about a purely random population of stars: we know that in the past they were all bound to the same dwarf galaxy; and we know that they form a dynamically cold population moving on similar orbits. They also probe a different region of the matter distribution in a galaxy than many other methods of mass determination, as their orbits take them far beyond the typical extent of those for the bulk of stars. Although conclusive results from this information have yet to be obtained, significant progress has been made in developing the methodologies for determining both the global mass distribution of the Milky Way's dark matter halo and the amount of dark matter substructure within it. Methods for measuring the halo shape are divided into "predictive methods," which predict the tidal debris properties from the progenitor satellite's mass and orbit, given an assumed parent galaxy mass distribution; and "fundamental methods," which exploit properties fundamental to the nature of tidal debris as global potential constraints. Methods for quantifying the prevalence of dark matter subhalos within halos through the analysis of the gaps left in tidal streams after these substructures pass through them are reviewed., Comment: Chapter~7 from the volume "Tidal Streams in the Local Group and Beyond: Observations and Implications"; ed. Newberg, H.~J., \& Carlin, J.~L.\ 2016, Spr inger International Publishing, Astrophysics and Space Science Library, 420\\ IS BN 978-3-319-19335-9; DOI 10.1007/978-3-319-19336-6\_1

We investigate whether the recently discovered Phoenix stream may be part of a much longer structure that includes the previously discovered Hermus stream. Using a simple model of the Galaxy with a disk, bulge, and a spherical dark matter halo, we show that a nearly circular orbit, highly inclined with respect to the disk, can be found that fits the positions, orientations, and distances of both streams. While the two streams are somewhat misaligned in the sense that they do not occupy the same plane, nodal precession due to the Milky Way disk potential naturally brings the orbit into line with each stream in the course of half an orbit. We consequently consider a common origin for the two streams as plausible. Based on our best fitting orbit, we make predictions for the positions, distances, radial velocities, and proper motions along each stream. If our hypothesis is borne out by measurements, then at ~183 degrees (~235 degrees with respect to the Galactic center) and ~76 kpc in length, Phoenix-Hermus would become the longest cold stream yet found. This would make it a particularly valuable new probe of the shape and mass of the Galactic halo out to ~20 kpc., Comment: 8 pages, 4 figures, accepted for publication in ApJ Letters

We simulate the tidal disruption of a collisionless N-body globular star cluster in a total of 300 different orbits selected to have galactocentric radii between 10 and 30 kpc in four dark matter halos: (a) a spherical halo with no subhalos, (b) a spherical halo with subhalos, (c) a realistic halo with no subhalos, and (d) a realistic halo with subhalos. This allows us to isolate and study how the halo's (lack of) dynamical symmetry and substructures affect the dispersal of tidal debris. The realistic halos are constructed from the snapshot of the Via Lactea II simulation at redshift zero. We find that the overall halo's lack of dynamical symmetry disperses tidal debris to make the streams fluffier, consistent with previous studies of tidal debris of dwarf galaxies in larger orbits than ours in this study. On the other hand, subhalos in realistic potentials can locally enhance the densities along streams, making streams denser than their counterparts in smooth potentials. We show that many long and thin streams can survive in a realistic and lumpy halo for a Hubble time. This suggests that upcoming stellar surveys will likely uncover more thin streams which may contain density gaps that have been shown to be promising probes for dark matter substructures., Comment: 11 pages, 9 figures, ApJ accepted. Minor clarification in abstract

The GD-1 star stream is currently the best available for identifying density fluctuations, "gaps", along its length as a test of the LCDM prediction of large numbers of dark matter sub-halos orbiting in the halo. Density variations of some form are present, since the variance of the density along the stream is three times that expected from the empirically estimated variation in the filtered mean star counts. The density variations are characterized with filters that approximate the shape of sub-halo gravitationally induced stream gaps. The filters locate gaps and measures their amplitude, leading to a measurement of the distribution of gap widths. To gain understanding of the factors influencing the gap width distribution, a suite of collisionless n-body simulations for a GD-1 like orbit in a Milky Way-like potential provides a dynamically realistic statistical prediction of the gap distribution. The simulations show that every location in the stream has been disturbed to some degree by a sub-halo. The small gaps found via the filtering are largely noise. Larger gaps, those longer than 1 kpc, or 10 degrees for GD-1, are the source of the excess variance. The suite of stream simulations shows that sub-halos at the predicted inner halo abundance or possibly somewhat higher can produce the required large sale density variations., Comment: ApJ accepted

Using a shallow, two-color survey carried out with the Dark Energy Camera, we detect the southern, possibly trailing arm of the Orphan Stream. The stream is reliably detected to a declination of $-38^\circ$, bringing the total known length of the Orphan stream to $108^\circ$. We find a slight offset or "S" shape in the stream at $\delta \simeq -14^\circ$ that would be consistent with the transition from leading to trailing arms. This coincides with a moderate concentration of $137 \pm 25$ stars (to $g = 21.6$) that we consider a possible remnant of the Orphan progenitor. The position of this feature is in agreement with previous predictions., Comment: 4 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

The velocity, position, and action variable evolution of a tidal stream drawn out of a star cluster in a triaxial isochrone potential containing a sub-halo population reproduces many of the orbital effects of more general cosmological halos but allows easy calculation of orbital actions. We employ a spherical shell code which we show accurately reproduces the results of a tree gravity code for a collisionless star cluster. Streams from clusters on high eccentricity orbits, $e\gtrsim 0.6$, can spread out so much that the amount of material at high enough surface density to stand out on the sky may be only a few percent of the stream's total mass. Low eccentricity streams remain more spatially coherent, but sub-halos both broaden the stream and displace the centerline with details depending on the orbits allowed within the potential. Overall, the majority of stream particles have changes in their total actions of only 1-2\%, leaving the mean stream relatively undisturbed. A halo with 1\% of the mass in sub-halos typically spreads the velocity distribution about a factor of two wider than would be expected for a smooth halo. Strong density variations, "gaps", along with mean velocity offsets, are clearly detected in low eccentricity streams for even a 0.2\% sub-halo mass fraction. Around one hundred velocity measurements per kiloparsec of stream will enable tests for the presence of a local sub-halo density as small as 0.2-0.5\% of the local mass density, with about 1\% predicted for 30 kiloparsec orbital radii streams., Comment: ApJ submitted after revisions in response to referee

A sufficiently extended satellite in the tidal field of a host galaxy loses mass to create nearly symmetric leading and trailing tidal streams. We study the case in which tidal heating drives mass loss from a low mass satellite. The stream effectively has two dynamical components, a common angular momentum core superposed with episodic pulses with a broader angular momentum distribution. The pulses appear as spurs on the stream, oscillating above and below the stream centerline, stretching and blurring in configuration space as they move away from the cluster. Low orbital eccentricity streams are smoother and have less differential motion than high eccentricity streams. The tail of a high eccentricity stream can develop a fan of particles which wraps around at apocenter in a shell feature. We show that scaling the essentially stationary action-angle variables with the cube root of the satellite mass allows a low mass satellite stream to accurately predict the features in the stream from a satellite a thousand times more massive. As a practical astrophysical application, we demonstrate that narrow gaps in a moderate eccentricity stream, such as GD-1, blur out to 50\% contrast over approximately 6 radial periods. A high eccentricity stream, such as Pal~5, will blur small gaps in only radial 2 orbits as can be understood from the much larger dispersion of angular momentum in the stream., Comment: ApJ accepted