Chemical abundance ratio of an UDG

Extreme chemical abundance ratio suggesting an exotic origin for an ultra-diffuse galaxy

Ignacio Martín-Navarro, Aaron J. Romanowsky, Jean P. Brodie Anna Ferré-Mateu, Adebusola Alabi, Duncan A. Forbes, Margarita Sharina Alexa Villaume, Viraj Pandya, and David Martinez-Delgado
University of California Observatories, 1156 High Street, Santa Cruz, CA 95064, USA
Max-Planck Institut für Astronomie, Konigstuhl 17, D-69117 Heidelberg, Germany
Department of Physics & Astronomy, San José State University, One Washington Square, San Jose, CA 95192, USA
Centre for Astrophysics & Supercomputing, Swinburne University of Technology, Hawthorn VIC 3122, Australia
Institut de Ciencies del Cosmos (ICCUB),Universitat de Barcelona (IEEC-UB), E02028 Barcelona, Spain
Special Astrophysical Observatory, Russian Academy of Sciences, N. Arkhyz, KChR 369167, Russia
Astronomisches Rechen-Institut, Zentrum für Astronomie, Universität Heidelberg, Mönchhofstr. 12–14, 69120
Accepted XXX. Received YYY; in original form ZZZ

Ultra diffuse galaxies are a population of extended galaxies but with relatively low luminosities. The origin of these objects remains unclear, largely due to the observational challenges of the low surface brightness Universe. We present here a detailed stellar population analysis of a relatively-isolated UDG, DGSAT I, based on spectroscopic data from the Keck Cosmic Web Imager integral field unit. The star formation history of DGSAT I seems to be extended, with a mean luminosity-weighted age of 3 Gyr, in agreement with previous photometric studies. However, we find a very high [Mg/Fe] abundance ratio, which is extreme even in the context of the highly alpha-enhanced massive ellipticals and ultra-faint dwarfs. The [Mg/Fe]-enhancement of DGSAT I appear to be 10 times higher than the most magnesium-enhanced stellar systems discovered to date, and suggests that the chemical enrichment of this object was dominated by core-collapse supernovae. Intriguingly, this breaks the canonical relation between [Mg/Fe] and star formation time-scale. With a measured velocity dispersion of km s, DGSAT I also shows a high mass-to-light ratio, which indicates that it is highly dark matter-dominated. The metal-poor conditions of DGSAT I may have enhanced the formation of massive stars, while at the same time, additional mechanisms are needed to prevent iron-rich yields from being recycled into stars. These results suggest that some ultra-diffuse galaxies could have experienced chemical enrichment episodes similar to the first building blocks of galaxies.

galaxies: formation – galaxies: evolution – galaxies: abundances – galaxies: stellar content
pubyear: 2018pagerange: Extreme chemical abundance ratio suggesting an exotic origin for an ultra-diffuse galaxyA.3

1 Introduction

The formation history of a galaxy leaves a permanent footprint on its chemical composition. In particular, the mass ratios of alpha-elements such as magnesium to iron ([Mg/Fe]) trace the time-scales for star formation, as these elements are produced by stars with different life-times (e.g. Worthey et al., 1992; Thomas et al., 1999; Thomas et al., 2005; Pipino & Matteucci, 2004). Both giant elliptical galaxies (e.g. de La Rosa et al., 2011; Conroy et al., 2014; La Barbera et al., 2014; McDermid et al., 2015; Martín-Navarro et al., 2018) and ultra-faint dwarfs (e.g. Norris et al., 2010; Vargas et al., 2013; Frebel et al., 2014; Frebel et al., 2016; Ishigaki et al., 2014; Roederer et al., 2016) are enriched in [Mg/Fe] by factors of up to 4 compared to the Milky Way.

Recently, a population of galaxies with the sizes of giants but the luminosities of dwarfs has received considerable attention (e.g. Merritt et al., 2014; van Dokkum et al., 2015; Koda et al., 2015; Mihos et al., 2015; Yagi et al., 2016; Mancera Piña et al., 2018; Zaritsky et al., 2018). The nature and origins of these ultra-diffuse galaxies (UDGs) are under debate (Amorisco & Loeb, 2016; Di Cintio et al., 2017; Chan et al., 2018; Rong et al., 2017; Alabi et al., 2018; Bennet et al., 2018; Carleton et al., 2018; Ogiya, 2018), and it has been suggested that they may be analogs of the primordial stellar systems that built up the halos of more luminous galaxies (Peng & Lim, 2016), with high levels of alpha elements.

There are indications that some UDGs are unusually metal-poor but detailed spectroscopic studies are still scarce due to the low-surface brightness of these objects (Kadowaki et al., 2017; Gu et al., 2018; Ruiz-Lara et al., 2018; Ferré-Mateu et al., 2018; Emsellem et al., 2018; Fensch et al., 2018), which makes them observationally challenging even for 10-meter class telescopes. Based on photometric measurements, Pandya et al. (2018) measured the stellar population properties of two UDGs. While one of these UDGs, VCC 1287, turned out to be old and metal-poor as expected for its stellar mass, the other UDG in their sample, DGSAT I, showed bluer optical colors (see Figure 2 in Pandya et al., 2018), reflecting the presence of young and/or metal-poor stellar populations. However, the lack of spectroscopic data prevented a more detailed stellar population analysis.

This relatively blue object DGSAT I was first identified by Martínez-Delgado et al. (2016), and it is particularly interesting because it apparently resides in a low density environment, at least compared to the densities where UDGs are typically found (e.g. van der Burg et al., 2016; Román & Trujillo, 2017). DGSAT I is a UDG likely associated with the outskirts of the Pisces-Perseus supercluster Martínez-Delgado et al. (2016), with a central surface brightness of mag arcsec and a total -band luminosity of . It has an optical color of 1.0, and size of kpc in the -band at a distance of 78 Mpc (Martínez-Delgado et al., 2016; Pandya et al., 2018). The estimated stellar mass of DGSAT I is 4.8 M, and it shows little sign of cold gas or current star formation (Papastergis et al., 2017).

