Issue 
A&A
Volume 541, May 2012



Article Number  A64  
Number of page(s)  11  
Section  Galactic structure, stellar clusters and populations  
DOI  https://doi.org/10.1051/00046361/201117347  
Published online  30 April 2012 
The origin and orbit of the old, metalrich, open cluster NGC 6791
Insights from kinematics
^{1} ESO, Alonso de Cordova 3107, Casilla 19001, Santiago, Chile
email: ljilkova@eso.org; gcarraro@eso.org
^{2} Department of Theoretical Physics and Astrophysics, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
^{3} Dipartimento di Astronomia, Universitá di Padova, Vicolo Osservatorio 3, 35122 Padova, Italy
^{4} Astronomical Institute, Academy of Sciences of the Czech Republic, Boční II 1401/1a, 141 31 Prague, Czech Republic
email: bruno@ig.cas.cz
^{5} Astronomical Institute, Faculty of Mathematics and Physics, Charles University in Prague, Ke Karlovu 3, 121 16 Prague, Czech Republic
^{6} LeibnizInstitut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
email: iminchev@aip.de
Received: 25 May 2011
Accepted: 1 March 2012
Context.NGC 6791 is a unique stellar system among Galactic open clusters, which is at the same time one of the oldest open clusters and the most metal rich. It is located inside the solar circle, harbors a large population of binary stars, and possibly experienced prolonged star formation. The combination of all these properties is puzzling and poses the intriguing question of its origin.
Aims. One possible scenario is that the cluster formed close to the Galactic Center and later migrated outward to its current location. In this work we study the cluster’s orbit and investigate the possible migration processes that may have displaced NGC 6791 to its presentday position, under the assumption that it actually formed in the inner disk.
Methods. To this aim we performed integrations of NGC 6791’s orbit in a potential consistent with the main Milky Way parameters. In addition to analytical expressions for halo, bulge and disk, we also consider the effect of bar and spiral arm perturbations, which are expected to be very important for the disk dynamical evolution, especially inside the solar circle. Starting from stateofthe art initial conditions for NGC 6791, we calculated 1000 orbits back in time for about 1 Gyr turning different nonaxisymmetric components of the global potential on and off. We then compared statistical estimates of the cluster’s recent orbital parameters with the orbital parameters of 10^{4} testparticles originating close to the Galactic Center (with initial galocentric radii in the range of 3–5 kpc) and undergoing radial migration during 8 Gyr of forward integration.
Results. We find that a model that incorporates a strong bar and spiral arm perturbations can indeed be responsible for the migration of NGC 6791 from the inner disk (galocentric radii of 3–5 kpc) to its presentday location. Such a model can provide orbital parameters that are close enough to the observed ones. However, the probability of this scenario as it results from our investigations is very low.
Key words: Galaxy: kinematics and dynamics / open clusters and associations: general / Galaxy: structure / open clusters and associations: individual: NGC 6791
© ESO, 2012
1. Introduction
Star clusters – globular and open – have long been recognized as ideal benchmarks for stellar evolution theories primarily because they are made of a large number of stars that share the same age, distance, and metallicity. This allows comparisons with stellar models in a statistical sense. The only exception has routinely been ω Cen, which has been known for quite some time to display a significant spread in color in its red giant branch (RGB) and inhomogeneities in chemical elements (Piotto et al. 2005; Villanova et al. 2010, and references therein).
However, in the past decades, signatures of multiple stellar populations – with different age and metallicity – have been detected in virtually all Milky Way (MW) globulars by means of very accurate photometry and spectroscopy (see Piotto 2009, for a review on the subject). Typical photometric signatures are the presence of multiple distinct main sequences (MS) or RGBs in the same cluster color magnitude diagram (CMD), while spectroscopic signatures are normally in the form of spreads in the abundance of iron and other elements much larger than observational uncertainties (Milone et al. 2009).
The rich, intermediateage star clusters of the Magellanic Clouds (MCs) do not seem to be an exception (Milone et al. 2009), although in this case signatures of prolonged star formation are more common than evidences of discrete episodes of star formation. Indeed, in the CMD of MCs clusters the MS is much broader than expected by photometric errors and binary stars, but does not separate into discrete sequences unlike in Galactic globulars.
As of now, the only MW open cluster where evidence of prolonged star formation has been found is NGC 6791 . The upper MS and turnoff region of this cluster are much wider than expected from photometric errors and binary stars alone. By separating the stars lying inside the cluster core and outside it, Twarog et al. (2011) found that inner stars occupy the red edge of the upper MS, while outer stars occupy the blue edge. The most viable explanation for this occurrence is an age spread of ~1 Gyr between the inner and the outer regions of the cluster, in the sense that inner stars formed before outer stars. Star formation therefore started probably in the center of the cluster and continued outward during quite a long time.
NGC 6791 is unique in the MW because it is one of the most massive old open clusters with an age of ~8 Gyr, and is extremely metalrich ([Fe/H] ~+0.40, Carraro et al. 2006, Origlia et al. 2006). No other open cluster in the MW is that rich in metals. This combination of age and metallicity is unique in the whole population of Galactic open clusters. While the RGB of NGC 6791 is broad in color (Janes 1984), no statistically significant evidence for a metallicity spread within the cluster – which would support a selfenrichment scenario – has been reported so far. Within observational errors, stars in NGC 6791 show the same metal abundance. As mentioned by Twarog et al. (2011), since the CMDbased age spread estimated for populous clusters is usually less than the 0.7 Gyr found for NGC 419 by Rubele et al. (2010), the CMD of NGC 6791 may have resembled that of the Magellanic Cloud clusters 7 Gyr ago.
The cluster’s current location is at 4 kpc from the Sun, close to the solar circle, at 8 kpc from the Galactic Center (GC) and about 0.8 kpc above the plane (derived from Brogaard et al. 2011). The metallicity of this cluster is at odds with the radial abundance gradient defined by the other MW open clusters (Carraro et al. 2006).
This wealth of unique properties clearly poses the basic question of the origin of this system. So far, to answer to this question only arguments based on the observed chemical characteristics of the cluster have been used. If this system is so metal rich, it must have formed very close to the Galactic bulge and then been deposited where it is now by some migration process. Indeed, only in the inner regions of the MW can star formation and metal enrichment be fast and strong enough to generate a stellar system like NGC 6791.
Alternatively, one can argue that NGC 6791 is the leftover of a dwarf elliptical that was disrupted by the MW some time ago and was then engulfed by the Galaxy, losing most of its mass by tidal interaction. This seems to be supported by the many extreme blue horizontal branch stars discovered in the cluster (Kaluzny & Rucinski 1995; Liebert et al. 1994; Buson et al. 2006) and by the recent measure of a significant UVupturn in the spectral energy distribution of the cluster (Buzzoni et al. 2012). Such a UVupturn is only found in moderately UVenhanced ellipticals and has the same strength as in, e.g., NGC 4552 and NGC 4486, which emit ~2% of their light at a shorter wavelength than 2500 Å. This feature does not have any counterpart in the population of Galactic open clusters.
The study of the NGC 6791 Galactic orbit did not reveal much to clarify the scenario so far, since depending on the quality of presentday kinematic data, simple backward calculation could produce different orbits – a regular orbit with different eccentricities e ~ 0.5 or ~0.3 by Bedin et al. (2006), and Wu et al. (2009), respectively, both reaching galocentric radii of ~10 kpc, and an orbit by Carraro et al. (2006) extending to more than 20 kpc in the Galactic disk.
In this paper we present a numerical study of the NGC 6791 Galactic orbit that includes the combined effect of the Galactic bar and spiral arms into an analytical MW model. At the same time, we aim at investigating in detail the scenario in which NGC 6791 formed close to the GC, where star formation might have been rapid and efficient enough to generate it. This scenario requires some radial migration process to move the cluster’s orbit outward in the Galactic disk.
Several different dynamical mechanisms triggering radial mixing are discussed in the recent literature:

