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.

We aim to investigate the connections existing between the density profiles of the stellar populations used to define a gravitationally bound stellar system and their star formation history: we do this by developing a general framework accounting for both classical stellar population theory and classical stellar dynamics. We extend the work of Pasetto et al. (2012) on a single composite-stellar population (CSP) to multiple CSPs, including also a phase-space description of the CSP concept. In this framework, we use the concept of distribution function to define the CSP in terms of mass, metallicity, and phase-space in a suitable space of existence $\mathbb{E}$ of the CSP. We introduce the concept of foliation of $\mathbb{E}$ to describe formally any CSP as sum of disjointed Simple Stellar Populations (SSP), with the aim to offer a more general formal setting to cast the equations of stellar populations theory and stellar dynamics theory. In doing so, we allow the CSP to be object of dissipation processes thus developing its dynamics in a general non-Hamiltonian framework. Furthermore, we investigate the necessary and sufficient condition to realize a multiple CSP consistent with its mass-metallicity and phase-space distribution function over its temporal evolution, for a collisionless CSP. Finally, analytical and numerical examples show the potential of the result obtained.

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

Aims. We describe the photometric content of the second data release of the Gaia project (Gaia DR2) and its validation along with the quality of the data. Methods. The validation was mainly carried out using an internal analysis of the photometry. External comparisons were also made, but were limited by the precision and systematics that may be present in the external catalogues used. Results. In addition to the photometric quality assessment, we present the best estimates of the three photometric passbands. Various colour-colour transformations are also derived to enable the users to convert between the Gaia and commonly used passbands. Conclusions. The internal analysis of the data shows that the photometric calibrations can reach a precision as low as 2 mmag on individual CCD measurements. Other tests show that systematic effects are present in the data at the 10 mmag level.

We present a new Galaxy population synthesis Model (GalMod). GalMod is a star-count model featuring an asymmetric bar/bulge as well as spiral arms and related extinction. The model, initially introduced in Pasetto et al. (2016b), has been here completed with a central bar, a new bulge description, new disk vertical profiles and several new bolometric corrections. The model can generate synthetic mock catalogs of visible portions of the Milky Way (MW), external galaxies like M31, or N-body simulation initial conditions. At any given time, e.g., a chosen age of the Galaxy, the model contains a sum of discrete stellar populations, namely bulge/bar, disk, halo. These populations are in turn the sum of different components: the disk is the sum of spiral arms, thin disks, a thick disk, and various gas components, while the halo is the sum of a stellar component, a hot coronal gas, and a dark matter component. The Galactic potential is computed from these population density profiles and used to generate detailed kinematics by considering up to the first four moments of the collisionless Boltzmann equation. The same density profiles are then used to define the observed color-magnitude diagrams in a user-defined field of view from an arbitrary solar location. Several photometric systems have been included and made available on-line and no limits on the size of the field of view are imposed thus allowing full sky simulations, too. Finally, we model the extinction adopting a dust model with advanced ray-tracing solutions. The model's web page (and tutorial) can be accessed at www.GalMod.org.

The second Gaia data release is based on 22 months of mission data with an average of 0.9 billion individual CCD observations per day. A data volume of this size and granularity requires a robust and reliable but still flexible system to achieve the demanding accuracy and precision constraints that Gaia is capable of delivering. The internal Gaia photometric system was initialised using an iterative process that is solely based on Gaia data. A set of calibrations was derived for the entire Gaia DR2 baseline and then used to produce the final mean source photometry. The photometric catalogue contains 2.5 billion sources comprised of three different grades depending on the availability of colour information and the procedure used to calibrate them: 1.5 billion gold, 144 million silver, and 0.9 billion bronze. These figures reflect the results of the photometric processing; the content of the data release will be different due to the validation and data quality filters applied during the catalogue preparation. The photometric processing pipeline, PhotPipe, implements all the processing and calibration workflows in terms of Map/Reduce jobs based on the Hadoop platform. This is the first example of a processing system for a large astrophysical survey project to make use of these technologies. The improvements in the generation of the integrated G-band fluxes, in the attitude modelling, in the cross-matching, and and in the identification of spurious detections led to a much cleaner input stream for the photometric processing. This, combined with the improvements in the definition of the internal photometric system and calibration flow, produced high-quality photometry. Hadoop proved to be an excellent platform choice for the implementation of PhotPipe in terms of overall performance, scalability, downtime, and manpower required for operations and maintenance.