The low-density environment where DGSAT I resides offers an opportunity to study the secular evolution of this type of object, and therefore, to understand different formation paths for UDGs without confounding environmental effects. In this paper, we make use of integral field spectroscopic data from the Keck Cosmic Web Imager (KCWI) to study in detail the stellar population properties of DGSAT I, as a representative example of a relatively isolated UDG. The layout is as follows: in Section 2 we present the data. In Sections 3 and 4 we describe the kinematic and stellar populations analysis, respectively. The main results of this paper are discussed in Section 5, and finally in Section 6 we summarize our main conclusions.

2 Data

Spectroscopic data were obtained through the Keck Cosmic Web Imager (Morrissey et al., 2018) integral field unit at Keck II Telescope, on 2017 Sept. 18,19 (PI Duncan Forbes) and Oct. 18,19 (PI Jean Brodie). We used the BM grating with the medium slicer, which gave us a spectral resolution of over a field-of-view of 16.5" 20.4" and a plate scale of 0.3 arcsec/spaxel. We covered a spectral range from to in the rest frame. The total usable integration time was 7.25 hours (on source). We show in Figure 1 the -band image of DGSAT I (Martínez-Delgado et al., 2016), with the KCWI field-of-view superimposed.

Figure 1: Archival Subaru/Suprime-Cam -band image of DGSAT I with the KCWI field-of-view marked in white. The galaxy photometric centre is marked with a cross. For reference, the yellow horizontal line has a length of 2 kpc.

Data were reduced using the kderp KCWI pipeline111 It performs a standard data reduction process that covers a basic CCD reduction (bias and over-scan subtraction, gain-correction, trimming and cosmic ray removal), flat-fielding, geometric solution, profile correction and atmospheric correction. Additionally, it performs a (relative) flux-calibration using standard stars.

Background objects were masked before collapsing the KCWI field-of-view into a single spectrum: in particular, the two brightest stars visible in Figure 1, along with a more distant galaxy only detected in emission (the extended object close to the center). In order to ensure a clean extraction of the DGSAT I spectrum, the masking process was done after collapsing all the KCWI exposures into a single image. In this way we maximized the signal-to-noise of the detection image, making sure we detected all possible objects contributing to the final spectrum. For safety, all masks extend well beyond the sources, where no more contaminated light is detected. This process was done manually in order to prevent any contamination in the spectrum of DGSAT I. Note also that the effective field of view of KCWI is slightly narrower than the nominal 16"x20" and thus, the brightest stars potentially affecting our results are barely included. No signs of additional objects were detected in the collapsed 2D cube. Due to the extremely low surface brightness of DGSAT I, we collapsed all the spaxels within the KCWI effective field of view (an aperture equivalent to  kpc assuming a 0.38 kpc/arcsec scale) into a single spectrum, reaching a signal-to-noise ratio of 30 Å.

3 Kinematics

Given the low surface brightness of DGSAT I, we carried out a more careful sky subtraction than that provided by the KCWI pipeline. We made use of pPXF (Cappellari & Emsellem, 2004; Cappellari, 2017) to fit simultaneously both sky emission and galaxy stellar absorption. This technique, particularly well suited for sky-dominated data (e.g. Weijmans et al., 2009), allowed us to minimize the residuals in our final spectra. Since a proper sky subtraction is necessary for a reliable analysis, in Appendix A.2 we tested the robustness of our results against changes in the sky subtraction process, finding no significant differences.

We fed pPXF with the PEGASE-HR stellar population synthesis models (Le Borgne et al., 2004; Prugniel & Soubiran, 2001), which have a nominal resolution of R=10000, in order to measure the potentially low velocity dispersion of DGSAT I. We found a systemic velocity km sand velocity dispersion km s. The latter implies a mass of M within a deprojected half-light radius of 5.8 kpc (Wolf et al., 2010), and an -band mass-to-light ratio of . DGSAT I is thus highly dark matter dominated in its central regions. The dynamical mass for DGSAT Iwas calculated at the half-light radius: a standard approach that is insensitive to orbital anisotropy. The specific formula used is (Wolf et al., 2010), where is in km s, is circularized and in kpc, and is in .

In order to test the robustness of this result, we measured and in three radial apertures (see Appendix A.2). We consistently recover a velocity dispersion value of 55 km s, with a standard deviation of 7 km s among these three radial apertures, in agreement with the integrated measurement quoted in the main text. Since DGSAT I was not centered in the field-of-view of KCWI (see Figure 1), the radial measurements of are sensitive to the rotation curve of this galaxy. We find slightly higher mean velocities in the outskirts ( km s) than in the center ( km s). However, the large 1- uncertainty, typically around 8 km s, prevents us from any confident claim of rotation in DGSAT I. The signal-to-noise of the radial apertures was lower () than that of the integrated spectrum. Although it is enough to confidently measure both and , it makes the radial analysis even more difficult. Its near-circular shape (Martínez-Delgado et al., 2016) further complicates the identification of a clear rotation axis.

4 Stellar population analysis

4.1 Star formation history

We derived the star formation history (SFH) of DGSAT I by fitting our spectrum with STECKMAP (Ocvirk et al., 2006a, b) plus, as for the kinematics, the PEGASE-HR set of templates. In short, STECKMAP finds the best-fitting linear combination of single stellar populations models that reproduce a given spectrum. Recovering the SFH from an observed spectrum is an ill-defined problem, and therefore regularization is needed. The level of regularization imposed by STECKMAP is tuned via the dimensionless and parameters, which effectively behave as Gaussian priors on the SFH and on the age–metallicity relation, respectively(Ocvirk et al., 2006b). It is important to note here that age information is not only encoded in Balmer lines, but is also distributed along the whole optical wavelength range (Ocvirk et al., 2006a). We assumed , imposed by the signal-to-noise of the data (see e.g. Sánchez-Blázquez et al., 2011) However, there is not a unique way of varying the amount of regularization and thus, the recovered SFH of DGSAT I may slightly depend on the adopted values for and . Appendix A.1 demonstrates how our conclusions do not depend on the assumed and values.

