The Pristine Dwarf-Galaxy survey II

The Pristine Dwarf-Galaxy survey – II. In-depth observational study of the faint Milky Way satellite Sagittarius II.

Nicolas Longeard, Nicolas Martin, Else Starkenburg, Rodrigo A. Ibata, Michelle L. M. Collins, Benjamin P. M. Laevens, Dougal Mackey, R. Michael Rich , David S. Aguado, Anke Arentsen, Pascale Jablonka, Jonay I. González Hernández, Julio F. Navarro, Rubén Sánchez-Janssen
Université de Strasbourg, CNRS, Observatoire astronomique de Strasbourg, UMR 7550, F-67000 Strasbourg, France
Max-Planck-Institut für Astronomy, Königstuhl 17, D-69117, Heidelberg, Germany
Leibniz Institute for Astrophysics Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
Department of Physics, University of Surrey, Guildford, GU2 7XH, Surrey, UK
Institute of Astrophysics, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile
Department of Astronomy, Yale University, New Haven, CT 06520, USA
Research School of Astronomy and Astrophysics, Australian National University, Canberra, ACT 2611, Australia
University of California Los Angeles, Department of Physics & Astronomy, Los Angeles, CA, USA
Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
GEPI, Observatoire de Paris, PSL Research University, CNRS, Place Jules Janssen, 92190 Meudon, France
Laboratoire d’astrophysique, École Polytechnique Fédérale de Lausanne (EPFL), Observatoire, 1290 Versoix, Switzerland
Instituto de Astrofisica de Canarias, Via Lactea, 38205 La Laguna, Tenerife, Spain
Universidad de La Laguna, Departamento de Astrofisica, 38206 La Laguna, Tenerife, Spain
Dept. of Physics and Astronomy, University of Victoria, P.O. Box 3055, STN CSC, Victoria BC V8W 3P6, Canada
NRC Herzberg Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada
Royal Observatory Edinburgh, Blackford Hill, Edinburgh, EH9 3HJ, UK
July 20, 2019
Abstract

We present an extensive study of the Sagittarius II (Sgr II) stellar system using MegaCam and photometry, narrow-band, metallicity-sensitive Calcium H&K doublet photometry, augmented with Keck II/DEIMOS multi-object spectroscopy. We are able to derive and refine the Sgr II structural and stellar properties: the colour-magnitude diagram implies Sgr II is old (12.0 0.5) Gyr and metal-poor. The CaHK photometry confirms the metal-poor nature of the satellite ([Fe/H] dex) and suggests that Sgr II hosts more than one single stellar population ( dex). From the deep spectroscopic data, the velocity dispersion of the system is found to be after excluding two potential binary stars. Using the Ca infrared triplet measured from our highest signal-to-noise spectra, we are able to confirm the metallicity and dispersion inferred from the Pristine photometric metallicities: ([Fe/H] dex, dex). Sgr II’s metallicity and absolute magnitude (M mag) place the system on the luminosity-metallicity relation of the Milky Way dwarf galaxies despite its small size. The low, but resolved metallicity and velocity dispersions paint the picture of a slightly dark matter-dominated satellite. Furthermore, using the Gaia Data Release 2, we constrain the orbit of the satellite and find an apocenter of and a pericenter of . The orbit of Sgr II is consistent with the trailing arm of the Sgr stream and indicate that it is possibly a satellite of the Sgr dSph that was tidally stripped from the dwarf’s influence.

keywords:
galaxy: Dwarf – Local Group – object: Sagittarius II, Sagittarius stream

1 Introduction

During the history of the Universe, structures such as galaxies formed hierarchically. Therefore, dwarf galaxies (DGs) are particularly old and metal-poor systems and targets of choice to study the history of the local universe. They are systems spanning a wide range of masses and luminosity. Bright dwarf galaxies such as Sculptor (Shapley, 1938b), Draco (Wilson, 1955), or Sextans (Irwin et al., 1990) have been known for decades (Mateo, 1998), but the extensive search for these systems over the last twenty years revealed fainter and fainter systems (Martin et al. 2006, Belokurov et al. 2007, Zucker et al. 2006b). Still, our knowledge of the Milky Way satellites remains incomplete. The recent discoveries of several of those faint galaxy candidates with (Willman et al., 2005a; Belokurov et al., 2007; Drlica-Wagner et al., 2015; Laevens et al., 2015b; Luque et al., 2016) are promising as they might well bring new perspectives to near-field cosmology (Bullock & Boylan-Kolchin, 2017).

The study of these nearby small-scale structures can allow one to explore various problematics in astrophysics, from the faint-end of the galaxy luminosity function (Koposov et al., 2009) to the validity of cosmological models. Therefore, DGs are important cosmological probes (Pawlowski et al. 2017; Tulin & Yu 2017) as the comparison of their observed properties with the predictions made by our current CDM model leads to some discrepancies that we have to understand in order to constrain and refine our cosmological models. DGs are also thought to be among the most dark matter (DM hereafter) dominated systems in the universe (Wolf et al., 2010) and could be useful for the detection of the elusive DM particle through self-annihilation processes (Bertone, Hooper & Silk 2005; Geringer-Sameth, Koushiappas & Walker 2015). However, using faint dwarf galaxies as cosmological probes can be challenging as their exceptional faintness comes with observational challenges. The overall properties and/or even the very nature of the recently discovered systems can sometimes be puzzling as the distinction between galaxy and globular cluster is difficult to make (Conn et al. 2018, Longeard et al. 2018). Therefore, only the combined efforts of deep photometric surveys, such as the Dark Energy Survey (The Dark Energy Survey Collaboration, 2005, DES), the Panoramic Survey Telescope and Rapid Response System (Chambers et al., 2016, PS1), or the Sloan Digital Sky Survey (York et al., 2000, SDSS), and spectroscopic observations can hope to improve our understanding on the faint-end of the luminosity function and the history of the Milky Way (MW).

In this context, we present here the study of the Milky Way satellite Sgr II, discovered by Laevens et al. (2015a, hereafter L15) in PS1, where it was identified as an old (12.5 Gyr) and metal-poor ([Fe/H] = -2.20 dex) dwarf-galaxy candidate. L15 noticed that Sagittarius II had a peculiar location on the sky: its position and distance were found to be consistent with the predictions of models for the Sagittarius stream (Law & Majewski, 2010). They concluded that this satellite might actually have been a satellite of the bright Sagittarius dwarf galaxy discovered by Ibata, Gilmore & Irwin (1994). However, spectroscopic observations were still needed at the time to dynamically tie the stream and Sagittarius II, as well as confirming the galaxy nature of the satellite. Sgr II was also studied by Mutlu-Pakdil et al. (2018, M18) with Magellan/MegaCam photometry, who confirmed the structural properties inferred from L15. Furthermore, using both blue horizontal branch stars (BHBs) and a CMD-fitting technique, they found an old ( Gyr), metal-poor ( dex) stellar population, with an alpha abundance ratio of /Fe dex, and a distance modulus mag. Moreover, they found a half-light radius of pc, and an absolute magnitude of mag. Based on all these properties, M18 concluded that the system was likely a globular cluster, and compared the satellite to several extended clusters of M31 associated to known streams, in the same way that Sgr II is suspected to belong to the Sgr stream. However, M18 emphasised the importance of a spectroscopic study to confirm their conclusion.

