Very few examples are known of red supergiant runaways, all of them descending from the more massive O-type precursors, but none from the lower mass B-type precursors, although runaway statistics among B-type stars suggest that K-type runaways must be relatively numerous. We study HD 137071, a star that has been considered so far as a normal K-type red giant. Its parallax measured by Gaia and the derived luminosity suggest that it is actually a supergiant, whereas its derived distance to the galactic plane and its spatial velocity of 54.1 km s$^{-1}$ with respect to the local standard of rest suggest that it is also a runaway star. However, intrinsic limitations in determining the trigonometric parallaxes of cool supergiants, even in the Gaia era, require accurate spectral classifications for confirmation. We reliably classify HD 137071 as a K4II star establishing its membership to the extreme Population I, which is in agreement with the luminosity derived using the Gaia DR2 parallax measurement. Kinematical data from the Gaia DR2 catalog confirm its high spatial velocity and its runaway nature. Combining the spectral classification with astrometric information, a state-of-the-art galactic potential model, and evolutionary models for high-mass stars we trace the motion of HD 137071 back to the proximities of the galactic plane and speculate on which of the two proposed mechanisms for the production of runaway stars may be responsible for the high velocity of HD 137071. The available data favor the formation of HD 137071 in a massive binary system where the more massive companion underwent a supernova explosion about 32 Myr ago.

There are few warp kinematic models of the Galaxy able to characterise structure and kinematics. These models are necessary to study the lopsidedness of the warp and the twisting of the line-of-nodes of the stellar warp, already seen in gas and dust. We use the \Gaia~Data Release 2 astrometric data up to $G=20$mag to characterise the structure of the Galactic warp, the vertical motions and the dependency on the age. We use two populations up to galactocentric distances of $16$kpc, a young (OB-type) and an old (Red Giant Branch, RGB). We use the nGC3 PCM and LonKin methods based on the Gaia observables, together with 2D projections of the positions and proper motions in the Galactic plane. We confirm the age dependency of the Galactic warp, both in positions and kinematics, being the height of the Galactic warp of about $0.2$kpc for the OB sample and of $1.$kpc for the RGB at a galactocentric distance of $14$kpc. Both methods find that the onset radius is $12\sim 13$kpc for the OB sample and $10\sim 11$kpc for the RGB. From the RGB sample, we find from galactocentric distances larger than $10$kpc the line-of-nodes twists away from the Sun-anticentre line towards galactic azimuths $\sim 180-200^{\circ}$ increasing with radius, though possibly influenced by extinction. The RGB sample reveals a slightly lopsided stellar warp with $\sim 250$pc between the up and down sides. The line of maximum of proper motions in latitude is systematically offset from the line-of-nodes estimated from the spatial data, which our models predict as a kinematic signature of lopsidedness. We also show a prominent wave-like pattern of a bending mode different in the OB and RGB, and substructures that might not be related to the Galactic warp nor to a bending mode. GDR2 triggers the need for complex kinematic models, flexible enough to combine both wave-like patterns and an S-shaped lopsided warp.[abridged]

We use Gaia DR2 magnitudes, colours and parallaxes for stars with G<12 to explore a 15-dimensional space that includes simultaneously the initial mass function (IMF) and a non-parametric star formation history (SFH) for the Galactic disc. This inference is performed by combining the Besancon Galaxy Model fast approximate simulations (BGM FASt) and an approximate Bayesian computation algorithm. We find in Gaia DR2 data an imprint of a star formation burst 2-3 Gyr ago, in the Galactic thin disc domain, and a present star formation rate (SFR) of about 1 Msun. Our results show a decreasing trend of the SFR from 9-10 Gyr to 6-7 Gyr ago. This is consistent with the cosmological star formation quenching observed at redshifts z<1.8. This decreasing trend is followed by a SFR enhancement starting at about 5 Gyr ago and continuing until about 1 Gyr ago which is detected with high statistical significance by discarding the null hypothesis of an exponential SFH with a p-value=0.002. We estimate, from our best fit model, that about the 50% of the mass used to generate stars, along the thin disc life, was expended in the period from 5 to 1 Gyr ago. The timescale and the amount of stellar mass generated during the SFR enhancement event lead us to hypothesise that its origin, currently under investigation, is not intrinsic to the disc. Thus, an external perturbation is needed for its explanation. Additionally we find, for the thin disc, a slope of the IMF of $\alpha_3 \approx 2$ for masses M>1.53 Msun and $\alpha_2 \approx 1.3$ for the mass range between 0,5 and 1,53 Msun. This is the first time that we consider a non-parametric SFH for the thin disc in the Besancon Galaxy Model. This new step, together with the capabilities of the Gaia DR2 parallaxes to break degeneracies between different stellar populations, allow us to better constrain the SFH and the IMF.

Context. Open clusters are very good tracers of the evolution of the Galactic disc. Thanks to Gaia, their kinematics can be investigated with an unprecedented precision and accuracy. Aims. The distribution of open clusters in the 6D phase space is revisited with Gaia DR2. Methods. The weighted mean radial velocity of open clusters was determined, using the most probable members available from a previous astrometric investigation that also provided mean parallaxes and proper motions. Those parameters, all derived from Gaia DR2 only, were combined to provide the 6D phase space information of 861 clusters. The velocity distribution of nearby clusters was investigated, as well as the spatial and velocity distributions of the whole sample as a function of age. A high quality subsample was used to investigate some possible pairs and groups of clusters sharing the same Galactic position and velocity. Results. For the high quality sample that has 406 clusters, the median uncertainty of the weighted mean radial velocity is 0.5 km/s. The accuracy, assessed by comparison to ground-based high resolution spectroscopy, is better than 1 km/s. Open clusters nicely follow the velocity distribution of field stars in the close Solar neighbourhood previously revealed by Gaia DR2. As expected, the vertical distribution of young clusters is very flat but the novelty is the high precision to which this can be seen. The dispersion of vertical velocities of young clusters is at the level of 5 km/s. Clusters older than 1 Gyr span distances to the Galactic plane up to 1 kpc with a vertical velocity dispersion of 14 km/s, typical of the thin disc. Five pairs of clusters and one group with five members are possibly physically related. Other binary candidates previously identified turn out to be chance alignment.