The best-fitting model and residuals are shown in Figure 2, where the most prominent absorption features are clearly noticeable. DGSAT I shows an extended SFH, with a mass-weighted age of 8.1 0.4 Gyr, and luminosity-weighted age of 2.8 0.5 Gyr. It took Gyr for DGSAT I to form 50% of its stellar mass. This result is in agreement with the photometric study of Pandya et al. (2018), where they estimated a luminosity-weighted age of 3 Gyr.

Figure 2: The KCWI (rest-frame) spectrum of DGSAT I is shown in black, along with the best-fitting model derived with STECKMAP (red), and the residuals (blue). The spectrum was normalized so its median value was equal to one. The most prominent spectral features are clearly visible (vertical grey bands).

In addition, we also measured the SFH of DGSAT I using a regularized pPXF solution (Cappellari, 2017). This minimization scheme is different than the STECKMAP implementation, and it allows for multiple metallicities per age bin. Thus, a metallicity spread at a given age should be captured by pPXF. With this additional method, we recover a luminosity-weighted age of 3.5 Gyr, and a mass-weighted value of 7.9 Gyr, in excellent agreement with the STECKMAP analysis. Similar results were found using STARLIGHT(Cid Fernandes et al., 2005). In order to quantify how much an extended SFH improves the fit over a single stellar population model, we compare the residuals to those expected from a Gaussian distribution with mean (i.e. age = 3.5 Gyr) and variance . A Kolmogorov–Smirnov test indicates that an extended star formation is rejected at a level of , while the single-age population is strongly disfavored at a level. This test, along with the consistency among different algorithms, shows that DGSAT I experienced an extended SFH. The cumulative SFH of DGSAT I as measured from its optical spectrum is shown in Fig. 3.

Figure 3: Cumulative SFH of DGSAT I. Horizontal dashed lines indicate when this object had already formed 95%, 75%, and 50% of its stellar mass, respectively. The width of the line corresponds to the estimated 1 uncertainty.

Qualitatively, the young luminosity-weighted ages measured in DGSAT I using both STECKMAP, STARLIGHT, and pPXF suggest a long-lasting star formation activity in this object. Since even low-mass halos in the local Universe seem to host a fraction of very old stars (Weisz et al., 2011), it is expected that the star formation in DGSAT I roughly lasted from the formation of these first stars in the early Universe until, at least, 3 Gyr ago.

4.2 Chemical composition

DGSAT I’s optical spectrum (Figure 4) reveals the peculiar chemical composition of this object. When an -enhanced ([Mg/Fe] = +0.4), metal-poor ([M/H] = ) stellar population model (red line in Fig. 2) is subtracted from the observed spectrum, the residuals associated with iron-sensitive and magnesium-sensitive features behave in the opposite manner. Magnesium lines are under-predicted by the model (negative residuals), while iron lines are over-predicted (positive residuals). For comparison, Figure 4 also shows the residuals of metal-poor, -enhanced globular clusters (Usher et al., 2017). The observed spectrum of DGSAT I and that of the globular clusters were divided by a constant factor so their median is equal to one. Note that best-fitting models account for possible systematics in the (relative) flux calibration of the data (see Appendix A.3).

Figure 4: Empirical assessment of the abundance pattern of DGSAT I. In black we show the residuals between the observed spectrum of DGSAT I and a metal-poor ([M/H] = ), -enhanced ([Mg/Fe] = ) stellar population model at rest-frame. Vertical colored regions indicate the index band-pass definition. While magnesium lines (in blue) are under-predicted by the model, iron lines (in red) behave in the opposite way, showing that the stellar populations of DGSAT I are more -enhanced than the model itself. For comparison, we present the same ratio but for metal-poor globular clusters in the Milky Way (magenta) (Usher et al., 2017). In this case, the observed spectrum is well represented by the stellar population model, without any additional Mg-enhancement. The typical 1- uncertainty for DGSAT I spectrally binned data is shown as an error bar. Data and models were convolved with a 2.5 Å-width boxcar kernel.

In order to quantify the chemical properties of DGSAT I we followed the standard line-strength analysis. We partially broke the age–metallicity degeneracy by fixing the age to those (luminosity- and mass-weighted) obtained with STECKMAP. The best-fitting metallicity was then derived by comparing the equivalent width of metallicity-sensitive lines in our spectrum to the PEGASE-HR model predictions (at fixed age), as show in Fig. 5. The typical error in the line-strength measurements was Å, which leads to the uncertainties shown in Fig. 5.

Figure 5: Line-strength analysis. The best-fitting metallicities based on magnesium lines are shown with blue symbols, while metallicities measured from iron lines are shown in red. Open and filled symbols result from assuming the mass-weighted (8.1 Gyr) and the luminosity-weighted (2.8 Gyr) ages measured with STECKMAP, respectively. For clarity, ages are vertically shifted according to the 1- uncertainties (vertical error bars). Independently of the assumed age, magnesium-sensitive lines indicate a much higher metallicity than that derived from iron-sensitive features, clearly revealing the abnormal [Mg/Fe] enhancement of DGSAT I. Note that the stellar population model at these metallicities is already -enhanced.

We found a metallicity value of [M/H] = dex assuming a mass-weighted age of 8.1 Gyr, as expected for a low- galaxy. However, metallicity measurements based on magnesium spectral features ([M/H]) are higher than the metallicities estimated from iron lines ([M/H]), as shown in Figure 5. Note that PEGASE-HR stellar population models are based on solar neighborhood stars, which are magnesium-rich in the low-metallicity regime. Hence, given the low-metallicity of DGSAT I, iron- and magnesium-based metallicities were in practice measured with respect to stellar population models with [Mg/Fe] . This implies that the stellar populations of DGSAT I have an [Mg/Fe] enhancement larger than the metal-poor stars in the solar neighborhood, i.e., larger than +0.4 dex. In the Appendix (Fig. 11) we test the robustness of the measured [Mg/Fe] against flux calibration issues.

