The large astrometric and photometric survey performed by the Gaia mission allows for a panoptic view of the Galactic disc and in its stellar cluster population. Hundreds of clusters were only discovered after the latest G data release (DR2) and have yet to be characterised. Here we make use of the deep and homogeneous Gaia photometry down to G=18 to estimate the distance, age, and interstellar reddening for about 2000 clusters identified with Gaia~DR2 astrometry. We use these objects to study the structure and evolution of the Galactic disc. We rely on a set of objects with well-determined parameters in the literature to train an artificial neural network to estimate parameters from the Gaia photometry of cluster members and their mean parallax. We obtain reliable parameters for 1867 clusters. Our new homogeneous catalogue confirms the relative lack of old clusters in the inner disc (with a few notable exceptions). We also quantify and discuss the variation of scale height with cluster age, and detect the Galactic warp in the distribution of old clusters. This work results in a large and homogenous cluster catalogue. However, the present sample is still unable to trace the Outer spiral arm of the Milky Way, which indicates that the outer disc cluster census might still be incomplete.

The evolution of the Milky Way disk, which contains most of the stars in the Galaxy, is affected by several phenomena. For example, the bar and the spiral arms of the Milky Way induce radial migration of stars and can trap or scatter stars close to orbital resonances. External perturbations from satellite galaxies can also have a role, causing dynamical heating of the Galaxy, ring-like structures in the disk and correlations between different components of the stellar velocity. These perturbations can also cause 'phase wrapping' signatures in the disk, such as arched velocity structures in the motions of stars in the Galactic plane. Some manifestations of these dynamical processes have already been detected, including kinematic substructure in samples of nearby stars, density asymmetries and velocities across the Galactic disk that differ from the axisymmetric and equilibrium expectations, especially in the vertical direction, and signatures of incomplete phase mixing in the disk. Here we report an analysis of the motions of six million stars in the Milky Way disk. We show that the phase-space distribution contains different substructures with various morphologies, such as snail shells and ridges, when spatial and velocity coordinates are combined. We infer that the disk must have been perturbed between 300 million and 900 million years ago, consistent with estimates of the previous pericentric passage of the Sagittarius dwarf galaxy. Our findings show that the Galactic disk is dynamically young and that modelling it as time-independent and axisymmetric is incorrect.

We study the dynamics near the unstable Lagrangian points in galactic bar models using dynamical system tools in order to determine the global morphology of a barred galaxy. We aim at the case of non-autonomous models, in particular with secular evolution, by allowing the bar pattern speed to decrease with time. We extend the concept of manifolds widely used in the autonomous problem to the Lagrangian Coherent Structures (LCS), widely used in fluid dynamics, which behave similar to the invariant manifolds driving the motion. After adapting the LCS computation code to the galactic dynamics problem, we apply it to both the autonomous and non-autonomous problems, relating the results with the manifolds and identifying the objects that best describe the motion in the non-autonomous case. We see that the strainlines coincide with the first intersection of the stable manifold when applied to the autonomous case, while, when the secular model is used, the strainlines still show the regions of maximal repulsion associated to both the corresponding stable manifolds and regions with a steep change of energy. The global morphology of the galaxy predicted by the autonomous problem remains unchanged.

Most galaxies have a warped shape when they are seen from an edge-on point of view. The reason for this curious form is not completely known so far and in this work we apply dynamical system tools to contribute to its explanation. Starting from a simple, but realistic, model formed by a bar and a disc, we study the effect produced by a small misalignment between the angular momentum of the system and its angular velocity. To this end, a precession model is developed and considered, assuming that the bar behaves like a rigid body. After checking that the periodic orbits inside the bar keep being the skeleton of the inner system, even after inflicting a precession to the potential, we compute the invariant manifolds of the unstable periodic orbits departing from the equilibrium points at the ends of the bar to get evidences of their warped shapes. As it is well known, the invariant manifolds associated with these periodic orbits drive the arms and rings of barred galaxies and constitute the skeleton of these building blocks. Looking at them from a side-on viewpoint, we find that these manifolds present warped shapes as those recognized in observations. Lastly, test particle simulations have been performed to determine how the stars are affected by the applied precession, confirming this way the theoretical results obtained.

