We present ten new ultra-cool dwarfs in seven wide binary systems discovered using $\textit{Gaia}$ DR2 data, identified as part of our $\textit{Gaia}$ Ultra-Cool Dwarf Sample project. The seven systems presented here include an L1 companion to the G5 IV star HD 164507, an L1: companion to the V478 Lyr AB system, an L2 companion to the metal-poor K5 V star CD-28 8692, an M9 V companion to the young variable K0 V star LT UMa, and three low-mass binaries consisting of late Ms and early Ls. The HD 164507, CD-28 8692, V478 Lyr, and LT UMa systems are particularly important benchmarks, because the primaries are well characterised and offer excellent constraints on the atmospheric parameters and ages of the companions. We find that the M8 V star 2MASS J23253550+4608163 is $\sim$2.5 mag overluminous compared to M dwarfs of similar spectral type, but at the same time it does not exhibit obvious peculiarities in its near-infrared spectrum. Its overluminosity cannot be explained by unresolved binarity alone. Finally, we present an L1+L2 system with a projected physical separation of 959 au, making this the widest L+L binary currently known.

Corona-Australis is one of the nearest regions to the Sun with recent and ongoing star formation, but the current picture of its stellar (and substellar) content is not complete yet. We take advantage of the second data release of the Gaia space mission to revisit the stellar census and search for additional members of the young stellar association in Corona-Australis. We applied a probabilistic method to infer membership probabilities based on a multidimensional astrometric and photometric data set over a field of 128 deg$^{2}$ around the dark clouds of the region. We identify 313 high-probability candidate members to the Corona-Australis association, 262 of which had never been reported as members before. Our sample of members covers the magnitude range between $G\gtrsim5$ mag and $G\lesssim20$ mag, and it reveals the existence of two kinematically and spatially distinct subgroups. There is a distributed `off-cloud' population of stars located in the north of the dark clouds that is twice as numerous as the historically known `on-cloud' population that is concentrated around the densest cores. By comparing the location of the stars in the HR-diagram with evolutionary models, we show that these two populations are younger than 10 Myr. Based on their infrared excess emission, we identify 28 Class II and 215 Class III stars among the sources with available infrared photometry, and we conclude that the frequency of Class II stars (i.e. `disc-bearing' stars) in the on-cloud region is twice as large as compared to the off-cloud population. The distance derived for the Corona-Australis region based on this updated census is $d=149.4^{+0.4}_{-0.4}$ pc, which exceeds previous estimates by about 20 pc.In this paper we provide the most complete census of stars in Corona-Australis available to date that can be confirmed with Gaia data.

Aims:We take advantage of the second data release of the Gaia space mission and the state-of-the-art astrometry delivered from very long baseline interferometry observations to revisit the structure and kinematics of the nearby Taurus star-forming region. Methods: We apply a hierarchical clustering algorithm for partitioning the stars in our sample into groups (i.e., clusters) that are associated with the various molecular clouds of the complex, and derive the distance and spatial velocity of individual stars and their corresponding molecular clouds. Results: We show that the molecular clouds are located at different distances and confirm the existence of important depth effects in this region reported in previous studies. For example, we find that the L 1495 molecular cloud is located at $d=129.9^{+0.4}_{-0.3}$ pc, while the filamentary structure connected to it (in the plane of the sky) is at $d=160.0^{+1.2}_{-1.2}$ pc. We report B 215 and L 1558 as the closest ($d=128.5^{+1.6}_{-1.6}$ pc) and most remote ($d=198.1^{+2.5}_{-2.5}$ pc) substructures of the complex, respectively. The median inter-cloud distance is 25 pc and the relative motion of the subgroups is on the order of a few km/s. We find no clear evidence for expansion (or contraction) of the Taurus complex, but signs of the potential effects of a global rotation. Finally, we compare the radial velocity of the stars with the velocity of the underlying $^{13}$CO molecular gas and report a mean difference of $0.04\pm0.12$ km/s (with r.m.s. of 0.63 km/s) confirming that the stars and the gas are tightly coupled.

Context. The study of star formation is extremely challenging due to the lack of complete and clean samples of young, nearby clusters, and star forming regions. The recent Gaia DR2 catalogue complemented with the deep, ground based COSMIC DANCe catalogue offers a new database of unprecedented accuracy to revisit the membership of clusters and star forming regions. The 30 Myr open cluster IC 4665 is one of the few well-known clusters of this age and it is an excellent target where to test evolutionary models and study planetary formation. Aims. We aim to provide a comprehensive membership analysis of IC 4665 and to study the following properties: empirical isochrones, distance, magnitude distribution, present-day system mass function, and spatial distribution. Methods. We use the Gaia DR2 catalogue together with the DANCe catalogue to look for members using a probabilistic model of the distribution of the observable quantities in both the cluster and background populations. Results. We obtain a final list of 819 candidate members which cover a 12.4 magnitude range (7 < J < 19.4). We find that 50% are new candidates, and we estimate a conservative contamination rate of 20%. This unique sample of members allows us to obtain a present-day system mass function in the range of 0.02-6 Msun, which reveals a number of details not seen in previous studies. In addition, they favour a spherically symmetric spatial distribution for this young open cluster. Conclusions. Our membership analysis represents a significant increase in the quantity and quality (low-contamination) with respect to previous studies. As such, it offers an excellent opportunity to revisit other fundamental parameters such as the age.