Both the empirical analysis of the DGSAT I spectrum and the difference between magnesium- and iron-based metallicities highlight the abnormal chemical enrichment of this UDG. A zero-order estimate of its [Mg/Fe] ratio can be done using the proxy (Vazdekis et al., 2015), and leads to a value of [Mg/Fe] dex. This estimate should be interpreted as a lower limit, since it does not take into account the decreasing sensitivity to [Mg/Fe] with decreasing metallicity (Vazdekis et al., 2015). To fully quantify the [Mg/Fe] ratio of DGSAT I, we made use of the new set of MILES stellar population models (Vazdekis et al., 2010; Vazdekis et al., 2015). This set of models consistently treat the effect of the [/Fe] ratio across a wide range of metallicities, from [M/H] = to [M/H] = +0.26. The MILES resolution is, however, lower than our data. We therefore extrapolated the line-strength MILES predictions to match the resolution of the KCWI spectrograph. This correction, typically times smaller than the error in the equivalent widths, is negligible compared to the intrinsic uncertainty of the data.

The [Mg/Fe] measurement was done by fitting the average of magnesium and iron lines equivalent widths, as shown in Fig. 6, assuming a luminosity weighted age of 3.5 Gyr for the MILES -enhanced models. We fixed the slope of the stellar initial mass function to that of the Milky Way (Kroupa, 2002). The measured value for DGSAT I corresponds to an iron metallicity of [Fe/H] , and an extreme magnesium enhancement of [Mg/Fe] . These values should however be interpreted with caution given the extreme extrapolations of the models predictions. The relative comparison shown in Fig. 5 is a more robust evidence of the peculiar abundance pattern of DGSAT I. The agreement between the results obtained using PEGASE-HR and MILES models indicates that the extreme [Mg/Fe] ratio of DGSAT I is not related to systematics in the stellar population modeling.

Figure 6: The index–index grid based on the MILES models shows that the stellar populations of DGSAT I are metal-poor ([M/H]) and extremely [Mg/Fe] enhanced ([Mg/Fe] ). The lower sensitivity of the [Mg/Fe] at low metallicities is also clear, as the same difference in [Mg/Fe] implies a weaker variation in the line-strength indices than in the metal-rich regime.

5 Discussion

The stellar population properties of DGSAT I appear to contradict one another. Under a canonical interpretation, a high [Mg/Fe] ratio is produced in short-lasting star formation events (Thomas et al., 2005), which cease before Type Ia supernovae (SN Ia) explode and pollute the inter-stellar medium (ISM) with iron-rich ejecta. However, the SFH of DGSAT I is very extended, lasting for more than 10 Gyr, which leaves plenty of time for SN Ia to enrich the ISM with iron.

Reconciling a long-lasting SFH with such a high [Mg/Fe] ratio suggests exotic processes regulating the chemical evolution of DGSAT I. Since magnesium is mainly produced in massive stars, it is natural to relate [Mg/Fe] enhancement to a stellar initial mass function (IMF) dominated by massive stars, i.e. a top-heavy IMF. Observationally, there is evidence suggesting that metal-poor galactic systems, as measured in this UDG, tend to favor the formation of massive stars (Worthey et al., 1992; Geha et al., 2013; Martín-Navarro et al., 2015). Hence, a relatively flat IMF would naturally alleviate the tension between the formation time-scale and the chemical composition of DGSAT I, although it would also speed up its chemical enrichment. The metal-poor environment within DGSAT I may have also favored the formation of peculiar SNe whose ejecta are overly magnesium-enhanced (Aoki et al., 2007).

The lack of iron, however, is harder to reconcile with an extended SFH. It appears that DGSAT I efficiently retained the magnesium-rich yields of SN II, but it was unable to recycle iron-rich SN Ia stellar ejecta. Such a selective mass loss could be consistent with a scenario where the large sizes of UDGs are due to supernovae-driven gas outflows, expanding both stellar and dark matter components (Di Cintio et al., 2017; Chan et al., 2018). Under this scenario, it would be possible that magnesium-rich yields were produced in a deeper potential well, before the galaxy expansion. In contrast, SN Ia would have exploded in a much shallower halo, facilitating the escape of their iron-rich ejecta. Moreover, the cooling time of stellar ejecta increases with decreasing density and metallicity (Martizzi et al., 2015), further delaying the recycling of iron-rich ejecta after the initial galaxy expansion.

Alternatively, it could be argued that the low stellar density of DGSAT I may suppress the formation of binary stars and therefore, of SNe Ia. However, at fixed metallicity, the [Mg/Fe] ratio does not seem to increase with decreasing stellar density (Shetrone et al., 2001; Pritzl et al., 2005), suggesting that iron-rich stellar ejecta are recycled even in very low-mass galaxies (Vargas et al., 2013). More extreme processes to explain the lack of iron in DGSAT I such as filamentary (Mg-rich) gas accretion from the cosmic web could also be possible (Sánchez Almeida, 2017), although it would require a constant cold accretion over a period of Gyr. Note also that this process is not unique to UDGs and would also apply to other (dwarf) galaxies where the chemical composition is known to be less extreme.