In this work, we show a thorough analysis of the stellar, structural and orbital properties of Sagittarius II, using deep broadband photometry from the Canadian-France-Hawaii Telescope (CFHT) MegaCam (MC hereafter) imager in the context of Pristine. The Pristine survey uses a narrow-band filter centred on the metallicity-sensitive Ca H&K doublet (Starkenburg et al. 2017) to identify metal-poor stars and estimate their metallicity using pure photometry. Keck II/DEIMOS spectroscopy are additionally used to constrain the system’s metallicity and kinematics. Finally, combined with the Gaia Data Release 2, we constrain the orbital properties of the satellite.

2 Observations and Data

2.1 Photometry

Our photometry consists of deep broadband and observations as well as narrow-band observations with the CaHK filter centred on the metallicity-sensitive Calcium HK doublet. It was observed using the wide-field imager MegaCam on the CFHT (Boulade et al., 2003). The CaHK photometry is part of a larger survey called Pristine (Starkenburg et al., 2017). The ongoing Pristine and Pristine dwarf-galaxy surveys aim to probe the northern part of the MW halo, the MW bulge, as well as many of the northern faint dwarf galaxies and dwarf-galaxy candidates. A study of the MW satellite Draco II, from a very similar data set to that of the current paper, is presented in Longeard et al. (2018, hereafter L18).

Observations were conducted in service mode by the CFHT staff during the night of July, 2nd, 2016 under good seeing conditions (). We refer to reader to L18 for the details of the Megacam data reduction. The star/galaxy separation is done using the Cambridge Astronomical Survey Unit (Irwin & Lewis, 2001) pipeline flags, which also indicate saturated sources. The MegaCam photometry is calibrated onto the PS1 photometric system with the same procedure as in L18: a cross-identification of all unsaturated point sources with photometric uncertainties below 0.05 mag in both catalogs is performed. The difference (respectively ) is expressed as a function of the colour . We then fit a third-order polynomial to translate MC photometry into PS1 through a 3 clipping procedure. The coefficients of the polynomials to transform (,) into (,) in this work are different from those in L18. We define and obtain:

with . The calibration coefficients are: , , for the band and , , for the band. All uncertainties on the polynomials coefficients are propagated into the photometric uncertainties.

All stars saturated in the MC photometry, filtered during the calibration process, are taken directly from PS1 and added to the final catalog, for a total of 83,355 stars. This catalog is finally dereddened using the dust map from Schlegel, Finkbeiner & Davis (1998) and the extinction coefficients from Schlafly & Finkbeiner (2011). In the rest of the text, we use the combined catalogue and the PS or MC subscripts are dropped.

2.2 Spectroscopy

Figure 1: Spectrum of a star in our spectroscopic data set, corresponding to a S/N . This spectrum is already processed and sky-substracted.

Spectroscopic follow-up observations of Sgr II were obtained with Keck/DEIMOS (Faber et al., 2003): “mask 1” was observed on the 2015-09-18, “mask 2” on the 2015-09-08 and “mask 3” on the 2015-09-12 . Mask 3 is a re-observation of mask 1, while mask 2 was chosen perpendicular to the other two in order to probe potential Sgr II members further away in the South/North direction (Figure 2). We refer the reader to Ibata et al. (2011) for a detailed description of our DEIMOS data reduction procedure. All stars with a signal-to-noise ratio (S/N) below 3 or with a velocity uncertainty greater than 15 km s were discarded. Following the procedure described in (Simon & Geha, 2007) and using the 51 stars observed at least twice, we assess the systematics in our sample and find a negligible bias of and a systematic uncertainty floor of . The heliocentric velocities of each star observed more than once are merged into one single measurement by taking the mean of all available quantities weighted by the inverse of their respective velocity uncertainties. The same procedure was followed for the equivalent widths of the Ca triplet. For illustrative purposes, an example of spectrum of S/N 90 is shown in Figure 1.

Finally, the existence of binaries in the sample is investigated for all stars with multiple velocity measurements. To do so, we define the quantity such that

with the heliocentric velocity of a star in the mask l (resp. for ), and the uncertainty on this measurement (resp. for ). If is greater than 2.5, the star is considered as a possible binary and flagged accordingly. Two stars, observed in mask 1 and 3, are identified as binaries trough this procedure, with differences in velocities of and .

3 Broadband photometry analysis

We present the one square degree field centred on Sgr II together with the spatial distribution of stars observed with spectroscopy in Figure 2.

Figure 2: Left panel: Spatial distribution of Sgr II-like stars, centred on (). The red contour defines the two half-light radii () of the satellite. Right panel: Close-up on the central region, with stars observed spectroscopically are colour-coded according to their heliocentric velocities. Filled circles represent very likely Sgr II members, while filled triangles represent Sgr II HB stars in the spectroscopic data set.

The colour-magnitude diagram (CMD) of all stars within two half-light radii () of the system is shown in the left panel of Figure 3, along with the spectroscopically observed stars. The CMD of the same areal coverage but selected in the outskirts of the MegaCam field of view is represented as a comparison in the middle panel. Sgr II’s main sequence (MS) and main sequence turn-off (MSTO) are very well defined thanks to the depth of the MegaCam data, and corresponds to an old ( 12 Gyr) and very metal-poor ([Fe/H] ) population. A few blue stragglers can be seen in the satellite. Sgr II also hosts a few blue horizontal branch stars around .

The BHB stars are useful as they are good distance tracers (Deason, Belokurov & Evans, 2011, D11). First, all 12 BHBs’ and are calibrated onto the SDSS photometry according to the colour equations of Tonry et al. (2012). The median absolute magnitude of the BHBs is obtained using equation (7) of D11 and yields mag, yielding a median distance modulus of mag or kpc.

Figure 3: Left panel: CMD within two of Sgr II. Its old ( 10 Gyr) and metal-poor ([Fe/H] -2.2) stellar population clearly stands out. One can notice the presence of a few blue stragglers in the system around , as well as Sgr II horizontal branch. Stars observed with spectroscopy are shown with coloured circles. The filled ones represent the confirmed spectroscopic members. Filled triangles show the location, in the CMD, of HB stars in the spectroscopic data set. The favoured isochrone for Sgr II, obtained in section 3, is shown as a solid red line (, dex, /Fe dex, mag). Middle panel: CMD of the field for an equivalent area. Right panel: Photometric uncertainties for the band. The grey contours show the mask used to select the Sgr II-like population showed in the spatial distribution of Figure 2.

3.1 Structural and CMD fitting

We use our MC photometry to refine the structural properties of Sgr II previously studied by L15 and M18 and determine its main stellar properties through a CMD and spatial distribution fitting procedure. The formalism of this analysis is detailed in Martin et al. (2016a) and L18. Though the main steps will be briefly detailed below, we refer the reader to these two references for more details. Six structural parameters will be inferred from our analysis: the centroid offsets along the X and Y axes with respect to the centre coordinates of the literature, and , in arcminutes, the ellipticity 111The ellipticity is defined as 1 - , with and the major and minor axis. , the half-light radius , the position angle , and the number of stars of the satellite. These structural properties are gathered in a parameter set noted .

We then define the CMD parameters derived by our CMD fitting procedure: the age of the satellite , the systemic metallicity , the [/Fe] abundance ratio, the distance modulus , and the fraction of Sgr II stars with respect to the total number of stars in the field chosen for the anlaysis. We regroup these properties into the set

For a given star , we consider its following properties: its and magnitudes, and its position offset from the center coordinates of the literature, and . These four properties are gathered into one set .