The European Space Agency Gaia satellite was launched into orbit around L2 in December 2013. This ambitious mission has strict requirements on residual systematic errors resulting from instrumental corrections in order to meet a design goal of sub-10 microarcsecond astrometry. During the design and build phase of the science instruments, various critical calibrations were studied in detail to ensure that this goal could be met in orbit. In particular, it was determined that the video-chain offsets on the analogue side of the analogue-to-digital conversion electronics exhibited instabilities that could not be mitigated fully by modifications to the flight hardware. We provide a detailed description of the behaviour of the electronic offset levels on microsecond timescales, identifying various systematic effects that are known collectively as offset non-uniformities. The effects manifest themselves as transient perturbations on the gross zero-point electronic offset level that is routinely monitored as part of the overall calibration process. Using in-orbit special calibration sequences along with simple parametric models, we show how the effects can be calibrated, and how these calibrations are applied to the science data. While the calibration part of the process is relatively straightforward, the application of the calibrations during science data processing requires a detailed on-ground reconstruction of the readout timing of each charge-coupled device (CCD) sample on each device in order to predict correctly the highly time-dependent nature of the corrections. We demonstrate the effectiveness of our offset non-uniformity models in mitigating the effects in Gaia data. We demonstrate for all CCDs and operating instrument and modes on board Gaia that the video-chain noise-limited performance is recovered in the vast majority of science samples.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

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

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

Aims. The photometric validation of the Gaia DR1 release of the ESA Gaia mission is described and the quality of the data shown. Methods. This is carried out via an internal analysis of the photometry using the most constant sources. Comparisons with external photometric catalogues are also made, but are limited by the accuracies and systematics present in these catalogues. An analysis of the quoted errors is also described. Investigations of the calibration coefficients reveal some of the systematic effects that affect the fluxes. Results. The analysis of the constant sources shows that the early-stage photometric calibrations can reach an accuracy as low as 3 mmag.

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.

Context. Gaia is an ESA cornerstone mission launched on 19 December 2013 aiming to obtain the most complete and precise 3D map of our Galaxy by observing more than one billion sources. This paper is part of a series of documents explaining the data processing and its results for Gaia Data Release 1, focussing on the G band photometry. Aims. This paper describes the calibration model of the Gaia photometric passband for Gaia Data Release 1. Methods. The overall principle of splitting the process into internal and external calibrations is outlined. In the internal calibration, a self-consistent photometric system is generated. Then, the external calibration provides the link to the absolute photometric flux scales. Results. The Gaia photometric calibration pipeline explained here was applied to the first data release with good results. Details are given of the various calibration elements including the mathematical formulation of the models used and of the extraction and preparation of the required input parameters (e.g. colour terms). The external calibration in this first release provides the absolute zero point and photometric transformations from the Gaia G passband to other common photometric systems. Conclusions. This paper describes the photometric calibration implemented for the first Gaia data release and the instrumental effects taken into account. For this first release no aperture losses, radiation damage, and other second-order effects have not yet been implemented in the calibration.

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

We report the discovery and characterisation of a deeply eclipsing AM CVn-system, Gaia14aae (= ASSASN-14cn). Gaia14aae was identified independently by the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014) and by the Gaia Science Alerts project, during two separate outbursts. A third outburst is seen in archival Pan-STARRS-1 (PS1; Schlafly et al. 2012; Tonry et al. 2012; Magnier et al. 2013) and ASAS-SN data. Spectroscopy reveals a hot, hydrogen-deficient spectrum with clear double-peaked emission lines, consistent with an accreting double degenerate classification. We use follow-up photometry to constrain the orbital parameters of the system. We find an orbital period of 49.71 min, which places Gaia14aae at the long period extremum of the outbursting AM CVn period distribution. Gaia14aae is dominated by the light from its accreting white dwarf. Assuming an orbital inclination of 90 degrees for the binary system, the contact phases of the white dwarf lead to lower limits of 0.78 M solar and 0.015 M solar on the masses of the accretor and donor respectively and a lower limit on the mass ratio of 0.019. Gaia14aae is only the third eclipsing AM CVn star known, and the first in which the WD is totally eclipsed. Using a helium WD model, we estimate the accretor's effective temperature to be 12900+-200 K. The three out-burst events occurred within 4 months of each other, while no other outburst activity is seen in the previous 8 years of Catalina Real-time Transient Survey (CRTS; Drake et al. 2009), Pan-STARRS-1 and ASAS-SN data. This suggests that these events might be rebrightenings of the first outburst rather than individual events.

I present multiband optical and UV Hubble Space Telescope photometry of the two Galactic globular clusters NGC6388 and NGC6441, in order to investigate the nature of the physical mechanism(s) responsible for the existence of an extended blue tail and of a slope in the horizontal branch. Further evidence that the horizontal branch tilt cannot be interpreted as an effect of differential reddening is provided, while I show that a possible solution of the puzzle is to assume that a small fraction of the stellar population in the two clusters is strongly helium enriched.

