Gaia is a cornerstone mission in the science programme of the European Space Agency (ESA). The spacecraft construction was approved in 2006, following a study in which the original interferometric concept was changed to a direct-imaging approach. Both the spacecraft and the payload were built by European industry. The involvement of the scientific community focusses on data processing for which the international Gaia Data Processing and Analysis Consortium (DPAC) was selected in 2007. Gaia was launched on 19 December 2013 and arrived at its operating point, the second Lagrange point of the Sun-Earth-Moon system, a few weeks later. The commissioning of the spacecraft and payload was completed on 19 July 2014. The nominal five-year mission started with four weeks of special, ecliptic-pole scanning and subsequently transferred into full-sky scanning mode. We recall the scientific goals of Gaia and give a description of the as-built spacecraft that is currently (mid-2016) being operated to achieve these goals. We pay special attention to the payload module, the performance of which is closely related to the scientific performance of the mission. We provide a summary of the commissioning activities and findings, followed by a description of the routine operational mode. We summarise scientific performance estimates on the basis of in-orbit operations. Several intermediate Gaia data releases are planned and the data can be retrieved from the Gaia Archive, which is available through the Gaia home page.

Context. At about 1000 days after the launch of Gaia we present the first Gaia data release, Gaia DR1, consisting of astrometry and photometry for over 1 billion sources brighter than magnitude 20.7. Aims. A summary of Gaia DR1 is presented along with illustrations of the scientific quality of the data, followed by a discussion of the limitations due to the preliminary nature of this release. Methods. The raw data collected by Gaia during the first 14 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium (DPAC) and turned into an astrometric and photometric catalogue. Results. Gaia DR1 consists of three components: a primary astrometric data set which contains the positions, parallaxes, and mean proper motions for about 2 million of the brightest stars in common with the Hipparcos and Tycho-2 catalogues -a realisation of the Tycho-Gaia Astrometric Solution (TGAS) -and a secondary astrometric data set containing the positions for an additional 1.1 billion sources. The second component is the photometric data set, consisting of mean G-band magnitudes for all sources. The G-band light curves and the characteristics of ∼3000 Cepheid and RR Lyrae stars, observed at high cadence around the south ecliptic pole, form the third component. For the primary astrometric data set the typical uncertainty is about 0.3 mas for the positions and parallaxes, and about 1 mas yr −1 for the proper motions. A systematic component of ∼0.3 mas should be added to the parallax uncertainties. For the subset of ∼94 000 Hipparcos stars in the primary data set, the proper motions are much more precise at about 0.06 mas yr −1 . For the secondary astrometric data set, the typical uncertainty of the positions is ∼10 mas. The median uncertainties on the mean G-band magnitudes range from the mmag level to ∼0.03 mag over the magnitude range 5 to 20.7. Conclusions. Gaia DR1 is an important milestone ahead of the next Gaia data release, which will feature five-parameter astrometry for all sources. Extensive validation shows that Gaia DR1 represents a major advance in the mapping of the heavens and the availability of basic stellar data that underpin observational astrophysics. Nevertheless, the very preliminary nature of this first Gaia data release does lead to a number of important limitations to the data quality which should be carefully considered before drawing conclusions from the data.

We present the results of a systematic search for massive black hole binaries in the Sloan Digital Sky Survey spectroscopic database. We focus on bound binaries, under the assumption that one of the black holes is active. In this framework, the broad lines associated to the accreting black hole are expected to show systematic velocity shifts with respect to the narrow lines, which trace the rest-frame of the galaxy. For a sample of 54 586 quasars and 3 929 galaxies at redshifts 0.1 < z < 1.5 we bruteforce model each spectrum as a mixture of two quasars at two different redshifts. The spectral model is a data-driven dimensionality reduction of the SDSS quasar spectra based on a matrix factorization. We identified 32 objects with peculiar spectra. Nine of them can be interpreted as black hole binaries. This doubles the number of known black hole binary candidates. We also report on the discovery of a new class of extreme double-peaked emitters with exceptionally broad and faint Balmer lines. For all the interesting sources, we present detailed analysis of the spectra, and discuss possible interpretations.

Theoretically, bound binaries of massive black holes are expected as the natural outcome of mergers of massive galaxies. From the observational side, however, massive black hole binaries remain elusive. Velocity shifts between narrow and broad emission lines in quasar spectra are considered a promising observational tool to search for spatially unresolved, dynamically bound binaries. In this series of papers we investigate the nature of such candidates through analyses of their spectra, images and multi-wavelength spectral energy distributions. Here we investigate the properties of the optical spectra, including the evolution of the broad line profiles, of all the sources identified in our previous study. We find a diverse phenomenology of broad and narrow line luminosities, widths, shapes, ionization conditions and time variability, which we can broadly ascribe to 4 classes based on the shape of the broad line profiles: 1) Objects with bell-shaped broad lines with big velocity shifts (>1000 km s −1 ) compared to their narrow lines show a variety of broad line widths and luminosities, modest flux variations over a few years, and no significant change in the broad line peak wavelength. 2) Objects with double-peaked broad emission lines tend to show very luminous and broadened lines, and little time variability. 3) Objects with asymmetric broad emission lines show a broad range of broad line luminosities and significant variability of the line profiles. 4) The remaining sources tend to show moderate to low broad line luminosities, and can be ascribed to diverse phenomena. We discuss the implications of our findings in the context of massive black hole binary searches.

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 100 K in effective temperature (T eff ), 0.1 mag in extinction (A 0 ), 0.2 dex in metallicity ([Fe/H]), and 0.25 dex in surface gravity (log g). 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 December 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.