Fig. 7 shows how DGSAT I compares to the general population of galaxies. Similarly to some other UDGs in denser environments (Beasley et al., 2016; van Dokkum et al., 2016; Toloba et al., 2018), DGSAT I is heavily dark-matter dominated, deviating strongly from the scaling relation of normal galaxies. Its abnormal chemical composition is also clear from Fig. 7, where DGSAT I exhibits a [Mg/Fe] 1.0 dex higher than any other galaxy known so far. Interestingly, there is tentative evidence of a trend within the UDG population, with more metal-poor objects being more [Mg/Fe] enhanced, similarly to other low-mass galaxies (Conselice et al., 2003; Penny & Conselice, 2008; Ferré-Mateu et al., 2018) and star in the Milky Way. Only a handful of resolved Milky Way stars (Aoki et al., 2007) present [Mg/Fe] values similar to that measured in DGSAT I. The similarity with these stars, likely descendants of the first core-collapse SNe (Audouze & Silk, 1995), further supports the pristine character of DGSAT I. However, it is not yet understood how its stellar populations have evolved unaltered after 10 Gyr of star formation. Neither the metallicity of DGSAT I, nor its mass, nor its size alone would be enough to explain how this galaxy had an extended SFH while preserving a high [Mg/Fe] ratio. Therefore, the evolutionary history of DGSAT I is likely to be the differentiating factor. If so, it would put strong constraints on the formation of low-density cores within low-mass galaxies via supernovae feedback (e.g. Governato et al., 2010; Pontzen & Governato, 2012), as this process should be then reflected in their chemical abundance patterns.

Figure 7: Scaling relations. In the upper panel we show the dynamical mass-to-light ratio vs. the dynamical mass, both within a half-light radius. DGSAT I and other UDGs (red points) are outliers relative to other objects (grey), with an excess of dark matter. In the bottom panel we show the [Mg/Fe] ratio vs. iron metallicity, for field stars and globular clusters (Pritzl et al., 2005) in the Milky Way and in satellite galaxies (Shetrone et al., 2001; Vargas et al., 2013; Aoki et al., 2007). A massive early-type galaxy at high redshift is also shown (Kriek et al., 2016), along with Coma UDGs (Ferré-Mateu et al., 2018). They follow the Milky Way relation of increasing [Mg/Fe] at lower metallicities, with two UDGs resembling average halo stars. DGSAT I is more extreme, with lower metallicity and higher [Mg/Fe].

6 Summary and conclusions

DGSAT I is a relatively isolated UDG with a measured velocity dispersion of km s, indicating a high dark matter-to-stellar mass ratio. Moreover, DGSAT I exhibits an extended SFH, in agreement with its blue colors. However, this long-lasting star formation time-scale is in apparent contradiction with its chemical properties, as DGSAT I hosts extremely magnesium-enhanced stellar populations. Reconciling such a high [Mg/Fe] abundance ratio, which seems to be a common feature of at least some UDGs in denser environments (Ruiz-Lara et al., 2018; Ferré-Mateu et al., 2018), with an extended SFH requires an inefficient recycling process of Fe-rich SN Ia stellar ejecta, either due to the low density of DGSAT I, or as a consequence of the particular formation path of this object.

Regardless of the exact mechanism(s) regulating the chemical evolution of DGSAT I, its stellar populations may provide a window into the early Universe. The lack of iron is an indication of undeveloped chemical evolution, suggesting that the star formation conditions were similar to those of primeval galaxies. Spectroscopic observations of larger samples of UDGs (Ruiz-Lara et al., 2018; Ferré-Mateu et al., 2018) are needed in order to further investigate the connections between star formation, abundance patterns and formation pathways for this class of exotic objects (Amorisco & Loeb, 2016; Di Cintio et al., 2017; Chan et al., 2018), in particular for isolated UDGs.


