© ESO, 2012

1. Introduction

The signs that the Milky Way is a barred galaxy have increased following the many discoveries about the morphology of the central regions of the Galaxy, most coming from detailed near-infrared (NIR) star counts. A bar-like distribution in the stars of our Galaxy is observed from the asymmetries in the infrared surface brightness maps (e.g. Dwek et al. 1995) and in source counts (e.g. Stanek et al. 1994; Hammersley et al. 1994), which both systematically show more stars at positive Galactic longitudes within  < 30° and close to the Galactic plane compared to the negative side.

While some papers refer to a thick triaxial structure (bulge or thick bar) with a semi-major axis of around 2.5 kpc and a position angle (PA) of 15°−30° with respect to the Sun-Galactic centre direction (e.g. López-Corredoira et al. 2005; Habing et al. 2006; Rattenbury et al. 2007; Vanhollebeke et al. 2009), other contributions suggest that there is also a long thin bar, the in-plane bar, with 4 kpc in semi-major axis length and a position angle of around 45° (e.g. López-Corredoira et al. 2007; Cabrera-Lavers et al. 2007; Vallenari et al. 2008; Churchwell et al. 2009). This long bar has a tip in the positive longitude at the beginning of the Scutum arm (Dame et al. 1986). As realized by Sevenster et al. (1999), the angle is around 25 deg when low latitudes are excluded from the fit, and around 45° in the plane regions within −15° < l < 30°.

Rather than two structures, Martinez-Valpuesta & Gerhard (2011) suggested that there is a scenario of a single structure with a twisted major axis. We find this proposal acceptable from a purely morphological point of view. However, bars form as a result of an instability in differentially rotating discs (Sellwood 1981), whereas bulges are a primordial galactic component. Therefore, the scenario bulge+bar is not the same thing as a single structure, and one should observe other differences apart from the morphology. We think that the particular features of this new proposal of a “single twisted bulge/bar” scenario leaves some observational facts unexplained, whereas the model of a misaligned bulge+long bar successfully explains them all (Lopez-Corredoira et al. 2011). Furthermore, there are two distinct populations in the thick bulge. The metal-rich population presents bar-like kinematics, while the metal-poor population shows kinematics corresponding to an old spheroid or a thick disc, so the two main scenarios for the bulge formation co-exist within the Milky Way bulge (Babusiaux et al. 2010).

We pay attention here to the long bar. We explore it in the negative longitudes (−20° < l < −8°, |b| < 2°) with the data of 2MASS+VISTA-VVV in the NIR. For the positive longitudes there are available studies including deep NIR surveys (e.g. Cabrera-Lavers et al. 2007, 2008). Four aspects are interesting to explore in this region:

• 1.

  The stellar population of the bar itself, which can be followed by a tracer, such as the red clump population, as already carried out in the positive Galactic longitudes (Hammersley et al. 2000; Cabrera-Lavers et al. 2007, 2008).

• 2.

  The tip of the long bar in the negative longitudes, which is tentatively placed at l ≈ −14° (López-Corredoira et al. 2001). Finding there a massive star formation region would be new supporting evidence for the presence of a long bar in our Galaxy, as it was for the tip of the bar in the positive longitudes (López-Corredoira et al. 1999; Negueruela et al. 2011; Green et al. 2011; Nguyen Luong et al. 2011). Such regions can form because of the concentrations of shocked gas where a galactic bar meets a spiral arm, as observed at the ends of the bars of face-on external galaxies.

• 3.

  The dust lane before the long bar in the negative Galactic longitudes (Calbet et al. 1996; Marshall et al. 2008).

• 4.

  The stellar density of the inner disc around and behind the long bar+bulge, and particularly the possibility that the in-plane regions have a deficit of stars, due to a truncation in this density or to the effect of a flare in the disc (López-Corredoira et al. 2004).

All these features will be analysed using NIR colour–magnitude diagrams of the stellar populations, deriving the extinction, and separating within them the red clump stars (old population) and the red supergiants or very bright red giants (young population).

This paper is linked to Amôres et al. (2012), where the stellar counts produced with the same datasets used here in the form of CMDs are analysed. After a brief introduction of the datasets that will be used, a detailed description of the analysis tools we used follows in Sect. 3. In Sect. 4 the main results of the paper are given, leading to a summary of the more relevant conclusions in Sect. 5.

