Showing 1–50 of 237 results for author: Carlberg, R

GD-1 Stellar Stream and Cocoon in the DESI Early Data Release

Authors: Monica Valluri, Parker Fagrelius, Sergey. E. Koposov, Ting S. Li, Oleg Y. Gnedin, Eric F. Bell, Raymond G. Carlberg, Andrew P. Cooper, Jessia N. Aguilar, Carlos Allende Prieto, Vasily Belokurov, Leandro Beraldo e Silva, David Brooks, Amanda Byström, Todd Claybaugh, Kyle Dawson, Arjun Dey, Peter Doel, Jaime E. Forero-Romero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Klaus Honscheid, T . Kisner, Anthony Kremin, A. Lambert , et al. (27 additional authors not shown)

Abstract: 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 compo… ▽ More 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. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ, 23 pages, 13 figures 4 tables

Star Stream Velocity Distributions in CDM and WDM Galactic Halos

Authors: Raymond G. Carlberg, Adrian Jenkins, Carlos S. Frenk, Andrew P. Cooper

Abstract: 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,… ▽ More 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. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: AAS Journals submitted

Euclid preparation. Detecting globular clusters in the Euclid survey

Authors: Euclid Collaboration, K. Voggel, A. Lançon, T. Saifollahi, S. S. Larsen, M. Cantiello, M. Rejkuba, J. -C. Cuillandre, P. Hudelot, A. A. Nucita, M. Urbano, E. Romelli, M. A. Raj, M. Schirmer, C. Tortora, Abdurro'uf, F. Annibali, M. Baes, P. Boldrini, R. Cabanac, D. Carollo, C. J. Conselice, P. -A. Duc, A. M. N. Ferguson, L. K. Hunt , et al. (247 additional authors not shown)

Abstract: Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid… ▽ More Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics. △ Less

Submitted 29 May, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A

Euclid: Early Release Observations -- Unveiling the morphology of two Milky Way globular clusters out to their periphery

Authors: D. Massari, E. Dalessandro, D. Erkal, E. Balbinot, J. Bovy, I. McDonald, A. M. N. Ferguson, S. S. Larsen, A. Lançon, F. Annibali, B. Goldman, P. B. Kuzma, K. Voggel, T. Saifollahi, J. -C. Cuillandre, M. Schirmer, M. Kluge, B. Altieri, A. Amara, S. Andreon, N. Auricchio, M. Baldi, A. Balestra, S. Bardelli, A. Basset , et al. (136 additional authors not shown)

Abstract: As part of the Euclid Early Release Observations (ERO) programme, we analyse deep, wide-field imaging from the VIS and NISP instruments of two Milky Way globular clusters (GCs), namely NGC 6254 (M10) and NGC 6397, to look for observational evidence of their dynamical interaction with the Milky Way. We search for such an interaction in the form of structural and morphological features in the cluste… ▽ More As part of the Euclid Early Release Observations (ERO) programme, we analyse deep, wide-field imaging from the VIS and NISP instruments of two Milky Way globular clusters (GCs), namely NGC 6254 (M10) and NGC 6397, to look for observational evidence of their dynamical interaction with the Milky Way. We search for such an interaction in the form of structural and morphological features in the clusters' outermost regions, which are suggestive of the development of tidal tails on scales larger than those sampled by the ERO programme. Our multi-band photometric analysis results in deep and well-behaved colour-magnitude diagrams that, in turn, enable an accurate membership selection. The surface brightness profiles built from these samples of member stars are the deepest ever obtained for these two Milky Way GCs, reaching down to $\sim30.0$ mag~arcsec$^{-2}$, which is about $1.5$ mag arcsec$^{-2}$ below the current limit. The investigation of the two-dimensional density map of NGC 6254 reveals an elongated morphology of the cluster peripheries in the direction and with the amplitude predicted by $N$-body simulations of the cluster's dynamical evolution, at high statistical significance. We interpret this as strong evidence for the first detection of tidally induced morphological distortion around this cluster. The density map of NGC 6397 reveals a slightly elliptical morphology, in agreement with previous studies, which requires further investigation on larger scales to be properly interpreted. This ERO project thus demonstrates the power of Euclid in studying the outer regions of GCs at an unprecedented level of detail, thanks to the combination of large field of view, high spatial resolution, and depth enabled by the telescope. Our results highlight the future Euclid survey as the ideal data set to investigate GC tidal tails and stellar streams. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 15 pages, 18 figures. Paper accepted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

The Pristine Inner Galaxy Survey (PIGS) IX. The largest detailed chemical analysis of very metal-poor stars in the Sagittarius dwarf galaxy

Authors: Federico Sestito, Sara Vitali, Paula Jofre, Kim A. Venn, David S. Aguado, Claudia Aguilera-Gómez, Anke Ardern-Arentsen, Danielle de Brito Silva, Raymond Carlberg, Camilla J. L. Eldridge, Felipe Gran, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Nicolas F. Martin, Tadafumi Matsuno, Samuel Rusterucci, Else Starkenburg, Akshara Viswanathan

Abstract: The most metal-poor stars provide valuable insights into the early chemical enrichment history of a system, carrying the chemical imprints of the first generations of supernovae. The most metal-poor region of the Sagittarius dwarf galaxy remains inadequately observed and characterised. To date, only $\sim4$ stars with [Fe/H]~$<-2.0$ have been chemically analysed with high-resolution spectroscopy.… ▽ More The most metal-poor stars provide valuable insights into the early chemical enrichment history of a system, carrying the chemical imprints of the first generations of supernovae. The most metal-poor region of the Sagittarius dwarf galaxy remains inadequately observed and characterised. To date, only $\sim4$ stars with [Fe/H]~$<-2.0$ have been chemically analysed with high-resolution spectroscopy. In this study, we present the most extensive chemical abundance analysis of 12 low-metallicity stars with metallicities down to [Fe/H]~$=-3.26$ and located in the main body of Sagittarius. These targets, selected from the Pristine Inner Galaxy Survey, were observed using the MIKE high-resolution spectrograph at the {\it Magellan-Clay} telescope, which allowed us to measure up to 17 chemical species. The chemical composition of these stars reflects the imprint of a variety of type~II supernovae (SNe~II). A combination of low- to intermediate-mass high-energy SNe and hypernovae ($\sim10-70\msun$) is required to account for the abundance patterns of the lighter elements up to the Fe-peak. The trend of the heavy elements suggests the involvement of compact binary merger events and fast-rotating (up to $\sim300\kms$) intermediate-mass to massive metal-poor stars ($\sim25-120\msun$) that are the sources of rapid and slow processes, respectively. Additionally, asymptotic giant branch stars contribute to a wide dispersion of [Ba/Mg] and [Ba/Eu]. The absence of an $α-$knee in our data indicates that type Ia supernovae did not contribute in the very metal-poor region ([Fe/H]~$\leq-2.0$). However, they might have started to pollute the interstellar medium at [Fe/H]~$>-2.0$, given the relatively low [Co/Fe] in this metallicity region. △ Less

Submitted 2 July, 2024; v1 submitted 30 April, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A. New plot on [Eu/Mg]. Some refs updated

