The recent discoveries of gravitational wave events and in one case also its electromagnetic (EM) counterpart allow us to study the Universe in a novel way. The increased sensitivity of the LIGO and Virgo detectors has opened the possibility for regular detections of EM transient events from mergers of stellar remnants. Gravitational wave sources are expected to have sky localisation up to a few hundred square degrees, thus Gaia as an all-sky multi-epoch photometric survey has the potential to be a good tool to search for the EM counterparts. In this paper we study the possibility of detecting EM counterparts to gravitational wave sources using the Gaia Science Alerts system. We develop an extension to current used algorithms to find transients and test its capabilities in discovering candidate transients on a sample of events from the observation periods O1 and O2 of LIGO and Virgo. For the gravitational wave events from the current run O3 we expect that about 16 (25) per cent should fall in sky regions observed by Gaia 7 (10) days after gravitational wave. The new algorithm will provide about 10 candidates per day from the whole sky.

The Gaia DR2 sample of short-timescale variable candidates results from the investigation of the first 22 months of Gaia photometry for a subsample of sources at the Gaia faint end. For this exercise, we limited ourselves to the case of suspected rapid periodic variability. Our study combines fast-variability detection through variogram analysis, high-frequency search by means of least-squares periodograms, and empirical selection based on the investigation of specific sources seen through the Gaia eyes (e.g. known variables or visually identified objects with peculiar features in their light curves). The progressive definition and validation of this selection criterion also benefited from supplementary ground-based photometric monitoring of a few preliminary candidates, performed at the Flemish Mercator telescope (Canary Islands, Spain) between August and November 2017. We publish a list of 3,018 short-timescale variable candidates, spread throughout the sky, with a false-positive rate up to 10-20% in the Magellanic Clouds, and a more significant but justifiable contamination from longer-period variables between 19% and 50%, depending on the area of the sky. Although its completeness is limited to about 0.05%, this first sample of Gaia short-timescale variables recovers some very interesting known short-period variables, such as post-common envelope binaries or cataclysmic variables, and brings to light some fascinating, newly discovered variable sources. In the perspective of future Gaia data releases, several improvements of the short-timescale variability processing are considered, by enhancing the existing variogram and period-search algorithms or by classifying the identified candidates. Nonetheless, the encouraging outcome of our Gaia DR2 analysis demonstrates the power of this mission for such fast-variability studies, and opens great perspectives for this domain of astrophysics.

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.

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.

We present the discovery of stellar tidal tails around the Large and the Small Magellanic Clouds in the Gaia DR1 data. In between the Clouds, their tidal arms are stretched towards each other to form an almost continuous stellar bridge. Our analysis relies on the exquisite quality of the Gaia's photometric catalogue to build detailed star-count maps of the Clouds. We demonstrate that the Gaia DR1 data can be used to detect variable stars across the whole sky, and in particular, RR Lyrae stars in and around the LMC and the SMC. Additionally, we use a combination of Gaia and Gale to follow the distribution of Young Main Sequence stars in the Magellanic System. Viewed by Gaia, the Clouds show unmistakable signs of interaction. Around the LMC, clumps of RR Lyrae are observable as far as ~20 degrees, in agreement with the most recent map of Mira-like stars reported in Deason et al (2016). The SMC's outer stellar density contours show a characteristic S-shape, symptomatic of the on-set of tidal stripping. Beyond several degrees from the center of the dwarf, the Gaia RR Lyrae stars trace the Cloud's trailing arm, extending towards the LMC. This stellar tidal tail mapped with RR Lyrae is not aligned with the gaseous Magellanic Bridge, and is shifted by some ~5 degrees from the Young Main Sequence bridge. We use the offset between the bridges to put constraints on the density of the hot gaseous corona of the Milky Way.

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.

Most stars are formed in a cluster or association, where the number density of stars can be high. This means that a large fraction of initially-single stars will undergo close encounters with other stars and/or exchange into binaries. We describe how such close encounters and exchange encounters can affect the properties of a planetary system around a single star. We define a singleton as a single star which has never suffered close encounters with other stars or spent time within a binary system. It may be that planetary systems similar to our own solar system can only survive around singletons. Close encounters or the presence of a stellar companion will perturb the planetary system, often leaving planets on tighter and more eccentric orbits. Thus planetary systems which initially resembled our own solar system may later more closely resemble some of the observed exoplanet systems.

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N-body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction are not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the solar system, the planets orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet-planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extra-solar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of solar-system-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the solar system around single stars, dynamical encounters in young stellar clusters may produce at least some of the extra-solar planetary systems observed in the solar neighbourhood.