1.
the resonant scattering by transient spiral structure was firstlyintroduced by Sellwood & Binney (2002). Starscan move inward and outward at the corotation radius whilepreserving their circular orbit. During the evolution of the disk theangular velocity (the corotation radius) changes and the spiralscan induce radial migration within the whole disk. Theimportance of the transient spiral structure radial migrationphenomenon was confirmed by Roškar et al.(2008; see also Roškar et al. 2011).

2.
Minchev & Famaey (2010) described another diffusion mechanism resulting from the resonance overlap of multiple rotating patterns, the Galactic bar, and spiral arms. Using testparticle simulations, these authors found that the overlap mechanism can be up to an order of magnitude more effective than the one due to transient spiral structure (for the model consistent with the MW it takes ~3 Gyr to produce the same radial mixing as in ~9 Gyr by the transients, see also Minchev & Famaey 2011). This migration mechanism is a strong function of the strength of MW bar and spiral arms. However, this dependency is nonlinear. The characteristics of the resonance overlap mechanism were later confirmed by Minchev et al. (2011) via a fully selfconsistent, treeSPH Nbody simulations. Among other results, these authors also found that the presence of a gaseous component makes the migration even more efficient.

3.
Small satellites have been found to effectively mix the outer disks of galaxies. Quillen et al. (2009) used disklike distributed testparticles perturbed by a lowmass satellite (of mass a few times 10^{9} M_{⊙}). These authors find that the disk stars from the outer disk can move to low eccentricity orbits inside their birth radii. Bird et al. (2012) compared radial migration in several fully selfconsistent Nbody simulations both with and without interaction with satellites. Among other conclusions, these authors showed that the radial migration induced by satellites is a distinct mechanism that does not happen in isolated systems.
During the past decades it has become widely accepted that the MW has a central bar. Evidence for this comes from observations of the inner Galaxy (e.g., NIR photometry – Blitz & Spergel 1991; Dwek et al. 1995; Binney et al. 1997; gas kinematics – Englmaier & Gerhard 1999; and microlensing – Rattenbury et al. 2007) and dynamical effects near the solar circle (Dehnen 2000; Minchev et al. 2007, 2010). It is, therefore, very important to consider the MW’s bar in orbital integration studies, especially in the inner disk, but also near the Sun.
While for the Sellwood & Binney (2002) mechanism to work shortlived transient spirals are required, in barred galaxies, such as the MW, spirals are more likely coupled with the bar, and are thus longerlived (Binney & Tremaine 2008; Quillen et al. 2011). Therefore, in this paper we will concentrate on the migration induced by spiralbar interaction, rather than transient spirals. The migration induced by orbiting satellites will also be unimportant for the inner disk.
In Sect. 2 we describe the adopted analytical Galactic model and introduce spiral arm and bar perturbation. The initial conditions and the results of the backward orbit integration are presented in Sect. 3. We concentrate on radial migration in Sect. 4 and, finally, discuss the results of this work in Sect. 5.
2. Galactic models and orbit integrator
To compute any Galactic orbit, it is necessary to adopt a model for the MW gravitational potential. The potential should reproduce the observed mass density of the MW. While the general picture of the today’s MW structure is more or less well established and supported by observations, a more detailed description is still under lively discussion, see e.g., Schönrich et al. (2010) or Coşkunoǧlu et al. (2011) for the recent analysis of the velocity of the local standard of rest (LSR), or Reid et al. (2009) for the uptodate description of MW rotation curve. These questions naturally bring difficultly quantifiable uncertainties for the construction of any MW model.
We chose to use a relatively simple model composed of several components. The axisymmetric timeindependent part was modeled by a bulge, disk and halo. Additionally, we included a nonaxisymmetric component – bar and spiral arms – that can act separately or together. For each component, an analytic expression for the density is provided, from which the corresponding potential and force can be derived.
The model is described in the following subsection. The main idea of our study is to qualitatively estimate the viability of the migration induced by the bar and spiral arms resonance overlap mechanism relative to the migration in simpler potentials, either purely axisymmetric or containing just one type of nonaxisymmetric perturbation.
2.1. Axisymmetric model
For the timeindependent axisymmetric component of the Galactic potential we used the model introduced by Flynn et al. (1996) that was also used by Gardner & Flynn (2010)^{1}. This model has three components – bulge, disk, and dark halo. The bulge is modeled as a superposition of two Plummer spheres (Plummer 1911). The disk is modeled as a superposition of three MiyamotoNagai disks (Miyamoto & Nagai 1975). Finally, the dark halo is modeled as a spherical logarithmic potential. For clarity, we give the equation for each component of the axisymmetric potential: where r is the galocentric radius (r^{2} = R^{2} + z^{2}), G is the gravitational constant, M_{C1} and M_{C2}, r_{C1} and r_{C2}, are the masses and scale radii of the bulge components, respectively. M_{Dn} and a_{n} are masses and scalelengths of the three disk components (n goes from 1 to 3). The parameter b is related to the scaleheight of the disk and is the same for all three components of the disk. The halo has core radius r_{0} and V_{H} is the asymptotic circular velocity at large radii (relative to r_{0}).
Flynn et al. (1996) chose the parameters of their model to reproduce the current MW observations at that time. We modified the values of some parameters to achieve a model consistent with more recent MW observations. As for the recent observations of the rotation curve at the solar galocentric distance R_{⊙}, Reid et al. (2009) estimated R_{⊙} = (8.4 ± 0.6) kpc and the solar circular rotation speed Θ_{⊙} = (254 ± 16) km s^{1}. To derive these values, they used the solar motion with respect to the LSR, v_{ ⊙ ,LSR}, as determined by Dehnen & Binney (1998) from the Hipparcos stars. The same Hipparcos data set was later reexamined by Schönrich et al. (2010), who obtained slightly different values for the v_{ ⊙ ,LSR} components. Namely, they found the velocity component in the direction of the Galactic rotation to be about 7 km s^{1} higher than Dehnen & Binney (1998). Considering these more recent results, we changed the original Flynn et al. (1996) values for the halo circular velocity V_{H} and the disk parameters (masses and scale radii) to obtain Θ_{⊙} = 243 km s^{1} at R_{⊙} = 8.4 kpc. These values are consistent with the Θ_{⊙} derived from the apparent proper motion of Sgr A* (Reid & Brunthaler 2004) using R_{⊙} = 8.4 kpc and the Schönrich et al. (2010) solar motion.
The superposition of three MiyamotoNagai disks builds up a disk with exponential falloff of the radial surface density on a wide range of radii. Similar as in Gardner & Flynn (2010), we chose the disk parameters to obtain the scalelength of the exponential falloff of the disk consistent with recent observations. Our disk gives a radial surface density profile approximately corresponding to the exponential falloff with the scalelength of 2.6 kpc for radii of 5–18 kpc, which is consistent with recent observations (Jurić et al. 2008).
Values for all parameters of the axisymmetrical component are given in Table 1.
Parameters of the axisymmetrical component of today’s MW model.
The above described axisymmetric potential was also used in the models including the nonaxisymmetric components. However, when we included the bar perturbation, the mass M_{C2} of one of the bulge components was reduced and the bar replaced most of the bulge mass. We kept the axisymmetric components unaltered when we turn the spiral arms perturbation on (see Sect. 2.3).
2.2. Galactic bar
For the bar model we adopted the Ferrers potential (see paragraph 2.5.3 (a) of Binney & Tremaine 2008,and we chose n = 2) of an inhomogeneous triaxial ellipsoid (Pfenniger 1984). The bar model is defined by the bar length, axes ratios, mass, and the angular velocity. For all parameters (listed in Table 2) we adopted the same values as reported by Pichardo et al. (2004) and refer the reader to their original paper for a justification of the various choices.
2.3. Spiral arms
The parameters of the spiral structure are much more uncertain than for the bar. Most observational studies suggest twoarmed or fourarmed spiral pattern, or their superposition (see for example Vallée 2005, 2008; Allen et al. 2008; Grosbol et al. 2011).
We modeled the gravitational potential perturbation induced by spiral arms as described by Cox & Gómez (2002). We used a twoarmed spiral pattern with their basic sinusoidal arms. The perturbation has the shape of a logarithmic spiral. The shape of the spiral arms model is defined by five parameters – radial scale length, scale height, pitch angle, amplitude of density perturbation (mass density of the arms in the Galactic plane at the solar radius), and angular velocity. The vertical structure of the arms’ mass density has the form sech^{2}(z/z_{s}), where z_{s} is a scale height. We chose the values of parameters consistent with a current picture of the MW spiral structure – see Table 2.
Drimmel & Spergel (2001) present a threedimensional model for the MW fitted to the COBE farinfrared and nearinfrared emission. They found that twoarmed structure dominates the nearinfrared emission and the arminterarm density contrast of their modeled stellar spiral arms has a value of about 1.8 at the solar radius (see Fig. 15 in Drimmel & Spergel 2001). Our choice of parameters for Cox & Gómez (2002) gives 2.0 for the same ratio. The value is also consistent with the observations of external galaxies with a weaker spiral structure (see, e.g., Block et al. 2004).
In addition to the parameters defining the shape of spiral arms, we needed to specify the initial orientation of the arms. We followed the result of the study of the COBE/DIRBE Kband emission made by Drimmel (2000), which is consistent with twoarmed spiral model of the MW.
Fig. 1 Maps of strength of nonaxisymmetric components in the combined (bar and spiral arms) today’s MW model (noted as MW1 in Sect. 4). Color scales map the ratio between radial (left) and tangential (right) force resulting from the nonmonopole part of the gravitational potential to the radial force due to the monopole part at the same radius. Maps are plotted for the Galactic plane. ⊙ symbol marks the position of the Sun. 

