The present paper is the culminating one of a series aimed to contribute to the understanding of the kinematic structures of the solar neighbourhood (SN), explaining the origin of the Local Arm and relating the moving groups with the spiral-arms resonances in the disk. With a model for the Galactic potential, with the Sun inside the spiral corotation resonance (CR), we integrate the 2D orbits of test particles distributed in birthplaces along the main spiral arms, the Local Arm, and in the axisymmetric disk. A comparison of the resulting U-V plane of the SN with that provided by Gaia DR2 confirms our previous conclusion that the moving groups of Coma Berenices, Pleiades, and Hyades are associated with the CR, and that the Hercules stream is formed by the bulk of high-order inner Lindblad resonances. The kinematic structures result from stellar orbits trapped by the spiral resonances in a timespan of ~ 1 Gyr, indicating the long-living nature of the spiral structure and challenging recent arguments in favor of short-lived structures originated from incomplete phase mixing in the Galactic disk. As a byproduct, our simulations give some insight into the birthplaces of the stars presently located in the SN; the majority of stars of the main moving groups and the Hercules stream were likely born in the Local Arm, while stars of the Sirius group possibly originated from the outer segment of the Sagittarius-Carina arm. We also propose the spiral resonances as the dynamical origin for the diagonal ridges in the Galactic distribution of rotation velocities.

In this work we revisit the issue of the rotation speed of the spiral arms and the location of the corotation radius of our Galaxy. This research was performed using homogeneous data set of young open clusters (age < 50 Myr) determined from Gaia DR2 data. The stellar astrometric membership were determined using proper motions and parallaxes, taking into account the full covariance matrix. The distance, age, reddening and metallicity of the clusters were determined by our non subjective multidimensional global optimization tool to fit theoretical isochrones to Gaia DR2 photometric data. The rotation speed of the arms is obtained from the relation between age and angular distance of the birthplace of the clusters to the present-day position of the arms. Using the clusters belonging to the Sagittarius-Carina, Local and Perseus arms, and adopting the Galactic parameters $R_0$ = 8.3 kpc and $V_0$ = 240 km\,s$^{-1}$, we determine a pattern speed of $28.2 \pm 2.1$ km\,s$^{-1}$\,kpc$^{-1}$, with no difference between the arms. This implies that the corotation radius is $R_c = 8.51 \pm 0.64$ kpc, close to the solar Galactic orbit ($R_c/R_0 = 1.02\pm0.07$).

Using the Gaia data release 2 (DR2), we analyzed the distribution of stars in the close vicinity of the Sun in the full 3D position-velocity space. We have found no evidence of incomplete phase mixing in the vertical direction of the disk, which could be originated by some external events. We show that the vertical phase space spiral $Z$-$V_z$ is produced by the well-known moving groups (MGs), mainly by Coma-Berenices, Pleiades-Hyades and Sirius, when the statistical characteristics (mean, median, or mode) of the azimuthal velocity $V_\varphi$ are used to analyze the distribution in the vertical position-velocity plane. This result does not invoke external perturbations and is independent on the internal dynamical mechanisms that originate the MGs. Our conclusions counterbalance current arguments in favor of short-lived (between 300 and 900 Myr) structures in the solar neighborhood. Contrarily, they support the hypothesis of a longer formation time scale (around a few Gyr) for the MGs.

