Highlight
Free Access
Issue
A&A
Volume 587, March 2016
Article Number L6
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201527500
Published online 01 March 2016

© ESO, 2016

1. Introduction

Core helium-burning red clump (RC) stars are useful distance indicators, because their magnitude changes slowly and smoothly with age and metallicity, and is predicted very accurately by stellar evolution models (e.g., Salaris & Girardi 2002). The RC is also a bright and distinct feature in the colour-magnitude diagram (CMD) of intermediate-to-old stellar populations, hence easy to identify. This makes it a powerful tracer of the Milky Way (MW) bulge structure. In particular, RC stars have been extensively used to trace the bar position, orientation and scale parameters (see the review of Gonzalez & Gadotti 2015, for a complete list of references), as well as to uncover the presence of a peanut (X-shape) in the outer bulge (McWilliam & Zoccali 2010; Nataf et al. 2010; Saito et al. 2011).

The infrared photometry from the VISTA Variables in the Vía Láctea (VVV) ESO public survey (Minniti et al. 2010) has been used by Wegg & Gerhard (2013) to constrain the first 3D density model of the bulge. Their model was rather poorly constrained close to the plane because the DR1 catalogues (Saito et al. 2012) used in that study are based on aperture photometry, which suffers strongly from incompleteness because of the high level of crowding in the region | b | ≤ 0.5°.

Here we present the first stellar density map that reaches all the way to the centre of the MW bulge, using RC star counts from new VVV catalogs that are based on PSF-fitting photometry. We directly compare the stellar density with kinematics, and derive the first fully empirical estimate of the stellar mass of the bulge.

2. Observation and data reduction

Our data are part of the VVV survey that covers ~180 sq. deg of MW bulge between | l | ≤ 10° and | b | ≤ 4.5° contiguously. A detailed description of the global survey, and of the multiband observations completed within the first two years (used in this study), can be found in Minniti et al. (2010) and Saito et al. (2012), respectively.

Unlike most previous studies that are based on DR11 VVV data, we performed PSF-fitting photometry with the DoPHOT code (Schechter et al. 1993; Alonso-García et al. 2012) on each of the 16 chips in the images, both in J and Ks, and for this study we retrieved the stacked pawprint images from the CASU pipeline2, following Alonso-García et al. (2015).

The photometric catalogues of the individual chips were calibrated and astrometrized to the VISTA system. The six pawprint catalogs that cover a VISTA tile were combined into tile catalogs and the magnitudes of stars, which were detected multiple times, weight-averaged.

The photometric completeness was estimated via artificial star experiments. For each tile, stars were randomly added to chip#12 only, assuming that the stellar density in this chip was representative of the whole tile. A more detailed description of the PSF photometry and completeness experiments for the whole VVV survey area is the subject of a forthcoming paper (Alonso-Garcia et al., in prep.). DoPHOT allowed us to obtain photometry with a typical error < 0.05 down to J ~ 18 and Ks ~ 17 in most of the bulge area. As expected, in the faint (Ks> 17) magnitude regime, the completeness of the catalogs can vary significantly across the bulge area, being less complete for tiles closer to the plane, where both crowding and extinction are more severe. However, in the magnitude range of interest for the present work, i.e., K0 ≤ 14, the derived catalogs are more than 70% complete, with the exception of a few fields (see Fig. 1).

thumbnail Fig. 1

Completeness map of the RC peak magnitude (K0RC\hbox{$K_0^{\rm RC}$}) across the bulge.

3. Analysis and results

We aim to investigate the 3D structure and symmetry of the inner bulge in detail, by tracing the stellar density from RC star counts. To build an RC-based stellar density map, we needed to address two main issues that relate to the observing strategy that was adopted by VVV when surveying the bulge area. The first is that there is a non-negligible overlap (~10%) between adjacent tiles. We therefore cross-correlated adjacent tiles to eliminate double detections. The second issue is that the centers of the tiles are not symmetric with respect to the plane (b = 0°), nor to the bulge projected minor axis (l = 0°). Hence we resampled the bulge area analysed here with a new uniform grid of tiles, hereafter called fields, each with a 1° × 1° area that is centred at integer values of both latitude and longitude between −9.5° ≤ l ≤ 10.5° and | b | ≤ 4.5°. This allows us to directly and consistently compare fields that are located symmetrically around the bulge projected axes without the need for interpolations on stellar density or completeness.