Open with DEXTER 
2.4. Strength of rotating patterns
We estimated the contribution of the nonaxisymmetric components by computing the force resulting from the bar or/and spiral arms’ gravitational potentials in the disk plane, following the approach of Combes & Sanders (1981) (see also Buta & Block 2001). Figure 1 shows maps of the ratio between the radial, resp. the tangential, force due to the nonmonopole part of the combined (bar and spiral arms) potential to the radial force due to the monopole part at the same radius.
The strength of the bar is often characterized by the bar strength parameter, which is defined as the maximal value of the ratio between the bar’s tangential force to the monopole radial force at the same radius. Our choice of bar model and parameters gives a bar strength parameter of 0.4. The model of spiral arms gives a value of 0.05 (using the tangential force from the spiral arms instead of the bar in the previous definition).
2.5. Orbit calculation
Once a model of gravitational potential is adopted, the corresponding equations of motion can be constructed. We used Hamilton’s equations, which give a set of firstorder differential equations that we integrated using a BulirschStoer integrator with adaptive timesteps (Press et al. 1992). The relative change in the Jacobi constant for the models with a single rotating pattern is estimated to be on the order of 10^{10}.
To see and understand the different orbits caused by the bar or/and spiral arms we calculated orbits in four versions of the MW model – axisymmetric potential, potential with the bar but not spiral arms, with spiral arms and no bar and, finally, a combined model with both bar and spiral arms.
3. Recent orbit from backward integration
To investigate the shape and properties of the recent orbit of NGC 6791, we integrated its orbit backward for 1 Gyr, starting from the cluster’s current position and velocity vectors. These initial vectors were obtained from the most recent observational estimates we found in the literature: equatorial celestial coordinates (α and δ for J2000.0 equinox), distance to the Sun d_{⊙}, radial velocity v_{r}, and PM.
Input observational data used for the calculation of the initial conditions – equatorial coordinates α, δ (J2000), radial velocity v_{r}, PM μ_{α}cosδ, μ_{δ}, and distance to the Sun d_{⊙}; together with their references.
3.1. Initial conditions
As mentioned above, NGC 6791 has been the subject of intense study in the past. However, we rely in this study on the most recent values for the aforementioned parameters, which we list in Table 3 together with their source. So that they can be used as initial conditions for orbit integration, the input quantities reported in Table 3 need to be transformed into the appropriate coordinate system.
We integrated the orbits in a cartesian galocentric righthanded coordinates system: the x axis points from the GC outward, in the direction of the Sun; the y axis points in antidirection of the LSR motion; and the z axis toward the North Galactic Pole. The position of the Sun is given by vector (x,y,z)_{⊙} = (8.4,0,0) kpc in this coordinate system. The transformation of the observational data to the Cartesian coordinate system centered on the Sun was performed following step by step the Johnson & Soderblom (1987) algorithm updated to the International Celestial Reference System (the coordinates of the North Galactic Pole α = 192.85948°, δ = 27.12825°). Adopting the solar motion with respect to the LSR (v_{ ⊙ ,LSR}) from Schönrich et al. (2010): (U,V,W)_{⊙} = (11.1,12.24,7.25) km s^{1} (righthanded system, with U in the direction to the GC and V in the Galactic rotation direction) – the velocity vector with respect to LSR was obtained. Finally, using the Sun galocentric distance of 8.4 kpc and the LSR rotation velocity of 243 km s^{1} (see Sect. 2.1), the position and velocity vectors were transformed into the coordinate system of integrations (denoted as x, y, z, v_{x}, v_{y}, v_{z}).
Fig. 2 Recent portion of Galactic orbit of NGC 6791 – orbit projections into the Galactic and meridional plane are plotted in the left and right column, respectively. Orbit in axisymmetric potential is plotted with the blue line. The red line shows the combined model (bar and spiral arms). The initial conditions (see Table 4) are given by mean observational input data (Table 3), the black dot marks the initial position. Integration was backward for time of 1 Gyr. 