The Sgr II radial density, , is modelled by an exponential radial profile while the foreground contamination is assumed constant over the field of view, i.e.

(1)

with the elliptical radius, which can be expressed using the projected sky coordinates as

(2)

For the -th star, the spatial likelihood can then simply be written as

(3)

where in the area of the sky over which the analysis is conducted.

The spatial likelihood of the Milky Way foreground contamination is assumed flat, which yields

(4)

The CMD likelihood function is built from the sum of two models: one for the foreground, , constructed empirically from the field CMD, and one to describe the Sgr II population taken as a single star population, and called . The foreground contamination model is built by selecting all stars outside 5 of the system centre and binning their distribution on the CMD. This distribution is smoothed by a gaussian kernel in both colour and magnitude of a width of 0.1 in an attempt to limit the effects of shot noise. is generated using a range of Darmouth isochrones (Dotter et al., 2008). We choose isochrones of different , , [/Fe], and distance modulus . The priors on each parameters are reported in Table 1. To build the PDF of a given stellar population, we simulate the CMD of a population of several million stars, based on its isochrone, luminosity function and photometric uncertainty at a given (,). To avoid aliasing effects, especially at the bright end of our models where the photometric uncertainties are unrealistically low, we add 0.01 in quadrature to the photometric uncertainties. Finally, each PDF is degraded to the completeness of the data estimated by following the method of Martin et al. (2016a). The 50 per cent completeness is reached at and mag.

The structural and CMD parameters are gathered into a single set . At the distance of Sgr II, the tip of the Red Giant Branch (RGB) tip is expected to be at . Furthermore, misidentified background galaxies start to pollute our photometry below . The fit does not take into account the horizontal branch stars as these are poorly modelled by the theoretical stellar population models. Therefore, the analysis is restricted in a specific CMD box defined with the following cuts: and . CMD and spatial properties are fitted at the same time through our own Markov Chain Monte Carlo (Hastings, 1970, MCMC) algorithm by maximising the likelihood of the following model:

(5)

with

(6)
(7)

Finally, the distance to Sgr II is constrained using a Gaussian prior based on the distance modulus derived from the median absolute magnitude of the BHBs in the first paragraph of section 3 ( mag). A Gaussian prior on the metallicity of the satellite is also applied and comes directly from the combination of the spectroscopic and CaHK metallicity measurements detailed in the sections 4 and 5 respectively ( dex). The inference of each parameter of is summed up in Table 1, and the 2D PDFs are shown in Figure 4.

Parameter Unit Prior Favoured model Uncertainties
Right ascension (ICRS) degrees 298.166284027
Declination (ICRS) degrees
degrees 18.93203283
degrees
arcmin 1.7
pc 35.5
degrees [0,180] 103
0.0 at the % CL
Distance modulus mag [18.90,19.45] 19.32
Distance kpc 73.1
Gyr [9,13.5] 12.0 0.5
dex [-2.4,-1.5] 0.03
dex [-0.2,0.6] 0.0
log(Luminosity)
mag -5.7 0.1
mag arcsec
0
mas yr
mas yr
Apocenter kpc
Pericenter kpc
0.44
0.4
160.3
km s kpc 6292
km s 17159
Table 1: Properties of Sgr II.

Figure 4: One and two-dimensional PDFs of the structural and CMD properties of Sgr II, inferred using the method described in section 3.1. Red lines correspond to the favoured inference, using both the distance prior based on BHBs and the metallicity prior from the spectroscopy, while the black contours show the case without any prior. These contours are defined as the usual 1, 2 and 3 confidence intervals in the 2D, gaussian case.

Figure 5: CMD within two half-light radii of Sgr II overplotted with the fiducials of MW globular clusters from Bernard et al. (2014). Stars with a Sgr II membership probability greater than ten per cent are shown as red circles. Spectroscopically confirmed members of Sgr II are shown as yellow diamonds. The fiducials are separated in metallicity bins, from the most metal-poor in dark blue to the most metal-rich available in orange, and shifted to the favoured distance modulus found for Sgr II.

The best-fit isochrone is shown as the red PDFs in Figure 4. Sgr II is found to be significantly old and metal-poor with an age of Gyr population along with a systemic metallicity of dex. Furthermore, the alpha abundance of this isochrone is solar ([/Fe] = 0.0), though we caution the reader about reading too much into this parameter given the roughness of the [/Fe] abundance ratio grid. Finally, the favoured distance modulus is mag, and corresponds to a physical distance of kpc. We compare these results by performing the fit without the BHBs or the spectroscopic metallicity priors. For this case, the PDFs are shown in black in Figure 4. The inferences of all the parameters are compatible: the stellar population is here found to be older (), more metal-poor (, at the lower edge of the metallicity grid) and with a distance modulus of = mag. All structural properties are perfectly compatible with L15 and M18: Sgr II is consistent with being spherical ( at the 95% CL) and has a size of arcminutes, translating into a physical size of 35.5 pc. All the main properties of Sgr II are summarised in Table 1.

Finally, we investigate the presence of RR Lyrae in the field by cross-identifying the PS1 RR Lyrae catalog of Sesar et al. (2017) with our photometry. Three RRLyrae are found in the vicinity of Sgr II. Two of these have similar distance modulii, as inferred from Sesar et al. (2017) (18.73 and 18.85 mag). However, the resulting distances are discrepant from both our BHB and CMD fitting analyses by 0.5 mag (roughly 10 kpc in physical distance). To confirm the distance modulus of Sgr II, we compare the CMD of the satellite with fiducials of MW globular clusters in PS1 (Bernard et al., 2014) in Figure 5. In this plot, all fiducials are deredenned and their distance modulii are corrected to correspond to our favoured model for Sgr II ( mag). The most metal-poor fiducials are undoubtedly the ones that better reproduce the features of the CMD of Sgr II, including the HB. Shifting the fiducials to the distance of the two RR Lyrae in the field would allow Sgr II features to be best reproduced by the more metal-rich fiducials in green. However, Sgr II’s metallicity is clearly lower than that as shown by our CMD fitting as well as the spectroscopic and CaHK analyses respectively in section 4 and 5. Therefore, the distance implied by the two RR Lyrae is clearly not that of Sgr II. A plausible origin for these two stars might just be the Sgr stream, as shown in section 6.

3.2 Luminosity

Figure 6: PDFs of the Sgr II luminosity and absolute magnitude. The favoured luminosity of the satellite is , corresponding to an absolute magnitude of mag.

The luminosity, absolute magnitude, and surface brightness of Sgr II are derived using the formalism of Martin et al. (2016a). The first step consists of drawing a set of parameters denoted from the final multi-dimensional PDF obtained through the analysis of section 3.1. These parameters are the number of stars N, an age , the metallicity , the alpha abundance ratio [], and the distance modulus . A CMD of the -th stellar population is then simulated; for each simulated star, we ensure that its location in the colour-magnitude diagram falls in the CMD box used to perform the structural and CMD fit ( and ). Furthermore, a completeness test is performed: the completeness of the survey is first estimated at the colour and magnitude of the simulated star. Then, two random numbers and between 0 and 1 are drawn: if the completenesses of the star in both and is greater than these numbers, it is flagged. When the number of flagged stars reaches , the fluxes of all simulated stars, flagged or not, are summed, which gives the luminosity of the satellite for the -th iteration. This procedure is repeated a thousand times in order to have visually pleasant PDFs.