The Pristine Inner Galaxy Survey (PIGS) VIII: Characterising the orbital properties of the ancient, very metal-poor inner Milky Way

Authors: Anke Ardern-Arentsen, Giacomo Monari, Anna B. A. Queiroz, Else Starkenburg, Nicolas F. Martin, Cristina Chiappini, David S. Aguado, Vasily Belokurov, Ray Carlberg, Stephanie Monty, GyuChul Myeong, Mathias Schultheis, Federico Sestito, Kim A. Venn, Sara Vitali, Zhen Yuan, Hanyuan Zhang, Sven Buder, Geraint F. Lewis, William H. Oliver, Zhen Wan, Daniel B. Zucker

Abstract: The oldest stars in the Milky Way (born in the first few billion years) are expected to have a high density in the inner few kpc, spatially overlapping with the Galactic bulge. We use spectroscopic data from the Pristine Inner Galaxy Survey (PIGS) to study the dynamical properties of ancient, metal-poor inner Galaxy stars. We compute distances using StarHorse, and orbital properties in a barred Ga… ▽ More The oldest stars in the Milky Way (born in the first few billion years) are expected to have a high density in the inner few kpc, spatially overlapping with the Galactic bulge. We use spectroscopic data from the Pristine Inner Galaxy Survey (PIGS) to study the dynamical properties of ancient, metal-poor inner Galaxy stars. We compute distances using StarHorse, and orbital properties in a barred Galactic potential. With this paper, we release the spectroscopic AAT/PIGS catalogue (13 235 stars). We find that most PIGS stars have orbits typical for a pressure-supported population. The fraction of stars confined to the inner Galaxy decreases with decreasing metallicity, but many very metal-poor stars (VMP, [Fe/H] < -2.0) stay confined (~ 60% stay within 5 kpc). The azimuthal velocity v$_φ$ also decreases between [Fe/H] = -1.0 and -2.0, but is constant for VMP stars (at ~ 40 km/s). The carbon-enhanced metal-poor (CEMP) stars in PIGS appear to have similar orbital properties compared to normal VMP stars. Our results suggest a possible transition between two spheroidal components - a more metal-rich, more concentrated, faster rotating component, and a more metal-poor, more extended and slower/non-rotating component. We propose that the former may be connected to pre-disc in-situ stars (or those born in large building blocks), whereas the latter may be dominated by contributions from smaller galaxies. This is an exciting era where large metal-poor samples, such as in this work (as well as upcoming surveys, e.g., 4MOST), shed light on the earliest evolution of our Galaxy. △ Less

Submitted 16 April, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: Accepted to MNRAS, new: Figures 4, 8 and 9 and Table 2 - Figure 9 shows two inner halo components

Gaia's brightest very metal-poor (VMP) stars: A metallicity catalogue of a thousand VMP stars from Gaia RVS spectra

Authors: Akshara Viswanathan, Else Starkenburg, Tadafumi Matsuno, Kim A. Venn, Nicolas F. Martin, Nicolas Longeard, Anke Ardern-Arentsen, Raymond G. Carlberg, Sebastien Fabbro, Georges Kordopatis, Martin Montelius, Federico Sestito, Zhen Yuan

Abstract: 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 ar… ▽ More 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. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 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

Journal ref: A&A 683, L11 (2024)

The Pristine survey -- XXIII. Data Release 1 and an all-sky metallicity catalogue based on Gaia DR3 BP/RP spectro-photometry

Authors: Nicolas F. Martin, Else Starkenburg, Zhen Yuan, Morgan Fouesneau, Anke Arentsen, Francesca De Angeli, Felipe Gran, Martin Montelius, René Andrae, Michele Bellazzini, Paolo Montegriffo, Anna F. Esselink, Hanyuan Zhang, Kim A. Venn, Akshara Viswanathan, David S. Aguado, Giuseppina Battaglia, Manuel Bayer, Piercarlo Bonifacio, Elisabetta Caffau, Patrick Côté, Raymond Carlberg, Sébastien Fabbro, Emma Fernández Alvar, Jonay I. González Hernández , et al. (15 additional authors not shown)

Abstract: We use the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way stars that has been mapping more than 6,500 deg^2 of the northern sky with the CFHT since 2015. These synthetic magnitudes are used for an ab… ▽ More We use the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way stars that has been mapping more than 6,500 deg^2 of the northern sky with the CFHT since 2015. These synthetic magnitudes are used for an absolute re-calibration of the deeper Pristine photometry and, combined with broadband Gaia information, synthetic and Pristine CaHK magnitudes are used to estimate photometric metallicities over the whole sky. The resulting metallicity catalogue is accurate down to [Fe/H]~-3.5 and is particularly suited for the exploration of the metal-poor Milky Way ([Fe/H]<-1.0). We make available here the catalogue of synthetic CaHK_syn magnitudes for all stars with BP/RP information in Gaia DR3, as well as an associated catalogue of more than ~30 million photometric metallicities for high S/N FGK stars. This paper further provides the first public DR of the Pristine catalogue in the form of higher quality recalibrated Pristine CaHK magnitudes and photometric metallicities for all stars in common with the BP/RP information in Gaia DR3. We demonstrate that, when available, the much deeper Pristine data greatly enhances the quality of the derived metallicities, in particular at the faint end of the catalogue (G_BP>16). Combined, both catalogues include more than 2 million metal-poor star candidates as well as more than 200,000 and ~8,000 very and extremely metal-poor candidates. Finally, we show that these metallicity catalogues can be used efficiently, among other applications, for Galactic archaeology, to hunt for the most metal-poor stars, and to study how the structure of the Milky Way varies with metallicity, from the flat distribution of disk stars to the spheroid-shaped metal-poor halo. (Shortened) △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 30 pages, 24 figures, submitted to A&A. First two authors are co-first author. The CaHK photometry catalogue and the two photometric metallicity catalogues are available, before acceptance, as large compressed csv files at: https://seafile.unistra.fr/d/ee0c0f05719d4368bcbb/

The stability of some galaxy disks is still perplexing

Authors: J A Sellwood, R G Carlberg

Abstract: The problem of how some disk galaxies avoid forming bars remains unsolved. Many galaxy models having reasonable properties continue to manifest vigorous instabilities that rapidly form strong bars and no widely-accepted idea has yet been advanced to account for how some disk galaxies manage to avoid this instability. It is encouraging that not all galaxies formed in recent cosmological simulations… ▽ More The problem of how some disk galaxies avoid forming bars remains unsolved. Many galaxy models having reasonable properties continue to manifest vigorous instabilities that rapidly form strong bars and no widely-accepted idea has yet been advanced to account for how some disk galaxies manage to avoid this instability. It is encouraging that not all galaxies formed in recent cosmological simulations possess bars, but the dynamical explanation for this result is unclear. The unstable mode that creates a bar is understood as a standing wave in a cavity that reflects off the disk center and the corotation radius, with amplification at corotation. Here we use simulations to address one further idea that may perhaps inhibit the feedback loop and therefore contribute to stability, which is to make the disk center dynamically hot and/or to taper away mass from the inner disk, which could be masked by a bulge. Unfortunately, we find that neither strategy makes much difference to the global stability of the disk in the models we have tried. While deep density cutouts do indeed prevent feedback through the center, they still reflect incoming waves and thereby provoke a slightly different instability that again leads to a strong bar. △ Less

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