We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)

Aims: This paper presents a new homogeneous catalogue of blue straggler stars (BSS) in Galactic open clusters. Methods: Photometric data for 216 clusters were collected from the literature and 2782 BSS candidates were extracted from 76 of them. Results: We found that the anticorrelation of BSS frequency vs. total magnitude identified in similar studies conducted on Galactic globular clusters extends to the open cluster regime: clusters with smaller total magnitude tend to have higher BSS frequencies. Moreover, a clear correlation between the BSS frequency (obtained normalising the total number of BSS either to the total cluster mass or, for the older clusters, to the total number of clump stars) and the age of the clusters was found. A simple model is developed here to try to explain this last and new result. The model allows us to ascertain the important effect played by mass loss in the evolution of open clusters.

We use two stellar populations in the globular cluster 47 Tucanae to trace its dynamical history: blue stragglers and low mass main sequence stars. We assumed that the blue stragglers were formed through stellar collisions in all regions of the cluster. We find that in the core of the cluster, models of collisional blue stragglers agree well with the observations as long as blue stragglers are still continuing to form and the mass function in the cluster is extremely biased towards massive stars (x=-8 where a Salpeter mass function has x=+1.35). We show that such an extreme mass function is supported by direct measurements of the luminosity function of main sequence stars in the centre of the cluster. In the middle region of our dataset (25'' to 130'' from the cluster centre), blue straggler formation seems to have stopped about half a Gyr ago. In the outskirts of the cluster, our models are least successful at reproducing the blue straggler data. Taken at face value, they indicate that blue straggler formation has been insignificant over the past billion years, and that a Salpeter mass function applies. However, it is more likely that the dominant formation mechanism in this part of the cluster is not the collisional one, and that our models are not appropriate for this region of the cluster. We conclude that blue stragglers can be used as tracers of dynamics in globular clusters, despite our incomplete understanding of how and where they were formed.

(Abridged) We derive the structure of the Galactic stellar Warp and Flare using 2MASS RC and RGB stars, selected at mean heliocentric distances of 3, 7 and 17 kpc. Our results are: (i) a clear stellar warp signature is derived for the 3 selected rings; (ii) the derived stellar warp is consistent (both in amplitude and phase-angle) with that for the Galactic interstellar dust and HI gas; (iii) the Sun seems not to fall on the line of nodes. The stellar warp phase-angle orientation (+15 degrees) is close to the orientation angle of the Galactic bar and this produces an asymmetric warp for the inner rings; (iv) a Northern/Southern warp symmetry is observed only for the ring at 17 kpc; (v) treating a mixture of thin and thick disk populations we trace the disk flaring and derive a constant scale-height (~0.65 kpc) within R(GC)~15 kpc. Further out, the disk flaring increase gradually reaching a mean scale-height of ~1.5 kpc at R(GC)~23 kpc; and (vi) these results provide further robust evidence that there is no disk radial truncation at R(GC)~14 kpc. In the particular case of the Canis Major over-density we confirm its coincidence with the Southern stellar maximum warp occurring near l=240. We present evidence to conclude that all observed parameters (e.g. number density, radial velocities, proper motion etc) of CMa are consistent with it being a normal Milky Way outer-disk population, thereby leaving no justification for a more complex interpretations of its origin. The present analysis does not provide a conclusive test of the structure or origin of the Monoceros Ring. Nevertheless, we show that a warped flared Milky Way contributes significantly at the locations of the Monoceros Ring.

The interpretation of globular cluster horizontal branch (HB) morphology is a classical problem that can significantly blur our understanding of stellar populations. In this paper, we present a new multivariate analysis connecting the effective temperature extent of the HB with other cluster parameters. The work is based on Hubble Space Telescope photometry of 54 Galactic globular clusters. The present study reveals an important role of the total mass of the globular cluster on its HB morphology. More massive clusters tend to have HBs more extended to higher temperatures. For a set of three input variables including the temperature extension of the HB, [Fe/H] and M_V, the first two eigenvectors account for the 90% of the total sample variance. Possible effects of cluster self-pollution on HB morphology, eventually stronger in more massive clusters, could explain the results here derived.

