Astrometric detection involves a precise measurement of stellar positions, and is widely regarded as the leading concept presently ready to find earth-mass planets in temperate orbits around nearby sun-like stars. The TOLIMAN space telescope[39] is a low-cost, agile mission concept dedicated to narrow-angle astrometric monitoring of bright binary stars. In particular the mission will be optimised to search for habitable-zone planets around Alpha Centauri AB. If the separation between these two stars can be monitored with sufficient precision, tiny perturbations due to the gravitational tug from an unseen planet can be witnessed and, given the configuration of the optical system, the scale of the shifts in the image plane are about one millionth of a pixel. Image registration at this level of precision has never been demonstrated (to our knowledge) in any setting within science. In this paper we demonstrate that a Deep Convolutional Auto-Encoder is able to retrieve such a signal from simplified simulations of the TOLIMAN data and we present the full experimental pipeline to recreate out experiments from the simulations to the signal analysis. In future works, all the more realistic sources of noise and systematic effects present in the real-world system will be injected into the simulations.

The Commission on Science and Information Technology (CTCI) of the Brazilian Astronomical Society (SAB) is tasked with assisting the Society on issues of astronomical data management, from its handling and the management of data centres and networks, to technical aspects of the archiving, storage and dissemination of data. In this paper we present a summary of the results of a survey recently conducted by the Commission to diagnose the status of several data-related issues within the Brazilian astronomical community, as well as some proposals derived therefrom.

The large astrometric and photometric survey performed by the Gaia mission allows for a panoptic view of the Galactic disc and in its stellar cluster population. Hundreds of clusters were only discovered after the latest G data release (DR2) and have yet to be characterised. Here we make use of the deep and homogeneous Gaia photometry down to G=18 to estimate the distance, age, and interstellar reddening for about 2000 clusters identified with Gaia~DR2 astrometry. We use these objects to study the structure and evolution of the Galactic disc. We rely on a set of objects with well-determined parameters in the literature to train an artificial neural network to estimate parameters from the Gaia photometry of cluster members and their mean parallax. We obtain reliable parameters for 1867 clusters. Our new homogeneous catalogue confirms the relative lack of old clusters in the inner disc (with a few notable exceptions). We also quantify and discuss the variation of scale height with cluster age, and detect the Galactic warp in the distribution of old clusters. This work results in a large and homogenous cluster catalogue. However, the present sample is still unable to trace the Outer spiral arm of the Milky Way, which indicates that the outer disc cluster census might still be incomplete.

Cosmic voids and their corresponding redshift-aggregated projections of mass densities, known as troughs, play an important role in our attempt to model the large-scale structure of the Universe. Understanding these structures leads to tests comparing the standard model with alternative cosmologies, constraints on the dark energy equation of state, and provides evidence to differentiate among gravitational theories. In this paper, we extend the subspace-constrained mean shift algorithm, a recently introduced method to estimate density ridges, and apply it to 2D weak-lensing mass density maps from the Dark Energy Survey Y1 data release to identify curvilinear filamentary structures. We compare the obtained ridges with previous approaches to extract trough structure in the same data, and apply curvelets as an alternative wavelet-based method to constrain densities. We then invoke the Wasserstein distance between noisy and noiseless simulations to validate the denoising capabilities of our method. Our results demonstrate the viability of ridge estimation as a precursor for denoising weak lensing quantities to recover the large-scale structure, paving the way for a more versatile and effective search for troughs.

The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

Open clusters (OCs) are popular tracers of the structure and evolutionary history of the Galactic disk. The OC population is often considered to be complete within 1.8 kpc of the Sun. The recent Gaia Data Release 2 (DR2) allows the latter claim to be challenged. We perform a systematic search for new OCs in the direction of Perseus using precise and accurate astrometry from Gaia DR2. We implement a coarse-to-fine search method. First, we exploit spatial proximity using a fast density-aware partitioning of the sky via a k-d tree in the spatial domain of Galactic coordinates, (l, b). Secondly, we employ a Gaussian mixture model in the proper motion space to quickly tag fields around OC candidates. Thirdly, we apply an unsupervised membership assignment method, UPMASK, to scrutinise the candidates. We visually inspect colour-magnitude diagrams to validate the detected objects. Finally, we perform a diagnostic to quantify the significance of each identified overdensity in proper motion and in parallax space We report the discovery of 41 new stellar clusters. This represents an increment of at least 20% of the previously known OC population in this volume of the Milky Way. We also report on the clear identification of NGC 886, an object previously considered an asterism. This letter challenges the previous claim of a near-complete sample of open clusters up to 1.8 kpc. Our results reveal that this claim requires revision, and a complete census of nearby open clusters is yet to be found.