The 1D marginalised PDFs of Sgr II’s luminosity and absolute magnitude are represented in Figure 6. The final favoured parameters are reported in Table 1. The luminosity of the satellite is found to be . This measurement is in agreement with both L15 and M18 ( ). Finally, we obtain a surface brightness of mag arcsec.

4 Narrow-band CaHK photometry analysis

Figure 7: Normalised distribution of Pristine photometric metallicities for all stars within 2r (solid red line). The same histogram is also shown for all field stars, i.e. stars outside 5r (black dashed line). Sgr II clearly peaks at [Fe/H] dex, while no such overdensity exists for the field distribution.

The CaHK photometry can be used to estimate the metallicity of Sgr II and its metallicity dispersion. The photometric metallicity of each star from (,,) can be estimated using the model detailed in Starkenburg et al. (2017). Pristine observations are shallower than our broadband and photometry (right panel of Figure 3) and therefore can only be used to estimate the photometric metallicity [Fe/H] of stars in our field down to mag. Starkenburg et al. (2017) show that the Pristine metallicities are slightly biased low as we go toward the metal-poor end of the calibration model. Therefore, we repeat the procedure presented in L18 and we first correct for this effect by binning in metallicity the sample used by Starkenburg et al. (2017), which provides both the SDSS spectroscopic metallicity and the Pristine photometric metallicity for several thousands stars. For each bin, the median difference between the SDSS and Pristine metallicities is computed. This procedure yields the bias as a function of the photometric metallicity, which is used to correct our whole Sgr II metallicity catalog. All stars with [Fe/H] or [Fe/H] are discarded as our Pristine model is not reliable for those stars (Youakim et al., 2017). This choice does not affect the analysis as the systemic metallicity of Sgr II is well within this range. Stars with a large uncertainty in the photometry () are rejected. All remaining stars within 2 are selected and their photometric metallicity distribution function (MDF) is shown in Figure 7 with the solid red line. The distribution of all field stars within is shown as the black dashed line for comparison.

Sgr II stars in red stand out clearly in Figure 7 as they form a pronounced peak around [Fe/H] dex that does not exist in the MDF of the field stars in black. To derive Sgr II’s metallicity properties, we assume that the population encapsulated inside 2r (corresponding to 206 stars) in Figure 7 can be modelled as the sum of the foreground MDF and a normally distributed photometric metallicity population associated with Sgr II stars. This assumption seems legitimate as the metallicity distribution at the metal-rich end of the red histogram in Figure 7 overlaps well with the black distribution, thus implying that the underlying foreground contamination MDF is comparable over the field of view. The Sgr II stellar population metallicity distribution is assumed to be normally distributed, with a mean [Fe/H] and a standard deviation of with being the intrinsic metallicity dispersion of Sgr II and the uncertainty on the photometric metallicity of the -th star. The likelihood model for the MW contamination stars is built by interpolating the [Fe/H] MDF of all stars outside 5r. This model is then smoothed by a gaussian kernel of 0.2 dex to account for poor statistics in some metallicity bins. The fit is performed through a MCMC algorithm, and we marginalise over the foreground contamination model. At each iteration, we randomly draw a photometric metallicity for all stars in the contamination subsample, according to their individual photometric metallicity uncertainties. Then, the procedure to build the foreground contamination model described above is repeated. In doing so, the analysis takes into account the overall uncertainty of the contamination MDF.

The 39, 88 and 95 % volume intervals on the final 2D posterior PDF, corresponding to the , and confidence levels for the 2D gaussian case, are shown in black solid line in Figure 8. We measure a significant, non-zero metallicity spread in Sgr II, with dex and to be particularly metal-poor ([Fe/H] dex), in agreement with the stellar population inferred through the CMD fitting. To ensure that this inference is not caused by any systematic effect, the same analysis was done in L18 with two metal-poor globular clusters in the Pristine footprint, M92 ans M15. The systemic metallicities of both clusters were found to be compatible with their previous estimates using spectroscopic data. Furthermore, their metallicity dispersions were unresolved, in agreement with previous studies, showing that the technique does not seem to be affected by a systematic effect.

Figure 8: Two-dimensional joint PDFs of the systemic metallicity and dispersion for Sgr II using the photometric CaHK metallicities (black) and the individual spectroscopic metallicities of member stars (grey). These two independent measurements are combined to give the final PDF shown as the dashed red line. The contours represent the 39, 88 and 95 % volume intervals. The associated one-dimensional marginalised PDFs for all cases are shown in the upper and right panels. Both methods are in agreement and shows that Sgr II has a small but measurable metallicity dispersion.

5 Spectroscopic analysis

5.1 Velocity properties

Figure 9: Heliocentric velocity histograms for the three spectroscopic samples. The fourth panel shows the merging of all samples. The grey histograms show the number of stars with a non-reliable photometric metallicity measurement in our sample that therefore cannot be filtered out by our technique. The peak of Sgr II stars around is pronounced and the disc contamination, from 0 to , is also quite populated. The last panel shows the final spectroscopic catalog filtered from “metal-rich” stars using photometric metallicities based on our CaHK photometry. For stars with reliable photometric metallicities, the ones with [Fe/H], i.e. compatible with Sgr II metallicity properties measured in section 3 and 5, are selected, while the others are discarded. Stars with mediocre quality CaHK measurement or [Fe/H] uncertainties are not discarded as their [Fe/H] is not reliable. The disc and halo contamination is notably reduced, and one star located in the Sgr II velocity peak is identified as a contaminant and removed from the sample.

Figure 10: Bottom right panel: Radial distances to the center of Sgr II with respect to the heliocentric velocities. Open circles represent stars that are metal-poor using their photometric metallicities, or stars with non-reliable photometric metallicity measurements. Black-filled dots represent member stars. Bottom left panel: Spectroscopic metallicities with respect to the heliocentric velocities of all stars from the final spectroscopic data set with S/N 12. Top panel: Two-dimensional joint PDF of the systemic spectroscopic metallicity and metallicity dispersion. The contours represent the local 39%, 88% and 95% volume intervals. Sgr II comes out as a very metal-poor satellite, with [Fe/H] dex, and seems chemically enriched: dex.

Figure 11: Pristine colour-colour diagram: the colour is represented on the x-axis, while the metallicity information is carried by the colour on the y-axis. Grey dots stand for all field stars in a range between 5 and 8 half-light radii to Sgr II. Most of them are halo and disc stars and form a stellar locus of more metal-rich stars ([Fe/H] 1 or above). Stars observed with spectroscopy are represented with circles colour-coded according to their heliocentric velocities. Among those, filled circles show the stars identified as spectroscopic members, while filled triangles stand for the HB stars in the spectroscopic data set. Above the grey stellar locus are located stars that become more and more metal-poor as we go towards the upper part of the diagram. Two iso-metallicity sequences are shown in red and green dashed lines, corresponding respectively to a photometric metallicity of [Fe/H] -1.6 and [Fe/H] -4.0. As expected, most of the stars in cyan, with a radial velocity compatible with Sgr II, are located in the metal-poor region enclaved by those two sequences. Hence, only stars within this region are selected for the final spectroscopic sample. Furthermore, we add a criteria on and discard all stars between in order to discard potential foreground white dwarfs.

Figure 12: Marginalised PDFs for the systemic velocity (left) and its associated dispersion (right) for Sgr II. The satellite is found to be dynamically cold: the velocity dispersion is only marginally resolved: , reaching 6.5 at the 95% confidence level.