(abridged) Accurate reddenings for Globular Clusters could be obtained by comparing the colour-temperature obtained using temperatures from reddening-free indicator (Halpha), with that given by standard colour-temperature calibrations. From a single 1300 seconds exposure with FLAMES at VLT2 we obtained spectra centred on Halpha (R=6000, 5<S/N<50) for 120 stars near the turn-off of NGC6752. This set of spectra was used to derive effective temperatures from fittings of Halpha profiles with typical errors of about +/- 200 K and reddening estimates with individual errors of 0.05 mag. Averaging all individual reddenings, a high precision reddening estimate has been obtained for the cluster: E(B-V)=0.046 +/- 0.005. The same exposure provided UVES spectra of seven stars near the red giant branch bump (R=40000, 20<S/N<40) that provided Fe abundances with internal errors of 0.026 dex, and with average metallicity [Fe/H]= -1.48 +/- 0.01 +/- 0.06 dex (random + systematic). Abundances were obtained for several other elements, allowing e.g. an accurate estimate of the ratio between the alpha-elements and Fe ([alpha/Fe]=+0.27 +/- 0.01). The O-Na anticorrelation is evident from our UVES data, in agreement with past results. This analysis shows the power of FLAMES for analysis of globular clusters: the accurate reddenings and metal abundances obtained by a procedure like that described here, combined with distance determinations from cluster dynamics or main sequence fitting, and high quality colour-magnitude diagrams, could allow derivation of ages with errors below 1 Gyr for individual globular clusters.

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 distance modulus and reddening determinations for 72 Galactic globular clusters from the homogeneous photometric database of Piotto et al. (2002), calibrated to the HST flight F439W and F555W bands. The distances have been determined by comparison with theoretical absolute magnitudes of the ZAHB. For low and intermediate metallicity clusters, we have estimated the apparent Zero Age Horizontal Branch (ZAHB) magnitude from the RR Lyrae level. For metal rich clusters, the ZAHB magnitude was obtained from the fainter envelope of the red HB. Reddenings have been estimated by comparison of the HST colour-magnitude diagrams (CMD) with ground CMDs of low reddening template clusters. The homogeneity of both the photometric data and the adopted methodological approach allowed us to obtain highly accurate relative cluster distances and reddenings. Our results are also compared with recent compilations in the literature.

Proper-motion, star counts and photometric catalog simulations are used to explain the detected stellar over-density in the region of Canis Major (CMa), claimed to be the core of a disrupted dwarf galaxy (Martin et al. 2004, Bellazzini et al. 2003), as due to the Galactic warp and flare in the external disk. We compare the kinematics of CMa M-giant selected sample with surrounding Galactic disk stars in the UCAC2 catalog and find no peculiar proper motion signature: CMa stars mimic thick disk kinematics. Moreover, when taking into account the Galactic warp and flare of the disk, 2MASS star count profiles reproduce the CMa stellar over-density. This star count analysis is confirmed by direct comparison with synthetic color-magnitude diagrams simulated with the Besancon models (Robin et al. 2003) that include the warp and flare of the disk. The presented evidence casts doubt on the identification of the CMa over-density as the core of a disrupted Milky Way satellite. This however does not make clear the origin of over-densities responsible for the ring structure in the anticenter direction of the Galactic halo (Newberg et al. 2002; Yanny et al. 2003; Zaggia et al. 2004, in preparation).

In a previous paper we reported on a discontinuity in the extreme horizontal branch (EHB) of the Galactic globular cluster NGC6752, which we called the second U-jump. This feature was attributed to a combination of post zero-age horizontal branch evolution and diffusion effects. In this follow-up study we analyze other EHB clusters and show that the second U-jump is a common feature among EHB clusters reaching T_{eff}\ge 23,000K, and that its onset in different clusters converges around T_{eff}\sim 21,000\pm3,000K. We also present near-ultraviolet diagrams of \omega Cen and NGC2808, the only two objects with spectroscopically confirmed ``blue hook'' stars (T_{eff}\ge 35,000K). We confirm predictions of a photometric discontinuity separating late from early-helium flashers. Moreover, we present empirical evidence that the second U-jump population might be mainly composed by early-helium flashers. Lastly, we revisit the discussion on the ubiquitous nature of the gaps and jumps so far identified in the blue HB tails, suggesting a possible discrete nature of the distribution in temperature of the HB stars.

