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.

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.

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.

RW Aur is a young binary star that experienced a deep dimming in 2010-11 in component A and a second even deeper dimming from summer 2014 to summer 2016. We present new unresolved multi-band photometry during the 2014-16 eclipse, new emission line spectroscopy before and during the dimming, archive infrared photometry between 2014-15, as well as an overview of literature data. Spectral observations were carried out with the Fibre-fed RObotic Dual-beam Optical Spectrograph on the Liverpool Telescope. Photometric monitoring was done with the Las Cumbres Observatory Global Telescope Network and James Gregory Telescope. Our photometry shows that RW Aur dropped in brightness to R = 12.5 in March 2016. In addition to the long-term dimming trend, RW Aur is variable on time scales as short as hours. The short-term variation is most likely due to an unstable accretion flow. This, combined with the presence of accretion-related emission lines in the spectra suggest that accretion flows in the binary system are at least partially visible during the eclipse. The equivalent width of [O I] increases by a factor of ten in 2014, coinciding with the dimming event, confirming previous reports. The blue-shifted part of the $H\alpha$ profile is suppressed during the eclipse. In combination with the increase in mid-infrared brightness during the eclipse reported in the literature and seen in WISE archival data, and constraints on the geometry of the disk around RW Aur A we arrive at the conclusion that the obscuring screen is part of a wind emanating from the inner disk.

Given that a majority of stars form in multiple systems, in order to fully understand the star- and planet-formation processes we must seek to understand them in multiple stellar systems. With this in mind, we present an analysis of the enigmatic binary T-Tauri system VV Corona Australis, in which both components host discs, but only one is visible at optical wavelengths. We seek to understand the peculiarities of this system by searching for a model for the binary which explains all the available continuum observations of the system. We present new mid-infrared interferometry and near-infrared spectroscopy along with archival millimetre-wave observations, which resolve the binary at 1.3mm for the first time. We compute a grid of pre-main-sequence radiative transfer models and calculate their posterior probabilities given the observed spectral energy distributions and mid-infrared interferometric visibilities of the binary components, beginning with the assumption that the only differences between the two components are their inclination and position angles. Our best-fitting solution corresponds to a relatively low luminosity T-tauri binary, with each component's disc having a large scale height and viewed at moderate inclination ($\sim 50^\circ$), with the infrared companion inclined by $\sim 5^\circ$ degrees more than the primary. Comparing the results of our model to evolutionary models suggests stellar masses $\sim 1.7$M$_\odot$ and an age for the system of 3.5Myr, towards the upper end of previous estimates. Combining these results with accretion indicators from near-IR spectroscopy, we determine an accretion rate of 4.0$\times 10^{-8}$M$_\odot$ yr$^{-1}$ for the primary. We suggest that future observations of VV~CrA and similar systems should prioritise high angular resolution sub-mm and near-IR imaging of the discs and high resolution optical/NIR spectroscopy of the central stars.

There is increasing evidence that episodic accretion is a common phenomenon in Young Stellar Objects (YSOs). Recently, the source HOPS 383 in Orion was reported to have a $\times 35$ mid-infrared -- and bolometric -- luminosity increase between 2004 and 2008, constituting the first clear example of a class 0 YSO (a protostar) with a large accretion burst. The usual assumption that in YSOs accretion and ejection follow each other in time needs to be tested. Radio jets at centimeter wavelengths are often the only way of tracing the jets from embedded protostars. We searched the Very Large Array archive for the available observations of the radio counterpart of HOPS 383. The data show that the radio flux of HOPS 383 varies only mildly from January 1998 to December 2014, staying at the level of $\sim 200$ to 300 $\mu$Jy in the X band ($\sim 9$ GHz), with a typical uncertainty of 10 to 20 $\mu$Jy in each measurement. We interpret the absence of a radio burst as suggesting that accretion and ejection enhancements do not follow each other in time, at least not within timescales shorter than a few years. Time monitoring of more objects and specific predictions from simulations are needed to clarify the details of the connection betwen accretion and jets/winds in YSOs.

Variability in emission lines is a characteristic feature in young stars and can be used as a tool to study the physics of the accretion process. Here we present a study of H{\alpha} variability in 15 T Tauri and Herbig Ae stars (K7-B2) over a wide range of time windows, from minutes, to hours, to days, and years. We assess the variability using linewidth measurements and the time series of line profiles. All objects show gradual, slow profile changes on time-scales of days. In addition, in three cases there is evidence for rapid variations in H{\alpha} with typical time-scales of 10 min, which occurs in 10% of the total covered observing time. The mean accretion-rate changes, inferred from the line fluxes,are 0.01-0.07 dex for time-scales of < 1 hour, 0.04-0.4 dex for time-scales of days, and 0.13-0.52 dex for time-scales of years. In Costigan et al. 2012 we derived an upper limit finding that the intermediate (days) variability dominated over longer (years) variability. Here our new results, based on much higher cadence observations, also provide a lower limit to accretion-rate variability on similar time-scales (days), thereby constraining the accretion rate variability physics in a much more definitive way. A plausible explanation for the gradual variations over days is an asymmetric accretion flow resulting in a rotational modulation of the accretion-related emission, although other interpretations are possible as well. In conjunction with our previous work, we find that the time-scales and the extent of the variability is similar for objects ranging in mass from ~ 0.1 to ~ 5 Msol. This confirms that a single mode of accretion is at work from T Tauri to Herbig Ae stars - across a wide range of stellar masses.

