The unprecedented astrometry from Gaia DR2 provides us with an opportunity to study in detail molecular clouds in the solar neighbourhood. Extracting the wealth of information in these data remains a challenge, however. We have further improved our Gaussian Processes-based, three-dimensional dust mapping technique to allow us to study molecular clouds in more detail. These improvements include a significantly better scaling of the computational cost with the number of stars, and taking into account distance uncertainties to individual stars. Using Gaia DR2 astrometry together with 2MASS and WISE photometry for 30 000 stars, we infer the distribution of dust out to 600 pc in the direction of the Orion A molecular cloud. We identify a bubble-like structure in front of Orion A, centred at a distance of about 350 pc from the Sun. The main Orion A structure is visible at slightly larger distances, and we clearly see a tail extending over 100 pc that is curved and slightly inclined to the line-of-sight. The location of our foreground structure coincides with 5-10 Myr old stellar populations, suggesting a star formation episode that predates that of the Orion Nebula Cluster itself. We identify also the main structure of the Orion B molecular cloud, and in addition discover a background component to this at a distance of about 460 pc from the Sun. Finally, we associate our dust components at different distances with the plane-of-the-sky magnetic field orientation as mapped by Planck. This provides valuable information for modelling the magnetic field in 3D around star forming regions.

Determining the mechanisms that drive the evolution of protoplanetary disks is a necessary step to understand how planets form. Here we measured the mass accretion rate for young stellar objects at age >5 Myr, a critical test for the current models of disk evolution. We present the analysis of the spectra of 36 targets in the ~5-10 Myr old Upper Scorpius region for which disk masses were measured with ALMA. We find that the mass accretion rates in this sample of old but still survived disks are similarly high as those of the younger (<3 Myr old) star-forming regions of Lupus and Cha I, when considering the dependence on stellar and disk mass. In particular, several disks show high mass accretion rates >10^-9 Msun/yr while having low disk masses. Furthermore, the median values of the measured mass accretion rates in the disk mass ranges where our sample is complete at a level ~60-80% are compatible in these three regions. At the same time, the spread of mass accretion rates at any given disk mass is still >0.9 dex even at age>5 Myr. These results are in contrast with simple models of viscous evolution, which would predict that the values of the mass accretion rate diminish with time, and a tighter correlation with disk mass at age>5 Myr. Similarly, simple models of internal photoevaporation cannot reproduce the observed mass accretion rates, while external photoevaporation might explain the low disk masses and high accretion rates. A partial possible solution to the discrepancy with the viscous models is that the gas-to-dust ratio of the disks at >5 Myr is significantly different and higher than the canonical 100, as suggested by some dust and gas disk evolution models. The results shown here require the presence of several inter-playing processes, such as detailed dust evolution, external photoevaporation and possibly MHD winds, to explain the secular evolution of protoplanetary disks.

We present a study of the three dimensional structure, kinematics, and age distribution of the Orion OB association, based on the second data release of the Gaia satellite (Gaia DR2). Our goal is to obtain a complete picture of the star formation history of the Orion complex and to relate our findings to theories of sequential and triggered star formation. We select the Orion population with simple photometric criteria, and we construct a three dimensional map in galactic Cartesian coordinates to study the physical arrangement of the stellar clusters in the Orion region. The map shows structures that extend for roughly $150 \, \mathrm{pc}$ along the line of sight, divided in multiple sub-clusters. We separate different groups by using the density based clustering algorithm DBSCAN. We study the kinematic properties of all the groups found by DBSCAN first by inspecting their proper motion distribution, and then by applying a kinematic modelling code based on an iterative maximum likelihood approach, which we use to derive their mean velocity, velocity dispersion and isotropic expansion. By using an isochrone fitting procedure we provide ages and extinction values for all the groups. We confirm the presence of an old population ($\sim 15$ Myr) towards the 25 Ori region, and we find that groups with ages of $12-15 \, \mathrm{Myr}$ are present also towards the Belt region. We notice the presence of a population of $\sim 10$ Myr also in front of the Orion A molecular cloud. Our findings suggest that star formation in Orion does not follow a simple sequential scenario, but instead consists of multiple events, which caused kinematic and physical sub-structure. To fully explain the detailed sequence of events, specific simulations and further radial velocity data are needed.