We acknowledge support from the National Science Foundation grants AST-1616598, AST-1616710, Marie Curie Global Fellowship, and from grants AYA2016-77237-C3-1-P and AYA2014-56795-P from the Spanish Ministry of Economy and Competitiveness (MINECO). AJR was supported as a Research Corporation for Science Advancement Cottrell Scholar. DAF and AFM thank the ARC for financial support via DP160101608. also acknowledges financial support through the Postdoctoral Junior Leader Fellowship Programme from ”la Caixa” Banking Foundation (LCF/BQ/LI18/11630007). MES is supported by the RFBR grant no. 18-02-00167 A. 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. DMD acknowledges support by Sonderforschungsbereich (SFB) 881 “The Milky Way System” of the German Research Foundation (DFG) (sub-project A2). We would like to thank P. van Dokkum, R. Leaman, G. van de Ven, A. L. Frebel and S. Penny for their useful comments. 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.


  • Alabi et al. (2018) Alabi A., et al., 2018, MNRAS, 479, 3308
  • Amorisco & Loeb (2016) Amorisco N. C., Loeb A., 2016, MNRAS, 459, L51
  • Aoki et al. (2007) Aoki W., Beers T. C., Christlieb N., Norris J. E., Ryan S. G., Tsangarides S., 2007, ApJ, 655, 492
  • Audouze & Silk (1995) Audouze J., Silk J., 1995, ApJ, 451, L49
  • Beasley et al. (2016) Beasley M. A., Romanowsky A. J., Pota V., Navarro I. M., Martinez Delgado D., Neyer F., Deich A. L., 2016, ApJ, 819, L20
  • Bennet et al. (2018) Bennet P., Sand D. J., Zaritsky D., Crnojević D., Spekkens K., Karunakaran A., 2018, ApJ, 866, L11
  • Cappellari (2017) Cappellari M., 2017, MNRAS, 466, 798
  • Cappellari & Emsellem (2004) Cappellari M., Emsellem E., 2004, PASP, 116, 138
  • Carleton et al. (2018) Carleton T., Errani R., Cooper M., Kaplinghat M., Peñarrubia J., 2018, preprint, (arXiv:1805.06896)
  • Cenarro et al. (2008) Cenarro A. J., Cervantes J. L., Beasley M. A., Marín-Franch A., Vazdekis A., 2008, ApJ, 689, L29
  • Chan et al. (2018) Chan T. K., Kereš D., Wetzel A., Hopkins P. F., Faucher-Giguère C.-A., El-Badry K., Garrison-Kimmel S., Boylan-Kolchin M., 2018, MNRAS, 478, 906
  • Cid Fernandes et al. (2005) Cid Fernandes R., Mateus A., Sodré L., Stasińska G., Gomes J. M., 2005, MNRAS, 358, 363
  • Coelho et al. (2009) Coelho P., Mendes de Oliveira C., Cid Fernandes R., 2009, MNRAS, 396, 624
  • Conroy et al. (2014) Conroy C., Graves G. J., van Dokkum P. G., 2014, ApJ, 780, 33
  • Conroy et al. (2018) Conroy C., Villaume A., van Dokkum P. G., Lind K., 2018, ApJ, 854, 139
  • Conselice et al. (2003) Conselice C. J., Gallagher III J. S., Wyse R. F. G., 2003, AJ, 125, 66
  • Di Cintio et al. (2017) Di Cintio A., Brook C. B., Dutton A. A., Macciò A. V., Obreja A., Dekel A., 2017, MNRAS, 466, L1
  • Emsellem et al. (2018) Emsellem E., et al., 2018, arXiv e-prints, p. arXiv:1812.07345
  • Faber et al. (1985) Faber S. M., Friel E. D., Burstein D., Gaskell C. M., 1985, ApJS, 57, 711
  • Fensch et al. (2018) Fensch J., et al., 2018, arXiv e-prints,
  • Ferré-Mateu et al. (2018) Ferré-Mateu A., et al., 2018, MNRAS, 479, 4891
  • Frebel et al. (2014) Frebel A., Simon J. D., Kirby E. N., 2014, ApJ, 786, 74
  • Frebel et al. (2016) Frebel A., Norris J. E., Gilmore G., Wyse R. F. G., 2016, ApJ, 826, 110
  • Geha et al. (2013) Geha M., et al., 2013, ApJ, 771, 29
  • Governato et al. (2010) Governato F., et al., 2010, Nature, 463, 203
  • Gu et al. (2018) Gu M., et al., 2018, ApJ, 859, 37
  • Ishigaki et al. (2014) Ishigaki M. N., Aoki W., Arimoto N., Okamoto S., 2014, A&A, 562, A146
  • Kadowaki et al. (2017) Kadowaki J., Zaritsky D., Donnerstein R. L., 2017, ApJ, 838, L21
  • Koda et al. (2015) Koda J., Yagi M., Yamanoi H., Komiyama Y., 2015, ApJ, 807, L2
  • Koleva et al. (2008) Koleva M., Prugniel P., Ocvirk P., Le Borgne D., Soubiran C., 2008, MNRAS, 385, 1998
  • Kriek et al. (2016) Kriek M., et al., 2016, Nature, 540, 248
  • Kroupa (2002) Kroupa P., 2002, Science, 295, 82
  • La Barbera et al. (2014) La Barbera F., Pasquali A., Ferreras I., Gallazzi A., de Carvalho R. R., de la Rosa I. G., 2014, MNRAS, 445, 1977
  • Le Borgne et al. (2004) Le Borgne D., Rocca-Volmerange B., Prugniel P., Lançon A., Fioc M., Soubiran C., 2004, A&A, 425, 881
  • Mancera Piña et al. (2018) Mancera Piña P. E., Peletier R. F., Aguerri J. A. L., Venhola A., Trager S., Choque Challapa N., 2018, MNRAS, 481, 4381
  • Martín-Navarro et al. (2015) Martín-Navarro I., et al., 2015, ApJ, 806, L31
  • Martín-Navarro et al. (2018) Martín-Navarro I., Vazdekis A., Falcón-Barroso J., La Barbera F., Yıldırım A., van de Ven G., 2018, MNRAS, 475, 3700
  • Martínez-Delgado et al. (2016) Martínez-Delgado D., et al., 2016, AJ, 151, 96
  • Martizzi et al. (2015) Martizzi D., Faucher-Giguère C.-A., Quataert E., 2015, MNRAS, 450, 504
  • McDermid et al. (2015) McDermid R. M., et al., 2015, MNRAS, 448, 3484
  • Merritt et al. (2014) Merritt A., van Dokkum P., Abraham R., 2014, ApJ, 787, L37
  • Mihos et al. (2015) Mihos J. C., et al., 2015, ApJ, 809, L21
  • Morrissey et al. (2018) Morrissey P., et al., 2018, ApJ, 864, 93
  • Norris et al. (2010) Norris J. E., Yong D., Gilmore G., Wyse R. F. G., 2010, ApJ, 711, 350
  • Ocvirk (2010) Ocvirk P., 2010, ApJ, 709, 88
  • Ocvirk et al. (2006a) Ocvirk P., Pichon C., Lançon A., Thiébaut E., 2006a, MNRAS, 365, 46
  • Ocvirk et al. (2006b) Ocvirk P., Pichon C., Lançon A., Thiébaut E., 2006b, MNRAS, 365, 74
  • Ogiya (2018) Ogiya G., 2018, MNRAS, 480, L106
  • Pandya et al. (2018) Pandya V., et al., 2018, ApJ, 858, 29
  • Papastergis et al. (2017) Papastergis E., Adams E. A. K., Romanowsky A. J., 2017, A&A, 601, L10
  • Peng & Lim (2016) Peng E. W., Lim S., 2016, ApJ, 822, L31
  • Penny & Conselice (2008) Penny S. J., Conselice C. J., 2008, MNRAS, 383, 247
  • Pipino & Matteucci (2004) Pipino A., Matteucci F., 2004, MNRAS, 347, 968
  • Pontzen & Governato (2012) Pontzen A., Governato F., 2012, MNRAS, 421, 3464
  • Pritzl et al. (2005) Pritzl B. J., Venn K. A., Irwin M., 2005, AJ, 130, 2140
  • Prugniel & Soubiran (2001) Prugniel P., Soubiran C., 2001, A&A, 369, 1048
  • Roederer et al. (2016) Roederer I. U., et al., 2016, AJ, 151, 82
  • Román & Trujillo (2017) Román J., Trujillo I., 2017, MNRAS, 468, 703
  • Rong et al. (2017) Rong Y., Guo Q., Gao L., Liao S., Xie L., Puzia T. H., Sun S., Pan J., 2017, MNRAS, 470, 4231
  • Ruiz-Lara et al. (2018) Ruiz-Lara T., et al., 2018, MNRAS, 478, 2034
  • Sánchez Almeida (2017) Sánchez Almeida J., 2017, in Fox A., Davé R., eds, Astrophysics and Space Science Library Vol. 430, Gas Accretion onto Galaxies. p. 67
  • Sánchez-Blázquez et al. (2011) Sánchez-Blázquez P., Ocvirk P., Gibson B. K., Pérez I., Peletier R. F., 2011, MNRAS, 415, 709
  • Shetrone et al. (2001) Shetrone M. D., Côté P., Sargent W. L. W., 2001, ApJ, 548, 592
  • Thomas et al. (1999) Thomas D., Greggio L., Bender R., 1999, MNRAS, 302, 537
  • Thomas et al. (2005) Thomas D., Maraston C., Bender R., Mendes de Oliveira C., 2005, ApJ, 621, 673
  • Toloba et al. (2018) Toloba E., et al., 2018, ApJ, 856, L31
  • Usher et al. (2017) Usher C., et al., 2017, MNRAS, 468, 3828
  • Vargas et al. (2013) Vargas L. C., Geha M., Kirby E. N., Simon J. D., 2013, ApJ, 767, 134
  • Vazdekis et al. (2010) Vazdekis A., Sánchez-Blázquez P., Falcón-Barroso J., Cenarro A. J., Beasley M. A., Cardiel N., Gorgas J., Peletier R. F., 2010, MNRAS, 404, 1639
  • Vazdekis et al. (2015) Vazdekis A., et al., 2015, MNRAS, 449, 1177
  • Weijmans et al. (2009) Weijmans A.-M., et al., 2009, MNRAS, 398, 561
  • Weisz et al. (2011) Weisz D. R., et al., 2011, ApJ, 739, 5
  • Wolf et al. (2010) Wolf J., Martinez G. D., Bullock J. S., Kaplinghat M., Geha M., Muñoz R. R., Simon J. D., Avedo F. F., 2010, MNRAS, 406, 1220
  • Worthey et al. (1992) Worthey G., Faber S. M., Gonzalez J. J., 1992, ApJ, 398, 69
  • Yagi et al. (2016) Yagi M., Koda J., Komiyama Y., Yamanoi H., 2016, ApJS, 225, 11
  • Zaritsky et al. (2018) Zaritsky D., et al., 2018, preprint, (arXiv:1811.04098)
  • de La Rosa et al. (2011) de La Rosa I. G., La Barbera F., Ferreras I., de Carvalho R. R., 2011, MNRAS, 418, L74
  • van Dokkum et al. (2015) van Dokkum P. G., Abraham R., Merritt A., Zhang J., Geha M., Conroy C., 2015, ApJ, 798, L45
  • van Dokkum et al. (2016) van Dokkum P., et al., 2016, ApJ, 828, L6
  • van der Burg et al. (2016) van der Burg R. F. J., Muzzin A., Hoekstra H., 2016, A&A, 590, A20

