# Near-infrared spectroscopic observations of massive young stellar object candidates in the central molecular zone

G. Nandakumar, M. Schultheis, A. Feldmeier-Krause, R. Schödel, N. Neumayer, F. Matteucci, N. Ryde, A. Rojas-Arriagada, A. Tej

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

G. Nandakumar, M. Schultheis, A. Feldmeier-Krause, R. Schödel, N. Neumayer, et al., Near-infrared spectroscopic observations of massive young stellar object candidates in the central molecular zone. Astronomy and Astrophysics - A&A, 2018, 609, pp.A109. 10.1051/0004-6361/201731918. hal-02307858

HAL Id: hal-02307858

https://hal.science/hal-02307858

Submitted on 3 Sep 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Near-infrared spectroscopic observations of massive young stellar object candidates in the central molecular zone*

G. Nandakumar<fnref n="1" />, M. Schultheis<fnref n="1" />, A. Feldmeier-Krause<fnref n="2" />, R. Schödel<fnref n="3" />, N. Neumayer<fnref n="4" />, F. Matteucci<fnref n="5" />, N. Ryde<fnref n="6" />, A. Rojas-Arriagada<fnref n="7" />,<fnref n="8" />, and A. Tej<fnref n="9" />

Received 8 September 2017 / Accepted 11 November 2017

## ABSTRACT

**Context.** The central molecular zone (CMZ) is a ~200 pc region around the Galactic centre. The study of star formation in the central part of the Milky Way is of great interest as it provides a template for the closest galactic nuclei.

**Aims.** We present a spectroscopic follow-up of photometrically selected young stellar object (YSO) candidates in the CMZ of the Galactic centre. Our goal is to quantify the contamination of this YSO sample by reddened giant stars with circumstellar envelopes and to determine the star formation rate (SFR) in the CMZ.

**Methods.** We obtained KMOS low-resolution near-infrared spectra (R ~ 4000) between 2.0 and 2.5 µm of sources, many of which have been previously identified by mid-infrared photometric criteria as massive YSOs in the Galactic centre. Our final sample consists of 91 stars with good signal-to-noise ratio. We separated YSOs from cool late-type stars based on spectral features of CO and Bry at 2.3 µm and 2.16 µm, respectively. We made use of spectral energy distribution (SED) model fits to the observed photometric data points from 1.25 to 24 µm to estimate approximate masses for the YSOs.

**Results.** Using the spectroscopically identified YSOs in our sample, we confirm that existing colour-colour diagrams and colour-magnitude diagrams are unable to efficiently separate YSOs and cool late-type stars. In addition, we define a new colour-colour criterion that separates YSOs from cool late-type stars in the H-Ks vs. H-[8.0] diagram. We use this new criterion to identify YSO candidates in the |*l*| < 1°5, |*b*| < 0°5 region and use model SED fits to estimate their approximate masses. By assuming an appropriate initial mass function (IMF) and extrapolating the stellar IMF down to lower masses, we determine a SFR of ~0.046 ± 0.026 M☉ yr⁻<fnref n="1" /> assuming an average age of 0.75 ± 0.25 Myr for the YSOs. This value is lower than estimates found using the YSO counting method in the literature.

**Conclusions.** Our SFR estimate in the CMZ agrees with the previous estimates from various methods and reaffirms that star formation in the CMZ is proceeding at a lower rate than predicted by various star forming models.

**Key words.** Galaxy: center – Galaxy: stellar content – stars: formation – stars: massive – stars: pre-main sequence – Hertzsprung-Russell and C-M diagrams

## 1. Introduction

The central molecular zone (CMZ) is the innermost ~200 pc region of the Milky Way, covering about −0°7 < l < 1°8 in longitude and −0°3 < b < 0°2 in latitude. It is a giant molecular cloud complex with an asymmetric distribution of molecular clouds (see e.g. Morris & Serabyn 1996; Martin et al. 2004; Oka et al. 2005). The understanding of the physical processes occurring in the CMZ of our Galaxy is crucial for the insight in the formation or evolution of our own Milky Way.

This prodigious reservoir of molecular gas is in an active region of star formation, where there is evidence of starburst activity in the last 100 000 yr (Yusef-Zadeh et al. 2009; hereafter

YHA09). The gas pressure and temperature are higher in the CMZ than in the Galactic disk; these conditions favour a larger Jeans mass for star formation and an initial mass function (IMF) biased towards more massive stars (see Serabyn & Morris 1996; Fatuzzo & Melia 2009). Thus it is essential to understand the modes of star formation and star formation history in the unique environment of the CMZ, both to gain insight into our own Milky Way and to provide a template for circumnuclear star formation in the closest galactic nuclei.

Longmore et al. (2013a) carried out a detailed study of the variations in star formation across the Galactic plane using observational tracers of dense gas (NH₃(1, 1)) and star formation activity (masers, HII regions). These authors showed that there is a deficiency in star formation activity tracers in the CMZ given the large reservoir of dense gas available. On the other hand, they

<footnote n="1">¹ Laboratoire Lagrange, Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Blvd de l'Observatoire, 06304 Nice, France e-mail: mathias.schultheis@oca.eu</footnote>

<footnote n="2">² The Department of Astronomy and Astrophysics, The University of Chicago, 5640 S. Ellis Ave., Chicago, IL 60637, USA</footnote>

<footnote n="3">³ Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía s/n, 18008 Granada, Spain</footnote>

<footnote n="4">⁴ Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="5">⁵ Dipartimento di Fisica, Università di Trieste, 34127 Trieste, Italy</footnote>

<footnote n="6">⁶ Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 43, 221 00 Lund, Sweden</footnote>

<footnote n="7">⁷ Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Santiago, Chile</footnote>

<footnote n="8">⁸ Millennium Institute of Astrophysics, Av. Vicuña Mackenna 4860, 782-0436 Macul, Santiago, Chile</footnote>

<footnote n="9">⁹ Indian Institute of Space Science and Technology, Thiruvananthapuram 695547, India</footnote>

<footnote>* Based on observations collected at the European Southern Observatory, Chile, programme number 097.C-0208(A).</footnote>
found that various star formation models predict much higher values of star formation rates (SFR). Barnes et al. (2017) determined the average SFR across the CMZ using a variety of extra-galactic luminosity-SFR conversion and found it to be comparable to previous measurements made from young stellar object (YSO) counting and free-free emission. Thus, they ruled out systematic uncertainties in the SFR measurements as the reason for low star formation in the CMZ.

The central few parsecs of the Milky Way host a massive young population of stars. There is strong evidence that young stars in the central parsec formed in situ (see Genzel et al. 2010, and references therein). Recent spectroscopic observations have provided further strong evidence to support this. Feldmeier-Krause et al. (2015) observed the central >4 pc² of the Galactic centre and identified >100 early-type young stars by spectral classification. These authors found that early-type stars are centrally concentrated favouring the in situ formation of early-type stars. Støstad et al. (2015) mapped a smaller area within 0.28–0.92 pc from SgrA* and found a break in the distribution of young stars at 0.52 pc. These authors concluded that this break possibly indicated an outer edge to the young stellar cluster in the Galactic centre, which is expected in the case of in situ star formation.

Until recently, most studies of YSOs in the CMZ have been based on infrared photometry (Felli et al. 2002; Schuller et al. 2006; YHA09). The YSO phase of a massive star is a relatively brief phase in which the star is surrounded by dense envelopes of gas and dust (Zinnecker & Yorke 2007). These objects are best identified by their point-source infrared radiation as well as excess flux values in the mid-infrared bands. Young stellar objects are classified into three classes or stages depending on their spectral index (Lada 1987) and spectral energy distribution (SED) models (Robitaille et al. 2006). Adams et al. (1987) identified Class I YSOs as protostars with infalling envelopes, Class II YSOs as stars with disks, and Class III YSOs as those stars having the SEDs of stellar photospheres. Analogous to this classification, Robitaille et al. (2006) defined three evolutionary stages based on their derived SED model properties: Stage I objects are young protostars embedded in an opaque infalling envelope, Stage II objects are stars surrounded by an opaque disk and dispersed envelope, and Stage III objects are stars with an optically thin disk. Felli et al. (2002) searched for YSO candidates using the mid-infrared excess derived from Infrared Space Observatory (ISO) photometry at 7 and 15 µm and found a strong concentration of YSO candidates in the inner Galaxy. Schuller et al. (2006) refined the ISO mid-infrared colour criteria and argued that slightly extended mid-infrared sources were more likely to be YSOs than point-like mid-infrared sources. The authors of YHA09 identified YSO candidates with Spitzer photometry at 3.6–24 µm. Their SED fitting techniques associated most of their YSO candidates with Stage I objects; they concluded that a recent starburst took place in the CMZ.