The catalogs were corrected for extinction by applying the reddening map of Gonzalez et al. (2011b) and adopting the Nishiyama et al. (2005) extinction law. From the de-reddened catalog of each field, stars lying on the red giant branch (RGB) were selected using a variable color cut aimed at excluding stars in the bright (K0 ≤ 15) main sequence of the foreground disk, with a typical color (JK)0 ~ 0.3. We then used RGB stars to derive the K0 luminosity function (LF), corrected this for incompleteness, and measured the RC mean magnitude (K0RC\hbox{$K^{\rm RC}_0$}) by following the prescription of Stanek et al. (1997). As such, the RGB LF was fitted with a second-order polynomial, together with a Gaussian to account for the peak of RC stars. In addition, as done by Nataf et al. (2011), Gonzalez et al. (2011b) and Wegg & Gerhard (2013), a second Gaussian was included in the fit to account for the presence of a second peak, which is less populated and fainter than the RC (Fig. 2). The position and relative intensity of the second peak is broadly consistent with the expected position of the RGB bump (Sweigart et al. 1990; Wegg & Gerhard 2013). However, as Gonzalez et al. (2011b), we noticed that the magnitude difference between the bulge RC and this faint peak is not constant close to the Galactic plane, as one would expect if it were a stellar evolutionary feature. The present data allowed us to map both the magnitude and the density of this second peak with respect to the bulge RC, and will be the subject of a forthcoming paper.

thumbnail Fig. 2

K0 LF (black histogram) for representative fields along the Galactic plane, after subtracting the contribution from the RGB stars. The red solid line shows the best fit to each LF, while the dashed lines identify the individual Gaussian fits to the RC (blue) and to the additional fainter peak (gray).

3.1. The RC density map

thumbnail Fig. 3

Density map in the longitude-latitude plane based on RC star counts. Star counts have been normalized to the maximum (Max). Solid contours are isodensity curves, linearly spaced by 0.1 × Max deg-2. The fields observed by GIBS are shown as circles, whose colors follow the measured velocity dispersion (σ in km s-1) in that given field.

By using the total number of RC stars detected in all the fields, normalized to the their maximum value across the whole area, we derive the stellar density map shown in Fig. 3. As expected, the map clearly shows a boxy/peanut bulge, with an increase of star counts towards the center. The coloured circles on top of the grayscale map show the GIBS fields (Zoccali et al. 2014), where radial velocities were measured for 200400 individual RC stars per field. The symbols are colour-coded according to the radial velocity dispersion measured in each of them. It is clear that the velocity dispersion follows the density of RC stars, and hence, the mass density. Therefore, the present map confirms the presence of a high stellar density peak in the central ~150 pc, although the peak is more axisymmetric than the one extrapolated using the GIBS fields (see their Fig. 11).

In general, the vertical extent of the isodensity contours is larger at positive longitude, as expected for a bar that has its close side at positive longitude. This effect was previously noted by Alard (2001) using 2MASS photometry, although limited to | b | ≥ 2°. The high quality of the VVV PSF-fitting photometry allows us to map the stellar density of the innermost bulge regions (| b | ≤ 1.5°) for the first time. The maximum density of stars is found in the region | l | ≤ 1° and | b | ≤ 0.5°, and, as shown by the 0.8 and 0.7 isodensity curves, appears slightly asymmetric with respect to the bulge’s minor axis. The central asymmetry becomes more noticeable as we progressively move outwards. In fact, as shown in Fig. 4, in the region between | l | ≤ 5°, | b | ≤ 2° we find that the counts are systematically higher at b< 0. The additional overdensity along the Galactic plane at l ~ −7° is also present in the GLIMPSE 4.5 μm map (Churchwell et al. 2009). Finally, as discussed before, the absolute density of the central field and the field at (−7°,0°) should be taken as a lower limit because these are the only two fields for which the completeness corrections at K0RC\hbox{$K_0^{\rm RC}$} are larger than 50% (Fig. 1).

thumbnail Fig. 4

Stellar density profile across the longitude, at fixed latitudes. Filled symbols with dotted lines show star counts at negative latitudes, while empty symbols and solid lines refers to positive latitudes.

3.2. A fully empirical estimate of the bulge stellar mass

In what follows we present a fully empirical estimate of the stellar mass of the bulge within the VVV area. We use the empirical initial mass function (IMF) measured by Zoccali et al. (2000) (recently confirmed by Calamida et al. 2015) in a HST/NICMOS field at (0,−6°), integrated to yield the mass contained in that field, and scale it to the whole bulge using the ratio of RC stars. The field of view (FoV) of NICMOS is too small to sample evolved stars, but a complete LF in a larger area (~8′ × 8′; from NTT/SOFI), including the NICMOS field, is presented in Zoccali et al. (2003). The scaling factor between NICMOS and SOFI (×609) is simply the ratio of the FoV areas, given that they sample the same bulge region. From this scaling factor, we can derive the mass contained in the SOFI field. Finally, to obtain the total stellar mass in the VVV area, we measure the scaling factor between the number of RGB+RC stars in the whole VVV area and in the SOFI field.