In this first paper we simulate the population of disc Red Clump stars to be observed by Gaia. We generate a set of test particles and we evolve it in a 3D barred Milky Way like galactic potential. We assign physical properties of the Red Clump trace population and a realistic 3D interstellar extinction model. We add Gaia observational constraints and an error model according to the pre-commissioning scientific performance assessments. We present and analyse two mock catalogues, offered to the community, that are an excellent test bed for testing tools being developed for the future scientific exploitation of Gaia data. The first catalogue contains stars up to Gaia G 20, while the second is the subset containing Gaia radial velocity data with a maximum error of sigmaVr=10 kms. Here we present first attempts to characterise the density structure of the Galactic bar in the Gaia space of observables. The Gaia large errors in parallax and the high interstellar extinction in the inner parts of the Galactic disc prevent us to model the bar overdensity. This result suggests the need to combine Gaia and IR data to undertake such studies. We find that IR photometric distances for this Gaia sample allow us to recover the Galactic bar orientation angle with an accuracy of approximately 5 degrees.

We explore the possibility of detecting and characterizing the warp of the stellar disc of our Galaxy using synthetic Gaia data. The availability of proper motions and, for the brightest stars radial velocities, adds a new dimension to this study. A family of Great Circle Cell Counts (GC3) methods is used. They are ideally suited to find the tilt and twist of a collection of rings, which allow us to detect and measure the warp parameters. To test them, we use random realizations of test particles which evolve in a realistic Galactic potential warped adiabatically to various final configurations. In some cases a twist is introduced additionally. The Gaia selection function, its errors model and a realistic 3D extinction map are applied to mimic three tracer populations: OB, A and Red Clump stars. We show how the use of kinematics improves the accuracy in the recovery of the warp parameters. The OB stars are demonstrated to be the best tracers determining the tilt angle with accuracy better than $\sim0.5$ up to Galactocentric distance of $\sim16$ kpc. Using data with good astrometric quality, the same accuracy is obtained for A type stars up to $\sim13$ kpc and for Red Clump up to the expected stellar cut-off. Using OB stars the twist angle is recovered to within $< 3^\circ$ for all distances.

We study the bar-driven dynamics in the inner part of the Milky Way by using invariant manifolds. This theory has been successfully applied to describe the morphology and kinematics of rings and spirals in external galaxies, and now, for the first time, we apply it to the Milky Way. We compute the orbits confined by the invariant manifolds of the unstable periodic orbits located at the ends of the bar. We discuss whether the COBE/DIRBE bar and the Long bar compose a single bar or two independent bars and perform a number of comparisons which, taken together, argue strongly in favour of the former. More specifically, we favour the possibility that the so-called COBE/DIRBE bar is the boxy/peanut bulge of a bar whose outer thin parts are the so-called Long bar. This possibility is in good agreement both with observations of external galaxies, with orbital structure theory and with simulations. We then analyse in detail the morphology and kinematics given by five representative Galactic potentials. Two have a Ferrers bar, two have a quadrupole bar and the last one a composite bar. We first consider only the COBE/DIRBE bar and then extend it to include the effect of the Long bar. We find that the large-scale structure given by the manifolds describes an inner ring, whose size is similar to the near and far 3-kpc arm, and an outer ring, whose properties resemble those of the Galactic Molecular Ring. We also analyse the kinematics of these two structures, under the different galactic potentials, and find they reproduce the main over-densities found in the galactic longitude-velocity CO diagram. Finally, we consider for what model parameters, the global morphology of the manifolds may reproduce the two outer spiral arms. We conclude that this would necessitate either more massive and more rapidly rotating bars, or including in the potential an extra component describing the spiral arms.

