Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://lnfm1.sai.msu.ru/~anna/article/AL09p609.pdf Äàòà èçìåíåíèÿ: Mon Aug 31 13:20:44 2009 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 22:54:29 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: m 63

ISSN 1063-7737, Astronomy Letters, 2009, Vol. 35, No. 9, pp. 609­624. c Pleiades Publishing, Inc., 2009. Published in Russian in Pis'ma v Astronomicheski Zhurnal, 2009, Vol. 35, No. 9, pp. 676­692. i

Kinematics of the Outer Pseudorings and the Spiral Structure of the Galaxy*
A. M. Mel'nik
1

1**

and P. Rautiainen2

Sternberg Astronomical Institute, Universitetskii pr., 13, Moscow, 119992 Russia 2 Astronomy Division, Department of Physical Sciences, University of Oulu, P.O. Box 3000, FIN-90014 Oulun yliopisto, Finland
Received November 13, 2008

Abstract--The kinematics of the outer rings and pseudorings is determined by two processes: the resonance tuning and the gas outflow. The resonance kinematics is clearly observed in the pure rings while the kinematics of the outflow is manifested itself in the pseudorings. The direction of systematical motions in the pure rings depends on the position angle of a point with respect to the bar major axis and on the class of the outer ring. The direction of the radial and azimuthal components of the residual velocities of young stars in the Perseus, Carina, and Sagittarius regions can be explained by the presence of the outer pseudoring of class R1 R2 in the Galaxy. We present models which reproduce the values and directions of the residual velocities of OB-associations in the Perseus and Sagittarius regions and also model reproducing the directions of the residual velocities in the Perseus, Sagittarius, and Carina regions. The kinematics of the Sagittarius region accurately defines the solar position angle with respect to the bar elongation, b = 4 5 ± 5 . PACS numbers : 98.10.+z; 98.20.-d; 98.20.Af; 98.35.-a; 98.35.Ac; 98.35.Hj; 98.35.Jk; 98.35.Pr DOI: 10.1134/S1063773709090047 Key words: Galaxy (Milky Way), spiral structure, kinematics and dynamics, resonances.

1. INTRODUCTION

1.1. The Galactic Spiral Structure
The spiral structure is clearly observed in many external galaxies viewed face-on, but in the Galaxy we are faced with a difficulty of the distance determination for the indicators of spiral structure. On the other hand, in the Galaxy we can study the field of space velocities using both line-of-sight velocities and proper motions, whereas with other galaxies we are almost completely limited to line-of-sight velocities. One of the best tracers of the spiral structure are H II regions--gas clouds which are ionized by young hot stars. These can be seen as bright radio objects throughout the disc of the Galaxy (Y.M. Georgelin and Y.P. Georgelin 1976; Russeil 2003; and other papers). Though the distances for the distant H II regions, r > 4 kpc, are determined from observations badly (usually they are derived from the kinematical models) the distribution of H II regions
**

displays the most general features of the Galactic spiral structure which can be formulated as follows. 1. The Galactic spiral structure is most pronounced at the Galactocentric distances R = 5­9 kpc. 2. The pitch angle of the Galactic spirals is quite small, i < 20 . 3. Many researchers believe the Galaxy is 4-armed. In the vicinity of 3 kpc from the Sun the optical data has revealed the existence of three fragments of the spiral structure: Sagittarius-Carina, Cygnus-Orion, and Perseus ones. The characteristic of these regions is the intense star formation manifested as increased concentration of young clusters and OB-associations (Humphreys 1979; and other papers). The Sagittarius-Carina and Perseus arm-fragments are often thought to be part of the global spiral structure (Y.P. Georgelin and Y.M. Georgelin ´ 1976; Efremov 1998; Russeil 2003; Vallee 2005; and other papers). The Cygnus-Orion fragment is usually regarded as the local arm. This conception appears mostly because of its location between two global arms.
609

The article was translated by the authors. E-mail: anna@sai.msu.ru


610

MEL'NIK, RAUTIAINEN

()

X

Y

1 kpc 30 km s­1

(b)

X

Perseus Local System

Y
Carina Cygnus

Sagittarius

The density-wave theory which had been already developed by the end of the 60s (Lin et al. 1969; Lin 1970; Roberts 1969) afforded an opportunity to determine the location of spiral arms with respect to the corotation radius (CR) through the direction of gas streaming motions in spiral arms. The direction of streaming motions in the Perseus region strongly suggests that its location is inside the CR. The density-wave theory describes the kinematics of young stars and gas in the Perseus, Cygnus, and Carina regions quite well (Mel'nik 2003). Nevertheless its application to the whole 3-kpc solar neighborhood encounters some difficulties. Young stars in the Carina and Sagittarius regions, through which the spiral arm is traditionally drawn, have different systematical motions and cannot be fragments of the same density-wave spiral arm. Moreover, the young stars with the systematical non-circular motions directed towards the Galactic center (the Perseus, Cygnus, and Carina regions) which could belong to the same, perhaps patchy, spiral arm fall nicely on a large-scale structure which appears to be the leading spiral arm (Mel'nik 2006). In the present paper we intend to build the dynamical models that reproduce the kinematics of young stars in the 3-kpc solar neighborhood. We'll show that model of the Galaxy with an outer ring of class R1 R2 can explain the kinematics of young stars in some regions.

1 kpc
Galactic center

30 km s­1

Fig. 1. (a) The residual velocities of OB-associations projected onto the galactic plane. The union of OB-associations into regions of intense star formation. (b) The average velocities VR and V of OB-associations in the regions of intense star formation. The X axis is directed away from the galactic center, the Y axis is in the direction of the Galactic rotation. The Sun is at the origin.