In practice, the LF for the VVV catalogues was restricted to the RGB at the RC level (K0RC±1.5\hbox{$K_0^{\rm RC}\pm1.5$} mag) where the completeness is high (and thus the correction small) across the whole bulge area. Because the RC magnitude varies across the fields, as a result of the bar orientation, we measured the RC position in each field individually and counted stars within ±1.5 mag of the peak. This procedure enabled us to obtain a spatial mass profile in 1° bins (Fig. 5). In the small region, where the RC is split because of the X-shape of the bulge (| l | < 4°, for b< −5°), we use the mean value of the two RCs. This is accurate enough for our purpose, which is just to find a normalization factor between the RGB in the SOFI and in the whole VVV area. Furthermore, we can assume symmetry around the Galactic plane to extend our mass estimate to the ~squared region, – hereafter VVVΛ| l | < 10° and | b | < 9.5° (~2.8 × 2.7 kpc, assuming a distance to the Galactic center of 8 kpc). In other words, the number of RGB+RC stars for b> 4.5°, outside the VVV area, is assumed to be equal to that at b< −4.5°. Finally, since the LF from the SOFI field was decontaminated of disk stars (amounting to 18% within ±1.5 mag around the RC), the number of RGB+RC stars in VVVΛ was also reduced by this quantity. The ratio between the total number of RGB+RC stars within ±1.5 mag from the RC in VVVΛ vs. SOFI is 58624.

In Zoccali et al. (2000) the IMF was integrated using the observational data for main-sequence stars (from the turnoff down to 0.15 M) and different assumptions for the IMF above 1 M and below 0.15 M. From recent determinations of the IMF of massive stars (see Bastian et al. 2010, for a review), we assume a Salpeter IMF (α = −2.35) above 1 M and a power law with a slope of α = −0.3 in the substellar regime. With the same initial-to-final-mass relation as in Zoccali et al. (2000), we found 570 M in the NICMOS field, hence 570 × 609 × 58 624 = 2.0 × 1010M in VVVΛ. This is the mass in stars and remnants of the Galactic bulge in the region | b | < 9.5°, | l | < 10°. It assumes a uniform disk contamination (at the level of 18% close to the RC) across the whole bulge area. While this assumption is accurate for fields close to (0,−6°), it is likely an overestimate for regions closer to the Galactic plane, and an underestimate for outer regions3. Currently, an accurate estimate of the disk-to-bulge fraction for the inner few kpc of the Galaxy is not yet available. However, by adopting a smoothly varying disk contamination from 20% at b = ± 9° down to 5% at b = ± 2°, and then back to 20% at b = 0° where the thin disk might become significant again (e.g., Dékány et al. 2015), the total mass would only increase by 7%. Other systematic error sources, such as the slope of the IMF above 1 M and below 0.08 M, or the initial-to-final mass, are estimated as follows. A generous ±0.3 dex variation of the assumed IMF slope changes the total mass by less than 12%. We further assume that small gradients on the bulge’s dominant stellar populations do not change the number ratio between RC and main-sequence stars by more than a few percent. In summary, a conservative total uncertainty amounts to 15%, which yields a total bulge mass within VVVΛ of 2.0 ± 0.3 × 1010M, including disk contamination, extrapolation of the IMF, mass of the remnants, and stellar population gradients.

thumbnail Fig. 5

Stellar mass profile across latitude in each field for fixed longitudes. Filled symbols with dotted lines represent b< 0, empty symbols with solid lines to b> 0 fields. The lower part of the plot has an expanded y-axis, for a better display of the high negative latitude fields. The number counts at b = 0° are rather noisy, especially at negative longitudes, owing to the higher extinction, hence higher incompleteness in those fields. The purple line includes a mild smoothing from which the central point was excluded, to avoid flattening the peak.

4. Conclusions

