Aims. We aim at creating a Bayesian model to infer the coefficients of PL or PLZ relations that propagates uncertainties in the observables in a rigorous and well founded way. Methods. We propose a directed acyclic graph to encode the conditional probabilities of the inference model that will allow us to infer probability distributions for the PL and PL(Z) relations. We evaluate the model with several semi-synthetic data sets and apply it to a sample of 200 fundamental mode and first overtone mode RR Lyrae stars for which Gaia DR1 parallaxes and literature Ks-band mean magnitudes are available. We define and test several hyperprior probabilities to verify their adequacy and check the sensitivity of the solution with respect to the prior choice. Results. The main conclusion of this work is the absolute necessity of incorporating the existing correlations between the observed variables (periods, metallicities and parallaxes) in the form of model priors in order to avoid systematically biased results, especially in the case of non-negligible uncertainties in the parallaxes. The tests with the semi-synthetic data based on the data set used in Gaia Collaboration et al. (2017) reveal the significant impact that the existing correlations between parallax, metallicity and periods have on the inferred parameters. The relation coefficients obtained here have been superseded by those presented in Muraveva et al. (2018a), that incorporates the findings of this work and the more recent Gaia DR2 measurements.

We present results from the analysis of 401 RR Lyrae stars (RRLs) belonging to the field of the Milky Way (MW). For a fraction of them multi-band ($V$, $K_{\rm s}$, $W1$) photometry, metal abundances, extinction values and pulsation periods are available in the literature and accurate trigonometric parallaxes measured by the Gaia mission alongside Gaia $G$-band time-series photometry have become available with the Gaia second data release (DR2) on 2018 April 25. Using a Bayesian fitting approach we derive new near-, mid-infrared period-absolute magnitude-metallicity ($PMZ$) relations and new absolute magnitude-metallicity relations in the visual ($M_V - {\rm [Fe/H]}$) and $G$ bands ($M_G - {\rm [Fe/H]}$), based on the Gaia DR2 parallaxes. We find the dependence of luminosity on metallicity to be higher than usually found in the literature, irrespective of the passband considered. Running the adopted Bayesian model on a simulated dataset we show that the high metallicity dependence is not caused by the method, but likely arises from the actual distribution of the data and the presence of a zero-point offset in the Gaia parallaxes. We infer a zero-point offset of $-0.057$ mas, with the Gaia DR2 parallaxes being systematically smaller. We find the RR Lyrae absolute magnitude in the $V$, $G$, $K_{\rm s}$ and $W1$ bands at metallicity of [Fe/H]=$-1.5$ dex and period of P = 0.5238 days, based on Gaia DR2 parallaxes to be $M_V = 0.66\pm0.06$ mag, $M_G = 0.63\pm0.08$ mag, $M_{K_{\rm s}} = -0.37\pm0.11$ mag and $M_{W1} = -0.41\pm0.11$ mag, respectively.

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.

The second Gaia data release (GDR2) provides precise five-parameter astrometric data (positions, proper motions and parallaxes) for an unprecedented amount of sources (more than $1.3$ billion, mostly stars). The use of this wealth of astrometric data comes with a specific challenge: how does one properly infer from these data the astrophysical parameters of interest? The main - but not only - focus of this paper is the issue of the estimation of distances from parallaxes, possibly combined with other information. We start with a critical review of the methods traditionally used to obtain distances from parallaxes and their shortcomings. Then we provide guidelines on how to use parallaxes more efficiently to estimate distances by using Bayesian methods. In particular also we show that negative parallaxes, or parallaxes with relatively larger uncertainties still contain valuable information. Finally, we provide examples that show more generally how to use astrometric data for parameter estimation, including the combination of proper motions and parallaxes and the handling of covariances in the uncertainties. The paper contains examples based on simulated Gaia data to illustrate the problems and the solutions proposed. Furthermore, the developments and methods proposed in the paper are linked to a set of tutorials included in the Gaia archive documentation that provide practical examples and a good starting point for the application of the recommendations to actual problems. In all cases the source code for the analysis methods is provided. Our main recommendation is to always treat the derivation of (astro-) physical parameters from astrometric data, in particular when parallaxes are involved, as an inference problem which should preferably be handled with a full Bayesian approach.

21 st century astrophysicists are confronted with the herculean task of distilling the maximum scientific return from extremely expensive and complex space- or ground-based instrumental projects. This paper concentrates in the mining of the time series catalog produced by the European Space Agency Gaia mission, launched in December 2013. We tackle in particular the problem of inferring the true distribution of the variability properties of Cepheid stars in the Milky Way satellite galaxy known as the Large Magellanic Cloud (LMC). Classical Cepheid stars are the first step in the so-called distance ladder: a series of techniques to measure cosmological distances and decipher the structure and evolution of our Universe. In this work we attempt to unbias the catalog by modelling the aliasing phenomenon that distorts the true distribution of periods. We have represented the problem by a 2-level generative Bayesian graphical model and used a Markov chain Monte Carlo (MCMC) algorithm for inference (classification and regression). Our results with synthetic data show that the system successfully removes systematic biases and is able to infer the true hyperparameters of the frequency and magnitude distributions.

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.

The ESA Gaia mission provides a unique time-domain survey for more than one billion sources brighter than G=20.7 mag. Gaia offers the unprecedented opportunity to study variability phenomena in the Universe thanks to multi-epoch G-magnitude photometry in addition to astrometry, blue and red spectro-photometry, and spectroscopy. Within the Gaia Consortium, Coordination Unit 7 has the responsibility to detect variable objects, classify them, derive characteristic parameters for specific variability classes, and provide global descriptions of variable phenomena. We describe the variability processing and analysis that we plan to apply to the successive data releases, and we present its application to the G-band photometry results of the first 14 months of Gaia operations that comprises 28 days of Ecliptic Pole Scanning Law and 13 months of Nominal Scanning Law. Out of the 694 million, all-sky, sources that have calibrated G-band photometry in this first stage of the mission, about 2.3 million sources that have at least 20 observations are located within 38 degrees from the South Ecliptic Pole. We detect about 14% of them as variable candidates, among which the automated classification identified 9347 Cepheid and RR Lyrae candidates. Additional visual inspections and selection criteria led to the publication of 3194 Cepheid and RR Lyrae stars, described in Clementini et al. (2016). Under the restrictive conditions for DR1, the completenesses of Cepheids and RR Lyrae stars are estimated at 67% and 58%, respectively, numbers that will significantly increase with subsequent Gaia data releases. Data processing within the Gaia Consortium is iterative, the quality of the data and the results being improved at each iteration. The results presented in this article show a glimpse of the exceptional harvest that is to be expected from the Gaia mission for variability phenomena. [abridged]