The study of the Galactic structure in the neutral hydrogen revealed the existence of regions with systematical non-circular motions which were later identified with spiral arms (Burton, 1971). The investigation of the Perseus region in OB-associations, red supergiants, molecular and neutral hydrogen showed the presence of systematical motions which were interpreted as motions directed towards the Galactic center (Burton and Bania 1974; Humphreys 1976; Avedisova and Palous 1989; Brand and Blitz 1993; and other papers). Recent studies agree with this result (Mel'nik et al. 2001; Sitnik 2003; and other papers).

1.2. The Kinematics of Young Stars within 3 kpc of the Sun OB-associations are the most suitable objects for the kinematical investigations in the wide solar neighborhood. These loose groups of high-luminosity stars have quite reliable distances (the average accuracy is about 15%). This good accuracy appears due to young clusters which often enter OB-associations (Garmany and Stencel 1992). But unlike young clusters OB-associations contain a sufficient number of stars with known line-of-sight velocities and proper motions. On average the space velocities of OB-associations are determined over 11 stars. The stellar proper motions were taken from the Hipparcos catalog (1997). The electronic version of the catalog of line-of-sight velocities and proper motions of OB-associations is available at http://lnfm1.sai.msu.ru/anna/page3.html (for more details, see Mel'nik et al. 2001). Figure 1a shows the residual velocities of OB-associations in a projection onto the Galactic plane. The residual velocities were calculated as the differences between the observed heliocentric velocities and the velocities due to the circular rotation law and the solar motion relative to the centroid
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS Table 1. Average residual velocities of OB-associations Region Sagittarius Carina Cygnus Local System Perseus R,kpc VR ,kms 5.6 6.5 6.9 7.4 8.4
-1

611

V ,km s

-1

l,deg 8­23

r,kpc

Associations

+11 ± 3 -6 ± 2 -7 ± 3 +5 ± 3 -7 ± 2

-1 ± 1 +5 ± 3 -11 ± 2 +0 ± 3 -5 ± 2

1.3­1.9 Sgr OB1, OB7, OB4, Ser OB1, OB2, Sct OB2, OB3

286­315 1.5­2.1 Car OB1, OB2, Cru OB1, Cen OB1, Coll 228, Tr 16, Hogg 16, NGC 3766, 5606 73­78 1.0­1.8 Cyg OB1, OB3, OB8, OB9 0­360 0.1­0.6 Per OB2, Mon OB1, Ori OB1, Vela OB2, Coll 121, 140, Sco OB2 104­135 1.8­2.8 Per OB1, NGC 457, Cas OB8, OB7, OB6, OB5, OB4, OB2, OB1, Cep OB1

of OB-associations. The parameters of the circular Galactic rotation law and the components of the solar motion were derived from the sample of OB-associations located within 3 kpc from the Sun (Mel'nik et al. 2001). In such a large region systematical non-circular motions of OB-associations can be regarded as random deviations from the rotation curve. The obtained parameters describe the rotation of the low-dispersed Galactic subsystem. Calculation of the residual velocities with respect to the rotation curve derived from the same objects yields minimal residual velocities, any other rotation curve would produce larger, on average, deviations. Besides, the use of the self-consistent distance scale also decreases the deviations from the rotation curve (Sitnik and Mel'nik 1996; Dambis et al. 2001). The derived rotation curve is nearly flat and corresponds to the linear velocity at the solar distance of 0 = 215 km s-1 , the value of remains constant inside 3 kpc from the Sun with the accuracy ±3%. We adopted the Galactocentric distance of the Sun to be R0 = 7.1 kpc (Dambis et al. 1995; Glushkova et al. 1998; and other papers). The distances for OB-associations from the catalog by Blaha and Humphreys (1989) were multiplied by a factor of 0.8 to be consistent with the so-called short distance scale for classical Cepheids (Berdnikov et al. 2000). Figure 1a also shows the union of OB-associations into regions of intense star formation which practically coincide with the stellar-gas complexes united by Efremov and Sitnik (1988). Figure 1b represents the average residual velocities of OB-associations in regions of intense star formation in projection onto the radial, VR , and azimuthal, V , directions. It is a generalization of Fig. 1a. The average residual velocities and their random errors are given in Table 1, which also contains the average Galactocentric distance R, the interval of galactic longitudes l, the range of heliocentric distances r , and names of OB-association the region includes.
ASTRONOMY LETTERS Vol. 35 No. 9 2009

The study of the Hipparcos catalog (1997) shows that systematical errors in positions of bright stars (mV < 9m ) don't exceed 0.0001 arcsec (Kovalevsky 2002). Stars of the Blaha and Humphreys' catalog (1989) having the known proper motion are bright . enough, their average magnitude equals mV = 7m3 . in the region 0 < r < 3 kpc and mV = 7m8 in the region 1.5 < r < 3 kpc. Since the mission of Hipparcos (1997) continued several years (37 months), we can suppose that systematical errors of proper motions of stars of OB-associations don't exceed 0.0001 arcsec yr-1 , that corresponds to the velocity of 1 km s-1 at the distance r = 2 kpc. Besides, the contribution of systematical errors decreases after averaging of proper motions over a large area: most of OB-associations occupy on the sky more than 10 square degrees. The velocity VR in the Perseus, Cygnus and Carina regions is directed towards the Galactic center and is about VR = -7 km s-1 . In the Sagittarius region it is directed away from the Galactic center, VR = +11 km s-1 . This means that Sagittarius and Carina regions cannot belong to the same trailing density-wave spiral arm: the Sagittarius region that is closer to the Galactic center has kinematics typical for a spiral arm outside the CR, whereas the Carina region has kinematics typical for its being inside the CR. Thus, the simple trailing density-wave spiral arm is not applicable to the Galaxy. The velocity V is conspicuous only in the Carina, Cygnus, and Perseus regions. In the Carina region the velocity V is directed in the sense of Galactic rotation (V = +5 km s-1 ) while in the Perseus (V = -5 km s-1 ) and Cygnus (V = -11 km s-1 ) regions it is in the opposite sense.


612

MEL'NIK, RAUTIAINEN

1.3. Morphology and Modelling of the Outer Pseudorings
The essential characteristic of the galaxies with the outer rings and pseudorings--incomplete rings made up of spiral arms--is the presence of the bar (Buta 1995; Buta and Combes 1996). Two main classes of the outer rings and pseudorings have been identified: the R1 rings (R1 pseudorings) elongated perpendicular to the bar and the R2 rings (R2 pseudorings) elongated parallel to the bar. In addition, there is a combined morphological type R1 R2 which shows elements of both classes. The R2 rings have elliptical shape, but the R1 rings are often "dimpled" near the bar ends. There is also a lot of outer rings/pseudorings that cannot be classified into previous classes because their morphological characteristics are unclear or the inclination prevents detailed classification (Buta 1995; Buta and Crocker 1991; Buta et al. 2007). The outer rings are typically observed in earlytype galaxies. For the galaxies in the lower red-shift range the frequency of the outer rings is found to be about 10% of all types of spiral galaxies. But for the early-type sample it increases to 20% (Buta and Combes 1996). A study by Buta (1995) shows the following distribution among outer ring types: 18% (R1 ), 37% (R1 ), under 1% (R2 ), 35% (R2 ), and 9% (R1 R2 ). The small fraction of the complete R2 rings may be due to the selection effects--they lack conspicuous features, for example, "dimples," and so the definite classification can be difficult due to orientation uncertainties. The test particle simulations (Schwarz 1981; Byrd et al. 1994; Rautiainen and Salo 1999) and N-body simulations (Rautiainen and Salo 2000) show that the outer pseudorings are typically located in the region of the Outer Lindblad Resonance (OLR) and are connected with two main families of periodic orbits. Main families of stable periodic orbits are followed by most non-periodic orbits after the introduction of the bar potential (Contopoulos and Papayannopulos 1980). The R1 rings are supported by periodic orbits lying inside the OLR and elongated perpendicular to the bar, while the R2 rings are supported by orbits situated outside the OLR and elongated along the bar. As convincingly shown by Schwarz (1981) the pseudorings appear before the pure rings. According to Schwarz (1981, 1984), the type of an outer ring is determined by two factors: the bar's strength and the initial distribution of gas particles. Under the stronger bar forcing the pure R1 ring no longer forms, instead, particles move outside through the OLR and form the pseudoring R2 . However, this idea contradicts observations: the R2 pseudorings are more frequently

observed in SAB galaxies than in SB ones (Buta 1995). Simulations demonstrate that such factors as the strong bar, a large radius of the initial particle distribution, and a large model time favor the formation of the R2 component. Byrd et al. (1994) find that the R1 component appears quickly and the R2 component forms slower. All their models yield the R1 pseudoring at the earlier time and the R2 or the combined R1 R2 pseudoring at the later time. In some galaxies with the combined morphology the R1 component can be seen in infrared, but the R2 component is usually prominent only in blue. Byrd et al. (1994) explain this fact by the age difference between two components. N -body simulations suggest another explanation. Rautiainen and Salo (2000) think that at least some stellar R1 rings are not remnants of previous ring-shaped star formation episode but forms due to self-gravity in the stellar subsystem. Galactic disks in some simulations demonstrate the presence of the slow modes near the radius of the outer ring (Rautiainen and Salo 1999, 2000). We like to emphasize that the present paper explains the kinematics of young stars without invoking the slow modes.

1.4. The Bar, Rotation Curve, and a Prototype of the Galaxy
Our Galaxy certainly has the bar. The gas kinematics in the central region, infrared photometry, star counts, and other modern tests confirm this fact (Weiner and Sellwood 1999; Benjamin et al. 2005; Englmaier and Gerhard 2006; Habing et al. 2006; Cabrera-Lavers et al. 2007; and other papers). There is ample evidence suggesting that the major axis of the bar is oriented in the direction b = 15 ­45 in such a way that the end of the bar closest to the Sun lies in the first quadrant. However, the angular speed of the bar b and its length are determined from observations badly. Some researchers believe that the CR of the bar lies at the distance range R = 3­4 kpc (Englmaier and Gerhard 2006; Habing et al. 2006; and references therein), whereas others suggest that the Galaxy has a longer bar with the major axis of a = 4­5 kpc (Weiner and Sellwood 1999; Benjamin et al., 2005; Cabrera-Lavers et al. 2007; and references therein). Analysis of orbits in barred galaxies shows that a bar cannot reach beyond its CR--outside this radius the main family of periodic orbits becomes oriented perpendicular to the bar, thus unable to support it (Contopoulos and Papayannopulos 1980). This sets an upper limit for the angular speed of the bar. However, the lower limit is less clear. In general the CR
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS

613

300 250

(a) Model 1

80 60

(b) Model 1

200 150 100 20 50 0 300 250 60 , km/s/kps , km/s 200 150 100 50 0 250 200 150 40 100 50 0 2 4 6 8 10 12 R, kpc
Fig. 2. Left panel: the rotation curves in models 1­3. The continuous lines show the total curve (km s-1 ) while the dashed, dotted and dash-dotted lines show the contribution of the bulge, bar, and halo, respectively. Right panel: the curves of the angular speeds. The continuous lines indicate the angular speeds and ± /2 while the dashed lines-- ± /4 (km s-1 kpc-1 ). The straight dotted line represents the angular speed of the bar bar in each case. The resonance distances are determined by its intersections with the curves of angular speeds. Distances are given in kpc.

40

0 Model 2 80 Model 2

40 20 0

Model 3

80 60

Model 3

20

0

2

4

6

8

10

12

is believed to lie within 1.0­1.4 Rbar , although higher values for RCR /Rbar has also been suggested for some galaxies (for example, Rautiainen et al. 2005). For our study the more important thing is the location of the OLR which is supposed to lie in the solar neighborhood (Kalnajs 1991; and other papers). For a flat rotation curve, the location of the OLR between the Sagittarius (R = 5.6 kpc) and Perseus (R = 8.4 kpc) regions corresponds to b lying within the limits b = 43­65 km s-1 kpc-1 and the CR within the range RCR = 3.3­5.0 kpc.
ASTRONOMY LETTERS Vol. 35 No. 9 2009

Observations don't give an unambiguous answer on the question about the form of the Galactic rotation curve in the central region, R < 3 kpc. The lineof-sight velocities at the tangential points indicate a peak (R 0.3 kpc) and a local minimum (R 3 kpc) on the rotation curve, though the depth of the minimum is less than 50 km s-1 (Burton and Gordon 1978). However, the apparent peak and minimum can be due to the perturbation of the circular velocities by the bar (Englmaier and Gerhard 2006; and other papers). In the outer region the Galactic rotation curve


614 Table 2. The potential components Model Bulge Mass [10 rb [kpc] Halo rc [kpc] vmax [km s Bar Mass [10 a, b [kpc]
10 -1 10

MEL'NIK, RAUTIAINEN

1 M] 1.22 0.31 8.0 ] 251.6 1.82 3.82, 1.20

2 1.21 0.33 4.15 209.5 0.99 3.98, 1.25

3 0.72 0.34 4.20 217.3 1.08 4.03, 1.26

included a disc component, because a flat rotation curve can be achieved without it. The bulge potential is a Plummer sphere (see, for example, Binney and Tremaine 2008), which defines the slope of the rotation curve in the inner region. The outer part is dominated by a halo whose rotation curve is 2 2 v 2 (r ) = vmax r 2 /(r 2 + rc ), (1) where vmax is the asymptotic maximum on the halo rotation curve and rc is a core radius. This should not be considered as a pure halo component, because, strictly speaking, we do not make specific assumptions of the halo-disc mass ratio. We are making two-dimensional simulations without self-gravity: the particles are feeling the axisymmetric potential (essentially the rotation curve and its slope, see Binney and Tremaine 2008) and the non-axisymmetric bar perturbation. The disc in this case can be thought to consist both of the bar and part of the halo. This approach is common in studies of ring formation or orbits in barred galaxies. The gravitational effects of the rings and spiral arms are omitted in this stage of study. The bar is modelled as a Ferrers ellipsoid, whose density distribution is 0 (1 - m2 )n , 0, m > 1, m1 (2)

M]

is nearly flat. The line-of-sight velocities for H IIregions, molecular clouds, and Cepheids suggest this idea (Brand and Blitz 1993; Dambis et al. 1995; Russeil 2003; and other papers). Our hypothesis is that the Galactic rotation curve is nearly flat at the distance range R = 1­12 kpc, though it can have a small peak and a local minimum in the central region. As will be shown below, the kinematics of young stars in the Perseus region indicates the existence of the R2 ring, while the velocities in the Sagittarius region suggest the presence of the R1 ring in the Galaxy. We suggest that the Galaxy has the combined morphology R1 R2 . Buta and Crocker (1991) regard the galaxies ESO 509-98 and ESO 507-16 as typical examples of the R1 R2 morphology. Here are some other examples of galaxies with the R1 R2 morphology which can be also supposed as possible prototypes of the Galaxy: ESO 245-1, NGC 1079, NGC 1211, NGC 3081, NGC 5101, NGC 5701, NGC 6782, and NGC 7098. 2. MODELS The simulation program we use was developed by Salo (1991). It has been used in both self-gravitating simulations (Rautiainen and Salo 1999) and models where the gravitational potential has been derived from infrared observations (for example, Rautiainen et al. 2005). Here we use analytical expressions for the potential components, as was also done by Byrd et al. (1994) in simulations with an earlier version of the same code. We constructed a series of models which have essentially flat rotation curve in the region corresponding to the solar neighborhood. The galactic potential in all models consists of three components: a bulge, halo and a bar (Fig. 2). We have not additionally

where m equals m2 = x2 /a2 + y 2 /b2 and in our case n = 1. The Ferrers ellipsoid is often used in gasdynamical simulations and in orbital analysis (for example, Athanassoula 1992; Romero-Gomez et al. 2007). The parameters of the previously mentioned gravitation potential components were scaled so that the linear velocity of the rotation curve at the solar distance equals observational one (215 km s-1 ). We achieved the best fit between the model and observed velocities using small variations in the scale and in the solar position angle. In particular, we suppose that the Sagittarius region (5.6 kpc) is lying near the point D of the R1 ring where velocities V equal zero (see Section 3.2). In simulation units the values of the bar's semi-axes, a and b, are the same in all models, but the fine-tuning introduces small differences in the length scales (Table 2). One time unit in our simulations corresponds to the physical time of T = 62­75 Myr. The angular speed of the bar has the value bar = 52­55 km s-1 kpc-1 which corresponds to the bar rotation period T = 112­118 Myr. In our models, the non-axisymmetric perturbation is gradually turned on during four bar rotation periods. This was done to avoid possible transient effects. The initial stage assumes only small deviations from the circular motion. One should especially note that
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS

615

10 5 0 ­5 ­10 10 5 0 ­5 ­10 10 5 0 ­5 ­10 10 5 0 ­5 ­10 ­10

Model 1, T = 7.5

Model 2, T = 7.5

Model 3, T = 7.5

Model 1, T = 11.25

Model 2, T = 11.25

Model 3, T = 11.25

Model 1, T = 15

Model 2, T = 15

Model 3, T = 15

Model 1, T = 25

Model 2, T = 25

Model 3, T = 25

­5

0

5

10 ­10

­5

0

5

10 ­10

­5

0

5

10

Fig. 3. The distribution of OB-particles in models 1­3at T = 7.5, 11.25, 15.0, and 25.0. ASTRONOMY LETTERS Vol. 35 No. 9 2009


616

MEL'NIK, RAUTIAINEN

Y VR

(a) V

Y

(b)

X

X

Fig. 4. The distribution of particles (gas + OB) with the negative and positive residual velocities in model 1 at the moment T = 15. The bar is oriented along the Y axis, the galaxy rotates clockwise, a division on the X and Y axis corresponds to 1 kpc. (a) Represents the radial velocities, while (b) the azimuthal ones. Only 20% of particles are shown. Particles with the positive residual velocities (VR > 5 km s-1 or V > 5 km s-1 ) are indicated by black circles, while particles with the negative ones (VR < -5 km s-1 or V < -5 km s-1 )-- by gray circles. Particles with the velocities close to zero (-5 < VR < +5 km s-1 or -5 < V < +5 km s-1 ) are shown by points.

the bar mass is included in the models from the beginning: initially only its m = 0 component has a nonzero value. The gas subsystem is modelled by 50 000 massless test particles (the initial surface density is uniform in the occupied region) that can collide with each other inelastically. Initially, the gas disk is cold, its velocity dispersion corresponds to only few km s-1 . The velocity dispersion rises during simulation, but does not exceed the typical values seen in the galaxies. To detect regions corresponding to OB-associations, we used a simple recipe of the star formation: in each collision there is a small probability (typically 0.1) that the gas cloud becomes a massless "OB-association." These associations have a limited lifetime before returning as gas particles. With the adopted units of lengths and velocities it corresponds to about four million years. Our analysis includes both gas particles and OB-particles. Though their behavior is quite similar, there are also some differences: the collision-based star formation recipe favors regions where orbits are crossing. Thus, the density of OB-particles is not directly proportional to the density of gas clouds. Figure 3 shows the distribution of OB-particles in models 1­3 at the different time steps. The evolution of models can be briefly described as follows. When the bar is turned on a nuclear ring and an inner ring are quickly forming near the inner Lindblad resonance and the inner 4/1-resonance, respectively (see, for example, Buta and Combes, 1996). These rings also quickly disappear, only model 3 keeps a well-defined

inner ring by the moment T = 7.5. As for the outer region, all models form an outer pseudoring R2 by the time T = 7.5 and nearly pure rings of class R1 R2 by the moment T = 15. In its early evolution the R2 pseudoring is not oriented strictly parallel to the bar, its major axis is slightly leading with respect to the bar (T = 11.5, models 1 and 3), but by the time T = 15 it adjusts to parallel or almost parallel alignment with the bar. In a long simulation, the R1 and R2 rings acquire a more circular form (T = 25). In the case of the strong bar (model 1) the evolution is faster in general. The evolution in the outer region depends also on the initial gas particle distribution. If the gas disc is very extended, an outer pseudoring of class R2 can form almost as fast as the R1 component, but if its extent is small (Rcut < ROLR ), the R2 component forms more slowly. A very small radius of the initial particle distribution can delay the appearance of the outer pseudorings and even inhibit the formation of an R2 pseudoring. 3. KINEMATICS OF THE OUTER RINGS AND PSEUDORINGS 3.1. The Resonance Kinematics The modelling suggests that the kinematics of the outer rings and pseudorings is determined by two processes: the resonance tuning and the gas outflow. The resonance kinematics is clearly observed in pure rings while the kinematics of the gas outflow is manifested itself in the pseudorings.
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS

617

30 20 10 0 ­10 ­20 ­30 15 10 5 0 ­5 ­10 ­15

V, km s­1

VR, km s­1

­ 100

0 , deg

100

Fig. 5. The oscillations of the velocities, VR and V , of particles (gas + OB) located in the narrow annulus R = 7.38 ± 0.16 kpc in model 2 at T = 15.

Figure 4 exhibits the distribution of gas particles and OB-particles with the positive and negative residual velocities in projection onto the radial and azimuthal directions at the time moment T = 15 (8 bar rotation periods). Model 1 was chosen for illustration, other models display a very similar velocity distribution. The residual velocities were calculated as differences between the model velocities and the velocities due to the rotation curve. The alternation of the quadrants with the negative and positive residual velocities is clearly seen in the outer rings. Besides, particles located at the same azimuthal angle but forming either the R2 or the R1 ring have the opposite residual velocities. The oscillation of the azimuthal velocity, V , is shifted by the angle /4 with respect to the radial-velocity oscillation. To study oscillations of the residual velocities we selected particles (gas + OB) located in the narrow annulus of R = 7.38 ± 0.16 kpc in different models at different moments. Figure 5 shows the velocity oscillations of selected particles in model 2 at T = 15. The position of the selected particles corresponds to the average radius of the R2 ring. The velocity profiles made at the same radius R but for different models and moments resemble each other, the main difference is seen in velocity amplitudes. Table 3 represents the amplitudes of velocity oscillations, fR and f , calculated for selected particles in different models at different moments. The highest velocity amplitudes are observed in model 1 while the amplitudes in models 2 and 3 are lower. In each model the lowest velocity amplitudes are observed at T = 25. We combined the samples of gas particles and OB-particles because their kinematics is very similar. That is not surprising because OB-associations in
ASTRONOMY LETTERS Vol. 35 No. 9 2009

our models have the velocities of their parent gas clouds.

3.2. The Orbital Kinematics in the Resonance Region The resonance between the epicyclic and orbital motions adjusts the epicyclic motions of gas clouds in accordance with orbital rotation. That creates systematical non-circular motions of gas clouds whose direction depends on the position angle of a point with respect to the bar major axis and on the class of the outer ring. Figure 6 shows the directions of the residual velocities at the different points of the outer pseudorings. It schematically represents the bar and two main families of periodic orbits. The galaxy rotates clockwise. Motions are considered in the reference frame rotating with the angular speed of the bar. In this frame stars located near the OLR rotate with the speed ( - b ) counterclockwise.
Table 3. Amplitudes of the velocity oscillations at the distance R = 7.38 kpc Model T = 7.5 fR , f ,km s T = 15 fR , f ,km s T = 25 fR , f ,km s
-1 -1 -1

1 24, 16 31, 15 25, 12

2 22, 12 23, 12 12, 6

3 24, 13 25, 12 15, 8


618

MEL'NIK, RAUTIAINEN



3 4 D C B

2

5E F 6 G

A1

H 8

7

Fig. 6. Orbital kinematics in the resonance region. It shows the bar and two main families of periodic orbits in the region of the OLR. The galaxy rotates clockwise. Motions are considered in the reference frame rotating with the speed of the bar. The orbit going through the points 1-2-3-4-5-6-7-8 represents the family of periodic orbits elongated along the bar. The orbit A-B-C-D-E-F-G-H denotes the family of periodic orbits elongated perpendicular to the bar. The epicycles drawing on the border of the picture show the position of the particle on the epicyclic orbit at the points 1­8. The vectors show the directions of the additional (residual) velocities due to the epicyclic motion. The dashed lines indicate the orbital segments with the negative radial velocities, while the solid lines--those with the positive velocities VR . The likely position of the Sun is shown by the specificsymbol.

Let us consider the family of periodic orbits elongated along the bar and located outside the OLR. It is represented by the orbit passing through the points 1­2­3­..­8 (Fig. 6). This family is the backbone of the R2 ring. The points 1 and 5 are the peri-centers of the orbit, but the points 3 and 7 are the apocenters. The epicycles drawing along the border of the figure demonstrate the position of a particle on the epicyclic orbit at the points 1­8. The arrows show the additional velocities caused by the epicyclic motion which can be regarded as the residual velocities. Let us consider the projection of the residual velocities onto the radial and azimuthal directions. At the points 2 and 6 the velocity VR is directed away from the galactic center and achieves its extremal positive value, while at the points 4 and 8 it is directed toward the galactic center and exhibits its extremal negative value. Stars have the negative radial velocities on the orbital segments 3­5 and 7­1 which are marked by dashed lines. Azimuthal residual velocity, V , equals zero at the points 2, 4, 6, and 8; it achieves its extremal

negative value at the points 3 and 7 and its extremal positive amount at the points 1 and 5. Another family of periodic orbits oriented perpendicular to the bar is represented by the orbit passing through the points A­B­C­..­H (Fig. 6). The "dimples" near the points C and G appear due to a relatively large size of the epicycle in comparison with the orbital radius. This family is the backbone of the R1 ring. The radial velocity, VR , achieves its extremal positive value at the points D and H and the extremal negative one at the points B and F. Points B, D, F, and H are located not at the middles of the corresponding segments but a bit closer to the bar's ends. The azimuthal velocity, V , is zero at the points B, D, F, and H; it achieves its extremal positive value at the points C and G and the extremal negative one at the points A and E. The directions of the residual velocities in models 1­3 at the moment T = 15 (Fig. 4) are in a good agreement with the resonance kinematics shown in Fig. 6. Probably, the kinematics of the pure rings
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS

619

Y VR

(a) V

Y

(b)

X

X

Fig. 7. The distribution of gas + OB particles with the negative and positive residual velocities in model 1 at T = 7.5 radial velocities (a), azimuthal velocities (b). Particles with the positive residual velocities are indicated by black circles, while particles with the negative ones--by gray circles. Particles with the residual velocities close to zero are shown by points. For more details see the caption to Fig. 4.

(T = 15, T = 25) is determined by the resonant orbits only. In the epicyclic approximation the elongation of the outer rings can be explained through the subtraction and addition of the size of the epicycle to the average radius of the ring. In this case the outer rings must demonstrate the following tendency: the smaller the elongation of the ring, the smaller the residual velocities. Models 1­3 confirm this dependence. The outer rings become more round by the time T = 25 (Fig. 3), that is accompanied with the decreasing velocity amplitudes (Table 3).

surface density of particles (gas + OB) in models 1­3 at different moments. At the initial moment the surface density is constant in the interval R = 1­10 kpc in all models. At the moment T = 15 all density profiles have a well-defined maximum corresponding to the location of the R2 ring (R = 7­8 kpc). It grows faster in model 1 than in models 2 and 3 what indicates the more intense gas outflow in model 1. In all models the maximum corresponding to the R2 ring is shifted 0.5 kpc outside the OLR, while the maximum corresponding to the R1 ring is lying in the vicinity of the outer 4/1 resonance. 4. COMPARISON WITH THE GALAXY

3.3. Kinematics of the Gas Outflow
The main kinematical feature of the outer pseudorings is the preponderance of positive radial velocities against negative ones. All models considered demonstrate the gas outflow at early moments of their evolution. Figure 7 represents the distribution of particles (gas + OB) with the positive and negative residual velocities in model 1 at T = 7.5 (4 bar rotation periods), other models have a very similar velocity distribution. Figure 7a clearly shows that particles with the positive radial velocities concentrate to the spiral arms. Figure 7b demonstrates another kinematical feature of the pseudorings: particles located in the inner parts of the spiral arms (R < ROLR ) have positive azimuthal velocities V , while those forming the outer parts of the spiral arms (R > ROLR ) have only negative ones. Schwarz (1981, 1984) finds that the gas outflow to the periphery plays an important role in the formation of the outer rings. Figure 8 shows the profiles of the
ASTRONOMY LETTERS Vol. 35 No. 9 2009

4.1. Positions of the Stellar-Gas Complexes with Respect to the Outer Rings
Let us suppose that the kinematics of young stars in the solar neighborhood is determined by their position with respect to the outer rings. Observations suggest that the Sun goes b = 15 ­45 behind the bar elongation, so all complexes studied must be located near the segments 3­5 and C­E of the outer rings (Fig. 6). The Perseus, Cygnus, and Carina regions having the negative velocities VR must belong to the ring/pseudoring R2 , but the Sagittarius region with the positive velocities VR --to the ring/pseudoring R1 . There is some problem with the Perseus (R = 8.4 kpc) and Carina (R = 6.5 kpc) regions: they have very different Galactocentric distances R and they both cannot accurately lie on the same ring. Otherwise the ring R2 must be highly elongated:


620

MEL'NIK, RAUTIAINEN

T = 7.5 200

m1 m2 m3

Table 4. Model 1. Average residual velocities Region T = 7.5 Sagittarius Carina Cygnus Local System Perseus T = 15.0 Sagittarius Carina Cygnus Local System Perseus VR , km s- 21 18 ­23 ­15 ­4 24 ­4 ­24 ­24 ­4
1

R , km s- 10 14 3 4 1 6 20 4 3 5

1

V , km s- 10 8 7 3 ­1 0 10 1 ­5 ­5

1

, km s- 12 11 2 2 2 5 9 3 3 7

1

n 164 322 33 63 191 41 91 18 159 207

100

T = 15 200 , kpc­2

Table 5. Model 2. Average residual velocities

100

Region T = 7.5 Sagittarius Carina Cygnus Local System Perseus T = 15.0 Sagittarius Carina Cygnus Local System Perseus

VR , km s- 11 12 ­16 ­13 ­4 10 9 ­16 ­20 ­7

1

R , km s- 3 4 2 2 2 3 6 1 2 6

1

V , km s- 0 0 11 6 ­2 2 1 8 1 ­5

1

, km s- 3 7 2 4 2 3 8 3 3 4

1

n 87 183 49 50 197 98 149 5 127 187

T = 25 200

100
Table 6. Model 3. Average residual velocities

0

2

4

6 R, kpc

8

10

12

Region T = 7.5 Sagittarius Carina Cygnus Local System Perseus T = 15.0 Sagittarius Carina Cygnus Local System Perseus

VR , km s- 13 15 ­17 ­13 ­4 13 7 ­18 ­21 ­8

1

R , km s- 4 6 2 2 1 4 9 5 2 7

1

V , km s- ­1 1 10 6 ­2 1 2 6 1 ­6

1

, km s- 4 9 2 3 1 4 11 2 3 5

1

n 105 246 31 43 179 106 138 12 144 191

Fig. 8. Profiles of the surface density of particles (gas + OB) in models 1­3 at different moments. At the initial moment the surface density is constant in the interval R = 1­10 kpc.

the ratio q of the minor and major axis must be of q 0.6. However, the observations suggest that the rings/pseudorings R2 are not so much elongated: their average axis ratio is q = 0.9 (Buta and Combes, 1996). So the Perseus region must lie a bit outside the R2 ring, while the Carina region--a bit inside it. The kinematics of the Sagittarius region is very important for the definition of a type of the Galactic morphology. The nearly zero residual azimuthal velocities, V = -1 km s-1 , here suggest the presence

of the R1 ring in the Galaxy but not of the spiral arms. Particles located in the inner parts of the spiral arms have only positive velocities V (Fig. 7b) and cannot
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS

621

(a) Y Model 1 VR Model 1 V

(b) Y

X Y Model 2 VR Model 2 V Y

X

X Y Model 3 VR Model 3 V Y

X

X

X

Fig. 9. The boundaries of the stellar-gas complexes overlaid above the distribution of particles with the positive and negative residual velocities in models 1­3 at the moment T = 15. The arrows represent the observed residual velocities. Particles with the positive residual velocities (VR > 5 km s-1 or V > 5 km s-1 ) are indicated by black-gray circles, while particles with the negative ones (VR < -5 km s-1 or V < -5 km s-1 )--by light-gray circles. For more details see the caption to Fig. 4.

reproduce V = -1 km s-1 . Thus, the Galaxy must also include the R1 ring.

velocity dispersions, R and , and the number n of particles (gas + OB) which appear to fall within the complex at a certain time moment. Table 7 represents the model and observed residual velocities for each region. The interval of the model velocities demonstrates the minimal and maximal velocity obtained in models 1­3at T = 15. The Perseus region (R = 8.4 kpc) has the negative values of the average velocities, VR and V , in all models (Table 7). Though the amplitude of the radial velocities at the average distance of the R2 ring (R = 7.4 kpc) is fairly large, fR = 23­31 km s-1 (Table 3, T = 15), its value drops to fR = 4­7 km s-1 at the distance R = 8.4 kpc. The decrease in the velocity amplitude is due to the fact that the Perseus region

4.2. Comparison between Models and Observations To compare models with observations we chose the moment T = 15 when the resonance kinematics is best-defined. Figure 9 shows the boundaries of the stellar-gas complexes overlaid above the distribution of particles with the positive and negative residual velocities at the moment T = 15. The solar position angle is supposed to be b = 45 in models 1 and 3 and b = 40 in model 2 (explanation will be given below). Tables 4­6 exhibit the average residual velocities, VR and V , of particles located within the boundaries of the stellar-gas complexes in models 1­3 at the moments T = 7.5 and 15.0. They also contain the
ASTRONOMY LETTERS Vol. 35 No. 9 2009


622

MEL'NIK, RAUTIAINEN

Table 7. Comparison between models and observations Region Sagittarius Carina Cygnus Local system Perseus VR
mod

,kms

-1

VR

obs

,km s

-1

V

mod

,kms

-1

V

obs

,km s

-1

(+10, +24) (-4, +9) (-24, -16) (-24, -20) (-4, -7)

+11 ± 3 -6 ± 2 -7 ± 3 +5 ± 3 -7 ± 2

(-1, +2) (+1, +10) (+1, +8) (-5, +1) (-5, -6)

-1 ± 2 +5 ± 3 -11 ± 2 0±3 -5 ± 2

is situated a bit outside the R2 ring and has less elongated orbit. Altogether, our models reproduce well the kinematics in the Perseus region: the difference between the model and observed velocities don't exceed 3km s-1 . The Sagittarius region (R = 5.6 kpc) lies near the point D of the R1 ring due to the proper scaling and choice of the solar position angle (Fig. 6). Particles

10 5 0 ­5 ­10 V, km s­1

Sagittarius

m1 m2 m3

5 0 ­5 ­10 ­15 20 30

Perseus

m1 m2 m3

40

b, deg

50

60

70

Fig. 10. (a) The dependence of the model azimuthal velocity V (km s-1 ) in the Sagittarius and Perseus regions (T = 15, models 1­3) on the solar position angle b (deg).

located near the point D have the positive velocity VR and nearly zero velocity V . Model 1 yields a very large value of the radial velocity, VR = 24 km s-1 , but models 2 and 3 create the velocities in a reasonable range, VR = 10­13 km s-1 . In models 2 and 3 the difference between the model and observed velocities don't exceed 3 km s-1 . Interestingly, that at the moment T = 7.5 model 1 produces the large positive azimuthal velocity in the Sagittarius region (V = +10 km s-1 ), but by the moment T = 15 its value drops to V = 0 km s-1 (Table 4). Probably, the large velocity V in model 1 is due to the intense gas outflow at T = 7.5. Models 2 and 3 do not show such feature (Tables 5, 6). The Carina region (R = 6.5 kpc) consists of two groups of particles with different velocities (Fig. 9). One group belongs to the ring R1 , but another--to the ring R2 . The mixture of two streams produces the larger values of the velocity dispersions, R and , here in comparison with the other regions (Tables 4­ 6, T = 15). The velocities of particles in the ring R2 are consistent with observations, so the increase of their relative number in the Carina region brings approaching between the model and observed velocities. All models reproduce the direction of the observed azimuthal residual velocity, but only model 1 reproduces the direction of the radial one. Model 1 has more elongated ring R2 , what causes more particles of the ring R2 to fall within the boundaries of the Carina region. Note that a small rescaling of models can significantly change the average velocity of particles in the Carina region. The Cygnus region (R = 6.9 kpc) is located between two outer rings where only a few particles are located. Although the direction of the model radial velocity agrees with observations, its absolute value is too high. As for the azimuthal velocities, none of our models can reproduce the observed negative azimuthal velocity here (Table 7). The Local System (R = 7.4 kpc) lies in the vicinity of the R2 ring in all models. It contains more particles at the moment T = 15 than at T = 7.5. In the context of morphology the moment T = 7.5 is more
ASTRONOMY LETTERS Vol. 35 No. 9 2009


KINEMATICS OF THE OUTER PSEUDORINGS

623

suitable for a comparison with observations, because at this moment the break in the R2 ring is fairly wide yet (Fig. 3). Our models reproduce well the nearly zero azimuthal velocity, but they cannot create the observed positive radial velocity here (Table 7). Thus, models 2 and 3 reproduce well the values and directions of the residual velocities in the Perseus and Sagittarius regions, but model 1 reproduces the directions of the residual velocities in the Perseus, Sagittarius, and Carina regions. Probably, the kinematics of young stars in the Perseus, Sagittarius, and Carina regions is determined mostly by the resonance, while other processes are dominant in the Cygnus region and in the Local System. The value of the model azimuthal velocity V in the Sagittarius and Perseus regions is especially sensitive to the choice of the solar position angle b . An increase of the angle b causes a decrease of the velocity V in the Sagittarius region and an opposite effect in the Perseus one. Figure 10 exhibits the variations of the average velocity V in the Sagittarius and Perseus regions with the change of the angle b . It is seen that all models reproduce the value of V = -1 km s-1 in the Sagittarius region under b = 45 ­48 . But in this range model 2 creates the small absolute values of the velocity V in the Perseus region, |V | = 3 km s-1 . To get the larger value of |V | in Perseus we adopted b = 40 for model 2. 5. CONCLUSIONS The kinematics of the outer rings and pseudorings is determined by two processes: the resonance tuning and the gas outflow. The resonance kinematics is clearly observed in the pure rings while the kinematics of the outflow is manifested itself in the pseudorings. The resonance between the epicyclic and orbital motion in the reference frame rotating with the bar speed adjusts the epicyclic motions of particles in accordance with the bar rotation. This adjustment creates the systematical non-circular motions, whose direction depends on the position angle of a point with respect to the bar elongation and on the class of the outer ring. Models of the Galaxy with the R1 R2 pseudoring reproduce well the radial and azimuthal components of the residual velocities of OB-associations in the Perseus and Sagittarius regions: the difference between the model and observed velocities does not exceed 3 km s-1 (models 2 and 3, T = 15). The kinematics in the Perseus region indicates the presence of the R2 ring in the Galaxy, while the velocities in the Sagittarius region suggest the existence of the R1 ring. The azimuthal velocities in the Sagittarius region accurately defines the solar position angle with respect to the bar elongation, b = 45 ± 5 . Besides,
ASTRONOMY LETTERS Vol. 35 No. 9 2009

model 1 reproduces the directions of the radial and azimuthal residual velocities in the Perseus, Sagittarius, and Carina regions. Probably, the kinematics of young stars in the Perseus, Sagittarius, and Carina regions is determined mostly by the resonant orbits. Our models have nearly flat rotation curve. The ring R1 R2 is forming after several bar rotation periods. The R1 ring lies in the vicinity of the outer 4/1-resonance, while the R2 pseudoring is shifted 0.5 kpc outside the OLR. We found that the mixed R1 R2 morphology appears when the bar is not very strong: in the case of a strong bar the R1 ring disappears. The strong bar and the initial distribution of particles in the disk of a large radius accelerate the formation of the R2 rings. The major semi-axis of the bar in our models has the value of a = 3.8­4.0 kpc which amounts 55% of the solar Galactocentric distance, R0 = 7.1 kpc, adopted here. The relatively large size of the bar and the large value of the solar position angle, b = 45 ± 5 , agree well with studies in which authors find the presence of a long bar in the Galaxy (Weiner and Sellwood, 1999; Benjamin et al., 2005; Cabrera-Lavers et al., 2007; and other papers). The model of the Galaxy with the R1 R2 pseudoring can explain some large-scale morphological features of the Galactic spiral structure. The so-called Carina spiral arm (Russeil, 2003) falls nicely onto the segment of the R2 ring. Note that two outer rings which are stretched perpendicular to each other can look like a 4-armed spiral pattern, especially if the ascending parts of the rings are brighter than the descending ones. But this spiral pattern appears not as a result of the spiral perturbation of the disk potential but due to the existence of specific orbits in the barred galaxies. ACKNOWLEDGMENTS We wish to thank Heikki Salo for providing his simulation program. This work was partly supported by the Russian Foundation for Basic Research (project no. 06-02-16077) and the Council for the Program of Support for Leading Scientific Schools (project no. NSh-433.2008.2). REFERENCES
1. E. Athanassoula, Mon. Not. R. Astron. Soc. 259, 328 (1992). 2. V. S. Avedisova and J. Palous, BAICz 40, 42 (1989). 3. R. A. Benjamin, E. Churchwell, B. L. Babler, et al., Astrophys. J. 630, L149 (2005). 4. L. N. Berdnikov, A. K. Dambis, and O. V. Vozyakova, Astron. Astrophys., Suppl. Ser. 143, 211 (2000).


624

MEL'NIK, RAUTIAINEN 27. R. M. Humphreys, Astrophys. J. 206, 114 (1976). 28. R. M. Humphreys, in IAU Symp. No. 84: The Large-Scale Characteristics of the Galaxy, Ed. by W. B. Burton (Reidel, Dordrecht, 1979) p. 93. 29. A. J. Kalnajs, in Dynamics of Disc Galaxies, Ed. ¨ ¨ by B. Sundelius (Goteborgs Univ., Goteborg, 1991), p. 323. 30. J. Kovalevsky, Modern Astrometry (Berlin, New York: Springer, 2002). 31. C. C. Lin, in IAU Symp. No. 38: The Spiral Structure of our Galaxy, Ed. by W. Becker and G. Contopoulos (Reidel, Dordrecht, 1970) p. 377. 32. C. C. Lin, C. Yuan, and F. H. Shu, Astrophys. J. 155, 721 (1969). 33. A. M. Mel'nik, Pis'ma Astron. Zh. 29, 349 (2003) [Astron. Lett. 29, 304 (2003)]. 34. A. M. Mel'nik, Pis'ma Astron. Zh. 32, 9 (2006) [Astron. Lett. 32, 7 (2006)]. 35. A. M. Mel'nik, A. K. Dambis, and A. S. Rastorguev, Pis'ma Astron. Zh. 27, 611 (2001) [Astron. Lett. 27, 521 (2001)]. 36. P. Rautiainen and H. Salo, Astron. Astrophys. 348, 737 (1999). 37. P. Rautiainen and H. Salo, Astron. Astrophys. 362, 465 (2000). 38. P. Rautiainen, H. Salo, and E. Laurikainen, Astrophys. J. 631, L129 (2005). 39. W. W. Roberts, Astrophys. J. 158, 123 (1969). ´ 40. M. Romero-Gomez, E. Athanassoula, J. J. Masde´ mont, and C. Garci´ a-Gomez, Astron. Astrophys. 472, 63 (2007). 41. D. Russeil, Astron. Astrophys. 397, 133 (2003). 42. H. Salo, Astron. Astrophys. 243, 118 (1991). 43. M. P. Schwarz, Astrophys. J. 247, 77 (1981). 44. M. P. Schwarz, Mon. Not. R. Astron. Soc. 209, 93 (1984). 45. T. G. Sitnik, Pis'ma Astron. Zh. 29, 356 (2003) [Astron. Lett. 29, 311 (2003)]. 46. T. G. Sitnik and A. M. Mel'nik, Pis'ma Astron. Zh. 22, 471 (1996) [Astron. Lett. 22, 422 (1996)]. 47. The Hipparcos and Tycho Catalogs, ESA SP-1200 (1997). ´ 48. J. P. Vallee, Astron. J. 130, 569 (2005). 49. B. J. Weiner and J. A. Sellwood, Astrophys. J. 524, 112 (1999).

5. J. Binney and S. Tremaine, Galactic Dynamics, (Princeton Univ. Press, Princeton, NJ, 2008). 6. C. Blaha and R. M. Humphreys, Astron. J. 98, 1598 (1989). 7. J. Brand and L. Blitz, Astron. Astrophys. 275, 67 (1993). 8. W. B. Burton, Astron. Astrophys. 10, 76 (1971). 9. W. B. Burton and T. M. Bania, Astron. Astrophys. 33, 425 (1974). 10. W. B. Burton and M. A. Gordon, Astron. Astrophys. 63, 7 (1978). 11. R. Buta, Astrophys. J, Suppl. Ser. 96, 39 (1995). 12. R. Buta and F. Combes, Fund. Cosmic Physics 17,95 (1996). 13. R. Buta, H. G. Corwin, and S. C. Odewahn, The de Vaucouleurs Atlas of Galaxies (Cambridge Univ. Press, 2007). 14. R. Buta and D. A. Crocker, Astron. J. 102, 1715 (1991). 15. G. Byrd, P. Rautiainen, H. Salo, R. Buta, and D. A. Crocker, Astron. J. 108, 476 (1994). 16. A. Cabrera-Lavers, P. L. Hammersley, C. Gonzalez-Fernandez, et al. Astron. Astrophys. 465, 825 (2007). 17. G. Contopoulos and Th. Papayannopoulos, Astron. Astrophys. 92, 33 (1980). 18. A. K. Dambis, A. M. Mel'nik, and A. S. Rastorguev, Pis'ma Astron. Zh. 21, 331 (1995) [Astron. Lett. 21, 291 (1995)]. 19. A. K. Dambis, A. M. Mel'nik, and A. S. Rastorguev, Pis'ma Astron. Zh. 27, 68 (2001) [Astron. Lett. 27, 58 (2001)]. 20. Y. N. Efremov and T. G. Sitnik, Pis'ma Astron. Zh. 14, 817 (1988) [Astron. Lett. 14, 347 (1988)]. 21. Y. N. Efremov, Astron. Astrophys. Trans. 15, 3 (1998). 22. P. Englmaier and O. Gerhard, Celestial Mechanics and Dynamical Astronomy, 94, 369 (2006). 23. C. D. Garmany and R. E. Stencel, Astron. Astrophys., Suppl. Ser. 94, 211 (1992). 24. Y. M. Georgelin and Y. P. Georgelin, Astron. Astrophys. 49, 57 (1976). 25. E. V. Glushkova, A. K. Dambis, A. M. Mel'nik, and A. S. Rastorguev, Astron. Astrophys. 329, 514 (1998). 26. H. J. Habing, M. N. Sevenster, M. Messineo, G. van de Ven, and K. Kuijken, Astron. Astrophys. 458, 151 (2006).

ASTRONOMY LETTERS

Vol. 35 No. 9

2009