More than half a million of the 1.69 billion sources in Gaia Data Release 2 (DR2) are published with photometric time series that exhibit light variations during the 22 months of observation. An all-sky classification of common high-amplitude pulsators (Cepheids, long-period variables, Delta Scuti / SX Phoenicis, and RR Lyrae stars) is provided for stars with brightness variations greater than 0.1 mag in G band. A semi-supervised classification approach was employed, firstly training multi-stage random forest classifiers with sources of known types in the literature, followed by a preliminary classification of the Gaia data and a second training phase that included a selection of the first classification results to improve the representation of some classes, before the improved classifiers were applied to the Gaia data. Dedicated validation classifiers were used to reduce the level of contamination in the published results. A relevant fraction of objects were not yet sufficiently sampled for reliable Fourier series decomposition, consequently classifiers were based on features derived from statistics of photometric time series in the G, BP, and RP bands, as well as from some astrometric parameters. The published classification results include 195,780 RR Lyrae stars, 150,757 long-period variables, 8550 Cepheids, and 8882 Delta Scuti / SX Phoenicis stars. All of these results represent candidates whose completeness and contamination are described as a function of variability type and classification reliability. Results are expressed in terms of class labels and classification scores, which are available in the vari_classifier_result table of the Gaia archive.

We identify and investigate known late M, L and T dwarfs in the Gaia second data release. This sample is being used as a training set in the Gaia data processing chain of the ultra cool dwarfs work package. We find 695 objects in the optical spectral range M8 to T6 with accurate Gaia coordinates, proper motions, and parallaxes which we combine with published spectral types and photometry from large area optical and infrared sky surveys. We find that 100 objects are in 47 multiple systems, of which 27 systems are published and 20 are new. These will be useful benchmark systems and we discuss the requirements to produce a complete catalog of multiple systems with an ultra cool dwarf component. We examine the magnitudes in the Gaia passbands and find that the GBP magnitudes are unreliable and should not be used for these objects. We examine progressively redder colour magnitude diagrams and see a notable increase in the main sequence scatter and a bivariate main sequence for old and young objects. We provide an absolute magnitude spectral sub-type calibration for G and GRP passbands along with linear fits over the range M8 L8 for other passbands.

Context. Ruprecht 147 is the oldest (2.5 Gyr) open cluster in the solar vicinity (< 300 pc), making it an important target for stellar evolution studies and exoplanet searches. Aims. Derive a census of members and the luminosity, mass, and spatial distributions of the cluster. Methods. We use an astro-photometric data set including all available information from the literature together with our own observations. We process the data with an updated version of an existent membership selection methodology. Results. We identify 259 high-probability candidate members, including 58 previously unreported. All these candidates cover the luminosity interval between G > 6 mag to i< 21 mag. The cluster luminosity and mass distributions are derived with an unprecedented level of details allowing us to recognize, among other features, the Wielen dip. The mass distribution in the low-mass regime drops sharply at 0.4 $M_{\odot}$ even though our data are sensitive to stellar masses down to 0.1 $M_{\odot}$, suggesting that most very-low-mass members left the cluster as the result of its dynamical evolution. In addition, the cluster is highly elongated (ellipticity $\sim$ 0.5) towards the galactic plane, and mass segregated. Conclusions. Our combined Gaia+DANCe data set allows us to obtain an extended list of cluster candidate members, and to derive luminosity, mass and projected spatial distributions in the oldest open cluster of the solar vicinity.

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.

The Gaia second Data Release (DR2) presents a first mapping of full-sky RR Lyrae stars and Cepheids observed by the spacecraft during the initial 22 months of science operations. The Specific Object Study (SOS) pipeline, developed to validate and fully characterise Cepheids and RR Lyrae stars (SOS Cep&RRL) observed by Gaia, has been presented in the documentation and papers accompanying the Gaia first Data Release. Here we describe how the SOS pipeline was modified to allow for processing the Gaia multiband (G, G_BP and G_RP) time series photometry of all-sky candidate variables and produce specific results for confirmed RR Lyrae stars and Cepheids that are published in the DR2 catalogue. The SOS Cep&RRL processing uses tools such as the period-amplitude and the period-luminosity relations in the G band. For the analysis of the Gaia DR2 candidates we also used tools based on the G_BP and G_RP photometry, such as the period-Wesenheit relation in (G,G_RP). Multiband time series photometry and characterisation by the SOS Cep&RRL pipeline are published in Gaia DR2 for 150,359 such variables (9,575 classified as Cepheids and 140,784 as RR Lyrae stars) distributed all over the sky. The sample includes variables in 87 globular clusters and 14 dwarf galaxies. To the best of our knowledge, as of 25 April 2018, variability of 50,570 of these sources (350 Cepheids and 50,220 RR Lyrae stars) is not known in the literature, hence likely they are new discoveries by Gaia. An estimate of the interstellar absorption is published for 54,272 fundamental-mode RR Lyrae stars from a relation based on the G-band amplitude and the pulsation period. Metallicities derived from the Fourier parameters of the light curves are also released for 64,932 RR Lyrae stars and 3,738 fundamental-mode classical Cepheids with period shorter than 6.3 days.

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.