We study the three dimensional arrangement of young stars in the solar neighbourhood using the second release of the Gaia mission (Gaia DR2) and we provide a new, original view of the spatial configuration of the star forming regions within 500 pc from the Sun. By smoothing the star distribution through a gaussian filter, we construct three dimensional density maps for early-type stars (upper-main sequence, UMS) and pre-main sequence (PMS) sources. The PMS and the UMS samples are selected through a combination of photometric and astrometric criteria. A side product of the analysis is a three dimensional, G-band extinction map, which we use to correct our colour-magnitude diagram for extinction and reddening. Both density maps show three prominent structures, Scorpius-Centaurus, Orion, and Vela. The PMS map shows a plethora of lower mass star forming regions, such as Taurus, Perseus, Cepheus, Cassiopeia, and Lacerta, which are less visible in the UMS map, due to the lack of large numbers of bright, early-type stars. We report the finding of a candidate new open cluster towards $l, b \sim 218.5^{\circ}, -2^{\circ}$, which could be related to the Orion star forming complex. We estimate ages for the PMS sample and we study the distribution of PMS stars as a function of their age. We find that younger stars cluster in dense, compact clumps, and are surrounded by older sources, whose distribution is instead more diffuse. The youngest groups that we find are mainly located in Scorpius-Centaurus, Orion, Vela, and Taurus. Cepheus, Cassiopeia, and Lacerta are instead more evolved and less numerous. Finally, we find that the three dimensional density maps show no evidence for the existence of the ring-like structure which is usually referred to as the Gould Belt.

We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures.

We present a study of the membership of the Hyades open cluster, derive kinematically-modelled parallaxes of its members, and study the colour-absolute magnitude diagram of the cluster. We use Gaia DR1 Tycho-Gaia Astrometric Solution (TGAS) data complemented by Hipparcos-2 data for bright stars not contained in TGAS. We supplement the astrometric data with radial velocities collected from a dozen literature sources. By assuming that all cluster members move with the mean cluster velocity to within the velocity dispersion, we use the observed and the expected motions of the stars to determine individual cluster membership probabilities. We subsequently derive improved parallaxes through maximum-likelihood kinematic modelling of the cluster. This method has an iterative component to deal with 'outliers', caused for instance by double stars or escaping members. Our method extends an existing method and supports the mixed presence of stars with and without radial velocities. We find 251 candidate members, 200 of which have a literature radial velocity, and 70 of which are new candidate members with TGAS astrometry. The cluster is roughly spherical in its centre but significantly flattened at larger radii. The observed colour-absolute magnitude diagram shows a clear binary sequence. The kinematically-modelled parallaxes that we derive are a factor ~1.7 / 2.9 more precise than the TGAS / Hipparcos-2 values and allow to derive an extremely sharp main sequence. This sequence shows evidence for fine-detailed structure which is elegantly explained by the full spectrum turbulence model of convection.

Context. The origin and dynamical evolution of star clusters is an important topic in stellar astrophysics. Several models have been proposed to understand the formation of bound and unbound clusters and their evolution, and these can be tested by examining the kinematical and dynamical properties of clusters over a wide range of ages and masses. Aims. We use the Gaia-ESO Survey products to study four open clusters (IC 2602, IC 2391, IC 4665, and NGC 2547) that lie in the age range between 20 and 50 Myr. Methods. We employ the gravity index $\gamma$ and the equivalent width of the lithium line at 6708 $\AA$, together with effective temperature $\rm{T_{eff}}$, and the metallicity of the stars in order to discard observed contaminant stars. Then, we derive the cluster radial velocity dispersions $\sigma_c$, the total cluster mass $\rm{M}_{tot}$, and the half mass radius $r_{hm}$. Using the $Gaia$-DR1 TGAS catalogue, we independently derive the intrinsic velocity dispersion of the clusters from the astrometric parameters of cluster members. Results. The intrinsic radial velocity dispersions derived by the spectroscopic data are larger than those derived from the TGAS data, possibly due to the different masses of the considered stars. Using $\rm{M}_{tot}$ and $r_{hm}$ we derive the virial velocity dispersion $\sigma_{vir}$ and we find that three out of four clusters are supervirial. This result is in agreement with the hypothesis that these clusters are dispersing, as predicted by the "residual gas expulsion" scenario. However, recent simulations show that the virial ratio of young star clusters may be overestimated if it is determined using the global velocity dispersion, since the clusters are not fully relaxed.