Comments: Added reference to Smercina et al (2023). Appeared in ApJ 958:182

Sub-Halo Spreading of Thin Tidal Star Streams

Authors: Raymond G. Carlberg, Hayley Agler

Abstract: 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 st… ▽ More 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. △ Less

Submitted 3 July, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: accepted for publication in the ApJ

Spiral instabilities: Mode saturation and decay

Authors: J A Sellwood, R G Carlberg

Abstract: This paper continues a series reporting different aspects of the behaviour of disc galaxy simulations that support spiral instabilities. The focus in this paper is to demonstrate how linear spiral instabilities saturate and decay, and how the properties of the disc affect the limiting amplitude of the spirals. Once again, we employ idealized models that each possess a single instability that we fo… ▽ More This paper continues a series reporting different aspects of the behaviour of disc galaxy simulations that support spiral instabilities. The focus in this paper is to demonstrate how linear spiral instabilities saturate and decay, and how the properties of the disc affect the limiting amplitude of the spirals. Once again, we employ idealized models that each possess a single instability that we follow until it has run its course. Remarkably, we find a tight correlation between the growth rate of the mode and its limiting amplitude, albeit from only six simulations. We show that non-linear orbit deflections near corotation cause the mode to saturate, and that the more time available in a slowly-growing mode creates the critical deflections at lower amplitude. We also find that scattering at the inner Lindblad resonance is insignificant until after the mode has saturated. Our objective in this series of papers, which we believe we have now achieved, has been to develop a convincing and well-documented account of the physical behaviour of the spiral patterns that have been observed in simulations by others, and by ourselves, for many decades. Understanding the simulations is an important step towards the greater objective, which is to find observational evidence from galaxies that could confront the identified mechanism. △ Less

Submitted 1 October, 2022; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 13 pages, 13 figures, to appear in MNRAS, corrected Fig 13 only

The Pristine Inner Galaxy Survey (PIGS) V: a chemo-dynamical investigation of the early assembly of the Milky Way with the most metal-poor stars in the bulge

Authors: Federico Sestito, Kim A. Venn, Anke Arentsen, David Aguado, Collin L. Kielty, Carmela Lardo, Nicolas F. Martin, Julio F. Navarro, Else Starkenburg, Fletcher Waller, Raymond G. Carlberg, Patrick François, Jonay I. González Hernández, Georges Kordopatis, Sara Vitali, Zhen Yuan

Abstract: 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… ▽ More 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. △ Less

Submitted 11 November, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 22 pages + appendix, 14 Figures + 3 in the appendix, Accepted for publication in MNRAS

The origin and evolution of the normal Type Ia SN 2018aoz with infant-phase reddening and excess emission

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago GonzÁlez-GaitÁn, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, LluÍs Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg, S. Bradley Cenko , et al. (18 additional authors not shown)

Abstract: 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 int… ▽ More 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$α$ 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. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

The Pristine Inner Galaxy Survey (PIGS) IV: A photometric metallicity analysis of the Sagittarius dwarf spheroidal galaxy

Authors: Sara Vitali, Anke Arentsen, Else Starkenburg, Paula Jofré, Nicolas F. Martin, David S. Aguado, Raymond Carlberg, Jonay I. González Hernández, Rodrigo Ibata, Georges Kordopatis, Khyati Malhan, Pau Ramos, Federico Sestito, Zhen Yuan, Sven Buder, Geraint F. Lewis, Zhen Wan, Daniel B. Zucker

Abstract: We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of… ▽ More We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of spectroscopic metallicities allows us to cover an unprecedented large area ($\sim 100$ square degrees) of the dwarf galaxy, and to study the spatial distribution of its members as function of metallicity with little selection effects. Our study compares the spatial distributions of a metal-poor population of 9719 stars with [Fe/H] $< -1.3$ and a metal rich one of 30115 stars with [Fe/H] $> -1.0$. The photometric Sgr sample also allows us to assemble the largest sample of 1150 very metal-poor Sgr candidates ([Fe/H] $< -2.0$). By investigating and fitting the spatial properties of the metal-rich and metal-poor population, we find a negative metallicity gradient which extends up to 12 degrees from the Sgr center (or $\sim 5.5$ kpc at the distance of Sgr), the limit of our footprint. We conclude that the relative number of metal-poor stars increases in the outer areas of the galaxy, while the central region is dominated by metal-rich stars. These finding suggest an outside-in formation process and are an indication of the extended formation history of Sgr, which has been affected by the tidal interaction between Sgr and the Milky Way. △ Less

Submitted 4 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: The paper was accepted on 20 September by MNRAS

Constraints on galaxy formation from the cosmic-infrared-background$\,-\,$optical-imaging cross-correlation using $\boldsymbol{Herschel}$ and UNIONS

Authors: Seunghwan Lim, Ryley Hill, Douglas Scott, Ludovic van Waerbeke, Jean-Charles Cuillandre, Raymond G. Carlberg, Nora Elisa Chisari, Andrej Dvornik, Thomas Erben, Stephen Gwyn, Alan W. McConnachie, Marc-Antoine Miville-Deschênes, Angus H. Wright, Pierre-Alain Duc