Gaia DR2 provides a unique all-sky catalogue of 550'737 variable stars, of which 151'761 are long-period variable (LPV) candidates with G variability amplitudes larger than 0.2 mag (5-95% quantile range). About one-fifth of the LPV candidates are Mira candidates, the majority of the rest are semi-regular variable candidates. For each source, G, BP , and RP photometric time-series are published, together with some LPV-specific attributes for the subset of 89'617 candidates with periods in G longer than 60 days. We describe this first Gaia catalogue of LPV candidates, and present various validation checks. Various samples of LPVs were used to validate the catalogue: a sample of well-studied very bright LPVs with light curves from the AAVSO that are partly contemporaneous with Gaia light curves, a sample of Gaia LPV candidates with good parallaxes, the ASAS_SN catalogue of LPVs, and the OGLE catalogues of LPVs towards the Magellanic Clouds and the Galactic bulge. The analyses of these samples show a good agreement between Gaia DR2 and literature periods. The same is globally true for bolometric corrections of M-type stars. The main contaminant of our DR2 catalogue comes from young stellar objects (YSOs) in the solar vicinity (within ~1 kpc), although their number in the whole catalogue is only at the percent level. A cautionary note is provided about parallax-dependent LPV attributes published in the catalogue. This first Gaia catalogue of LPVs approximately doubles the number of known LPVs with amplitudes larger than 0.2 mag, despite the conservative candidate selection criteria that prioritise low contamination over high completeness, and despite the limited DR2 time coverage compared to the long periods characteristic of LPVs. It also contains a small set of YSO candidates, which offers the serendipitous opportunity to study these objects at an early stage of the Gaia data releases.

The Gaia Data Release 2 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability 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.

Aims. We develop, test and characterise of a new statistical tool (intelligent system) for the sifting and analysis of nearby young open cluster (NYOC) populations. Methods. Using a Bayesian formalism, this statistical tool is able to obtain the posterior distributions of parameters governing the cluster model. It also uses hierarchical bayesian models to establish weakly informative priors, and incorporates the treatment of missing values and non-homogeneous (heteroscedastic) observational uncertainties. Results. From simulations, we estimate that this statistical tool renders kinematic (proper motion) and photometric (luminosity) distributions of the cluster population with a contamination rate of $5.8 \pm 0.2$ %. The luminosity distributions and present day mass function agree with the ones found by Bouy et al. (2015b) on the completeness interval of the survey. At the probability threshold of maximum accuracy, the classifier recovers $\sim$ 90% of Bouy et al. (2015b) candidate members and finds 10% of new ones. Conclusions. A new statistical tool for the analysis of NYOC is introduced, tested and characterised. Its comprehensive modelling of the data properties allows it to get rid of the biases present in previous works. In particular, those resulting from the use of only completely observed (non-missing) data and the assumption of homoskedastic uncertainties. Also, its Bayesian framework allows it to properly propagate observational uncertainties into membership probabilities and cluster velocity and luminosity distributions. Our results are in a general agreement with those from the literature, although we provide the most up-to-date and extended list of candidate members of the Pleiades cluster.

Methods. We compute Bayesian evidences and Bayes Factors for a set of variations of the classical radial models by King (1962), Elson et al. (1987) and Lauer et al. (1995). The variations incorporate different degrees of model freedom and complexity, amongst which we include biaxial (elliptical) symmetry, and luminosity segregation. As a by-product of the model comparison, we obtain posterior distributions and maximum a posteriori estimates for each set of model parameters. Results. We find that the model comparison results depend on the spatial extent of the region used for the analysis. For a circle of 11.5 parsecs around the cluster centre (the most homogeneous and complete region), we find no compelling reason to abandon Kings model, although the Generalised King model, introduced in this work, has slightly better fitting properties. Furthermore, we find strong evidence against radially symmetric models when compared to the elliptic extensions. Finally, we find that including mass segregation in the form of luminosity segregation in the J band, is strongly supported in all our models. Conclusions. We have put the question of the projected spatial distribution of the Pleiades cluster on a solid probabilistic framework, and inferred its properties using the most exhaustive and least contaminated list of Pleiades candidate members available to date. Our results suggest however that this sample may still lack about 20% of the expected number of cluster members. Therefore, this study should be revised when the completeness and homogeneity of the data can be extended beyond the 11.5 parsecs limit. Such study will allow a more precise determination of the Pleiades spatial distribution, its tidal radius, ellipticity, number of objects and total mass.

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.

Gaia Data Release 1 contains parallaxes for more than 700 Galactic Cepheids and RR Lyrae stars, computed as part of the Tycho-Gaia Astrometric Solution (TGAS). We have used TGAS parallaxes, along with literature ($V, I, J, {K_\mathrm{s}}, W_1$) photometry and spectroscopy, to calibrate the zero point of the Period-Luminosity and Period-Wesenheit relations of classical and type II Cepheids, and the near-infrared Period-Luminosity, Period-Luminosity-Metallicity and optical Luminosity-Metallicity relations of RR Lyrae stars. In this contribution we briefly summarise results obtained by fitting these basic relations adopting different techniques that operate either in parallax or distance (absolute magnitude) space.

In this article we present an overview of the ESA Gaia mission and of the unprecedented impact that Gaia will have on the field of variable star research. We summarise the contents and impact of the first Gaia data release on the description of variability phenomena, with particular emphasis on pulsating star research. The Tycho-Gaia astrometric solution, although limited to 2.1 million stars, has been used in many studies related to pulsating stars. Furthermore a set of 3,194 Cepheids and RR Lyrae stars with their times series have been released. Finally we present the plans for the ongoing study of variable phenomena with Gaia and highlight some of the possible impacts of the second data release on variable, and specifically, pulsating stars.