Appendix A Robustness of the stellar population analysis

a.1 Star formation histories

A key point in our analysis is the fact that DGSAT I exhibits an extended SFH. To investigate to what extent our choice of the regularization parameters affects the measured SFH, we repeated the STECKMAP analysis but assuming . In Fig. 8 we show the cumulative mass distribution of DGSAT I with () and without () regularization. In both cases, the star formation lasts for more than 10 Gyr, appearing even more extended in the case without regularization. Hence, regularization is not responsible for the long-lasting SFH of DGSAT I. As an additional test, we also measured the SFH from previously published optical spectra (Martínez-Delgado et al., 2016). Although in this case the spectral coverage was not enough to measure the abundance pattern of DGSAT I, we obtained ages and metallicities in agreement with the KCWI data.

Figure 8: The extended SFH of DGSAT I. In green we show the fiducial cumulative mass function for DGSAT I, assuming . Horizontal dashed lines mark when 50%, 75% and 95% of the stellar mass was formed. In blue we show the cumulative mass function but leaving the STECKMAP solution almost un-regularized (), which tends to favor even more extended star formation histories.Our choice of the regularization parameters is therefore not spuriously delaying the mass assembly of this object.

Finally, it is worth noting that a significant population of blue horizontal branch (BHB) stars in metal-poor objects could artificially bias measured ages towards younger values (Cenarro et al., 2008; Conroy et al., 2018), which may affect the recovered SFH of DGSAT I. This bias is introduced because BHB stars do not follow canonical isochrones, being hotter than expected for their luminosity and thus, they are not properly modelled by standard stellar population synthesis models. Effectively, this means that they can be misinterpreted as a young (hot) component. However, it has been shown that, using STECKMAP, the net effect of BHB stars is to resemble a burst of recent (age 0.1 Gyr) star formation (Ocvirk, 2010). Similar results (i.e. a burst of young stars rather than a more extended star formation) have been found for similar algorithms (Coelho et al., 2009), which, in practice, could be used to identify BHB stars in unresolved stellar populations (Koleva et al., 2008). Hence, our measurement of an extended SFH in DGSAT I is not likely due to an inaccurate modeling of the horizontal branch morphology.

a.2 Line-strength analysis

Because of the unusually high [Mg/Fe] measurement, we have thoroughly looked for systematic effects on the data. With a relatively low signal-to-noise of 30 per Å, metallicity and abundance pattern measurements based on individual spectral features are unreliable. However, as shown in Fig. 5, all magnesium and iron line residuals have opposite behaviors. This is a crucial point, since different line-strengths have different band-passes and pseudo-continuum definitions. The fact that magnesium-based metallicities are dex higher than iron-based ones, independently of the index choice, shows that our [Mg/Fe] estimation is not driven by data systematics such as flux calibration issues.

We made two extra tests to further assess the robustness of our [Mg/Fe] measurement. First, we randomly selected 20 of the 36 individual exposures of DGSAT I and then measured the stellar population properties only over this subset of the data. We repeated this process five times, making sure that our results are not driven by anomalous frames. However, we find that, irrespectively of the sub-sample of frames, magnesium-based and iron-based metallicities systematically depart from each other, as shown in Fig. 9. We also repeated this test but stacking consecutive exposures, and found the same result.