Abstract: 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 infrared background and cosmic optical background fluctuations. With a combined sky area of $91\,{\rm deg}^2$ and a minimum resolved scale of $18\,$arcsec, the cros… ▽ More 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 infrared background and cosmic optical background fluctuations. With a combined sky area of $91\,{\rm deg}^2$ and a minimum resolved scale of $18\,$arcsec, 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 22\,σ$ ($\gtrsim 18\,σ$ 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 approximately a factor of 10 higher than the extragalactic signal, with the contamination being estimated using a linear regression from several external survey maps; however, for the other fields, the contamination is typically 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. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 35 pages, 27 figures, submitted to MNRAS, comments are welcome

The Dark Matter Halo of M54

Authors: Raymond G. Carlberg, Carl J. Grillmair

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 28 June, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: ApJ accepted

Infant-phase reddening by surface Fe-peak elements in a normal Type Ia Supernova

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago Gonzalez-Gaitan, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, Lluis Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Iair Arcavi, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg , et al. (19 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 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 Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

Authors: Khyati Malhan, Rodrigo A. Ibata, Sanjib Sharma, Benoit Famaey, Michele Bellazzini, Raymond G. Carlberg, Richard D'Souza, Zhen Yuan, Nicolas F. Martin, Guillaume F. Thomas

Abstract: 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, stella… ▽ More 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. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: Published in ApJ. Key results in Table 3. The paper contains 38 pages, 17 figures, 4 Tables

The Pristine survey -- XVI. The metallicity of 26 stellar streams around the Milky Way detected with the STREAMFINDER in Gaia EDR3

Authors: Nicolas F. Martin, Rodrigo A. Ibata, Else Starkenburg, Zhen Yuan, Khyati Malhan, Michele Bellazzini, Akshara Viswanathan, David Aguado, Anke Arentsen, Piercarlo Bonifacio, Ray Carlberg, Jonay I. González Hernández, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Alan W. McConnachie, Julio Navarro, Rubén Sánchez-Janssen, Federico Sestito, Guillaume F. Thomas, Kim A. Venn, Sara Vitali, Karina T. Voggel

Abstract: We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by the application of the STREAMFINDER algorithm to the Gaia EDR3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, suppo… ▽ More We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by the application of the STREAMFINDER algorithm to the Gaia EDR3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, supporting the reliability of STREAMFINDER in finding streams and the power of Pristine to measure precise metallicities. We further present 6 candidate structures with coherent phase-space and metallicity signals that are very likely streams. The majority of studied streams are very metal-poor (14 structures with [Fe/H]<-2.0) and include 3 systems with [Fe/H]<-2.9 (C-11, C-19, and C-20). These streams could be the closest debris of low-luminosity dwarf galaxies or may have originated from globular clusters of significantly lower metallicity than any known current Milky Way globular cluster. Our study shows that the promise of the Gaia data for Galactic Archeology studies can be substantially strengthened by quality photometric metallicities, allowing us to peer back into the earliest epochs of the formation of our Galaxy and its stellar halo constituents. △ Less

Submitted 25 August, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 25 pages, 29 figures, accepted for publication in MNRAS. Updated version has only minor changes compared to initial submission. The metallicity of most streams have marginally changed in the updated version

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

Authors: Nicolas F. Martin, Kim A. Venn, David S. Aguado, Else Starkenburg, Jonay I. González Hernández, Rodrigo A. Ibata, Piercarlo Bonifacio, Elisabetta Caffau, Federico Sestito, Anke Arentsen, Carlos Allende Prieto, Raymond G. Carlberg, Sébastien Fabbro, Morgan Fouesneau, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Khyati Malhan, Lyudmila I. Mashonkina, Alan W. McConnachie, Julio F. Navarro, Rubén Sánchez Janssen, Guillaume F. Thomas, Zhen Yuan , et al. (1 additional authors not shown)

Abstract: 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 s… ▽ More 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. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 14 pages, 6 figures, 5 tables. Nature, accepted

Ram Pressure Candidates in UNIONS

Authors: Ian D. Roberts, Laura C. Parker, Stephen Gwyn, Michael J. Hudson, Raymond Carlberg, Alan McConnachie, Jean-Charles Cuillandre, Kenneth C. Chambers, Pierre-Alain Duc, Hisanori Furusawa, Raphael Gavazzi, Vanessa Hill, Mark E. Huber, Rodrigo Ibata, Martin Kilbinger, Simona Mei, Yannick Mellier, Satoshi Miyazaki, Masamune Oguri, Richard J. Wainscoat

Abstract: We present a search for disturbed, candidate ram pressure stripping galaxies across more than 50 spectroscopically selected SDSS groups and clusters. Forty-eight ram pressure candidates are visually identified in these systems using high quality UNIONS imaging from the Canada-France Hawaii Telescope, covering ~6200 and ~2800 square degrees in the u- and r-bands respectively. Ram pressure candidate… ▽ More We present a search for disturbed, candidate ram pressure stripping galaxies across more than 50 spectroscopically selected SDSS groups and clusters. Forty-eight ram pressure candidates are visually identified in these systems using high quality UNIONS imaging from the Canada-France Hawaii Telescope, covering ~6200 and ~2800 square degrees in the u- and r-bands respectively. Ram pressure candidates are found in groups and clusters spanning a wide range in halo mass and include ~30 ram pressure candidates in the group regime ($M_h < 10^{14}$). The observed frequency of ram pressure candidates shows substantial scatter with group/cluster mass, but on average is larger in clusters ($M_h > 10^{14}\,M_\odot$) than groups ($M_h < 10^{14}\,M_\odot$) by a factor of ~2. We find that ram pressure candidates are most commonly low-mass galaxies and have enhanced star formation rates relative to star-forming field galaxies. The enhancement in star formation is largely independent of galaxy mass and strongest for galaxies in clusters. As a result of the large survey footprint and excellent image quality from UNIONS, we are able to identify disturbed galaxies, potentially affected by ram pressure stripping, across a wide range of host environment. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 12 pages, 9 figures, 1 appendix, accepted for publication in MNRAS

The Pristine survey XIV: chemical analysis of two ultra-metal-poor stars

Authors: C. Lardo, L. Mashonkina, P. Jablonka, P. Bonifacio, E. Caffau, D. S. Aguado, J. I. González Hernández, F. Sestito, C. L. Kielty, K. A. Venn, V. Hill, E. Starkenburg, N. F. Martin, T. Sitnova, A. Arentsen, R. G. Carlberg, J. F. Navarro, G. Kordopatis

Abstract: Elemental abundances of the most metal-poor stars reflect the conditions in the early Galaxy and the properties of the first stars. We present a spectroscopic follow-up of two ultra metal-poor stars ([Fe/H]<-4.0) identified by the survey {\em Pristine}: Pristine 221.8781+9.7844 and Pristine 237.8588+12.5660 (hereafter Pr 221 and Pr 237, respectively). Combining data with earlier observations, we f… ▽ More Elemental abundances of the most metal-poor stars reflect the conditions in the early Galaxy and the properties of the first stars. We present a spectroscopic follow-up of two ultra metal-poor stars ([Fe/H]<-4.0) identified by the survey {\em Pristine}: Pristine 221.8781+9.7844 and Pristine 237.8588+12.5660 (hereafter Pr 221 and Pr 237, respectively). Combining data with earlier observations, we find a radial velocity of -149.25 $\pm$ 0.27 and -3.18 $\pm$ 0.19 km/s for Pr 221 and Pr 237, respectively, with no evidence of variability between 2018 and 2020. From a one-dimensional (1D) local thermodynamic equilibrium (LTE) analysis, we measure [Fe/H]$_{\rm LTE}$=-4.79 $\pm$ 0.14 for Pr 221 and [Fe/H]$_{\rm LTE}$=-4.22 $\pm$ 0.12 for Pr 237, in good agreement with previous studies. Abundances of Li, Na, Mg, Al, Si, Ca, Ti, Fe, and Sr were derived based on the non-LTE (NLTE) line formation calculations. When NLTE effects are included, we measure slightly higher metallicities: [Fe/H]$_{\rm NLTE}$=-4.40 $\pm$ 0.13 and [Fe/H]$_{\rm NLTE}$=-3.93 $\pm$ 0.12, for Pr 221 and Pr 237, respectively. Analysis of the G-band yields [C/Fe]$_{\rm 1D-LTE} \leq$ +2.3 and [C/Fe]$_{\rm 1D-LTE} \leq$ +2.0 for Pr 221 and Pr 237. Both stars belong to the low-carbon band. Upper limits on nitrogen abundances are also derived. Abundances for other elements exhibit good agreement with those of stars with similar parameters. Finally, to get insight into the properties of their progenitors, we compare NLTE abundances to theoretical yields of zero-metallicity supernovae. This suggests that the supernovae progenitors had masses ranging from 10.6 to 14.4 M$_{\odot}$ and low-energy explosions with 0.3-1.2 $\times$ 10$^{51}$ erg. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 17 pages, 11 figures, accepted for publication in MNRAS

Euclid preparation: I. The Euclid Wide Survey

Authors: R. Scaramella, J. Amiaux, Y. Mellier, C. Burigana, C. S. Carvalho, J. -C. Cuillandre, A. Da Silva, A. Derosa, J. Dinis, E. Maiorano, M. Maris, I. Tereno, R. Laureijs, T. Boenke, G. Buenadicha, X. Dupac, L. M. Gaspar Venancio, P. Gómez-Álvarez, J. Hoar, J. Lorenzo Alvarez, G. D. Racca, G. Saavedra-Criado, J. Schwartz, R. Vavrek, M. Schirmer , et al. (216 additional authors not shown)

Abstract: Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced s… ▽ More Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 43 pages, 51 figures, submitted to A&A

Journal ref: A&A 662, A112 (2022)

The Pristine Dwarf-Galaxy survey -- IV. Probing the outskirts of the dwarf galaxy Boötes I

Authors: Nicolas Longeard, Pascale Jablonka, Anke Arentsen, Guillaume F. Thomas, David S. Aguado, Raymond G. Carlberg, Romain Lucchesi, Khyati Malhan, Nicolas Martin, Alan W. McConnachie, Julio F. Navarro, Rubén Sánchez-Janssen, Federico Sestito, Else Starkenburg, Zhen Yuan

Abstract: 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'… ▽ More 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. △ Less

Submitted 22 June, 2022; v1 submitted 22 July, 2021; originally announced July 2021.

The Pristine survey XIII: Uncovering the very metal-poor tail of the thin disc

Authors: Emma Fernández-Alvar, Georges Kordopatis, Vanessa Hill, Else Starkenburg, Akshara Viswanathan, Nicolas F. Martin, Guillaume F. Thomas, Julio F. Navarro, Khyati Malhan, Federico Sestito, Jonay I. González Hernández, Raymond G. Carlberg

Abstract: 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_φ$) is similar to the tangential velocity in the galactic longitude direction ($V_{\ell}$). This… ▽ More 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_φ$) is similar to the tangential velocity in the galactic longitude direction ($V_{\ell}$). This allows us to estimate $V_φ$ 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. △ Less