We identify and investigate known ultracool stars and brown dwarfs that are being observed or indirectly constrained by the Gaia mission. These objects will be the core of the Gaia ultracool dwarf sample composed of all dwarfs later than M7 that Gaia will provide direct or indirect information on. We match known L and T dwarfs to the Gaia first data release, the Two Micron All Sky Survey and the Wide-field Infrared Survey Explorer AllWISE survey and examine the Gaia and infrared colours, along with proper motions, to improve spectral typing, identify outliers and find mismatches. There are 321 L and T dwarfs observed directly in the Gaia first data release, of which 10 are later than L7. This represents 45 % of all the known LT dwarfs with estimated Gaia G magnitudes brighter than 20.3 mag. We determine proper motions for the 321 objects from Gaia and the Two Micron All Sky Survey positions. Combining the Gaia and infrared magnitudes provides useful diagnostic diagrams for the determination of L and T dwarf physical parameters. We then search the Tycho-Gaia astrometric solution Gaia first data release subset to find any objects with common proper motions to known L and T dwarfs and a high probability of being related. We find 15 new candidate common proper motion systems.

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]

We present an overview of the Specific Objects Study (SOS) pipeline developed within the Coordination Unit 7 (CU7) of the Gaia Data Processing and Analysis Consortium (DPAC), the coordination unit charged with the processing and analysis of variable sources observed by Gaia, to validate and fully characterise Cepheids and RR Lyrae stars observed by the spacecraft. We describe how the SOS for Cepheids and RR Lyrae stars (SOS Cep&RRL) was specifically tailored to analyse Gaia's G-band photometric time-series with a South Ecliptic Pole (SEP) footprint, which covers an external region of the Large Magellanic Cloud (LMC). G-band time-series photometry and characterization by the SOS Cep&RRL pipeline (mean magnitude and pulsation characteristics) are published in Gaia Data Release 1 (Gaia DR1) for a total sample of 3,194 variable stars, 599 Cepheids and 2,595 RR Lyrae stars, of which 386 (43 Cepheids and 343 RR Lyrae stars) are new discoveries by Gaia. All 3,194 stars are distributed over an area extending 38 degrees on either side from a point offset from the centre of the LMC by about 3 degrees to the north and 4 degrees to the east. The vast majority, but not all, are located within the LMC. The published sample also includes a few bright RR Lyrae stars that trace the outer halo of the Milky Way in front of the LMC.

This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster.

Stellar clusters are open windows to understand stellar evolution. Specifically, the change with time and the dependence on mass of different stellar properties. As such, they are our laboratories where different theories can be tested. We try to understand the origin of the connection between lithium depletion in F, G and K stars, rotation and activity, in particular in the Pleiades open cluster. We have collected all the relevant data in the literature, including information regarding rotation period, binarity and activity, and cross-matched with proper motions, multi-wavelength photometry and membership probability from the DANCe database. In order to avoid biases, only Pleiades single members with probabilities larger than p=0.75 have been included in the discussion. Results. The analysis confirms that there is a strong link between activity, rotation and the lithium equivalent width excess, specially for the range Lum(bol) = 0.5-0.2 Lsun (about K2-K7 spectral types or 0.75-0.95 Msun). It is not possible to disentangle these effects but we cannot exclude that the observed lithium overabundance is partially an observational effect due to enhanced activity, due to a large coverage by stellar spots induced by high rotation rates. Since a bona fide lithium enhancement is present in young, fast rotators, both activity and rotation should play a role in the lithium problem.

We present a high-precision proper motion study of 873 X-ray and spectroscopically selected stars in the massive OB association Cygnus OB2 as part of the DANCe project. These were calculated from images spanning a 15 year baseline and have typical precisions < 1 mas/yr. We calculate the velocity dispersion in the two axes to be $\sigma_\alpha(c) = 13.0^{+0.8}_{-0.7}$ and $\sigma_\delta(c) = 9.1^{+0.5}_{-0.5}$ km/s, using a 2-component, 2-dimensional model that takes into account the uncertainties on the measurements. This gives a 3-dimensional velocity dispersion of $\sigma_{3D} = 17.8 \pm 0.6$ km/s implying a virial mass significantly larger than the observed stellar mass, confirming that the association is gravitationally unbound. The association appears to be dynamically unevolved, as evidenced by considerable kinematic substructure, non-isotropic velocity dispersions and a lack of energy equipartition. The proper motions show no evidence for a global expansion pattern, with approximately the same amount of kinetic energy in expansion as there is in contraction, which argues against the association being an expanded star cluster disrupted by process such as residual gas expulsion or tidal heating. The kinematic substructures, which appear to be close to virial equilibrium and have typical masses of 40-400 M$_\odot$, also do not appear to have been affected by the expulsion of the residual gas. We conclude that Cyg OB2 was most likely born highly substructured and globally unbound, with the individual subgroups born in (or close to) virial equilibrium, and that the OB association has not experienced significant dynamical evolution since then.

One of the fundamental goals of modern Astronomy is to estimate the physical parameters of galaxies from images in different spectral bands. We present a hierarchical Bayesian model for obtaining age maps from images in the \Ha\ line (taken with Taurus Tunable Filter (TTF)), ultraviolet band (far UV or FUV, from GALEX) and infrared bands (24, 70 and 160 microns ($\mu$m), from Spitzer). As shown in S\'anchez-Gil et al. (2011), we present the burst ages for young stellar populations in the nearby and nearly face on galaxy M74. As it is shown in the previous work, the \Ha\ to FUV flux ratio gives a good relative indicator of very recent star formation history (SFH). As a nascent star-forming region evolves, the \Ha\ line emission declines earlier than the UV continuum, leading to a decrease in the \Ha\/FUV ratio. Through a specific star-forming galaxy model (Starburst 99, SB99), we can obtain the corresponding theoretical ratio \Ha\ / FUV to compare with our observed flux ratios, and thus to estimate the ages of the observed regions. Due to the nature of the problem, it is necessary to propose a model of high complexity to take into account the mean uncertainties, and the interrelationship between parameters when the \Ha\ / FUV flux ratio mentioned above is obtained. To address the complexity of the model, we propose a Bayesian hierarchical model, where a joint probability distribution is defined to determine the parameters (age, metallicity, IMF), from the observed data, in this case the observed flux ratios \Ha\ / FUV. The joint distribution of the parameters is described through an i.i.d. (independent and identically distributed random variables), generated through MCMC (Markov Chain Monte Carlo) techniques.