The systemic velocity and velocity dispersion are derived using the deep spectroscopic observations of the system. The heliocentric velocity distributions of each mask separately are summarised in the top three panels of Figure 9, and are combined to obtain one single velocity distribution shown in the fourth panel of the same figure. We present the radial distance of each star with respect to their radial velocities in Figure 10. The full dataset is detailed in Table 2.

The velocity peak of Sgr II stands out at around km s while contaminating MW stars are distributed sparsely all over the velocity space, and in particular can be located in the vicinity of the Sgr II velocity peak. Because of the small number of stars in the Sgr II population, the velocity properties can be polluted by the contamination. Ideally, those stars would have to be identified and discarded from the spectroscopic sample. Sgr II is an old and metal-poor system as suggested by its CMD, and confirmed in section 3, 4 and 5, whereas the contaminating foreground is expected to be more metal-rich overall. Therefore, the contamination stars could be discarded based on their metallicities. Even though the individual [Fe/H] are only accessible for the brightest stars in our sample in our spectroscopic sample, the Pristine photometry can be used here.

The Pristine colour-colour diagram is shown in Figure 11. Field stars, i.e. a randomly selected sample of all stars outside five half-light radii, are represented in small black dots and form a clear stellar locus. This diagram is constructed so that the individual metallicity of a given star decreases from the bottom right to the top left, while more metal-rich stars occupy the lower part of this colour-colour space. Stars observed with spectroscopy are colour-coded according to their heliocentric velocities, provided they pass the following CaHK photometry and [Fe/H] quality cuts: and . Stars that do not match these criteria are not discarded from the final spectroscopic sample, whether they are more metal-rich or not in the model, as our goal here is to select only stars for which the CaHK photometry is reliable enough to ensure their metal-poor nature. In section 4, we found that Sagittarius II has a systemic metallicity of [Fe/H] dex and has a resolved metallicity dispersion. Therefore, within the subsample of stars that passed the CaHK photometry cuts discussed above, we choose to select only stars with , as a Sgr II-like system would likely have a star formation history too short to produce significantly more metal-rich stars. The region of the diagram that corresponds to such a metallicity cut is represented by the two iso-metallicity sequences in green and red dashed lines in Figure 11. Two cuts in are also applied in order to discard potential white dwarfs and metal-rich stars. The final spectroscopic velocity distribution is shown in the last panel of Figure 9. A significant number of MW stars with a reliable Pristine photometric metallicity measurements are cleaned out from the catalog as their metallicities are too high to be members of Sgr II, even if the satellite has a metallicity spread. In particular, one star in the immediate vicinity of the Sgr II velocity peak is identified as a more metal-rich contaminants using this technique ( dex) and therefore discarded.

The resulting velocity distribution is assumed to be the sum of three normally distributed populations: one for Sgr II stars, and two others corresponding to the MW foreground disc and halo stars. Each individual likelihood is weighted by its CMD and structural probability membership determined in section 5. One can write the individual likelihood of the -th star as

(8)

with and the individual velocity uncertainty of the k-th star, the intrinsic velocity dispersion, the systematic threshold derived in section 2.2. is the systemic velocity of Sgr II. and are the fractions of stars respectively in the MW disc and halo populations. is defined as , with the intrinsic velocity dispersion of the disc population. The corresponding quantity for the halo population is written , while is the systemic velocity of the disc population in the sample (resp. for the halo population). is the usual one-dimensional normal distribution. We run a MCMC analysis and show the resulting marginalised 1D PDFs in Figure 12. At each iteration of the MCMC, the systematic threshold is randomly drawn from its PDF. The favoured systemic velocity is km s. The velocity dispersion of Sgr II is km s, reaching 6.5 km s at the 95% confidence interval, thus showing that Sagittarius II is a dynamically cold satellite. A similar analysis was performed for the inner (r 1 arcmin) and outer (r 1 arcmin) regions and no statistical difference in terms of velocity dispersion was found.

5.2 Metallicity properties

To infer the metallicity properties of Sgr II from the spectrocopy, we create a subsample constituted of stars brighter than and a S/N ratio above 12 from our final spectroscopic sample, for a total of 26 stars. The spectroscopic metallicity is estimated using the calibration from Starkenburg et al. (2010) based on the Ca triplet. This method is originally calibrated for RGB stars, however, Leaman et al. (2013) showed that it can applied to stars up to two magnitudes fainter (see also Carrera et al. 2013). The resulting sample consists of six member stars, for which the individual spectroscopic metallicities are reported in Table 2 under “Fe/H]”. The distribution of spectroscopic metallicities with respect to the radial velocity is shown in the bottom left panel of Figure 10, and shows the existence of a clump of stars at around dex at the velocity of Sgr II.

To derive the systemic metallicity and metallicity dispersion of Sgr II, we assume that the spectroscopic metallicity of Sgr II stars are normally distributed and weigh each star with its CMD and structural probability membership, giving the following likelihood function

(9)

with , being the individual uncertainty on the spectroscopic metallicity of the -th star, and the intrisic metallicity dispersion of Sagittarius II. The 39, 88 and and 95 volume intervals are represented by black solid lines on the resulting 2D probability distribution functions (PDF) in the top left panel of Figure 10. Sgr II is confirmed to be metal-poor, with [Fe/H] = -2.23 0.05 dex. Moreover, we find a metallicity dispersion of dex. This spread in metallicity is driven by the two brightest stars identified as members of Sgr II, for which the spectroscopic metallicity is accurately measured. They have respectively a spectroscopic metallicity of dex and dex. Furthermore, since they are among the stars that were observed multiple times in our catalog, it is possible to infer their individual spectroscopic metallicities using the Ca triplet equivalent widths of each run separately. For both stars, the metallicities obtained from each spectroscopic run in which they were observed are consistent with one another, suggesting that their final are not driven by one spurious equivalent widths measurement in one of the three spectroscopic samples. In addition with being consistent with the CMD of Sgr II and its systemic velocity, the two stars are also remarkably compatible with the satellite’s proper motion inferred in section 6. Taken all together, we favour the fact that these two stars are indeed members of Sgr II, and there is more than one stellar population in the system.

The two independent measurements of the metallicity and dispersion of the satellite, using the observations on the one hand and the spectra on the other, are perfectly compatible. The results of both methods are then combined into one single measurement by performing the product of the two 2D joint PDFs. We show the corresponding 39, 88 and 95 % volume intervals in red thick line in Figure 8. This final measurement yields a systemic metallicity of [Fe/H] = -2.28 0.03 dex and a metallicity dispersion of dex.

6 Gaia DR2 proper motions and orbit

Figure 13: Distribution of the proper motions of MW contamination stars in small grey dots, and our Sgr II-like population shown with dots colour-coded according to their proper motion membership probability, derived from a gaussian mixture model. The seven Sgr II members with a proper motion measurement in Gaia are represented with diamonds and the twelve HB stars with triangles. The systemic proper motion of Sgr II ( mas.yr, mas.yr) is represented with a large red dot.

Figure 14: Projections of the orbit of Sgr II on the X-Y, X-Z and Y-Z planes backwards and forwards over 2.0 Gyr. Twenty-one orbits are shown here: the one based on the favoured position, distance, radial velocity and PMs of the satellite (as the thick, darker grey line), and twenty others using random realisations of those parameters (as thin, slightly transparent grey lines). The red circle is the current position of Sgr II, while the magenta triangle is the one of the Sgr dSph. A N-body simulation of the trailing arm of the Sgr stream (Law & Majewski, 2010) is shown in green. The MW disk is shown in black, with a chosen radius of 15 kpc.