Context. Open clusters are very good tracers of the evolution of the Galactic disc. Thanks to Gaia, their kinematics can be investigated with an unprecedented precision and accuracy. Aims. The distribution of open clusters in the 6D phase space is revisited with Gaia DR2. Methods. The weighted mean radial velocity of open clusters was determined, using the most probable members available from a previous astrometric investigation that also provided mean parallaxes and proper motions. Those parameters, all derived from Gaia DR2 only, were combined to provide the 6D phase space information of 861 clusters. The velocity distribution of nearby clusters was investigated, as well as the spatial and velocity distributions of the whole sample as a function of age. A high quality subsample was used to investigate some possible pairs and groups of clusters sharing the same Galactic position and velocity. Results. For the high quality sample that has 406 clusters, the median uncertainty of the weighted mean radial velocity is 0.5 km/s. The accuracy, assessed by comparison to ground-based high resolution spectroscopy, is better than 1 km/s. Open clusters nicely follow the velocity distribution of field stars in the close Solar neighbourhood previously revealed by Gaia DR2. As expected, the vertical distribution of young clusters is very flat but the novelty is the high precision to which this can be seen. The dispersion of vertical velocities of young clusters is at the level of 5 km/s. Clusters older than 1 Gyr span distances to the Galactic plane up to 1 kpc with a vertical velocity dispersion of 14 km/s, typical of the thin disc. Five pairs of clusters and one group with five members are possibly physically related. Other binary candidates previously identified turn out to be chance alignment.

Aims: In this work, we aim to provide a reliable list of gravitational lens (GL) candidates based on a search performed over the entire Gaia Data Release 2 (Gaia DR2). We also show that the sole astrometric and photometric informations coming from the Gaia satellite yield sufficient insights for supervised learning methods to automatically identify GL candidates with an efficiency that is comparable to methods based on image processing. Methods: We simulated 106,623,188 lens systems composed of more than two images, based on a regular grid of parameters characterizing a non-singular isothermal ellipsoid lens model in the presence of an external shear. These simulations are used as an input for training and testing our supervised learning models consisting of Extremely Randomized Trees. The latter are finally used to assign to each of the 2,129,659 clusters of celestial objects a discriminant value that reflects the ability of our simulations to match the observed relative positions and fluxes from each cluster. Once complemented with additional constraints, these discriminant values allowed us to identify GL candidates out of the list of clusters. Results: We report the discovery of 15 new quadruply-imaged lens candidates with angular separations less than 6" and assess the performance of our approach by recovering 12 out of the 13 known quadruply-imaged systems with all their components detected in Gaia DR2 with a misclassification rate of fortuitous clusters of stars as lens systems that is below one percent. Similarly, the identification capability of our method regarding quadruply-imaged systems where three images are detected in Gaia DR2 is assessed by recovering 10 out of the 13 known quadruply-imaged systems having one of their constituting images discarded. The associated misclassification rate varying then between 5.8% and 20%, depending on the image we decided to remove.

We report a framework for spectroscopic follow-up design for optimizing supernova photometric classification. The strategy accounts for the unavoidable mismatch between spectroscopic and photometric samples, and can be used even in the beginning of a new survey -- without any initial training set. The framework falls under the umbrella of active learning (AL), a class of algorithms that aims to minimize labelling costs by identifying a few, carefully chosen, objects which have high potential in improving the classifier predictions. As a proof of concept, we use the simulated data released after the Supernova Photometric Classification Challenge (SNPCC) and a random forest classifier. Our results show that, using only 12\% the number of training objects in the SNPCC spectroscopic sample, this approach is able to double purity results. Moreover, in order to take into account multiple spectroscopic observations in the same night, we propose a semi-supervised batch-mode AL algorithm which selects a set of $N=5$ most informative objects at each night. In comparison with the initial state using the traditional approach, our method achieves 2.3 times higher purity and comparable figure of merit results after only 180 days of observation, or 800 queries (73% of the SNPCC spectroscopic sample size). Such results were obtained using the same amount of spectroscopic time necessary to observe the original SNPCC spectroscopic sample, showing that this type of strategy is feasible with current available spectroscopic resources. The code used in this work is available in the COINtoolbox: https://github.com/COINtoolbox/ActSNClass .

