Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://lnfm1.sai.msu.ru/~anna/article/aa14646-10.pdf Äàòà èçìåíåíèÿ: Tue Jun 25 16:42:56 2013 Äàòà èíäåêñèðîâàíèÿ: Fri Feb 28 01:06:06 2014 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï ï

A&A 519, A70 (2010) DOI: 10.1051/0004-6361/201014646
c ESO 2010

Astronomy & Astrophysics

N -body simulations in reconstruction of the kinematics of young stars in the Galaxy
P. Rautiainen1 and A. M. Mel'nik2
1 2

Department of Physics/Astronomy Division, University of Oulu, PO Box 3000, 90014 Oulun yliopisto, Finland e-mail: pertti.rautiainen@oulu.fi Sternberg Astronomical Institute, 13, Universitetskii pr., Moscow 119992, Russia

Received 31 March 2010 / Accepted 21 May 2010
ABSTRACT

Aims. We try to determine the Galactic structure by comparing the observed and modeled velocities of OB-associations in the 3 kpc solar neighborhood. Methods. We made N -body simulations with a rotating stellar bar. The galactic disk in our model includes gas and stellar subsystems. The velocities of gas particles averaged over large time intervals (8 bar rotation periods) are compared with the observed velocities of the OB-associations. Results. Our models reproduce the directions of the radial and azimuthal components of the observed residual velocities in the Perseus and Sagittarius regions and in the Local system. The mean difference between the model and observed velocities is V = 3.3 km s-1 . The optimal value of the solar position angle b providing the best agreement between the model and observed velocities is b = 45 ± 5 , in good accordance with several recent estimates. The self-gravitating stellar subsystem forms a bar, an outer ring of subclass R1 , and slower spiral modes. Their combined gravitational perturbation leads to time-dependent morphology in the gas subsystem, which forms outer rings with elements of the R1 - and R2 -morphology. The success of N -body simulations in the Local System is likely due to the gravity of the stellar R1 -ring, which is omitted in models with analytical bars.
Key words. Galaxy: structure ­ Galaxy: kinematics and dynamics

1. Introduction
The consensus since the 1990s has been that the Milky Way is a barred galaxy (see, e.g. Blitz & Spergel 1991; Blitz et al. 1993). The estimate for the size of the large-scale bar has grown from initial Rbar 2-3 kpc to current estimates Rbar = 3-5 kpc (Habing et al. 2006; Cabrera-Lavers et al. 2007, 2008). The position angle of the bar is thought to be in the range 15 -45 (Blitz et al. 1993; Kuijken 1996; Weiner & Sellwood 1999; Benjamin et al. 2005; Englmaier & Gerhard 2006; Cabrera-Lavers et al. 2007; Minchev et al. 2010). The differences in the position angle estimates may indicate that the innermost structure is actually a triaxial bulge (Cabrera-Lavers et al. 2008). On the other hand, this ambiguity may be partly caused by our unfavorable viewing angle near the disk plane, which also hinders study of other aspects of Galactic morphology. The suggested configurations for the spiral morphology of the Galaxy include models or sketches containing from two to six spiral arms (see e.g. VallÈe 2005; VallÈe 2008, and references therein). A case has also been suggested where a two-armed structure dominates in the old stellar population, whereas the gas and young stellar population exhibits a four-armed structure (LÈpine et al. 2001; Churchwell et al. 2009). In addition to spiral arms, there may be an inner ring or pseudoring surrounding the bar, which manifests itself as the so-called 3-kpc arm(s) (Dame & Thaddeus 2008; Churchwell et al. 2009). Also, speculations about a nuclear ring with a major axis of about 1.5 kpc have been made (Rodriguez-Fernandez & Combes 2008). Different kinds of rings ­ nuclear rings, inner rings and outer rings ­ are often seen in the disks of spiral galaxies, especially if there is