We have used the VVV J and Ks single-epoch data to derive accurate PSF-fitting photometry in the bulge region between −9.5° ≤ l ≤ 10.5° and | b | ≤ 4.5°. The photometric accuracy of the catalogs, which include modest incompleteness corrections (the result of crowding and extinction) allowed us to build the first density profile of the inner bulge towards the Galactic plane by means of RC star counts. We found a good agreement between the measured RC density distribution and the velocity dispersion map provided by the GIBS survey. In particular, the central high σ-peak at | l | ≤ 2° and | b | ≤ 3° in the derived velocity dispersion curve is well matched by a high stellar density peak in the central ~150 pc. This tight link between the stellar density profile and the velocity dispersion map confirms that the flattening of the K0RC\hbox{$K^{\rm RC}_0$} trend with longitude, for | l | ≤ 4°, (Nishiyama et al. 2005; Gonzalez et al. 2011a) is more likely the result of an inner compact and axisymmetric spheroid rather than the presence of a nuclear bar (Gerhard & Martinez-Valpuesta 2012).

We further exploit the depth and high quality of the new VVV photometry that also extends to the innermost bulge

regions, and derive the first purely empirical estimate of the bulge stellar mass by scaling the observed RC stellar density map with the observed bulge LF from Zoccali et al. (2000), Zoccali et al. (2003). We measure the mass in star and stellar remnants of the Galactic bulge as 2.0 ± 0.3 × 1010M. Previous bulge mass estimates, all model-dependent, span a considerable range from 6.1 × 109 (Robin et al. 2012) to 1.6 × 1010M (Nesti & Salucci 2013; Sofue 2013; Portail et al. 2015). A possible reason why our result, although compatible within the quoted errors, is larger than the mass estimated by dynamical models (e.g. Portail et al. 2015) is because our sampled volume is not limited along the line of sight. In addition, the kinematics used so far to constrain the models, do not extend to the innermost bulge regions (| l | ≤ 1° and b = −1°, −2°), where a velocity dispersion peak was found (Zoccali et al. 2014).


1

Catalogs obtained from aperture photometry on individual tiles.

3

The disk contamination in three fields at different latitudes was estimated in Zoccali et al. (2008) using the Besançon model. At (0, −6°) it is very close to 20%, in Baade’s Window it is significantly smaller, while at (0, −12°) is larger (cf., their Table 4).

Acknowledgments

The authors gratefully acknowledge the use of data from the ESO Public Survey program ID 179.B-2002, taken with the VISTA telescope, and data products from the Cambridge Astronomical Survey Unit. Support for M.Z. and D.M. is provided by the BASAL CATA Center for Astrophysics and Associated Technologies through grant PFB-06, and the Ministry for the Economy, Development, and Tourism’s Programa Iniciativa Cientifica Milenio through grant IC120009, awarded to Millenium Institute of Astrophysics (MAS). D.M. and M.Z. acknowledge support from FONDECYT regular grant No. 1130196 and 1150345, respectively. J.A.G. acknowledges support by the FIC-R Fund, allocated to Project 30321072. M.H. acknowledges financial support from the Chilean BASAL Centro de Excelencia en Astrofisica y Tecnologias Afines (CATA), grant PFB-06/2007.