Open with DEXTER 
Fig. 3 Distributions of cluster’s recent orbital parameters. The distributions for perigalacticon R_{p}, apogalacticon R_{a}, vertical height of orbit z_{max}, and eccentricity e are plotted in columns. Different colors show results for different MW models: axisymmetric model – light blue, model with bar – black, model with spiral arms – yellow, and red for the bar and spiral arms combined model. 

Open with DEXTER 
The input observational data come naturally with some associated uncertainties, which translate into uncertainties in the initial condition vectors. To take the observational uncertainties into account, we followed Dinescu et al. (1999) (see also Wu et al. 2009; Vande Putte et al. 2010), and calculated a set of 1000 initial condition realizations. The initial velocity and position vectors were extracted from normally distributed values of radial velocity, PM, and distance to the Sun. The dispersions of a given distribution was taken as the observational uncertainty of the input data from Table 3. We integrated an orbit for each initial condition vector backward for 1 Gyr. The initial positions and velocities (x, y, z, v_{x}, v_{y}, v_{z}) derived from the mean observational input data are given in Table 4.
3.2. Orbits and orbital parameters
Galactic orbits are routinely characterized by a set of socalled orbital parameters: the perigalacticon R_{p}, the apogalacticon R_{a}, the vertical amplitude z_{max} defined as the maximal distance from the Galactic plane and, finally, the eccentricity e defined as (R_{a} − R_{p})/(R_{a} + R_{p}). Depending on the type of orbit, values of orbital parameters can differ for individual revolutions.
For irregular orbits orbital parameters change with time as the orbit evolves. This is particularly the case if the orbit is near a resonance for a single pattern or near a resonance overlap for bar + spiral arms. Since chaotic behavior of orbits is expected in these regions, it is impossible to recover/predict an orbit on longer timescales. Therefore, to see the current orbital history, we calculated the orbital parameters for the most recent revolution (defined by azimuthal change of 2π) for each set of initial conditions. The distributions of these orbital parameters are shown in Fig. 3.
To highlight and quantify the influence of the various nonaxisymmetric components, we integrated orbits in all four MW models described in Sect. 2. The orbits (backward integration for time of 1 Gyr) in axisymmetric and combined models given by the mean values of observational input data (see Table 4 for corresponding initial conditions) are plotted in Fig. 2 – projection into the Galactic and meridional plane on the left and right, respectively.
The orbits in models including only bar or spiral arms differ only slightly (around 10%) from the results of the axisymmetric model (see distributions of orbital parameters in Fig. 3). With increasing model complexity there is a trend of having larger apogalacticon and lower perigalacticon, which, in turn, yields a higher eccentricity. This is never higher than 0.33, however, which seems to make highly eccentric orbits for NGC 6791 implausible.
According to these calculations, NGC 6791 is not expected to move much outside the solar circle, at odds with previous suggestions (see next Sect. 3.3). The orbit is presently located close to its perigalacticon and this remains inside the solar circle on average, but never gets very close to the bulge. However, a scenario in which NGC 6791 formed close to the bulge and then migrated outward deserves a closer look, which we shall present in Sect. 4.
3.3. Comparison with previous orbit integration
To our knowledge, the orbit of NGC 6791 was studied three times – by Carraro et al. (2006), Bedin et al. (2006), and by Wu et al. (2009).
Carraro et al. (2006) used complex potential derived from an Nbody gasdynamical model (Fux 1997, 1999). They obtained an eccentric orbit (e = 0.59) extending beyond the solar circle. However, the initial conditions were derived from absolute PM and radial velocity values available at that moment and admittedly less precise than the ones we are using here, and therefore it cannot be compared directly with the results presented here.
Bedin et al. (2006) used a purely axisymmetric Allen & Santillan (1991) model with the most uptodate PM and radial velocity (the same values as we use here). However, reviewing the details of the initial conditions setup, we identified a problem (misdefinition of the x axis orientation), which prevents a proper comparison with our results.
Wu et al. (2009) calculated the orbit using the Allen & Santillan (1991) and two other axisymmetric potential models. As Carraro et al. (2006), they derived the initial conditions from less precise data than we use here, but still their orbital parameters for the Allen & Santillan (1991) potential are within errorbars consistent with results given in this study.
In any case, the present study supersedes any previous one, in the sense that statistical analysis is performed for the orbital parameters adopting different Galactic model versions.
4. Forward orbit integration and radial migration
Parameters of initial condition distributions for the forward orbits integration.
The inner disk, close to the bulge, is a highdensity region where star formation is very efficient and metal enrichment fast (Bensby et al. 2010). It is therefore conceivable to imagine that NGC 6791 could have formed there and then diffused outward. To estimate the probability and efficiency of radial migration we carried out a forward orbits integration. More specifically, we integrated forward in time a number of orbits with an initial location close to the GC. By comparing their orbital properties with the recent NGC 6791 orbital parameters obtained in Sect. 3, we attempted to estimate the probability that radial migration moved the NGC 6791 orbit outward from the inner Galactic regions.
We followed 10^{4} initially closetocircular orbits distributed between 3 and 5 kpc from the GC. Initial positions of orbits were equidistantly separated in azimuth. Then, the radius was generated randomly for each azimuth, from a uniform distribution between 3–5 kpc (see Table 5 for a detailed description of initial condition distributions in phasespace).
To quantify the relative influence of different nonaxisymmetric components, we followed the orbits in all four flavors of the Galactic potential model described in Sect. 2 (axisymmetric, barred, spiral, combined). When included, the amplitude of the nonaxisymmetric components grows simultaneously in time for 0.4 Gyr from 0 to its maximum value with the timedependency described by Dehnen (2000, given by his Eq. (4)). The gradual growth of nonaxisymmetric perturbations is introduced to allow a progressive adaptation of initially nearly circular orbits. The growth of a real bar and real spiral arms is certainly much more complex and not simultaneous, but modeling these time dependencies is beyond the scope of our model. We tested that our results are not significantly sensitive, apart from the first few hundreds of Myr, to the above gradual growth when compared to orbit integrations in which the bar and/or spiral arms were switched on abruptly to their full strength from the beginning.
According to Minchev & Famaey (2010) 3 Gyr should be sufficient for disk radial mixing. On the other hand, the age of the NGC 6791 cluster is estimated to be ~8 Gyr as mentioned in Sect. 1. Therefore, we integrated the orbits forward in time for 8 Gyr. Finally, the orbital parameters were calculated for individual rotations of each orbit as in Sect. 3.2.
The properties of the bar and spiral arms that we used for the integrations of the NGC 6791’s orbit (as listed in Table 2) are based on recent observations of the MW (see Sect. 2). We indicate this presentday MW model as MW1. From Nbody simulations and observations of external galaxies at various redshifts, we can infer the evolution of bar and spiral arms with time. For example, owing to the expected increase of stellar velocity dispersion with time, the disk becomes more stable and thus its spiral structure weakens. On the other hand, gas inflow can cool down the disk and the structure can strengthen again. Central mass concentrations (e.g., gas inflow) can decrease bar amplitudes (e.g., Athanassoula et al. 2005; Debattista et al. 2006) and even destroy bars, reforming later from the newly formed stars, but with a decreased strength and size (e.g., Bournaud & Combes 2002; Combes 2008).
Therefore, although we do not know the MW bar and spiral structure history, we can safely assume that those were much stronger in the past. For this reason we also followed orbits in a potential including stronger and faster spiral arms, as well as a stronger bar. This model is denoted as MW2. In Table 6 we list the values of the parameters for which MW1 and MW2 differ. We decided to show an example of a model that is deliberately tuned to favor the radial migration caused by stronger bar and spiral arms and also by the choice of the rotating speed of the spiral arms pattern, bringing its resonances closer to the GC. We note that the parameter space (bar and spiral arms shape, extent, strength, and pattern speed) within the framework of our test particle integrations in analytically expressed potentials is very large and its systematic exploration was not the goal of the paper.
Fig. 4 Time evolution of the apogalacticon R_{a} against perigalacticon R_{p} plane for the MW1 model. A mean probability of finding an orbit within a given bin of R_{a} and R_{p} (where the mean is calculated over the indicated 1 Gyr time interval) is mapped. The bin size is 0.1 × 0.1 kpc. The probability for every bin is calculated as a sum of orbits with the bin’s R_{a} and R_{p} values, where the contribution of each orbit is divided by the total number of revolutions per orbit within the time interval (so the total contribution of each orbit to each map is equal to 1). Finally, the probability is normalized to the total number of orbits 10^{4} and converted to %. Each column shows orbits for models including different components – axisymmetric model, model with bar, model with spiral arms, and finally bar and spiral arms combined model. Orbits are mapped within time intervals of 1 Gyr (0–1, 3–4, and 7–8 Gyr) as indicated in the left upper corner of each plot. The NGC 6791’s values of R_{a} and R_{p} with their three sigma errorbars are plotted with a red cross in each plot. The four lines in each of the plots show the lines of constant eccentricity for values 0.0 (circular orbits), 0.25, 0.5, and 0.75, respectively, with decreasing slope. 

