Issue |
A&A
Volume 681, January 2024
|
|
---|---|---|
Article Number | L13 | |
Number of page(s) | 7 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202348554 | |
Published online | 15 January 2024 |
Letter to the Editor
New KiDS in town
Sextans II: A new stellar system on the outskirts of the Milky Way
1
INAF-Osservatorio Astronomico di Capodimonte, Via Moiariello, 16, 80131 Napoli, Italy
e-mail: massimiliano.gatto@inaf.it
2
INAF-Osservatorio di Astrofisica e Scienza dello Spazio, Via Gobetti, 93/3, 40129 Bologna, Italy
3
Physics Department, University of Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
4
INFN-Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
5
Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
6
Ruhr University Bochum, Faculty of Physics and Astronomy, Astronomical Institute (AIRUB), German Centre for Cosmological Lensing, 44780 Bochum, Germany
7
Kapteyn Astronomical Institute, University of Groningen, PO Box 800, 9700 AV Groningen, The Netherlands
8
School of Physics and Astronomy, Sun Yat-sen University, Guangzhou, 519082 Zhuhai Campus, PR China
9
Department of Physics and Astronomy, University of Victoria, Victoria, BC V8P 1A1, Canada
Received:
10
November
2023
Accepted:
22
December
2023
We report the discovery of a significant and compact over-density of old and metal-poor stars in the fourth data release of the KiDS survey (DR4). The discovery is confirmed by deeper HSC-SSP data revealing the old main sequence turn-off of a stellar system located at a distance from the sun of D⊙ = 145-13+14 kpc in the direction of the Sextans constellation. The system has an absolute integrated magnitude (MV = −3.9 ± 0.4), half-light radius (rh = 193-46+61 pc), and ellipticity (e = 0.46-0.15+0.11) values that are typical of ultra-faint dwarf galaxies (UFDs). The central surface brightness is near the lower limits of known local dwarf galaxies of a similar integrated luminosity, as expected for stellar systems that have escaped detection until now. The distance of the newly found system suggests that it is likely to be a satellite of our own Milky Way and we have thus tentatively named it Sextans II (KiDS-UFD-1).
Key words: surveys / Hertzsprung-Russell and C-M diagrams / Galaxy: halo / galaxies: clusters: general / galaxies: dwarf / galaxies: photometry
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
The past two decades have seen a surge of interest in the research of Milky Way (MW) satellites thanks to the advent of digital wide-field deep panchromatic photometric surveys. The application of high-performance over-density detection techniques to large data sets has led to a substantial increase in the number of known stellar systems inhabiting the MW halo (e.g. Belokurov et al. 2007; Walsh et al. 2009; Koposov et al. 2015; Torrealba et al. 2016; Cerny et al. 2021, and references therein). The quest for MW satellite galaxies carried out with modern facilities has opened up the study of the faint end of the galaxy luminosity function to so-called ultra-faint dwarf (UFD) galaxies (see Belokurov 2013, for discussion and references). These objects provide valuable insights into the mass assembly history of the MW. In addition, the UFDs are the most dark matter-dominated objects in the Universe, as well as the faintest, oldest, and least chemically evolved galaxies known so far (see e.g. the detailed review by Simon 2019, and references therein). Therefore, they represent fossil evidence that is useful for probing the very early stages of the Universe, such as the epoch of the reionisation (e.g. Bovill & Ricotti 2011; Wheeler et al. 2015). Because of their pristine nature, UFDs are also ideal systems for testing and improving our understanding of the sub-branches of the stellar evolutionary models, such as the synthesis of heavy elements (e.g. Ji et al. 2016).
The majority of faint galaxies discovered in recent years were identified through a few large sky surveys that offer unprecedented photometric depth, including Sloan Digital Sky Survey (SDSS; Willman et al. 2005a,b; Belokurov et al. 2006, 2007, 2009; Belokurov et al. 2010; Zucker et al. 2006b,a), Dark Energy Survey (DES; Bechtol et al. 2015; Drlica-Wagner et al. 2015; Kim & Jerjen 2015; Koposov et al. 2015; Luque et al. 2016), PanSTARRS (Laevens et al. 2015a,b), ATLAS (Torrealba et al. 2016), DELVE (Mau et al. 2020; Cerny et al. 2021, 2023), DESI (Martínez-Delgado et al. 2022; Sand et al. 2022), MagLiteS (Torrealba et al. 2018), HSC-SSP (Homma et al. 2016, 2018, 2019), and UNIONS (Smith et al. 2023a,b). In this study, we expand the search for unknown stellar systems to Kilo-Degree Survey (KiDS; de Jong et al. 2013), which has not yet been exploited for this purpose.
KiDS is an extensive multi-band photometric survey conducted using the VLT Survey Telescope (VST; Capaccioli & Schipani 2011; Capaccioli et al. 2012). It has observed about 1350 square degrees of the sky in the u, g, r, and i bands. KiDS was conceived as a cosmological survey, mainly with an aim to use weak gravitational lensing of galaxies to study the assembly of large-scale structures. However, as a side-product of extra-galactic studies, KiDS provides a deep and accurate catalogue of stars, making it suitable in the search for low-surface brightness stellar systems in the Local Group. In this Letter, we present a new UFD galaxy uncovered in a large-scale search for low-luminosity stellar over-densities over the whole KiDS Data Release 4 (DR4) area.
2. KiDS data
We used the most recent KiDS data release (DR4, Kuijken et al. 2019), which includes 1006 of imaging and catalogue data with limiting 5σ AB magnitudes of 24.23 ± 0.12, 25.12 ± 0.14, 25.02 ± 0.13, and 23.68 ± 0.27 mag in the u, g, r, and i filters, respectively. These data are ≳2 mag deeper than SDSS and PanSTARRS DR1. The data also include photometry in the Z, Y, J, H, Ks infrared (IR) bands obtained with the VISTA Kilo-degree INfrared Galaxy survey (VIKING; Edge et al. 2013) carried out on the VISTA telescope. KiDS covers two different stripes in the sky, namely, a northern and a southern stripe (KiDS-N and KiDS-S, hereafter; see Table 2 in Kuijken et al. 2019, for the exact boundaries).
From the KiDS DR4 catalogue, we selected sources identified as stars, namely objects that are flagged as 1, 4, or 5 in the SG2DPHOT column (see Appendix A.1.2 in Kuijken et al. 2019). In this way, we obtained a list of 4 422 730 stars in KiDS-N and 3 116 050 stars in KiDS-S1. Successively, we filtered out sources with low-quality photometry by imposing a maximum tolerated photometric uncertainty in both the g and r bands of 0.2 mag. Nonetheless, the catalogues at our disposal still include distant galaxies that have been erroneously classified as stars. To enhance the purity of our catalogues, we took advantage of the multi-band photometry provided by KiDS, utilising colour-colour diagrams to further purge unresolved galaxies. Specifically, we focused on the g − r versus g − i plot, where stars occupy a distinctive, narrow sequence (see for example Fig. 18 in Kuijken et al. 2019). To this end, we selected only sources falling within the boundaries of a polygon defined by the following coordinates (expressed in mag): (g − r,r − i) = (1.25, 0.25);(1.25, 0.8);(−0.65, −0.05);(−0.65, −0.6). It is important to note that we also opted to exclude red stars, characterised by g − r > 1.25 mag, as this colour range is dominated by M-type dwarf stars located within the MW disk. After the colour-colour selection, we are left with 3 558 843 stars, 2 127 949 of which are located within the footprint of KiDS-N, and the remaining 1 430 894 stars are in KiDS-S.
3. Method
We tackled the task of unveiling unknown faint MW satellites by identifying the over-densities of stars in the sky relative to the local average stellar background. To achieve this, we utilised the algorithm and the procedure described in Gatto et al. (2020). In addition, to augment the chances of detecting faint stellar systems projected against the dominant MW stellar population, we also followed the approach of Walsh et al. (2009), which has been widely used in recent years for the discovery of low-luminosity galaxies and globular clusters in the Milky Way halo (e.g. Bechtol et al. 2015; Koposov et al. 2015; Torrealba et al. 2016; Homma et al. 2019; Cerny et al. 2021, and references therein). This method consists of a preliminary isochrone-based matched filter on the catalogue of stars. In particular, prior to running the searching algorithm, we selected stars whose positions in the g − r, r colour-magnitude diagram (CMD) were within 0.05 mag in colour from a given isochrone representative of an old and metal-poor stellar population2, which is characteristic of the faint galaxies expected to inhabit the MW halo. We opted to work with the suite of PARSEC isochrones (Bressan et al. 2012) and the isochrone filter was carried out by adopting three different ages, namely, log(t/yr) = 10.0, 10.06, 10.12 dex, and three different metallicities, or [Fe/H] ∼ −2.2, −1.5, −1.2 dex. Additionally, we varied the distance modulus of the isochrone in a range between 16 and 23.6 mag (corresponding to a range of physical distances between ∼16 and ∼520 kpc) with steps of 0.4 mag. Therefore, we scanned the KiDS catalogue with 180 different isochrone filters.
For each isochrone, we analysed the resulting filtered map as follows. We divided the map into 4° ×4° regions, and then build density maps using the star coordinates as input parameters, moothed through a kernel density estimation (KDE) technique adopting a Gaussian kernel. Subsequently, the algorithm counts the number of stars (Nstars) in each pixel of the sub-field and measures the number of standard deviations above a local mean:
where Nbkg and σbkg represent the average of the local stellar background and its standard deviation, respectively3. In this work, we use a pixel size of 30″× 30″and a bandwidth of the Gaussian kernel function of 1′4. As a sanity check, we conducted a validation of the algorithm by assessing its capability to accurately identify known stellar systems residing within the KiDS footprint, including Galactic globular clusters and dwarf galaxies. In this way, our algorithm recovered Leo V with Spxl = 21.5, by applying an isochrone filter5 with log t = 10.06 dex, [Fe/H] ∼ −2.2, and a distance modulus = 21.2 mag, corresponding to a distance of approximately 175 kpc. This determination is closely aligned with the previously estimated distance to the dwarf, reported as 173 ± 5 kpc by Medina et al. (2017).
4. Results
Utilising the methodology outlined in the previous section, we identified a highly promising over-density in the Sextans constellation. Our detection revealed a Spxl = 10.6 stellar over-density by applying an isochrone filter with log(t/yr) = 10.12 dex, [Fe/H] ∼ −1.5, and a distance modulus (m − M)0 = 20.4 magnitudes. The left panel of Fig. 1 shows the density map of stars in the region of the newfound galaxy candidate, for the aforementioned isochrone selection filter. The over-density stands out clearly. The subsequent three panels of the figure illustrate the CMDs of stars within a 5′ radius from the centre of the over-density, the CMD of a representative local field covering an equivalent area, and the Hess diagram derived through the subtraction of the over-density and field CMDs, from left to right, respectively. The inspection of these CMDs unveils a well-sampled red giant branch (RGB) within the inner 5′, that is absent in the local field CMD, suggesting the possible detection of a stellar system. Three likely candidate horizontal branch (HB) stars, highlighted as cyan squares in the figure, are found at (g − r, r)≃(0, 21 − 21.5) mag6. Unfortunately, the photometric depth of the KiDS data is not sufficient to reach the main sequence turn-off (MSTO) of the candidate7, whose detection would be crucial in order to confirm it as a genuine stellar system and to better assess its properties. In the next section, we describe the successful analysis of the candidate with deeper photometry.
![]() |
Fig. 1. Density map of a region with radius 1° centered on the galaxy candidate. The density map is obtained using stars from the KiDS catalog, and filtered by an isochrone with log(t/yr) = 10.12 dex, [Fe/H] ∼ −1.5 and a distance modulus = 20.4 mag (first panel). We adopted a Gaussian kernel with bandwidth = 2.0′ to smooth the figure. CMD of stars located within 5′ from the over-density center. A red isochrone with log(t/yr) = 10.12 dex; [Fe/H] ∼ −1.5 dex, and m − M = 20.4 mag is superimposed on the figure (second panel). Cyan squares indicate the position of stars compatible with the HB evolutionary stage. CMD of stars within a local field having the same area covered by the CMD built for the dwarf galaxy candidate (third panel). Hess diagram obtained by subtracting the CMDs displayed in the second and third panels (fourth panel). |
4.1. Deeper photometry in HSC-SSP
The Hyper Suprime Cam Subaru Strategic Program (HSC-SSP; Aihara et al. 2018a,b) is a wide-field optical-near infrared survey conducted using the 8.2-meter Subaru Telescope. The HSC-SSP is designed to cover approximately 1400 square degrees of the Northern sky across five distinct filters, namely, g, r, i, z, y. The survey 5σ magnitude limits for point-sources are impressively deep: 26.5 mag in the g band, 26.1 mag in the r band, and 25.9 mag in the i band (Aihara et al. 2018b), surpassing the limiting magnitude of the KiDS survey by more than one mag in each passband. HSC-SSP has already demonstrated its power for the detection of low-surface-brightness stellar systems (Homma et al. 2016, 2018, 2019). In particular, Homma et al. (2019) reported the discovery of Bootes IV, one of the faintest UFDs identified to date, characterised by surface brightness levels of μ ≥ 32 mag arcsec−2. Since our candidate lies within the HSC-SSP footprint, this survey provides an excellent means to follow it up.
From the HSC-SSP catalogue provided by the public third data release (Aihara et al. 2022), which included data for slightly less than 700 square degrees, we selected point-like sources within the region of the newly discovered over-density by means of the ’extendedness’ parameter in i-band, which is the filter with the best median seeing (i.e. ∼0.6″, Aihara et al. 2022). To further refine the star-galaxy separation, we used a colour-colour diagram, as done above for the KiDS data. In this case, we adopted (g − r, z − y) diagram that is particularly efficient in this sense. Consequently, we selected all point-like sources located within the polygon defined by the criteria: (g − r, z − y) = (−0.3, −0.12);(−0.3, 0.08);(1.2, 0.25);(1.2, 0). The following sections of this manuscript are based on the analysis of this catalogue.
4.2. Analysis of the CMD
In Fig. 2, we show the results of the analysis of the deeper HSC-SSP catalogue. The top-left panel shows the density map of the stars filtered by the best-fit isochrone (indicated by a red solid line in the bottom row of the same figure), fully confirming the detection of a strong and compact stellar over-density at the expected location. In contrast, the right panel of the same figure illustrates the density map for sources catalogued as galaxies based on the ‘extendedness’ parameter and, subsequently, filtered by the same isochrone. Notably, no discernible over-density of galaxies in the same region is observed.
![]() |
Fig. 2. Top row: density map of stars included in the HSC-SSP catalog, and filtered by an isochrone characterised by log(t/yr) = 10.12 dex; [Fe/H] ∼ −2.2 dex, and m − M = 20.8 mag (left). The depicted region encompasses 80′ in radius centered on the dwarf galaxy. The figure has been smoothed with a Gaussian kernel with bandwidth = 2′. Density map of sources included in the HSC-SSP catalog and classified as galaxies, filtered by the same isochrone parameters adopted for the stars (right). Bottom row: CMD of the stars located within 2 × rh from the galaxy center (panel a). A red (green) isochrone with log(t/yr) = 10.00 (log(t/yr) = 10.12); [Fe/H] ∼ −2.2 dex, and m − M = 20.8 mag is overlaid to the figure. CMD of stars within a local field of the same area as that used to construct the CMD in the panel to the left (panel b). Hess diagram obtained by subtracting the CMDs displayed in panels a and b (panel c). CMD of the stars located within 4 × rh from the galaxy center (panel d). An inset positioned at the top of the panel offers a zoomed-in view of the HB stage, capturing a narrow colour interval and highlighting the presence of four closely-aligned HB stars. |
Panel a of the second row of Fig. 2 displays the CMD8 of all stars located within twice the half-light radius of the structure (for the determination of rh, please refer to Sect. 4.3). The sub-giant branch (SGB) and the MSTO associated with the RGB already detected in KiDS data are unequivocally revealed. Notably, these features are not present in the CMD of the control field (same area, located 20′ from the centre of the over-density; see also the Hess diagram obtained by subtracting the two former CMDs). This compelling evidence strongly suggests that our candidate is a genuine stellar system, likely a faint local dwarf galaxy, which has hitherto gone unreported in the existing literature. In what follows, we refer to it as Sextans II (KiDS-UFD-1)9.
An in-depth analysis of the CMD also discloses the presence of four stars likely belonging to the HB evolutionary phase. Specifically, two stars with colours of −0.06 mag and −0.08 mag reside within the half-light radius of the galaxy. The third star, slightly exceeding the rh, exhibits g − r = −0.09 mag, while the fourth one exhibits g − r = 0.03 and is located at ∼8.3′ from the centre of the system (≃1.8rh, see Sect. 4.3). The absolute magnitude of the HB is a valuable distance indicator and we used it as the strongest constraint on the distance modulus. Through a visual isochrone fitting, we re-determined the parameters that provide the best match to the observed CMD: log(t/yr)≃10.12 dex; [Fe/H] ≃ −2.2 dex; and m − M = 20.8 ± 0.2 mag. This is the red solid line superimposed on all panels in the second row of Fig. 2.
The resulting distance modulus corresponds to a physical distance of approximately 150 kpc, well within the range of Galactic satellites. The green isochrone overlaid on panel a in the second row of Fig. 2 deviates from the best-matching isochrone only in terms of age. Specifically, it is characterised by log(t/yr) = 10.00 dex, showing that decreasing the age of the stellar model by ∼ ± 2 − 3 Gyr does not significantly affect the absolute magnitude of the HB (hence, our distance estimate), as expected.
Finally, panel d in the second row of Fig. 2 presents the CMD of Sextans II, including all stars confined within 4 × rh. This larger region enables a more comprehensive search for additional HB stars that might be associated with the system, albeit at the cost of higher contamination from MW stars within the CMD. We identified a potential fifth HB member at ∼15.9′ from the galaxy centre, corresponding to 3.4 × rh. This star has a colour of −0.09 mag, similar to the other three HB stars, and indeed it is barely discernible in the main panel. To enhance the visibility of these four closely colour-matched HB stars, we provide a zoomed-in view of the HB evolutionary phase in an inset at the top of panel d, confined to a narrow range of colours encompassing those of the four stars.
4.3. Structural parameters
To estimate the structural parameters, such as the centre, half-light radius, ellipticity, and position angle of the newfound stellar system, we adopted the maximum likelihood approach presented in Martin et al. (2008). In particular, Martin et al. (2008) model the two-dimensional (2D) position of the stars in the sky in the surroundings of the galaxy with the combination of the exponential profile of the galaxy and of a uniform background.
We adopted a Markov chain Monte Carlo (MCMC) technique to contemporaneously estimate the free parameters of our model: x0; y0; rh; e; θ; Nstars, employing only the stars located within the best-fit isochrone mask. Our MCMC implementation utilised 1000 different walkers, each executing a chain composed of 10000 steps. We discarded the first 1000 steps as burn-in, resulting in a robust estimation of the model parameters. The MCMC analysis was conducted using the emcee (Foreman-Mackey et al. 2013) Python library. The best-fit parameters and their uncertainties were determined from the median, 16th, and 84th quantiles of the posterior distributions. In Fig. 3, we present a corner plot displaying the posterior distribution of the six free parameters, while Table 1 lists the median values and uncertainties of these parameters.
![]() |
Fig. 3. Posterior probability distributions for each of the six free parameters of the model, derived from the MCMC approach described in Sect. 4.3. |
Properties of Sextans II (KiDS-UFD-1).
The left panel of Fig. 4 illustrates the radial density profile (RDP) of the system, where r represents the elliptical radius, as defined in Eq. (4) in Martin et al. (2008). The red solid line denotes the best fit of the model for the 2D position of the stars in the sky derived by employing the median values of the six parameters obtained through MCMC analysis. It is interesting to note that, both the derived half-light radius, pc (assuming D⊙ = 145 kpc), and ellipticity
are well within the range typical for local UFDs (see the latest version of the McConnachie 2012, catalogue). The right panel of Fig. 4 displays the 2D relative position of the filtered stars with respect to the galaxy centre, along with ellipses aligned according to the position angle determined by the MCMC procedure. These ellipses are defined by semi-major axes of 2rh (dashed ellipse) and 4rh (solid ellipse). Additionally, we indicate the relative position of the five stars in the HB evolutionary stage as cyan diamonds. As already stated, four HB stars are situated within 2rh, while the fifth one is at almost 4rh.
![]() |
Fig. 4. Radial density profile of KiDS-UFD-1 (left). Here, rell is the elliptical radius defined in Eq. (4) in Martin et al. (2008). The red solid line refers to the best-fit model with the parameters derived via the MCMC approach described in Sect. 4.3, while the gray lines represent 1000 random extractions from the posterior distributions of the free parameters of the model. Relative position of all stars filtered by the best-fit isochrone with respect to the galaxy center (right). The red dashed and solid ellipses have ellipticity and position angle derived from the MCMC procedure and semi-major axes of 2 × rh and 4 × rh, respectively. Cyan diamonds indicate the relative position of the five HB stars. |
4.4. Integrated magnitude and stellar mass
In this section, we outline the procedure employed to estimate the absolute magnitude in the V-band for the purpose of comparing the properties of Sextans II with other recently discovered UFDs. Our approach follows the methodology reported in Martin et al. (2008), which is based on Monte Carlo sampling of a synthetic CMD representative of Sextans II. Specifically, we built a synthetic stellar population with log(t/yr) = 10.12 dex; [Fe/H] ∼ −2.2 dex, placed at a distance modulus 20.8 mag. We considered a total of Nstars = 52 ± 12, that is the value obtained from the MCMC approach by adopting a cut at r = 24 mag, which is a conservative threshold adopted to mitigate potential issues of completeness in the central regions of the galaxy that might bias the estimation of its absolute magnitude. We generated 1000 realisations of the CMD, populating it with a Kroupa (2002) initial mass function (IMF), and subsequently computed the average absolute magnitude in the g and r SDSS filters, as well as in the Johnson V band10. Additionally, we determined the stellar mass of the galaxy. The median values, along with the 16th and 84th quantiles, were adopted to represent the best values and associated uncertainties, respectively. These results are presented in the final rows of Table 1.
Figure 5 depicts the MV–rh plane for galaxies and candidate galaxies reported in catalog by McConnachie (2012; version 2021), as well as for galactic globular clusters included in the catalog by Harris (1996; version 2010). The position of Sextans II within this plot clearly establishes it as one of the UFDs with the lowest surface brightness discovered to date, although it is important to note that some extremely faint low-surface brightness galaxies identified in recent years may not be represented in this figure.
![]() |
Fig. 5. Absolute integrated V-band magnitude as a function of the log of the half-light radius for local low-luminosity stellar systems. Galactic globular clusters, represented as open grey triangles, are from the 2010 version of the Harris (1996) catalogue. Confirmed and candidate dwarf galaxies, represented as filled circles colour-coded according to their heliocentric distance, are taken from the 2021 version of the catalog by McConnachie (2012). The same colour-coding is applied to the filled square marking the position of Sextans II, that is also highlighted with a red contour. The stellar systems with MV and rh values that are most similar to Sextans II have been labelled. The diagonal lines are loci of constant surface brightness, from at μ = 24.0 mag arcsec−2 (top) down to μ = 34.0 mag arcsec−2 (bottom), in steps of μ = 2.0 mag arcsec−2. |
5. Summary and conclusions
In this Letter we report on the discovery of an old and metal-poor stellar system located in the outskirts of the Milky Way, Sextans II (KiDS-UFD-1). This discovery resulted from an extensive search for low-surface brightness stellar systems in the MW halo within the KiDS survey catalogues. The KiDS detection, limited by photometric depth to the RGB and HB of the system, was confirmed with deeper HSC-SSP photometry which clearly unveiled the old MS of Sextans II down to about 1 mag below the MSTO. The isochrone fitting of the CMD revealed that the newly discovered system is ancient (t ≥ 10 Gyr) and metal-poor ([Fe/H] ≤ −1.5) and is located at a distance of D⊙ ≃ 145 kpc from the Sun and RGC ≃ 128.8 kpc from the center of the Galaxy. The structural parameters of Sextans II are typical of UFDs in the Local Group, strongly suggesting that the new system is a faint spheroidal satellite of the MW at the low end of the surface brightness distribution of local dwarfs. The most plausible alternative hypothesis is a disrupting globular cluster. However, even the most extended known globular clusters in the MW have rh ≤ 30.0 pc, more than five times smaller than what we determined for Sextans II; hence the UDF hypothesis is, by far, the most likely. The final word on the nature of the system can be provided only by a proper spectroscopic follow up of a reasonable sample of member stars, which may be challenging, given the magnitude range spanned by candidate RGB members11. The chemical composition can reveal the presence of a metallicity spread, while the kinematics can reveal the amount of non-baryonic dark matter required to keep the system at the dynamical equilibrium, both features that are generally used to discriminate between dwarf galaxies and star clusters (Willman & Strader 2012).
The research of unknown stellar systems is still ongoing, and Sextans II might not be the only newfound stellar system detected in the KiDS DR4 catalog. Moreover, this work underscores the importance of extending this research to the forthcoming fifth release of the KiDS survey, that will provide multi-band photometric data covering the entire 1350 square degrees of the KiDS footprint. This discovery serves as a reminder that the census of the Milky Way’s satellite galaxy population is far from complete. The upcoming wide and deep surveys, such as the Legacy Survey of Space and Time (LSST) set to be conducted at the imminent Vera Rubin Observatory, will likely contribute significantly to expanding our knowledge of Milky Way satellites, much like the impactful contributions of the SDSS and PanSTARRS surveys before it.
Magnitudes provided by KiDS are corrected for interstellar extinction using the Schlegel et al. (1998) reddening maps re-calibrated after Schlafly & Finkbeiner (2011), in combination with the Fitzpatrick’s law (Fitzpatrick 1999).
Please refer to Gatto et al. (2020), in particular their Sect. 3.1, for a detailed description of the detection algorithm and definitions of all the involved quantities.
Although the 5σ AB magnitudes reach r ≃ 25 mag, the cuts applied to the catalogue, described in Sect. 2, limited our analysis to r ≃ 23.5 − 24 mag in the CMD, preventing us from probing fainter magnitudes.
It is important to highlight that the HSC filter system is slightly different from the SDSS filter system. To work with the SDSS filter system available in the PARSEC library, we adopted the correction reported in Homma et al. (2016).
Concurrently and independently to our research, Homma et al. (2023) also reported the discovery of Sextans II by using the HSC-SSP catalogue.
To convert the g SDSS filter magnitude into the Johnson V band magnitude we applied the conversion relationship described in Jester et al. (2005).
In the recently published FASHI catalogue of HI sources (Zhang et al. 2024) we found a source at cz = 10 443.8 km s−1 (FASHI 20230004672, estimated distance ∼140 Mpc) that coincides in the sky with Sextans II, has a characteristic size (3.9′) comparable with the rh of Sextans II, and a mass compatible with a galaxy group. While we are confident that our colour cuts get rid of most of the contamination by unresolved galaxies, only spectroscopic follow-up will ascertain if there is significant contamination of Sextans II from galaxies possibly associated with this distant HI source.
Acknowledgments
We warmly thank the anonymous referee for the suggestions that helped us to improve the manuscript. Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme IDs 177.A-3016, 177.A-3017, 177.A-3018 and 179.A-2004, and on data products produced by the KiDS consortium. The KiDS production team acknowledges support from: Deutsche Forschungsgemeinschaft, ERC, NOVA and NWO-M grants; Target; the University of Padova, and the University Federico II (Naples). The Hyper Suprime-Cam (HSC) collaboration includes the astronomical communities of Japan and Taiwan, and Princeton University. The HSC instrumentation and software were developed by the National Astronomical Observatory of Japan (NAOJ), the Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), the University of Tokyo, the High Energy Accelerator Research Organization (KEK), the Academia Sinica Institute for Astronomy and Astrophysics in Taiwan (ASIAA), and Princeton University. Funding was contributed by the FIRST program from the Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. This paper is based on data collected at the Subaru Telescope and retrieved from the HSC data archive system, which is operated by the Subaru Telescope and Astronomy Data Center (ADC) at NAOJ. Data analysis was in part carried out with the cooperation of Center for Computational Astrophysics (CfCA), NAOJ. We are honored and grateful for the opportunity of observing the Universe from Maunakea, which has the cultural, historical and natural significance in Hawaii. M.G. acknowledges the INAF AstroFIt grant 1.05.11. C.T. acknowledges the INAF grant 2022 LEMON. JTAdJ is supported by the Netherlands Organisation for Scientific Research (NWO) through grant 621.016.402. H. Hildebrandt is supported by a DFG Heisenberg grant (Hi 1495/5-1), the DFG Collaborative Research Center SFB1491, as well as an ERC Consolidator Grant (No. 770935). K.K. acknowledges support from the Royal Society and Imperial College. S.Z. acknowledges the support from the Deutsche Forschungsgemeinschaft (DFG) SFB1491. N.R.N. acknowledges financial support from the National Science Foundation of China, Research Fund for Excellent International Scholars (grant n. 12150710511), and from the research grant from China Manned Space Project n. CMS-CSST-2021-A01.
References
- Aihara, H., Arimoto, N., Armstrong, R., et al. 2018a, PASJ, 70, S4 [NASA ADS] [Google Scholar]
- Aihara, H., Armstrong, R., Bickerton, S., et al. 2018b, PASJ, 70, S8 [NASA ADS] [Google Scholar]
- Aihara, H., AlSayyad, Y., Ando, M., et al. 2022, PASJ, 74, 247 [NASA ADS] [CrossRef] [Google Scholar]
- Bechtol, K., Drlica-Wagner, A., Balbinot, E., et al. 2015, ApJ, 807, 50 [NASA ADS] [CrossRef] [Google Scholar]
- Belokurov, V. 2013, New Astron. Rev., 57, 100 [CrossRef] [Google Scholar]
- Belokurov, V., Zucker, D. B., Evans, N. W., et al. 2006, ApJ, 647, L111 [NASA ADS] [CrossRef] [Google Scholar]
- Belokurov, V., Zucker, D. B., Evans, N. W., et al. 2007, ApJ, 654, 897 [Google Scholar]
- Belokurov, V., Walker, M. G., Evans, N. W., et al. 2009, MNRAS, 397, 1748 [CrossRef] [Google Scholar]
- Belokurov, V., Walker, M. G., Evans, N. W., et al. 2010, ApJ, 712, L103 [NASA ADS] [CrossRef] [Google Scholar]
- Bovill, M. S., & Ricotti, M. 2011, ApJ, 741, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [Google Scholar]
- Capaccioli, M., & Schipani, P. 2011, The Messenger, 146, 2 [NASA ADS] [Google Scholar]
- Capaccioli, M., Schipani, P., de Paris, G., et al. 2012, Science from the Next Generation Imaging and Spectroscopic Surveys, 1 [Google Scholar]
- Cerny, W., Pace, A. B., Drlica-Wagner, A., et al. 2021, ApJ, 910, 18 [CrossRef] [Google Scholar]
- Cerny, W., Martínez-Vázquez, C. E., Drlica-Wagner, A., et al. 2023, ApJ, 953, 1 [NASA ADS] [CrossRef] [Google Scholar]
- de Jong, J. T. A., Verdoes Kleijn, G. A., Kuijken, K. H., & Valentijn, E. A. 2013, Exp. Astron., 35, 25 [Google Scholar]
- Drlica-Wagner, A., Bechtol, K., Rykoff, E. S., et al. 2015, ApJ, 813, 109 [Google Scholar]
- Edge, A., Sutherland, W., Kuijken, K., et al. 2013, The Messenger, 154, 32 [NASA ADS] [Google Scholar]
- Fitzpatrick, E. L. 1999, PASP, 111, 63 [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
- Gatto, M., Ripepi, V., Bellazzini, M., et al. 2020, MNRAS, 499, 4114 [Google Scholar]
- Harris, W. E. 1996, AJ, 112, 1487 [Google Scholar]
- Homma, D., Chiba, M., Okamoto, S., et al. 2016, ApJ, 832, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Homma, D., Chiba, M., Okamoto, S., et al. 2018, PASJ, 70, S18 [NASA ADS] [CrossRef] [Google Scholar]
- Homma, D., Chiba, M., Komiyama, Y., et al. 2019, PASJ, 71, 94 [NASA ADS] [CrossRef] [Google Scholar]
- Homma, D., Chiba, M., Komiyama, Y., et al. 2023, PASJ, submitted, [arXiv:2311.05439] [Google Scholar]
- Jester, S., Schneider, D. P., Richards, G. T., et al. 2005, AJ, 130, 873 [Google Scholar]
- Ji, A. P., Frebel, A., Chiti, A., & Simon, J. D. 2016, Nature, 531, 610 [NASA ADS] [CrossRef] [Google Scholar]
- Kim, D., & Jerjen, H. 2015, ApJ, 808, L39 [NASA ADS] [CrossRef] [Google Scholar]
- Koposov, S. E., Belokurov, V., Torrealba, G., & Evans, N. W. 2015, ApJ, 805, 130 [NASA ADS] [CrossRef] [Google Scholar]
- Kroupa, P. 2002, Science, 295, 82 [Google Scholar]
- Kuijken, K., Heymans, C., Dvornik, A., et al. 2019, A&A, 625, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Laevens, B. P. M., Martin, N. F., Bernard, E. J., et al. 2015a, ApJ, 813, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Laevens, B. P. M., Martin, N. F., Ibata, R. A., et al. 2015b, ApJ, 802, L18 [NASA ADS] [CrossRef] [Google Scholar]
- Luque, E., Queiroz, A., Santiago, B., et al. 2016, MNRAS, 458, 603 [NASA ADS] [CrossRef] [Google Scholar]
- Martin, N. F., de Jong, J. T. A., & Rix, H.-W. 2008, ApJ, 684, 1075 [NASA ADS] [CrossRef] [Google Scholar]
- Martínez-Delgado, D., Karim, N., Charles, E. J. E., et al. 2022, MNRAS, 509, 16 [Google Scholar]
- Mau, S., Cerny, W., Pace, A. B., et al. 2020, ApJ, 890, 136 [NASA ADS] [CrossRef] [Google Scholar]
- McConnachie, A. W. 2012, AJ, 144, 4 [Google Scholar]
- Medina, G. E., Muñoz, R. R., Vivas, A. K., et al. 2017, ApJ, 845, L10 [NASA ADS] [CrossRef] [Google Scholar]
- Sand, D. J., Mutlu-Pakdil, B., Jones, M. G., et al. 2022, ApJ, 935, L17 [CrossRef] [Google Scholar]
- Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [Google Scholar]
- Simon, J. D. 2019, ARA&A, 57, 375 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, S. E. T., Cerny, W., Hayes, C. R., et al. 2023a, ApJ, submitted, [arXiv:2311.10147] [Google Scholar]
- Smith, S. E. T., Jensen, J., Roediger, J., et al. 2023b, AJ, 166, 76 [CrossRef] [Google Scholar]
- Torrealba, G., Koposov, S. E., Belokurov, V., & Irwin, M. 2016, MNRAS, 459, 2370 [NASA ADS] [CrossRef] [Google Scholar]
- Torrealba, G., Belokurov, V., Koposov, S. E., et al. 2018, MNRAS, 475, 5085 [NASA ADS] [CrossRef] [Google Scholar]
- Walsh, S. M., Willman, B., & Jerjen, H. 2009, AJ, 137, 450 [NASA ADS] [CrossRef] [Google Scholar]
- Wheeler, C., Oñorbe, J., Bullock, J. S., et al. 2015, MNRAS, 453, 1305 [NASA ADS] [CrossRef] [Google Scholar]
- Willman, B., & Strader, J. 2012, AJ, 144, 76 [Google Scholar]
- Willman, B., Blanton, M. R., West, A. A., et al. 2005a, AJ, 129, 2692 [NASA ADS] [CrossRef] [Google Scholar]
- Willman, B., Dalcanton, J. J., Martinez-Delgado, D., et al. 2005b, ApJ, 626, L85 [Google Scholar]
- Zhang, C.-P., Zhu, M., Jiang, P., et al. 2024, Sci. China Phys. Mech. Astron., 67, 219511 [NASA ADS] [CrossRef] [Google Scholar]
- Zucker, D. B., Belokurov, V., Evans, N. W., et al. 2006a, ApJ, 650, L41 [Google Scholar]
- Zucker, D. B., Belokurov, V., Evans, N. W., et al. 2006b, ApJ, 643, L103 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
![]() |
Fig. 1. Density map of a region with radius 1° centered on the galaxy candidate. The density map is obtained using stars from the KiDS catalog, and filtered by an isochrone with log(t/yr) = 10.12 dex, [Fe/H] ∼ −1.5 and a distance modulus = 20.4 mag (first panel). We adopted a Gaussian kernel with bandwidth = 2.0′ to smooth the figure. CMD of stars located within 5′ from the over-density center. A red isochrone with log(t/yr) = 10.12 dex; [Fe/H] ∼ −1.5 dex, and m − M = 20.4 mag is superimposed on the figure (second panel). Cyan squares indicate the position of stars compatible with the HB evolutionary stage. CMD of stars within a local field having the same area covered by the CMD built for the dwarf galaxy candidate (third panel). Hess diagram obtained by subtracting the CMDs displayed in the second and third panels (fourth panel). |
In the text |
![]() |
Fig. 2. Top row: density map of stars included in the HSC-SSP catalog, and filtered by an isochrone characterised by log(t/yr) = 10.12 dex; [Fe/H] ∼ −2.2 dex, and m − M = 20.8 mag (left). The depicted region encompasses 80′ in radius centered on the dwarf galaxy. The figure has been smoothed with a Gaussian kernel with bandwidth = 2′. Density map of sources included in the HSC-SSP catalog and classified as galaxies, filtered by the same isochrone parameters adopted for the stars (right). Bottom row: CMD of the stars located within 2 × rh from the galaxy center (panel a). A red (green) isochrone with log(t/yr) = 10.00 (log(t/yr) = 10.12); [Fe/H] ∼ −2.2 dex, and m − M = 20.8 mag is overlaid to the figure. CMD of stars within a local field of the same area as that used to construct the CMD in the panel to the left (panel b). Hess diagram obtained by subtracting the CMDs displayed in panels a and b (panel c). CMD of the stars located within 4 × rh from the galaxy center (panel d). An inset positioned at the top of the panel offers a zoomed-in view of the HB stage, capturing a narrow colour interval and highlighting the presence of four closely-aligned HB stars. |
In the text |
![]() |
Fig. 3. Posterior probability distributions for each of the six free parameters of the model, derived from the MCMC approach described in Sect. 4.3. |
In the text |
![]() |
Fig. 4. Radial density profile of KiDS-UFD-1 (left). Here, rell is the elliptical radius defined in Eq. (4) in Martin et al. (2008). The red solid line refers to the best-fit model with the parameters derived via the MCMC approach described in Sect. 4.3, while the gray lines represent 1000 random extractions from the posterior distributions of the free parameters of the model. Relative position of all stars filtered by the best-fit isochrone with respect to the galaxy center (right). The red dashed and solid ellipses have ellipticity and position angle derived from the MCMC procedure and semi-major axes of 2 × rh and 4 × rh, respectively. Cyan diamonds indicate the relative position of the five HB stars. |
In the text |
![]() |
Fig. 5. Absolute integrated V-band magnitude as a function of the log of the half-light radius for local low-luminosity stellar systems. Galactic globular clusters, represented as open grey triangles, are from the 2010 version of the Harris (1996) catalogue. Confirmed and candidate dwarf galaxies, represented as filled circles colour-coded according to their heliocentric distance, are taken from the 2021 version of the catalog by McConnachie (2012). The same colour-coding is applied to the filled square marking the position of Sextans II, that is also highlighted with a red contour. The stellar systems with MV and rh values that are most similar to Sextans II have been labelled. The diagonal lines are loci of constant surface brightness, from at μ = 24.0 mag arcsec−2 (top) down to μ = 34.0 mag arcsec−2 (bottom), in steps of μ = 2.0 mag arcsec−2. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.