To infer the orbit of Sgr II, we first build a sample of Sgr II-like population based on the mask shown in the right panel of Figure 3. The proper motions of those stars are retrieved from the Gaia Data Release 2 (Gaia Collaboration et al., 2018). All member stars identified with spectroscopy and bright enough to have a proper motion measurement in Gaia are naturally present in this sample. Furthermore, the Gaia DR2 data are also cross-matched with the potential HB stars within two half-light radii of the satellite. Twelve HB stars have a proper motion measurement in Gaia and are added to the sample shown in Figure 13.

The inference of the Sgr II proper motion is performed with a gaussian mixture model. We assume that the sample can be modelled by the sum of two bivariate gaussians: one for the Sgr II population and another for the foreground MW contamination. The sets of parameters inferred from the analysis are composed of the proper motions in both directions, their dispersions and correlation , for Sgr II (,,,,) and for the contamination (,,,,). The proper motion properties of the -th star are defined as with the uncertainty on the proper motion in the RA direction (respectively for DEC). The individual likelihood is

(10)

where is a two-dimensional gaussian and the fraction of Sgr II stars in the sample.

The gaussian mixture model gives a systemic proper motion of mas yr and mas yr for Sgr II. These proper motions take into account the systematic error on the one for dSph derived by Gaia Collaboration et al. (2018b). We also inferred the proper motion of the system using the HB and spectroscopic member stars only, and found a compatible result with (,) = (,) mas yr.

Our estimate is discrepant from the one of Massari & Helmi (2018) who find a proper motion of (,) = (,) mas yr. They rely on the convergence of the astrometric parameters through a clipping procedure, with an initial guess on those parameters based on the potential HB stars of Sgr II. However, our measurement based only on HB and spectroscopic member stars gives credit to the proper motion found in this work, and disfavours the estimate of the work of Massari & Helmi (2018), which might be biased by the foreground contamination.

The orbit of the satellite can then be inferred using the GALPY package (Bovy, 2015). The MW potential chosen to integrate the orbit is a modified “MWPotential14” constituted of three main components: a power-law, exponentially cut-off bulge, a Miyamoto-Nagai Disc, and a NFW DM halo with a virial mass of . Further details about this MW potential model can be found in Bovy (2015). We integrate 2000 orbits backwards and forwards, each time by randomly drawing a position, distance, radial velocity and proper motions from their respective PDFs, and extract for each realisation the pericenter, apocenter and ellipticity. Each orbit is shot over 2 Gyr. The favoured orbit (i.e. the favoured position, distance, radial velocity and PMs) is shown in Figure 14 in the X-Y, X-Z and Y-Z planes, along with the Sgr stream. Twenty other random realisations of Sgr II orbits are also shown in grey, partially transparent lines.

The analysis yields a pericenter of , an apocenter of and an orbital ellipticity of . Moreover, Figure 14 shows that the orbit of Sgr II is compatible with the trailing arm of the Sgr stream, despite being slightly tilted from it, especially in the Y direction.

7 Discussion

Figure 15: Comparison of Sgr II with other GCs and dwarf galaxies of the Milky Way. Squares represent dwarf galaxies while circles represent globular clusters, and the diamond corresponds to Dra II. Triangles stand for recently discovered dwarf-galaxy candidates that await confirmation. Hollow markers correspond to systems for which no metallicity dispersion measurement can be found in the literature. The solid line in the top-left panel corresponds to the luminosity-metallicity relation of Kirby et al. (2013) for dwarf spheroidals and dwarf irregulars. Dashed lines represent the RMS about this relation, also taken from Kirby et al. (2013). Among the 123 globular clusters presented here, the properties of 116 were extracted from Harris (1996) catalog, revised in 2010. For the remaining ones (Kim 1, Kim 2, Kim 3, Laevens 1, Balbinot 1, Munoz 1 and SMASH 1) parameters of the discovery publications were used (Kim & Jerjen (2015), Kim et al. (2015), Kim et al. (2016), Laevens et al. (2014), Balbinot et al. (2013), Muñoz et al. (2012) and Martin et al. (2016b)). Globular cluster metallicity spread measurements are taken from Willman & Strader (2012) and references therein: Carretta et al. (2006, 2007, 2009, 2011), Cohen et al. (2010), Gratton et al. (2007), Johnson & Pilachowski (2010), and Marino et al. (2011). McConnachie (2012) and Willman & Strader (2012) are used to compile the properties of the dwarf galaxies represented here. The 18 dwarf galaxies represented here are: Bootes I (Belokurov et al., 2006; Norris et al., 2010), Canes Venatici I (Zucker et al., 2006b), Canes Venatici II (Sakamoto & Hasegawa, 2006), Coma Berinices, Hercules, Leo IV and Segue I (Belokurov et al., 2007), Draco and Ursa Minor (Wilson, 1955), Fornax (Shapley, 1938b), Leo I and Leo II (Harrington & Wilson, 1950), Pisces II (Belokurov et al., 2010), Sculptor (Shapley, 1938a), Sextans (Irwin et al., 1990), Ursa Major I (Willman et al., 2005b), Ursa Major II (Zucker et al., 2006a), Willman I (Willman et al., 2005a). Their metallicity and metallicity spreads were drawn from Kirby et al. (2008), Kirby et al. (2010), Norris et al. (2010), Willman et al. (2011). The dwarf galaxy candidates discovered recently and shown on this figure are Bootes II (Koch & Rich, 2014), DES1 (Luque et al., 2016; Conn et al., 2018), Eridanus III (Bechtol et al., 2015; Conn et al., 2018; Koposov et al., 2015b), Hyades II (Martin et al., 2015), Pegasus III (Kim & Jerjen, 2015), Reticulum II and Horologium I (Koposov et al., 2015a), Segue II (Belokurov et al., 2009), and the most significant candidates of Drlica-Wagner et al. (2015): Gru II, Tuc III, and Tuc IV.

We used deep MegaCam broadband photometry, the narrow-band Pristine survey and DEIMOS spectroscopy to conduct a thorough study of the Milky Way satellite Sgr II. By performing a CMD and structural analysis, the satellite is found to have a size of 35.5 pc, and is located at kpc based on the combination of BHB stars distances and a CMD fitting procedure. The favoured stellar population is old ( Gyr) and metal-poor. Using our spectroscopic catalog, we are able to find the systemic velocity of Sgr II to be . The velocity dispersion yields km s and is at the 95% confidence interval. From this spectroscopic analysis, 22 stars are identified as members of the satellite and reported in Table 2. The individual photometric metallicities provided by the Pristine survey are used to show that Sgr II is a very metal-poor system, with [Fe/H] dex, and has a resolved metallicity dispersion: dex. These two chemical properties are perfectly supported by our spectroscopic analysis of six RGB stars. We applied the Ca triplet calibration from Starkenburg et al. (2010) to derive Sgr II spectroscopic metallicity and its associated dispersion: [Fe/H] dex and dex. Combining the CaHK and spectroscopic measurements, we obtain refined estimates of both parameters: [Fe/H] = dex and dex. Finally, using the Gaia DR2 data, the proper motion of Sgr II is inferred to be (,) = (,) mas yr. This yields an apocenter and pericenter of and kpc respectively.