In a series of papers, we propose a theory to explain the formation and properties of rings and spirals in barred galaxies. The building blocks of these structures are orbits guided by the manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In this paper, the last of the series, we present more comparisons of our theoretical results to observations and also give new predictions for further comparisons. Our theory provides the right building blocks for the rectangular-like bar outline and for ansae. We consider how our results can be used to give estimates for the pattern speed values, as well as their effect on abundance gradients in barred galaxies. We present the kinematics along the manifold loci, to allow comparisons with the observed kinematics along the ring and spiral loci. We consider gaseous arms and their relations to stellar ones. We discuss several theoretical aspects and stress that the orbits that constitute the building blocks of the spirals and rings are chaotic. They are, nevertheless, spatially well confined by the manifolds and are thus able to outline the relevant structures. Such chaos can be termed `confined chaos' and can play a very important role in understanding the formation and evolution of galaxy structures and in galactic dynamics in general. This work, in agreement with several others, argues convincingly that galactic dynamic studies should not be limited to the study of regular motions and orbits.

In this series of papers, we propose a theory to explain the formation and properties of rings and spirals in barred galaxies. The building blocks of these structures are orbits guided by the manifolds emanating from the unstable Lagrangian points located near the ends of the bar. In this paper we focus on a comparison of the morphology of observed and of theoretical spirals and rings and we also give some predictions for further comparisons. Our theory can account for spirals as well as both inner and outer rings. The model outer rings have the observed $R_1$, $R_1'$, $R_2$, $R_2'$ and $R_1R_2$ morphologies, including the dimples near the direction of the bar major axis. We explain why the vast majority of spirals in barred galaxies are two armed and trailing, and discuss what it would take for higher multiplicity arms to form. We show that the shapes of observed and theoretical spirals agree and we predict that stronger non-axisymmetric forcings at and somewhat beyond corotation will drive more open spirals. We compare the ratio of ring diameters in theory and in observations and predict that more elliptical rings will correspond to stronger forcings. We find that the model potential may influence strongly the numerical values of these ratios.

Two-armed, grand design spirals and inner and outer rings in barred galaxies can be due to orbits guided by the manifolds emanating from the vicinity of the L1 and L2 Lagrangian points, located at the ends of the bar. We first summarise the necessary theoretical background and in particular we describe the dynamics around the unstable equilibrium points in barred galaxy models, and the corresponding homoclinic and heteroclinic orbits. We then discuss two specific morphologies and the circulation of material within the corresponding manifolds. We also discuss the case where mass concentrations at the end of the bar can stabilise the L1 and L2 and the relevance of this work to the gas concentrations in spirals and rings.

In this paper we present building blocks which can explain the formation and properties both of spirals and of inner and outer rings in barred galaxies. We first briefly summarise the main results of the full theoretical description we have given elsewhere, presenting them in a more physical way, aimed to an understanding without the requirement of extended knowledge of dynamical systems or of orbital structure. We introduce in this manner the notion of manifolds, which can be thought of as tubes guiding the orbits. The dynamics of these manifolds can govern the properties of spirals and of inner and outer rings in barred galaxies. We find that the bar strength affects how unstable the L1 and L2 Lagrangian points are, the motion within the 5A5A5Amanifold tubes and the time necessary for particles in a manifold to make a complete turn around the galactic centre. We also show that the strength of the bar, or, to be more precise, of the non-axisymmetric forcing at and somewhat beyond the corotation region, determines the resulting morphology. Thus, less strong bars give rise to R1 rings or pseudorings, while stronger bars drive R2, R1R2 and spiral morphologies. We examine the morphology as a function of the main parameters of the bar and present descriptive two dimensional plots to that avail. We also derive how the manifold morphologies and properties are modified if the L1 and L2 Lagrangian points become stable. Finally, we discuss how dissipation affects the manifold properties and compare the manifolds in gas-like and in stellar cases. Comparison with observations, as well as clear predictions to be tested by observations will be given in an accompanying paper.

We study the role of the unstable equilibrium points in the transfer of matter in a galaxy using the potential of a rotating triaxial system. In particular, we study the neighbourhood of these points for energy levels and for main model parameters where the zero velocity curves just open and form a bottleneck in the region. For these energies, the transfer of matter from the inner to the outer parts and vice versa starts being possible. We study how the dynamics around the unstable equilibrium points is driven, by performing a partial normal form scheme and by computing the invariant manifolds of periodic orbits and quasi-periodic orbits using the reduced Hamiltonian. In particular, we compute some homoclinic and heteroclinic orbits playing a crucial role. Our results also show that in slow rotating and/or axisymmetric systems the hyperbolic character of the equilibrium points is cancelled, so that no transfer of matter is possible through the bottleneck.

We propose a theory to explain the formation of spiral arms and of all types of outer rings in barred galaxies, extending and applying the technique used in celestial mechanics to compute transfer orbits. Thus, our theory is based on the chaotic orbital motion driven by the invariant manifolds associated to the periodic orbits around the hyperbolic equilibrium points. In particular, spiral arms and outer rings are related to the presence of heteroclinic or homoclinic orbits. Thus, R1 rings are associated to the presence of heteroclinic orbits, while R1R2 rings are associated to the presence of homoclinic orbits. Spiral arms and R2 rings, however, appear when there exist neither heteroclinic nor homoclinic orbits. We examine the parameter space of three realistic, yet simple, barred galaxy models and discuss the formation of the different morphologies according to the properties of the galaxy model. The different morphologies arise from differences in the dynamical parameters of the galaxy.

We propose a theory to explain the formation of both spirals and rings in barred galaxies using a common dynamical framework. It is based on the orbital motion driven by the unstable equilibrium points of the rotating bar potential. Thus, spirals, rings and pseudo-rings are related to the invariant manifolds associated to the periodic orbits around these equilibrium points. We examine the parameter space of three barred galaxy models and discuss the formation of the different morphological structures according to the properties of the bar model. We also study the influence of the shape of the rotation curve in the outer parts, by making families of models with rising, flat or falling rotation curves in the outer parts. The differences between spiral and ringed structures arise from differences in the dynamical parameters of the host galaxies.

We propose a new theory to explain the formation of spiral arms and of all types of outer rings in barred galaxies. We have extended and applied the technique used in celestial mechanics to compute transfer orbits. Thus, our theory is based on the chaotic orbital motion driven by the invariant manifolds associated to the periodic orbits around the hyperbolic equilibrium points. In particular, spiral arms and outer rings are related to the presence of heteroclinic or homoclinic orbits. Thus, R1 rings are associated to the presence of heteroclinic orbits, while R1R2 rings are associated to the presence of homoclinic orbits. Spiral arms and R2 rings, however, appear when there exist neither heteroclinic nor homoclinic orbits. We examine the parameter space of three realistic, yet simple, barred galaxy models and discuss the formation of the different morphologies according to the properties of the galaxy model. The different morphologies arise from differences in the dynamical parameters of the galaxy.

In this and in a previous paper (Romero-Gomez et al. 2006) we propose a theory to explain the formation of both spirals and rings in barred galaxies using a common dynamical framework. It is based on the orbital motion driven by the unstable equilibrium points of the rotating bar potential. Thus, spirals, rings and pseudo-rings are related to the invariant manifolds associated to the periodic orbits around these equilibrium points. We examine the parameter space of three barred galaxy models and discuss the formation of the different morphological structures according to the properties of the bar model. We also study the influence of the shape of the rotation curve in the outer parts, by making families of models with rising, flat, or falling rotation curves in the outer parts. The differences between spiral and ringed structures arise from differences in the dynamical parameters of the host galaxies. The results presented here will be discussed and compared with observations in a forthcoming paper.

We propose a new theory for the formation of rR_1 ring structures, i.e. for ring structures with both an inner and an outer ring, the latter having the form of ``8''. We propose that these rings are formed by material from the stable and unstable invariant manifolds associated with the Lyapunov orbits around the equilibrium points of a barred galaxy. We discuss the shape and velocity structure of the rings thus formed and argue that they are in agreement with the observed properties of rR_1 structures.