The CMZ, however, suffers from very large and spatially variable interstellar extinction ($A_V$ = 20–40 mag; see e.g. Schultheis et al. 2009). The significant foreground extinction causes evolved stars with circumstellar envelopes, such as mass-losing asymptotic giant branch (AGB) stars, to have infrared colours similar to those of YSOs. Schultheis et al. (2003) demonstrated that near-infrared spectra are a powerful tool to distinguish YSOs from reddened AGB stars. These authors found that YSO samples in the CMZ selected by photometric colour criteria are heavily contaminated by AGB stars, red giants, and even supergiants (see their Fig. 5). By contrast, they showed that moderate-resolution spectra in the *H* and *K* bands

delineate YSOs from evolved stars by the absence or presence (respectively) of CO absorption at ~2.3 µm. Detectable YSOs at the distance of the Galactic centre (~8 kpc; see e.g. Boehle et al. 2016; Gillessen et al. 2017) are all massive, and thus never show 2.3 µm CO absorption; instead, they are featureless around 2.3 µm or show 2.3 µm CO in emission (Geballe & Persson 1987; Carr 1989; Hanson et al. 1997; Bik et al. 2006).

A recent improvement in YSO selection in the CMZ came from using Spitzer/IRS spectra to select YSOs. An et al. (2009; 2011, hereafter An11), presented Spitzer/IRS 5–35 µm spectra of 107 YSO candidates selected from 3.6–8.0 µm Spitzer photometry (Ramírez et al. 2008). The authors of An11 identified massive YSOs in the CMZ by the presence of gas-phase absorption from $C_2H_2$ (13.7 µm), HCN (14.0 µm), and $CO_2$ (15.0 µm) as well as strong and broad 15.2 µm $CO_2$ ice absorption. They found that mid-infrared spectra confirm only 33% of YSO candidates selected by their photometric criteria, and confirm 57% of YSO candidates selected photometrically by YHA09. Immer et al. (2012) analysed 5–40 µm Spitzer/IRS spectra of 57 YSO candidates selected from 7 and 15 µm ISO colours and spatial extent at 15 µm (Schuller et al. 2006). These authors identified 25% of their sources as YSOs and an additional 37% as H II regions. There is disagreement in the YSO classification even among the common sources in An11 and Immer et al. (2012) samples, suggesting uncertainties in spectroscopic YSO classification schemes as well. Koepferl et al. (2015), using radiative transfer models and realistic synthetic observations (with HYPERION, Robitaille 2011; and the FluxCompensator, Koepferl & Robitaille 2017), re-examined the YHA09 YSO sample and showed that embedded main sequence stars contaminate the YHA09 sample. These recent studies demonstrate significant contamination of photometrically selected YSO candidates samples by non-YSOs, which has important implications for CMZ star formation rates derived from photometry (YHA09).

In this paper, we present moderate-resolution 2.0–2.5 µm spectroscopic follow-up observations of a sample of 91 photometrically identified YSO candidates in the CMZ using *K*-band Multi Object Spectrograph (KMOS; Sharples et al. 2013) at VLT-UT1 (Anto). Our goal is to distinguish YSOs from evolved late-type stars by their near-infrared spectra. We discuss and show the contaminating evolved late-type stars in different colour-magnitude (CMD) and colour-colour (CCD) diagrams and define a new colour-colour criterion to distinguish these stars using our spectroscopically identified YSO sample. We estimate the SFR in the CMZ based on YSO counting and on SED fitting techniques.

## 2. Sample selection, observations, and data reduction

### 2.1. Sample selection

We selected the sample for our observation from the photometric catalogue of SIRIUS (Nishiyama et al. 2006) and point-source catalogue of Spitzer IRAC survey of the Galactic centre (Ramírez et al. 2008). The $JHK_S$ photometry from the SIRIUS catalogue has average $10\sigma$ limiting magnitudes of $J = 17.1$, $H = 16.6$ and $K_S = 15.6$ mag, while the 3.6, 4.5, 5.8, and 8.0 µm bands from the Spitzer IRAC catalogue has confusion limits of 12.4, 12.1, 11.7, and 11.2 mag, respectively. We divided our sample into three categories of high, medium, and low priorities; the highest priority (priority 1) was given to those sources in our sample that are photometrically identified
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Field distribution of our observations overlaid on the 3.6 µm Spitzer image (Stolovy et al. 2006). The white circles represent the 7/2 diameter fields that have been observed. The numbers are assigned to identify the fields, details of which are given in Table 1. Shown in red are the locations of prominent sources in the CMZ such as Sgr A*, Sgr B2, Sgr C, Quintuplet, and Arches clusters.</figcaption></figure>

YSO candidates in YHA09. We divided the rest of the sample into medium (priority 2) and low (priority 3) priorities using the following criteria that select sources exhibiting excess emission in mid-infrared regimes:

<table><thead><tr><th>Medium priority</th><th>Low priority</th></tr></thead><tbody><tr><td>[3.6] - [4.5] > 0.5</td><td>[3.6] - [8.0] > 2</td></tr><tr><td>[4.5] - [5.8] > 0.5</td><td>K<sub>S</sub> < 17</td></tr><tr><td>[5.8] - [8.0] > 1</td><td></td></tr><tr><td>K<sub>S</sub> < 17.</td><td></td></tr></tbody></table>

The KMOS consists of 24 integral field units (IFU) that can be arranged in a 7/2 diameter field per configuration, and it is crucial to prevent the 24 IFUs on the 24 pick-off arms from blocking each other. We use the KMOS ARM Allocator (KARMA), which assigns the maximum number of highest priority targets to the 24 pick-off arms, followed by medium and low priority targets thereby leaving as few arms as possible unallocated.

## 2.2. KMOS observations

Our spectroscopic observations were carried out with KMOS at VLT-UT1 (Antu) on June 23, 2016. Each of the 24 IFUs in KMOS has a field of view of 2''×8 × 2''×8. The spectral resolution of KMOS is R ~ 4300 with the wavelength range covering 1925–2500 µm. We prepared 22 fields with unique IFU configuration covering a significant part of the CMZ. But owing to bad weather conditions, only eight fields could be observed with an integration time of 900 s for each field in the nod to sky mode. The observations were carried out under photometric conditions with seeing ~0''×8. The observed field positions are shown in Fig. 1 and the field details with the number of different priority sources in each field are given in Table 1.

The fields are in regions with high stellar density as can be seen in Fig. 1. Hence, it was not possible to observe sky in each IFU by dithering or nodding to a new position within the respective field. So we observed a dark cloud G359.94+0.17

($\alpha$ ~ 266°2, $\delta$ ~ -28°9; Dutra & Bica 2001) in free dither mode with 2''×8 dither in between two 900 s exposures to carry out a proper sky subtraction. The sky observation was carried out after every two field observations. B and A type stars were observed for telluric corrections after every sky offset.

## 2.3. Data reduction

We used the ESO KMOS Recipe Flexible Execution Workbench (Reflex; Freudling et al. 2013) for data reduction. It organises the science and associated sky and calibration data together based on the calibration source type and their proximity in time to science observations. This is followed by dark level correction of frames, flat fielding, wavelength calibration, spatial illumination correction, telluric correction, sky subtraction, and cube reconstruction of the science data by dedicated pipeline recipes (or stages).

We made use of an IDL routine to remove the Bry absorption line from each telluric spectrum by fitting the Bry line with a Lorentz profile and subtracting it from the telluric spectrum. This routine also removes the stellar continuum by dividing it by a blackbody spectrum. We also removed cosmic rays from the final reconstructed object cube with a 3D version of L.A. Cosmic (van Dokkum 2001).