also a large-scale bar (Buta & Combes 1996). Thus, the presence of an outer ring in the Galaxy may also be considered plausible (Kalnajs 1991). Since the outer rings have an elliptic form, the broken outer rings (pseudorings) resemble two tightly wound spiral arms. Nevertheless their connection with the density-wave spiral arms is not very obvious because their formation does not need the spiral-shaped perturbation in the stellar disk. The main ingredient for their formation is a rotating bar. Both test particle simulations (Schwarz 1981; Byrd et al. 1994; Bagley et al. 2009) with an analytical bar and N -body simulations (Rautiainen & Salo 1999, 2000), where the bar forms in the disk by instability, show that the outer rings and pseudorings are typically located in the region of the outer Lindblad resonance (OLR). Two main classes of the outer rings and pseudorings have been identified: the R1 -rings and R1 -pseudorings elongated perpendicular to the bar and the R2 -rings and R2 -pseudorings elongated parallel to the bar. In addition, there is a combined morphological type R1 R2 that shows elements of both classes (Buta 1986; Buta & Crocker 1991; Buta 1995; Buta & Combes 1996; Buta et al. 2007). Schwarz (1981) connected two main types of the outer rings with two main families of periodic orbits existing near the OLR of the bar (Contopoulos & Papayannopoulos 1980; Contopoulos & Grosbol 1989). The stability of orbits enables gas clouds to follow them for a long time period. The R1 -rings are supported by x1 (2)-orbits (using the nomenclature of Contopoulos & Grosbol 1989) lying inside the OLR and elongated perpendicular to the bar, while the R2 -rings are supported by x1 (1)-orbits situated a bit outside the OLR and elongated along the bar. There Page 1 of 13

Article published by EDP Sciences


A&A 519, A70 (2010)

is also another conception of the ring formation. Romero-GÑmez et al. (2007) show that Lyapunov periodic orbits around L1 and L2 equilibrium points can lead to the formation of the spiral arms and the outer rings. They associate the spiral arms emanating from the bar's tips with the unstable manifolds of Lyapunov orbits. This approach can be useful for explaining of the motion of gas particles as well (Athanassoula et al. 2009). Besides the bar the galactic disks often contain spiral arms, which modify the shape of the gravitational perturbation. In the simplest case, the pattern speeds of the bar and spiral arms are the same. In many studies this assumption has been used for constructing the gravitational potential from near-IR observations (which represent the old stellar population better than the visual wavelengths). Several galaxies with outer rings have been modeled by this method, and findings are in good accordance with studies made by using analytical bars: the outer rings tend to be located near the OLR (Salo et al. 1999), although in some cases they can be completely confined within the outer 4/1-resonance, (Treuthardt et al. 2008). A real galactic disk provides further complications, which can be studied by N -body models, where the bars and spiral arms are made of self-gravitating particles. In particular, there can often be one or more modes rotating more slowly than the bar (Sellwood & Sparke 1988; Masset & Tagger 1997; Rautiainen & Salo 1999). Even if there is an apparent connection between the ends of the bar and the spiral arms, it is no guarantee that the pattern speeds are equal ­ the break between the components may be seen only for a short time before the connection reappears (see Fig. 2 in Sellwood & Sparke 1988). Sometimes the bar mode can contain a considerable spiral part that forms the observed spiral, together with the slower modes (Rautiainen & Salo 1999). The multiple modes can also introduce cyclic or semi-cyclic variations in the outer spiral morphology: outer rings of different types can appear and disappear temporarily (Rautiainen & Salo 2000). In Mel'nik & Rautiainen (2009, hereafter Paper I), we considered models with analytical bars. In this case the motion of gas particles is determined only by the bar. We found that the resonance between the epicyclic motion and the orbital motion creates systematical noncircular motions that depend on the position angle of a point with respect to the bar elongation and on the class of the outer ring. The resonance kinematics typical of the outer ring of subclass R1 R2 reproduces the observed velocities in the Perseus and Sagittarius regions well. In Paper I we also suggested that the two-component outer ring could be misinterpreted as a four-armed spiral. In some galaxies with the combined R1 R2 -morphology, the R1 -component can also be seen in the near infrared, but the R2 -component is usually prominent only in blue (Byrd et al. 1994). This could explain the ambiguity of the number of spiral arms in the Galaxy. N -body simulations confirm that the R1 -rings can be forming in the self-gravitating stellar subsystem, while the R2 -rings usually exist only in the gas component (Rautiainen & Salo 2000). In the present paper we study the effect of multiple modes and their influence on the kinematics and distribution of gas particles. We construct N -body models to study the influence of self-gravity in the stellar component on the kinematics of gas particles. We compare the model velocities of gas particles with the observed velocities of OB-associations in the neighborhood 3 kpc from the Sun. This paper has the following structure. Observational data are considered in Sect. 2. Section 3 is devoted to models and describes the essential model parameters, the evolution of the
Page 2 of 13

