In a cosmological setting, the disc of a galaxy is expected to continuously experience gravitational torques and perturbations from a variety of sources, which can cause the disc to wobble, flare and warp. Specifically, the study of galactic warps and their dynamical nature can potentially reveal key information on the formation history of galaxies and the mass distribution of their halos. Our Milky Way presents a unique case study for galactic warps, thanks to the detailed knowledge of its stellar distribution and kinematics. Using a simple model of how the warp's orientation is changing with time, we here measure the precession rate of the Milky Way's warp using 12 million giant stars from Gaia Data Release 2, finding that it is precessing at $10.86 \pm 0.03_{stat} \pm 3.20_{sys}$ km/s/kpc in the direction of Galactic rotation, about one third the angular velocity at the Sun's position in the Galaxy. The direction and magnitude of the warp's precession rate favour the scenario that the warp is the results of a recent or ongoing encounter with a satellite galaxy, rather than the relic of the ancient assembly history of the Galaxy.

We present rotation period measurements for 107 M dwarfs in the mass range $0.15-0.70 M_\odot$ observed within the context of the APACHE photometric survey. We measure rotation periods in the range 0.5-190 days, with the distribution peaking at $\sim$ 30 days. We revise the stellar masses and radii for our sample of rotators by exploiting the Gaia DR2 data. For $\sim 20\%$ of the sample, we compare the photometric rotation periods with those derived from different spectroscopic indicators, finding good correspondence in most cases. We compare our rotation periods distribution to the one obtained by the Kepler survey in the same mass range, and to that derived by the MEarth survey for stars in the mass range $0.07-0.25 M_\odot$. The APACHE and Kepler periods distributions are in good agreement, confirming the reliability of our results, while the APACHE distribution is consistent with the MEarth result only for the older/slow rotators, and in the overlapping mass range of the two surveys. Combining the APACHE/Kepler distribution with the MEarth distribution, we highlight that the rotation period increases with decreasing stellar mass, in agreement with previous work. Our findings also suggest that the spin-down time scale, from fast to slow rotators, changes crossing the fully convective limit at $\approx0.3 M_\odot$ for M dwarfs. The catalogue of 107 rotating M dwarfs presented here is particularly timely, as the stars are prime targets for the potential identification of transiting small planets with TESS and amenable to high-precision mass determination and further atmospheric characterization measurements.

Flat rotation curves in disk galaxies represent the main evidence for large amounts of surrounding dark matter. Despite of the difficulty in identifying the dark matter contribution to the total mass density in our Galaxy, stellar kinematics, as tracer of gravitational potential, is the most reliable observable for gauging different matter components. This work tests the flatness of the MW rotation curve with a simple general relativistic model suitable to represent the geometry of a disk as a stationary axisymmetric dust metric at a sufficiently large distance from a central body. Circular velocities of unprecedented accuracy were derived from the Gaia DR2 data for a carefully selected sample of disk stars. We then fit these velocities to both the classical, i.e. including a dark matter halo, rotation curve model and a relativistic analogue, as derived form the solution of Einstein's equation. The GR-compliant MW rotational curve model results statistically indistinguishable from its state-of-the-art DM analogue. This supports our ansatz that a stationary and axisymmetric galaxy-scale metric could "fill the gap" in a baryons-only Milky Way, suggestive of star orbits dragged along the background geometry. We confirmed that geometry is a manifestation of gravity according to the Einstein theory, in particular the weak gravitational effect due to the off-diagonal term of the metric could mimic for a "DM-like" effect in the observed flatness of the MW rotation curve. In the context of Local Cosmology, our findings are suggestive of a Galaxy phase-space as the exterior gravitational field of a Kerr-like source (inner rotating bulge) without the need of extra-matter.

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.

Using Gaia DR2 astrometry, we map the kinematic signature of the Galactic stellar warp out to a distance of 7 kpc from the Sun. Combining Gaia DR2 and 2MASS photometry, we identify, via a probabilistic approach, 599 494 upper main sequence stars and 12 616 068 giants without the need for individual extinction estimates. The spatial distribution of the upper main sequence stars clearly shows segments of the nearest spiral arms. The large-scale kinematics of both the upper main sequence and giant populations show a clear signature of the warp of the Milky Way, apparent as a gradient of 5-6 km/s in the vertical velocities from 8 to 14 kpc in Galactic radius. The presence of the signal in both samples, which have different typical ages, suggests that the warp is a gravitationally induced phenomenon.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $\sim$5~kpc to $\sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)

The nature and origin of the Galactic warp represent one of the open questions posed by Galactic evolution. Thanks to Gaia high precision absolute astrometry, steps towards the understanding of the warp's dynamical nature can be made. Indeed, proper motions for long-lived stable warp are expected to show measurable trends in the component vertical to the galactic plane. Within this context, we search for the kinematic warp signal in the first Gaia data release (DR1). By analyzing distant spectroscopically-identified OB stars in the Hipparcos subset in Gaia DR1, we find that the kinematic trends cannot be explained by a simple model of a long-lived warp. We therefore discuss possible scenarios for the interpretation of the obtained results. We also present current work in progress to select a larger sample of OB star candidates from the Tycho-Gaia Astrometric Solution (TGAS) subsample in DR1, and delineate the points that we will be addressing in the near future.

The mechanism responsible for the warp of our Galaxy, as well as its dynamical nature, continues to remain unknown. With the advent of high precision astrometry, new horizons have been opened for detecting the kinematics associated with the warp and constraining possible warp formation scenarios for the Milky Way. The aim of this contribution is to establish whether the first Gaia data release (DR1) shows significant evidence of the kinematic signature expected from a long-lived Galactic warp in the kinematics of distant OB stars. As the first paper in a series, we present our approach for analyzing the proper motions and apply it to the sub-sample of Hipparcos stars. We select a sample of 989 distant spectroscopically-identified OB stars from the New Reduction of Hipparcos, of which 758 are also in Gaia DR1, covering distances from 0.5 to 3 kpc from the Sun. We develop a model of the spatial distribution and kinematics of the OB stars from which we produce the probability distribution functions of the proper motions, with and without the systematic motions expected from a long-lived warp. A likelihood analysis is used to compare the expectations of the models with the observed proper motions from both Hipparcos and Gaia DR1. We find that the proper motions of the nearby OB stars are consistent with the signature of a warp, while those of the more distant stars (parallax<1 mas) are not. The kinematics of our sample of young OB stars suggests that systematic vertical motions in the disk cannot be explained by a simple model of a stable long-lived warp. The Galactic warp may either be a transient feature, or additional phenomena are acting on the gaseous component of the Milky Way, causing systematic vertical motions that are masking the expected warp signal. A larger and deeper sample of stars with Gaia astrometry will be needed to constrain the dynamical nature of the Galactic warp.