We describe the methodologies that, taking advantage of Gaia-DR1 and the Gaia-ESO Survey data, enable the comparison of observed open star cluster sequences with stellar evolutionary models. The final, long-term goal is the exploitation of open clusters as age calibrators. We perform a homogeneous analysis of eight open clusters using the Gaia-DR1 TGAS catalogue for bright members, and information from the Gaia-ESO Survey for fainter stars. Cluster membership probabilities for the Gaia-ESO Survey targets are derived based on several spectroscopic tracers. The Gaia-ESO Survey also provides the cluster chemical composition. We obtain cluster parallaxes using two methods. The first one relies on the astrometric selection of a sample of bona fide members, while the other one fits the parallax distribution of a larger sample of TGAS sources. Ages and reddening values are recovered through a Bayesian analysis using the 2MASS magnitudes and three sets of standard models. Lithium depletion boundary (LDB) ages are also determined using literature observations and the same models employed for the Bayesian analysis. For all but one cluster, parallaxes derived by us agree with those presented in Gaia Collaboration et al. (2017), while a discrepancy is found for NGC 2516; we provide evidence supporting our own determination. Inferred cluster ages are robust against models and are generally consistent with literature values. The systematic parallax errors inherent in the Gaia DR1 data presently limit the precision of our results. Nevertheless, we have been able to place these eight clusters onto the same age scale for the first time, with good agreement between isochronal and LDB ages where there is overlap. Our approach appears promising and demonstrates the potential of combining Gaia and ground-based spectroscopic datasets.

We use the first data release of the Gaia mission to explore the three dimensional arrangement and the age ordering of the many stellar groups towards the Orion OB association, aiming at a new classification and characterization of the stellar population. We make use of the parallaxes and proper motions provided in the Tycho Gaia Astrometric Solution (TGAS) sub-set of the Gaia catalogue, and of the combination of Gaia and 2MASS photometry. In TGAS we find evidence for the presence of a young population, at a parallax $\varpi \sim 2.65 \, \mathrm{mas}$, loosely distributed around some known clusters: 25 Ori, $\epsilon$ Ori and $\sigma$ Ori, and NGC 1980 ($\iota$ Ori). The low mass counterpart of this population is visible in the color-magnitude diagrams constructed by combining Gaia and 2MASS photometry. We study the density distribution of the young sources in the sky. We find the same groups as in TGAS, and also some other density enhancements that might be related to the recently discovered Orion X group, the Orion dust ring, and to the $\lambda$ Ori complex. We estimate the ages of this population and we infer the presence of an age gradient going from 25 Ori (13-15 Myr) to the ONC (1-2 Myr). We confirm this age ordering by repeating the Bayesian fit using the Pan-STARRS1 data. The estimated ages towards the NGC 1980 cluster span a broad range of values. This can either be due to the presence of two populations coming from two different episodes of star formation or to a large spread along the line of sight of the same population. Our results form the first step towards using the Gaia data to unravel the complex star formation history of the Orion region in terms of the different star formation episodes, their duration, and their effects on the surrounding interstellar medium.

Understanding how young stars and their circumstellar disks form and evolve is key to explain how planets form. The evolution of the star and the disk is regulated by different processes, both internal to the system or related to their environment. The former include accretion of material onto the central star, wind emission, and photoevaporation of the disk due to high-energy radiation from the central star. These are best studied spectroscopically, and the distance to the star is a key parameter in all these studies. Here we present new estimates of the distance to a complex of nearby star-forming clouds obtained combining TGAS distances with measurement of extinction on the line of sight. Furthermore, we show how we plan to study the effects of the environment on the evolution of disks with Gaia, using a kinematic modelling code we have developed to model young star-forming regions.

We present optical depth and temperature maps of the Perseus molecular cloud, obtained combining dust emission data from the Herschel and Planck satellites and 2MASS/NIR dust extinction maps. The maps have a resolution of 36 arcsec in the Herschel regions, and of 5 arcmin elsewhere. The dynamic range of the optical depth map ranges from $1\times10^{-2}\, \mathrm{mag}$ up to $20 \,\mathrm{mag}$ in the equivalent K band extinction. We also evaluate the ratio between the $2.2 \,\mathrm{\mu m}$ extinction coefficient and the $850 \,\mathrm{\mu m}$ opacity. The value we obtain is close to the one found in the Orion B molecular cloud. We show that the cumulative and the differential area function of the data (which is proportional to the probability distribution function of the cloud column density) follow power laws with index respectively $\simeq -2$, and $\simeq -3$. We use WISE data to improve current YSO catalogues based mostly on \emph{Spitzer} data and we build an up-to-date selection of Class~I/0 objects. Using this selection, we evaluate the local Schmidt law, $\Sigma_{\mathrm{YSO}} \propto \Sigma_{\mathrm{gas}}^{\beta}$, showing that $\beta=2.4 \pm 0.6$. Finally, we show that the area-extinction relation is important for determining the star formation rate in the cloud, which is in agreement with other recent works.