stellar and gas components: formation of the bar and the interplay between the bar and slower spiral modes. In Sect. 3 we also analyze the general features of the gas morphology. Section 4 is devoted to the comparison between the observed and modeled kinematics. Both momentary and average velocities of gas particles are considered. The influence of the bar position angle b on the model velocities is also investigated in Sect. 4, as are the evolutionary aspects of kinematics. Section 5 consists of conclusions and discussion.

2. Observational data
We have compared the mean residual velocities of OB-associations in the regions of intense star formation with those of gas particles in our models. These regions practically coincide with the stellar-gas complexes identified by Efremov & Sitnik (1988). The residual velocities characterize the non-circular motions in the galactic disks. They are calculated as differences between the observed heliocentric velocities (corrected for the motion to the apex) and the velocities due to the circular rotation law. We used the list of OB-associations by Blaha & Humphreys (1989), the line-ofsight velocities (Barbier-Brossat & Figon 2000), and proper motions (Hipparcos 1997; van Leeuwen 2007) to calculate their median velocities along the galactic radius-vector, VR , and in the azimuthal direction, V . Figure 1 shows the residual velocities of OB-associations in the regions of intense star formation. It also indicates the grouping of OB-associations into stellar-gas complexes. For each complex we calculated the mean residual velocities of OB-associations, which are listed in Table 1. Positive radial residual velocities VR are directed away from the Galactic center, and the positive azimuthal residual velocities V are in the sense of Galactic rotation. Table 1 also contains the rms errors of the mean velocities, the mean Galactocentric distances R of OB-associations in the complexes, the corresponding intervals of galactic longitudes l and heliocentric distances r, and names of OB-associations the region includes (see also Mel'nik & Dambis 2009). The Galactic rotation curve derived from an analysis of the kinematics of OB-associations is nearly flat in the 3-kpc solar neighborhood and corresponds to the linear velocity at the solar distance of 0 = 220 km s-1 (Mel'nik et al. 2001; Mel'nik & Dambis 2009). The nearly flat form of the Galactic rotation curve was found in many other studies (Burton & Gordon 1978; Clemems 1985; Brand & Blitz 1993; Pont 1994; Dambis et al. 1995; Russeil 2003; Bobylev et al. 2007). We adopted the Galactocentric distance of the Sun to be R0 = 7.5 kpc (Rastorguev et al. 1994; Dambis et al. 1995; Glushkova et al. 1998, and other papers), which is consistent with the so-called short distance scale for classical Cepheids (Berdnikov et al. 2000).

3. Models
3.1. The model parameters

We made several N -body models, which satisfy "broad observational constraints": the rotation curve is essentially flat and the size of the bar is acceptable. From these models we have chosen our best-fitting case, which we describe here in more detail. The rotation curve of our best-fitting model is illustrated in Fig. 2. In the beginning, the rotation curve is slightly falling in the solar neighborhood, but the mass rearrangement in the disk during the bar formation makes it rise slightly. We scaled


P. Rautiainen and A. M. Mel'nik: N -body simulations of the Galaxy

Fig. 1. a) The residual velocities of OB-associations projected on to the galactic plane. It also shows the grouping of OB-associations into regions of intense star formation. b) The mean VR - and V -velocities of OB-associations in the stellar-gas complexes. The X -axis is directed away from the galactic center, and the Y -axis is in the direction of the galactic rotation. The Sun is at the origin.

the simulation units to correspond to our preferred values of the solar distance from the Galactic center and the local circular velocity. This also gives the scales for masses and time units. However, in the following discussion we will use simulation time units, one corresponding to approximately 100 million years, and the full length of the simulation is 6 Gyr. The bulge and halo components are analytical, whereas the stellar disk is self-gravitating. The bulge is represented by a Plummer sphere, mass Mbulge = 1.17 â 1010 M , and scale length Rbulge = 0.61 kpc. The dark halo was included as a component giving a halo rotation curve of form V ( R) = V
max