2. Data

2.1. VISTA-VVV

The ESO public near-IR variability survey VVV is scanning the Milky Way bulge and an adjacent section of the mid-plane with high star formation activity (Minniti et al. 2010). The survey will take 1929 hours of observations with the four metre VISTA telescope, across an area of approximately 550   deg in the inner galaxy. The final products will be a deep IR atlas in five passbands (ZYJHK_S) and multi-epoch photometry in the K_S-band to build a catalogue of about a billion point sources, including about a million variable sources. The main goal of the VVV survey is to enable the construction of a 3D map of the inner Milky Way using well-understood distance indicators, such as RR Lyrae, Cepheids and red clump giants.

The first dataset, including observations from year 2010, has just been released. Saito et al. (2012) discuss the observations, reductions, photometry, calibrations, and catalogues. For most regions of interest the limiting magnitude is m_KS ~ 18.0, and the image scale is 0.34′′/pixel. A more detailed description of the VVV data used can be found in Amôres et al. (2012).

2.2. 2MASS

Although we are interested in the inner Milky Way, best viewed with VVV, we also study brighter populations, some of which might be too bright for this survey. We supplement the VVV data using 2MASS All-Sky Release (Skrutskie et al. 2006), accessing its point source catalogue through the Gator interface¹.

3. Data analysis

3.1. Combining 2MASS and VVV

We have made intensive use of the colour–magnitude diagram (CMD) instead of a source-by-source analysis, so a detailed combination of both datasets is not needed. Instead, for a more subtle approach we opted to adopt hard thresholds in K_S and J: above these thresholds we rely on 2MASS photometry, and below them VVV is the data source.

Determining these thresholds requires careful calculation of the completeness of 2MASS. We restricted ourselves to −20° < l < −8°, and in these regions completeness is dominated by source confusion, owing to the high stellar density. For each of the fields we built a histogram in magnitude, fit the brighter part with a second degree polynomial, and checked when the observed counts fall below 80% of the predicted value offered by this fit, adopting this value as our completeness limit. As we want to have a homogeneous set of CMDs, we selected a single set of thresholds for all our fields: (J,K_S) = (14,12.5).

As we are comparing data from two different surveys, with different instrumental configurations and calibration procedures, we also checked that both photometries are consistent. As both catalogues overlap in magnitude, for each field we crossmatched the common stars and checked for differences in magnitude. The typical values are around 0.05   mag in J and 0.02   mag in K_S, although in some fields subject to severe reddening these values can be higher, reaching differences up to 0.4   mag. In all our subsequent analysis, these systematic differences have been added to VVV magnitudes.

We also checked the behaviour of the magnitude difference against colour and found no evident dependence, but in the redder part of the CMD (with J − K_S > 4.5), there are few common stars, so this effect could be masked.

3.2. The infrared colour magnitude diagram

Since the advent of large-scale infrared surveys, the infrared CMD, usually of (J − K_S) versus K_S has become one of the standard tools for studying the structure of the Milky Way. Beyond the penetrating powers of this wavelength range that allow us to reach the inner parts of our Galaxy, in these diagrams it is relatively easy to separate the two main visible populations, dwarfs and giants. Although the morphology of the CMD changes greatly in different areas of the Galaxy, since they are heavily affected by interstellar extinction, in-plane fields offer a similar aspect, exemplified in Fig. 1.

This CMD was obtained for a field centred at (l,b) = (−11°,0°). The main features that dominate the figure are outlined in the right-hand panel: in blue and occupying the lower left-hand region of the plot, the dwarfs of the disc are the more abundant population, separated from the giants due to the extinction and their intrinsic brightness. Even if caused by distance, dwarfs and giants could appear to us with the same apparent magnitudes, since giants are intrinsically brighter they would have to be farther away and be the subject of heavier interstellar extinction, since we are dealing with in-plane fields. This implies that they will appear shifted towards the red, occupying then the top right-hand part of the diagram.