Sgr II is in perfect agreement with the luminosity-metallicity relations for dwarf galaxies (Kirby et al., 2013), as it is shown in the top left panel of Figure 15. Sgr II is however somewhat of an outlier in the plane (bottom-left panel of Figure 15), which led M18 to conclude that Sgr II is a globular cluster. However, the locus of dwarf galaxies in this plane becomes uncertain at low luminosities. The satellite is still more extended than the vast majority of MW globular clusters as shown in the bottom left panel of Figure 15, although two of them have a comparable size: Crater (Belokurov et al., 2014; Laevens et al., 2014) and Terzan 5 (Terzan, 1968). These two extended clusters do not, however, share the same metallicity properties as Sgr II: Terzan 5 is a bulge cluster with [Fe/H] and Crater is more metal-rich with a systemic metallicity of [Fe/H] (Weisz et al., 2016). Our two estimates of the metallicity dispersion of Sgr II both yield similar results and suggest that the satellite was able to retain its gas and form successive generation of stars, thus suggesting the presence of a dark matter halo (Willman & Strader, 2012). However, this result is driven by two bright RGB stars that have discrepant metallicity measurements. If one of the two were misidentified as a Sgr II member, the claim of a metallicity dispersion would be weaker.

The question of the dynamical mass of Sgr II remains open. We can use the relation of Walker et al. (2009) to estimate the expected velocity dispersion of a purely baryonic system. Assuming a mass-to-light ratio of 2 for an old and metal-poor stellar population (McLaughlin & van der Marel, 2005), Sgr II would have a velocity dispersion of , which is not incompatible with our inference of . Nevertheless, taken at face value, our velocity dispersion measurement implies that Sgr II has a dynamical mass-to-light ratio of and favours a slightly DM-dominated system under the usual assumption of dynamical equilibrium and sphericity. If this is confirmed, it would mean that Sgr II inhabits one of the lowest mass DM subhalo. Alternatively, this result could be driven by the compactness of the satellite, whose stars only probe the inner parts of the subhalo.

Taken together, these two pieces of evidence (marginally resolved metallicity dispersion and plausibly non-baryonic M/L) would indicate that Sgr II is more likely a dwarf galaxy rather than a cluster.

Before the submission of this work, a spectroscopic study of Sgr II was presented at the AAS iPoster222<https://aas233-aas.ipostersessions.com/default.aspx?s=E7-10-7C-92-5D-B1-84-24-1F-B5-07-1A-BF-2E-10-65> session (Simon et al., 2019). Using Magellan/IMACS spectroscopy, they found a systemic velocity and metallicity compatible to the ones in this work: and dex. Their velocity dispersion is also consistent with ours: . Finally, the proper motion they derive for Sgr II ((,) = (,) mas yr) is also compatible with our work. However, they estimate a very low metallicity dispersion, with dex at the 95 per cent confidence limit. Therefore, they conclude that the satellite is a globular cluster. Once the two data sets are made public, a thorough investigation is needed to understand the source of this discrepancy. Anyhow, it illustrates the difficulty of studying and understanding such faint systems.

Independently of the nature of Sgr II, the orbit we infer for the satellite is compatible with the trailing arm of the Sagittarius stream according to the model by Law & Majewski (2010) (Figure 14). However, we note that the agreement between the two orbits is not perfect and, in particular, that the position of Sgr II today and its favoured movement in the Y-Z galactocentric plane are slightly offset from the plane of the Sgr stream. Three hypotheses can be formulated to explain this discrepancy:

  • The fact that the Sgr stream and the Sgr II orbits are compatible is purely coincidental.

  • Sgr II is linked to the stream, and the discrepancy between Sgr II and the stream in the Y direction, if real, could be explained by the fact that Sgr dSph satellites were stripped first and with a different energy than that of stars represented in the simulation.

  • Sgr II is linked to the stream and is also representative of its behaviour around the MW. No model is able to match all the observational constraints existing for the Sgr stream (Fardal et al., 2019). The observed difference in the orbital plane of Sgr and Sgr II could suggest that the behaviour of the distant Sgr stream wrap that Sgr II would be associated to is not perfectly described by the Law & Majewski (2010) simulation.

If either the second or the third scenario is the valid one, it would mean that Sgr II is a new, exciting example of satellite of a satellite. Similarly to the Magellanic Clouds, the Sgr dSph would then have brought its own cohort of satellites that have now been deposited in the MW halo. Moreover, it would also bring some precious insights on the orbit of the Sgr stream in regions where it is poorly constrained.

Acknowledgments

We gratefully thank the CFHT staff for performing the observations in queue mode, for their reactivity in adapting the schedule, and for answering our questions during the data-reduction process. We thank Nina Hernitschek for granting us access to the catalogue of Pan-STARRS variability catalogue.

ES, KY, and AA gratefully acknowledge funding by the Emmy Noether program from the Deutsche Forschungsgemeinschaft (DFG). This work has been published under the framework of the IdEx Unistra and benefits from a funding from the state managed by the French National Research Agency as part of the investments for the future program. NFM, RI, and NL gratefully acknowledge support from the French National Research Agency (ANR) funded project “Pristine” (ANR-18-CE31-0017) along with funding from CNRS/INSU through the Programme National Galaxies et Cosmologie and through the CNRS grant PICS07708. The authors thank the International Space Science Institute, Berne, Switzerland for providing financial support and meeting facilities to the international team “Pristine”. JIGH acknowledges financial support from the Spanish Ministry project MINECO AYA2017-86389-P, and from the Spanish MINECO under the 2013 Ramón y Cajal program MINECO RYC-2013-14875.

Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, the Institut National des Sciences de l’Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii.

Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. Furthermore, the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain.

The Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE).

This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