The DANCe survey provides photometric and astrometric (position and proper motion) measurements for approximately 2 millions unique sources in a region encompassing $\approx$80deg$^{2}$ centered around the Pleiades cluster. We aim at deriving a complete census of the Pleiades, and measure the mass and luminosity function of the cluster. Using the probabilistic selection method described in Sarro+2014, we identify high probability members in the DANCe ($i\ge$14mag) and Tycho-2 ($V\lesssim$12mag) catalogues, and study the properties of the cluster over the corresponding luminosity range. We find a total of 2109 high probability members, of which 812 are new, making it the most extensive and complete census of the cluster to date. The luminosity and mass functions of the cluster are computed from the most massive members down to $\approx$0.025M$_{\odot}$. The size, sensitivity and quality of the sample result in the most precise luminosity and mass functions observed to date for a cluster. Our census supersedes previous studies of the Pleiades cluster populations, both in terms of sensitivity and accuracy.

We present the variability processing and analysis that is foreseen for the Gaia mission within Coordination Unit 7 (CU7) of the Gaia Data Processing and Analysis Consortium (DPAC). A top level description of the tasks is given.

We aim at identifying the cluster's members by deriving membership probabilities for the sources within 1 degree of the cluster's center, going further away than equivalent previous studies. We measure accurate proper motions and multi-wavelength (optical and near-infrared) photometry using ground based archival images of the cluster. We use these measurements to compute membership probabilities. The list of candidate members from Barrado+2001 is used as training set to identify the cluster's locus in a multi-dimensional space made of proper motions, luminosities and colors. The final catalog includes 338892 sources with multi-wavelength photometry. Approximately half (194452) were detected at more than two epochs and we measured their proper motion and used it to derive membership probability. A total of 4349 candidate members with membership probabilities greater than 50% are found in this sample in the luminosity range between 10 and 22mag. The slow proper motion of the cluster and the overlap of its sequence with the field and background sequences in almost all color-magnitude and color-color diagrams complicate the analysis and the contamination level is expected to be significant. Our study nevertheless provides a coherent and quantitative membership analysis of Messier 35 based on a large fraction of the best ground-based data sets obtained over the past 18 years. As such, it represents a valuable input for follow-up studies using in particular the Kepler K2 photometric time series.

Following the recent discovery of a large population of young stars in front of the Orion Nebula, we carried out an observational campaign with the DECam wide-field camera covering ~10~deg^2 centered on NGC 1980 to confirm, probe the extent of, and characterize this foreground population of pre-main-sequence stars. We confirm the presence of a large foreground population towards the Orion A cloud. This population contains several distinct subgroups, including NGC1980 and NGC1981, and stretches across several degrees in front of the Orion A cloud. By comparing the location of their sequence in various color-magnitude diagrams with other clusters, we found a distance and an age of 380pc and 5~10Myr, in good agreement with previous estimates. Our final sample includes 2123 candidate members and is complete from below the hydrogen-burning limit to about 0.3Msun, where the data start to be limited by saturation. Extrapolating the mass function to the high masses, we estimate a total number of ~2600 members in the surveyed region. We confirm the presence of a rich, contiguous, and essentially coeval population of about 2600 foreground stars in front of the Orion A cloud, loosely clustered around NGC1980, NGC1981, and a new group in the foreground of the OMC-2/3. For the area of the cloud surveyed, this result implies that there are more young stars in the foreground population than young stars inside the cloud. Assuming a normal initial mass function, we estimate that between one to a few supernovae must have exploded in the foreground population in the past few million years, close to the surface of Orion A, which might be responsible, together with stellar winds, for the structure and star formation activity in these clouds. This long-overlooked foreground stellar population is of great significance, calling for a revision of the star formation history in this region of the Galaxy.

We present a new technique designed to take full advantage of the high dimensionality (photometric, astrometric, temporal) of the DANCe survey to derive self-consistent and robust membership probabilities of the Pleiades cluster. We aim at developing a methodology to infer membership probabilities to the Pleiades cluster from the DANCe multidimensional astro-photometric data set in a consistent way throughout the entire derivation. The determination of the membership probabilities has to be applicable to censored data and must incorporate the measurement uncertainties into the inference procedure. We use Bayes' theorem and a curvilinear forward model for the likelihood of the measurements of cluster members in the colour-magnitude space, to infer posterior membership probabilities. The distribution of the cluster members proper motions and the distribution of contaminants in the full multidimensional astro-photometric space is modelled with a mixture-of-Gaussians likelihood. We analyse several representation spaces composed of the proper motions plus a subset of the available magnitudes and colour indices. We select two prominent representation spaces composed of variables selected using feature relevance determination techniques based in Random Forests, and analyse the resulting samples of high probability candidates. We consistently find lists of high probability (p > 0.9975) candidates with $\approx$ 1000 sources, 4 to 5 times more than obtained in the most recent astro-photometric studies of the cluster. The methodology presented here is ready for application in data sets that include more dimensions, such as radial and/or rotational velocities, spectral indices and variability.