The effect of distance (that moves stars towards dimmer apparent magnitudes) and extinction (that makes stars appear dimmer and redder) results in an ideal population of stars with the same intrinsic magnitude and colour and evenly distributed along the line-of-sight that would appear as a diagonal line (not necessarily straight) in this CMD. The natural spread in magnitude of any real population changes this line into a strip of a certain width, and the uneven stellar density of the Galaxy translates into overdensities in the (J − K_S,K_S) space. Because the later the type of a giant the brighter and redder it becomes, the lower envelope of this region is composed by late type G giants, but the bulk of it will correspond to giants of spectral types around K2, also known as red clump giants (RCG), which are the most abundant giant population. The redder tail of Fig. 1 would be composed almost exclusively of M giants. The relative abundance of RCGs is high enough so that they appear as a distinct strip on the CMD while later type giants, much less abundant, are visible towards redder and brighter magnitudes.

To build the CMDs used in this work (including the one from Fig. 1), we always use 0.5° × 0.5° fields centered on the nominal coordinates. As the source density is high, instead of plotting star by star, we built a density diagram using tiles of 0.1 mag in colour and magnitude, as can be seen in Fig. 1. Also, the variation range of N, the number of stars in each one of these tiles, is wide, so that to appreciate all the substructure present in these CMD it is mandatory to plot instead some function that compresses this range. After a lengthy trial and error process we opted for the function (1)where N is the number of stars per tile and I the plotted intensity. This selection combines the compression factor of a logarithmic intensity plot but without losing the ability to show density variations of  ~100 stars per tile, changes that in a pure logarithmic plot would be invisible. As our analysis of these density variations was purely qualitative, we opted to use a function that maximizes the visible information instead of one with immediate physical translation, such as a square root or a logarithm.

The resulting CMDs, used in subsequent sections, can be seen in Appendix A.

3.3. The red clump method

As said above, the ubiquity of the RCGs makes the so-called red clump method one of the standard tools for analysing CMDs in the infrared. Its use and application is discussed in depth in López-Corredoira et al. (2002, hereafter LC02), but in this section we offer a brief outline of its application to the 2MASS+VVV composite CMDs, as the results of this method is used thoroughly through the text.

These RCGs have a very narrow luminosity function that it is also calibrated well (LC02), so their intrinsic colour and luminosity are well known. It is possible then to transform their position in the (J − K_S,K_S) plane to (d,A_K), once we assume an extinction law (in the form of the relative extinction coefficient A_J/A_K). A more detailed description of the method and the errors involved can be found in Cabrera-Lavers et al. (2008) and references therein. Following this study, we use (J − K_S)₀ = 0.70 ± 0.05, M_K = −1.62 ± 0.03 and A_J/A_K = 2.518.

To perform this calculation there is a previous step that needs to be taken: the calculation of the (J − K_S,K_S) curve that represents the distribution of RCGs over the CMD. The simplest way to do this is to cut the CMD into several strips in K_S and then fit the (J − K_S) distribution of each strip with the sum of a second degree polynomial and a Gaussian; the μ of this Gaussian gives the fiducial point that represents the locus of the RCG for that strip (LC02). The problem here is that, in the presence of high interstellar extinction (see for example the CMDs for b = 0° in Fig. A.1), the RCG strip almost appears horizontal, and so projecting the data of each strip over the (J − K_S) is not optimal, because it mixes stars with very different colours but with roughly the same K_S magnitude, resulting in a loss of information.

To overcome this, as a first step in our analysis, a rough fit of the RCG strip is performed by eye, and this (J − K_S,K_S)₀ curve is used as the entry point of an automated process that, instead of performing cuts over K_S and projecting them over the (J − K_S) axis, defines a local coordinate system (X,Y) that is perpendicular to this (J − K_S,K_S)₀ curve at each point. As these new axes take the reddening vector into account, stars with similar magnitudes but different reddening appear separated when projecting over either X or Y, so a much more refined (J − K_S,K_S) curve can be derived.

Once we have the variation in (J − K_S) with K_S for the RCGs, we can extract this population from the CMD. Although this can be optimized for each line-of-sight, and because our study requires the comparison of data at different Galactic longitudes, we consider that everything within 0.5   mag in colour and magnitude from this curve is a RCG.

3.4. Photometric search for massive stars

One of the best tracers of stellar formation regions are red supergiants (RSG). Since these stars are very bright in the infrared (with a conservative limit of M_K < −8, Wood et al. 1983), they shine through the obscuration present in the plane. One of the expected effects of a galactic bar is a triggering of star formation at its tips and, depending on the bar strength, over its dust lanes (see, for example, Sheth et al. 2002). As such, the study of the distribution of RSGs with Galactic longitude could give us clues on the structure of the bar.