Open with DEXTER 
Fig. 5 Time evolution of the apogalacticon R_{a} against perigalacticon R_{p} plane for the MW2 model, see Fig. 4 for a detailed description of the plot. 

Open with DEXTER 
In detail, the bar semimajor axis of the MW2 model was chosen to be longer, of 5 kpc (but the axis ratio for other two axes was kept as for MW1 shape of bar) and the bar itself was more massive – so that on a time scale of 0.4 Gyr it replaced the whole mass of the bulge. As a consequence, the strength parameter of this bar is 0.5. As for spiral arms, the spiral pattern was set to have a higher amplitude and higher pitch angle. The angular velocity of the bar in this case was decreased to 45 km s^{1} kpc^{1}. This brings the bar’s corotation radius to 5.1 kpc and so the ratio of bar’s corotation radius to its semimajor axis is consistent with the value of 0.9–1.3 found in most barred galaxies (see paragraph 6.5.1 in Binney & Tremaine 2008,and references therein). On the other hand, the spiral arms angular velocity is increased to 30 km s^{1} kpc^{1}. The resonances of spiral arms are therefore closer to the GC where the orbits are originally distributed. See Table 7 for resonance radii in both models.
Before starting to comment more on the results, we recall the reader that we aimed to search for similarities between the presentday orbital parameters of NGC 6791 and the orbital parameters of orbits that started close to the GC and then migrated outward. If there were similarities, we could claim that the scenario in which NGC 6791 could form in the very inner disk is plausible also from a dynamical point of view, and not supported only by the observed chemical characteristics of the cluster.
The outcome of these experiments is illustrated in Figs. 4–9. Figures 4 and 5 map the evolution of mean probability of orbit with a given combination of R_{a} and R_{p}. For the direct comparison with the recent orbital parameters of NGC 6791 we also indicate its values with the red cross in the maps. Figures 6 and 7 map the mean probability in the vertical amplitude z_{max} against the R_{a} and R_{p} plane for the MW2 model, respectively. As mentioned above, model MW2 was chosen to obtain an enhanced migration. Comparing Figs. 4 and 5, we can see that this is indeed the case, illustrated by the higher values of apogalacticon seen in Fig. 5. Moreover, some orbits have rotations with similar orbital parameters as we derived for NGC 6791 (the position of the red cross in the R_{a} against R_{p} plane is within the area populated by some of the forward integrated orbits).
From Figs. 6 and 7 we see that only the MW2 combined model can produce orbits that also show a vertical amplitude similar to NGC 6791. Only very few orbits (169 out of 10^{4}) can at least once per the total integration time of 8 Gyr reproduce the R_{a} and R_{p} values similar to NGC 6791, and only some of these orbits (13 out of the 169) meet all three orbital parameters together. To investigate these orbits in more detail, we display all three orbital parameters together in Fig. 9 for the time interval of 7–8 Gyr of the integration.
The experiment with forward integrated orbits illustrated in Figs. 4 to 9 helps us to provide the following considerations:

only model MW2 that includes bar and spiral arms produces the apo, perigalacticons, and eventually also the vertical amplitude similar to the backward integration of the current orbit (see the maps for the combined model in the rightmost column of Figs. 5–7);

Fig. 8 shows an example of an orbit that migrated close to the recent NGC 6791 one. Initially, the orbit is close to circular with the radius of 4 kpc, and later starts to oscillate between very eccentric and close to circular again but has a radius of about 7 kpc;

from our experiment, we estimate the probability for the MW2 model including bar and spiral arms to produce the present orbit of NGC 6791 restricted to 2D Galactic plane to be about 0.40% (during the total integration time of 8 Gyr). This probability is calculated in the same way as for maps in Figs. 4–7: as a sum of orbits reaching R_{a} and R_{p} within 3σ intervals of NGC 6791, where the contribution of each orbit is divided by the total number of revolutions per orbit within the time interval of 1 Gyr (so the total contribution of each orbit within the time interval is equal to 1). Finally, the probability is normalized to the total number of orbits 10^{4} and we sum the probabilities for time intervals of 1 Gyr within the whole integration time of 8 Gyr. If the vertical amplitude is taken into account, the probability has an even lower value of 0.03%.
5. Summary and conclusions
Fig. 6 Time evolution of the apogalacticon R_{a} against vertical amplitude z_{max} plane for the MW2 model. Similar to Fig. 4, results for the models including different components are plotted in columns; time increases with rows. 