The UV excess shown by elliptical galaxies in their spectra is believed to be caused by evolved low-mass stars, in particular sdB stars. The stellar system most similar to the ellipticals for age and metallicity, in which it is possible to resolve these stars, is the bulge of our Galaxy. sdB star candidates were observed in the color magnitude diagram of a bulge region by Zoccali et al. (2003). The follow-up spectroscopic analysis of these stars confirmed that most of these stars are bulge sdBs, while some candidates turned out to be disk sdBs or cool stars. Both spectroscopic and photometric data and a spectral library are used to construct the integrated spectrum of the observed bulge region from the UV to the optical: the stars in the color magnitude diagram are associated to the library spectra, on the basis of their evolutionary status and temperature. The total integrated spectrum is obtained as the sum of the spectra associated to the color magnitude diagram. The comparison of the obtained integrated spectrum with old single stellar population synthetic spectra calculated by Bruzual & Charlot (2003) agrees with age and metallicity of the bulge found by previous work. The bulge integrated spectrum shows only a very weak UV excess, but a too strict selection of the sample of the sdB star candidates in the color magnitude diagram and the exclusion of post-Asymptotic Giant Branch stars could have influenced the result.

In this paper we present multiband optical and UV Hubble Space Telescope photometry of the two Galactic globular clusters NGC 6388 and NGC 6441. We have obtained the first UV color-magnitude diagrams for NGC 6388 and NGC 6441. These diagrams confirm previous results, obtained in optical bands, about the presence of a sizeable stellar population of extremely hot horizontal branch stars. At least in NGC 6388, we find a clear indication that at the hot end of the horizontal branch the distribution of stars forms a hook-like feature, closely resembling those observed in NGC 2808 and Omega Centauri. We briefly review the theoretical scenarios that have been suggested for interpreting this observational feature. We also investigate the tilted horizontal branch morphology and provide further evidence that supports early suggestions that this feature cannot be interpreted as an effect of differential reddening. We show that a possible solution of the puzzle is to assume that a small fraction - ranging between 10-20% - of the stellar population in the two clusters is strongly helium enriched (Y=0.40 in NGC 6388 and Y=0.35 in NGC 6441). The occurrence of a spread in the He abundance between the canonical value (Y=0.26) and the quoted upper limits can significantly help in explaining the whole morphology of the horizontal branch and the pulsational properties of the variable stars in the target clusters.

Recent observations and theories suggest that extreme horizontal branch (EHB) stars and their progeny should be the cause of the UV excess seen in the spectra of many elliptical galaxies. Since the Galactic Bulge is the closest representation of an old, metal-rich spheroid in which we are able to study the EHB scenario in detail, we obtained spectra of bulge EHB star candidates and we confirm their status as hot evolved stars. It is the first time that such stars are unambiguously identified in the Galactic Bulge.

We present accurate relative ages for a sample of 55 Galactic globular clusters. The ages have been obtained by measuring the difference between the horizontal branch and the turnoff in two, internally photometrically homogeneous databases. The mutual consistency of the two data sets has been assessed by comparing the ages of 16 globular clusters in common between the two databases. We have also investigated the consistency of our relative age determination within the recent stellar model framework. All clusters with [Fe/H]<-1.7 are found to be old, and coeval, with the possible exception of two objects, which are marginally younger. The age dispersion for the metal poor clusters is 0.6 Gyr (rms), consistent with a null age dispersion. Intermediate metallicity clusters (-1.7<[Fe/H]<-0.8) are on average 1.5 Gyr younger than the metal poor ones, with an age dispersion of 1.0 Gyr (rms), and a total age range of ~3 Gyr. About 15% of the intermediate metallicity clusters are coeval with the oldest clusters. All the clusters with [Fe/H]>-0.8 are ~1 Gyr younger than the most metal poor ones, with a relatively small age dispersion, though the metal rich sample is still too small to allow firmer conclusions. There is no correlation of the cluster age with the Galactocentric distance. We briefly discuss the implication of these observational results for the formation history of the Galaxy.

We present color-magnitude diagrams (CMD) of the globular cluster NGC 6388 based on HST multiband photometry (F255W, F336W, F439W, F555W). In this paper we focus our attention on the peculiar horizontal branch of this cluster. After a careful reddening correction, we fitted the observed CMD with theoretical models. For the first time we demonstrated that the HB of this very metal rich globular cluster extends beyond Te > 30.000 K, showing clear evidence of a population of blue hook (D'Cruz et al.2000, Brown et al. 2001) stars. Moreover, we could demonstrate that the HB tilt (slope) is not a reddening effect and is present in all the photometric bands.