The aim of this work is to contribute to the understanding of the stellar velocity distribution in the solar neighborhood (SN). We propose that the structures on the $U-V$ planes, known as the moving groups, can be mainly explained by the spiral arms perturbations. The applied model of the Galactic disk and spiral arms, with the parameters defined by observational data and with pattern speed $\Omega_p=$28.0 km s$^{-1}$ kpc$^{-1}$, is the same that allowed us to explain the origin of the Local Arm and the Sun's orbit trapped inside the corotation resonance (CR). We show that the $U-V$ picture of the SN consists of the main component, associated with the CR, and the inner and outer structures, which we could associate with the Hercules and Sirius streams, respectively. The Coma-Berenices and Hyades-Pleiades groups and the Sun itself belong to the main part. The substructures of Hercules are formed mainly by the nearby 8/1, 12/1, and even 6/1 inner Lindblad resonances, while Sirius is shaped by the bulk of overlapping outer Lindblad resonances, -8/1, -12/1, -16/1, which are stuck to the CR. This richness in resonances only exists near corotation, which should be of the spiral arms, not of the Galactic bar, whose stable corotation zone is far away from the Sun. The model's predictions of the velocity distribution match qualitatively and quantitatively the distribution provided by Gaia DR2.

Observational data indicate that the Milky Way is a barred spiral galaxy. Computation facilities and availability of data from Galactic surveys stimulate the appearance of models of the Galactic structure. More efforts to build dynamical models containing both spiral arms and the central bar/bulge are needed. We expand the study of the stellar dynamics in the Galaxy by adding the bar/bulge component to a model with spiral arms introduced in our previous paper. The model is tested by applying it to the solar neighborhood, where observational data are more precise. We model analytically the potential of the Galaxy to derive the force field in its equatorial plane. The model comprises an axisymmetric disc derived from the observable rotation curve, four arms with Gaussian-shaped groove profiles, and a classical elongated/oblate ellipsoidal bar/bulge structure. The parameters describing the bar/bulge are constrained by observations and the stellar dynamics, and their possible limits are determined. A basic model results in a bar of 2.9 kpc in length, with a mass of the order of a few 10$^9M_\odot$. The size and orientation of the bar are also restricted by the position of masers with VLBI distances. The bar's rotation speed is constrained to $\Omega_{\rm bar}<50$ km s$^{-1}$ kpc$^{-1}$ taking into account the allowed mass range. We conclude that our basic model is compatible with observations and with the dynamical constraints. The model explains simultaneously the bulk of the main moving groups, associated here with the spiral corotation resonance, and the Hercules stream, associated with several inner high-order spiral resonances; in particular, with the 8/1 resonance. From the dynamical constraints on the bar's angular speed, it is unlikely that the bar's OLR lies near the solar circle; moreover, its proximity would compromise the stability of the Local Arm structure.

The Local arm of the Milky Way, a short spiral feature near the Sun whose existence is known for decades, was recently observed in detail with different tracers. Many efforts have been dedicated to elaborate plausible hypotheses concerning the origin of the main spiral arms of the Galaxy; however, up to now, no specific mechanism for the origin of the Local arm was proposed. Here we explain, for the first time, the Local arm as an outcome of the spiral corotation resonance, which traps arm tracers and the Sun inside it. We show that the majority of maser sources belonging to the Local arm, together with the Sun, evolve inside the corotation resonance, never crossing the main spiral arms but instead oscillating in the region between them. This peculiar behavior of the Sun could have numerous consequences to our understanding of the local kinematics of stars, the Galactic Habitable Zone, and the Solar System evolution.

Context: Resonances in the stellar orbital motion under perturbations from spiral arms structure play an important role in the evolution of the disks of spiral galaxies. The epicyclic approximation allows the determination of the corresponding resonant radii on the equatorial plane (for nearly circular orbits), but is not suitable in general. Aims: We expand the study of resonant orbits by analysing stellar motions perturbed by spiral arms with Gaussian-shaped profiles without any restriction on the stellar orbital configurations, and we expand the concept of Lindblad (epicyclic) resonances for orbits with large radial excursions. Methods: We define a representative plane of initial conditions, which covers the whole phase space of the system. Dynamical maps on representative planes are constructed numerically, in order to characterize the phase-space structure and identify the precise location of resonances. The study is complemented by the construction of dynamical power spectra, which provide the identification of fundamental oscillatory patterns in the stellar motion. Results: Our approach allows a precise description of the resonance chains in the whole phase space, giving a broader view of the dynamics of the system when compared to the classical epicyclic approach, even for objects in retrograde motion. The analysis of the solar neighbourhood shows that, depending on the current azimuthal phase of the Sun with respect to the spiral arms, a star with solar kinematic parameters may evolve either inside the stable co-rotation resonance or in a chaotic zone. Conclusions: Our approach contributes to quantifying the domains of resonant orbits and the degree of chaos in the whole Galactic phase-space structure. It may serve as a starting point to apply these techniques to the investigation of clumps in the distribution of stars in the Galaxy, such as kinematic moving groups.