are omitted, the velocity dispersion of the test particles rises much higher into the range 25-50 km s-1 . The model used in the kinematical analysis contains 40 000 gas particles initially distributed as a uniform disk with an outer radius of 9.2 kpc.
3.2. Evolution of the stellar component

R
2 c

R2 + R

,

(1)

where Vmax = 210 km s-1 is the asymptotic maximum on the halo contribution to the rotation curve and Rc = 7.6 kpc the core radius. The N -body models are two-dimensional, and the gravitational potential due to self-gravitating particles is calculated by using a logarithmic polar grid (108 radial and 144 azimuthal cells). The N -body code we used has been written by Salo (for more details on the code, see Salo 1991; Salo & Laurikainen 2000). The value of the gravitation softening is about 0.2 kpc on the adopted length scale. The mass of the disk Mdisk = 3.51 â 1010 M . The disk is composed of 8 million gravitating stellar particles, whose initial distribution is an exponential disk reaching about 10 scale lengths. The disk and halo have nearly equal contribution to the rotation curve at the solar distance. The initial scale length of the disk was about 2 kpc, but after the bar formation, it forms a twin profile disk: the inner profile becomes steeper and the outer profile shallower, and the exponential scale length corresponds to about 3 kpc outside the bar region. The initial value of the Toomre-parameter QT was 1.75. The gas disk was modeled by inelastically colliding test particles as was done in Paper I. The initial velocity dispersion of the gas disk was low, about 2 km s-1 , but it reached typical values in the range 5-15 km s-1 during the simulation. If collisions

The inner regions quickly develop a small spiral (at T 2.5), which then evolves to a clear bar (T 5). Its original pattern speed b is about 80 km s-1 kpc-1 , meaning that when it forms it does not have an Inner Lindblad Resonance (ILR). In its early phase the bar slows down quite quickly (b 60 km s-1 kpc-1 at T = 10), but the deceleration rate soon settles down: b 54 km s-1 kpc-1 at T = 20 and b 47 km s-1 kpc-1 at T = 55. In this model the bar's slowing down is accompanied by its growth, and the bar can always be considered dynamically fast (see e.g. Debattista & Sellwood 2000). Using the same method to determine the bar length as Rautiainen et al. (2008) (a modification of one used by Erwin 2005), we get Rbar = 4.0 ± 0.6 kpc at T = 55 and RCR /Rbar = 1.2 ± 0.2. There is no secondary bar in this model. The amplitude spectra of the relative density perturbations (see e.g. Masset & Tagger 1997; Rautiainen & Salo 1999) (Fig. 3) show that the bar mode is not the only one in the disk, but there are also slower modes. The strongest of these modes, hereafter the S1 mode, has an overlap of resonance radii with the bar: the corotation radius of the bar is approximately the same as the inner 4/1-resonance radius of the slower mode (at T = 55 the RCR of the bar and the inner 4/1 resonance radius of the S1 mode are both about 4.6 kpc). This resonance overlap does not seem to be a coincidence: when the amplitude spectra from different time intervals are compared, one can see that both the bar and the S1 modes slow down so that the resonance overlap remains (see Fig. 3). Furthermore, this resonance overlap was the most common case in the simulations of Rautiainen & Salo (1999). Also, the S1 mode has a strong m = 1 signal and a maximum near its corotation at 7.1 kpc. The bar mode is also seen as a strong signal in the m = 4 spectrum, but only inside CR ­ the
Page 3 of 13


A&A 519, A70 (2010)
250 250

200

200

V [km s-1]

100

V [km s-1] 2 4 6 R [kpc] 8 10 12

150

150

100

50 0 0

50 0 0

2

4

6 R [kpc]

8

10

12