In the near end of the bar, around l ~ 28° several clusters of these stars have already been detected (i.e. Negueruela et al. 2011; González-Fernández & Negueruela 2012). In these clusters, RSGs will have apparent magnitudes 7 < m_k < 4. If we assume that the bar forms a 45° angle with the Sun-Galactic centre line and has a 4   kpc semimajor axis (Cabrera-Lavers et al. 2008), this would put the far end of the bar at  ~11.2 kpc, while the near one sits at  ~5.9 kpc. This difference, in distance modulus, is of  ~1.4   mag, so the same RSG that we see at the near end of the bar would have magnitudes in the range 9 < m_k < 6, allowing for a  ~0.5   mag in extra extinction due to the presence of the dust lane of the bar.

It is worth noting that, although we have selected a φ = 45° bar, other possible position angles have been postulated (for example, Babusiaux & Gilmore 2005, with φ = 22° ± 5.5°), but always with φ < 45°. This implies that for these bar configurations, the far end should be farther away than for a 45° angle bar. The difference in brightness for a given RSG between both configurations is only of 0.1   mag (corresponding to 0.6   kpc), and our selection criteria for RSG is therefore not biased against these models.

For this reason, the possible RSGs situated at the far end of the bar would be very difficult to separate from the population of late type giants (they also have very similar intrinsic colours) of the disc. One possible way to separate this population is to use the reddening free pseudo-colour Q, defined as Q = (J − H) − 1.8 × (H − K_S) (Negueruela & Schurch 2007). Generally speaking, the later (and intrinsically redder) a star is, the higher value of Q it will have, with late type giants (types beyond K0) having Q > 0.41. RSGs, because their atmospheres are very extended, complex, and often present mass loss, will not deredden properly, and so they roughly occupy the range 0.1 < Q < 0.4 (Messineo et al. 2012).

Another applicable threshold affects the colour of the candidates. RSGs are intrinsically red sources, with intrinsic colours (J − K_S)₀ > 0.7. Since we are also looking at sources under heavy extinction (as is expected for the inner Milky Way), we can safely assume that any RSG will have (J − K_S) > 1.3. This step eliminates late type dwarfs from the disc that will fulfill most of the previous criteria but, as they are located on the local disc will have lower colour excesses.

Summing up, the criteria that we used to select RSG candidates over 2MASS data are

• 1.

  0.1 ≤ Q ≤ 0.4,

• 2.

  m_K ≤ 9,

• 3.

  (J − K_S) ≥ 1.3.

Negueruela et al. (2011) performed a spectroscopic survey around the massive cluster RSGC3, confirming the nature of stars selected based on these criteria, and find that the initial sample is entirely composed of RSGs and very evolved giants, with stellar types later than M5. A 38% fraction of the objects turned out to be class I sources (supergiants), while 28% were class II (bright giants), and the remaining 34% identified as late class III giants.

4. Results

4.1. Structure of the Galactic bar

One of the most striking facts when analysing the CMDs produced with 2MASS+VVV for fields with |l| < 0° and b = 0° is that there is a sudden drop in the number of giants around l = −13° (see Fig. 2). Since the total number of sources per square degree would be dominated by the dwarf population of the disc and we are more interested in the innermost region of the Milky Way, we focus our study on the giants, that reach beyond the disc within the limiting magnitude of VVV.

Using the method described in Sect. 3.3, we can extract the RCGs from a CMD and estimate their distance. This allows us to reconstruct the in-plane 2D distribution of these stars, as can be seen in Fig. 3. Beyond l = −8°, where extinction and source density affect our sample, the behaviour of the number of giants with (l,d) is rather smooth, showing traces of a rather straight structure that ends beyond l ~ −12°, where there is an appreciable drop in the number of visible giants. This could easily be associated with a bar-like structure, but before doing this there are some issues that need to be tackled. This drop in the number of giants visible in Fig. 2 could be due to an increment in the extinction along the line-of-sight, but this is unlikely: in this same figure it can be seen that the morphology of the RCG strip is similar in all the lines of sight, with higher latitudes and lower extinction. This implies a similar (d,A_K) behaviour for all of them. Using the method described in Sect. 3.3, we can estimate the variation of E(J − K_S) with d, plotted in Fig. 4. Comparing it with Fig. 2 we see that, although the extinction along the line-of-sight is almost equal for both of them, the number of visible giants is appreciably lower for l = −13°.

