Issue |
A&A
Volume 660, April 2022
|
|
---|---|---|
Article Number | L3 | |
Number of page(s) | 5 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202243263 | |
Published online | 01 April 2022 |
Letter to the Editor
Distance to the Brick cloud using stellar kinematics
1
Instituto de Astrofísica de Andalucía (CSIC), University of Granada, Glorieta de la astronomía s/n, 18008 Granada, Spain
e-mail: amartinez@iaa.es
2
Max-Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
Received:
3
February
2022
Accepted:
21
March
2022
Context. The central molecular zone at the Galactic center is currently being studied intensively to understand how star formation proceeds under the extreme conditions of a galactic nucleus. Knowing the position of molecular clouds along the line of sight toward the Galactic center has had important implications in our understanding of the physics of the gas and star formation in the central molecular zone. It was recently claimed that the dense molecular cloud G0.253 + 0.016 (the Brick) has a distance of ∼7.20 kpc from the Sun. That would place it outside of the central molecular zone, and therefore of the nuclear stellar disk, but still inside the Bulge.
Aims. Theoretical considerations as well as observational studies show that stars that belong to the nuclear stellar disk have different kinematics from those that belong to the inner Bulge. Therefore, we aim to constrain the distance to the Brick by studying the proper motions of the stars in the area.
Methods. We used ESO HAWK-I/VLT imaging data from epochs 2015 and 2019 to compute proper motions on the Brick and in a nearby comparison field free of dark clouds.
Results. The stellar population seen toward the nuclear stellar disk shows the following three kinematic components: (1) Bulge stars with an isotropic velocity dispersion of ∼3.5 micro-arc second per year; (2) eastward moving stars on the near side of the nuclear stellar disk; and (3) westward moving stars on the far side of the nuclear stellar disk. We clearly see all three components toward the comparison field. However, toward the Brick, which blocks the light from stars behind it, we can only see kinematic components (1) and (2).
Conclusions. While the Brick blocks the light from stars on the far side of the nuclear stellar disk, the detection of a significant component of eastward streaming stars implies that the Brick must be located inside the nuclear stellar disk and, therefore, that it forms part of the central molecular zone.
Key words: proper motions / infrared: general / stars: kinematics and dynamics / Galaxy: center / techniques: high angular resolution / Galaxy: structure
© ESO 2022
1. Introduction
In the Galactic center, at a distance of ∼8 kpc (McNamara et al. 2000; Ghez et al. 2008; Genzel et al. 2010), we can find the central molecular zone (CMZ), a flattened structure of gas and dust within a radius of ∼200 pc (Kruijssen et al. 2014; Tress et al. 2020). Occupying a similar space as the CMZ, there is the so-called nuclear stellar disk (NSD), a rotating structure (Schönrich et al. 2015; Shahzamanian et al. 2022), with a radius of ∼150 pc and a scale height of ∼45 pc (Launhardt et al. 2002; Gallego-Cano et al. 2020; Sormani et al. 2020).
Although most stellar mass in the NSD probably formed ≳8 Gyr ago, it has been one of the most active star forming regions of the Milky Way in the past 100 Myr (Matsunaga et al. 2011; Nogueras-Lara et al. 2020). Nevertheless, the current star formation rate appears to be about a factor of ten lower than expected, considering the high gas density and dense molecular clouds present in the area (Longmore et al. 2013).
In order to fully understand the star formation history in the Galactic center and the evolution of the Galactic Bulge, it is crucial to better understand the structures of the NSD and CMZ (Sormani et al. 2022). One important question to address that would help us move toward a better understanding of this region is where the dense molecular clouds are located along the line of sight (Kruijssen et al. 2015).
In a recent publication by Zoccali et al. (2021), the molecular cloud G0.53+0.016 (the Brick) was suggested to be at a distance of 7.2 ± 0.20 kpc from the Sun, that is to say outside the NSD, but still inside the Bulge. This result is based on the photometric analysis of the Brick stars, which they found brighter than they should be if they were placed at the Galactic center.
This claim was recently challenged by Nogueras-Lara et al. (2021a), who found that the Brick presents a stellar population compatible with that of the NSD (Nogueras-Lara et al. 2020). They concluded that it is located at a distance of 8.4 ± 0.5 kpc, which would place it inside of the CMZ and NSD.
In this Letter we cross-check the work of Nogueras-Lara et al. (2021a) with kinematic data, considering one of these two possibilities based on the assumption that the Brick is mostly opaque to near-infrared light (Longmore et al. 2012) and, therefore, we can only detect stars in front of the cloud. The first option being that the Brick is outside of the NSD and inside the Bulge, in which case we would expect to find only stars with typical Bulge proper motions, in other words zero mean proper motion and isotropic velocity dispersion (Kunder et al. 2012; Shahzamanian et al. 2022). Alternatively, the second option being that the Brick is inside the NSD, in which case we would expect to find two different kinematic groups: one belonging to the Bulge population and the other to the NSD itself, that is to say with stars on the near side of the disk moving eastward and those on the far side moving westward (Schönrich et al. 2015; Shahzamanian et al. 2022).
2. Data and methods
We combined imaging data from the wide-field near-infrared camera HAWK-I/VLT from two epochs, 2015 and 20191. The 2015 data, henceforth D15, correspond to the GALACTICNUCLEUS survey (GNS, Nogueras-Lara et al. 2018, 2019), acquired with fast photometry mode and reduced with the speckle holography algorithm (Schödel et al. 2013) to provide a 0.2″ homogeneous angular resolution (Nogueras-Lara et al. 2019). The 2019 data, henceforth D19, were acquired with the new ground layer adaptive optics assisted by laser (GRAAL), with an average angular resolution of 0.4″ (Zoccali et al. 2021) for H and Ks bands.
The images from both data sets overlap on a small region on the Brick, which is outlined in Fig. 1, along with the comparison field. Both fields have the same dimensions of 0.027° by 0.034° (3.7 pc by 4.7 pc), centered on the coordinates shown in Table 1, at a projected distance of ∼0.22° / 30 pc and ∼0.28° / 39 pc, respectively, to the Galactic east of Sgr A*.
![]() |
Fig. 1. Pointings from GNS (outlined in black) and D19 (in white) overlaid on a 4.5 μm Spitzer/IRAC image (Schödel et al. 2014). The GNS pointings cover smaller areas because of detector windowing. The overlapping areas used for proper motion calculations are shown in blue (on the Brick) and green (comparison field). |
Comparison and Brick field coordinates.
2.1. D15
The acquisition and reduction of D15 is described in Nogueras-Lara et al. (2018, 2019). Due to saturation issues with the Ks band for stars brighter than ∼11.5 mag, we used only H band data for the calculation of proper motions.
2.2. D19
We used two pointings of D19 (outlined in white in Fig. 1). These data consist of 15 jittered exposures (DIT = 10 s, NDIT = 1) and 32 jittered exposures (DIT = 10 s, NDIT = 3), respectively. We applied standard data reduction (sky subtraction, flat fielding, and bad pixel correction). The sky was produced from the jittered images by stacking the different frames and selecting the pixels with the smallest values. This procedure was necessary because the extreme source density toward the target does not result in a homogeneous sky when applying a median. For the photometry and astrometry, we used StarFinder (Diolaiti et al. 2000). In order to compute the mean position and uncertainties of the detected stars, we aligned all jittered frames employing the first one as a reference, using a linear polynomial and considered stars in different frames to be the same ones if they were less than 1 pixel apart. Due to the long exposure time of these images, nonlinearity and saturation effects become increasingly important for stars brighter than H = 14, resulting in increased astrometric uncertainties (Fig. 2). Therefore, for the fine alignment of the frames, we considered only stars with magnitudes 14 ≤ H ≤ 16, which have the smallest position uncertainties according to Fig. 2. We only accepted stars detected in all frames. We used the mean and the uncertainty of the mean of the positions in all frames. We only used areas covered by at least four pointings. We calibrated D19 photometrically via common stars with the SIRIUS/IRSF catalog of Nishiyama et al. (2006), as was done for the GNS data (Nogueras-Lara et al. 2019).
![]() |
Fig. 2. Relative astrometric uncertainty as a function of H magnitude for all stars in the Brick and comparison fields that are common to D15 (red dots) and D19 (black dots). |
2.3. Proper motions
We aligned D19 with D15 using common stars. We applied a color cut H − Ks > 1.3 to remove foreground stars (see Nogueras-Lara et al. 2020). For an estimation of the displacement and rotation angle between both data sets, we made a first alignment by finding similar three-point asterisms with the astroalign package in python (Beroiz et al. 2020) and a second one with a polynomial fit (IDL polywarp), considering two stars to be the same if they were less than 1 pixel apart. This procedure was iterated until the number of common stars remained stable. We computed the alignment uncertainties using a Jackknife resampling approach; we found that a polynomial of degree 2 resulted in the most accurate solution.
After aligning the lists, we computed the velocities for each star by subtracting the positions of common stars and divided them by the time baseline (∼4 year). We computed the uncertainties quadratically, propagating the errors in position for each star and the alignment uncertainty for the D19 stars. From that point on, we considered only stars with proper motion uncertainties below 2.0 mas yr−1 and with an absolute magnitude difference with GNS smaller than 0.5 magnitudes (Fig. 3).
![]() |
Fig. 3. Uncertainties for vx and vy as a function of H magnitude for stars with colors H − Ks > 1.3 on the Brick field. We considered only stars with velocity uncertainties smaller than 2 mas yr−1 (red dashed line). We did not consider stars with an absolute difference in magnitude with GNS bigger than 0.5 (red dots). We applied the same criteria in both fields. |
3. Results
We studied stellar kinematics in two fields: the Brick field and the Comparison field to Galactic west of the Brick on the NSD (Fig. 1).
3.1. Comparison field
Figure 4 (top left) shows the proper motion distribution perpendicular to the Galactic plane. We computed the binning of the data with the python function numpy.histogram_bin_edges (Harris et al. 2020), and considered the uncertainty for the height of each bin as the square root of the number of the stars in that bin. We fit it with different Gaussian models, using the package dynesty (Speagle 2020) in python, which estimates Bayesian posterior probabilities and evidences. We tried doing this with one, two, and three-Gaussians models. To estimate the goodness of each fit, we used the logarithm of posterior probability, log Z (see Shahzamanian et al. 2022). We found that the best fit is achieved with two Gaussians. We show the best fit Gaussian model along with the proper motion distribution in Fig. 4 (top left) and their parameters in Table 2; uncertainties correspond to a 1σ spread of each parameter’s posterior distribution. We interpret the red Gaussian in Fig. 4 (top left), with the broader velocity dispersion, to represent the stars from the Bulge, in agreement with Kunder et al. (2012) and Shahzamanian et al. (2022), and the narrow ones to represent the stars from NSD.
![]() |
Fig. 4. Proper motion distribution for the comparison field (top row) and the Brick field (bottom row). Left-hand plots correspond to the perpendicular component, where the red Gaussian represents the Bulge stars and the black one corresponds to the NSD stars. Right-hand plots correspond to the parallel component, where the red Gaussian represents the Bulge stars and the yellow and black ones represent the NSD stars moving eastward and westward, respectively. |
The proper motions parallel to the Galactic plane show a broader distribution than the ones perpendicular to the Galactic plane (Fig. 4, top right), which is in agreement with Shahzamanian et al. (2022). We followed the procedure explained in the previous paragraph with perpendicular components. We found the best fit corresponding to a three-Gaussians model. The mean Gaussians and the proper motion distribution are shown in Fig. 4 (top right) and their parameters are found in Table 2. We interpret the three-Gaussians model in Fig. 4 (top right) as follows: bulge stars, with about zero proper motion and high velocity dispersion (Soto et al. 2012; Shahzamanian et al. 2022); and eastward and westward moving stars on the near and far side of the NSD (Schönrich et al. 2015; Shahzamanian et al. 2022).
The lower amplitude of the Gaussian, representing the westward moving stars, is probably a consequence of these stars experiencing higher extinction than the eastward moving stars since the former probably belong to the far side of the NSD. To check this, we de-reddened the stars in the comparison field using a star-by-star approach assuming that the intrinsic color H − Ks is similar for all the stars analyzed (see Nogueras-Lara et al. 2021b). In Fig. 5 we show mean Ks extinction as a function of μl. On average, the extinction of eastward moving stars is smaller than for the westward moving stars, which justifies why we detected fewer westward moving stars (see also Fig. 10 of Shahzamanian et al. 2022).
![]() |
Fig. 5. Extinction as a function of μl for the comparison field. Squares and bars represent the mean extinctions and their statistical uncertainties for stars in bins of 3 mas yr−1 width. Statistical uncertainties for μl are smaller than the width of the squares. |
To compute the mean proper motion values of the NSD and Bulges stars, since we do not consider an absolute frame of reference, our basic assumption is that the mean proper motion of all stars is zero. This assumption does not appear to be valid in the comparison field since we detected more stars from the near side of the NSD than from the far side, that is, there is a net eastward proper motion present in the field. This makes our reference frame drift eastward and this results in an underestimation of eastward motion and an overestimation of westward motions. If we correct this by assuming the same mean eastward and westward velocities, then we see that the near and far sides of the NSD rotate with approximately 2.27 mas yr−1 which, within uncertainties, agrees with measurements and theoretical expectations (Shahzamanian et al. 2022; Schönrich et al. 2015; Sormani et al. 2022). The mean velocities ( and
) and the velocities’ dispersion (σb and σl ) for the Bulge population are in agreement with the values found in the literature (Soto et al. 2012; Kunder et al. 2012; Shahzamanian et al. 2022) and also theoretical studies (Sormani et al. 2022).
3.2. Brick field
In the Brick field, we proceeded as we did in the comparison field for the perpendicular component. Also, in this case, two Gaussians optimally describe the data (Fig. 4, bottom left and Table 3). These results agree with the values found for the comparison field for the perpendicular components.
For the parallel component, the three-Gaussians model does not converge to make a single solution. The best fit to the distribution corresponds to a two-Gaussians model (Fig. 4, bottom row and Table 3). The absence of a third Gaussian in the model could mean that the mean opacity on the Brick is too high to observe stars through it, blocking the light from the stars in the far side of the NSD, which agrees with the results found by Nogueras-Lara et al. (2021a). This is not surprising because the Brick is an extremely dense mid-infrared dark cloud (Longmore et al. 2012). Another possibility for the presence of only two Gaussians is that the number of detected sources through the Brick is not large enough to affect the distribution significantly.
As in the comparison field, in this case, with no sources moving westward (or just a small number of them), the reference frame drifted eastward, and in a more pronounced way because of the extreme extinction caused by the Brick. As a consequence, the mean velocity for the Bulge Gaussian is not around zero any more and the velocity of eastward moving stars is underestimated. Since we know that the mean velocity of Bulge stars should be around zero, we can compensate for the drifting effect by requiring the mean motion of the Bulge stars to be zero. That leaves us with a mean velocity for the parallel component of the NSD of μeast = 1.57 + 0.84 = 2.41 mas yr−1. This value agrees, within the uncertainties, with our results in the comparison field and with the rotation velocities for the disk found by Shahzamanian et al. (2022) and Schönrich et al. (2015).
4. Conclusions
We analyzed the stellar kinematics for two different regions, one on the dense molecular cloud G0.253+0.016 (the Brick) and another on a Comparison field westward from it, delimited with blue and green in Fig. 1. We found that, in both cases, the models that best describe the parallel and perpendicular components of the proper motion distributions are compatible with the presence of two different populations.
In the comparison field, for the perpendicular component of the velocity, we can optimally fit the proper motion distribution with two Gaussians (Fig. 4, top left), representing a Bulge and NSD population. We interpret the broader Gaussian to characterize the stars from the Bulge and the narrow one to correspond to the population of the NSD. Three Gaussians optimally describe the proper motion distribution for the parallel component: with one corresponding to the Bulge population and the other two corresponding to the NSD population, moving eastward and westward, respectively (Fig. 4, top right). In the Brick field, as for the comparison field, we clearly identify two different populations, which are represented with two Gaussians for the perpendicular component (Fig. 4, bottom left) and, because of extinction, also two Gaussians for the parallel component (Fig. 4, bottom right).
The velocity dispersion values for the Bulge populations in both fields are in agreement with the ones found by Soto et al. (2012), Kunder et al. (2012), and Shahzamanian et al. (2022). With respect to the parallel component, the values obtained in both fields are compatible with each other and also with the rotation velocity for the NSD found by Schönrich et al. (2015) and Shahzamanian et al. (2022).
If the Brick were located at a distance from the Sun of ∼7.20 kpc, as claimed by Zoccali et al. (2021), it would be placed well outside the NSD, but still inside the Bulge. Consequently, we would not expect to find traces of two different populations in the distribution, but only of one. Instead, we can clearly distinguish two different populations. Furthermore, the parameters of these distributions are in agreement with those found by other authors for the NSD and Bulge populations. These are strong indicators suggesting that the Brick is embedded in the NSD and the CMZ, and hence that it is located at the Galactic center.
This conclusion is supported by different arguments in other studies. For example, the turbulent velocity dispersion found for the gas in the Brick cloud does not match that found in other molecular clouds in the Galactic plane, which are significantly lower (see e.g., Henshaw et al. 2016, 2019).
Acknowledgments
We acknowledge financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). We acknowledge financial support from national project PGC2018-095049-B-C21 (MCIU/AEI/FEDER, UE). FN-L gratefully acknowledges support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project-ID 138713538 – SFB 881 (“The Milky Way System”, subproject B8), and the sponsorship provided by the Federal Ministry for Education and Research of Germany through the Alexander von Humboldt Foundation.
References
- Beroiz, M., Cabral, J., & Sanchez, B. 2020, Astron. Comput., 32, 100384 [NASA ADS] [CrossRef] [Google Scholar]
- Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, in Adaptive Optical Systems Technology, ed. P. L. Wizinowich, International Society for Optics and Photonics (SPIE), 4007, 879 [NASA ADS] [CrossRef] [Google Scholar]
- Gallego-Cano, E., Schödel, R., Nogueras-Lara, F., et al. 2020, A&A, 634, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Genzel, R., Eisenhauer, F., & Gillessen, S. 2010, Rev. Mod. Phys., 82, 3121 [Google Scholar]
- Ghez, A. M., Salim, S., Weinberg, N. N., et al. 2008, ApJ, 689, 1044 [Google Scholar]
- Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357 [NASA ADS] [CrossRef] [Google Scholar]
- Henshaw, J. D., Longmore, S. N., Kruijssen, J. M. D., et al. 2016, MNRAS, 457, 2675 [Google Scholar]
- Henshaw, J. D., Ginsburg, A., Haworth, T. J., et al. 2019, MNRAS, 485, 2457 [Google Scholar]
- Kruijssen, J. M. D., Longmore, S. N., Elmegreen, B. G., et al. 2014, MNRAS, 440, 3370 [Google Scholar]
- Kruijssen, J. M. D., Dale, J. E., & Longmore, S. N. 2015, MNRAS, 447, 1059 [Google Scholar]
- Kunder, A., Koch, A., Rich, R. M., et al. 2012, AJ, 143, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Launhardt, R., Zylka, R., & Mezger, P. G. 2002, A&A, 384, 112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Longmore, S. N., Rathborne, J., Bastian, N., et al. 2012, ApJ, 746, 117 [Google Scholar]
- Longmore, S. N., Bally, J., Testi, L., et al. 2013, MNRAS, 429, 987 [NASA ADS] [CrossRef] [Google Scholar]
- Matsunaga, N., Kawadu, T., Nishiyama, S., et al. 2011, Nature, 477, 188 [Google Scholar]
- McNamara, D. H., Madsen, J. B., Barnes, J., & Ericksen, B. F. 2000, PASP, 112, 202 [NASA ADS] [CrossRef] [Google Scholar]
- Nishiyama, S., Nagata, T., Kusakabe, N., et al. 2006, ApJ, 638, 839 [Google Scholar]
- Nogueras-Lara, F., Gallego-Calvente, A. T., Dong, H., et al. 2018, A&A, 610, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2019, A&A, 631, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2020, Nat. Astron., 4, 377 [Google Scholar]
- Nogueras-Lara, F., Schödel, R., Neumayer, N., & Schultheis, M. 2021a, A&A, 647, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nogueras-Lara, F., Schödel, R., & Neumayer, N. 2021b, A&A, 653, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schödel, R., Yelda, S., Ghez, A., et al. 2013, MNRAS, 429, 1367 [Google Scholar]
- Schödel, R., Feldmeier, A., Kunneriath, D., et al. 2014, A&A, 566, A47 [Google Scholar]
- Schönrich, R., Aumer, M., & Sale, S. E. 2015, ApJ, 812, L21 [Google Scholar]
- Shahzamanian, B., Schoedel, R., Nogueras-Lara, F., et al. 2022, A&A, in press https://doi.org/10.1051/0004-6361/202142687 [Google Scholar]
- Sormani, M. C., Magorrian, J., Nogueras-Lara, F., et al. 2020, MNRAS, 499, 7 [Google Scholar]
- Sormani, M. C., Sanders, J. L., Fritz, T. K., et al. 2022, MNRAS, 512, 1857 [CrossRef] [Google Scholar]
- Soto, M., Kuijken, K., & Rich, R. M. 2012, A&A, 540, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Speagle, J. S. 2020, MNRAS, 493, 3132 [Google Scholar]
- Tress, R. G., Sormani, M. C., Glover, S. C. O., et al. 2020, MNRAS, 499, 4455 [NASA ADS] [CrossRef] [Google Scholar]
- Zoccali, M., Valenti, E., Surot, F., et al. 2021, MNRAS, 502, 1246 [Google Scholar]
All Tables
All Figures
![]() |
Fig. 1. Pointings from GNS (outlined in black) and D19 (in white) overlaid on a 4.5 μm Spitzer/IRAC image (Schödel et al. 2014). The GNS pointings cover smaller areas because of detector windowing. The overlapping areas used for proper motion calculations are shown in blue (on the Brick) and green (comparison field). |
In the text |
![]() |
Fig. 2. Relative astrometric uncertainty as a function of H magnitude for all stars in the Brick and comparison fields that are common to D15 (red dots) and D19 (black dots). |
In the text |
![]() |
Fig. 3. Uncertainties for vx and vy as a function of H magnitude for stars with colors H − Ks > 1.3 on the Brick field. We considered only stars with velocity uncertainties smaller than 2 mas yr−1 (red dashed line). We did not consider stars with an absolute difference in magnitude with GNS bigger than 0.5 (red dots). We applied the same criteria in both fields. |
In the text |
![]() |
Fig. 4. Proper motion distribution for the comparison field (top row) and the Brick field (bottom row). Left-hand plots correspond to the perpendicular component, where the red Gaussian represents the Bulge stars and the black one corresponds to the NSD stars. Right-hand plots correspond to the parallel component, where the red Gaussian represents the Bulge stars and the yellow and black ones represent the NSD stars moving eastward and westward, respectively. |
In the text |
![]() |
Fig. 5. Extinction as a function of μl for the comparison field. Squares and bars represent the mean extinctions and their statistical uncertainties for stars in bins of 3 mas yr−1 width. Statistical uncertainties for μl are smaller than the width of the squares. |
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.