Aims. Galaxy mass models based on simple and analytical functions for the density and potential pairs have been widely proposed in the literature. Disk models constrained by kinematic data alone give information on the global disk structure only very near the Galactic plane. We attempt to circumvent this issue by constructing disk mass models whose three-dimensional structures are constrained by a recent Galactic star counts model in the near-infrared and also by observations of the hydrogen distribution in the disk. Our main aim is to provide models for the gravitational potential of the Galaxy that are fully analytical but also with a more realistic description of the density distribution in the disk component. Methods. From the disk model directly based on the observations (here divided into the thin and thick stellar disks and the HI and H$_2$ disks subcomponents), we produce fitted mass models by combining three Miyamoto-Nagai disk profiles of any "model order" (1, 2, or 3) for each disk subcomponent. The Miyamoto-Nagai disks are combined with models for the bulge and "dark halo" components and the total set of parameters is adjusted by observational kinematic constraints. A model which includes a ring density structure in the disk, beyond the solar Galactic radius, is also investigated. Results. The Galactic mass models return very good matches to the imposed observational constraints. In particular, the model with the ring density structure provides a greater contribution of the disk to the rotational support inside the solar circle. The gravitational potential models and their associated force-fields are described in analytically closed forms, and in addition, they are also compatible with our best knowledge of the stellar and gas distributions in the disk component. The gravitational potential models are suited for investigations of orbits in the Galactic disk.

We analyse the secular effects of a long-lived Galactic spiral structure on the stellar orbits with mean radii close to the corotation resonance. By test-particle simulations and different spiral potential models with parameters constrained on observations, we verified the formation of a minimum with amplitude ~ 30% - 40% of the background disk stellar density at corotation. Such minimum is formed by the secular angular momentum transfer between stars and the spiral density wave on both sides of corotation. We demonstrate that the secular loss (gain) of angular momentum and decrease (increase) of mean orbital radius of stars just inside (outside) corotation can counterbalance the opposite trend of exchange of angular momentum shown by stars orbiting the librational points L_4/5 at the corotation circle. Such secular processes actually allow steady spiral waves to promote radial migration across corotation. We propose some observational evidences for the minimum stellar density in the Galactic disk, such as its direct relation to the minimum in the observed rotation curve of the Galaxy at the radius R ~ 9 kpc (for R_0 = 7.5 kpc), as well as its association with a minimum in the distribution of Galactic radii of a sample of open clusters older than 1 Gyr. The closeness of the solar orbit radius to the corotation resonance implies that the solar orbit lies inside a ring of minimum surface density (stellar + gas). This also implies in a correction to larger values for the estimated total mass of the Galactic disk, and consequently, a greater contribution of the disk component to the inner rotation curve of the Galaxy.