The density distribution of Fig. 3 could also be an effect of the variations in the limiting magnitude with Galactic longitude. A simple exponential density distribution, when observed with the instrumental setup of VVV, would show a maximum in the number of stars near completeness, for when we approach this limit we start to gradually lose more visible stars, modulating the subjacent stellar distribution. As hinted before, in off-plane fields these effects are less severe, and in Fig. 5 we can see how this bar-like structure is visible, well before the completeness limit is reached.

This structure corresponds very well to the extension of the Long Bar, as has been observed previously in integrated light (Hammersley et al. 2000) and stellar counts (López-Corredoira et al. 2001), to the fourth quadrant. Assuming that this bar is symmetric with respect to the Galactic centre, the positions of its tips provide a strong constraint on its size and orientation. We place its far end at l = −13° (as suggested by Figs. 3, 9, and 5), and the nearby end in l = 27° (Cabrera-Lavers et al. 2008). Although these measurements are not very precise, the error should not exceed  ± 1°. Assuming that R_GC = 8.0 ± 0.5   kpc, we obtain an angle of 41° ± 5° and a total length of 8.2 ± 0.6   kpc. This is close to the values calculated by López-Corredoira et al. (2001); Cabrera-Lavers et al. (2007, 2008), but clearly different from the proposed geometry of Babusiaux & Gilmore (2005).

The structure in Figs. 3 and 5 seems to have a wider angle, close to 45° but the range in Galactic longitudes is small, and these differences in orientation are difficult to judge. Also, the derivation of d for this figure is affected by extinction, which could warp the visible structure if there is variation between different lines of sight. The properties derived using only the angular positions of the bar ends are very robust, not affected by any of these effects, the distance to the Galactic centre being the only assumption made.

Another interesting feature that, for example, can be observed in the CMDs of Fig. 2 and the extinction curves from Fig. 4 is that there is a huge amount of interstellar material concentrated in a narrow strip around 9   kpc at l = −13° (although this distance estimate varies with the assumed extinction law, see Sect. 4.3). This wall of extinction makes the RCGs appear as an almost horizontal strip on the CMD, but its effect is constrained to the plane. In Fig. 6 we plot the diagrams of a vertical cut at l = −12.5° (near the end of the bar), and we can see clearly how the effect of interstellar extinction decreases from b = 0.0° to b = −0.5°. This feature of the interstellar medium can be clearly followed in the plane from l = −13° to l = −8° (see Figs. A.1 and A.2).

In these longitudes it is easy to see that this structure is too far away to be mistaken for the 3-kpc arm/molecular ring (contrary to what might happen in the first quadrant), so this sudden increase in the extinction (hence, the density of interstellar material) is a clear clue to the existence of a leading dust lane in the Long Bar of the Milky Way. A similar structure was previously detected in the inner Milky Way over extinction maps and molecular emission (Marshall et al. 2008), although at lower |l|. Dust lanes are a common companion to stellar bars in other galaxies (for example, Comerón et al. 2009), so the presence of such a structure in our Galaxy is a hint of the existence of a bar up until l ~ −12°.

As we move outwards, this dust lane disappears, but the interstellar extinction grows more or less continuosly for l < −14.5°. As can be seen in Fig. 7, the variation in E(J − K_S) with d is quite different from that of the dust lane. Instead of a sudden increase in the amount of interstellar material, the variation of Fig. 7 is much smoother. This matches very well what is expected of a near-tangent line-of-sight to a gas rich structure such as an arm (see Sect. 4.4).

4.2. Structure of the inner disc

A careful look at Fig. 5 presents us with another interesting feature: the depletion of giants present from −18° ≤ l ≤ −12° (particularly evident at b = 1.5°). This truncation of the stellar density of the disc is visible in every Galactic latitude except for |b| ≤ 0.5°, due to completeness (see Fig. 3). This is one of the strongest visible effects of a bar; its typical potential (Binney & Tremaine 2008) has two maxima in the direction of its minor axis (labelled L₄ and L₅ in Fig. 8, following the notation of Binney & Tremaine 2008). These maxima sweep away stars and gas outwards, effectively emptying the disc inside the corotation radius. In Fig. 5 this appears clearly in the first plot, and there are traces of the same effect in the third, but masked mostly by the effects of extinction. The method described in Sect. 3.3 assumes that the extinction varies smoothly along the line-of-sight, so the clumpiness of the real distribution of A_K with d translates into a cumulative blur of the features in Figs. 3 and 5.