We report JVLA 8-10 GHz ($\lambda$=3.0-3.7 cm) monitoring observations toward the YSO cluster R Coronae Australis (R\,CrA), taken in 2012, from March 15 to September 12. These observations were planned to measure the radio flux variabilities in timescales from 0.5 hours to several days, to tens of days, and up to $\sim$200 days. We found that among the YSOs detectable in individual epochs, in general, the most reddened objects in the \textit{Spitzer} observations show the highest mean 3.5 cm Stokes \textit{I} emission, and the lowest fractional variabilities on $<$200-day timescales. The brightest radio flux emitters in our observations are the two reddest sources IRS7W and IRS7E. In addition, by comparing with observations taken in 1996-1998 and 2005, we found that the radio fluxes of these two sources have increased by a factor $\sim$1.5. The mean 3.5-cm fluxes of the three Class I/II sources IRSI, IRS2, and IRS6 appear to be correlated with their accretion rates derived by a previous near infrared line survey. The weakly accreting Class I/II YSOs, or those in later evolutionary stages, present radio flux variability on $<$0.5-hour timescales. Some YSOs were detected only during occasional flaring events. The source R\,CrA went below our detection limit during a few fading events.

Characterization of the properties of young brown dwarfs are important to constraining the formation of objects at the extreme low-mass end of the IMF. While young brown dwarfs share many properties with solar-mass T Tauri stars, differences may be used as tests of how the physics of accretion/outflow and disk chemistry/dissipation depend on the mass of the central object. This article summarizes the presentations and discussions during the splinter session on 'Disks, accretion and outflows of brown dwarfs' held at the CoolStars17 conference in Barcelona in June 2012. Recent results in the field of brown dwarf disks and outflows include the determination of brown dwarf disk masses and geometries based on Herschel far-IR photometry (70-160 um), accretion properties based on X-Shooter spectra, and new outflow detections in the very low-mass regime.

We present the results of a variability study of accreting young stellar objects in the Chameleon I star-forming region which is based on ~300 high resolution optical spectra from the multi-object fibre spectrograph FLAMES/GIRAFFE at the ESO/VLT. Twenty five objects with spectral types from G2-M5.75 were observed 12 times over the course of 15 months. Using the emission lines Ha (6562.81 A) and Ca II (8662.1 A) as accretion indicators we found 10 accreting and 15 non-accreting objects. We derived accretion rates for all accretors in the sample using the Ha equivalent width, Ha 10% width and the CaII equivalent width. The mean amplitude of variations in derived accretion rate from Ha equivalent width was ~ 0.37 dex, from Ca II equivalent width ~0.83 dex and from Ha 10% width ~1.11 dex. Based on the large amplitude of variations in accretion rates derived from the Ha 10% width with respect to the other diagnostics, we do not consider it to be a reliable accretion rate estimator. Taking the variations in Ha equivalent width and CaII equivalent width accretion rates to be closer to the true value, they suggest that the spread which has been found around the accretion rate to stellar mass relation is not due to the variability of individual objects on time-scales of weeks to ~1 year. From these variations we can also infer that the accretion rates are stable within < 0.37 dex over time-scales of less than 15 months. A major portion of the accretion variability was found to occur on less than the shortest time-scales in our observations, 8-25 days, which is comparable with the rotation periods of these young stellar objects. This could be an indication that what we are probing is spatial structure in the accretion flows, and also suggests that observations on time-scales of ~a couple of weeks are sufficient to limit the total extent of accretion rate variations in typical young stars.

FU Tau A is a young very low mass object in the Taurus star forming region which was previously found to have strong X-ray emission and to be anomalously bright for its spectral type. In this study we discuss these characteristics using new information from quasi-simultaneous photometric and spectroscopic monitoring. From photometric time series obtained with the 2.2m telescope on Calar Alto we measure a period of ~4d for FU Tau A, most likely the rotation period. The short-term variations over a few days are consistent with the rotational modulation of the flux by cool, magnetically induced spots. In contrast, the photometric variability on timescales of weeks and years can only be explained by the presence of hot spots, presumably caused by accretion. The hot spot properties are thus variable on timescales exceeding the rotation period, maybe due to long-term changes in the accretion rate or geometry. The new constraints from the analysis of the variability confirm that FU Tau A is affected by magnetically induced spots and excess luminosity from accretion. However, accretion is not sufficient to explain its anomalous position in the HR diagram. In addition, suppressed convection due to magnetic activity and/or an early evolutionary stage need to be invoked to fully account for the observed properties. These factors cause considerable problems in estimating the mass of FU Tau A and other objects in this mass/age regime, to the extent that it appears questionable if it is feasible to derive the Initial Mass Function for young low-mass stars and brown dwarfs.