RA (deg) DEC (deg) (mas.yr) (mas.yr) S/N [Fe/H] [Fe/H] Member
298.20599167 21.98790000 18.53 0.01 17.99 0.01 19.44 0.01 16.6 1.3 2.949 0.506 -6.861 0.274 29.2 0.00 N
298.13157500 21.98582778 19.00 0.01 18.54 0.01 19.77 0.01 27.5 1.5 -5.194 0.738 -4.88 0.384 22.5 0.00 N
298.13137917 21.98273333 19.16 0.01 18.50 0.01 20.13 0.01 -39.4 2.2 0.336 0.69 -7.317 0.38 25.8 1.39 0.10 0.00 N
298.15507917 21.98054167 19.34 0.01 18.66 0.01 20.52 0.02 12.5 1.7 1.159 0.889 -8.058 0.507 21.2 0.00 N
298.15097500 21.95283333 19.46 0.01 18.90 0.01 20.39 0.02 140.9 1.7 -6.681 1.002 0.739 0.503 21.5 1.02 0.15 0.00 N
298.19072083 21.96758056 19.88 0.01 19.54 0.01 20.19 0.01 -118.1 6.2 -0.867 1.551 1.216 0.81 14.4 3.61 0.17 -2.85 0.23 0.00 N
298.18688333 21.97713889 20.43 0.01 19.87 0.01 21.41 0.03 -107.3 2.8 -5.917 2.251 -10.725 1.175 12.6 -1.26 0.12 0.00 N
298.18320417 21.96361944 21.28 0.01 20.73 0.01 22.22 0.06 -15.4 3.5 5.6 0.30 N
298.13665417 21.97423889 21.41 0.01 20.90 0.01 22.18 0.06 -7.8 7.4 4.7 1.40 0.30 0.04 N
298.19440000 21.99070556 21.60 0.01 21.22 0.02 22.17 0.06 -67.3 3.8 3.9 1.45 0.35 0.00 N
298.19622500 21.99368889 22.40 0.02 21.66 0.02 23.62 0.18 -2.0 13.4 2.8 0.00 N
298.16959583 22.17463056 17.71 0.01 16.94 0.01 19.06 0.01 89.7 1.1 -3.409 0.228 -5.191 0.172 31.4 0.00 N
298.18001250 22.07175000 17.90 0.01 17.05 0.01 19.32 0.01 0.2 0.9 1.896 0.225 -7.431 0.133 33.5 0.00 N
298.16810000 22.18830000 18.27 0.01 17.65 0.01 19.27 0.01 105.1 1.3 -3.679 0.287 -3.852 0.18 37.2 1.01 0.11 0.00 N
298.19661667 22.14635278 18.73 0.01 18.31 0.01 19.38 0.01 29.2 1.8 -0.708 0.441 -4.393 0.262 21.1 1.30 0.13 0.00 N
298.17483750 22.16606389 18.82 0.01 18.33 0.01 19.65 0.01 -54.0 1.3 1.647 0.507 -3.121 0.37 27.6 0.00 N
298.19189583 22.19686944 18.90 0.01 18.34 0.01 19.79 0.01 65.8 4.0 -1.137 0.558 -3.422 0.38 36.6 1.15 0.12 0.00 N
298.19827500 22.14498333 18.94 0.01 18.35 0.01 19.77 0.01 89.3 1.3 -1.161 0.478 -7.535 0.281 25.3 1.54 0.13 0.00 N
298.16145833 22.08265556 18.98 0.01 18.20 0.01 19.83 0.01 -182.8 0.9 -0.27 0.46 -0.825 0.269 26.3 2.62 0.13 -2.27 0.04 0.98 Y
298.16120000 22.00828611 19.04 0.01 18.41 0.01 20.13 0.01 -85.1 1.4 3.636 0.626 -8.935 0.343 26.7 0.00 N
298.15403750 22.11108333 19.17 0.01 18.71 0.01 19.92 0.01 51.4 1.9 1.085 0.593 -4.926 0.335 17.5 1.02 0.15 0.00 N
298.15852917 22.05846944 19.42 0.01 19.59 0.01 19.67 0.01 -167.9 9.5 0.865 1.014 -0.73 0.571 5.0 -0.43 0.29 0.00 HB
298.19296667 22.02218889 19.44 0.01 18.93 0.01 20.21 0.01 122.6 1.8 -0.935 0.985 -1.78 0.535 17.8 1.39 0.13 0.00 N
298.15097500 22.07938611 19.56 0.01 19.17 0.01 19.94 0.01 -174.0 2.1 -1.205 0.805 -0.816 0.463 16.9 3.16 0.29 -2.57 0.14 0.00 N
298.17763750 22.04601111 19.62 0.01 19.83 0.01 19.63 0.01 -135.3 4.2 -4.001 1.206 -1.43 0.651 11.1 -0.92 0.28 0.00 HB
298.17252500 22.07410833 19.66 0.01 19.90 0.01 19.67 0.01 -76.1 10.3 -0.958 1.284 -2.143 0.663 11.7 -0.05 0.23 0.00 HB
298.12763750 22.17288611 19.75 0.01 19.37 0.01 20.20 0.01 -15.5 2.6 2.824 0.831 -6.75 0.504 15.0 2.38 0.29 -1.77 0.16 0.00 N
298.20524167 22.02750556 19.83 0.01 19.01 0.01 21.05 0.02 201.1 3.4 -3.185 1.126 -15.463 0.607 17.8 1.27 0.13 0.00 B
298.18212917 22.05709444 19.85 0.01 19.04 0.01 21.08 0.03 28.4 1.8 -0.533 1.135 -6.337 0.594 14.9 1.24 0.12 0.03 N
298.14820833 22.00248056 19.89 0.01 19.52 0.01 20.39 0.02 163.6 1.9 0.014 1.404 -4.599 0.742 14.3 1.82 0.23 -2.5 0.12 0.00 N
298.16396667 22.06349722 19.94 0.01 19.53 0.01 20.46 0.02 41.4 2.0 3.398 1.562 -8.363 0.947 14.0 2.04 0.25 -2.01 0.11 0.00 N
298.16217500 22.05441111 19.96 0.01 19.24 0.01 20.72 0.02 -176.0 1.5 -1.281 1.022 -1.531 0.539 15.8 2.65 0.16 -2.24 0.08 0.99 Y
298.14762500 22.18983611 20.06 0.01 19.67 0.01 20.63 0.02 -227.3 3.3 -1.419 1.191 -5.42 0.742 11.3 1.57 0.18 0.00 N
298.19400833 22.08638333 20.38 0.01 19.70 0.01 21.09 0.03 -174.9 2.3 -0.305 1.333 0.776 0.766 13.2 2.59 0.23 -2.09 0.11 0.99 Y
298.19723333 22.12452500 20.51 0.01 19.94 0.01 21.45 0.03 -323.9 2.6 -10.391 1.957 -3.304 1.107 12.0 0.01 N
298.14902500 22.02358611 20.91 0.01 20.49 0.01 21.45 0.03 -71.1 3.6 7.2 2.05 0.31 0.00 N
298.18112083 22.06077222 20.93 0.01 20.29 0.01 21.60 0.04 -179.5 2.4 9.7 2.54 0.28 1.00 Y
298.17317500 22.11583611 21.10 0.01 20.51 0.01 21.83 0.04 -173.9 3.8 6.2 2.06 0.22 0.97 B
298.14954583 22.10703056 21.24 0.01 20.62 0.01 21.89 0.05 -177.4 3.7 6.2 2.69 0.45 0.99 Y
298.18245417 22.10563889 21.27 0.01 20.67 0.01 21.89 0.05 -176.3 6.0 5.6 2.67 0.36 0.99 Y
298.15165833 22.14928333 21.54 0.01 21.13 0.02 22.07 0.05 10.6 10.6 4.6 1.82 0.51 0.00 N
298.16334583 22.14321667 21.60 0.01 21.16 0.02 22.27 0.06 71.0 4.9 4.0 1.31 0.36 0.00 N
298.16188333 22.05321389 21.56 0.01 20.99 0.01 22.31 0.07 -177.5 7.8 5.7 1.81 0.37 1.00 Y
298.18229167 22.04275278 21.61 0.01 20.99 0.01 22.22 0.06 -177.1 3.9 5.5 2.96 0.47 0.99 Y
298.19592917 22.13295000 21.65 0.01 21.03 0.01 22.88 0.10 22.0 6.6 -0.3 0.90 N
298.19628750 22.18246389 21.79 0.01 21.30 0.02 22.57 0.08 -285.7 12.8 4.2 1.33 0.37 0.02 N
298.13948750 22.17620000 21.81 0.01 21.41 0.02 22.37 0.07 174.5 11.6 2.6 1.49 0.43 0.00 N
298.15003333 22.01741667 21.77 0.01 21.01 0.01 22.59 0.08 -276.5 3.6 4.5 2.68 0.56 0.00 N
298.14941250 22.17936667 21.82 0.01 21.21 0.02 22.92 0.11 -59.4 5.6 4.1 0.57 N
298.16383333 22.18616667 21.90 0.01 21.11 0.02 23.39 0.16 -11.5 5.2 3.3 0.00 N
298.16461667 22.09165278 21.93 0.01 21.41 0.02 22.58 0.08 -172.1 4.5 3.7