This truncation is not only visible in the stellar content, but it also affects interstellar matter. As can be seen in Fig. 4, from  ~4.5 to  ~9   kpc there is little variation in E(J − K_S) with d. As the colour excess is proportional to A_KS and this quantity, in turn, is a function of the integral of the volume of interstellar matter along the line-of-sight, a flat curve implies that there is little to no underlying gas, matching the behaviour of the giants in Fig. 5.

Even so, the photometric depth of VVV is high enough to probe the far disc for the first time. In some of the CMDs of Appendix A, the trace of the RCG runs vertically for m_K ≥ 12 (particularly evident for (l,b) = (−12.5, −1.0)). This trace shows a secondary maximum, visible on several lines of sight (as, for example, in Fig. 6). These maxima appear around 13   kpc, and are clearly visible in the reconstructed distribution of Fig. 5 (more evident in the middle and bottom panels).

This increase in stellar density shows us where the disc regains its density beyond corotation. This effect has also been noted in Gonzalez et al. (2011) at |b| = 1°, although the authors interpret these maxima as traces for an arm at 11.2   kpc. This seems unlikely, since a one-degree separation from the plane implies a height of 200   pc at 11   kpc (300   pc for a 1.5° line-of-sight), and arms are structures enclosed in the thin disc, so have very small height-scales. On top of this, the effective volume of an arm that a line-of-sight crossing it probes is very small unless at a tangent point. Since these structures are mostly composed of young populations, they are not traced very well with RCGs (a late population, with τ ~ 2   Gyr), and it is very unlikely that they could appear clearly at any line-of-sight, as in Fig. 5.

4.3. The behaviour of interstellar extinction

Many of the previous results depend, or at least are affected by, extinction. In the derivation of Figs. 3 and 5 we have assumed a simple extinction law that remains constant with d, l, and b and that yields A_J/A_K = 2.518. As can be seen in Fig. 4, beyond 9   kpc, the slope of the (d,E(J − K_S)) curve becomes negative. This is physically impossible, and is a direct effect of this choice of extinction law, for the derivation of d depends on the calculation of A_K and, as follows: This implies: Once we isolate a given population of stars that we will use as standard candles, from two measurables, (J − K_S) and m_KS, we therefore derive d, E(J − K_S), and A_KS. To complete the recipe we need to know, apart from the aforementioned A_J/A_KS, the luminosity function of the population under study, through its (J − K_S)₀ and M_KS. While a variation in any of these three constants would change Fig. 4, we know that:

• An offset in (J − K_S)₀ or M_KS would, in exchange, offset the curves in Fig. 4, but not change their slope or overall shape.

• For the case of the RCGs, there is a known dependency of (J − K_S)₀ with metallicity (Cabrera-Lavers et al. 2007), but we can see in Fig. 4 that, since the extinction increases in such a quick manner, most of the variation (particularly the part in which the curve acquires negative slope) takes place in less than a kiloparsec, so the metallicity spread is expected to be low. For this reason, while the overall shape of the d vs. E(J − K_S) curve might change due to this effect, this variation would affect the slope of the nearby part of the graph with respect to the points that are more distant from the observer, but it will not change the slope measured for d > 8   kpc.

Taking these points into consideration, the only remaining possibility to explain this negative slope is our choice of a value of A_J/A_KS not representative of the observed physical conditions. We can tinker with this value in order to bring the curves in Fig. 4 into coherence with the expected behaviour of extinction with distance (that, by definition, needs to be a monotonically increasing function), obtaining a lower limit for the observed value of this ratio. We find that so far away from the typical values of Rieke & Lebofsky (1985) or Cardelli et al. (1989), both close to A_J/A_KS = 2.5. This points to an interstellar material that is more transparent towards the red than predicted by the extinction law, more in concordance with the results from studies as Nishiyama et al. (2009) and Stead & Hoare (2009). This effect is not exclusive to the lines of sight used to produce Fig. 4; heavy extinction transforms the RCG strip in a CMD into a nearly horizontal band, so similar behaviors are observed in all the fields with −8.0° < l < −13° and b = 0.0° in Figs. A.1 and A.2.