SN 1572 (Tycho Brahe's supernova) clearly belongs to the Ia (thermonuclear) type. It was produced by the explosion of a white dwarf in a binary system. Its remnant has been the first of this type to be explored in search of a possible surviving companion, the mass donor that brought the white dwarf to the point of explosion. A high peculiar motion with respect to the stars at the same location in the Galaxy, mainly due to the orbital velocity at the time of the explosion, is a basic criterion for the detection of such companions. Radial velocities from the spectra of the stars close to the geometrical center of Tycho's supernova remnant, plus proper motions of the same stars, obtained by astrometry with the {\it Hubble Space Telescope}, have been used so far. In addition, a detailed chemical analysis of the atmospheres of a sample of candidate stars had been made. However, the distances to the stars, remained uncertain. Now, the Second {\it Gaia} Data Release (DR2) provides unprecedent accurate distances and new proper motions for the stars can be compared with those made from the {\it HST}. We consider the Galactic orbits that the candidate stars to SN companion would have in the future. We do this to explore kinematic peculiarity. We also locate a representative sample of candidate stars in the Toomre diagram. Using the new data, we reevaluate here the status of the candidates suggested thus far, as well as the larger sample of the stars seen in the central region of the remnant.

Context: The Initial Mass Function (IMF) plays a crucial role on galaxy evolution and its implications on star formation theory make it a milestone for the next decade. It is in the intermediate and high mass ranges where the uncertainties of the IMF are larger. This is a major subject of debate and analysis both for Galactic and extragalactic science. Aims: Our goal is to constrain the IMF of the Galactic thin disc population using both Galactic Classical Cepheids and Tycho-2 data. Methods: For the first time the Besan\c{c}on Galaxy Model (BGM) has been used to characterise the Galactic population of the Classical Cepheids. We have modified the age configuration in the youngest populations of the BGM thin disc model to avoid artificial discontinuities in the age distribution of the simulated Cepheids. Three statistical methods, optimized for different mass ranges, have been developed and applied to search for the best IMF that fits the observations. This strategy allows us to quantify variations in the Star Formation History (SFH), the stellar density at Sun position and the thin disc radial scale length. A rigorous treatment of unresolved multiple stellar systems has been undertaken adopting a spatial resolution according to the catalogues used. Results: For intermediate masses, our study favours a composite field-star IMF slope of $\alpha=3.2$ for the local thin disc, excluding flatter values such as the Salpeter IMF ($\alpha=2.35$). Moreover, a constant Star Formation History is definitively excluded, the three statistical methods considered here show that it is inconsistent with the observational data. Conclusions: Using field stars and Galactic Classical Cepheids, we have found, above $1M_\odot$, an IMF steeper than the canonical stellar IMF of associations and young clusters. This result is consistent with the predictions of the Integrated Galactic IMF.

We develop a new theoretical framework to generate Besan\c{c}on Galaxy Model fast approximate simulations (BGM FASt) to address fundamental questions of the Galactic structure and evolution performing multi-parameter inference. As a first application of our strategy we simultaneously infer the IMF, the star formation history and the stellar mass density in the Solar Neighbourhood. The BGM FASt strategy is based on a reweighing scheme, that uses a specific pre-sampled simulation, and on the assumption that the distribution function of the generated stars in the Galaxy can be described by an analytical expression. To validate BGM FASt we execute a set of tests. Finally, we use BGM FASt with an approximate Bayesian computation algorithm to obtain the posterior PDF of the inferred parameters, by comparing synthetic versus Tycho-2 colour-magnitude diagrams. Results: The validation shows a very good agreement between BGM FASt and the standard BGM, with BGM FASt being $\approx 10^4$ times faster. By analysing Tycho-2 data we obtain a thin disc star formation history decreasing in time and a present rate of $1.2 \pm 0.2 M_\odot/yr$. The resulting total stellar mass density in the Solar Neighbourhood is $0.051_{-0.005}^{+0.002} M_\odot/pc^3$ and the local dark matter density is $0.012 \pm 0.001 M_\odot/pc^3$. For the composite IMF we obtain a slope of $\alpha_2={2.1}_{-0.3}^{+0.1}$ in the mass range between $0.5 M_\odot$ and $1.53M_\odot$. The results of the slope at the high mass range are trustable up to $4M_\odot$ and highly depend on the choice of the extinction map (obtaining $\alpha_3={2.9}_{-0.2}^{+0.2}$ and $\alpha_3={3.7}_{-0.2}^{+0.2}$ respectively, for two different extinction maps). Systematic uncertainties are not included. Conclusions: The good performance of BGM FASt demonstrates that it is a very valuable tool to perform multi-parameter inference using Gaia data releases.

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 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.

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.

Gaia DR2 has delivered full-sky 6-D measurements for millions of stars, and the quest to understand the dynamics of our Galaxy has entered a new phase. Our aim is to reveal and characterize the kinematic sub-structure of the different Galactic neighbourhoods, to form a picture of their spatial evolution that can be used to infer the Galactic potential, its evolution and its components. We take ~5 million stars in the Galactic disk from the Gaia DR2 catalogue and build the velocity distribution of many different Galactic Neighbourhoods distributed along 5 kpc in Galactic radius and azimuth. We decompose their distribution of stars in the V_R-V_phi plane with the wavelet transformation and asses the statistical significance of the structures found. We detect many kinematic sub-structures (arches and more rounded groups) that diminish their azimuthal velocity as a function of Galactic radius in a continuous way, connecting volumes up to 3 kpc apart in some cases. The decrease rate is, on average, of ~23 km/s/kpc. In azimuth, the kinematic sub-structures present much smaller variations. We also observe a duality in their behaviour: some conserve their vertical angular momentum with radius (e.g., Hercules), while some seem to have nearly constant kinetic energy (e.g., Sirius). These two trends are consistent with the approximate predictions of resonances and of phase mixing, respectively. Besides, the overall spatial evolution of Hercules is consistent with being related to the Outer Lindblad Resonance of the Bar. We also detect structures without apparent counterpart in the vicinity of the Sun. The various trends observed and their continuity with radius and azimuth allows for future work to deeply explore the parameter space of the models. Also, the characterization of extrasolar moving groups opens the opportunity to expand our understanding of the Galaxy beyond the Solar Neighbourhood.

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 Young Local Associations constitute an excellent sample for the study of a variety of astrophysical topics, especially the star formation process in low-density environments. Data from the Gaia mission allows us to undertake studies of the YLAs with unprecedented accuracy. We determine the dynamical age and place of birth of a set of associations in a uniform and dynamically consistent manner. There are nine YLAs in our sample $\epsilon$ Chamaeleontis, TW Hydrae, $\beta$ Pictoris, Octans, Tucana-Horologium, Columba, Carina, Argus and AB Doradus. We designed a method for deriving the dynamical age of the YLAs based on the orbital integration. It involves a strategy to account for the effect of observational errors and we tested it using mock YLAs. Finally, we applied it to our set of nine YLAs with astrometry from the first Gaia data release and complementary on-ground radial velocities from the literature. Our orbital analysis yields a first estimate of the dynamical age of 3$^{+9}_{-0}$ Myr, 13$^{+7}_{-0}$ Myr and 5$^{+23}_{-0}$ Myr for $\epsilon$ Cha, $\beta$ Pict and Tuc-Hor, respectively. For four other associations (Oct, Col, Car and Arg), we provide a lower limit for the dynamical age. Our rigorous error treatment indicates that TW Hya and AB Dor deserve further study. Our dynamical ages are compatible spectroscopic and isochrone fitting ages obtained elsewhere. From the orbital analysis, we suggest a scenario with two episodes of star formation: one ~40 Myr ago in the first quadrant that gave birth to $\epsilon$ Cha, TW Hya and $\beta$ Pic, and another 5-15 Myr ago close to the Sun that formed Tuc-Hor, Col, and Car. Future Gaia data will provide the necessary accuracy to improve the present results, especially for the controversial age determinations, and additional evidence for the proposed scenario, once a complete census of YLAs and better membership can be obtained.

We use the Gaia-TGAS data to compare the transverse velocities in Galactic longitude (coming from proper motions and parallaxes) in the Milky Way disk for negative and positive longitudes as a function of distance. The transverse velocities are strongly asymmetric and deviate significantly from the expectations for an axisymmetric Galaxy. The value and sign of the asymmetry changes at spatial scales of several tens of degrees in Galactic longitude and about 0.5 kpc in distance. The asymmetry is statistically significant at 95% confidence level for 57% of the region probed, which extends up to ~1.2 kpc. A percentage of 24% of the region studied shows absolute differences at this confidence level larger than 5 km/s and 7% larger than 10 km/s. The asymmetry pattern shows mild variations in the vertical direction and with stellar type. A first qualitative comparison with spiral arm models indicates that the arms are unlikely to be the main source of the asymmetry. We briefly discuss alternative origins. This is the first time that global all-sky asymmetries are detected in the Milky Way kinematics, beyond the local neighbourhood, and with a purely astrometric sample.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

With the start of the Gaia era, the time has come to address the major challenge of deriving the star formation history and evolution of the disk of our MilkyWay. Here we review our present knowledge of the outer regions of the Milky Way disk population. Its stellar content, its structure and its dynamical and chemical evolution are summarized, focussing on our lack of understanding both from an observational and a theoretical viewpoint. We describe the unprecedented data that Gaia and the upcoming ground-based spectroscopic surveys will provide in the next decade. More in detail, we quantify the expect accuracy in position, velocity and astrophysical parameters of some of the key tracers of the stellar populations in the outer Galactic disk. Some insights on the future capability of these surveys to answer crucial and fundamental issues are discussed, such as the mechanisms driving the spiral arms and the warp formation. Our Galaxy, theMilkyWay, is our cosmological laboratory for understanding the process of formation and evolution of disk galaxies. What we learn in the next decades will be naturally transferred to the extragalactic domain.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

NGC 7067 is a young open cluster located in the direction between the first and the second Galactic quadrants and close to the Perseus spiral arm. This makes it useful for studies of the nature of the Milky Way spiral arms. Stromgren photometry taken with the Wide Field Camera at the Isaac Newton Telescope allowed us to compute individual physical parameters for the observed stars and hence to derive cluster's physical parameters. Spectra from the 1.93-m telescope at the Observatoire de Haute-Provence helped to check and improve the results. We obtained photometry for 1233 stars, individual physical parameters for 515 and spectra for 9 of them. The 139 selected cluster members lead to a cluster distance of 4.4+/-0.4 kpc, with an age below log10(t(yr))=7.3 and a present Mass of 1260+/-160Msun. The morphology of the data reveals that the centre of the cluster is at (ra,dec)=(21:24:13.69,+48:00:39.2) J2000, with a radius of 6.1arcsec. Stromgren and spectroscopic data allowed us to improve the previous parameters available for the cluster in the literature.

Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame.

We present the results of the short-term constancy monitoring of candidate Gaia Spectrophotometric Standard Stars (SPSS). We obtained time series of typically 1.24 hour - with sampling periods from 1-3 min to a few hours, depending on the case - to monitor the constancy of our candidate SPSS down to 10 mmag, as required for the calibration of Gaia photometric data. We monitored 162 out of a total of 212 SPSS candidates. The observing campaign started in 2006 and finished in 2015, using 143 observing nights on nine different instruments covering both hemispheres. Using differential photometry techniques, we built light curves with a typical precision of 4 mmag, depending on the data quality. As a result of our constancy assessment, 150 SPSS candidates were validated against short term variability, and only 12 were rejected because of variability including some widely used flux standards such as BD+174708, SA 105-448, 1740346, and HD 37725.

We introduce a new set of simulations of Milky Way-sized galaxies using the AMR code ART + hydrodynamics in a $\Lambda$CDM cosmogony. The simulation series is named GARROTXA and follow the formation of a halo/galaxy from z~$=$~60 to z~$=$~0. The final virial mass of the system is $\sim$7.4$\times$10$^{11}$M$_{\odot}$. Our results are as follows: (a) contrary to many previous studies, the circular velocity curve shows no central peak and overall agrees with recent MW observations. (b) Other quantities, such as M$\_{*}$(6$\times$10$^{10}$M$_\odot$) and R$_d$ (2.56 kpc), fall well inside the observational MW range. (c) We measure the disk-to-total ratio kinematically and find that D/T=0.42. (d) The cold gas fraction and star formation rate (SFR) at z=0, on the other hand, fall short from the values estimated for the Milky Way. As a first scientific exploitation of the simulation series, we study the spatial distribution of the hot X-ray luminous gas. We have found that most of this X-ray emitting gas is in a halo-like distribution accounting for an important fraction but not all of the missing baryons. An important amount of hot gas is also present in filaments. In all our models there is not a massive disk-like hot gas distribution dominating the column density. Our analysis of hot gas mock observations reveals that the homogeneity assumption leads to an overestimation of the total mass by factors 3 to 5 or to an underestimation by factors $0.7-0.1$, depending on the used observational method. Finally, we confirm a clear correlation between the total hot gas mass and the dark matter halo mass of galactic systems.

We present a method to identify Ultra Faint Dwarf Galaxy (UFDG) candidates in the halo of the Milky Way using the future Gaia catalogue and we explore its detection limits and completeness. The method is based on the Wavelet Transform and searches for over-densities in the combined space of sky coordinates and proper motions, using kinematics in the search for the first time. We test the method with a Gaia mock catalogue that has the Gaia Universe Model Snapshot (GUMS) as a background, and use a library of around 30 000 UFDGs simulated as Plummer spheres with a single stellar population. For the UFDGs we use a wide range of structural and orbital parameters that go beyond the range spanned by real systems, where some UFDGs may remain undetected. We characterize the detection limits as function of the number of observable stars by Gaia in the UFDGs with respect to that of the background and their apparent sizes in the sky and proper motion planes. We find that the addition of proper motions in the search improves considerably the detections compared to a photometric survey at the same magnitude limit. Our experiments suggest that Gaia will be able to detect UFDGs that are similar to some of the known UFDGs even if the limit of Gaia is around 2 magnitudes brighter than that of SDSS, with the advantage of having a full-sky catalogue. We also see that Gaia could even find some UFDGs that have lower surface brightness than the SDSS limit.

The Gaia SpectroPhotometric Standard Stars (SPSS) survey started in 2006, it was awarded almost 450 observing nights, and accumulated almost 100,000 raw data frames, with both photometric and spectroscopic observations. Such large observational effort requires careful, homogeneous, and automated data reduction and quality control procedures. In this paper, we quantitatively evaluate instrumental effects that might have a significant (i.e.,$\geq$1%) impact on the Gaia SPSS flux calibration. The measurements involve six different instruments, monitored over the eight years of observations dedicated to the Gaia flux standards campaigns: DOLORES@TNG in La Palma, EFOSC2@NTT and ROSS@REM in La Silla, CAFOS@2.2m in Calar Alto, BFOSC@Cassini in Loiano, and LaRuca@1.5m in San Pedro Martir. We examine and quantitatively evaluate the following effects: CCD linearity and shutter times, calibration frames stability, lamp flexures, second order contamination, light polarization, and fringing. We present methods to correct for the relevant effects, which can be applied to a wide range of observational projects at similar instruments.

Aims. The main goal of this study is to detect the stellar overdensity associated with the Perseus arm in the anticenter direction. Methods. We used the physical parameters derived from Str\"omgren photometric data to compute the surface density distribution as a function of galactocentric distance for different samples of intermediate young stars. The radial distribution of the interstellar absorption has also been derived. Results. We detected the Perseus arm stellar overdensity at 1.6+-0.2 kpc from the Sun with a significance of 4-5{\sigma} and a surface density amplitude of around 10%, slightly depending on the sample used. Values for the radial scale length of the Galactic disk have been simultaneously fitted obtaining values in the range [2.9,3.5] kpc for the population of the B4-A1 stars. Moreover, the interstellar visual absorption distribution is congruent with a dust layer in front of the Perseus arm. Conclusions. This is the first time that the presence of the Perseus arm stellar overdensity has been detected through individual star counts, and its location matches a variation in the dust distribution. The offset between the dust lane and the overdensity indicates that the Perseus arm is placed inside the co-rotation radius of the Milky Way spiral pattern.

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 derive the vertical velocities of disk stars in the range of Galactocentric radii of R=5-16 kpc within 2 kpc in height from the Galactic plane. This kinematic information is connected to dynamical aspects in the formation and evolution of the Milky Way, such as the passage of satellites and vertical resonance and determines whether the warp is a long-lived or a transient feature. We used the proper motions of the PPMXL survey, correcting of systematic errors with the reference of quasars. From the color-magnitude diagram K versus (J-K) we selected the standard candles corresponding to red clump giants and used the information of their proper motions to build a map of the vertical motions of our Galaxy. We derived the kinematics of the warp both analytically and through a particle simulation to fit these data. Complementarily, we also carried out the same analysis with red clump giants spectroscopically selected with APOGEE data, and we predict the improvements in accuracy that will be reached with future Gaia data. A simple model of warp with the height of the disk z_w(R,phi)=gamma (R-R_sun) sin(phi-phi_w) fits the vertical motions if d(gamma)/dt/gamma=-34+/-17 Gyr^{-1}; the contribution to d(gamma)/dt comes from the southern warp and is negligible in the north. The vertical motion in the warp apparently indicates that the main S-shaped structure of the warp is a long-lived feature, whereas the perturbation that produces an irregularity in the southern part is most likely a transient phenomenon. With the use of the Gaia end-of-mission products together with spectroscopically classified red clump giants, the precision in vertical motions can be increased by an order of magnitude at least.

Aims. The aim is to derive accurate stellar distances and extinctions for young stars of our survey in the Galactic anticenter direction using the Str\"omgren photometric system. This will allow a detailed mapping of the stellar density and absorption toward the Perseus arm. Methods. We developed a new method for deriving physical parameters from Str\"omgren photometry and also implemented and tested it. This is a model-based method that uses the most recent available stellar atmospheric models and evolutionary tracks to interpolate in a 3D grid of the unreddened indexes [m1], [c1] and Hbeta. Distances derived from both this method and the classical pre-Hipparcos calibrations were tested against Hipparcos parallaxes and found to be accurate. Results. Systematic trends in stellar photometric distances derived from empirical calibrations were detected and quantified. Furthermore, a shift in the atmospheric grids in the range Teff=[7000,9000]K was detected and a correction is proposed. The two methods were used to compute distances and reddening for about 12000 OBA-type stars in our Str\"omgren anticenter survey. Data from the IPHAS and 2MASS catalogs were used to complement the detection of emission line stars and to break the degeneracy between early and late photometric regions. We note that photometric distances can differ by more than 20%, those derived from the empirical calibrations being smaller than those derived with the new method, which agree better with the Hipparcos data.

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 map the kinematics of stars in simulated galaxy disks with spiral arms using the velocity ellipsoid vertex deviation (l$_v$). We use test particle simulations, and for the first time, fully self-consistent high resolution N-body models. We compare our maps with the Tight Winding Approximation model analytical predictions. We see that for all barred models spiral arms rotate closely to a rigid body manner and the vertex deviation values correlate with the density peaks position bounded by overdense and underdense regions. In such cases, vertex deviation sign changes from negative to positive when crossing the spiral arms in the direction of disk rotation, in regions where the spiral arms are in between corotation (CR) and the Outer Lindblad Resonance (OLR). By contrast, when the arm sections are inside the CR and outside the OLR, l$_v$ changes from negative to positive.We propose that measurements of the vertex deviations pattern can be used to trace the position of the main resonances of the spiral arms. We propose that this technique might exploit future data from Gaia and APOGEE surveys. For unbarred N-body simulations with spiral arms corotating with disk material at all radii, our analysis suggests that no clear correlation exists between l$_v$ and density structures.

The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allows us to test hypotheses on the star formation history, star evolution, and chemical and dynamical evolution of the Galaxy. The new version of the Besancon Galaxy model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the initial mass function (IMF) and star formation rate (SFR) of the Galactic disc. We present a new strategy for the generation of thin disc stars, which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. The local dynamical self-consistency is maintained in this new scheme. We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. For the first time, a full sky comparison is performed between BGM and data. This strategy allows us to constrain the IMF slope at high masses, which is found to be close to 3.0 and excludes a shallower slope such as Salpeter's one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 Mo/pc^3 implying that there is no compelling evidence for the significant amount of dark matter in the disc. While the model is fitted to Tycho-2 data, which is a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

High resolution N-body simulations using different codes and initial condition techniques reveal two different behaviours for the rotation frequency of transient spiral arms like structures. Whereas unbarred disks present spiral arms nearly corotatingwith disk particles, strong barred models (bulged or bulge-less) quickly develop a bar-spiral structure dominant in density, with a pattern speed almost constant in radius. As the bar strength decreases the arm departs from bar rigid rotation and behaves similar to the unbarred case. In strong barred models we detect in the frequency space other subdominant and slower modes at large radii, in agreement with previous studies, however we also detect them in the configuration space. We propose that the distinctive behaviour of the dominant spiral modes can be exploited in order to constraint the nature of Galactic spiral arms by the astrometric survey GAIA and by 2-D spectroscopic surveys like CALIFA and MANGA in external galaxies.

Aims. The main purpose is to map the radial variation of the stellar space density for the young stellar population in the Galactic anticenter direction in order to understand the structure and location of the Perseus spiral arm. Methods. A uvbyHbeta Stromgren photometric survey covering 16sqrdeg in the anticenter direction was carried out using the Wide Field Camera at the Isaac Newton Telescope. This is the natural photometric system for identifying young stars and obtaining accurate estimates of individual distances and ages. The calibration to the standard system was undertaken using open clusters. Results. We present a main catalog of 35974 stars with all Stromgren indexes and a more extended one with 96980 stars with partial data. The central 8sqrdeg have a limiting magnitude of V<17m, while the outer region reaches V<15.5m. These large samples will permit us to analyze the stellar surface density variation associated to the Perseus arm also to study the properties of the stellar component and the interstellar extinction in the anticenter direction.

We describe two ground based observing campaigns aimed at building a grid of approximately 200 spectrophotometric standard stars (SPSS), with an internal ~1% precision and tied to Vega within ~3%, for the absolute flux calibration of data gathered by Gaia, the ESA astrometric mission. The criteria for the selection and a list of candidates are presented, together with a description of the survey strategy and the adopted data analysis methods. We also discuss a short list of notable rejected SPSS candidates and difficult cases, based on identification problems, literature discordant data, visual companions, and variability. In fact, all candidates are also monitored for constancy (within \pm5 mmag, approximately). In particular, we report on a CALSPEC standard, 1740346, that we found to be a delta Scuti variable during our short-term monitoring (1-2 h) campaign.

With the advent of the Gaia data, the unprecedented kinematic census of great part of the Milky Way disc will allow us to characterise the local kinematic groups and new groups in different disc neighbourhoods. First, we show here that the models predict a stellar kinematic response to the spiral arms and bar strongly dependent on disc position. For example, we find that the kinematic groups induced by the spiral arm models change significantly if one moves only ~ 0.6 kpc in galactocentric radius, but ~ 2 kpc in azimuth. There are more and stronger groups as one approaches the spiral arms. Depending on the spiral pattern speed, the kinematic imprints are more intense in nearby vicinities or far from the Sun. Secondly, we present a preliminary study of the kinematic groups observed by RAVE. This sample will allow us, for the first time, to study the dependence on Galactic position of the (thin and thick) disc moving groups. In the solar neighbourhood, we find the same kinematics groups as detected in previous surveys, but now with better statistics and over a larger spatial volume around the Sun. This indicates that these structures are indeed large scale kinematic features.

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.

We study the spiral arm influence on the solar neighbourhood stellar kinematics. As the nature of the Milky Way (MW) spiral arms is not completely determined, we study two models: the Tight-Winding Approximation (TWA) model, which represents a local approximation, and a model with self-consistent material arms named PERLAS. This is a mass distribution with more abrupt gravitational forces. We perform test particle simulations after tuning the two models to the observational range for the MW spiral arm properties. We explore the effects of the arm properties and find that a significant region of the allowed parameter space favours the appearance of kinematic groups. The velocity distribution is mostly sensitive to the relative spiral arm phase and pattern speed. In all cases the arms induce strong kinematic imprints for pattern speeds around 17 km/s/kpc (close to the 4:1 inner resonance) but no substructure is induced close to corotation. The groups change significantly if one moves only ~0.6 kpc in galactocentric radius, but ~2 kpc in azimuth. The appearance time of each group is different, ranging from 0 to more than 1 Gyr. Recent spiral arms can produce strong kinematic structures. The stellar response to the two potential models is significantly different near the Sun, both in density and kinematics. The PERLAS model triggers more substructure for a larger range of pattern speed values. The kinematic groups can be used to reduce the current uncertainty about the MW spiral structure and to test whether this follows the TWA. However, groups such as the observed ones in the solar vicinity can be reproduced by different parameter combinations. Data from velocity distributions at larger distances are needed for a definitive constraint.

We study the phase space available to the local stellar distribution using a Galactic potential consistent with several recent observational constraints. We find that the induced phase space structure has several observable consequences. The spiral arm contribution to the kinematic structure in the solar neighborhood may be as important as the one produced by the Galactic bar. We suggest that some of the stellar kinematic groups in the solar neighborhood, like the Hercules structure and the kinematic branches, can be created by the dynamical resonances of self-gravitating spiral arms and not exclusively by the Galactic bar. A structure coincident with the Arcturus kinematic group is developed when a hot stellar disk population is considered, which introduces a new perspective on the interpretation of its extragalactic origin. A bar-related resonant mechanism can modify this kinematic structure.We show that particles in the dark matter disk-like structure predicted by recent LCDM galaxy formation experiments, with similar kinematics to the thick disk, are affected by the same resonances, developing phase space structures or dark kinematic groups that are independent of the Galaxy assembly history and substructure abundance. We discuss the possibility of using the stellar phase space groups as constraints to non-axisymmetric models of the Milky Way structure.

Over the last decade, several groups of young (mainly low-mass) stars have been discovered in the solar neighbourhood (closer than ~100 pc), thanks to cross-correlation between X-ray, optical spectroscopy and kinematic data. These young local associations offer insights into the star formation process in low-density environments, shed light on the substellar domain, and could have played an important role in the recent history of the local interstellar medium. Ages estimates for these associations have been derived in the literature by several ways. In this work we have studied the kinematic evolution of young local associations and their relation to other young stellar groups and structures in the local interstellar medium, thus casting new light on recent star formation processes in the solar neighbourhood. We compiled the data published in the literature for young local associations, including the astrometric data from the new Hipparcos reduction. Using a realistic Galactic potential we integrated the orbits for these associations and the Sco-Cen complex back in time. Combining these data with the spatial structure of the Local Bubble and the spiral structure of the Galaxy, we propose a recent history of star formation in the solar neighbourhood. We suggest that both the Sco-Cen complex and young local associations originated as a result of the impact of the inner spiral arm shock wave against a giant molecular cloud. The core of the giant molecular cloud formed the Sco-Cen complex, and some small cloudlets in a halo around the giant molecular cloud formed young local associations several million years later. We also propose a supernova in young local associations a few million years ago as the most likely candidate to have reheated the Local Bubble to its present temperature.

Context. Recent studies have suggested that moving groups have a dynamic or "resonant" origin. Under this hypothesis, these kinematic structures become a powerful tool for studying the large-scale structure and dynamics of the Milky Way. Aims. We aim to characterize these structures in the U-V-age-[Fe/H] space and establish observational constraints that will allow us to study their origin and evolution. Methods. We apply multiscale techniques -wavelet denoising (WD)- to an extensive compendium of more than 24000 stars in the solar neighbourhood with the best available astrometric, photometric and spectroscopic data. Results. We confirm that the dominant structures in the U-V plane are the branches of Sirius, Coma Berenices, Hyades-Pleiades and Hercules, which are nearly equidistant in this kinematic plane and show a negative slope. The abrupt drops in the velocity distribution are characterized. We find a certain dependence of these kinematic structures on Galactic position with a significant change of contrast among substructures inside the branches. A large spread of ages is observed for all branches. The Hercules branch is detected in all subsamples with ages older than ~ 2 Gyr and the set of the other three branches is well established for stars > 400 Myr. The age-metallicity relation of each branch is examined and the relation between kinematics and metallicity is studied. Conclusions. Not all of these observational constraints are successfully explained by the recent models proposed for the formation of such kinematic structures. Simulations incorporating stellar ages and metallicities are essential for future studies. The comparison of the observed and simulated distributions obtained by WD will provide a physical interpretation of the existence of the branches in terms of local or large-scale dynamics. [Abridged]

We investigate the possible existence of an extended halo of early-type stars around Cygnus OB2, which is hinted at by near-infrared color-color diagrams, and its relationship to Cygnus OB2 itself, as well as to the nearby association Cygnus OB9 and to the star forming regions in the Cygnus X North complex. A total of 96 early-type stars are identified in the targeted region, which amounts to nearly half of the observed sample. Most of them have featureless near-infrared spectra as expected from OB stars at the available resolution. Another 18 stars that display Brackett emission lines can be divided between evolved massive stars (most likely Be stars) and Herbig Ae/Be stars based on their infrared excesses. A component associated with Cygnus OB9/NGC 6910 is clearly identified, as well as an enhancement in the surface density of early-type stars at Cygnus X North. We also find a field population, consisting largely of early B giants and supergiants, which is probably the same as identified in recent studies of the inner 1-degree circle around Cygnus OB2. The age and large extension of this population discards a direct relationship with Cygnus OB2 or any other particular association. Earlier claims of the possible large extent of Cygnus OB2 beyond its central, very massive aggregate seem to be dismissed by our findings. The existence of a nearly ubiquitous population of evolved stars with massive precursors suggests a massive star formation history in Cygnus having started long before the formation of the currently observed OB associations in the region.

Context: Over the last decade, several groups of young (mainly low-mass) stars have been discovered in the solar neighbourhood (closer than ~100 pc), thanks to cross-correlation between X-ray, optical spectroscopy and kinematic data. These young local associations -- including an important fraction whose members are Hipparcos stars -- offer insights into the star formation process in low-density environments, shed light on the substellar domain, and could have played an important role in the recent history of the local interstellar medium. Aims: To study the kinematic evolution of young local associations and their relation to other young stellar groups and structures in the local interstellar medium, thus casting new light on recent star formation processes in the solar neighbourhood. Methods: We compiled the data published in the literature for young local associations. Using a realistic Galactic potential we integrated the orbits for these associations and the Sco-Cen complex back in time. Results: Combining these data with the spatial structure of the Local Bubble and the spiral structure of the Galaxy, we propose a recent history of star formation in the solar neighbourhood. We suggest that both the Sco-Cen complex and young local associations originated as a result of the impact of the inner spiral arm shock wave against a giant molecular cloud. The core of the giant molecular cloud formed the Sco-Cen complex, and some small cloudlets in a halo around the giant molecular cloud formed young local associations several million years later. We also propose a supernova in young local associations a few million years ago as the most likely candidate to have reheated the Local Bubble to its present temperature.

In this work we propose a scenario for the history of the recent star formation (during the last 20-30 Myr) in the nearest solar neighbourhood (~150 pc), from the study of the spatial and kinematic properties of the members of the so-called young local associations, the Sco-Cen complex and the Local Bubble, the most important structure observed in the local interstellar medium (ISM).

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the proven principles of ESA's Hipparcos mission to create an all-sky survey of about one billion stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. Through its massive measurement of stellar distances, motions and multi-colour photometry it will provide fundamental data necessary for unravelling the structure, formation and evolution of the Galaxy. This paper presents the design and performance of the broad- and medium-band set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands (extending from the ultraviolet to the far-red), the criteria, and the methodology on which this choice has been based are discussed in detail. We analyse the photometric capabilities for characterizing the luminosity, temperature, gravity and chemical composition of stars. We also discuss the automatic determination of these physical parameters for the large number of observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. Finally, the capability of the photometric system to deal with the main Gaia science case is outlined.

One of the most promising space missions of ESA is the astrometric satellite Gaia, which will provide very precise astrometry and multicolour photometry, for all 1.3 billion objects to V~20, and radial velocities with accuracies of a few km/s for most stars brighter than V ~ 17. Consequently, full homogeneous six-dimensional phase-space information for a huge number of stars will become available. Our Monte Carlo simulator has been used to estimate the number of white dwarfs potentially observable by Gaia. From this we assess which would be the white dwarf luminosity functions which Gaia will obtain and discuss in depth the scientific returns of Gaia in the specific field of white dwarf populations. Scientific attainable goals include, among others, a reliable determination of the age of the Galactic disk, a better knowledge of the halo of the Milky Way and the reconstruction of the star formation history of the Galactic disk. Our results also demonstrate the potential impact of a mission like Gaia in the current understanding of the white dwarf cooling theory.

We present the results of a narrow-band infrared imaging survey of a narrow strip (12' wide) around the galactic equator between 6 deg and 21 deg of galactic longitude aimed at detecting field stars with strong CO absorption, mainly late-type giants and supergiants. Our observations include follow-up low resolution spectroscopy (R = 980) of 191 selected candidates in the H and K bands. Most of these objects have photometric and spectroscopic characteristics consistent with them being red giants, and some display broad, strong absorption wings due to water vapor absorption between the H and K bands. We also identify in our sample 18 good supergiant candidates characterized by their lack of noticeable water absorption, strong CO bands in the H and K windows, and HK_S photometry suggestive of high intrinsic luminosity and extinctions reaching up to A_V ~40 mag. Another 9 additional candidates share the same features except for weak H2O absorption, which is also observed among some M supergiants in the solar neighbourhood. Interesting differences are noticed when comparing our stars to a local sample of late-type giants and supergiants, as well as to a sample of red giants in globular clusters of moderately subsolar metallicity and to a sample of bulge stars. (...) We propose that the systematic spectroscopic differences of our inner Galaxy stars are due to their higher metallicities that cause deeper mixing in their mantles, resulting in lower surface abundances of C and O and higher abundances of CN, which contribute to the strength of the CaI and NaI features at low resolution. Our results stress the limitations of using local stars as templates for the study of composite cool stellar populations such as central starbursts in galaxies (Abridged).

We report the discovery of a new Wolf-Rayet star in the direction of Cygnus. The star is strongly reddened but quite bright in the infrared, with J = 9.22, H = 8.08 and K = 7.09 (2MASS). On the basis of its H + K spectrum, we have classified WR 142a a WC8 star. We have estimated its properties using as a reference those of other WC8 stars in the solar neighbourhood as well as those of WR 135, whose near-infrared spectrum is remarkably similar. We thus obtain a foreground reddening of A(V) = 8.1 mag, M(J) = -4.3, log(L/Lo) = 5.0 - 5.2, R = 0.8 Ro, T = 125,000 K, M = 7.9 - 9.7 Mo, and a mass loss of (1.4 - 2.3)e-05 Mo/yr. The derived distance modulus, DM = 11.2 +/- 0.7 mag, places it in a region occupied by several OB associations in the Cygnus arm, and particularly in the outskirts of both Cygnus OB2 and Cygnus OB9. The position in the sky alone does not allow us to unambiguously assign the star to either association, but based on the much richer massive star content of Cygnus OB2 membership in this latter association appears to be more likely.

The young star velocity field is analysed by means of a galactic model which takes into account solar motion, differential galactic rotation and spiral arm kinematics. We use two samples of Hipparcos data, one containing O- and B-type stars and another one composed of Cepheid variable stars. The robustness of our method is tested through careful kinematic simulations. Our results show a galactic rotation curve with a classical value of $A$ Oort constant for the O and B star sample ($A^{\mathrm{OB}} =$ 13.7-13.8 km s$^{-1}$ kpc $^{-1}$) and a higher value for Cepheids ($A^{\mathrm{Cep}} =$ 14.9-16.9 km s$^{-1}$ kpc $^{-1}$, depending on the cosmic distance scale chosen). The second-order term is found to be small, compatible with a zero value. The study of the residuals shows the need for a $K$-term up to a heliocentric distance of 4 kpc, obtaining a value $K = -$(1-3) km s$^{-1}$ kpc$^{-1}$. The results obtained for the spiral structure from O and B stars and Cepheids show good agreement. The Sun is located relatively near the minimum of the spiral perturbation potential ($\psi_\odot =$ 284-20$\degr$) and very near the corotation circle. The angular rotation velocity of the spiral pattern was found to be $\Omega_{\mathrm{p}} \approx 30$ km s$^{-1}$ kpc $^{-1}$.

Two samples of O- and B-type stars and Cepheids with Hipparcos data have been used to characterize galactic rotation and spiral arm kinematics in the solar neighbourhood. An extensive set of simulations has been performed in order to assess the capabilities of the method and its sensitivity to sample errors and biases in the adopted parameters.

The ESA HIPPARCOS satellite has provided astrometry of unprecedented accuracy, allowing to reassess, improve and refine the pre-HIPPARCOS luminosity calibrations. We review the "classical" absolute magnitude calibrations with the Stroemgren-Crawford intermediate-band photometric system. A small zero point correction of about 2-4% seems necessary, as well as to refine the dependences on metallicity and projected rotational velocity. The need of a rigorous statistical treatment of the extremely precise HIPPARCOS to derive definite dependences of the luminosity on physical stellar parameters is emphasized.

A sample of O- and B-type stars with Hipparcos astrometric data, ages computed from Str\"omgren photometry and radial velocities, has been used to characterize the structure, age and kinematics of the Gould Belt system. The local spiral structure of our galaxy is determined from this sample, and also from a sample of Hipparcos Cepheid stars. The Gould Belt, with an orientation with respect to the galactic plane of $i_{\mathrm{G}} = 16$-$22^\circ$ and $\Omega_{\mathrm{G}} = 275$-$295^\circ$, extends up to a distance of 600 pc from the Sun. Roughly the 60-65% of the O and B stars younger than 60 Myr in the solar neighbourhood belong to this structure. Our results indicate that the kinematical behaviour of this system is complex, with an expansion motion in the solar neighbourhood ($R <$ 300 pc). In the frame of the Lin's theory, and analysing the O and B stars further than 600 pc and the Cepheids, we found a galactic spiral structure characterized by a 4-arm spiral pattern with the Sun located at $\psi_\odot = 350$-$355 \pm 30^\circ$ -- near the Sagittarius-Carina arm -- and outside the corotation circle. The angular rotation speed of the spiral pattern was found to be $\Omega_{\mathrm{p}} = 31$-$32 \pm 4$ km s$^{-1}$ kpc$^{-1}$.

The local velocity field of young stars is dominated by the galactic rotation, the kinematics of the Gould Belt and the nearest OB associations and open clusters, and the kinematics of the spiral structure. We re-examined here this local velocity field by using a large sample of nearby O and B stars from the Hipparcos Catalogue. The high quality astrometric data are complemented with a careful compilation of radial velocities and Str\"omgren photometry, which allows individual photometric distances and ages to be derived. The Gould Belt extends up to 600 pc from the Sun with an inclination with respect to the galactic plane of $i_{\mathrm{G}} = 16$-$22\degr$ and the ascending node placed at $\Omega_{\mathrm{G}} = 275$-$295\degr$. Approximately 60% of the stars younger than 60 Myr belong to this structure. The values found for the Oort constants when different samples selected by age or distance were used allowed us to interpret the systematic trends observed as signatures induced by the kinematic behaviour of the Gould Belt. The contribution of Sco-Cen and Ori OB1 complexes in the characterization of the expansion of the Gould Belt system is also discussed. We found that a positive $K$-term remains when these aggregates are excluded. From the kinematic behaviour of the stars and their spatial distribution we derive an age for the Gould Belt system in the interval 30-60 Myr.

The disruption of stellar systems, such as open clusters or stellar complexes, stands out as one of the most reasonable physical processes accounting for the young moving groups observed in the solar neighbourhood. In the present study we analyse some of the mechanisms that are important in the kinematic evolution of a group of unbound stars, such as the focusing phenomenon and its ability to recover the observed moving group's velocity dispersions, and the efficiency of disc heating and galactic differential rotation in disrupting unbound stellar systems. Our main tools used to perform this analysis are both the epicycle theory and the integration of the equations of motion using a realistic gravitational potential of the Galaxy. The study of the trajectories followed by stars in each of the Pleiades moving group substructures found by Asiain et al. (1999) allows us to determine their stellar spatial and velocity distribution evolution. The kinematic properties of these substructures are compared to those of a simulated stellar complex which has evolved under the influence of the galactic gravitational potential and the disc heating. We conclude that a constant diffusion coefficient compatible with the observational heating law is able to explain the velocity and spatial dispersions of the Pleiades moving group substructures that are younger than $\sim 1.5 \cdot 10^8$ yr.

The distance modulus of the Large Magellanic Cloud (LMC) is the first step on the determination of the cosmic distance scale. In this paper it is obtained from a new calibration of the classical Cepheid Period-Luminosity (PL) relation. Our new calibration is obtained by applying the LM method -- a parametric method designed to obtain luminosity calibrations from astrometric information -- to the Hipparcos data for Classical Cepheids. The resulting distance modulus is $18.35^{m}\pm 0.13$.

The LM method, designed to exploit the Hipparcos data to obtain luminosity calibrations, is applied to derive luminosity calibrations for RR Lyrae and classical Cepheids. From these calibrations the distance to the Large Magellanic Cloud (LMC) is estimated. The distance moduli provided by the two calibrations are in good agreement, giving a value of 18.3m, while several previous calibrations using Hipparcos data provided inconsistent results between both types of stars. This result suggests that the Hubble constant should have a value of H_0 = 79 km s-1 Mpc-1.

The application of the Stromgren photometric luminosity calibrations to different types of CP stars is reexamined in the light of the new Hipparcos data. A first attempt is made to use the LM statistical parallax method (Luri et al., 1996) -- based on the maximum likelihood principle -- to obtain a calibration of the absolute magnitude as a function of two Stromgren colour indices, thus reflecting effective temperature and evolution. Its application to a sample of Si stars and to a sample of normal main sequence stars in the same spectral range allows us to compare the calibrations obtained and to discuss the position of Si stars in the HR diagram. Additionally, a sample of bonafide, spectroscopically selected Am stars together with normal main sequence stars are used to derive a new absolute magnitude calibration for late A-type main sequence stars, taking into account the effects of evolution, metallicity and stellar rotation.

The position in the HR diagram and the kinematic characteristics of different kinds of CP stars of the upper main sequence are obtained using the LM method (Luri et al., 1996). Most of the CP stars are main sequence stars occupying the whole width of the sequence. From a kinematic point of view, they belong to the young disk population (ages < 1.5 Gyr). It has also been found that, on kinematic grounds, the behaviour of lambda Bootis stars is similar to the one observed for normal stars of the same spectral range. On the other hand, roAp and noAp stars show the same kinematic characteristics. The peculiar velocity distribution function has been decomposed into a sum of three dimensional gaussians and the presence of Pleiades, Sirius and Hyades moving groups has been clearly established. Finally, a small number of CP stars are found to be high-velocity objects.

We present new optical and infrared photometric observations and high resolution H$\alpha$ spectra of the periodic radio star \lsi. The optical photometric data set covers the time interval 1985-1993 and amounts to about a hundred nights. A period of $\sim$26 days is found in the V band. The infrared data also present evidence for a similar periodicity, but with higher amplitude of variation (0\rmag 2). The spectroscopic observations include 16 intermediate and high dispersion spectra of \lsi\ collected between January 1989 and February 1993. The H$\alpha$ emission line profile and its variations are analyzed. Several emission line parameters -- among them the H$\alpha$ EW and the width of the H$\alpha$ red hump -- change strongly at or close to radio maximum, and may exhibit periodic variability. We also observe a significant change in the peak separation. The H$\alpha$ profile of \lsi\ does not seem peculiar for a Be star. However, several of the observed variations of the H$\alpha$ profile can probably be associated with the presence of the compact, secondary star.