We extracted spectra from 173 data cubes using *kmos_extract_spec* recipe with the ESO Recipe Execution Tool (EsoRex). The *kmos_extract_spec* recipe extracts a spectrum from a data cube with the option of defining a mask manually or automatically by fitting a normalised profile to the image of the data cube. We identified multiple sources in 53 data cubes and extracted their spectra by defining the mask manually. Finally, we extracted nearly 250 spectra and used IRAF median filtering with seven pixels as filter size to smoothen the spectra. After discarding spectra with a signal-to-noise ratio (S/N) below 20 or having negative flux values, there were 91 spectra left in our sample with good S/N. Among these spectra, there are 15 from data cubes with multiple sources, and we selected the spectrum of the brightest source in the data cube.
<table><caption>**Table 1.** Details of the observed fields.</caption><thead><tr><th>Field No.</th><th><i>l</i> (°)</th><th><i>b</i> (°)</th><th>Priority 1</th><th>Priority 2</th><th>Priority 3</th><th>Sky</th><th>Random source</th></tr></thead><tbody><tr><td>1</td><td>359.2568</td><td>-0.0813</td><td>13</td><td>3</td><td>3</td><td>3</td><td>2</td></tr><tr><td>4</td><td>359.5854</td><td>-0.0351</td><td>13</td><td>5</td><td>4</td><td>1</td><td>1</td></tr><tr><td>5</td><td>359.5131</td><td>-0.0986</td><td>9</td><td>3</td><td>11</td><td>1</td><td>0</td></tr><tr><td>6</td><td>359.8000</td><td>-0.0500</td><td>15</td><td>1</td><td>7</td><td>1</td><td>0</td></tr><tr><td>8</td><td>0.0419</td><td>-0.0547</td><td>0</td><td>5</td><td>18</td><td>1</td><td>0</td></tr><tr><td>14</td><td>0.6674</td><td>-0.0527</td><td>12</td><td>9</td><td>3</td><td>1</td><td>0</td></tr><tr><td>15</td><td>0.0319</td><td>0.0506</td><td>3</td><td>9</td><td>12</td><td>2</td><td>0</td></tr><tr><td>18</td><td>359.1303</td><td>-0.0080</td><td>12</td><td>1</td><td>5</td><td>4</td><td>2</td></tr></tbody></table>

Notes. The (*l*, *b*) of field centres and the number of high, medium, and low priority sources, and the manually allocated sky and random sources for the free arms are listed.

### 3. Classification

In this section, we classify our targets as YSOs and late-type stars based on their spectra. Using this classification, we evaluate different photometric YSO classification criteria and suggest a new criterion to distinguish YSOs and late-type stars.

#### 3.1. Spectroscopic classification

We classified our spectra mainly based on the presence or absence of the ¹²CO (2, 0)-band at 2.3 µm. CO absorptions bands are typically found in late-type G, K, M giants and AGB stars. In addition, we also used Bry, found in emission, absorption, or with a P Cygni-type profile in massive YSOs (Cooper et al. 2013). We carried out background subtraction during data reduction to make sure that the Bry emission lines are intrinsic to the source. Still, we expect contamination from the OB main sequence or post-main sequence Wolf Rayet stars, some of whose spectra also show Bry in emission attributed to their stellar wind (Mauerhan et al. 2010a). CO band emission at 2.3 µm is also considered to be an indication for the presence of a dense circumstellar disk and hence a YSO signature (Bik et al. 2006). Some spectra show a featureless continuum at 2.0–2.5 µm; these could be either YSOs (Greene & Lada 1996) or dusty late-type carbon-type (WC) Wolf-Rayet stars (Mauerhan et al. 2010a).

We measured the equivalent width (EW) of the ¹²CO (2, 0) band at 2.3 µm using the same CO band and continuum points as in Ramírez et al. (2000). In addition, we measured the EW of Bry line at 2.16 µm. Figure 2 shows the EW(CO) vs. EW(Bry) plot of our 91 sources. We find that there are two separate groups of stars with a very evident gap, which we approximately denote by the dashed line at EW(CO) = 10 Å. Positive values of EW indicate that the line is in absorption while negative values indicate it is in emission. Thus we classify the stars to the right of the dashed line as cool, late-type stars and those to the left as YSOs. All cool, late-type stars (represented by filled black circles) lie very close to EW(Bry) = 0 Å indicating the absence of this particular feature in their spectra. The majority of stars we classify as YSOs also show no Bry feature, while approximately five stars show Bry in emission (EW(Bry) < 1 Å) and only one star shows Bry in absorption (EW(Bry) > 1 Å), which also shows CO in emission (EW(CO) = -10 Å).

Based on the above-mentioned classification scheme, there are 23 spectroscopically identified YSOs in our sample. Figures 3 and 4 show the reduced normalised spectra of YSOs and cool late-type stars, respectively. We searched for previously identified sources in the SIMBAD database with a search

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Equivalent widths measured for the ¹²CO (2, 0) line and Bry line for YSOs (blue stars) and late-type stars (black dots). We separate the two populations approximately using the dashed line at EW(Bry) = 10 Å. Different symbols and colours represent the classification of the SIMBAD matches to our sources by searching within 2′′0 (see Sect. 3.1).</figcaption></figure>

radius of 2′′0 from each of the 23 YSOs to check their status in the literature. Table 2 lists all 23 sources with their Spitzer Space Telescope (SST) Galactic centre (GC) number, equatorial coordinates in degrees, distance of the SIMBAD match from the source, source type along with corresponding references, and JHK_S magnitudes from SIRIUS catalogue. There are seven sources in common with An11, out of which only one is confirmed to be YSO and one is considered to be a possible YSO by An11. The remaining five sources have been classified as non-YSOs, with #517724 also classified as an OB super giant star in Mauerhan et al. (2010b) based on absorption lines of Bry at 2.1661 µm, NIII at 2.115 µm, and He I at 2.058, 2.113, and 2.1647 µm. Two other sources (#528828 and #599826) have been classified as possible long period variable stars in Matsunaga et al. (2009) and Glass et al. (2001) based on periods estimated using near-infrared observations, although no clear periodicity was found for them. Two sources (#531300 and #584563) have been classified as blue super giant stars based on weak Bry emission or absorption features, NIII and CIV contributions, and HeI absorption profiles at 2.058 µm in addition to detected Paschen-α (Pα) excess (Mauerhan et al. 2010a). The counterpart to the source #238110 has been classified as X-ray source in Muno et al. (2003). Three sources (#358063, #395315, and #520760) are in common with the YHA09 sample, out of which #520760 is classified also as a radio source in
<table><caption>**Table 2.** Details of the search for previously identified sources within 2′′0 of our spectroscopically identified YSOs in the SIMBAD database.</caption><thead><tr><th>SST GC No.</th><th>RA (°)</th><th>Dec (°)</th><th>Distance (in ′′)</th><th>Source type</th><th>J (mag)</th><th>H (mag)</th><th>K<sub>S</sub> (mag)</th><th>Field</th></tr></thead><tbody><tr><td>238110</td><td>265.93584</td><td>-29.67287</td><td>0.64</td><td>X-ray source<sup>a</sup></td><td>17.19</td><td>13.64</td><td>11.84</td><td>18</td></tr><tr><td>517724</td><td>266.40542</td><td>-28.89823</td><td>0.07</td><td>OB supergiant<sup>b</sup>, non-YSO<sup>c</sup></td><td>15.75</td><td>12.80</td><td>11.23</td><td>15</td></tr><tr><td>520760</td><td>266.41002</td><td>-28.89086</td><td>1.54</td><td>radio source<sup>d</sup>, YSO<sup>e</sup></td><td>...</td><td>15.59</td><td>13.17</td><td>15</td></tr><tr><td>524419</td><td>266.41569</td><td>-28.89559</td><td>0.50</td><td>non-YSO<sup>e</sup></td><td>...</td><td>15.86</td><td>14.02</td><td>15</td></tr><tr><td>525666</td><td>266.41759</td><td>-28.89117</td><td>0.04</td><td>non-YSO<sup>e</sup></td><td>...</td><td>14.49</td><td>12.74</td><td>15</td></tr><tr><td>531300</td><td>266.42633</td><td>-28.87979</td><td>0.17</td><td>Blue super giant<sup>b</sup></td><td>14.70</td><td>11.67</td><td>10.11</td><td>15</td></tr><tr><td>535007</td><td>266.43185</td><td>-28.87358</td><td>0.05</td><td>non-YSO<sup>e</sup></td><td>...</td><td>16.05</td><td>13.81</td><td>15</td></tr><tr><td>528828</td><td>266.42248</td><td>-28.90728</td><td>0.11</td><td>Long period variable star<sup>f</sup></td><td>...</td><td>...</td><td>13.69</td><td>8</td></tr><tr><td>541457</td><td>266.44147</td><td>-28.90716</td><td>...</td><td>...</td><td>...</td><td>14.28</td><td>12.53</td><td>8</td></tr><tr><td>553700</td><td>266.45999</td><td>-28.91312</td><td>...</td><td>...</td><td>15.11</td><td>13.09</td><td>12.27</td><td>8</td></tr><tr><td>563727</td><td>266.47539</td><td>-28.97628</td><td>...</td><td>...</td><td>...</td><td>14.75</td><td>12.97</td><td>8</td></tr><tr><td>567598</td><td>266.48123</td><td>-28.93786</td><td>...</td><td>...</td><td>...</td><td>14.71</td><td>12.92</td><td>8</td></tr><tr><td>609669</td><td>266.54470</td><td>-28.92073</td><td>...</td><td>...</td><td>...</td><td>16.24</td><td>14.17</td><td>8</td></tr><tr><td>610642</td><td>266.54618</td><td>-28.92802</td><td>0.05</td><td>Maybe YSO<sup>c</sup></td><td>...</td><td>15.26</td><td>12.69</td><td>8</td></tr><tr><td>584563</td><td>266.50690</td><td>-28.92095</td><td>0.17</td><td>Blue super giant<sup>b</sup></td><td>13.45</td><td>10.69</td><td>9.12</td><td>8</td></tr><tr><td>585974</td><td>266.50900</td><td>-28.95653</td><td>...</td><td>...</td><td>16.80</td><td>11.73</td><td>9.42</td><td>8</td></tr><tr><td>599826</td><td>266.52982</td><td>-28.93144</td><td>0.52</td><td>Long period variable star<sup>f</sup></td><td>...</td><td>...</td><td>13.87</td><td>8</td></tr><tr><td>373107</td><td>266.17890</td><td>-29.33702</td><td>...</td><td>...</td><td>16.48</td><td>13.59</td><td>11.94</td><td>4</td></tr><tr><td>388790</td><td>266.20389</td><td>-29.39521</td><td>0.40</td><td>non-YSO<sup>e</sup></td><td>...</td><td>14.12</td><td>12.05</td><td>5</td></tr><tr><td>352034</td><td>266.14518</td><td>-29.39360</td><td>...</td><td>...</td><td>...</td><td>12.52</td><td>10.95</td><td>5</td></tr><tr><td>395315</td><td>266.21438</td><td>-29.34097</td><td>0.02</td><td>YSO<sup>e</sup></td><td>...</td><td>15.95</td><td>13.65</td><td>4</td></tr><tr><td>405235</td><td>266.23023</td><td>-29.26057</td><td>0.28</td><td>non-YSO<sup>e</sup></td><td>...</td><td>14.77</td><td>12.97</td><td>4</td></tr><tr><td>358063</td><td>266.15475</td><td>-29.31017</td><td>0.14</td><td>YSO<sup>e</sup></td><td>15.59</td><td>13.42</td><td>11.79</td><td>4</td></tr></tbody></table>

