arXiv:1711.11299v1 [astro-ph.SR] 30 Nov 2017 © oraogteLclGopglxe htssansa format star sustain that mos the Group (e.g. of Local one the st as among massive out poor stands resolved promis- population metal-poor, a intermediate-age very makes Its find on) to now environment from ing SagDIG 1927-17.7; UKS 594-4, cetdXX eevdYY noiia omZZZ form original in YYY; Received XXX. Accepted iimGarcia Miriam MNRAS ⋆ eima ihrrdhfs( inter-g the higher of at measurements medium from and history formation star 2014 the 2010 asv tr nteSgtaisDafIrglrGalaxy Irregular Dwarf Sagittarius the in stars Massive 1 etn ( A Sextans nttt eAtoíiad aais nteiln fL Pal La of island Mucha the los on † de Canarias, de Roque Astrofísica del de Observatorio Instituto Spanish the in installed sSClk ( iro SMC-like their that is indicate studies subsequent but abundances, o trahwt ut-betsetorpsa –0 tele- IC 8–10m the ( galaxies to at WLM directed ar and spectrographs 3109 were stars NGC analyses multi-object massive quantitative with metal-poorer First scopes. but reach regime, at low-Z now the for dard Group. Local the of galaxies irregular dwarf metal mostly decreasing of environments e in stars enthusiastic massive an resolved s stimulated Cos first has the previous of This physics Universe. at the the towards feedback step and necessary a formation and sta epochs, star massive interpret metal-poor to very key of evolution the Understanding INTRODUCTION 1 1 / etod srbooí CI-NA.Cta eTreó A a Torrejón de Crtra. (CSIC-INTA). Astrobiología de Centro 0O 10 -al [email protected] E-mail: ae nosrain aewt h rnTlsoi Canaria Telescopio Gran the with made observations on Based 07TeAuthors The 2017 ∼ h aitru wr reua aay( Irregular Dwarf Sagittarius The ; ; h ml aelncCod(M)i h urn stan- current the is (SMC) Cloud Magellanic Small The conci 2012 McConnachie 1 rme ta.2011 al. et Tramper ⊙ / orte l 2015 al. et Bouret 0Z 30 000 ( afre l 2004 al. et Kaufer , ⊙ 1 aah ta.2016 al. et Camacho – ∼ 5 vrg ealct fteUies ttepa of peak the at Universe the of average 1 (2017) / Fe 5 ⊙ 1 nta ( instead ) ,wt ausifre rmtelocation the from inferred values with ), † ihecuaigylw1 low encouragingly with ) .Floige Following ). aa ikno 2014 Dickinson & Madau .Hwvr hsvlei tl a from far still is value this However, ). rslne l 2007 al. et Bresolin acae l 2014 al. et Garcia ihbt 1 both with ) h aitru wr reua aay(aDG 12 (SagDIG, Galaxy Irregular Dwarf Sagittarius The o ealct asv tr odtekyt nepe nume high- interpret galaxies, to starburst key re–ionization, the including hold verse stars massive metallicity Low ABSTRACT hs ol etems ea-ormsiesaso h Local the words: of Key stars massive metal-poor hig most from the be confirmation could Pending unve these candidate. observations supergiant The red stars. massive one hosts 10- Gr SagDIG with the that accessible with firms executed analogues spectroscopy for low-resolution quest presents the in landmark tant tla content stellar ff rswr ietdto directed were orts aka ff / ; g I,ESO dIG, Sgr tr:msie–Sas al–ye–Glxe:individual: Galaxies: – early–type Stars: – massive Stars: r odiscover to ort ma. / 0Fe 10 ; O 7 err tal. et Herrero ). oe tal. et Hosek avrk .280Treó eAdz(ard,Spain (Madrid), Ardoz de Torrejón 28850 4. km jalvir hso the of chos ⊙ content n ⊙ (GTC), s content, oxygen alactic asof tars 1613, rpit1Dcme 07Cmie sn NA L MNRAS using Compiled 2017 December 1 Preprint Fe- t sis rs mic and ars. ion e ne gemn ihteNtoa cec Foundation. Science National the Astr with in agreement Research for under Universities of Association the by ated dpe o h emsOII us( 1 runs OSIRIS team’s the for adapted were VPH R2000B and slit resolution arcsec achieve 1.2 to The used Garcia). M. GTC67 (PI programme (GTC), 13A Canarias spectrograp Telescopio OSIRIS Gran the 10.4-m with the executed were observations The REDUCTION AND DATA 2 Section in provided are work future for prospects some and stars pres massive paper young Section This of galaxy. stars. observations massive spectroscopic for first explored the been never has it 2002 H confirmed determined only abundance the oxygen the populations, youngest the 2009 k,12ace eigepsr iewr eoe oteslit the g to devoted of were hours time 2 exposure of seeing total arcsec 1.2 A sky, sources. UV–bright several enclose and 12 low remarkably also is fterdgatbac RB htrnefo [Fe from range that (RGB) branch giant red the of ntecneto h oa ru nSection in Group Local the of context the in a lge ntresS- n D7t rs h aay(Fig. galaxy the cross to SD-7 and SD-1 targets on aligned was agtietfiain pcrltpsadmmesi r as are membership Section and in types spectral identification, Target − aaiso h oa nvre iaigIZ1 ( Zw18 I rivaling Universe, Local the of galaxies .5( 2.45 iraf h aawr eue olwn standard following reduced were data The ,ti au seuvln o1 to equivalent is value ( abundance this photospheric ), Solar the by Scaled ). hl aDGi osdrdoeo h xrml metal-poor extremely the of one considered is SagDIG While sdsrbtdb h ainlOtclAtooyObservato Astronomy Optical National the by distributed is ecr ta.2014 al. et Beccari 3 h icvrdrdsprin addt sdiscussed is candidate supergiant red discovered The . 2 ecie h bevtosaddt reduction. data and observations the describes ii + einSgI-I# ( SagDIG-HII#3 region esitsproa and supernovae log(O ; + R aahnsve l 1999 al. et Karachentsev log(O nTlsoi aaista con- that Canarias Telescopio an ∼ / ospoessi h atUni- past the in processes rous ls eecps hspaper This telescopes. class m -eouinflo-pstudies, follow-up h-resolution H) Group. ldtreOAtp tr and stars OBA–type three iled 00at 1000 / 1O 21 / = H) .7 ersnsa impor- an represents 7.37) = ⊙ 7.37 . aca&Hreo2013 Herrero & Garcia ⋆ ∼ aDG–Galaxies: – SagDIG 0050Å h slit The 4000–5500Å. − 0 4 . 0 13 toe ta.1991 al. et Strobel ial,summary Finally, . . 11 ih ta.2015 al. et Filho A iraf T E γ nm (AURA) onomy ( tl l v3.0 file style X sln tal. et Asplund ain tal. et Saviane rybursts. –ray 1 .Regarding ). / H] procedures = − y oper- ry, 5 . sessed . to 2.0 nthe in . from ents at h rey 1 ). ), - ) ) 2 M. Garcia

in the LMC and the SMC, and the spectral classification atlas by Gray (2000)2. They were also compared against SMC and MW standards. Membership of the programme stars with counterpart in M14’s catalogue, that registered relative shifts from multi- HST visits, was ensured by the lack of proper motions. In addi- tion, the estimated absolute magnitudes were compared against the calibrated values for the spectral types assigned in this work. The consistency of the stellar radial velocities and SagDIG’s systemic velocity (−78.5 km s−1, McConnachie 2012) was not determinant because of the low spectral resolution of these observations(∼ 300 −1 −1 km s , hence typical vrad uncertainties of ∼ 100 km s ). The following subsections provide more details on the targets that were positively identified as SagDIG massive stars.

3.1 SD-1 Figure 1. SagDIG: GTC–OSIRIS acquisition image (1.2 arcsec seeing, filter–free), and programme stars. The HST–ACS images reveal that SD-1 is the blend of two stars ∼ 1.1 arcsec apart, SD-1-blue (M14-33466; B14-8903, V=20.190) and SD-1-red (B14-8909, V=21.774). SD-1 will be identified with The spectra were flux–calibrated against the flux standard GD190, the brightest star of the blend, SD-1-blue, from now on. achieving about 5 per cent accuracy for λ >4250Å. Finally, the The Balmer series is clearly detected in its spectrum (Fig. 2), spectra were set to the heliocentric frame and corrected by radial and the comparison against the sky spectrum confirms that the velocity. The latter were calculated from the Doppler shift mea- stellar component dominates Hγ and Hβ. The data also display sured at the core of the Balmer lines. The reduced spectra are shown He i+He ii 4026Å, He i 4387Å and He i 4471Å, but the detection of in Figs. 2-3. He ii lines is unclear. He ii 4541Å is absent and the He ii 4200Å and The programme faced the challenges related to observations He ii 5411Å absorptions overlap with sky features. of faint, distant stars: contamination by nearby sources and small The exception is He ii 4686Å, detected in emission. Its width number of counts over the background level. For instance the spec- together with the absence of H i and He i emission lines ren- trum of SD-1, with V∼20.2, registered a factor ∼2.5 less counts ders a nebular origin unlikely, and the comparison against the than the extracted sky spectrum. The resulting signal to noise ratio subtracted sky spectrum discards the feature as an artefact. The (SNR) is consequently poor even after careful sky subtraction and He ii 4686Å emission could be consistent with a late–O type, with some Solar features remained. These shortcomings discourage fu- the rest of He ii lines too weak for detection at the SNR of the spec- ture optical range spectroscopic observations of targets fainter than trum. The non–detection of Mg ii 4481Å and Si iii 4552Å supports V=20 mag in grey time. the late–O classification, although it could also be caused by the Targets SD-1 to 4 are found in an area of extended contin- low metallicity of the star. uum emission (detected in Far– imaging as well, e.g. The spectrum of SD-1 is reminiscent of star C1_31 in NGC 55 Momany et al. 2014) that may contaminate the extracted spectra. (Castro et al. 2008). C1_31 is mostly featureless except for the In order to minimize blends, the sources were extracted discarding Balmer lines and a prominent broad He ii 4686Å emission, with the wings of the spatial distribution of photons in the 2-D images. He i 4471Å clearly seen but no He ii 4541Å. It was classified as early–O supergiant. Similarly, SD-1 is classified as late–O I, with the He ii 4686Å emission and narrow Balmer lines supporting 3 SOURCE IDENTIFICATION, SPECTRAL the supergiant luminosity class. The estimated absolute magnitude =− CLASSIFICATION AND MEMBERSHIP MV 6.1 agrees with the calibrated value for late–O supergiants (Massey 1998). SagDIG is projected close to the Galactic Centre and the line of sight is severely contaminated with stars. Hence, ensur- ing that the programme targets belong to the galaxy was an issue. 3.2 SD-2 This section discusses jointly target identification, spectral classifi- SD-2 is the blend of M14-31722 (F606W=20.278) and M14-31570 cation, membership, and their consistency. The location of the pro- (F606W=21.839), 0.7 arcsec apart and unresolved. SD-2 is identi- gramme stars in the colour–magnitude diagram (CMD) is evaluated fied with the brightest source M14-31722. as a final check. Its spectrum shows the Balmer series, He i 4471Å, The programme targets were identified in the Beccari et al. He i 4387Å and He i 4922Å (Fig. 2). There may be some (2014, hereafter B14) and Momany et al. (2014, M14) catalogues Mg ii 4481Å, but no lines of Si iii, Si iv or He ii. The calibrated stel- by comparing the GTC acquisition images against HST–ACS lar flux corresponds to a colder object than SD-1. Since the width archival observations (GO-9820 P.I. Y. Momany). Identification of the Balmer lines is similar to SD-1, and the strength of the He i numbers and photometry are listed in Table 1. lines is almost at noise level, SD-2 is classified as a BA-supergiant. The spectral classification of blue stars, whose most promi- The estimated absolute magnitude MV =−5.1 is consistent with a nent spectral feature is the Balmer series, was carried out following blue supergiant located at SagDIG. classical criteria (see e.g. Garcia & Herrero 2013). Their spectra are shown in Fig. 2. The rest of the sample (Fig. 3) was classi- fied after Dorda (2016)’s scheme for cool stars (G–types and later) 2 https://ned.ipac.caltech.edu/level5/Gray/frames.html

MNRAS 000, 1–5 (2017) First massive stars in SagDIG 3

Figure 2. GTC–OSIRIS spectra of blue massive stars in SagDIG. The data have been rebinned to 3 pixels for clarity, and corrected by heliocentric and . The three spectra show a clear Balmer series.

Table 1. Observed stars: spectral types (SpT) from this work, radial velocities (vrad) used in Figs. 2–3, and cross-identification numbers and photometry from M14 and B14. MV was computed from aparent magnitudes, distance modulus (SagDIG stars: DM=25.10 Momany et al. (2005); foreground stars: DM=14.50), and Cox (2000)’s intrinsic colors to estimate E V − I and correct from extinction. The last column flags SagDIG membership.

ID SpT vrad M14 F475W F606W F814W B14 VIE V − I MV Memb. [km s−1] ID ID

SD-1 Late–O I -65 33466 20.189 20.125 20.024 8903 20.190 20.154 0.51 -6.1 yes SD-2 BA I -55 31722 20.355 20.278 20.120 ‡ 20.38 20.24 0.12 -5.1 yes SD-3 OB V -95 32939 21.408 21.480 21.515 8822 21.442 21.585 0.28 -4.3 yes SD-4 G0–2 I -100 30481 21.558 20.503 19.394 8805 20.731 19.473 0.42 -5.3 yes SD-5 G4 III -90 3019 20.887 20.055 19.144 ‡ 20.38 19.23 0.26 5.3 fg SD-6 G5 V -70 † 18.59 17.82 7783 18.875 17.912 0.07 4.2 fg SD-7 K4 V -50 † 18.76 17.74 8017 19.124 17.819 0.42 3.7 fg (†) Star without M14-counterpart. F606W and F814W magnitudes estimated from B14’s V- and I-magnitudes, using Rejkuba et al. (2005)’s transformations. (‡) Star without B14-counterpart. VI-magnitudes estimated from M14, using Rejkuba et al. (2005)’s transformations.

3.3 SD-3 the stellar component, and the profile is wider than SD-1’s, suggest- ing class III or V. Since the absolute magnitude agrees with Massey Target SD-3 is faint and off–slit, resulting in extremely low spectral (1998)’s calibration for O9–B0 V (M = −4.4 to −3.8), SD-3 is SNR. The sky contribution amounts to 83 per cent of the total flux V classified as an OB-dwarf. at blue wavelengths, increasing to 91 per cent at 4800Å. SD-3 overlaps with the red component of the SD-1 blend. However, there is a shift between SD-3’s strongest features and the 3.4 SD-4 same lines in SD-1 reflecting the different radial velocities of the two stars and confirming that the lines genuinely belong to SD-3. SD-4 (M14-30481, F606W=20.503) is contaminated by a nearby The spectrum of SD-3 (Fig. 2) exhibits the Balmer series and red star with similar magnitude (M14-30025, F606W=20.360). Ac- He i 4471Å. The He ii lines, but also Si iii 4552Å and Mg ii 4481Å, cording to their proper motions, both belong to SagDIG. are absent suggesting late–O or early–B spectral type. The spec- Luminosity class was constrained in the first place, following trum lacks any He ii 4686Å, but contamination by SD-1-red begins Dorda (2016). Because the 5167Å component is the strongest of to be significant in this region. The wings of Hγ are dominated by the magnesium triplet Mg i 5167, 5172, 5184Å, SD-4 was classified

MNRAS 000, 1–5 (2017) 4 M. Garcia

Figure 4. SagDIG’s colour–magnitude diagram built from M14’s catalogue (grey dots). Marigo et al. (2017)’s isochrones for Z=0.001 Figure 3. Red programme stars. The flux–calibrated spectra were smoothed (∼1/15 Z⊙) are included for comparison, with colour coding to 5 points, and corrected by heliocentric and radial velocity. In both pan- log ageyrs=6.0,6.49,7.0,7.49,8.0. Younger isochrones are not included nels the upper spectrum has been shifted by a constant for clarity (SD-4: because they overlap at F606W−F814W∼0 in the F814W=20–22 range. 1.5E−17 erg cm−2 s−1 Å−1, SD-6: 3E−17 erg cm−2 s−1 Å−1). The isochrones were shifted by SagDIG distance modulus m − M0=25.10 and reddened by the equivalent to E B − V=0.22 (AF606W = 0.646, AF814W = 0.410, after Sirianni et al. (2005)). a supergiant. The relative strength of Fe i 5247, 5250, 5255Å vs Ca i+Fe i+Ti i 5270Å, and Mn i 5433Å vs Mn i 5447Å support this luminosity class. The presence of the G–band and its relative strength compared the Balmer decrement of SagDIG-HII#3 was adopted (Lee et al. to Hγ and Ca i 4227Å suggests that the spectral type is slightly 2003; Tang et al. 2016). later than G0. The relative strength of Fe i 4325Å compared to SD-3 is located at the main sequence where isochrones of both Hγ and Fe i 4383Å, on the other hand, suggests spectral type log ageyrs=4.5–7.0 overlap, in accordance with its OB V spec- G0 or earlier. A G0–2 I type is finally assigned. tral type. SD-2 is found between the log ageyrs=7.0 and 7.49 It is noted that the estimated absolute magnitude M =−5.3 is V isochrones, as expected considering that BA–supergiants have slighlty fainter than the calibrated M =−6.4 G0 I. Nevertheless, V evolved off the main sequence and are slightly older than O–stars. SD-4’s membership to SagDIG is supported by the lack of detec- The location of SD-4 is consistent with a RSG with similar age tion of proper motions by M14. Since red supergiants (RSG) stars log ageyrs=7.49. may experience spectral variability (see Sect. 4), SD-4 is proposed SD-1 is found between the log ageyrs=7.0–7.49 isochrones, as a RSG-candidate in SagDIG. however, it is more reddened than the other programme stars. If its E B − V=0.38 spectroscopic colour excess is taken at face value, the isochrones would need an additional shift of 3.5 Target stars in the colour–magnitude diagram ∆ F606W − F814W ∼0.2 mag and ∆F814W ∼0.3 mag to account The location of the target stars in the CMD is consistent with the for the increased extinction. This would bring SD-1 close to the spectral types assigned in this work (Fig. 4). The three OBA stars log ageyrs=6.49 isochrone (3 Myr), that better agrees with the are found in the upper part of the blue plume. The red stars SD-5, expected age of an O–type supergiant. SD-6 and SD-7 are enclosed in the vertical sequence of Galactic The unknown amount of internal reddening in SagDIG hin- foreground contamination. SD-4, on the other hand, is above the ders the interpretation of the CMD (SD-1,2,3 would be mistaken as tip of the RGB consistently with the star being a red-supergiant. intermediate−mass stars if E B − V f g was used) and may hide a The loci of the youngest stars in the CMD is also con- significant population of young and massive stars. Extinction can sistent with the theoretical isochrones (Z=0.001 or ∼1/15 Z⊙ explain the paucity of stars in the upper blue plume, which is at Marigo et al. 2017) but require a larger correction than the reported odds with the extended emission of the galaxy in the UV. These E B − V f g=0.12 foreground extinction (e.g. Momany et al. results may be considered a reminder that the total galactic stellar 2005). Internal reddening can indeed be expected, since SD-1 to content may be underestimated from the luminosity function, and 4 are located in a concentration of neutral (Young & Lo advices against selecting massive star candidates from their loca- 1997, also B14). The E B − V≃0.22 colour excess calculated from tion in classical CMDs only.

MNRAS 000, 1–5 (2017) First massive stars in SagDIG 5

4 THE MOST METAL-POOR RSG TO DATE ACKNOWLEDGMENTS

Red supergiants are He–burning descendants of Mini.30 M⊙ stars M. Garcia would like to thank F. Najarro for fruitful discussions and and the coolest stage in the life-cycle of massive stars. They are support from grants FIS2012-39162-C06-01, ESP2013-47809-C3- gathering increased interest because some RSGs (9

MNRAS 000, 1–5 (2017)