We discuss the main properties of the Galactic globular cluster (GC) blue straggler stars (BSS), as inferred from our new catalog containing nearly 3000 BSS. The catalog has been extracted from the photometrically homogeneous V vs. (B-V) color-magnitude diagrams (CMD) of 56 GCs, based on WFPC2 images of their central cores. In our analysis we used consistent relative distances based on the same photometry and calibration. The number of BSS has been normalized to obtain relative frequencies (F_{BSS}) and specific densities (N_S) using different stellar populations extracted from the CMD. The cluster F_{BSS} is significantly smaller than the relative frequency of field BSS. We find a significant anti-correlation between the BSS relative frequency in a cluster and its total absolute luminosity (mass). There is no statistically significant trend between the BSS frequency and the expected collision rate. F_{BSS} does not depend on other cluster parameters, apart from a mild dependence on the central density. PCC clusters act like normal clusters as far as the BSS frequency is concerned. We also show that the BSS luminosity function for the most luminous clusters is significantly different, with a brighter peak and extending to brighter luminosities than in the less luminous clusters. These results imply that the efficiency of BSS production mechanisms and their relative importance vary with the cluster mass.

Recent HST observations of a large sample of globular clusters reveal that every cluster contains between 40 and 400 blue stragglers. The population does not correlate with either stellar collision rate (as would be expected if all blue stragglers were formed via collisions) or total mass (as would be expected if all blue stragglers were formed via the unhindered evolution of a subset of the stellar population). In this paper, we support the idea that blue stragglers are made through both channels. The number produced via collisions tends to increase with cluster mass. In this paper we show how the current population produced from primordial binaries decreases with increasing cluster mass; exchange encounters with third, single, stars in the most massive clusters tend to reduce the fraction of binaries containing a primary close to the current turn-off mass. Rather their primaries tend to be somewhat more massive (~1-3 M_sun) and have evolved off the main sequence, filling their Roche lobes in the past, often converting their secondaries into blue stragglers (but more than 1 Gyr or so ago and thus they are no longer visible as blue stragglers). We show that this decline in the primordial blue straggler population is likely to be offset by the increase in the number of blue stragglers produced via collisions. The predicted total blue straggler population is therefore relatively independent of cluster mass, thus matching the observed population. This result does not depend on any particular assumed blue straggler lifetime.

We present observational estimates of the magnitude difference between the luminosity function RGB bump and the HB (Delta(bump-zahb)), and of star counts in the bump region (R_bump), for a sample of 54 Galactic globular clusters (GCs) observed by the HST. The large sample of stars resolved in each cluster, and the high photometric accuracy of the data allowed us to detect the bump also in a number of metal poor clusters. To reduce the photometric uncertainties, empirical values are compared with theoretical predictions obtained from a set of updated canonical stellar evolution models which have been transformed directly into the HST flight system. We found an overall qualitative agreement between theory and observations. Quantitative estimates of the confidence level are hampered by current uncertainties on the GCs metallicity scale, and by the strong dependence of Delta(bump-zahb) on the cluster metallicity. In case of the R_bump parameter, which is only weakly affected by the metallicity, we find a very good quantitative agreement between theoretical canonical models and observations. For our full cluster sample the average difference between predicted and observed R_bump values is practically negligible, and ranges from -0.002 to -0.028, depending on the employed metallicity scale. The observed dispersion around these values is entirely consistent with the observational errors on R_bump. As a comparison, the value of R_bump predicted by theory in case of spurious bump detections would be \~0.10 smaller than the observed ones. We have also tested the influence on the predicted Delta(bump-zahb) and R_bump values of an He-enriched component in the cluster stellar population. We find that the predicted Delta(bump-zahb) and R_bump values are only marginally affected.

In this paper we report on the most significant results from a statistical analysis of the main properties of globular cluster blue straggler stars (BSS) extracted from the HST snapshot database of photometrically homogeneous CMDs (Piotto et al. 2002). The BSS relative frequency presents a significant anticorrelation with the collisional rate and with the cluster total absolute luminosity.

We present the complete photometric database and the color-magnitude diagrams for 74 Galactic globular clusters observed with the HST/WFPC2 camera in the F439W and F555W bands. A detailed discussion of the various reduction steps is also presented, and of the procedures to transform instrumental magnitudes into both the HST F439W and F555W flight system and the standard Johnson B and V systems. We also describe the artificial star experiments which have been performed to derive the star count completeness in all the relevant branches of the color magnitude diagram. The entire photometric database and the completeness function will be made available on the Web immediately after the publication of the present paper.