Machine learning techniques offer a precious tool box for use within astronomy to solve problems involving so-called big data. They provide a means to make accurate predictions about a particular system without prior knowledge of the underlying physical processes of the data. In this article, and the companion papers of this series, we present the set of Generalized Linear Models (GLMs) as a fast alternative method for tackling general astronomical problems, including the ones related to the machine learning paradigm. To demonstrate the applicability of GLMs to inherently positive and continuous physical observables, we explore their use in estimating the photometric redshifts of galaxies from their multi-wavelength photometry. Using the gamma family with a log link function we predict redshifts from the PHoto-z Accuracy Testing simulated catalogue and a subset of the Sloan Digital Sky Survey from Data Release 10. We obtain fits that result in catastrophic outlier rates as low as ~1% for simulated and ~2% for real data. Moreover, we can easily obtain such levels of precision within a matter of seconds on a normal desktop computer and with training sets that contain merely thousands of galaxies. Our software is made publicly available as an user-friendly package developed in Python, R and via an interactive web application (https://cosmostatisticsinitiative.shinyapps.io/CosmoPhotoz). This software allows users to apply a set of GLMs to their own photometric catalogues and generates publication quality plots with minimum effort from the user. By facilitating their ease of use to the astronomical community, this paper series aims to make GLMs widely known and to encourage their implementation in future large-scale projects, such as the Large Synoptic Survey Telescope.

Astronomical observations of extended sources, such as cubes of integral field spectroscopy (IFS), encode auto-correlated spatial structures that cannot be optimally exploited by standard methodologies. This work introduces a novel technique to model IFS datasets, which treats the observed galaxy properties as realizations of an unobserved Gaussian Markov random field. The method is computationally efficient, resilient to the presence of low-signal-to-noise regions, and uses an alternative to Markov Chain Monte Carlo for fast Bayesian inference, the Integrated Nested Laplace Approximation (INLA). As a case study, we analyse 721 IFS data cubes of nearby galaxies from the CALIFA and PISCO surveys, for which we retrieve the maps of the following physical properties: age, metallicity, mass and extinction. The proposed Bayesian approach, built on a generative representation of the galaxy properties, enables the creation of synthetic images, recovery of areas with bad pixels, and an increased power to detect structures in datasets subject to substantial noise and/or sparsity of sampling. A snippet code to reproduce the analysis of this paper is available in the COIN toolbox, together with the field reconstructions of the CALIFA and PISCO samples.

We report the spectroscopic confirmation and modeling of the quadruply imaged quasar GRAL113100-441959, the first gravitational lens (GL) to be discovered mainly from astrometric considerations. Follow-up spectra obtained with Keck/LRIS reveal the lensing nature of this quadruply-imaged quasar with redshift $z_s = 1.090 \pm 0.002$, but show no evidence of the central lens galaxy. Using the image positions and $G$-band flux ratios provided by Gaia Data Release 2 as constraints, we model the system with a singular power-law elliptical mass distribution (SPEMD) plus external shear, to different levels of complexity. We show that relaxing the isothermal constraint of the SPEMD is not statistically significant, and thus we simplify the SPEMD to a singular isothermal ellipsoid to estimate the Einstein radius of the main lens galaxy $\theta_{\text{E}} = 0.\!\!^{\prime\prime}851$, the intensity and position angle of the external shear $(\gamma,\theta_{\gamma}) = (0.044,11.\!\!^{\circ}5)$, and we predict the lensing galaxy position to be $(x_{\text{gal}},y_{\text{gal}}) = (-0.\!\!^{\prime\prime}424,-0.\!\!^{\prime\prime}744)$ with respect to image A. We provide time delay predictions for pairs of images, assuming a plausible range of lens redshift values $z_l$ between $0.5$ and $0.9$. We finally examine the impact on time delays of the so-called Source Position Transformation, a family of degeneracies existing between different lens density profiles that reproduce most of the lensing observables equally well. We show that this effect contributes significantly to the time delay error budget and cannot be ignored during the modeling. This has implications for robust cosmography applications of lensed systems. GRAL113100-441959 is the first in a series of seven new spectroscopically confirmed GLs discovered from Gaia Data Release 2.

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.

Open clusters are convenient probes of the structure and history of the Galactic disk. They are also fundamental to stellar evolution studies. The second Gaia data release contains precise astrometry at the sub-milliarcsecond level and homogeneous photometry at the mmag level, that can be used to characterise a large number of clusters over the entire sky. In this study we aim to a establish list of members and derive mean parameters, in particular distances, for as many clusters as possible, making use of Gaia data alone. We compile a list of thousands of known or putative clusters from the literature. We then apply an unsupervised membership assignment code, UPMASK, to the Gaia DR2 data contained within the fields of those clusters. We obtained a list of members and cluster parameters for 1229 clusters. As expected, the youngest clusters are seen to be tightly distributed near the Galactic plane and to trace the spiral arms of the Milky Way, while older objects are more uniformly distributed, deviate further from the plane, and tend to be located at larger Galactocentric distances. Thanks to the quality of GaiaDR2 astrometry, the fully homogeneous parameters derived in this study are the most precise to date. Furthermore, we report on the serendipitous discovery of 60 new open clusters in the fields analysed during this study.

Thanks to its spatial resolution the ESA/Gaia space mission offers a unique opportunity to discover new multiply-imaged quasars and to study the already known lensed systems at sub-milliarcsecond astrometric precisions. In this paper, we address the detection of the known multiply-imaged quasars from the Gaia Data Release 2 and determine the astrometric and photometric properties of the individually detected images found in the Gaia DR2 catalogue. We have compiled an exhaustive list of quasar gravitational lenses from the literature to search for counterparts in the Gaia Data Release 2. We then analyze the astrometric and photometric properties of these Gaia's detections. To highlight the tremendous potential of Gaia at the sub-milliarcsecond level we finally perform a simple Bayesian modeling of the well-known gravitational lens system HE0435-1223, using Gaia Data Release 2 and HST astrometry. From 478 known multiply imaged quasars, 200 have at least one image found in the Gaia Data Release 2. Among the 41 known quadruply-imaged quasars of the list, 26 have at least one image in the Gaia Data Release 2, 12 of which are fully detected, 6 have three counterparts, 7 have two and 1 has only one. As expected, the modeling of HE0435-1223 shows that the model parameters are significantly better constrained when using Gaia astrometry compared to HST astrometry, in particular the relative positions of the background quasar source and the centroid of the deflector. The Gaia sub-milliarcsecond astrometry also significantly reduces the parameter correlations. Besides providing an up-to-date list of multiply imaged quasars and their detection in the Gaia DR2, this paper shows that more complex modeling scenarios will certainly benefit from Gaia sub-milliarcsecond astrometry.

Context. Strong gravitationally lensed quasars are among the most interesting and useful observable extragalactic phenomena. Because their study constitutes a unique tool in various fields of astronomy, they are highly sought, not without difficulty. Indeed, even in this era of all-sky surveys, their recognition remains a great challenge, with barely a few hundred currently known systems. Aims. In this work we aim to detect new strongly lensed quasar candidates in the recently published Gaia Data Release 2 (DR2), which is the highest spatial resolution astrometric and photometric all-sky survey, attaining effective resolutions from 0.4" to 2.2". Methods. We cross-matched a merged list of quasars and candidates with the Gaia DR2 and found 1,839,143 counterparts within 0.5". We then searched matches with more than two Gaia DR2 counterparts within 6". We further narrowed the resulting list using astrometry and photometry compatibility criteria between the Gaia DR2 counterparts. A supervised machine learning method, Extremely Randomized Trees, is finally adopted to assign to each remaining system a probability of being lensed. Results. We report the discovery of three quadruply-imaged quasar candidates that are fully detected in Gaia DR2. These are the most promising new quasar lens candidates from Gaia DR2 and a simple singular isothermal ellipsoid lens model is able to reproduce their image positions to within $\sim$1 mas. This letter demonstrates the gravitational lens discovery potential of Gaia.

The Tycho-Gaia Astrometric Solution (TGAS) subset of the first Gaia catalogue contains an unprecedented sample of proper motions and parallaxes for two million stars brighter than G~12 mag. We take advantage of the full astrometric solution available for those stars to identify the members of known open clusters and compute mean cluster parameters using either TGAS or UCAC4 proper motions, and TGAS parallaxes. We apply an unsupervised membership assignment procedure to select high probability cluster members, we use a Bayesian/MCMC technique to fit stellar isochrones to the observed 2MASS JHK$_s$ magnitudes of the member stars and derive cluster parameters (age, metallicity, extinction, distance modulus), and we combine TGAS data with spectroscopic radial velocities to compute full Galactic orbits. We obtain mean astrometric parameters (proper motions and parallaxes) for 128 clusters closer than about 2 kpc, and cluster parameters from isochrone fitting for 26 of them located inside 1 kpc from the Sun. We show the orbital parameters obtained from integrating 36 orbits in a Galactic potential.

We invoke a Gaussian mixture model (GMM) to jointly analyse two traditional emission-line classification schemes of galaxy ionization sources: the Baldwin-Phillips-Terlevich (BPT) and $\rm W_{H\alpha}$ vs. [NII]/H$\alpha$ (WHAN) diagrams, using spectroscopic data from the Sloan Digital Sky Survey Data Release 7 and SEAGal/STARLIGHT datasets. We apply a GMM to empirically define classes of galaxies in a three-dimensional space spanned by the $\log$ [OIII]/H$\beta$, $\log$ [NII]/H$\alpha$, and $\log$ EW(H${\alpha}$), optical parameters. The best-fit GMM based on several statistical criteria suggests a solution around four Gaussian components (GCs), which are capable to explain up to 97 per cent of the data variance. Using elements of information theory, we compare each GC to their respective astronomical counterpart. GC1 and GC4 are associated with star-forming galaxies, suggesting the need to define a new starburst subgroup. GC2 is associated with BPT's Active Galaxy Nuclei (AGN) class and WHAN's weak AGN class. GC3 is associated with BPT's composite class and WHAN's strong AGN class. Conversely, there is no statistical evidence -- based on four GCs -- for the existence of a Seyfert/LINER dichotomy in our sample. Notwithstanding, the inclusion of an additional GC5 unravels it. The GC5 appears associated to the LINER and Passive galaxies on the BPT and WHAN diagrams respectively. Subtleties aside, we demonstrate the potential of our methodology to recover/unravel different objects inside the wilderness of astronomical datasets, without lacking the ability to convey physically interpretable results. The probabilistic classifications from the GMM analysis are publicly available within the COINtoolbox (https://cointoolbox.github.io/GMM\_Catalogue/).

Context: The first Gaia data release (DR1) delivered a catalogue of astrometry and photometry for over a billion astronomical sources. Within the panoply of methods used for data exploration, visualisation is often the starting point and even the guiding reference for scientific thought. However, this is a volume of data that cannot be efficiently explored using traditional tools, techniques, and habits. Aims: We aim to provide a global visual exploration service for the Gaia archive, something that is not possible out of the box for most people. The service has two main goals. The first is to provide a software platform for interactive visual exploration of the archive contents, using common personal computers and mobile devices available to most users. The second aim is to produce intelligible and appealing visual representations of the enormous information content of the archive. Methods: The interactive exploration service follows a client-server design. The server runs close to the data, at the archive, and is responsible for hiding as far as possible the complexity and volume of the Gaia data from the client. This is achieved by serving visual detail on demand. Levels of detail are pre-computed using data aggregation and subsampling techniques. For DR1, the client is a web application that provides an interactive multi-panel visualisation workspace as well as a graphical user interface. Results: The Gaia archive Visualisation Service offers a web-based multi-panel interactive visualisation desktop in a browser tab. It currently provides highly configurable 1D histograms and 2D scatter plots of Gaia DR1 and the Tycho-Gaia Astrometric Solution (TGAS) with linked views. An innovative feature is the creation of ADQL queries from visually defined regions in plots. [abridged]

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

Two of the main problems encountered in the development and accurate validation of photometric redshift (photo-z) techniques are the lack of spectroscopic coverage in feature space (e.g. colours and magnitudes) and the mismatch between photometric error distributions associated with the spectroscopic and photometric samples. Although these issues are well known, there is currently no standard benchmark allowing a quantitative analysis of their impact on the final photo-z estimation. In this work, we present two galaxy catalogues, Teddy and Happy, built to enable a more demanding and realistic test of photo-z methods. Using photometry from the Sloan Digital Sky Survey and spectroscopy from a collection of sources, we constructed datasets which mimic the biases between the underlying probability distribution of the real spectroscopic and photometric sample. We demonstrate the potential of these catalogues by submitting them to the scrutiny of different photo-z methods, including machine learning (ML) and template fitting approaches. Beyond the expected bad results from most ML algorithms for cases with missing coverage in feature space, we were able to recognize the superiority of global models in the same situation and the general failure across all types of methods when incomplete coverage is convoluted with the presence of photometric errors - a data situation which photo-z methods were not trained to deal with up to now and which must be addressed by future large scale surveys. Our catalogues represent the first controlled environment allowing a straightforward implementation of such tests. The data are publicly available within the COINtoolbox (https://github.com/COINtoolbox/photoz_catalogues).

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

Context. This paper presents an overview of the photometric data that are part of the first Gaia data release. Aims. The principles of the processing and the main characteristics of the Gaia photometric data are presented. Methods. The calibration strategy is outlined briefly and the main properties of the resulting photometry are presented. Results. Relations with other broadband photometric systems are provided. The overall precision for the Gaia photometry is shown to be at the milli-magnitude level and has a clear potential to improve further in future releases.

We aim at performing a kinematic census of young stellar objects (YSOs) in the Rho Ophiuchi F core and partially in the E core of the L1688 dark cloud. We run a proper motion program at the ESO New Technology Telescope (NTT) with the Son of ISAAC (SOFI) instrument over nine years in the near-infrared. We complemented these observations with various public image databases to enlarge the time base of observations and the field of investigation to 0.5 deg X 0.5 deg. We derived positions and proper motions for 2213 objects. From these, 607 proper motions were derived from SOFI observations with a ~1.8 mas/yr accuracy while the remaining objects were measured only from auxiliary data with a mean precision of about ~3 mas/yr. We performed a kinematic analysis of the most accurate proper motions derived in this work, which allowed us to separate cluster members from field stars and to derive the mean properties of the cluster. From the kinematic analysis we derived a list of 68 members and 14 candidate members, comprising 26 new objects with a high membership probability. These new members are generally fainter than the known ones. We measured a mean proper motion of (mu_RA*, mu_DEC)=(-8.2, -24.3)+/-0.8 mas/yr for the L1688 dark cloud. A supervised classification was applied to photometric data of members to allocate a spectral energy distribution (SED) classification to the unclassified members.} We kinematically confirmed that the 56 members that were known from previous studies of the Rho Ophiuchi F cluster and that were also part of our survey are members of the cluster, and we added 26 new members. We defined the evolutionary status of the unclassified members of the cluster. We showed that a large part (23) of these new members are probably brown dwarfs, which multiplies the number of known substellar objects in the cluster by a factor of 3.3.

We developed a hierarchical Bayesian model (HBM) to investigate how the presence of Seyfert activity relates to their environment, herein represented by the galaxy cluster mass, $M_{200}$, and the normalized cluster-centric distance, $r/r_{200}$. We achieved this by constructing an unbiased sample of galaxies from the Sloan Digital Sky Survey, with morphological classifications provided by the Galaxy Zoo Project. A propensity score matching approach is introduced to control for the effects of confounding variables: stellar mass, galaxy colour, and star formation rate. The connection between Seyfert-activity and environmental properties in the de-biased sample is modelled within an HBM framework using the so-called logistic regression technique, suitable for the analysis of binary data (e.g., whether or not a galaxy hosts an AGN). Unlike standard ordinary least square fitting methods, our methodology naturally allows modelling the probability of Seyfert-AGN activity in galaxies on their natural scale, i.e. as a binary variable. Furthermore, we demonstrate how an HBM can incorporate information of each particular galaxy morphological type in a unified framework. In elliptical galaxies, our analysis indicates a strong correlation of Seyfert-AGN activity with $r/r_{200}$, and a weaker correlation with the mass of the host. In spiral galaxies these trends do not appear, suggesting that the link between Seyfert activity and the properties of spiral galaxies are independent of the environment.

Machine learning techniques offer a plethora of opportunities in tackling big data within the astronomical community. We present the set of Generalized Linear Models as a fast alternative for determining photometric redshifts of galaxies, a set of tools not commonly applied within astronomy, despite being widely used in other professions. With this technique, we achieve catastrophic outlier rates of the order of ~1%, that can be achieved in a matter of seconds on large datasets of size ~1,000,000. To make these techniques easily accessible to the astronomical community, we developed a set of libraries and tools that are publicly available.

(Abridged) Gaia aims to make a 3-dimensional map of 1,000 million stars in our Milky Way to unravel its kinematical, dynamical, and chemical structure and evolution. Gaia's on-board detection software discriminates stars from spurious objects like cosmic rays and Solar protons. For this, parametrised point-spread-function-shape criteria are used. This study aims to provide an optimum set of parameters for these filters. We developed an emulation of the on-board detection software, which has 20 free, so-called rejection parameters which govern the boundaries between stars on the one hand and sharp or extended events on the other hand. We evaluate the detection and rejection performance of the algorithm using catalogues of simulated single stars, double stars, cosmic rays, Solar protons, unresolved galaxies, and asteroids. We optimised the rejection parameters, improving - with respect to the functional baseline - the detection performance of single and double stars, while, at the same time, improving the rejection performance of cosmic rays and of Solar protons. We find that the minimum separation to resolve a close, equal-brightness double star is 0.23 arcsec in the along-scan and 0.70 arcsec in the across-scan direction, independent of the brightness of the primary. We find that, whereas the optimised rejection parameters have no significant impact on the detectability of de Vaucouleurs profiles, they do significantly improve the detection of exponential-disk profiles. We also find that the optimised rejection parameters provide detection gains for asteroids fainter than 20 mag and for fast-moving near-Earth objects fainter than 18 mag, albeit this gain comes at the expense of a modest detection-probability loss for bright, fast-moving near-Earth objects. The major side effect of the optimised parameters is that spurious ghosts in the wings of bright stars essentially pass unfiltered.

We report the first analytical expression purely constructed by a machine to determine photometric redshifts ($z_{\rm phot}$) of galaxies. A simple and reliable functional form is derived using $41,214$ galaxies from the Sloan Digital Sky Survey Data Release 10 (SDSS-DR10) spectroscopic sample. The method automatically dropped the $u$ and $z$ bands, relying only on $g$, $r$ and $i$ for the final solution. Applying this expression to other $1,417,181$ SDSS-DR10 galaxies, with measured spectroscopic redshifts ($z_{\rm spec}$), we achieved a mean $\langle (z_{\rm phot} - z_{\rm spec})/(1+z_{\rm spec})\rangle\lesssim 0.0086$ and a scatter $\sigma_{(z_{\rm phot} - z_{\rm spec})/(1+z_{\rm spec})}\lesssim 0.045$ when averaged up to $z \lesssim 1.0$. The method was also applied to the PHAT0 dataset, confirming the competitiveness of our results when faced with other methods from the literature. This is the first use of symbolic regression in cosmology, representing a leap forward in astronomy-data-mining connection.

Besides its major objective tuned to the detection of the stellar galactic population the Gaia mission experiment will also observe a large number of galaxies. In this work we intend to evaluate the number and the characteristics of the galaxies that will effectively pass the onboard selection algorithm of Gaia. The detection of objects in Gaia will be performed in a section of the focal plane known as the Sky Mapper. Taking into account the Video Processing Algorithm criterion of detection and considering the known light profiles of discs and bulges galaxies we assess the number and the type of extra-galactic objects that will be observed by Gaia. We show that the stellar disk population of galaxies will be very difficult to be observed. On the contrary the spheroidal component of elliptical galaxies and bulges having higher central surface brightness and steeper brightness profile will be more easy to be detected. We estimate that most of the 20 000 elliptical population of nearby galaxies inside the local region up to 170 Mpc are in condition to be observed by Gaia. A similar number of bulges could also be observed although the low luminosity bulges should escape detection. About two thirds of the more distant objects up to 600 Mpc could also be detected increasing the total sample to half a million objects including ellipticals and bulges. The angular size of the detected objects will never exceed 4.72 arcsec which is the size of the largest transmitted windows. An heterogeneous population of elliptical galaxies and bulges will be observable by Gaia. This nearby Universe sample of galaxies should constitute a very rich and interesting sample to study their structural properties and their distribution.

The nearby TW Hydrae Association (TWA) is currently a benchmark for the study of formation and evolution of young low-mass stars, circumstellar disks and the imaging detection of planetary companions. For such studies, it is crucial to evaluate the distance to group members in order to access their physical properties. Membership of several stars is strongly debated and age estimates vary from one author to another with doubts about coevality. We revisit the kinematic properties of the TWA in light of new trigonometric parallaxes and proper motions to derive the dynamical age of the association and physical parameters of kinematic members. Using observations performed with the NTT/ESO telescope we measured trigonometric parallaxes and proper motions for 13 stars in TWA. With the convergent point method we identify a co-moving group with 31 TWA stars. We deduce kinematic distances for 7 members of the moving group that lack trigonometric parallaxes. A traceback strategy is applied to the stellar space motions of a selection of 16 of the co-moving objects with accurate and reliable data yielding a dynamical age for the association of t = 7.5 +/- 0.7 Myr. Using our new parallaxes and photometry available in the literature we derive stellar ages and masses from theoretical evolutionary models. With new parallax and proper motion measurements from this work and current astrometric catalogs we provide an improved and accurate database for TWA stars to be used in kinematical analysis. We conclude that the dynamical age obtained via traceback strategy is consistent with previous age estimates for the TWA, and is also compatible with the average ages derived in the present paper from evolutionary models for pre-main sequence stars.

We develop a method for membership assignment in stellar clusters using only photometry and positions. The method, UPMASK, is aimed to be unsupervised, data driven, model free, and to rely on as few assumptions as possible. It is based on an iterative process, principal component analysis, clustering algorithm, and kernel density estimations. Moreover, it is able to take into account arbitrary error models. An implementation in R was tested on simulated clusters that covered a broad range of ages, masses, distances, reddenings, and also on real data of cluster fields. Running UPMASK on simulations showed that it effectively separates cluster and field populations. The overall spatial structure and distribution of cluster member stars in the colour-magnitude diagram were recovered under a broad variety of conditions. For a set of 360 simulations, the resulting true positive rates (a measurement of purity) and member recovery rates (a measurement of completeness) at the 90% membership probability level reached high values for a range of open cluster ages ($10^{7.1}-10^{9.5}$ yr), initial masses ($0.5-10\times10^3$M$_{\sun}$) and heliocentric distances ($0.5-4.0$ kpc). UPMASK was also tested on real data from the fields of the open cluster Haffner~16 and of the closely projected clusters Haffner~10 and Czernik~29. These tests showed that even for moderate variable extinction and cluster superposition, the method yielded useful cluster membership probabilities and provided some insight into their stellar contents. The UPMASK implementation will be available at the CRAN archive.

The ESA Gaia mission, to be launched during 2013, will observe billions of objects, among which many galaxies, during its scanning of the sky. This will provide a large space-based dataset with unprecedented spatial resolution. Because of its natural Galactic and Astrometric priority, Gaia's observational strategy was optimized for point sources. Nonetheless, it is expected that 10^6 sources will be extragalactic, and a large portion of them will be angularly small galaxies. Although the mission was designed for point sources, a dedicated analysis of the raw data will allow the recovery of morphology of those objects at a 0.2" level. This may constitute a unique all-sky survey of such galaxies. We describe the conceptual design of the method we created for performing the morphological analysis of these objects as well as first results obtained from data simulations of low-resolution, highly binned, satellite data. Based on the obtained results we conclude that it is possible to push the limits of the Gaia space mission by analyzing galaxy morphology. (Abridged)

We construct a theoretical model to predict the number of orphan afterglows (OA) from gamma-ray bursts (GRBs) triggered by primordial metal-free (Pop III) stars expected to be observed by the Gaia mission. In particular, we consider primordial metal-free stars that were affected by radiation from other stars (Pop III.2) as a possible target. We use a semi-analytical approach that includes all relevant feedback effects to construct cosmic star formation history and its connection with the cumulative number of GRBs. The OA events are generated using the Monte Carlo method, and realistic simulations of Gaia's scanning law are performed to derive the observation probability expectation. We show that Gaia can observe up to 2.28 $\pm$ 0.88 off-axis afterglows and 2.78 $\pm$ 1.41 on-axis during the five-year nominal mission. This implies that a nonnegligible percentage of afterglows that may be observed by Gaia ($\sim 10%$) could have Pop III stars as progenitors.

Context. This study has been developed in the framework of the computational simulations executed for the preparation of the ESA Gaia astrometric mission. Aims. We focus on describing the objects and characteristics that Gaia will potentially observe without taking into consideration instrumental effects (detection efficiency, observing errors). Methods. The theoretical Universe Model prepared for the Gaia simulation has been statistically analyzed at a given time. Ingredients of the model are described, giving most attention to the stellar content, the double and multiple stars, and variability. Results. In this simulation the errors have not been included yet. Hence we estimate the number of objects and their theoretical photometric, astrometric and spectroscopic characteristics in the case that they are perfectly detected. We show that Gaia will be able to potentially observe 1.1 billion of stars (single or part of multiple star systems) of which about 2% are variable stars, 3% have one or two exoplanets. At the extragalactic level, observations will be potentially composed by several millions of galaxies, half million to 1 million of quasars and about 50,000 supernovas that will occur during the 5 years of mission. The simulated catalogue will be made publicly available by the DPAC on the Gaia portal of the ESA web site http://www.rssd.esa.int/gaia/.

(abridged) A complete census of planetary systems around a volume-limited sample of solar-type stars (FGK dwarfs) in the Solar neighborhood with uniform sensitivity down to Earth-mass planets within their Habitable Zones out to several AUs would be a major milestone in extrasolar planets astrophysics. This fundamental goal can be achieved with a mission concept such as NEAT - the Nearby Earth Astrometric Telescope. NEAT is designed to carry out space-borne extremely-high-precision astrometric measurements sufficient to detect dynamical effects due to orbiting planets of mass even lower than Earth's around the nearest stars. Such a survey mission would provide the actual planetary masses and the full orbital geometry for all the components of the detected planetary systems down to the Earth-mass limit. The NEAT performance limits can be achieved by carrying out differential astrometry between the targets and a set of suitable reference stars in the field. The NEAT instrument design consists of an off-axis parabola single-mirror telescope, a detector with a large field of view made of small movable CCDs located around a fixed central CCD, and an interferometric calibration system originating from metrology fibers located at the primary mirror. The proposed mission architecture relies on the use of two satellites operating at L2 for 5 years, flying in formation and offering a capability of more than 20,000 reconfigurations (alternative option uses deployable boom). The NEAT primary science program will encompass an astrometric survey of our 200 closest F-, G- and K-type stellar neighbors, with an average of 50 visits. The remaining time might be allocated to improve the characterization of the architecture of selected planetary systems around nearby targets of specific interest (low-mass stars, young stars, etc.) discovered by Gaia, ground-based high-precision radial-velocity surveys.

We derive lists of proper-motions and kinematic membership probabilities for 49 open clusters and possible open clusters in the zone of the Bordeaux PM2000 proper motion catalogue ($+11^{\circ}\le\delta\le+18^{\circ}$). We test different parametrisations of the proper motion and position distribution functions and select the most successful one. In the light of those results, we analyse some objects individually. The segregation between cluster and field member stars, and the assignment of membership probabilities, is accomplished by applying a new and fully automated method based on both parametrisations of the proper motion and position distribution functions, and genetic algorithm optimization heuristics associated with a derivative-based hill climbing algorithm for the likelihood optimization. We present a catalogue comprising kinematic parameters and associated membership probability lists for 49 open clusters and possible open clusters in the Bordeaux PM2000 catalogue region. We note that this is the first determination of proper motions for five open clusters. We confirm the non-existence of two kinematic populations in the region of 15 previously suspected non-existent objects.

TWA22 was initially regarded as a member of the TW Hydrae association (TWA). In addition to being one of the youngest (~8Myr) and nearest (~20pc) stars to Earth, TWA22 has proven to be very interesting after being resolved as a tight, very low-mass binary. This binary can serve as a very useful dynamical calibrator for pre-main sequence evolutionary models. However, its membership in the TWA has been recently questioned despite due to the lack of accurate kinematic measurements. Based on proper motion, radial velocity, and trigonometric parallax measurements, we aim here to re-analyze the membership of TWA22 to young, nearby associations. Using the ESO NTT/SUSI2 telescope, we observed TWA22AB during 5 different observing runs over 1.2 years to measure its trigonometric parallax and proper motion. This is a part of a larger project measuring trigonometric parallaxes and proper motions of most known TWA members at a sub-milliarcsec level. HARPS at the ESO 3.6m telescope was also used to measure the system's radial velocity over 2 years. We report an absolute trigonometric parallax of TWA22AB, 57.0mas, corresponding to a distance 17.5pc from Earth. Measured proper motions of TWA22AB are -175.8mas/yr in right ascension and -21.3mas/yr in declination. Finally, from HARPS measurements, we obtain a radial velocity 14.8km/s.

In 2005, Scholz and collaborators (Scholz et al. 2005) discovered, in a proper motion survey, a young brown dwarf SSSPMJ1102-3431(SSSPMJ1102) of spectral type M8.5, probable member of the TW Hydrae Association (TWA) and possible companion of the T Tauri star TW Hya. The physical characterization of SSSPMJ1102 was based on the hypothesis that it forms a binary system with TW Hya. The recent discovery of a probable giant planet inside the TW Hya protoplanetary disk with a very short-period (Setiawan et al. 2008) and a disk around SSSPMJ1102 (Riaz and Gizis 2008) make it especially interesting and important to measure well the physical parameters of SSSPMJ1102. Trigonometric parallax and proper motion measurements of SSSPMJ1102 are necessary to test for TWA membership and, thus, to determine the mass and age of this young brown dwarf and the possibility that it forms a wide binary system with TW Hya. Two years of regular observations at the ESO NTT/SUSI2 telescope, have enabled us to determine the trigonometric parallax and proper motion of SSSPMJ1102. Our parallax and proper motion determination allow us to precisely describe the physical properties of this low mass object and to confirm its TWA membership. Our results are not incompatible with the hypothesis that SSSPMJ1102 is a binary companion of the star TW Hya.

The central region of our Galaxy is a very peculiar environment, containing magnetic fields in excess of 100 mG and gas densities reaching ~ 10^4cm^-3. This region was observed as a strong source of GeV and TeVs gammas, what suggests that a mechanism of proton-neutron conversion could be taking place therein. We propose that the Galactic Center must also be a source of EeV neutrons due to the conversion of ultra high energy cosmic ray protons into neutrons via p-p interactions inside this region. This scenario should be falsifiable by the Pierre Auger Observatory after a few years of full exposure.