**Notes.** For each source represented by their SST GC number, their equatorial coordinates in degrees, distance from the source, source type along with corresponding references, and $JHK_S$ magnitudes from SIRIUS catalogue are listed. For sources with no counterpart within 2′′0 and no valid photometry in J or H bands, we use “...” for the corresponding column. (a) Muno et al. (2003); (b) Mauerhan et al. (2010a); (c) An et al. (2011); (d) Dong et al. (2015); (e) Yusef-Zadeh et al. (2009); (f) Matsunaga et al. (2009), Glass et al. (2001).

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Normalised spectra of YSOs classified based on the absence of $^{12}$CO (2, 0) band absorption line and presence of Bry emission or absorption line. Dashed lines and shaded areas represent the central wavelengths and range of continuum used for measuring equivalent widths of Bry, $^{12}$CO (2, 0), $^{12}$CO (3, 1), and $^{12}$CO (4, 2) bands. The SST GC No. for each source is specified adjacent to the corresponding spectra.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Same as Fig. 2 but for cool late-type stars in fields 8 and 15, classified based on the presence of $^{12}$CO (2, 0) band absorption line.</figcaption></figure>

Dong et al. (2015). The remaining eight sources do not have any counterparts in the SIMBAD database within 2′′0.

In Fig. 2, we show the sources classified in the literature as X-ray sources, OB or blue supergiants, and non-YSOs using separate symbols. Regarding the classification of five sources as

