Gas Flow Models in the Milky Way Embedded Bars
Total Page:16
File Type:pdf, Size:1020Kb
Astronomy & Astrophysics manuscript no. paper_dyna_2 c ESO 2018 October 31, 2018 Gas flow models in the Milky Way embedded bars N. J. Rodriguez-Fernandez1 and F. Combes2 1 IRAM, 300 rue de la Piscine, 38406 St. Martin d'Heres, France e-mail: [email protected] 2 LERMA, Observatoire de Paris, 61 Av de l'Observatoire, 75014 Paris, France Received September 15, 1996; accepted March 16, 1997 ABSTRACT Context. The gas distribution and dynamics in the inner Galaxy present many unknowns as the origin of the asymmetry of the lv-diagram of the Central Molecular Zone (CMZ). On the other hand, there are recent evidences in the stellar component of the presence of a nuclear bar that could be slightly lopsided. Aims. Our goal is to characterize the nuclear bar observed in 2MASS maps and to study the gas dynamics in the inner Milky Way taking into account this secondary bar. Methods. We have derived a realistic mass distribution by fitting the 2MASS star counts map of Alard (2001) with a model including three components (disk, bulge and nuclear bar) and we have simulated the gas dynamics, in the deduced gravitational potential, using a sticky-particles code. Results. Our simulations of the gas dynamics reproduce successfully the main characteristics of the Milky Way for a bulge orientation of 20◦ − 35◦ with respect to the Sun-Galactic Center (GC) line and a pattern speed of 30-40 km s−1 kpc−1. In our models the Galactic Molecular Ring (GMR) is not an actual ring but the inner parts of the spiral arms, while the 3-kpc arm and its far side counterpart are lateral arms that contour the bar. Our simulations reproduce, for the first time, the parallelogram shape of the lv-diagram of the CMZ as the gas response to the nuclear bar. This bar should be oriented by an angle of ◦ ◦ 9 ∼ 60 − 75 with respect to the Sun-GC line and its mass amounts to (2 − 5:5) 10 M . We show that the observed asymmetry of the CMZ cannot be due to lopsidedness of the nuclear bar as suggested by the 2MASS maps. Conclusions. We do not find clear evidences of lopsidedness in the stellar potential. We propose that the observed asymmetry of the central gas layer can be due to the infalling of gas into the CMZ in the l=1.3◦-complex. Key words. Galaxy: structure { Galaxy: center { Galaxy: kinematics and dynamics { ISM: kinematics and dynamics { Methods: numerical 1. Introduction Gaussians functions. Since COBE, the direct evidence of the bar has increased based on photometric surveys and Studying the structure and the dynamics of the Milky Way star counts in the near-infrared (see for instance L´opez- galaxy is necessary to understand its formation and evolu- Corredoira et al. 2000, 2005; Babusiaux & Gilmore 2005, tion. For example, the large scale star formation rate, the and references therein). Although, there is a consensus arXiv:0806.4252v1 [astro-ph] 26 Jun 2008 gas transfer to the nucleus and the fueling of the central on the triaxial nature of the bulge, the actual shape super-massive black hole are strongly influenced by the and the inclination is still under debate with currently large scale dynamics of the Galaxy. no clear agreement among different authors. The inclina- At least three components are needed to explain the tion angles found range from ∼ 10◦ − 15◦ (Binney et al. rotation curve of the Galaxy: a disc, a bulge and a mas- 1991; Freudenreich 1998; Lopez-Corredoira et al. 1997) sive dark halo. In near-infrared images, the bulge is a tri- to ∼ 30◦ − 45◦ (Weiner & Sellwood 1999; Sevenster 1999; axial structure difficult to distinguish from a bar. The L´opez-Corredoira et al. 2001) with a number of works giv- first direct evidence of the barred nature of the Milky ing values in the range 15◦ − 30◦ (Mulder & Liem 1986; Way was found by Blitz & Spergel (1991). The analysis Binney et al. 1997; Fux 1999; Englmaier & Gerhard 1999; of COBE/DIRBE data (Dwek et al. 1995; Binney et al. Bissantz & Gerhard 2002; Bissantz et al. 2003; Babusiaux 1997; Freudenreich 1998) supported Blitz & Spergel re- & Gilmore 2005; L´opez-Corredoira et al. 2005). sults. Dwek et al. (1995) found that the peanut-shaped bulge seen in the COBE image is best fitted by boxy In addition to the bar-like boxy-bulge, some recent works have also proposed the existence of another bar that Send offprint requests to: N. J. Rodriguez-Fernandez would be longer and thinner than the bulge (Picaud et al. 2 N. J. Rodriguez-Fernandez and F. Combes: Gas flow models in the Milky Way embedded bars 2003; L´opez-Corredoira et al. 2001, 2007; Benjamin et al. To our knowledge, the works by Binney et al. (1991), 2005). As the bulge, the near side of this long and thin bar Jenkins & Binney (1994) and Lee et al. (1999) are the will also be located in the first Galactic quadrant. In con- only dynamical models that have tried to interpret the trast, the long bar would be oriented at an angle of ∼ 45◦ dynamics of the complex Galactic center region (the cen- with respect to the line of sight towards Galactic center. tral 500 pc). Other models of the Galactic dynamics have In the present paper we will not discuss this component been mainly devoted to larger scales, from the works of as the mass should be dominated by the bulge. Hereafter, Mulder & Liem (1986) and Wada et al (1994) to the more we will refer to the triaxial boxy bulge as \large bar" or recent works by Fux (1999), Englmaier & Gerhard (1999), simply \bulge". Weiner & Sellwood (1999) and Bissantz et al. (2003). The Due to the difficulty to measure the stellar veloci- bar pattern rotation speed for most of those models are ties, most of our knowledge of the Galactic dynamics in the range from 40 to 60 km s−1 kpc−1. Regarding the come from spectroscopic observations of the atomic and corotation radius, it is estimated to be in the range 3-4 molecular interstellar medium (Rougoor & Oort 1960; kpc by some works (Englmaier & Gerhard 1999; Bissantz Scoville 1972; Liszt & Burton 1978; Bally et al. 1987; et al. 2003; Habing et al. 2006) and in the range 4-5 kpc Dame et al. 2001). The gas exhibit large non-circular mo- by other works (Combes 1996; Sevenster 1999; Fux 1999). tions in the inner Galaxy. Among the possible explana- In contrast, Binney et al. (1991) placed the corotation at tions for these kinematics, de Vaucouleurs (1964) proposed 2.4 kpc and found a pattern speed of 80 km s−1 kpc−1. the gas response to a bar potential. This explanation has As pointed out by Combes (1996), the discrepancy be- also been discussed in other observational works as those tween Binney et al. estimation and other estimations can of Kerr (1967), Bania (1977) or Liszt & Burton (1978). be due to the fact that the molecular gas distribution in Nevertheless, it is the comparison of the observations with the Galactic center region used by Binney et al. could be numerical simulations of the gas dynamics what have al- the response to a small nuclear bar instead of the large lowed to better understand the structure of the Galaxy bar of the Galaxy. and to derive important dynamical information as the pat- Indeed, there are recent evidences of the presence of tern speed of the bar. The first attempt was that of Mulder a distinct structure in the nuclear region of the Galaxy. & Liem (1986) solving the gas-dynamical equations to get The star counts map made by Alard (2001) with 2MASS the gas flow in a weak bar for the whole Galaxy. data shows an excess of counts in the inner four degrees of A major step in our understanding of the gas kinemat- the Galaxy that resembles a small nuclear bar. This bar ics in the Galactic Center was the paper by Binney et al. could be lopsided with the barycenter displaced to nega- (1991). These authors made a comparison of the CO(1-0) tive longitudes. Nishiyama et al. (2005) have also found and HI longitude-velocity diagrams (hereafter lv-diagram) evidence of a structure different to the large bar or bulge of Bally et al. (1987) and Burton & Liszt (1978), with the in the inner degrees of the Galaxy. lv locus of closed stellar orbits in a bar potential (x1 and In the present paper, we want to investigate the gas re- x2 orbits). The lv-diagram of the molecular gas in the in- sponse to the mass distribution inferred from the 2MASS ner ∼ 4◦ of the Galaxy shows a parallelogram shape that data. In particular, we model the central mass distribution resembles the lv locus of the cusped x1 orbit for a bar with as a small nuclear bar. We are interested in the role of this inclination of ∼ 16◦. In spite of the globally successful in- nuclear bar to explain the differences among the results of terpretation of the CO and HI lv-diagrams, the Binney previous works (Binney et al. 1991; Fux 1999; Combes et al. (1991) model cannot account for the observed asym- 1996). We also want to study whether a lopsided nuclear metry of the CO parallelogram, which is lopsided towards bar could explain the apparent lopsidedness of the molec- positive longitudes while the modeled parallelograms are ular gas distribution in the nuclear region, whose com- centered at longitude 0◦. In addition the expected upper plex dynamics are far from being understood.