Submitted 1 October, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 18 pages, 14 figures. Accepted for publication in MNRAS

Testing for Dark Matter in the Outskirts of Globular Clusters

Authors: Raymond G. Carlberg, Carl J. Grillmair

Abstract: 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… ▽ More 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. △ Less

Submitted 7 September, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: AAS revised version submitted

Simulating Globular Clusters in Dark Matter Sub-Halos

Authors: Raymond G. Carlberg, Laura C. Keating

Abstract: 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 clust… ▽ More 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. △ Less

Submitted 29 August, 2021; v1 submitted 28 May, 2021; originally announced May 2021.

Comments: ApJ revised submission

The Pristine Inner Galaxy Survey (PIGS) III: carbon-enhanced metal-poor stars in the bulge

Authors: Anke Arentsen, Else Starkenburg, David S. Aguado, Nicolas F. Martin, Vinicius M. Placco, Raymond Carlberg, Jonay I. González Hernández, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Lyudmila I. Mashonkina, Julio F. Navarro, Kim A. Venn, Sven Buder, Geraint F. Lewis, Zhen Wan, Daniel B. Zucker

Abstract: The most metal-deficient stars hold important clues about the early build-up and chemical evolution of the Milky Way, and carbon-enhanced metal-poor (CEMP) stars are of special interest. However, little is known about CEMP stars in the Galactic bulge. In this paper, we use the large spectroscopic sample of metal-poor stars from the Pristine Inner Galaxy Survey (PIGS) to identify CEMP stars ([C/Fe]… ▽ More The most metal-deficient stars hold important clues about the early build-up and chemical evolution of the Milky Way, and carbon-enhanced metal-poor (CEMP) stars are of special interest. However, little is known about CEMP stars in the Galactic bulge. In this paper, we use the large spectroscopic sample of metal-poor stars from the Pristine Inner Galaxy Survey (PIGS) to identify CEMP stars ([C/Fe] > +0.7) in the bulge region and to derive a CEMP fraction. We identify 96 new CEMP stars in the inner Galaxy, of which 62 are very metal-poor ([Fe/H] < -2.0); this is more than a ten-fold increase compared to the seven previously known bulge CEMP stars. The cumulative fraction of CEMP stars in PIGS is $42^{\,+14\,}_{\,-13} \%$ for stars with [Fe/H] < -3.0, and decreases to $16^{\,+3\,}_{\,-3} \%$ for [Fe/H] < -2.5 and $5.7^{\,+0.6\,}_{\,-0.5} \%$ for [Fe/H] < -2.0. The PIGS inner Galaxy CEMP fraction for [Fe/H] < -3.0 is consistent with the halo fraction found in the literature, but at higher metallicities the PIGS fraction is substantially lower. While this can partly be attributed to a photometric selection bias, such bias is unlikely to fully explain the low CEMP fraction at higher metallicities. Considering the typical carbon excesses and metallicity ranges for halo CEMP-s and CEMP-no stars, our results point to a possible deficiency of both CEMP-s and CEMP-no stars (especially the more metal-rich) in the inner Galaxy. The former is potentially related to a difference in the binary fraction, whereas the latter may be the result of a fast chemical enrichment in the early building blocks of the inner Galaxy. △ Less

Submitted 7 May, 2021; originally announced May 2021.

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

The Likelihood of Undiscovered Globular Clusters in the Outskirts of the Milky Way

Authors: Jeremy J. Webb, Raymond G. Carlberg

Abstract: 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 e… ▽ More 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. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 13 pages, 9 figures, 4 tables, resubmitted to MNRAS for publication

Spiral instabilities: Linear and nonlinear effects

Authors: J. A. Sellwood, R. G. Carlberg

Abstract: We present a study of the spiral responses in a stable disc galaxy model to co-orbiting perturbing masses that are evenly spaced around rings. The amplitudes of the responses, or wakes, are proportional to the masses of the perturbations, and we find that the response to a low-mass ring disperses when it is removed -- behaviour that is predicted by linear theory. Higher mass rings cause nonlinear… ▽ More We present a study of the spiral responses in a stable disc galaxy model to co-orbiting perturbing masses that are evenly spaced around rings. The amplitudes of the responses, or wakes, are proportional to the masses of the perturbations, and we find that the response to a low-mass ring disperses when it is removed -- behaviour that is predicted by linear theory. Higher mass rings cause nonlinear changes through scattering at the major resonances, provoking instabilities that were absent before the scattering took place. The separate wake patterns from two rings orbiting at differing frequencies, produce a net response that is an apparently shearing spiral. When the rings have low mass, the evolution of the simulation is both qualitatively and quantitatively reproduced by linear superposition of the two separate responses. We argue that apparently shearing transient spirals in simulations result from the superposition of two or more steadily rotating patterns, each of which is best accounted for as a normal mode of the non-smooth disc. △ Less

Submitted 11 August, 2021; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: 14 pages, 19 figures and one animation. The animation may be viewed at http://www.physics.rutgers.edu/~sellwood/animated-fig.html

Journal ref: MNRAS, 500:5043 (2021)

The hidden past of M92: Detection and characterization of a newly formed 17° long stellar stream using the Canada-France Imaging Survey

Authors: Guillaume F. Thomas, Jaclyn Jensen, Alan McConnachie, Patrick Côté, Kim Venn, Nicolas Longeard, Raymond Carlberg, Scott Chapman, Jean-Charles Cuillandre, Benoit Famaey, Laura Ferrarese, Stephen Gwyn, François Hammer, Rodrigo A. Ibata, Khyati Malhan, Nicolas F. Martin, Simona Mei, Julio F. Navarro, Céline Reylé, Else Starkenburg

Abstract: We present an analysis of the structure, kinematics and orbit of a newly found stellar stream emanating from the globular cluster M92 (NGC 6341). This stream was discovered in an improved matched-filter map of the outer Galaxy, based on a "color-color-magnitude" diagram, created using photometry from the Canada-France Imaging Survey (CFIS) and the Pan-STARRS 1 3$π$ survey (PS1). We find the stream… ▽ More We present an analysis of the structure, kinematics and orbit of a newly found stellar stream emanating from the globular cluster M92 (NGC 6341). This stream was discovered in an improved matched-filter map of the outer Galaxy, based on a "color-color-magnitude" diagram, created using photometry from the Canada-France Imaging Survey (CFIS) and the Pan-STARRS 1 3$π$ survey (PS1). We find the stream to have a length of 17° (2.5 kpc at the distance of M92), a width dispersion of $0.29$°(42 pc) and a stellar mass of $[3.17 \pm{0.89}] \times 10^{4} $ M$_\odot$ ($10 \%$ of the stellar mass of the current main body of M92). We examine the kinematics of main sequence, red giant and blue horizontal branch stars belonging to the stream and that have proper motion measurements from the second data release of Gaia. N-body simulations suggest that the stream was likely formed very recently (during the last $\sim 500$ Myr) forcing us to question the orbital origin of this ancient, metal-poor globular cluster. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 14 pages, 8 figures. Accepted for publication in ApJ

The Pristine Inner Galaxy Survey (PIGS) II: Uncovering the most metal-poor populations in the inner Milky Way

Authors: Anke Arentsen, Else Starkenburg, Nicolas F. Martin, David S. Aguado, Daniel B. Zucker, Carlos Allende Prieto, Vanessa Hill, Kim. A. Venn, Raymond G. Carlberg, Jonay I. González Hernández, Lyudmila I. Mashonkina, Julio F. Navarro, Rubén Sánchez-Janssen, Mathias Schultheis, Guillaume F. Thomas, Kris Youakim, Geraint F. Lewis, Jeffrey D. Simpson, Zhen Wan, Roger E. Cohen, Doug Geisler, Julia E. O'Connell

Abstract: 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 Surv… ▽ More 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. △ Less

Submitted 15 June, 2020; originally announced June 2020.

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

Mass loss from massive globular clusters in tidal fields

Authors: Yohai Meiron, Jeremy J. Webb, Jongsuk Hong, Peter Berczik, Rainer Spurzem, Raymond G. Carlberg

Abstract: 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… ▽ More 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. △ Less

Submitted 24 March, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

The Pristine Dwarf-Galaxy survey -- III. Revealing the nature of the Milky Way globular cluster Sagittarius II

Authors: Nicolas Longeard, Nicolas Martin, Rodrigo Ibata, Else Starkenburg, Pascale Jablonka, David S. Aguado, Raymond G. Carlberg, Patrick Côté, Jonay I. González Hernández, Romain Lucchesi, Khyati Malhan, Julio F. Navarro, Rubén Sánchez-Janssen, Guillaume F. Thomas, Kim Venn, Alan W. McConnachie

Abstract: 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… ▽ More 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. △ Less

Submitted 12 May, 2020; originally announced May 2020.

A Deep CFHT Optical Search for a Counterpart to the Possible Neutron Star -- Black Hole Merger GW190814

Authors: Nicholas Vieira, John J. Ruan, Daryl Haggard, Maria R. Drout, Melania C. Nynka, Hope Boyce, Kristine Spekkens, Samar Safi-Harb, Raymond G. Carlberg, Rodrigo Fernández, Anthony L. Piro, Niloufar Afsariardchi, Dae-Sik Moon

Abstract: 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-Hawai… ▽ More 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. △ Less

Submitted 13 May, 2020; v1 submitted 20 March, 2020; originally announced March 2020.

Comments: 23 pages, 13 figures including appendices; accepted to ApJ with minor revisions

The Pristine survey VIII: The metallicity distribution function of the Milky Way halo down to the extremely metal-poor regime

Authors: K. Youakim, E. Starkenburg, N. F. Martin, G. Matijevic, D. S. Aguado, C. Allende Prieto, A. Arentsen, P. Bonifacio, R. G. Carlberg, J. I. González Hernández, V. Hill, G. Kordopatis, C. Lardo, J. F. Navarro, P. Jablonka, R. Sánchez Janssen, F. Sestito, G. F. Thomas

Abstract: The Pristine survey uses narrow-band photometry to derive precise metallicities down to the extremely metal-poor regime ([Fe/H] < -3), and currently consists of over 4 million FGK-type stars over a sky area of $\sim 2~500\, \mathrm{deg}^2$. We focus our analysis on a subsample of $\sim 80~000$ main sequence turnoff stars with heliocentric distances between 6 and 20 kpc, which we take to be a repre… ▽ More The Pristine survey uses narrow-band photometry to derive precise metallicities down to the extremely metal-poor regime ([Fe/H] < -3), and currently consists of over 4 million FGK-type stars over a sky area of $\sim 2~500\, \mathrm{deg}^2$. We focus our analysis on a subsample of $\sim 80~000$ main sequence turnoff stars with heliocentric distances between 6 and 20 kpc, which we take to be a representative sample of the inner halo. The resulting metallicity distribution function (MDF) has a peak at [Fe/H] = -1.6, and a slope of $Δ$(LogN)/$Δ[Fe/H] = 1.0 \pm 0.1$ in the metallicity range of -3.4 < [Fe/H] < -2.5. This agrees well with a simple closed-box chemical enrichment model in this range, but is shallower than previous spectroscopic MDFs presented in the literature, suggesting that there may be a larger proportion of metal-poor stars in the inner halo than previously reported. We identify the Monoceros/TriAnd/ACS/EBS/A13 structure in metallicity space in a low latitude field in the anticenter direction, and also discuss the possibility that the inner halo is dominated by a single, large merger event, but cannot strongly support or refute this idea with the current data. Finally, based on the MDF of field stars, we estimate the number of expected metal-poor globular clusters in the Milky Way halo to be 5.4 for [Fe/H] < -2.5 and 1.5 for [Fe/H] < -3, suggesting that the lack of low metallicity globular clusters in the Milky Way is not due simply to statistical undersampling. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 17 pages, 12 figures, 5 tables, accepted for publication in MNRAS

Dynamical Simulations of the First Globular Clusters

Authors: Raymond G. Carlberg

Abstract: 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 c… ▽ More 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. △ Less

Submitted 4 April, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: ApJ accepted, edited

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

Authors: Federico Sestito, Nicolas F. Martin, Else Starkenburg, Anke Arentsen, Rodrigo A. Ibata, Nicolas Longeard, Collin Kielty, Kristopher Youakim, Kim A. Venn, David S. Aguado, Raymond G. Carlberg, Jonay I. Gonzalez Hernandez, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Khyati Malhan, Julio F. Navarro, Ruben Sanchez-Janssen, Guillame Thomas, Eline Tolstoy, Thomas G. Wilson, Pedro Alonso Palicio, Spencer Bialek, Rafael Garcia-Dias, Romain Lucchesi , et al. (4 additional authors not shown)

Abstract: 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 g… ▽ More 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. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: Submitted to MNRAS letter, 5 pages with 1 figure. Appendices A and B with 3 figures

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

Authors: Nicholas J. Fantin, Patrick Côté, Alan W. McConnachie, Pierre Bergeron, Jean-Charles Cuillandre, Stephen D. J. Gwyn, Rodrigo A. Ibata, Guillaume F. Thomas, Raymond G. Carlberg, Sébastien Fabbro, Misha Haywood, Ariane Lançon, Geraint F. Lewis, Khyati Malhan, Nicolas F. Martin, Julio F. Navarro, Douglas Scott, Else Starkenburg

Abstract: 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… ▽ More 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. △ Less

Submitted 6 November, 2019; originally announced November 2019.

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

Authors: Khyati Malhan, Rodrigo A. Ibata, Raymond G. Carlberg, Michele Bellazzini, Benoit Famaey, Nicolas F. Martin

Abstract: 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 t… ▽ More 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. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Comments: 6 pages, 3 figures, 1 Table, accepted for publication in ApJ

LRP2020: Astrostatistics in Canada

Authors: Gwendolyn Eadie, Arash Bahramian, Pauline Barmby, Radu Craiu, Derek Bingham, Renée Hložek, JJ Kavelaars, David Stenning, Samantha Benincasa, Guillaume Thomas, Karun Thanjavur, Jo Bovy, Jan Cami, Ray Carlberg, Sam Lawler, Adrian Liu, Henry Ngo, Mubdi Rahman, Michael Rupen

Abstract: (Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summe… ▽ More (Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada. △ Less

Submitted 19 October, 2019; originally announced October 2019.

Comments: White paper E017 submitted to the Canadian Long Range Plan LRP2020

The Pristine survey -- VII. A cleaner view of the Galactic outer halo using blue horizontal branch stars

Authors: Else Starkenburg, Kris Youakim, Nicolas Martin, Guillaume Thomas, David S. Aguado, Anke Arentsen, Raymond G. Carlberg, Jonay I. Gonzalez Hernandez, Rodrigo Ibata, Nicolas Longeard, Alan W. McConnachie, Julio Navarro, Ruben Sanchez-Janssen, Kim A. Venn

Abstract: 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… ▽ More 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. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 14 pages, 11 figures. Accepted for publication in MNRAS

Pristine IX: CFHT ESPaDOnS Spectroscopic Analysis of 115 Bright Metal-Poor Candidate Stars

Authors: Kim Venn, Collin Kielty, Federico Sestito, Else Starkenburg, Nicolas Martin, David S. Aguado, Anke Arentsen, Piercarlo Bonifacio, Elisabetta Caffau, Vanessa Hill, Pascale Jablonka, Carmela Lardo, Lyudmilla Mashonkina, Julio F. Navarro, Chris Sneden, Guillaume Thomas, Kris Youakim, Jonay I. Gonzalez-Hernandez, Ruben Sanchez Janssen, Ray Carlberg, Khyati Malhan

Abstract: A chemo-dynamical analysis of 115 metal-poor candidate stars selected from the narrow-band Pristine photometric survey is presented based on CFHT high-resolution ESPaDoNS spectroscopy. We have discover 28 new bright (V < 15) stars with [Fe/H]<-2.5 and 5 with [Fe/H]<-3.0 for success rates of 40% (28/70) and 19% (5/27), respectively. A detailed model atmospheres analysis is carried out for the 28 ne… ▽ More A chemo-dynamical analysis of 115 metal-poor candidate stars selected from the narrow-band Pristine photometric survey is presented based on CFHT high-resolution ESPaDoNS spectroscopy. We have discover 28 new bright (V < 15) stars with [Fe/H]<-2.5 and 5 with [Fe/H]<-3.0 for success rates of 40% (28/70) and 19% (5/27), respectively. A detailed model atmospheres analysis is carried out for the 28 new metal-poor stars. Stellar parameters were determined from SDSS photometric colours, Gaia DR2 parallaxes, MESA/MIST stellar isochrones, and the initial Pristine survey metallicities, following a Bayesian inference method. Chemical abundances are determined for 10 elements (Na, Mg, Ca, Sc, Ti, Cr, Fe, Ni, Y, Ba). Most stars show chemical abundance patterns that are similar to the normal metal-poor stars in the Galactic halo; however, we also report the discoveries of a new r-process rich star, a new CEMP-s candidate with [Y/Ba]>0, and a metal-poor star with very low [Mg/Fe]. The kinematics and orbits for all of the highly probable metal-poor candidates are determined by combining our precision radial velocities with Gaia DR2 proper motions. Some stars show unusual kinematics for their chemistries, including planar orbits, unbound orbits, and highly elliptical orbits that plunge deeply into the Galactic bulge (Rperi < 0.5 kpc); also, eight stars have orbital energies and actions consistent with the Gaia-Enceladus accretion event. This paper contributes to our understanding of the complex chemo-dynamics of the metal-poor Galaxy, and increases the number of known bright metal-poor stars available for detailed nucleosynthetic studies. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: Submitted to MNRAS, 19 pages (22 figures), and 7 tables

The Pristine Inner Galaxy Survey (PIGS) I: Tracing the kinematics of metal-poor stars in the Galactic bulge

Authors: Anke Arentsen, Else Starkenburg, Nicolas F. Martin, Vanessa Hill, Rodrigo Ibata, Andrea Kunder, Mathias Schultheis, Kim A. Venn, Daniel B. Zucker, David Aguado, Ray Carlberg, Jonay I. González Hernández, Carmela Lardo, Nicolas Longeard, Khyati Malhan, Julio F. Navarro, Ruben Sánchez-Janssen, Federico Sestito, Guillaume Thomas, Kris Youakim, Geraint F. Lewis, Jeffrey D. Simpson, Zhen Wan

Abstract: Our Galaxy is known to contain a central boxy/peanut-shaped bulge, yet the importance of a classical, pressure-supported component within the central part of the Milky Way is still being debated. It should be most visible at low metallicity, a regime that has not yet been studied in detail. Using metallicity-sensitive narrow-band photometry, the Pristine Inner Galaxy Survey (PIGS) has collected a… ▽ More Our Galaxy is known to contain a central boxy/peanut-shaped bulge, yet the importance of a classical, pressure-supported component within the central part of the Milky Way is still being debated. It should be most visible at low metallicity, a regime that has not yet been studied in detail. Using metallicity-sensitive narrow-band photometry, the Pristine Inner Galaxy Survey (PIGS) has collected a large sample of metal-poor ([Fe/H] < -1.0) stars in the inner Galaxy to address this open question. We use PIGS to trace the metal-poor inner Galaxy kinematics as function of metallicity for the first time. We find that the rotational signal decreases with decreasing [Fe/H], until it becomes negligible for the most metal-poor stars. Additionally, the velocity dispersion increases with decreasing metallicity for -3.0 < [Fe/H] < -0.5, with a gradient of -44 $\pm$ 4 km$\,$s$^{-1}\,$dex$^{-1}$. These observations may signal a transition between Galactic components of different metallicities and kinematics, a different mapping onto the boxy/peanut-shaped bulge for former disk stars of different metallicities and/or the secular dynamical and gravitational influence of the bar on the pressure-supported component. Our results provide strong constraints on models that attempt to explain the properties of the inner Galaxy. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 5 pages + appendices, accepted to MNRAS Letters

Dwarfs or giants? Stellar metallicities and distances in the Canada-France-Imaging-Survey from $ugrizG$ multi-band photometry

Authors: Guillaume F. Thomas, Nicholaas Annau, Alan McConnachie, Sebastien Fabbro, Hossen Teimoorinia, Patrick Côté, Jean-Charles Cuillandre, Stephen Gwyn, Rodrigo A. Ibata, Else Starkenburg, Raymond Carlberg, Benoit Famaey, Nicholas Fantin, Laura Ferrarese, Vincent Hénault-Brunnet, Jaclyn Jensen, Ariane Lancon, Geraint F. Lewis, Nicolas F. Martin, Julio F. Navarro, Céline Reylé, Rubén Sánchez-Janssen

Abstract: We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]… ▽ More We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]$<-1.2$, the algorithm succeeds in identifying more than 70% of the giants in the training/test set, with a dwarf contamination fraction below 30% (with respect to the SDSS/SEGUE dataset). The photometric metallicity estimates have uncertainties better than 0.2 dex when compared with the spectroscopic measurements. The distances estimated by the algorithm are valid out to a distance of at least $\sim 80$ kpc without requiring any prior on the stellar distribution, and have fully independent uncertainities that take into account both random and systematic errors. These advances allow us to estimate these stellar parameters for approximately 12 million stars in the photometric dataset. This will enable studies involving the chemical mapping of the distant outer disc and the stellar halo, including their kinematics using the Gaia proper motions. This type of algorithm can be applied in the Southern hemisphere to the first release of LSST data, thus providing an almost complete view of the external components of our Galaxy out to at least $\sim 80$ kpc. Critical to the success of these efforts will be ensuring well-defined spectroscopic training sets that sample a broad range of stellar parameters with minimal biases. A catalogue containing the training/test set and all relevant parameters within the public footprint of CFIS is available online. △ Less

Submitted 9 October, 2019; v1 submitted 7 October, 2019; originally announced October 2019.

Comments: 27 pages. 21 Figures. Accepted for publication in AAS. Training/testing sample available at: http://www.canfar.net/storage/list/gthomas/dwarfOrGiants.fits

Characteristic radii of the Milky Way Globular Clusters

Authors: Andrés E. Piatti, Jeremy Webb, Raymond Carlberg

Abstract: We report on the extent of the effects of the Milky Way's gravitational field in shaping the structural parameters and internal dynamics of its globular cluster population. We make use of a homogeneous, up-to-date data set with kinematics, structural properties, current and initial masses of 156 globular clusters. In general, cluster radii increase as the Milky Way potential weakens; with the core… ▽ More We report on the extent of the effects of the Milky Way's gravitational field in shaping the structural parameters and internal dynamics of its globular cluster population. We make use of a homogeneous, up-to-date data set with kinematics, structural properties, current and initial masses of 156 globular clusters. In general, cluster radii increase as the Milky Way potential weakens; with the core and Jacobi radii being those which increase at the slowest and fastest rate respectively. We interpret this result as the innermost regions of globular clusters being less sensitive to changes in the tidal forces with the Galactocentric distance. The Milky Way's gravitational field also seems to have differentially accelerated the internal dynamical evolution of individual clusters, with those toward the bulge appearing dynamically older. Finally we find a sub-population consisting of both compact and extended globular clusters (as defined by their rh/rJ ratio) beyond 8 kpc that appear to have lost a large fraction of their initial mass lost via disruption. Moreover, we identify a third group with rh/rJ > 0.4, which have lost an even larger fraction of their initial mass by disruption. In both cases the high fraction of mass lost is likely due to their large orbital eccentricities and inclination angles, which lead to them experiencing more tidal shocks at perigalacticon and during disc crossings. Comparing the structural and orbital parameters of individual clusters allows for constraints to be placed on whether or not their evolution was relaxation or tidally dominated. △ Less

Submitted 4 September, 2019; originally announced September 2019.

Comments: 13 pages, 7 figures. Accepted for publication in MNRAS

Modelling the Effects of Dark Matter Substructure on Globular Cluster Evolution with the Tidal Approximation

Authors: Jeremy J. Webb, Jo Bovy, Raymond G. Carlberg, Mark Gieles

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 17 pages, 10 figures, Accepted for publication in MNRAS

Spiral instabilities: Mechanism for recurrence

Authors: J. A. Sellwood, Ray G. Carlberg

Abstract: We argue that self-excited instabilities are the cause of spiral patterns in simulations of unperturbed stellar discs. In previous papers, we have found that spiral patterns were caused by a few concurrent waves, which we claimed were modes. The superposition of a few steadily rotating waves inevitably causes the appearance of the disc to change continuously, and creates the kind of shearing spira… ▽ More We argue that self-excited instabilities are the cause of spiral patterns in simulations of unperturbed stellar discs. In previous papers, we have found that spiral patterns were caused by a few concurrent waves, which we claimed were modes. The superposition of a few steadily rotating waves inevitably causes the appearance of the disc to change continuously, and creates the kind of shearing spiral patterns that have been widely reported. Although we have found that individual modes last for relatively few rotations, spiral activity persists because fresh instabilities appear, which we suspected were excited by the changes to the disc caused by previous disturbances. Here we confirm our suspicion by demonstrating that scattering at either of the Lindblad resonances seeds a new groove-type instability. With this logical gap closed, our understanding of the behaviour in the simulations is almost complete. We believe that our robust mechanism is a major cause of spiral patterns in the old stellar discs of galaxies, including the Milky Way where we have previously reported evidence for resonance scattering in the recently released Gaia data. △ Less

Submitted 15 July, 2019; v1 submitted 10 June, 2019; originally announced June 2019.

Comments: 18 pages, 15 figures, minor revisions in response to a thoughtful referee report

Butterfly in a Cocoon, Understanding the origin and morphology of Globular Cluster Streams: The case of GD-1

Authors: Khyati Malhan, Rodrigo A. Ibata, Raymond G. Carlberg, Monica Valluri, Katherine Freese

Abstract: 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… ▽ More 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$σ$ 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. △ Less

Submitted 30 June, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: 11 pages, 10 figures, accepted for publication in The Astrophysical Journal

Report number: NORDITA-2019-022, LCTP-19-06