A change in the A_J/A_KS extinction ratio implies an expansion/contraction in the distance scale in Figs. 3 and 5, but will not change the overall shape of the structural features discussed. Because the analysis of these figures and the results drawn from are qualitative and the physical characteristics of the bar are calculated from the Galactic longitudes of its ends, a measure that is independent of extinction, our results are not affected by the choice of extinction law. Since a proper determination of the value of A_J/A_KS is beyond the scope of this work, we opt to retain the value from Rieke & Lebofsky (1985) for homogeneity with similar studies (for example, Cabrera-Lavers et al. 2007).

This distance scale factor introduced through A_J/A_KS can be appreciated when comparing Fig. 3 with 5: the bar overdensity seems to shift away from the observer at higher |b|. As increasing A_J/A_KS implies decreasing A_KS and a higher derived d, the error in this last parameter due to the difference in A_J/A_KS will be a monotonically increasing function of E(J − K_S), hence higher for in-plane fields.

4.4. Bar-driven stellar formation in the Milky Way

A structure such as the Long Bar is expected to affect its environment in several ways. One of them is the triggering of star formation where its tips interact with the Galactic disc. Since RSGs are massive stars, their lifespans are expected to be short, just a few tens of Myr, so they trace recent stellar formation. We built a sample of photometrically selected RSG candidates following the methods discussed in Sect. 3.4; namely, we selected stars with 0.1 ≤ Q ≤ 0.4, (J − K_S) > 1.3 and m_K < 9. Since these values are well within the completeness of 2MASS, we use this survey. Doing so lets us obtain candidates for all the inner Galaxy, from l = 40° to l = −30°.

The vertical distribution of RSG candidates is, as expected for a young population, concentrated in low b fields, and it can be roughly approximated with a Gaussian of σ = 2.4°. We find that in the region with −30° ≤ l ≤ 40° there are close to 30   000 candidates with |b| ≤ 1° (roughly a density of 200   deg^-2) and 60   000 in off-plane fields.

As reasoned in Cabrera-Lavers et al. (2007), the long bar is a very flat structure, so to locate traces of its associated stellar formation it is natural to restrict ourselves to the Galactic plane. Since several structures coexist in the region with |l| < 30°, we compared the in-plane (| b| ≤ 1°) distribution of candidate RSGs with off-plane (1° < |b| ≤ 90°) fields. The results can be seen in Fig. 9 (to facilitate comparison, the off-plane distribution has been normalized to have the same total number of objects as for in-plane fields). As can be seen, rather than being smooth, the in-plane distribution presents a series of overdensities not present in off-plane fields (or present in a much less conspicuous manner).

These overdensities (outside the Galactic center) give us hints to possible underlying structures. This distribution follows the 2.2   μm and 12   μm emission studied by Hammersley et al. (1994) closely. As the structure of the first quadrant of the Galaxy is much better studied, the nature of some of these peaks has already been established: the overdensity at l = 32° is related to the tangent line-of-sight towards the Scutum-Crux arm, and its counterpart in the fourth quadrant was already established by Hammersley et al. (1994) beyond l < −29° (hence outside our analysis). The next peak at l = 28° is known to be related to one of the more massive episodes of star formation in our Galaxy, containing several massive clusters with M_⊙ > 10⁴   M_⊙, and probably with a combined mass well in excess of 10⁵   M_⊙ (González-Fernández & Negueruela 2012; Negueruela et al. 2011). Since the wings of this peak are wide, it also contains the extremely massive clusters RSGC1 (Figer et al. 2006; Negueruela et al. 2010) and Stephenson 2 (Davies et al. 2007). All these clusters have been related to the interaction of the Long Bar and the Scutum-Crux Arm/molecular ring; this area also contains the more massive star forming region of the Milky Way, W43 (Nguyen Luong et al. 2011, although it does not harbour any RSG yet).

It is expected that a similar wealth of structures should occur at the far end of the bar, here identified at l = −13°, although the associated stellar formation has yet to be discovered. Some recent studies of massive clusters have started this search, but so far the only confirmed massive association detected is Mercer 81 (Davies et al. 2012). It is associated by the authors to the far tip of the bar, but sitting at l = −21.6° is likely to be part of the 3-kpc arm, with a tangent at l = −22° (LC01). This is a prominent feature in gas emission that was the first indication of our Galaxy being barred (de Vaucouleurs 1964). Although it was previously suspected that this arm has not been subject to recent stellar formation (Dame & Thaddeus 2008), in Fig. 9 it appears clearly delineated through its tangent line-of-sight, pointing towards an associated RSG population.