The Gaia satellite will survey the entire celestial sphere down to 20th magnitude, obtaining astrometry, photometry, and low resolution spectrophotometry on one billion astronomical sources, plus radial velocities for over one hundred million stars. Its main objective is to take a census of the stellar content of our Galaxy, with the goal of revealing its formation and evolution. Gaia's unique feature is the measurement of parallaxes and proper motions with hitherto unparalleled accuracy for many objects. As a survey, the physical properties of most of these objects are unknown. Here we describe the data analysis system put together by the Gaia consortium to classify these objects and to infer their astrophysical properties using the satellite's data. This system covers single stars, (unresolved) binary stars, quasars, and galaxies, all covering a wide parameter space. Multiple methods are used for many types of stars, producing multiple results for the end user according to different models and assumptions. Prior to its application to real Gaia data the accuracy of these methods cannot be assessed definitively. But as an example of the current performance, we can attain internal accuracies (RMS residuals) on F,G,K,M dwarfs and giants at G=15 (V=15-17) for a wide range of metallicites and interstellar extinctions of around 100K in effective temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]), and 0.25dex in surface gravity (logg). The accuracy is a strong function of the parameters themselves, varying by a factor of more than two up or down over this parameter range. After its launch in November 2013, Gaia will nominally observe for five years, during which the system we describe will continue to evolve in light of experience with the real data.

We present an automated classification of stars exhibiting periodic, non-periodic and irregular light variations. The Hipparcos catalogue of unsolved variables is employed to complement the training set of periodic variables of Dubath et al. with irregular and non-periodic representatives, leading to 3881 sources in total which describe 24 variability types. The attributes employed to characterize light-curve features are selected according to their relevance for classification. Classifier models are produced with random forests and a multistage methodology based on Bayesian networks, achieving overall misclassification rates under 12 per cent. Both classifiers are applied to predict variability types for 6051 Hipparcos variables associated with uncertain or missing types in the literature.