References

  1. Alard, C. 2001, A&A, 379, L44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alonso-García, J., Mateo, M., Sen, B., et al. 2012, AJ, 143, 70 [NASA ADS] [CrossRef] [Google Scholar]
  3. Alonso-García, J., Dékány, I., Catelan, M., et al. 2015, AJ, 149, 99 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bastian, N., Covey, K. R., & Meyer, M. R. 2010, ARA&A, 48, 339 [NASA ADS] [CrossRef] [Google Scholar]
  5. Calamida, A., Sahu, K. C., Casertano, S., et al. 2015, ApJ, 810, 8 [NASA ADS] [CrossRef] [Google Scholar]
  6. Churchwell, E., Babler, B. L., Meade, M. R., et al. 2009, PASP, 121, 213 [NASA ADS] [CrossRef] [Google Scholar]
  7. Dékány, I., Minniti, D., Majaess, D., et al. 2015, ApJ, 812, L29 [NASA ADS] [CrossRef] [Google Scholar]
  8. Gerhard, O., & Martinez-Valpuesta, I. 2012, ApJ, 744, L8 [NASA ADS] [CrossRef] [Google Scholar]
  9. Gonzalez, O. A., & Gadotti, D. A. 2015, in Galactic Bulges, eds. E. Laurikainen, R. Peletier, D. Gadotti (Springer Publishing), in press [arXiv:1503.07252] [Google Scholar]
  10. Gonzalez, O. A., Rejkuba, M., Minniti, D., et al. 2011a, A&A, 534, L14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Gonzalez, O. A., Rejkuba, M., Zoccali, M., Valenti, E., & Minniti, D. 2011b, A&A, 534, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. McWilliam, A., & Zoccali, M. 2010, ApJ, 724, 1491 [NASA ADS] [CrossRef] [Google Scholar]
  13. Minniti, D., Lucas, P. W., Emerson, J. P., et al. 2010, New A, 15, 433 [Google Scholar]
  14. Nataf, D. M., Udalski, A., Gould, A., Fouqué, P., & Stanek, K. Z. 2010, ApJ, 721, L28 [NASA ADS] [CrossRef] [Google Scholar]
  15. Nataf, D. M., Udalski, A., Gould, A., & Pinsonneault, M. H. 2011, ApJ, 730, 118 [NASA ADS] [CrossRef] [Google Scholar]
  16. Nesti, F., & Salucci, P. 2013, J. Cosmol. Astropart. Phys., 7, 16 [Google Scholar]
  17. Nishiyama, S., Nagata, T., Baba, D., et al. 2005, ApJ, 621, L105 [NASA ADS] [CrossRef] [Google Scholar]
  18. Portail, M., Wegg, C., Gerhard, O., & Martinez-Valpuesta, I. 2015, MNRAS, 448, 713 [NASA ADS] [CrossRef] [Google Scholar]
  19. Robin, A. C., Marshall, D. J., Schultheis, M., & Reylé, C. 2012, A&A, 538, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Saito, R. K., Zoccali, M., McWilliam, A., et al. 2011, AJ, 142, 76 [NASA ADS] [CrossRef] [Google Scholar]
  21. Saito, R. K., Hempel, M., Minniti, D., et al. 2012, A&A, 537, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Salaris, M., & Girardi, L. 2002, MNRAS, 337, 332 [NASA ADS] [CrossRef] [Google Scholar]
  23. Schechter, P. L., Mateo, M., & Saha, A. 1993, PASP, 105, 1342 [NASA ADS] [CrossRef] [Google Scholar]
  24. Sofue, Y. 2013, PASJ, 65, 118 [NASA ADS] [Google Scholar]
  25. Stanek, K. Z., Udalski, A., Szymański, M., et al. 1997, ApJ, 477, 163 [NASA ADS] [CrossRef] [Google Scholar]
  26. Sweigart, A. V., Greggio, L., & Renzini, A. 1990, ApJ, 364, 527 [NASA ADS] [CrossRef] [Google Scholar]
  27. Wegg, C., & Gerhard, O. 2013, MNRAS, 435, 1874 [Google Scholar]
  28. Zoccali, M., Cassisi, S., Frogel, J. A., et al. 2000, ApJ, 530, 418 [NASA ADS] [CrossRef] [Google Scholar]
  29. Zoccali, M., Renzini, A., Ortolani, S., et al. 2003, A&A, 399, 931 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Zoccali, M., Hill, V., Lecureur, A., et al. 2008, A&A, 486, 177 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Zoccali, M., Gonzalez, O. A., Vasquez, S., et al. 2014, A&A, 562, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Figures

thumbnail Fig. 1

Completeness map of the RC peak magnitude (K0RC\hbox{$K_0^{\rm RC}$}) across the bulge.

In the text
thumbnail Fig. 2

K0 LF (black histogram) for representative fields along the Galactic plane, after subtracting the contribution from the RGB stars. The red solid line shows the best fit to each LF, while the dashed lines identify the individual Gaussian fits to the RC (blue) and to the additional fainter peak (gray).

In the text
thumbnail Fig. 3

Density map in the longitude-latitude plane based on RC star counts. Star counts have been normalized to the maximum (Max). Solid contours are isodensity curves, linearly spaced by 0.1 × Max deg-2. The fields observed by GIBS are shown as circles, whose colors follow the measured velocity dispersion (σ in km s-1) in that given field.

In the text
thumbnail Fig. 4

Stellar density profile across the longitude, at fixed latitudes. Filled symbols with dotted lines show star counts at negative latitudes, while empty symbols and solid lines refers to positive latitudes.

In the text
thumbnail Fig. 5

Stellar mass profile across latitude in each field for fixed longitudes. Filled symbols with dotted lines represent b< 0, empty symbols with solid lines to b> 0 fields. The lower part of the plot has an expanded y-axis, for a better display of the high negative latitude fields. The number counts at b = 0° are rather noisy, especially at negative longitudes, owing to the higher extinction, hence higher incompleteness in those fields. The purple line includes a mild smoothing from which the central point was excluded, to avoid flattening the peak.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.