This arm occupies roughly the same space as the so-called molecular ring, a feature visible mostly in molecular emission (Jackson et al. 2006). Both structures are mostly young and will be poorly traced by RCGs. As a result, only their tangents are visible in our data and, it is visible in Fig. 5, is very difficult to distinguish between the two possible configurations, a closed ring or two tightly wound arms, by looking at both our CMDs and RSG distribution.

The exact stellar nature of the remaining overdensity in the near end, visible in Fig. 9 spanning from l ~ 20° to l ~ 15°, has not been established yet. It might be part of the Scutum-Crux arm, might be stellar formation of the bar itself ², or a part of the molecular ring or the 3-kpc arm (Hammersley et al. 1994) (or any combination of the three). It does not have a clear counterpart at the far end, and this could point to it being part of the 3-kpc arm. As can be seen in Fig. 8, the cross section offered by the far pass of the arm to the line-of-sight is much wider, so all the stellar formation that in the first quadrant is visible in 20° ≤ l ≤ 10° would appear, due to the effect of projection, around the tangential line-of-sight towards the arm in the fourth quadrant (at l = −22°). If there is any stellar formation along the dust lane, in Fig. 9 it would be visible with l ≤ −13° and so it would be mixed with that of the tip of the Long Bar and the increased density towards the Galactic centre. As such, it is almost impossible to distinguish between both hypotheses without a proper spectroscopic study of the individual sources themselves.

5. Conclusions

In this work we obtained combined CMDs of 2MASS and VVV, mapping for the first time the innermost Milky Way at negative Galactic longitudes. Using these data along with photometrically selected RSGs – used as proxies for recent stellar formation – we mapped the structures present at these latitudes, offering a fresh picture of the structure of our Galaxy.

The main conclusions drawn follow:

• Based on the analysis of the CMDs and the 2D distribution of the RCGs along several lines of sight, we showed there is a flat structure that runs up to l ~ −12°. It corresponds to the extension of the Long Bar into the fourth quadrant. With this new constraint, we concluded that it forms an angle of 41° ± 5° with the Sun-Galactic centre line.

• This Long Bar affects the structure of the disc. One of the more severe effects is the truncation of the stellar density with R_GC ≤ 5   kpc. Because it corresponds with two maxima of the bar potential, this area is devoid of stars and gas (outside the body of the bar itself).

• Using photometrically selected RSGs as a proxy for recent stellar formation, we demonstrated another of the predicted effects of a galactic bar, the triggering of formation regions at its ends, clearly visible at l ~ 27° and l ~  −12°. The tangents to the Scutum-Crux and the so-called 3-kpc arm are also visible in the distribution of these stars, showing that in both arms there are traces of recent massive stellar formation.

---

Acknowledgments

This research is partially supported by the Spanish Ministerio de Ciencia e Innovación (MCINN) under AYA2010-21697-C05-5 and the Consolider-Ingenio 2010 Program grant CSD2006-00070: First Science with the GTC (http://www.iac.es/consolider-ingenio-gtc). M.L.C. was supported by the grant AYA2007-67625-CO2-01 of the Spanish Science Ministry. Eduardo Amôres obtained financial support for this work from the Fundação para a Ciência e Tecnologia (FCT) under the grant SFRH/BPD/42239/2007 and CNPq (311838/2011-1). We gratefully acknowledge use of data from the ESO Public Survey programme ID 179.B-2002 taken with the VISTA telescope, data products from the Cambridge Astronomical Survey Unit, and funding from the FONDAP Centre for Astrophysics 15010003, the BASAL CATA Centre for Astrophysics and Associated Technologies PFB-06, the MILENIO Milky Way Millennium Nucleus from the Ministry of Economy’s ICM grant P07-021-F, and by Proyecto FONDECYT Regular No. 1090213 from CONICYT. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the NASA/ IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research has made use of Aladin.

References

Online material

Appendix A: VVV+2MASS CMDs of the inner Milky Way

All Figures