In order to understand the Barium abundance distribution in the Galactic disk based on Cepheids, one must first be aware of important effects of the corotation resonance, situated a little beyond the solar orbit. The thin disk of the Galaxy is divided in two regions that are separated by a barrier situated at that radius. Since the gas cannot get across that barrier, the chemical evolution is independent on the two sides of it. The barrier is caused by the opposite directions of flows of gas, on the two sides, in addition to a Cassini-like ring void of HI (caused itself by the flows). A step in the metallicity gradient developed at corotation, due to the difference in the average star formation rate on the two sides, and to this lack of communication between them. In connection with this, a proof that the spiral arms of our Galaxy are long-lived (a few billion years) is the existence of this step. When one studies the abundance gradients by means of stars which span a range of ages, like the Cepheids, one has to take into account that stars, contrary to the gas, have the possibility of crossing the corotation barrier. A few stars born on the high metallicity side are seen on the low metallicity one, and vice-versa. In the present work we re-discuss the data on Barium abundance in Cepheids as a function of Galactic radius, taking into account the scenario described above. The [Ba/H] ratio, plotted as a function of Galactic radius, apparently presents a distribution with two branches in the external region (beyond corotation). One can re-interpret the data and attribute the upper branch to the stars that were born on the high metallicity side. The lower branch, analyzed separately, indicates that the stars born beyond corotation have a rising Barium metallicity as a function of Galactic radius.

We propose a new, more realistic, description of the perturbed gravitational potential of spiral galaxies, with spiral arms having Gaussian-shaped groove profiles. We investigate the stable stellar orbits in galactic disks, using the new perturbed potential. The influence of the bulge mass on the stellar orbits in the inner regions of a disk is also investigated. The new description offers the advantage of easy control of the parameters of the Gaussian profile of its potential. We find a range of values for the perturbation amplitude from 400 to 800 km^2 s^{-2} kpc^{-1} which implies a maximum ratio of the tangential force to the axisymmetric force between 3% and 6%, approximately. Good self-consistency of arm shapes is obtained between the Inner Lindblad resonance (ILR) and the 4:1 resonance. Near the 4:1 resonance the response density starts to deviate from the imposed logarithmic spiral form. This creates bifurcations that appear as short arms. Therefore the deviation from a perfect logarithmic spiral in galaxies can be understood as a natural effect of the 4:1 resonance. Beyond the 4:1 resonance we find closed orbits which have similarities with the arms observed in our Galaxy. In regions near the center, in the presence of a massive bulge, elongated stellar orbits appear naturally, without imposing any bar-shaped potential, but only extending the spiral perturbation a little inward of the ILR. This suggests that a bar is formed with a half-size around 3 kpc by a mechanism similar to that of the spiral arms. The potential energy perturbation that we adopted represents an important step in the direction of self-consistency, compared to previous sine function descriptions of the potential. Our model produces a realistic description of the spiral structure, able to explain several details that were not yet understood.

The connection between some features of the metallicity gradient in the Galactic disk, best revealed by Open Clusters and Cepheids, and the spiral structure, is explored. The step-like abrupt decrease in metallicity at 8.5 kpc (with R_0= 7.5 kpc, or at 9.5 kpc if R_0 = 8.5 kpc is adopted) is well explained by the corotation ring-shaped gap in the density of gas, which isolates the internal and external regions of the disk one from the other. This solves a long standing problem of understanding the different chemical characteristics of the inner and outer parts of the disk. The time required to build up the metallicity difference between the two sides of the step is a measure of the minimal life-time of the present grand-design spiral pattern structure, of the order of 3 Gyr. The plateaux observed on each side of the step are interpreted in terms of the large scale radial motion of the stars and of the gas flow induced by the spiral structure. The star-formation rate revealed by the density of open clusters is maximum in the Galactic radial range from 6 to 12 kpc (with an exception of a narrow gap at corotation), coinciding with the region where the 4-arms mode is allowed to exist. We argue that most of the old open clusters situated at large galactocentric radii were born in this inner region where conditions more favorable to star-formation are found. The ratio of $\alpha$-elements to Fe of the sample of Cepheids does not vary appreciably with the Galactic radius, which reveals an homogeneous history of star formation. Different arguments are given showing that usual approximations of chemical evolution models, which assume fast mixing of metallicity in the azimuthal direction and ignore the existence of the spiral arms, are a poor ones.