Open with DEXTER 
Fig. 7 Time evolution of the perigalacticon R_{p} against vertical amplitude z_{max} plane for the MW2 model. See Fig. 4 for a detailed description. 

Open with DEXTER 
Fig. 8 One of the Galactic orbits and its orbital parameters evolutions calculated in the MW2 model including bar and spiral arms – upper left plot shows the orbit in the Galactic plane; upper right plot shows the timeevolution of the vertical amplitude; the bottom row shows evolutions of eccentricity, and peri and apogalacticon in the left and right plot, respectively. For one of its revolutions, this orbit fits the 3σ intervals for perigalacticon, apogalacticon, and vertical amplitude of the NGC 6791 orbit, these moments are marked with circle symbols in plots on the right. 

Open with DEXTER 
Fig. 9 Apogalacticon R_{a} against perigalacticon R_{p} plane for orbits integrated in the MW2 model during the time interval of 7–8 Gyr. The range of R_{a} aswellas of R_{p} is limited to the 3σ intervals of orbital parameters values of NGC 6791. Each revolution with this R_{a} and R_{p} is shown by a point. The color of each point maps the vertical amplitude z_{max}. Revolutions that match NGC 6791 for all three orbital parameters (within their 3σ intervals) are marked with squared boxes. 

Open with DEXTER 
We investigated a scenario in which the old, metalrich open cluster NGC 6791 formed in the inner disk, close to the bulge, and then moved outward because of perturbations induced by the bar and spiral arms. This scenario might indeed explain its unique properties, since star formation close to the bulge is strong and metal enrichment fast. NGC 6791’s large mass and compactness can explain why it has survived to the present time.
We conducted our investigation using a numerical tool that allowed us to integrate orbits in different Galactic potentials. We calculated the recent (1 Gyr back in time) Galactic orbit of NGC 6791 in various MW models – axisymmetric, including bar and/or spiral structure. The observational uncertainties were taken into account in Monte Carlo fashion, calculating a set of 1000 orbits with their initial conditions distributed according to PM, radial velocity, and distance observational uncertainties. We found that orbits in the axisymmetric model and models including only one rotating pattern differ only slightly. The orbit in the combined bar and spiral arms model is more eccentric with eccentricity up to ~0.3, and has a higher vertical amplitude of ~1 kpc.
To investigate whether these orbital parameters can be achieved as a consequence of migration processes from the inner disk, we also followed a set of forward integrations (for 8 Gyr) with initial positions closer to the GC (with galocentric radii of 3–5 kpc). To this aim, we considered two different realizations of the MW – presentday model: the same we used for the backward integrations (MW1), and models with rotating patterns more supporting the radial migration process (MW2).
The purpose of this analysis was not to reproduce the cluster dynamical history precisely, but to estimate how efficient the migration induced by the MW bar and spiral structure could be and how high the probability of the orbit realization within the limitations of our model is. We indeed found that our MW2 model that incorporates both the bar and spiral arm perturbations, can produce orbits with apogalacticon and perigalacticon similar to the actual values of NGC 6791. However, the probability of this scenario, as quantified from our experiment, is low, approximately 0.40%. Moreover, our migration scenario struggles to reproduce the vertical amplitude of NGC 6791. This might be an artifact of the Flynn et al. (1996) potential, which does not represent the vertical structure of the disk well over the studied range of galocentric distances, as well as of uncertainties in the modeled vertical structure of the bar and spiral arms. We are aware that this could also imply a problem with the migration scenario itself, a question that is unfortunately insolvable within the limitations of the current model.
Our study is unable to provide evidence for the inner origin of NGC 6791 and its later displacement due to radial migration caused by the bar and the spiral arms resonance overlap. This could be viewed as an argument in support of the external (accreted dwarf satellite) origin of NGC 6791, in line with the UV upturn argument mentioned in Sect. 1. On the other hand, our simple model gives a nonzero probability of migration from the inner galaxy, which would therefore deserve additional investigation with improved models of the MW and with a more systematically covering of model’s parameter space.
The recent results by Minchev & Famaey (2010) and Minchev et al. (2011) show that the addition of a second perturbation creates stochastic regions throughout the galactic disk and thus allows for a largescale migration without resorting to the effect of shortlived transients. We note that for transient spirals this effect would be amplified. We expect that the addition of a secondary spiral density wave, as found in Nbody simulations of barred systems (Quillen et al. 2011; Minchev et al., in prep.), that may possibly exist in the MW (Naoz & Shaviv 2007), will make the migration process more efficient. Therefore, the strengths of bar and spirals we used in this study may be the upper limits on what is required to push NGC 6791 outward.
Future studies should strive to incorporate more realistic dynamical models, e.g., considering the time evolution of the bar and spiral structure parameters, the effect of a gaseous component initially present in the disk, and continuous gas accretion.
Parameters for MW models used for the forward integrations.
Radii (in kpc) of the most important resonances for the Galactic bar and spiral arms MW models – corotation R_{CR}, inner Lindblad resonances R_{ILR}, and outer Lindblad resonances R_{OLR}.
We also tested the model of Allen & Santillan (1991), slightly modified to follow the more recent observations rotation curve as described below. The nature of the general results presented in Sects. 3 and 4 stays similar as in the case when the axisymmetric component is modeled as described in Sect. 2.1. In its simplicity the model introduced by Allen & Santillan (1991) does not follow the most updated knowledge of the MW mass distribution. However, it has been recently and extensively used by several teams dealing with orbits of MW clusters (e.g., Allen et al. 2006, 2008; Bellini et al. 2010; Wu et al. 2009) or with the solar neighborhood kinematics (Antoja et al. 2009, 2011).
Acknowledgments
This research has made use of the WEBDA database, operated at the Institute for Astronomy of the University of Vienna (http://www.univie.ac.at/webda/). L.J. acknowledges the support of twoyear ESO Ph.D. studentship, held in ESO, Santiago, as well as by grants No. 205/08/H005 (Czech Science Foundation) and MUNI/A/0968/2009 (Masaryk University in Brno). B.J. was supported by the grant LC06014Center for Theoretical Astrophysics (Czech Ministry of Education) and by the research plan AV0Z10030501 (Academy of Sciences of the Czech Republic). We thank the anonymous referee, whose valuable comments and suggestions helped to significantly improve the paper.
References
 Allen, C., & Santillan, A. 1991, Rev. Mex. Astron. Astrofis., 22, 255 [NASA ADS] [Google Scholar]
 Allen, C., Moreno, E., & Pichardo, B. 2006, ApJ, 652, 1150 [NASA ADS] [CrossRef] [Google Scholar]
 Allen, C., Moreno, E., & Pichardo, B. 2008, ApJ, 674, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Antoja, T., Valenzuela, O., Pichardo, B., et al. 2009, ApJ, 700, L78 [NASA ADS] [CrossRef] [Google Scholar]
 Antoja, T., Figueras, F., RomeroGómez, M., et al. 2011, MNRAS, 418, 1423 [NASA ADS] [CrossRef] [Google Scholar]
 Athanassoula, E., Lambert, J. C., & Dehnen, W. 2005, MNRAS, 363, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Bedin, L. R., Piotto, G., Carraro, G., King, I. R., & Anderson, J. 2006, A&A, 460, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bellini, A., Bedin, L. R., Pichardo, B., et al. 2010, A&A, 513, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bensby, T., Feltzing, S., Johnson, J. A., et al. 2010, A&A, 512, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Binney, J., & Tremaine, S. 2008, Galactic Dynamics, second edition (Princeton University Press) [Google Scholar]
 Binney, J., Gerhard, O., & Spergel, D. 1997, MNRAS, 288, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Bird, J. C., Kazantzidis, S., & Weinberg, D. H. 2012, MNRAS, 420, 913 [NASA ADS] [CrossRef] [Google Scholar]
 Blitz, L., & Spergel, D. N. 1991, ApJ, 379, 631 [NASA ADS] [CrossRef] [Google Scholar]
 Block, D. L., Buta, R., Knapen, J. H., et al. 2004, AJ, 128, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Bournaud, F., & Combes, F. 2002, A&A, 392, 83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brogaard, K., Bruntt, H., Grundahl, F., et al. 2011, A&A, 525, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buson, L. M., Bertone, E., Buzzoni, A., & Carraro, G. 2006, Baltic Astron., 15, 49 [NASA ADS] [Google Scholar]
 Buta, R., & Block, D. L. 2001, ApJ, 550, 243 [NASA ADS] [CrossRef] [Google Scholar]
 Buzzoni, A., Bertone, E., Carraro, G., & Buson, L. 2012, ApJ, 749, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Carraro, G., Villanova, S., Demarque, P., et al. 2006, ApJ, 643, 1151 [NASA ADS] [CrossRef] [Google Scholar]
 Combes, F. 2008 [arXiv:0811.0153] [Google Scholar]
 Combes, F., & Sanders, R. H. 1981, A&A, 96, 164 [NASA ADS] [Google Scholar]
 Coşkunoǧlu, B., Ak, S., Bilir, S., et al. 2011, MNRAS, 412, 1237 [NASA ADS] [Google Scholar]
 Cox, D. P., & Gómez, G. C. 2002, ApJS, 142, 261 [NASA ADS] [CrossRef] [Google Scholar]
 Debattista, V. P., Mayer, L., Carollo, C. M., et al. 2006, ApJ, 645, 209 [NASA ADS] [CrossRef] [Google Scholar]
 Dehnen, W. 2000, AJ, 119, 800 [NASA ADS] [CrossRef] [Google Scholar]
 Dehnen, W., & Binney, J. J. 1998, MNRAS, 298, 387 [NASA ADS] [CrossRef] [Google Scholar]
 Dinescu, D. I., Girard, T. M., & van Altena, W. F. 1999, AJ, 117, 1792 [NASA ADS] [CrossRef] [Google Scholar]
 Drimmel, R. 2000, A&A, 358, L13 [NASA ADS] [Google Scholar]
 Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181 [NASA ADS] [CrossRef] [Google Scholar]
 Dwek, E., Arendt, R. G., Hauser, M. G., et al. 1995, ApJ, 445, 716 [NASA ADS] [CrossRef] [Google Scholar]
 Englmaier, P., & Gerhard, O. 1999, MNRAS, 304, 512 [NASA ADS] [CrossRef] [Google Scholar]
 Flynn, C., SommerLarsen, J., & Christensen, P. R. 1996, MNRAS, 281, 1027 [NASA ADS] [CrossRef] [Google Scholar]
 Fux, R. 1997, A&A, 327, 983 [NASA ADS] [Google Scholar]
 Fux, R. 1999, A&A, 345, 787 [NASA ADS] [Google Scholar]
 Gardner, E., & Flynn, C. 2010, MNRAS, 405, 545 [NASA ADS] [Google Scholar]
 Grosbol, P., Carraro, G., & Beletski, Y. 2011, The Messenger, 03, 47 [NASA ADS] [Google Scholar]
 Janes, K. A. 1984, PASP, 96, 977 [NASA ADS] [CrossRef] [Google Scholar]
 Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [NASA ADS] [CrossRef] [Google Scholar]
 Jurić, M., Ivezić, Ž., Brooks, A., et al. 2008, ApJ, 673, 864 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Kaluzny, J., & Rucinski, S. M. 1995, A&AS, 114, 1 [NASA ADS] [Google Scholar]
 Liebert, J., Saffer, R. A., & Green, E. M. 1994, AJ, 107, 1408 [NASA ADS] [CrossRef] [Google Scholar]
 Milone, A. P., Bedin, L. R., Piotto, G., & Anderson, J. 2009, A&A, 497, 755 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Minchev, I., & Famaey, B. 2010, ApJ, 722, 112 [NASA ADS] [CrossRef] [Google Scholar]
 Minchev, I., & Famaey, B. 2011, EAS Pub. Ser., 45, 299 [CrossRef] [EDP Sciences] [Google Scholar]
 Minchev, I., Nordhaus, J., & Quillen, A. C. 2007, ApJ, 664, L31 [NASA ADS] [CrossRef] [Google Scholar]
 Minchev, I., Boily, C., Siebert, A., & Bienayme, O. 2010, MNRAS, 407, 2122 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Minchev, I., Famaey, B., Combes, F., et al. 2011, A&A, 527, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Miyamoto, M., & Nagai, R. 1975, PASJ, 27, 533 [NASA ADS] [Google Scholar]
 Naoz, S., & Shaviv, N. J. 2007, New A, 12, 410 [NASA ADS] [CrossRef] [Google Scholar]
 Origlia, L., Valenti, E., Rich, R. M., & Ferraro, F. R. 2006, ApJ, 646, 499 [NASA ADS] [CrossRef] [Google Scholar]
 Pfenniger, D. 1984, A&A, 134, 373 [NASA ADS] [Google Scholar]
 Pichardo, B., Martos, M., & Moreno, E. 2004, ApJ, 609, 144 [NASA ADS] [CrossRef] [Google Scholar]
 Piotto, G. 2009, in IAU Symp. 258, ed. E. E. Mamajek, D. R. Soderblom, & R. F. G. Wyse, 233 [Google Scholar]
 Piotto, G., Villanova, S., Bedin, L. R., et al. 2005, ApJ, 621, 777 [NASA ADS] [CrossRef] [Google Scholar]
 Plummer, H. C. 1911, MNRAS, 71, 460 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN, The art of scientific computing (Cambridge: University Press) [Google Scholar]
 Quillen, A. C., Minchev, I., BlandHawthorn, J., & Haywood, M. 2009, MNRAS, 397, 1599 [NASA ADS] [CrossRef] [Google Scholar]
 Quillen, A. C., Dougherty, J., Bagley, M. B., Minchev, I., & Comparetta, J. 2011, MNRAS, 417, 762 [NASA ADS] [CrossRef] [Google Scholar]
 Rattenbury, N. J., Mao, S., Sumi, T., & Smith, M. C. 2007, MNRAS, 378, 1064 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, M. J., & Brunthaler, A. 2004, ApJ, 616, 872 [NASA ADS] [CrossRef] [Google Scholar]
 Reid, M. J., Menten, K. M., Zheng, X. W., et al. 2009, ApJ, 700, 137 [NASA ADS] [CrossRef] [Google Scholar]
 Roškar, R., Debattista, V. P., Quinn, T. R., Stinson, G. S., & Wadsley, J. 2008, ApJ, 684, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Roškar, R., Debattista, V. P., Loebman, S. R., Ivezić, Ž., & Quinn, T. R. 2011 [arXiv:1101.1202] [Google Scholar]
 Rubele, S., Kerber, L., & Girardi, L. 2010, MNRAS, 403, 1156 [NASA ADS] [CrossRef] [Google Scholar]
 Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829 [NASA ADS] [CrossRef] [Google Scholar]
 Sellwood, J. A., & Binney, J. J. 2002, MNRAS, 336, 785 [NASA ADS] [CrossRef] [Google Scholar]
 Twarog, B. A., Carraro, G., & AnthonyTwarog, B. J. 2011, ApJ, 727, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Vallée, J. P. 2005, AJ, 130, 569 [NASA ADS] [CrossRef] [Google Scholar]
 Vallée, J. P. 2008, AJ, 135, 1301 [NASA ADS] [CrossRef] [Google Scholar]
 Van de Putte, D., Garnier, T. P., Ferreras, I., Mignani, R. P., & Cropper, M. 2010, MNRAS, 407, 2109 [NASA ADS] [CrossRef] [Google Scholar]
 Villanova, S., Carraro, G., Scarpa, R., & Marconi, G. 2010, New A, 15, 520 [NASA ADS] [CrossRef] [Google Scholar]
 Wu, Z., Zhou, X., Ma, J., & Du, C. 2009, MNRAS, 399, 2146 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Input observational data used for the calculation of the initial conditions – equatorial coordinates α, δ (J2000), radial velocity v_{r}, PM μ_{α}cosδ, μ_{δ}, and distance to the Sun d_{⊙}; together with their references.
Parameters of initial condition distributions for the forward orbits integration.
Radii (in kpc) of the most important resonances for the Galactic bar and spiral arms MW models – corotation R_{CR}, inner Lindblad resonances R_{ILR}, and outer Lindblad resonances R_{OLR}.
All Figures
Fig. 1 Maps of strength of nonaxisymmetric components in the combined (bar and spiral arms) today’s MW model (noted as MW1 in Sect. 4). Color scales map the ratio between radial (left) and tangential (right) force resulting from the nonmonopole part of the gravitational potential to the radial force due to the monopole part at the same radius. Maps are plotted for the Galactic plane. ⊙ symbol marks the position of the Sun. 

Open with DEXTER  
In the text 
Fig. 2 Recent portion of Galactic orbit of NGC 6791 – orbit projections into the Galactic and meridional plane are plotted in the left and right column, respectively. Orbit in axisymmetric potential is plotted with the blue line. The red line shows the combined model (bar and spiral arms). The initial conditions (see Table 4) are given by mean observational input data (Table 3), the black dot marks the initial position. Integration was backward for time of 1 Gyr. 

Open with DEXTER  
In the text 
Fig. 3 Distributions of cluster’s recent orbital parameters. The distributions for perigalacticon R_{p}, apogalacticon R_{a}, vertical height of orbit z_{max}, and eccentricity e are plotted in columns. Different colors show results for different MW models: axisymmetric model – light blue, model with bar – black, model with spiral arms – yellow, and red for the bar and spiral arms combined model. 

Open with DEXTER  
In the text 
Fig. 4 Time evolution of the apogalacticon R_{a} against perigalacticon R_{p} plane for the MW1 model. A mean probability of finding an orbit within a given bin of R_{a} and R_{p} (where the mean is calculated over the indicated 1 Gyr time interval) is mapped. The bin size is 0.1 × 0.1 kpc. The probability for every bin is calculated as a sum of orbits with the bin’s R_{a} and R_{p} values, where the contribution of each orbit is divided by the total number of revolutions per orbit within the time interval (so the total contribution of each orbit to each map is equal to 1). Finally, the probability is normalized to the total number of orbits 10^{4} and converted to %. Each column shows orbits for models including different components – axisymmetric model, model with bar, model with spiral arms, and finally bar and spiral arms combined model. Orbits are mapped within time intervals of 1 Gyr (0–1, 3–4, and 7–8 Gyr) as indicated in the left upper corner of each plot. The NGC 6791’s values of R_{a} and R_{p} with their three sigma errorbars are plotted with a red cross in each plot. The four lines in each of the plots show the lines of constant eccentricity for values 0.0 (circular orbits), 0.25, 0.5, and 0.75, respectively, with decreasing slope. 

Open with DEXTER  
In the text 
Fig. 5 Time evolution of the apogalacticon R_{a} against perigalacticon R_{p} plane for the MW2 model, see Fig. 4 for a detailed description of the plot. 

Open with DEXTER  
In the text 
Fig. 6 Time evolution of the apogalacticon R_{a} against vertical amplitude z_{max} plane for the MW2 model. Similar to Fig. 4, results for the models including different components are plotted in columns; time increases with rows. 

Open with DEXTER  
In the text 
Fig. 7 Time evolution of the perigalacticon R_{p} against vertical amplitude z_{max} plane for the MW2 model. See Fig. 4 for a detailed description. 

Open with DEXTER  
In the text 
Fig. 8 One of the Galactic orbits and its orbital parameters evolutions calculated in the MW2 model including bar and spiral arms – upper left plot shows the orbit in the Galactic plane; upper right plot shows the timeevolution of the vertical amplitude; the bottom row shows evolutions of eccentricity, and peri and apogalacticon in the left and right plot, respectively. For one of its revolutions, this orbit fits the 3σ intervals for perigalacticon, apogalacticon, and vertical amplitude of the NGC 6791 orbit, these moments are marked with circle symbols in plots on the right. 

Open with DEXTER  
In the text 
Fig. 9 Apogalacticon R_{a} against perigalacticon R_{p} plane for orbits integrated in the MW2 model during the time interval of 7–8 Gyr. The range of R_{a} aswellas of R_{p} is limited to the 3σ intervals of orbital parameters values of NGC 6791. Each revolution with this R_{a} and R_{p} is shown by a point. The color of each point maps the vertical amplitude z_{max}. Revolutions that match NGC 6791 for all three orbital parameters (within their 3σ intervals) are marked with squared boxes. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.