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.

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.

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.

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.

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.

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.

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.

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

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.