As a final test, we attempted a radial analysis of our IFU data, defining 3 circular apertures (annuli) at kpc, kpc, and kpc from the center of the galaxy. If the [Mg/Fe] over-enhancement of DGSAT I were to due to small-scale variations in the flux calibration, [Mg/Fe] measurements should abruptly change across the detector. However, we find a rather constant magnesium-enhancement over the whole KCWI field-of-view. There is only marginal evidence of a slightly higher [Mg/Fe] ratio in the outskirts of DGSAT I, which would be in agreement with a standard outside-in formation scenario (Pipino & Matteucci, 2004). Also, note that the metallicity does not significantly change across the galaxy, although the large uncertainties in the metallicity measurement prevent us from obtaining a reliable radial gradient. That indicates that our measurements are not driven by the off-center over-density in the North-West of DGSAT I (Martínez-Delgado et al., 2016).

Figure 9: Robustness of the [Mg/Fe] measurement. Differences between magnesium-based and iron-based metallicities (Z and Z, respectively) for the tests described above. In black we show our reference value. In blue we show the five bootstrapping realizations, where the measurement was done over 20 randomly selected frames out of the 36 individual exposures. In red, we show the radial profile, from the centre (R), to the outskirts (R kpc). Dark green symbol shows the measurement after removing the index continua (see details below), while the two light green points are the result after varying the second-order polynomial correction to the continuum. Low corresponds to a low-order polynomial, while high indicates a higher order correction. The agreement between all these measurements indicates that systematics in our data are not likely driving the observed abundance pattern. Error bars correspond to 1- uncertainties.

For clarity, Fig. 10 shows the three circular apertures mentioned above. Note that they are completely independent, i.e., they are not cumulative apertures. We also show the effective field of view after trimming the edges. The masked regions containing visible background/foreground objects are marked in red.

Figure 10: Radial apertures and masks. As in Fig. 2 but showing the effective field of view of KCWI (dashed line) and the masked background/foreground objects (red circles). The three radial annuli apertures described above (R, R, and R) are indicated with doted white circles.

a.3 Spectrophotometric calibration

A precise calibration of the KCWI data across all wavelengths is needed in order to robustly measure abundance pattern variations. Relative differences in the flux level could be artificially interpreted as changes in the measured equivalent widths, becoming a potential source of systematic errors in our analysis. In practice, we tried to minimize this effect by combining line-strength over a wide range of wavelengths. In Fig. 5 we show how all magnesium lines predict higher metallicities than iron lines. This would be hard to explain solely by flux calibration issues since magnesium-sensitive and iron-sensitive features are spread across the optical range. In order to mimic this effect, any flux calibration issue would have to precisely follow the location of the different lines, changing the relative flux response in such a way that all magnesium and iron lines behave in the same way. Moreover, the Mg 4780 is located 200Å away from the other lines, but it shows similar variations. Note also that while measuring the kinematics a multiplicative polynomial was applied to the data in order to match the (flux-calibrated) set of stellar population models. Before measuring any equivalent width, we corrected our data for this additional term. This second-order correction minimizes the effect of systematics introduced by, e.g., the spectral type of our standard star, as stellar population models have been accurately flux calibrated. Finally, it is worth noting that the central band-pass of the Mgb and Mg2 features is roughly the same, but the blue and red pseudo-continua, are significantly different (Faber et al., 1985). We intentionally included these two line-strengths in our analysis to account for flux calibration issues. Since both line indices, probing the same spectral feature but based on different continua show consistent metallicity measurements, our results seem to be robust again uncertainties in the spectrophotometry calibration.

The robustness of our analysis is also evident from Fig. 4. The residuals in this figure empirically show the peculiar [Mg/Fe] abundance ratio of DGSAT I after subtracting the best-fitting model from STECKMAP. This best-fitting model includes a continuum term that corrects for flux calibration issues, which in practice means that the differences shown in Fig. 4 are not likely a consequence of the data reduction/calibration.

Given the sensitivity of our analysis to potential flux calibration issues, we have done three additional tests to complement all the previous arguments. First, we have repeated the analysis shown in Fig. 5, but this time removing any information encoded in the continuum. The process was the following: for each line-strength index, we used blue and red pseudo-continua to fit a straight line. Then, we divided the spectrum around the index by this best-fitting straight line in order to remove the continuum, and we finally measured the equivalent width. While similar to the standard way of measuring line-strengths, this approach remove any short-scale continuum variation, making the analysis virtually insensitive to flux calibration issues. We did the same (independently) to the PEGASE-HR stellar population templates, and we then compared data and models. Finally, we have also increased and decreased the order of the polynomial used in the kinematics analysis (see above). By doing this we can assess the effect of this correction on the measured abundance pattern. The result of all these tests are presented in Fig 11, showing the robustness of our approach. It is worth mentioning that Z-Z in Fig 11 already corresponds to [Mg/Fe] = +0.4, i.e., that of inherited from the (metal-poor) stellar population models. Hence, the measured positive values indicate an even higher [Mg/Fe] value for DGSAT I.

Figure 11: Robustness of the [Mg/Fe] measurement: flux calibration. Same as in Fig. 5 but with equivalent widths measured after removing the continuum in data and models using the index pseudo-continua definitions. Removing any short-scale variation in the spectrum of DGSAT I does not change our main result, showing that the extreme [Mg/Fe] measurement is not related to systematic uncertainties in the spectrophotometric calibration. Error bars correspond to 1- uncertainties.

The separation between magnesium and iron lines changes but not significantly enough to alter our conclusions, further confirming that the extreme [Mg/Fe] value measured for DGSAT I is not driven by a poor flux calibration. The rightmost symbol in Fig. 9 shows the difference between magnesium-based and iron-based metallicities measured after removing the index continua, compared to all the other tests described above.

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description