We aimed to assess the accuracy of the Gaia teff and logg estimates as derived with current models and observations. We assessed the validity of several inference techniques for deriving the physical parameters of ultra-cool dwarf stars. We used synthetic spectra derived from ultra-cool dwarf models to construct (train) the regression models. We derived the intrinsic uncertainties of the best inference models and assessed their validity by comparing the estimated parameters with the values derived in the bibliography for a sample of ultra-cool dwarf stars observed from the ground. We estimated the total number of ultra-cool dwarfs per spectral subtype, and obtained values that can be summarised (in orders of magnitude) as 400000 objects in the M5-L0 range, 600 objects between L0 and L5, 30 objects between L5 and T0, and 10 objects between T0 and T8. A bright ultra-cool dwarf (with teff=2500 K and \logg=3.5 will be detected by Gaia out to approximately 220 pc, while for teff=1500 K (spectral type L5) and the same surface gravity, this maximum distance reduces to 10-20 pc. The RMSE of the prediction deduced from ground-based spectra of ultra-cool dwarfs simulated at the Gaia spectral range and resolution, and for a Gaia magnitude G=20 is 213 K and 266 K for the models based on k-nearest neighbours and Gaussian process regression, respectively. These are total errors in the sense that they include the internal and external errors, with the latter caused by the inability of the synthetic spectral models (used for the construction of the regression models) to exactly reproduce the observed spectra, and by the large uncertainties in the current calibrations of spectral types and effective temperatures.

A complete periodic star extraction and classification scheme is set up and tested with the Hipparcos catalogue. The efficiency of each step is derived by comparing the results with prior knowledge coming from the catalogue or from the literature. A combination of two variability criteria is applied in the first step to select 17 006 variability candidates from a complete sample of 115 152 stars. Our candidate sample turns out to include 10 406 known variables (i.e., 90% of the total of 11 597) and 6600 contaminating constant stars. A random forest classification is used in the second step to extract 1881 (82%) of the known periodic objects while removing entirely constant stars from the sample and limiting the contamination of non-periodic variables to 152 stars (7.5%). The confusion introduced by these 152 non-periodic variables is evaluated in the third step using the results of the Hipparcos periodic star classification presented in a previous study (Dubath et al. [1]).

We present an evaluation of the performance of an automated classification of the Hipparcos periodic variable stars into 26 types. The sub-sample with the most reliable variability types available in the literature is used to train supervised algorithms to characterize the type dependencies on a number of attributes. The most useful attributes evaluated with the random forest methodology include, in decreasing order of importance, the period, the amplitude, the V-I colour index, the absolute magnitude, the residual around the folded light-curve model, the magnitude distribution skewness and the amplitude of the second harmonic of the Fourier series model relative to that of the fundamental frequency. Random forests and a multi-stage scheme involving Bayesian network and Gaussian mixture methods lead to statistically equivalent results. In standard 10-fold cross-validation experiments, the rate of correct classification is between 90 and 100%, depending on the variability type. The main mis-classification cases, up to a rate of about 10%, arise due to confusion between SPB and ACV blue variables and between eclipsing binaries, ellipsoidal variables and other variability types. Our training set and the predicted types for the other Hipparcos periodic stars are available online.

We present a novel automated methodology to detect and classify periodic variable stars in a large database of photometric time series. The methods are based on multivariate Bayesian statistics and use a multi-stage approach. We applied our method to the ground-based data of the TrES Lyr1 field, which is also observed by the Kepler satellite, covering ~26000 stars. We found many eclipsing binaries as well as classical non-radial pulsators, such as slowly pulsating B stars, Gamma Doradus, Beta Cephei and Delta Scuti stars. Also a few classical radial pulsators were found.

Solar explosive events are commonly explained as small scale magnetic reconnection events, although unambiguous confirmation of this scenario remains elusive due to the lack of spatial resolution and of the statistical analysis of large enough samples of this type of events. In this work, we propose a sound statistical treatment of data cubes consisting of a temporal sequence of long slit spectra of the solar atmosphere. The analysis comprises all the stages from the explosive event detection to its characterization and the subsequent sample study. We have designed two complementary approaches based on the combination of standard statistical techniques (Robust Principal Component Analysis in one approach and wavelet decomposition and Independent Component Analysis in the second) in order to obtain least biased samples. These techniques are implemented in the spirit of letting the data speak for themselves. The analysis is carried out for two spectral lines: the C IV line at 1548.2 angstroms and the Ne VIII line at 770.4 angstroms. We find significant differences between the characteristics of the line profiles emitted in the proximities of two active regions, and in the quiet Sun, most visible in the relative importance of a separate population of red shifted profiles. We also find a higher frequency of explosive events near the active regions, and in the C IV line. The distribution of the explosive events characteristics is interpreted in the light of recent numerical simulations. Finally, we point out several regions of the parameter space where the reconnection model has to be refined in order to explain the observations.

The CoRoT (Convection, Rotation and planetary Transits) space mission provides a valuable opportunity to monitor stars with uninterrupted time sampling for up to 150 days at a time. The study of RR Lyrae stars, performed in the framework of the Additional Programmes belonging to the exoplanetary field, will particularly benefit from such dense, long-duration monitoring. The Blazhko effect in RR Lyrae stars is a long-standing, unsolved problem of stellar astrophysics. We used the CoRoT data of the new RR Lyrae variable CoRoT 101128793 (f0=2.119 c/d, P=0.4719296 d) to provide us with more detailed observational facts to understand the physical process behind the phenomenon. The CoRoT data were corrected for one jump and the long-term drift. We applied different period-finding techniques to the corrected timeseries to investigate amplitude and phase modulation. We detected 79 frequencies in the light curve of CoRoT 101128793. They have been identified as the main frequency f0, and its harmonics, two independent terms, the terms related to the Blazhko frequency, and several combination terms. A Blazhko frequency fB=0.056 c/d and a triplet structure around the fundamental radial mode and harmonics were detected, as well as a long-term variability of the Blazhko modulation. Indeed, the amplitude of the main oscillation is decreasing along the CoRoT survey. The Blazhko modulation is one of the smallest observed in RR Lyrae stars. Moreover, the additional modes f1=3.630 and f2=3.159 c/d are detected. Taking its ratio with the fundamental radial mode into account, the term f1 could be the identified as the second radial overtone. Detecting of these modes in horizontal branch stars is a new result obtained by CoRoT.

* Aims. We describe here the main functionalities of the LAEX (Laboratorio de Astrofisica Estelar y Exoplanetas/Laboratory for Stellar Astrophysics and Exoplanets) and NASA portals for CoRoT Public Data. The CoRoT archive at LAEX was opened to the community in January 2009 and is managed in the framework of the Spanish Virtual Observatory. NStED (NASA Star and Exoplanet Database) serves as the CoRoT portal for the US astronomical community. NStED is a general purpose stellar and exoplanet archive with the aim of providing support for NASA planet finding and characterisation goals, and the planning and support of NASA and other space missions. CoRoT data at LAEX and NStED can be accessed at http://sdc.laeff.inta.es/corotfa/ and http://nsted.ipac.caltech.edu,respectively. * Methods. Based on considerable experience with astronomical archives, the aforementioned archives are designed with the aim of delivering science-quality data in a simple and efficient way. * Results. LAEX and NStED not only provide access to CoRoT Public Data but furthermore serve a variety of observed and calculated astrophysical data. In particular, NStED provides scientifically validated information on stellar and planetary data related to the search for and characterization of extrasolar planets, and LAEX makes any information from Virtual Observatory services available to the astronomical community.

Context. Discovery of new variability classes in large surveys using multivariate statistics techniques such as clustering, relies heavily on the correct understanding of the distribution of known classes as point processes in parameter space. Aims. Our objective is to analyze the correspondence between the classical stellar variability types and the clusters found in the distribution of light curve parameters and colour indices of stars in the CoRoT exoplanet sample. The final aim is to help in the identification on new types of variability by first identifying the well known variables in the CoRoT sample. Methods. We apply unsupervised classification algorithms to identify clusters of variable stars from modes of the probability density distribution. We use reference variability databases (Hipparcos and OGLE) as a framework to calibrate the clustering methodology. Furthermore, we use the results from supervised classification methods to interpret the resulting clusters. Results.We interpret the clusters in the Hipparcos and OGLE LMC databases in terms of large-amplitude radial pulsators in the classical instability strip and of various types of eclipsing binaries. The Hipparcos data also provide clear distributions for low-amplitude nonradial pulsators. We show that the preselection of targets for the CoRoT exoplanet programme results in a completely different probability density landscape than the OGLE data, the interpretation of which involves mainly classes of low-amplitude variability in main-sequence stars. Our findings will be incorporated to improve the supervised classification used in the CoRoT catalogue production, once the existence of new classes or subtypes will be confirmed from complementary spectroscopic observations.

The classification of time series from photometric large scale surveys into variability types and the description of their properties is difficult for various reasons including but not limited to the irregular sampling, the usually few available photometric bands, and the diversity of variable objects. Furthermore, it can be seen that different physical processes may sometimes produce similar behavior which may end up to be represented as same models. In this article we will also be presenting our approach for processing the data resulting from the Gaia space mission. The approach may be classified into following three broader categories: supervised classification, unsupervised classifications, and "so-called" extractor methods i.e. algorithms that are specialized for particular type of sources. The whole process of classification- from classification attribute extraction to actual classification- is done in an automated manner.

A campaign on the open cluster NGC 5617 was organized in order to characterize the pulsations and to better understand the internal structure of its stars. The variability of the cluster members was never studied before. We present the observations taken and an up-to-date analysis of the obtained time series, especially of several SPB candidates we discovered.

We aim to extend and test the classifiers presented in a previous work against an independent dataset. We complement the assessment of the validity of the classifiers by applying them to the set of OGLE light curves treated as variable objects of unknown class. The results are compared to published classification results based on the so-called extractor methods.Two complementary analyses are carried out in parallel. In both cases, the original time series of OGLE observations of the Galactic bulge and Magellanic Clouds are processed in order to identify and characterize the frequency components. In the first approach, the classifiers are applied to the data and the results analyzed in terms of systematic errors and differences between the definition samples in the training set and in the extractor rules. In the second approach, the original classifiers are extended with colour information and, again, applied to OGLE light curves. We have constructed a classification system that can process huge amounts of time series in negligible time and provide reliable samples of the main variability classes. We have evaluated its strengths and weaknesses and provide potential users of the classifier with a detailed description of its characteristics to aid in the interpretation of classification results. Finally, we apply the classifiers to obtain object samples of classes not previously studied in the OGLE database and analyse the results. We pay specific attention to the B-stars in the samples, as their pulsations are strongly dependent on metallicity.

The fast classification of new variable stars is an important step in making them available for further research. Selection of science targets from large databases is much more efficient if they have been classified first. Defining the classes in terms of physical parameters is also important to get an unbiased statistical view on the variability mechanisms and the borders of instability strips. Our goal is twofold: provide an overview of the stellar variability classes that are presently known, in terms of some relevant stellar parameters; use the class descriptions obtained as the basis for an automated `supervised classification' of large databases. Such automated classification will compare and assign new objects to a set of pre-defined variability training classes. For every variability class, a literature search was performed to find as many well-known member stars as possible, or a considerable subset if too many were present. Next, we searched on-line and private databases for their light curves in the visible band and performed period analysis and harmonic fitting. The derived light curve parameters are used to describe the classes and define the training classifiers. We compared the performance of different classifiers in terms of percentage of correct identification, of confusion among classes and of computation time. We describe how well the classes can be separated using the proposed set of parameters and how future improvements can be made, based on new large databases such as the light curves to be assembled by the CoRoT and Kepler space missions.

The Variable Star One-shot Project (VSOP) is aimed at (1) providing the variability type and spectral type of all unstudied variable stars, (2) process, publish, and make the data available as automatically as possible, and (3) generate serendipitous discoveries. This first paper describes the project itself, the acquisition of the data, the dataflow, the spectroscopic analysis and the on-line availability of the fully calibrated and reduced data. We also present the results on the 221 stars observed during the first semester of the project. We used the high-resolution echelle spectrographs HARPS and FEROS in the ESO La Silla Observatory (Chile) to survey known variable stars. Once reduced by the dedicated pipelines, the radial velocities are determined from cross correlation with synthetic template spectra, and the spectral types are determined by an automatic minimum distance matching to synthetic spectra, with traditional manual spectral typing cross-checks. The variability types are determined by manually evaluating the available light curves and the spectroscopy. In the future, a new automatic classifier, currently being developed by members of the VSOP team, based on these spectroscopic data and on the photometric classifier developed for the COROT and Gaia space missions, will be used. We confirm or revise spectral types of 221 variable stars from the GCVS. We identify 26 previously unknown multiple systems, among them several visual binaries with spectroscopic binary individual components. We present new individual results for the multiple systems V349 Vel and BC Gru, for the composite spectrum star V4385 Sgr, for the T-Tauri star V1045 Sco, and for DM Boo which we re-classify as a BY Draconis variable. The complete data release can be accessed via the VSOP web site.

The Variable Star One-shot Project (VSOP) aimed at providing to the world-wide stellar community the necessary one-shot spectrum of unstudied variable stars, too often classified as such by an analysis of photometric data only. The VSOP has established an new kind of observational model, where all steps from observations to spectral analysis, are automatized (or are underway to be fully automatized). The project is centralized on a collaborative wiki website. The VSOP operational model is very successful, data is continously flowing and being analyszed, and VSOP is now a worldwide open collaboration of people with very different and complementary skills and expertise. The idea of a central wiki website has been extended by one of us to propose a new service to the whole astronomical community, called Wikimbad. Wikimbad is an open wiki website aimed at collecting, organizing and making publicly available all kind of reduced and published astronomical data. Its strengths and a comparison with the Virtual Observatory are discussed. See: http://vsop.sc.eso.org and http://wikimbad.org

In this work we present a system for the automatic classification of the light curves of eclipsing binaries. This system is based on a classification scheme that aims to separate eclipsing binary sistems according to their geometrical configuration in a modified version of the traditional classification scheme. The classification is performed by a Bayesian ensemble of neural networks trained with {\em Hipparcos} data of seven different categories including eccentric binary systems and two types of pulsating light curve morphologies.