Fig. 2. The rotation curve (solid line) of the N -body model at T = 0 (left) and at T = 55 (right). The contributions from the bulge (dash-dotted line), disk (dashed line) and halo (dotted line) are also indicated.
Stars, m=2, T=5-15
60 50 kpc-1] 40 30 20 10 0 0 kpc-1] 60 50 40 30 20 10 0 0

Stars, m=2, T=20-30

-1

p [km s

2

4

6 R [kpc]

8

10

12

p [km s

-1

2

4

6 R [kpc]

8

10

12

Stars, m=2, T=50-60
60 50 p [km s-1 kpc-1] 40 30 20 10 0 0 p [km s-1 kpc-1] 60 50 40 30 20 10 0 0

Stars, m=1, T=50-60

2

4

6 R [kpc]

8

10

12

2

4

6 R [kpc]

8

10

12

Gas, m=2, T=50-60
60 50 kpc-1] 40 30 20 10 0 0 kpc-1] 60 50 40 30 20 10 0 0

Gas, m=1, T=50-60

-1

p [km s

2

4

6 R [kpc]

8

10

12

p [km s

-1

2

4

6 R [kpc]

8

10

12

Fig. 3. The amplitude spectra of the relative density perturbations in the model disk. The frames show the amplitude spectra of the stellar or gas component at various times (indicated on the frame titles). The contour levels are 0.025, 0.05, 0.1, 0.2, 0.4, and 0.8, calculated with respect to the azimuthal average surface density at each radius. The continuous lines show the frequencies and ± /m, and the dashed curves indicate the frequencies ± /4inthe m = 2 amplitude spectrum.

Page 4 of 13


P. Rautiainen and A. M. Mel'nik: N -body simulations of the Galaxy

Bar
10 10

S1

0

0

-10 -10

0

10

-10 -10

0

10

Fig. 4. The reconstructed modes in the stellar component (see text) for T = 50-60 time interval. The enhanced density compared to the azimuthally averaged profile at each radius is shown. The shades of gray (darker corresponds to higher surface density) have been chosen to emphasize the features. The circles in the bar mode indicate ILR (1.4 kpc), CR (4.6 kpc), and OLR (8.1 kpc), whereas the inner 4/1 (4.6 kpc) and CR (7.1 kpc) are shown for the mode S1.

spiral part seems to be almost pure m = 2 mode. Altogether, the signals with m > 2 tend to be much weaker than features seen in m = 1 and m = 2 amplitude spectra. We have also tried to reconstruct the shapes of the modes seen in the amplitude spectra. This was done by averaging the surface density in coordinate frames rotating with the same angular velocities as the modes. No assumptions were made about the shapes of the modes. On the other hand, one should take these reconstructions with some caution, because the evolution of the two modes, the effect of slower (but weaker) modes, and short-lived waves may affect them. The results for the bar and the S1 mode at the time interval T = 50-60 are shown in Fig. 4. The mode p = 47 km s-1 kpc-1 clearly shows the bar and symmetrical spiral structure that forms an R1 outer ring or pseudoring. By the T = 50-60 interval, the density amplitude of the bar mode is about 15-20 per cent in the outer ring region, where the maxima and minima have roughly the same strength. On the other hand, by T = 50-60, the mode p = 31 km s-1 kpc-1 is clearly lopsided, which is not surprising considering the signal seen in the m = 1 amplitude spectrum. There is a minimum with an amplitude of about 30% and a maximum of about 15% at R 7 kpc, which corresponds to the CR of the S1 mode. Earlier, at T 20, the S1 mode does not have the m = 1 characteristic but exhibits a multiple-armed structure beyond its CR, accompanied by a clear signal in the m = 3 amplitude spectrum.
3.3. The morphological changes in the gas component

wound pair of spiral arms. On the broader sense, the overall Hubble stage of the model stays the same for several Gyr. Although the slow modes in the stellar component can be clearly seen outside the bar radius (about 4 kpc), they become pronounced in the gas from R 6 kpc. To study their effect on the gas morphology, we selected gas particles located at the annulus 7 < R < 10 kpc and calculated their number within every 5 -sector along . Such density profiles were built for 301 moments from the interval T = 30-60 (T 3-6 Gyr) with a step T = 0.1 (10 Myr). Earlier stages were not considered, because then the pattern speed of the bar was changing so fast that it complicated the analysis. At every moment the distribution of gas density along was approximated by one-fold (m = 1), twofold (m = 2), and four-fold (m = 4) sinusoidal wave: = 0 + Am cos(m + m ), (2)

The amplitude spectra for the gas component at the interval T = 50-60 are also shown in Fig. 3. Due to fewer particles, they include more noise, but otherwise they are quite similar. In addition to the bar mode, the S1 mode is also seen, but now it is more conspicuous in the m = 1 spectrum. The result of having several modes is the quite complicated evolution of the model (see Fig. 5): at different times, the morphology of the outer gaseous disk can be described as R1 R2 , R2 , R1 or just as open spiral arms, which can sometimes be followed over 400 degrees. There is no evolutionary trend between the morphological stages, since they all appear several times during the model time span. The shape of the inner ring also changes by being sometimes more elongated or even consisting of tightly

where is the gas density in a segment, 0 is the average density in the annulus, m and Am are the phase and amplitude of the corresponding sinusoidal approximation, respectively. Figure 6 demonstrates the motion of maxima in the distribution of gas particles along . We made the density profiles in the reference frame co-rotating with the bar, whose major axis is always oriented in the direction = 0 . Azimuthal angle is increasing in the sense of the galactic rotation, so the supposed position of the Sun is about = 315 . To illustrate the motion of density crests, we selected two intervals T = 35.5-37.5 and T = 52.5-54.5 with a high amplitude of density perturbation. These density profiles indicate the motion of density maxima in the opposite direction to that of galactic rotation (i.e. they actually rotate more slowly than the bar), which means an increase in the phase m of the sinusoidal wave (Eq. (2)). Figure 7 exhibits the variations in the phase m and amplitude Am of the sinusoidal wave at the time intervals T = 30-40, 40-50, and 50-60. The subscripts 1 and 2 are related to the oneand two-fold sinusoids. Rotation of the density maxima causes the sharp changes in the phase when it achieves the value of = 360, and at the new turn its value must fall to zero. These changes enable us to accurately calculate the mean values of the periods for the propagation of the sinusoidal waves, which appear to be P1 = 3.3 ± 0.4and P2 = 1.5 ± 0.4. Remember that we study the density oscillations in the reference frame co-rotating
Page 5 of 13


A&A 519, A70 (2010)

Fig. 5. The gas morphology at selected times. The bar is vertical in all frames, whose width is 20 kpc.

with the bar, so the period P of beating oscillations between the bar and slow modes is determined by the relation: Pm = 2 · m(b - sl ) (3)

The periods, P1 and P2 , appear to correspond to slow modes rotating with the pattern speeds = 28 ± 2 km s-1 kpc-1 and = 26 ± 6 km s-1 kpc-1 , respectively. It is more convenient to use simulation units here. The transformation coefficient between them and (km s-1 kpc-1 ) is k = 9.77, and the value of b is b = 4.8s.u. The m = 4 wave manifested itself in two density maxima separated by the angle 90 (Fig. 6, right panel). The analysis of phase motion of four-fold sinusoid reveals the period P4 = 0.81 ± 0.15, which corresponds to slow mode rotating with the speed sl = 28 ± 4km s-1 kpc-1 (Eq. (3)). Probably, it is mode = 28 ± 4 km s-1 kpc-1 that causes the strong variations in gas density with the periods P1 = 3.3, P2 = 1.5, and P4 = 0.8 when it works as m = 1, m = 2, and m = 4 density
Page 6 of 13

perturbations, respectively. This mode is well-defined in the gas and star power spectra made for the interval T = 50-60 (Fig. 3). Let us have a look at the amplitude variations (Fig. 7). The highest value of A2 equal A2 = 200 (particles per 5 -sector) is observed at the time T = 36.0 (left panel). On the other hand, A1 achieves its highest value of A1 = 220 at the time T = 56.5 (right panel). Amplitude A4 reaches its maximum value of A4 = 180 at the time interval T = 53-55. Thus, the highest values of the amplitudes A1 , A2 , and A4 are nearly the same. Figure 6 (left panel) indicates the growth of the amplitude of m = 2 perturbation under a specific orientation of the density clumps. The amplitude of the sinusoidal wave is at its maximum at the moments T = 36.0 and 37.5 when the density clumps are located near the bar's minor axis, = 90 and 270. This growth is also seen in Fig. 7 (left panel) for the interval T = 30-40: the amplitude A2 is at its maximum at the moments when 2 180 . This phase corresponds to the location of maxima of m = 2 sinusoid at = 90 and = 270 (Eq. (2)).


P. Rautiainen and A. M. Mel'nik: N -body simulations of the Galaxy
'
500 250 0 -250 -500 0 T=36.0 90 180 270 360 T=35.5

'
500 250 0 -250 -500 0 T=53.0 90 180 270 360 T=52.5

500 250 0 -250 -500

500 250 0 -250 -500

0 T=36.5

90

180

270

360

0 T=53.5

90

180

270

360

500 250 0 -250 -500

500 250 0 -250 -500

0 T=37.0

90

180

270

360

0 T=54.0

90

180

270

360

500 250 0 -250 -500

500 250 0 -250 -500

0 T=37.5

90

180

270

360

0 T=54.5

90

180

270

360

500 250 0 -250 -500

500 250 0 -250 -500

0

90

180

270

360

0

90

180

270

360





Fig. 6. The perturbation in the density of gas particles, = - 0 , located at the annulus 7 < R < 10 kpc along azimuthal angle built for different moments. It also shows its approximation by two-fold (left panel) and one-fold (right panel) sinusoids.

2

T=30 - 40

A2
180 120 60



1

T=30 - 40

A

1

270 180 90 30 32 T=40 - 50 34 36 38

270 180 90 30 32 T=40 - 50 34 36 38

180 120 60 40

40



2

A2
180 120 60



1

A

1

270 180 90 40 42 44 46 48

270 180 90 40 42 44 46 48

180 120 60 50

50



2

T=50 - 60

A2
180 120 60



1

T=50 - 60

A

1

270 180 90 50 52 54 56 58

270 180 90 50 52 54 56 58

180 120 60 60

60

T

T

Fig. 7. Variations in the phase (black curve) and amplitude A (gray curve) of the sinusoids that approximate the distribution of gas particles located at the distances 7 < R < 10 kpc along . Subscripts 1 and 2 are related to the one- and two-fold sinusoids, respectively.

Page 7 of 13


A&A 519, A70 (2010)
Momentary velocities VR
30 Sag 20 10 0 -10 -20 -30 50 30 Car 20 10 0 -10 -20 -30 50 30 Cyg 20 10 0 -10 -20 -30 50 30 LS 20 10 0 -10 -20 -30 50 30 Per 20 10 0 -10 -20 -30 50 52 54 56 58 60 52 54 56 58 60 52 54 56 58 60 52 54 56 58 60 52 54 56 58 60

Momentary velocities V
30 Sag 20 10 0 -10 -20 -30 50 30 Car 20 10 0 -10 -20 -30 50 30 Cyg 20 10 0 -10 -20 -30 50 30 LS 20 10 0 -10 -20 -30 50 30 Per 20 10 0 -10 -20 -30 50 52 54 56 58 60 52 54 56 58 60 52 54 56 58 60 52 54 56 58 60 52 54 56 58 60

Time

Time

Fig. 8. Variations in the mean velocities of gas particles located within the boundaries of the stellar-gas complexes. The left panel is related to the radial component VR and the right one to the azimuthal one V .

Our analysis revealed slight variations in the speed of the strongest slow mode, and they depend on its orientation with respect to the bar: Fig. 7 (left panel) shows that the tilt of the phase curve, 2 (t), is variable. We can see that the slow mode rotates a bit faster when 2 180 (density clumps are near the bar's minor axis) and more slowly when 2 = 0 or 360 (the clumps are near the bar's major axis). Probably, the variations in the speed of the slow mode are connected with the change in the form of the density crests due to tidal interaction between the bar mode (bar+R1 -ring) and the slow mode.

4. Kinematics of gas particles. Comparison with observations
4.1. Momentary and average velocities

We start our kinematical study with the interval T = 50-60 (5-6 Gyr in physical ti