non-YSOs by An11, the large pixel sizes of *Spitzer*/IRS spectra can lead to misidentification of sources in high stellar density regions such as in the CMZ. Also one non-YSO (#405235) shows Bry in absorption and CO in emission, indicating the presence of a dense circumstellar disk and considered to be a massive YSO feature though rarely seen (Bik et al. 2006; Cooper et al. 2013). Thus we stick with our classification scheme for this source. Since no clear periodicity was found for the two sources classified as long period variable stars, we assume that these are
<figure><img src="image_7.png" /><figcaption>Fig. 5. Panel a: [5.8] – [8.0] vs. [3.6] – [4.5] diagram used to identify different classes of YSOs based on the disk and envelope models of low mass YSOs as shown in Allen et al. (2004) and Megeath et al. (2004). Class II YSOs are expected to be concentrated in the small box, while Class I YSOs in the bigger box. Panel b: Same colour-colour diagram as (a) with the region shown by black polygon where Stage I YSOs are expected to lie (Robitaille et al. 2006). Panel c: [8.0] – [24] vs. [24] diagram showing the criteria used by YHA09 to choose their sample of possible YSO candidates (region to the right of the dashed line). The black arrows in (a) and (b) represent the extinction vector estimated for $A_K = 2$ mag (typical of the CMZ) using the $\frac{A_1}{A_K}$ relations from Nishiyama et al. (2009). In each diagram, the spectroscopically identified YSOs are shown using blue star symbols and the cool late-type stars using black filled circles. Other symbols and colours represent the classification of the SIMBAD matches to our sources by searching within $2''/0$ (see Sect. 3.1).</figcaption></figure>

YSOs as well. Two blue supergiants (#584563 and #531300) show Bry in emission, while no clear emission is seen for #517724 classified as OB supergiant. We made an approximate EW measurement of NIII at 2.115 µm and found that all three sources mentioned above as well as the radio source #238110 show a clear emission feature that is not seen in the rest of our YSOs. Thus we consider their classification as supergiants to be acceptable. Based on these measurements we conclude that none of the rest of our YSOs are O/B supergiants.

## 3.2. Classification using photometric criteria

Several previous studies have made use of colour-colour diagrams (CCDs) to define criteria to classify YSOs. We are trying with our spectroscopic sample of YSOs and non-YSOs to establish new and more reliable photometric criteria to distinguish YSOs from non-YSOs. Initially, we made use of some YSO classification criteria using CCDs implemented in the literature to check whether they are able to classify and separate our sample of spectroscopically identified YSOs from late-type stars. For this, we obtained the photometry of point sources at 3.6 µm, 4.5 µm, 5.8 µm, 8.0 µm (Ramírez et al. 2008), and 24 µm from the MIPSGAL catalogue (Gutermuth & Heyer 2015), and we merged these with our data by searching within radii of $2''/0$. We find 87 sources with valid 3.6 to 8.0 µm photometry and 28 sources with valid 24 µm photometry in our sample.

Figure 5 shows two CCDs and a colour-magnitude diagram (CMD) with respective criteria from the literature to classify YSOs. In Fig. 5a, we plot [5.8] – [8.0] vs. [3.6] – [4.5] with a small and big box representing the regions belonging to Class II and Class I YSOs, respectively. Those areas come from the disk and envelope models of low mass YSOs as shown in Allen et al. (2004) and Megeath et al. (2004). Figure 5b shows the same plot, but with a polygon used to define the region enclosing Stage I YSOs as defined in Robitaille et al. (2006). YHA09 used the CMD ([8.0] – [24] vs. [24]) in Fig. 5c to choose their sample of possible YSO candidates in the CMZ by considering all sources lying to the right of the dashed line as YSOs. We estimate and show the extinction vector (black arrow) for the two CCDs, assuming an $A_K$ of 2 mag (corresponding to $A_V$ = 30 mag; typical of the CMZ from Schultheis et al. 2009) using the $\frac{A_1}{A_K}$ relations from Nishiyama et al. (2009). It is clear from each plot that, even after taking the extinction into account, there is severe contamination from the late-type stars in the regions defined to

<figure><img src="image_8.png" /><figcaption>Fig. 6. $H - K_S$ vs. $H - [8.0]$ diagram. This plot shows a clear trend for cool late-type stars separating the YSOs in our sample. The dashed line represents our proposed criterion to separate YSOs and late-type stars. The black arrow represents the extinction vector estimated for $A_K = 2$ mag (typical of the CMZ) using the $\frac{A_1}{A_K}$ relations from Nishiyama et al. (2009).</figcaption></figure>

contain YSOs and that there is no clean colour-colour criterion visible.

Hence, we tested various combinations of colours and magnitudes to clearly separate YSOs from cool late-type stars in our sample. We see a clear linear trend followed by the late-type stars only in the CCD $H - K_S$ vs. $H - [8.0]$, as shown in Fig. 6, while the YSOs exhibit redder $H - [8.0]$ colours and are thus clearly separated. We defined a rough criterion to separate YSOs from cool late-type stars, taking into account the $H - K_S$ cut at 1.5 estimated to remove foreground sources (assuming $A_V$ = 30 mag and using extinction laws of Nishiyama et al. 2009). A similar $H - K_S$ cut was also suggested by Schödel et al. (2010) to remove foreground sources based on their study towards the central part of the Galaxy. We estimate the extinction vector again, as mentioned before, and show that it is almost parallel to the line separating YSOs from late-type stars. Thus the extinction does not greatly affect our criterion. We want to stress that our sample is small and this criterion needs to be confirmed by a larger sample,

$$ (H - [8.0]) = 2.75 \times (H - K_S) + 1.75; \quad 1.5 < (H - K_S) \le 5. \quad (1) $$
# 4. Mass estimates using SED model fits

In the above section, we spectroscopically identified 23 YSOs from among 91 sources in the CMZ. In this section, we construct their near- and mid-infrared SEDs and fit these sources using synthetic SED models for YSOs to constrain their stellar parameters (e.g. stellar radius $R_\star$, effective temperature $T_{\text{eff}}$, stellar mass $M_\star$, total stellar luminosity $L_\star$, and extinction in V-band $A_V$).

## 4.1. Robitaille et al. models

Robitaille et al. (2006) presented a set of approximately 20 000 radiative transfer models with corresponding SEDs assuming an accretion scenario with a central star surrounded by an accretion disk, infalling envelope, and bipolar cavities, i.e. YSOs. Robitaille et al. (2007) presented a tool to fit these YSO model SEDs to observations, providing a range in the parameter (e.g. stellar mass, total luminosity, extinction in V-band, envelope accretion rate, and age) space corresponding to a set of best fit models.

These models are largely in use to estimate approximate values of stellar parameters for a photometrically or spectroscopically identified sample of YSOs. YHA09 and An11 used these models to classify YSOs into various evolutionary stages based on the envelope infall rate and disk accretion rate of each source. YHA09 in turn estimate the SFR in the CMZ using masses they obtain from the SED fits.

Recently, Robitaille (2017), hereafter R17, introduced an improved set of SED models for YSOs that covers a much wider range of parameter space and excluding most of model-dependent parameters in addition to several other improvements. Unlike previous models, there are 18 different sets of models with increasing complexity that vary from a single central star to a star in an ambient medium surrounded by accreting disk, infalling envelope, and bipolar cavities as described in detail in R17. We used the latest R17 models to fit the SEDs of sources in our sample to estimate stellar parameter such as stellar radius, luminosity, or effective temperature. We used these parameters to determine the stellar masses of our YSO sample.

## 4.2. SED fits

We construct the SEDs for our sample using wavelengths ranging from 1.25–24 µm. As mentioned in Sect. 3.2, we used the *JHK*<sub>S</sub> photometry from the SIRIUS catalogue, 3.6 to 8.0 µm photometry from Ramírez et al. (2008) and 24 µm photometry from Gutermuth & Heyer (2015). In addition, we use the 15 µm photometry from the ISOGAL point source catalogue (Omont et al. 2003a,b) so that we constrain the SEDs over a large wavelength range. Searching within 2′/0 of YSOs in our sample, we find seven sources in the ISOGAL PSC out of which only two have valid 15 µm magnitudes. Within the same search radius, we find six YSOs with a match in the 24 µm catalogue, all of which have valid photometry.

Thus, among 23 spectroscopically identified YSOs in our sample, in addition to 1.25–8.0 µm magnitudes, six sources have only 24 µm magnitudes, two sources have only 15 µm magnitudes, and there are no sources with valid magnitude determined at both 15 µm as well as 24 µm. We find that the SED fitting by the model requires data points at $\lambda > 12$ µm to give reliable results. For that reason, we carried out SED fits using the SED fitting tool only for these eight sources using the above-mentioned set of photometry.

We assume the source distance to be in the range 7 kpc < $R < 9$ kpc from the Sun and interstellar extinction along the line of sight to the Galactic centre to be in the range 20 mag < $A_V < 50$ mag (Schultheis et al. 2009), ensuring that these sources belong to the CMZ. These assumptions ensure that the conditions at the Galactic centre are considered while fitting the model SEDs to our observed SEDs. We assume typical errors of 0.05 mag for *JHK*<sub>S</sub> photometry and 0.1 mag for 3.6 to 8.0 µm photometry, while ISOGAL and MIPSGAL catalogues provide typical errors of ~0.05 mag for 15 µm and 0.1 mag for 24 µm photometry, respectively. To include reasonable fitting results, we selected all SEDs that satisfy $\chi^2 - \chi^2_{\text{best}} < 5$ per data points for each source in all 18 model sets. For each model set, $\chi^2_{\text{best}}$ represents $\chi^2$ value of the best fit. Thus each source SED is fitted with 18 different model sets, each of which gives a best fit SED with a $\chi^2_{\text{best}}$ value. For each source, we selected the model set corresponding to the best fit SED with the lowest $\chi^2_{\text{best}}$ value as the one that best represents the evolutionary stage of the source. Figure 7 shows typical examples of SED fitting results for eight YSOs in our sample.

## 4.3. Fit parameters and mass

We chose the model set corresponding to the best fit SED with the lowest $\chi^2_{\text{best}}$ value as mentioned above and estimate mean values of $A_V$, $T_{\text{eff}}$, and stellar radius, $R_\star$ from all the fits satisfying the $\chi^2$ cut in the chosen model set. We estimated approximate values for the stellar luminosity, $L_\star$, using the Stefan-Boltzmann law from $T_{\text{eff}}$ and $R_\star$, assuming solar $T_{\text{eff}}$ to be 5772 K. To determine an approximate mass for each YSO, we used the pre-main sequence (PMS) tracks for stars with metallicity of $Z = 0.02$ and mass range of 0.8 $M_\odot$ to 60 $M_\odot$ from Bernasconi & Maeder (1996). The masses were sampled in a non-uniform manner with stellar tracks provided for 0.8, 1.0, 1.5, 2.0, 3.0, 5.0, 9.0, 15.0, 25.0, and 60.0 $M_\odot$.

We calculated the separation of each source from the stellar track for each mass in the $\log_{10}(L_\star) - \log_{10}(T_{\text{eff}})$ space, and assigned each source the mass corresponding to the track at the least separation. This exercise was carried out for all fits that satisfy the $\chi^2$ cut in the chosen model set and we estimated a mean mass from them. The standard deviation in mass from all SED fits can be used to make a rough estimate of the error. A zero error for mass is obtained when there is only one SED fit or when there is a single closest stellar track to source positions from all SED fits of the chosen model set. We assume their mass uncertainties to be limited by the mass sampling of PMS tracks. Figure 8 shows the $\log_{10}(L_\star)$ vs. $\log_{10}(T_{\text{eff}})$ plot with the stellar tracks and location of 8 spectroscopically identified YSOs. Mass estimates range from ~8 to 20 $M_\odot$, as expected for high mass YSOs.

We estimated the uncertainties for $T_{\text{eff}}$, $L_\star$, and $A_V$ similarly from standard deviation in their values from all SED fits of the chosen model set. Table 3 lists the main model fit parameters and the estimated masses for the YSOs. The $A_V$ values estimated by the SED models are mostly close to the lower end of our constraints (20 mag < $A_V < 50$ mag), which is not the expected case. So we used the extinction map of Schultheis et al. (2009) to estimate the foreground visual extinction close to the location of our sample, by searching within the radius corresponding to the pixel size of the extinction map (2'). The estimated $A_{V\text{map}}$ values are listed in the last column of Table 3. We find significant difference in $A_V$ from models and $A_{V\text{map}}$ from the extinction map (mean difference ~9.3 mag), suggesting that the models
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** SED fits for 8 YSOs in our sample with available photometric data up to 24 µm, performed by the SED fitter tool using different YSO SED models in R17. The best fit is shown using solid black line and grey lines show all other fits that satisfy the criteria: $\chi^2 - \chi^2_{best} < 5$ per data point.</figcaption></figure>

need to be improved. A similar disagreement between $A_V$ from Robitaille et al. (2007) models and $A_V$ from Schultheis et al. (2009) extinction map was estimated by An11 for their spectroscopically identified YSOs.

**5. Star formation rate in the CMZ**

One of the commonly used ways to estimate the SFR in the CMZ is the YSO counting method. In this method, masses of photometrically or spectroscopically confirmed YSOs in the region are estimated either from SED fits or, using the zero-age main sequence (ZAMS) luminosity-mass relations. By assuming an appropriate IMF and extrapolating the stellar IMF down to lower masses, the total embedded stellar population mass of the region can then be estimated. Owing to the low number of spectroscopically identified YSOs, it is not possible to apply this method to our spectroscopic sample.

However it is possible to use our photometric selection criterion (see Sect. 3.2) based on the $H - K_S$ vs. $H$-[8.0] diagram to obtain a much more complete sample of YSOs in the CMZ. For this, we used the photometric catalogue of SIRIUS towards the Galactic centre from which we selected our observed sample. We combined this sample with 3.6–8.0 µm photometry from Ramírez et al. (2008), 24 µm photometry from Gutermuth & Heyer (2015), and 15 µm photometry from ISO-GAL PSC. Within $|l| < 1^\circ.5$ and $|b| < 0^\circ.5$ we find 16 180 sources

<figure><img src="image_10.png" /><figcaption>**Fig. 8.** $\log_{10}(L_\star)$ vs. $\log_{10}(T_{\text{eff}})$ diagram showing the location of YSOs (blue stars) and the pre-main sequence (PMS) stellar tracks from Bernasconi & Maeder (1996) for different masses in black. The errors in $T_{\text{eff}}$ and $L_\star$ are also shown.</figcaption></figure>

with valid photometric magnitudes in H, $K_S$, and 8.0 µm bands and in either of the two bands: 15 µm or 24 µm. We then applied our criterion (see Eq. (1)), identifying 334 sources as YSOs. As
<table><caption>**Table 3.** Details of the SED fit parameters.</caption><thead><tr><th>SST GC No.</th><th>Model</th><th>N<sub>data</sub></th><th>N<sub>fits</sub></th><th>χ<sup>2</sup><sub>best</sub></th><th>⟨A<sub>V</sub>⟩ (mag)</th><th>⟨log<sub>10</sub>(L<sub>⋆</sub>)⟩(L<sub>☉</sub>)</th><th>⟨T<sub>eff</sub>⟩ (K)</th><th>⟨M<sub>⋆</sub>⟩(M<sub>☉</sub>)</th><th>A<sub>Vmap</sub> (mag)</th></tr></thead><tbody><tr><td>238110</td><td>spubhmi<sup>a</sup></td><td>8</td><td>2</td><td>96.7</td><td>20.5 ± 0.4</td><td>4.4 ± 0.0</td><td>8329 ± 936</td><td>20.0 ± 5.0</td><td>25.5 ± 1.1</td></tr><tr><td>352034</td><td>sp-s-i<sup>b</sup></td><td>7</td><td>329</td><td>1.1</td><td>22.8 ± 2.4</td><td>3.7 ± 0.4</td><td>6668 ± 2213</td><td>12.4 ± 3.8</td><td>45.0 ± 11.3</td></tr><tr><td>358063</td><td>sp-s-i<sup>b</sup></td><td>8</td><td>45</td><td>30.1</td><td>20.3 ± 0.3</td><td>3.7 ± 0.1</td><td>9091 ± 979</td><td>9.4 ± 1.5</td><td>32.3 ± 5.8</td></tr><tr><td>373107</td><td>sp-h-i<sup>c</sup></td><td>8</td><td>9</td><td>104.5</td><td>24.8 ± 0.6</td><td>4.1 ± 0.0</td><td>10524 ± 499</td><td>15.0*</td><td>31.7 ± 5.1</td></tr><tr><td>388790</td><td>spubsmi<sup>d</sup></td><td>7</td><td>13</td><td>0.9</td><td>22.2 ± 1.2</td><td>3.7 ± 0.2</td><td>6946 ± 1255</td><td>12.7 ± 2.9</td><td>43.0 ± 12.4</td></tr><tr><td>395315</td><td>spubsmi<sup>d</sup></td><td>5</td><td>91</td><td>1.9</td><td>24.2 ± 3.5</td><td>3.2 ± 0.2</td><td>9397 ± 3981</td><td>7.8 ± 2.0</td><td>34.3 ± 4.2</td></tr><tr><td>405235</td><td>spu-smi<sup>e</sup></td><td>7</td><td>1</td><td>32.7</td><td>20.0</td><td>3.9</td><td>6671</td><td>15.0<sup>+</sup></td><td>37.5 ± 2.5</td></tr><tr><td>609669</td><td>spubsmi<sup>d</sup></td><td>7</td><td>132</td><td>8.6</td><td>22.8 ± 2.7</td><td>3.3 ± 0.2</td><td>9932 ± 3436</td><td>8.2 ± 1.7</td><td>26.7 ± 2.6</td></tr></tbody></table>

**Notes.** For each source represented by their SST GC No., average values of parameters from the model set corresponding to the best fit SED with the lowest $\chi^2_{best}$ value are chosen and listed. The errors or uncertainties for parameters are roughly determined from the standard deviation in mass from all fits. The last column shows the $A_V$ values from the extinction map of Schultheis et al. (2009) within 2' of each source. (a) A complex model in which the central star with a disk, a variable disk inner radius, and bipolar cavities are enclosed in a rotationally flattened envelope structure surrounded by ambient interstellar medium. (b) Disks around a central star with a non-variable inner radius. No surrounding envelope or ambient interstellar medium. (c) Same as b except that the disk inner radius is variable. (d) Same as a except that the disk inner radius is set to the dust sublimation radius. (e) Same as d except that there are no bipolar cavities. (e) Only one closest stellar track to source position (see Fig. 8) from all SED fits of the chosen model set. (+) Only one SED fit to the observed photometry by the chosen model set.

<figure><img src="image_11.png" /></figure>

**Fig. 9.** Plots showing candidate YSOs in the CMZ selected using our criterion and their distribution in the (l, b) plane. Left: $H - K_S$ vs. $H - [8.0]$ diagram is used to select YSOs from the source sample constructed by combining NIR catalogue of Galactic centre using SIRIUS (Nagashima et al. 1999; Nagayama et al. 2003), 3.6–8.0 µm photometry from Ramírez et al. (2008), 24 µm photometry from Gutermuth & Heyer (2015) and 15 µm photometry from ISOGAL PSC. The dashed line represents our proposed criteria to separate YSOs and late-type stars. Right: spatial distribution of YSOs in the (l, b) plane colour coded with the number of YSOs in (l, b) bins of 0:05 each. The white patch close to (l, b) ~ (0:0, 0:0) is the observational artefact from the Gutermuth & Heyer (2015) and the ISOGAL PSC catalogues where data are lacking.

seen in Sect. 3.2, foreground sources were removed by default using this criterion. However, OB supergiants can still contaminate our YSO sample.

Figure 9 shows the selected YSOs in $H - K_S$ vs. $H - [8.0]$ diagram (left panel) and spatial distribution of YSOs in the (l, b) plane colour coded with the number of YSOs in (l, b) bins of 0:05 each. The white patch close to ~0:0 latitude and longitude is an observational artefact from the Gutermuth & Heyer (2015) and the ISOGAL PSC catalogues where data are lacking. As a result, the YSO count is higher at the negative longitudes than in positive longitudes, contrary to the fact that two-thirds of molecular gas is on positive longitudes (Bania 1977; Bally et al. 1988; Morris & Serabyn 1996; Oka et al. 2005, etc.). We performed the SED fitting using R17 models for these 334 sources and we determined approximate masses for these sources as mentioned in Sect. 4.3.

We chose ~190 sources with $\chi^2_{best} < 35$ (chosen based on the average value of $\chi^2_{best}$ among the 8 spectroscopically confirmed YSOs) to plot the mass distribution, which ranges from 2.7 $M_⊙$ to 35 $M_⊙$. The distribution peaks at ~8 $M_⊙$, emphasizing that the majority of YSOs are in the high mass range. Thus, we miss the low mass stars, and hence we adopted the Kroupa IMF (Kroupa 2001) to fit it to the peak of our distribution, extrapolated it to lower masses, and estimated the total embedded stellar population in the CMZ. The Kroupa IMF for different mass ranges is defined as follows:

$$ \zeta(M) = AM^{-2.3} \quad \text{for } 0.5 M_\odot \le M \le 120 M_\odot, \qquad (2) $$

$$ \zeta(M) = BM^{-1.3} \quad \text{for } 0.08 M_\odot \le M \le 0.5 M_\odot, \qquad (3) $$

$$ \zeta(M) = CM^{-0.3} \quad \text{for } 0.01 M_\odot \le M \le 0.08 M_\odot, \qquad (4) $$
<figure><img src="image_12.png" /><figcaption>Fig. 10. Mass distribution histogram (filled blue bars) of YSOs. These YSOs were selected using our new photometric colour–colour criteria from the SIRIUS catalogue within |l| < 1°5 and |b| < 0°5. Masses are determined approximately using the YSO parameters from SED fits and pre-main sequence tracks of Bernasconi & Maeder (1996). The black curve represents the Kroupa IMF (Kroupa 2001) fitted to the peak of the distribution. We estimate the mass of the underlying YSO population from the area under the curve from 0.01 M$_{\odot}$ to 120 M$_{\odot}$.</figcaption></figure>

where A, B, and C are scaling factors. We follow the method described in Immer et al. (2012) and fit our mass distribution histogram with a curve of the form as in Eq. (2) by non-linear least squares fitting routine. The fitting results in a value of A = 7339, which we used to obtain $\zeta(M)$ at M = 0.5 M$_{\odot}$ assuming a continuous IMF and thus estimate B = 14677 from Eq. (3). We carried out the same exercise to estimate C = 183 464. Figure 10 shows the mass distribution histogram and the fit we performed on the distribution. Finally, we estimate the total mass of YSOs to be ~35 000 M$_{\odot}$ in the CMZ using

$$M_{\text{tot}} = \int_{0.01}^{120} M\zeta(M)dM. \qquad (5)$$

Assuming all YSOs that constitute our sample have an average age of 0.75 ± 0.25 Myr, we estimated the average SFR to be ~0.046 M$_{\odot}$ yr$^{-1}$. If we assume a different IMF (e.g. Salpeter) and change the integration limits in the mass range in addition to including the mass uncertainties from the individual SED fitting (see Table 3) and the uncertainty in the assumed age, our estimated error in the derived SFR is of the order of ±0.026 M$_{\odot}$ yr$^{-1}$. We also changed our colour criterion by reducing the H-K$_{S}$ cut to 1.0 instead of 1.5 to account for the variability in extinction across CMZ and the estimated SFR is still within the uncertainty limit of +0.026 M$_{\odot}$ yr$^{-1}$.

Our SFR estimate is lower than values from previous studies of YHA09, An11, and Immer et al. (2012). YHA09 and Immer et al. (2012) applied the YSO counting method of photometrically identified YSOs to calculate SFR of ~0.14 M$_{\odot}$ yr$^{-1}$ (YSO lifetime ~ 0.1 Myr) and ~0.08 M$_{\odot}$ yr$^{-1}$ (YSO lifetime ~ 1 Myr), respectively. An11 carried out a spectroscopic identification of YSOs among sources in common with YHA09 and derived a value of 0.07 M$_{\odot}$ yr$^{-1}$ based on the 50% contamination they found. Based on the re-examination of YHA09 sample using radiative transfer models and realistic synthetic observations, Koepferl et al. (2015) estimated the SFR to be lower by a factor of three or more.

In addition to the YSO counting method, there have been studies that have employed the infrared luminosity-SFR relation,

free-free emission from the ionised gas (i.e. bremsstrahlung radiation) at cm wavelengths to estimate the mass of the underlying YSO population, column density threshold, and volumetric star forming relations to estimate and predict, in turn, the SFR in the CMZ (Longmore et al. 2013a; Barnes et al. 2017). Longmore et al. (2013a) estimated the SFR in the |l| < 1°0, |b| < 0°5 to be ~0.015 M$_{\odot}$ yr$^{-1}$ based on the free-free emission contribution to the 33 GHz flux using Wilkinson Microwave Anisotropy Probe (WMAP) data. But the predictions from the column density threshold and volumetric star formation relations exceed the observed SFR with estimates of 0.78 M$_{\odot}$ yr$^{-1}$ and 0.41 M$_{\odot}$ yr$^{-1}$, respectively. Given that predictions from these star formation relations and models are largely dependent on the mass of dense gas, it is important to make sure that different tracers of the dense gas are reliable probes. For example, Mills & Battersby (2017) have shown that HNCO might be a better cloud mass probe than HCN 1–0 in the Galactic centre environment.

Barnes et al. (2017) found an average global SFR of ~0.09 ± 0.02 M$_{\odot}$ yr$^{-1}$ in the same l, b range from the luminosity-SFR relations using 24 µm, 70 µm, and total infrared bolometric luminosity. Based on the observational evidence that the individual clouds and clusters are connected along a coherent velocity structure in position-position-velocity (PPV) space (Henshaw et al. 2016), Barnes et al. (2017) determined the SFR of individual clouds in the CMZ using the dynamical orbit model of Kruijssen et al. (2015) assuming that star formation within these clouds is tidally triggered at the pericentre of the orbit (Longmore et al. 2013b). These authors find the total SFR within these clouds to be 0.03 to 0.071 M$_{\odot}$ yr$^{-1}$. Table 4 lists the details of the SFR estimated in the CMZ using various methods based on past studies.

Thus, the SFR estimate in the CMZ from different methods (including our estimate) all point to a lower value than expected given the large reservoir of dense gas available. There are several physical explanations attributed to this dearth of star formation in the CMZ. Longmore et al. (2013a) suggested that the additional turbulent energy in the gas, as indicated by the larger internal cloud velocity dispersion, could be providing support against gravitational collapse. The other explanations include episodic star formation in the CMZ due to spiral instabilities, high turbulent pressure in the CMZ, and that the gas is not self-gravitating as discussed in detail in Kruijssen et al. (2014).

## 6. Summary

With the aim of estimating the SFR in the CMZ using spectroscopic identification of YSOs, we prepared a detailed observation of 22 fields using KMOS. From the 8 fields we observed, we extracted clean spectra for 91 sources. Based on the CO absorption found in cool, late-type stars and Bry emission seen in YSOs, we were able to clearly separate YSOs from cool, late-type stars in the EW(CO) vs. EW(Bry) diagram. We plotted our spectroscopically classified YSOs and late-type stars in the colour-colour and colour-magnitude diagrams used in the literature to classify YSOs. We found that various criteria used to classify YSOs in such diagrams were not able to remove contaminants. We suggest a new criterion in the H-K$_{S}$ vs. H-[8.0] colour-colour diagram wherein we see a clear separation of YSOs and late-type stars.

We used the new and improved version of SED models for YSOs in R17 to fit the observed photometry in the wavelength range of 1.25–24 µm for 8 YSOs. From the radii and temperatures we obtained from the SED fit, we estimated their masses to be greater than 8 M$_{\odot}$. Since we needed a bigger sample
<table><caption>**Table 4.** Details of the SFR estimated using different methods in the literature.</caption><thead><tr><th>Method</th><th>Region covered</th><th>SFR (M<sub>☉</sub> yr<sup>-1</sup>)</th><th>Reference</th></tr></thead><tbody><tr><td>YSO counting (photometric criterion)</td><td>|l| &lt; 1°3, |b| &lt; 0°17</td><td>0.14<sup>a</sup></td><td>YHA09</td></tr><tr><td>YSO counting (spectroscopic criterion)</td><td>...</td><td>0.07<sup>a</sup></td><td>An11</td></tr><tr><td>YSO counting (photometric criterion)</td><td>|l| &lt; 1°5, |b| &lt; 0°5</td><td>0.08<sup>b</sup></td><td>Immer et al. (2012)</td></tr><tr><td rowspan="2">Free-free emission – SFR</td><td>|l| &lt; 1°0, |b| &lt; 0°5</td><td>0.015</td><td>Longmore et al. (2013a)</td></tr><tr><td>|l| &lt; 1°0, |b| &lt; 1°0</td><td>0.06</td><td>Longmore et al. (2013a)</td></tr><tr><td>Column density threshold</td><td>|l| &lt; 1°0, |b| &lt; 0°5</td><td>0.78</td><td>Longmore et al. (2013a)</td></tr><tr><td>Volumetric SF relations</td><td>|l| &lt; 1°0, |b| &lt; 0°5</td><td>0.41</td><td>Longmore et al. (2013a)</td></tr><tr><td>Infrared luminosity–SFR</td><td>|l| &lt; 1°0, |b| &lt; 0°5</td><td>0.09±0.02</td><td>Barnes et al. (2017)</td></tr><tr><td>YSO counting</td><td>|l| &lt; 1°5, |b| &lt; 0°5</td><td>0.046 ± 0.026</td><td>This work</td></tr></tbody></table>

**Notes.** The method used to estimate the SFR, the region of the CMZ covered, estimated SFR, and corresponding references are listed.ª Assumed age of YSOs ~ 0.1 Myr.ᵇ Assumed age of YSOs ~ 1 Myr.

to estimate the SFR in the CMZ, we searched for sources within $|l| < 1°5$ and $|b| < 0°5$ with valid photometry in *H*, $K_S$ (IRSF catalogue), 8.0 µm (Ramírez et al. 2008) and 15 µm (ISOGAL PSC), or 24 µm (Gutermuth & Heyer 2015) bands. We identified 334 YSOs based on the criterion we defined in the *H* − $K_S$ vs. *H* − [8.0] diagram. We performed SED fits for these sources using R17 models resulting in 190 sources with a good fit, and their estimated masses range from 2.7 to 35 $M_☉$, peaking at ~8 $M_☉$. The total mass of YSOs in the CMZ was then estimated to be ~35 000 $M_☉$ by extrapolating to lower masses using a Kroupa IMF between 0.01 and 120 $M_☉$. Assuming an average age of 0.75 ± 0.25 Myr for YSOs, we estimate the SFR to be ~0.046 ± 0.026 $M_☉$ yr$^{-1}$, that is slightly lower than found in previous studies.

It is necessary to carry out follow-up spectroscopic infrared observations to obtain a statistically significant YSO sample to further constrain our colour–colour criterion to identify YSOs. This will help us accurately determine the SFR, which is an important ingredient in the chemical evolution models of the Galaxy, as well as to understand star formation of the Galactic centre and to serve as a template for circumnuclear star formation in the other galactic nuclei.

**Acknowledgements.** We wish to thank the anonymous referee for the extremely useful comments that greatly improved the quality of this paper. G.N. and M.S. acknowledges the Programme National de Cosmologie et Galaxies (PNCG) of CNRS/INSU, France, for financial support. The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement No. [614922]. This research made use of the SIMBAD database (operated at CDS, Strasbourg, France).

## References

Adams, F. C., Lada, C. J., & Shu, F. H. 1987, ApJ, 312, 788

Allen, L. E., Calvet, N., D’Alessio, P., et al. 2004, ApJS, 154, 363

An, D., Ramírez, S. V., Sellgren, K., et al. 2009, ApJ, 702, L128

An, D., Ramírez, S. V., Sellgren, K., et al. 2011, ApJ, 736, 133

Bally, J., Stark, A. A., Wilson, R. W., & Henkel, C. 1988, ApJ, 324, 223

Bania, T. M. 1977, ApJ, 216, 381

Barnes, A. T., Longmore, S. N., Battersby, C., et al. 2017, MNRAS, 469, 2263

Bernasconi, P. A., & Maeder, A. 1996, A&A, 307, 829

Bik, A., Kaper, L., & Waters, L. B. F. M. 2006, A&A, 455, 561

Boehle, A., Ghez, A. M., Schödel, R., et al. 2016, ApJ, 830, 17

Cooper, H. D. B., Lumsden, S. L., Oudmaijer, R. D., et al. 2013, MNRAS, 430, 1125

Dong, H., Mauerhan, J., Morris, M. R., Wang, Q. D., & Cotera, A. 2015, MNRAS, 446, 842

Dutra, C. M., & Bica, E. 2001, A&A, 376, 434

Fatuzzo, M., & Melia, F. 2009, PASP, 121, 585

Feldmeier-Krause, A., Neumayer, N., Schödel, R., et al. 2015, A&A, 584, A2

Felli, M., Testi, L., Schuller, F., & Omont, A. 2002, A&A, 392, 971

Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96

Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Rev. Mod. Phys., 82, 3121

Gillessen, S., Plewa, P. M., Eisenhauer, F., et al. 2017, ApJ, 837, 30

Glass, I. S., Matsumoto, S., Carter, B. S., & Sekiguchi, K. 2001, MNRAS, 321, 77

Greene, T. P., & Lada, C. J. 1996, AJ, 112, 2184

Guttermuth, R. A., & Heyer, M. 2015, AJ, 149, 64

Henshaw, J. D., Longmore, S. N., Kruijssen, J. M. D., et al. 2016, MNRAS, 457, 2675

Immer, K., Schuller, F., Omont, A., & Menten, K. M. 2012, A&A, 537, A121

Koeperfer, C. M., & Robitaille, T. P. 2017, ApJ, 849, 3

Koeperfer, C. M., Robitaille, T. P., Morales, E. F. E., & Johnston, K. G. 2015, ApJ, 799, 53

Kroupa, P. 2001, MNRAS, 322, 231

Kruijssen, J. M. D., Longmore, S. N., Elmegreen, B. G., et al. 2014, MNRAS, 440, 3370

Kruijssen, J. M. D., Dale, J. E., & Longmore, S. N. 2015, MNRAS, 447, 1059

Lada, C. J. 1987, in Star Forming Regions, eds. M. Peimbert, & J. Jugaku, IAU Symp., 115, 1

Longmore, S. N., Bally, J., Testi, L., et al. 2013a, MNRAS, 429, 987

Longmore, S. N., Kruijssen, J. M. D., Bally, J., et al. 2013b, MNRAS, 433, L15

Martin, C. L., Walsh, W. M., Xiao, K., et al. 2004, ApJS, 150, 239

Matsunaga, N., Kawadu, T., Nishiyama, S., et al. 2009, MNRAS, 399, 1709

Mauerhan, J. C., Cotera, A., Dong, H., et al. 2010a, ApJ, 725, 188

Mauerhan, J. C., Muno, M. P., Morris, M. R., Stolovy, S. R., & Cotera, A. 2010b, ApJ, 710, 706

Megeath, S. T., Allen, L. E., Guttermuth, R. A., et al. 2004, ApJS, 154, 367

Mills, E. A. C., & Battersby, C. 2017, ApJ, 835, 76

Morris, M., & Serabyn, E. 1996, ARA&A, 34, 645

Muno, M. P., Baganoff, F. K., Bautz, M. W., et al. 2003, ApJ, 599, 465

Nagashima, C., Nagayama, T., Nakajima, Y., et al. 1999, in Star Formation 1999, ed. T. Nakamoto, 397

Nagayama, T., Nagashima, C., Nakajima, Y., et al. 2003, in Instrument Design and Performance for Optical/Infrared Ground-based Telescopes, eds. M. Iye, & A. F. M. Moorwood, Proc. SPIE, 4841, 459

Nishiyama, S., Nagata, T., Kusakabe, N., et al. 2006, ApJ, 638, 839

Nishiyama, S., Tamura, M., Hatano, H., et al. 2009, ApJ, 696, 1407

Oka, T., Geballe, T. R., Goto, M., Usuda, T., & McCall, B. J. 2005, ApJ, 632, 882
Omont, A., Gilmore, G. F., Alard, C., et al. 2003a, A&A, 403, 975

Omont, A., Gilmore, G. F., Alard, C., et al. 2003b, VizieR Online Data Catalog: II/243

Ramírez, S. V., Sellgren, K., Carr, J. S., et al. 2000, ApJ, 537, 205

Ramírez, S. V., Arendt, R. G., Sellgren, K., et al. 2008, ApJS, 175, 147

Robitaille, T. P. 2011, A&A, 536, A79

Robitaille, T. P. 2017, A&A, 600, A11

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256

Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328

Schödel, R., Najarro, F., Muzic, K., & Eckart, A. 2010, A&A, 511, A18

Schuller, F., Omont, A., Glass, I. S., et al. 2006, A&A, 453, 535

Schultheis, M., Lançon, A., Omont, A., Schuller, F., & Ojha, D. K. 2003, A&A, 405, 531

Schultheis, M., Sellgren, K., Ramírez, S., et al. 2009, A&A, 495, 157

Serabyn, E., & Morris, M. 1996, Nature, 382, 602

Sharples, R., Bender, R., Agudo Berbel, A., et al. 2013, The Messenger, 151, 21

Stolovy, S., Ramirez, S., Arendt, R. G., et al. 2006, in J. Phys. Conf. Ser., 54, eds. R. Schödel, G. C. Bower, M. P. Muno, S. Nayakshin, & T. Ott, 176

Støstad, M., Do, T., Murray, N., et al. 2015, ApJ, 808, 106

van Dokkum, P. G. 2001, PASP, 113, 1420

Yusef-Zadeh, F., Hewitt, J. W., Arendt, R. G., et al. 2009, ApJ, 702, 178

Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481