ESA recently called for new "Science Ideas" to be investigated in terms of feasibility and technological developments -- for technologies not yet sufficiently mature. These ideas may in the future become candidates for M or L class missions within the ESA Science Program. With the launch of Gaia in December 2013, Europe entered a new era of space astrometry following in the footsteps of the very successful Hipparcos mission from the early 1990s. Gaia is the successor to Hipparcos, both of which operated in optical wavelengths, and Gaia is two orders of magnitude more accurate in the five astrometric parameters and is surveying four orders of magnitude more stars in a vast volume of the Milky Way. The combination of the Hipparcos/Tycho-2 catalogues with the first early Gaia data release will give improved proper motions over a long ~25 year baseline. The final Gaia solution will also establish a new optical reference frame by means of quasars, by linking the optical counterparts of radio (VLBI) sources defining the orientation of the reference frame, and by using the zero proper motion of quasars to determine a non-rotating frame. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. Traditionally, this problem is overcome by switching to the infra-red but this was not possible with Gaia's CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires that the CCDs operate in TDI mode, increasing their complexity.

Positions and proper motions of Gaia sources are expressed in a reference frame that ideally should be non-rotating relative to distant extragalactic objects, coincident with the International Celestial Reference System (ICRS), and consistent across all magnitudes. For sources fainter than 16th magnitude this is achieved thanks to Gaia's direct observations of quasars. At brighter magnitudes it is difficult to validate the quality of the reference frame due to the scarcity of comparison data. This paper examines the use of VLBI observations of radio stars to determine the spin and orientation of the bright reference frame of Gaia. Simultaneous estimation of the six spin and orientation parameters makes optimal use of VLBI data and makes it possible to include even single-epoch VLBI observations in the solution. The method is applied to Gaia Data Release 2 (DR2) using published VLBI data for 41 radio stars. Results for the 26 best-fitting sources indicate that the bright reference frame of Gaia DR2 is rotating relative to the faint quasars at a rate of about 0.1 mas/yr, significant at 2-sigma level. This supports a similar conclusion based on a comparison with stellar positions in the Hipparcos frame. The accuracy is currently limited by the small number of radio sources used, by uncertainties in the Gaia DR2 proper motions, and by the astrophysical nature of the radio stars. While the origin of the indicated rotation is understood and can be avoided in future data releases, it remains important to validate the bright reference frame of Gaia by independent observations. This can be achieved using VLBI astrometry, which may require re-observing the old sample of radio stars as well as measuring new objects. The unique historical value of positional measurements is stressed and VLBI observers are urged to ensure that relevant positional information is preserved for the future.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

Gaia is currently revolutionizing modern astronomy. However, much of the Galactic plane, center and the spiral arm regions are obscured by interstellar extinction, rendering them inaccessible because Gaia is an optical instrument. An all-sky near infrared (NIR) space observatory operating in the optical NIR, separated in time from the original Gaia would provide microarcsecond NIR astrometry and millimag photometry to penetrate obscured regions unraveling the internal dynamics of the Galaxy.

We report the discovery of the first likely black hole in a non-interacting binary system with a field red giant. By combining radial velocity measurements from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) with photometric variability data from the All-Sky Automated Survey for Supernovae (ASAS-SN), we identified the bright rapidly-rotating giant 2MASS J05215658+4359220 as a binary system with a massive unseen companion. Subsequent radial velocity measurements reveal a system with an orbital period of 83 days and near-zero eccentricity. The photometric variability period of the giant is consistent with the orbital period, indicative of star spots and tidal synchronization. Constraints on the giant's mass and radius from its luminosity, surface gravity, and temperature imply an unseen companion with mass of $3.3^{+2.8}_{-0.7}$ M$_\odot$, indicating a low-mass black hole or an exceedingly massive neutron star. Measurement of the astrometric binary motion by {\it Gaia} will further characterize the system. This discovery demonstrates the potential of massive spectroscopic surveys like APOGEE and all-sky, high-cadence photometric surveys like ASAS-SN to revolutionize our understanding of the compact object mass function, and to test theories of binary star evolution and the supernova mechanism.

A deep understanding of the Milky Way galaxy, its formation and evolution requires observations of huge numbers of stars. Stellar photometry, therefore, provides an economical method to obtain intrinsic stellar parameters. With the addition of distance information - deriving reliable ages from photometry is a possibility. We have developed a Bayesian method that generates 2D probability maps of a star's age and metallicity from photometry and parallax using isochrones. Our synthetic tests show that including a near-UV passband enables us to break the degeneracy between a star's age and metallicity for certain evolutionary stages. It is possible to find well-constrained ages and metallicities for turn-off and sub-giant stars with colours including a U band and a parallax with uncertainty less than ~20%. Metallicities alone are possible for the main sequence and giant branch. We find good agreement with the literature when we apply our method to the Gaia benchmark stars, particularly for turn-off and young stars. Further tests on the old open cluster NGC 188, however, reveal significant limitations in the stellar isochrones. The ages derived for the cluster stars vary with evolutionary stage, such that turn-off ages disagree with those on the sub-giant branch, and metallicities vary significantly throughout. Furthermore, the parameters vary appreciably depending on which colour combinations are used in the derivation. We identify the causes of these mismatches and show that improvements are needed in the modelling of giant branch stars and in the creation and calibration of synthetic near-UV photometry. Our results warn against applying isochrone fitting indiscriminately. In particular, the uncertainty on the stellar models should be quantitatively taken into account. Further efforts to improve the models will result in significant advancements in our ability to study the Galaxy.

In the era of large-scale surveys of stars in the Milky Way, stellar ages are crucial for studying the evolution of the Galaxy. But determining ages of field stars is notoriously difficult; therefore, we attempt to determine benchmark ages for the extensively studied Gaia benchmark stars which can be used for validation purposes. By searching the literature for age estimates from different methods and deriving new ages based on Bayesian isochrone fitting, we are able to put reliable limits on the ages of 16 out of the 33 benchmark stars. The giants with well-defined ages are all young, and an expansion of the sample to include older giants with asteroseismic ages would be beneficial. Some of the stars have surface parameters inconsistent with isochrones younger than 16 Gyr. Including $\alpha$-enhancement in the models when relevant resolves some of these cases, but others clearly highlight discrepancies between the models and observations. We test the impact of atomic diffusion on the age estimates by fitting to the actual surface metallicity of the models instead of the initial value and find that the effect is negligible except for a single turn-off star. Finally, we show that our ability to determine isochrone-based ages for large spectroscopic surveys largely mirrors our ability to determine ages for these benchmark stars, except for stars with $\log g \gtrsim 4.4$ dex since their location in the HR diagram is almost age insensitive. Hence, isochrone fitting does not constrain their ages given the typical uncertainties of spectroscopic stellar parameters.

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.

Context. Gaia DR2 offers unparalleled precision on stars' parallaxes and proper motions. This allows the prediction of microlensing events for which the lens stars (and any planets they possess) are nearby and may be well studied and characterised. Aims. We identify a number of potential microlensing events that will occur before the year 2035.5, 20 years from the Gaia DR2 reference epoch. Methods. We query Gaia DR2 for potential lenses within 100 pc, extract parallaxes and proper motions of the lenses and background sources, and identify potential lensing events. We estimate the lens masses from Priam effective temperatures, and use these to calculate peak magnifications and the size of the Einstein radii relative to the lens stars' habitable zones. Results. We identify 7 future events with a probability > 10% of an alignment within one Einstein radius. Of particular interest is DR2 5918299904067162240 (WISE J175839.20-583931.6), magnitude G = 14.9, which will lens a G = 13.9 background star in early 2030, with a median 23% net magnification. Other pairs are typically fainter, hampering characterisation of the lens (if the lens is faint) or the ability to accurately measure the magnification (if the source is much fainter than the lens). Of timely interest is DR2 4116504399886241792 (2MASS J17392440-2327071), which will lens a background star in July 2020, albeit with weak net magnification (0.03%). Median magnifications for the other 5 high-probability events range from 0.3% to 5.3%. The Einstein radii for these lenses are 1-10 times the radius of the habitable zone, allowing these lensing events to pick out cold planets around the ice line, and filling a gap between transit and current microlensing detections of planets around very low-mass stars. Conclusions. We provide a catalogue of the predicted events to aid future characterisation efforts... [abridged]

The goal of this paper is to demonstrate the outstanding quality of the second data release of the Gaia mission and its power for constraining many different aspects of the dynamics of the satellites of the Milky Way. We focus here on determining the proper motions of 75 Galactic globular clusters, nine dwarf spheroidal galaxies, one ultra-faint system, and the Large and Small Magellanic Clouds. Using data extracted from the Gaia archive, we derived the proper motions and parallaxes for these systems, as well as their uncertainties. We demonstrate that the errors, statistical and systematic, are relatively well understood. We integrated the orbits of these objects in three different Galactic potentials, and characterised their properties. We present the derived proper motions, space velocities, and characteristic orbital parameters in various tables to facilitate their use by the astronomical community. Our limited and straightforward analyses have allowed us for example to (i) determine absolute and very precise proper motions for globular clusters; (ii) detect clear rotation signatures in the proper motions of at least five globular clusters; (iii) show that the satellites of the Milky Way are all on high-inclination orbits, but that they do not share a single plane of motion; (iv) derive a lower limit for the mass of the Milky Way of 9.8^{+6.7}_{-2.7} x 10^{11} Msun based on the assumption that the Leo I dwarf spheroidal is bound; (v) derive a rotation curve for the Large Magellanic Cloud based solely on proper motions that is competitive with line-of-sight velocity curves, now using many orders of magnitude more sources; and (vi) unveil the dynamical effect of the bar on the motions of stars in the Large Magellanic Cloud. All these results highlight the incredible power of the Gaia astrometric mission, and in particular of its second data release.

Gaia Data Release 2 (Gaia DR2) contains results for 1693 million sources in the magnitude range 3 to 21 based on observations collected by the European Space Agency Gaia satellite during the first 22 months of its operational phase. We describe the input data, models, and processing used for the astrometric content of Gaia DR2, and the validation of these results performed within the astrometry task. Some 320 billion centroid positions from the pre-processed astrometric CCD observations were used to estimate the five astrometric parameters (positions, parallaxes, and proper motions) for 1332 million sources, and approximate positions at the reference epoch J2015.5 for an additional 361 million mostly faint sources. Special validation solutions were used to characterise the random and systematic errors in parallax and proper motion. For the sources with five-parameter astrometric solutions, the median uncertainty in parallax and position at the reference epoch J2015.5 is about 0.04 mas for bright (G<14 mag) sources, 0.1 mas at G=17 mag, and 0.7 mas at G=20 mag. In the proper motion components the corresponding uncertainties are 0.05, 0.2, and 1.2 mas/yr, respectively. The optical reference frame defined by Gaia DR2 is aligned with ICRS and is non-rotating with respect to the quasars to within 0.15 mas/yr. From the quasars and validation solutions we estimate that systematics in the parallaxes depending on position, magnitude, and colour are generally below 0.1 mas, but the parallaxes are on the whole too small by about 0.03 mas. Significant spatial correlations of up to 0.04 mas in parallax and 0.07 mas/yr in proper motion are seen on small (<1 deg) and intermediate (20 deg) angular scales. Important statistics and information for the users of the Gaia DR2 astrometry are given in the appendices.

The second release of Gaia data (Gaia DR2) contains the astrometric parameters for more than half a million quasars. This set defines a kinematically non-rotating reference frame in the optical domain referred to as the Gaia-CRF2. The Gaia-CRF2 is the first realisation of a non-rotating global optical reference frame that meets the ICRS prescriptions, meaning that it is built only on extragalactic sources. It consists of the positions of a sample of 556 869 sources in Gaia DR2, obtained from a positional cross-match with the ICRF3-prototype and AllWISE AGN catalogues. The sample constitutes a clean, dense, and homogeneous set of extragalactic point sources in the magnitude range G from 16 to 21 mag with accurately known optical positions. The median positional uncertainty is 0.12 mas for G < 18 mag and 0.5 mas at G = 20 mag. Large-scale systematics are estimated to be in the range 20 to 30 muas. The accuracy claims are supported by the parallaxes and proper motions of the quasars in Gaia DR2. The optical positions for a subset of 2820 sources in common with the ICRF3-prototype show very good overall agreement with the radio positions, but several tens of sources have significantly discrepant positions.

(abridged) We develop Bayesian methods and detection criteria for orbital fitting, and revise the detectability of exoplanets in light of the in-flight properties of Gaia. Limiting ourselves to one-planet systems as a first step of the development, we simulate Gaia data for exoplanet systems over a grid of S/N, orbital period, and eccentricity. The simulations are then fit using Markov chain Monte Carlo methods. We investigate the detection rate according to three information criteria and the delta chi^2. For the delta chi^2, the effective number of degrees of freedom depends on the mission length. We find that the choice of the Markov chain starting point can affect the quality of the results; we therefore consider two limit possibilities: an ideal case, and a very simple method that finds the starting point assuming circular orbits. Using Jeffreys' scale of evidence, the fraction of false positives passing a strong evidence criterion is < ~0.2% (0.6%) when considering a 5 yr (10 yr) mission and using the Akaike information criterion or the Watanabe-Akaike information criterion, and <0.02% (<0.06%) when using the Bayesian information criterion. We find that there is a 50% chance of detecting a planet with a minimum S/N=2.3 (1.7). This sets the maximum distance to which a planet is detectable to ~70 pc and ~3.5 pc for a Jupiter-mass and Neptune-mass planet, respectively, assuming a 10 yr mission, a 4 au semi-major axis, and a 1 M_sun star. The period is the orbital parameter that can be determined with the best accuracy, with a median relative difference between input and output periods of 4.2% (2.9%) assuming a 5 yr (10 yr) mission. The median accuracy of the semi-major axis of the orbit can be recovered with a median relative error of 7% (6%). The eccentricity can also be recovered with a median absolute accuracy of 0.07 (0.06).

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.

Determination of absolute parallaxes by means of a scanning astrometric satellite such as Hipparcos or Gaia relies on the short-term stability of the so-called basic angle between the two viewing directions. Uncalibrated variations of the basic angle may produce systematic errors in the computed parallaxes. We examine the coupling between a global parallax shift and specific variations of the basic angle, namely those related to the satellite attitude with respect to the Sun. The changes in observables produced by small perturbations of the basic angle, attitude, and parallaxes are calculated analytically. We then look for a combination of perturbations that has no net effect on the observables. In the approximation of infinitely small fields of view, it is shown that certain perturbations of the basic angle are observationally indistinguishable from a global shift of the parallaxes. If such perturbations exist, they cannot be calibrated from the astrometric observations but will produce a global parallax bias. Numerical simulations of the astrometric solution, using both direct and iterative methods, confirm this theoretical result. For a given amplitude of the basic angle perturbation, the parallax bias is smaller for a larger basic angle and a larger solar aspect angle. In both these respects Gaia has a more favourable geometry than Hipparcos. In the case of Gaia, internal metrology is used to monitor basic angle variations. Additionally, Gaia has the advantage of detecting numerous quasars, which can be used to verify the parallax zero point.

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.

As part of the data processing for Gaia Data Release~1 (Gaia DR1) a special astrometric solution was computed, the so-called auxiliary quasar solution. This gives positions for selected extragalactic objects, including radio sources in the second realisation of the International Celestial Reference Frame (ICRF2) that have optical counterparts bright enough to be observed with Gaia. A subset of these positions was used to align the positional reference frame of Gaia DR1 with the ICRF2. We describe the properties of the Gaia auxiliary quasar solution for a subset of sources matched to ICRF2, and compare their optical and radio positions at the sub-mas level. Their formal standard errors are better than 0.76~milliarcsec (mas) for 50% of the sources and better than 3.35~mas for 90%. Optical magnitudes are obtained in Gaia's unfiltered photometric G band. The comparison with the radio positions of the defining sources shows no systematic differences larger than a few tenths of a mas. The fraction of questionable solutions, not readily accounted for by the statistics, is less than 6%. Normalised differences have extended tails requiring case-by-case investigations for around 100 sources, but we have not seen any difference indisputably linked to an optical-radio offset in the sources.

The European Space Agency's Gaia satellite was launched into orbit around L2 in December 2013 with a payload containing 106 large-format scientific CCDs. The primary goal of the mission is to repeatedly obtain high-precision astrometric and photometric measurements of one thousand million stars over the course of five years. The scientific value of the down-linked data, and the operation of the onboard autonomous detection chain, relies on the high performance of the detectors. As Gaia slowly rotates and scans the sky, the CCDs are continuously operated in a mode where the line clock rate and the satellite rotation spin-rate are in synchronisation. Nominal mission operations began in July 2014 and the first data release is being prepared for release at the end of Summer 2016. In this paper we present an overview of the focal plane, the detector system, and strategies for on-orbit performance monitoring of the system. This is followed by a presentation of the performance results based on analysis of data acquired during a two-year window beginning at payload switch-on. Results for parameters such as readout noise and electronic offset behaviour are presented and we pay particular attention to the effects of the L2 radiation environment on the devices. The radiation-induced degradation in the charge transfer efficiency (CTE) in the (parallel) scan direction is clearly diagnosed; however, an extrapolation shows that charge transfer inefficiency (CTI) effects at end of mission will be approximately an order of magnitude less than predicted pre-flight. It is shown that the CTI in the serial register (horizontal direction) is still dominated by the traps inherent to the manufacturing process and that the radiation-induced degradation so far is only a few per cent. Finally, we summarise some of the detector effects discovered on-orbit which are still being investigated.

Gaia Data Release 1 (Gaia DR1) contains astrometric results for more than 1 billion stars brighter than magnitude 20.7 based on observations collected by the Gaia satellite during the first 14 months of its operational phase. We give a brief overview of the astrometric content of the data release and of the model assumptions, data processing, and validation of the results. For stars in common with the Hipparcos and Tycho-2 catalogues, complete astrometric single-star solutions are obtained by incorporating positional information from the earlier catalogues. For other stars only their positions are obtained by neglecting their proper motions and parallaxes. The results are validated by an analysis of the residuals, through special validation runs, and by comparison with external data. Results. For about two million of the brighter stars (down to magnitude ~11.5) we obtain positions, parallaxes, and proper motions to Hipparcos-type precision or better. For these stars, systematic errors depending e.g. on position and colour are at a level of 0.3 milliarcsecond (mas). For the remaining stars we obtain positions at epoch J2015.0 accurate to ~10 mas. Positions and proper motions are given in a reference frame that is aligned with the International Celestial Reference Frame (ICRF) to better than 0.1 mas at epoch J2015.0, and non-rotating with respect to ICRF to within 0.03 mas/yr. The Hipparcos reference frame is found to rotate with respect to the Gaia DR1 frame at a rate of 0.24 mas/yr. Based on less than a quarter of the nominal mission length and on very provisional and incomplete calibrations, the quality and completeness of the astrometric data in Gaia DR1 are far from what is expected for the final mission products. The results nevertheless represent a huge improvement in the available fundamental stellar data and practical definition of the optical reference frame.

The Gaia all-sky astrometric survey is challenged by several issues affecting the spacecraft stability. Amongst them, we find the focus evolution, straylight and basic angle variations Contrary to pre-launch expectations, the image quality is continuously evolving, during commissioning and the nominal mission. Payload decontaminations and wavefront sensor assisted refocuses have been carried out to recover optimum performance. Straylight and basic angle variations several orders of magnitude greater than foreseen were found and studied during commissioning by the Gaia scientists (payload experts). Building on their investigations, an ESA-Airbus DS working group was established during the early nominal mission and worked on a detailed root cause analysis. In parallel, Gaia scientists have also continued analysing the data, most notably comparing the BAM signal to global astrometric solutions, with remarkable agreement. In this contribution, a status review of these issues will be provided, with emphasis on the mitigation schemes and the lessons learned for future space missions where extreme stability is a key requirement.

Gaia's very accurate astrometric measurements will allow the International Celestial Reference Frame (ICRF) to be improved by a few orders of magnitude in the optical. Several sets of quasars are used to define a kinematically stable non-rotating reference frame with the barycentre of the Solar System as its origin. Gaia will also observe a large number of galaxies which could obtain accurate positions and proper motions although they are not point-like. The optical stability of the quasars is critical and we investigate how accurately the reference frame can be recovered. Various proper motion patterns are also present in the data, the best known is caused by the acceleration of the Solar System Barycentre, presumably, towards the Galactic centre. We review some other less-well-known effects that are not part of standard astrometric models. We model quasars and galaxies using realistic sky distributions, magnitudes and redshifts. Position variability is introduced using a Markov chain model. The reference frame is determined using the algorithm developed for the Gaia mission which also determines the acceleration of the Solar System. We also test a method to measure the velocity of the Solar System barycentre in a cosmological frame. We simulate the recovery of the reference frame and the acceleration of the Solar System and conclude that they are not significantly disturbed in the presence of quasar variability which is statistically averaged. However, the effect of a non-uniform sky distribution of the quasars can result in a correlation between the reference frame and acceleration which degrades the solution. Our results suggest that an attempt should be made to astrometrically determine the redshift dependent apparent drift of galaxies due to our velocity relative to the CMB, which in principle could allow the determination of the Hubble parameter.

Context. The Gaia project will determine positions, proper motions, and parallaxes for more than one billion stars in our Galaxy. It is known that Gaia's two telescopes are affected by a small but significant variation of the basic angle between them. Unless this variation is taken into account during data processing, e.g. using on-board metrology, it causes systematic errors in the astrometric parameters, in particular a shift of the parallax zero-point. Previously, we suggested an early reduction of Gaia data for the subset of Tycho-2 stars (Tycho-Gaia Astrometric Solution; TGAS). Aims. We aim to investigate whether quasars can be used to independently verify the parallax zero-point already in early data reductions. This is not trivially possible as the observation interval is too short to disentangle parallax and proper motion for the quasar subset. Methods. We repeat TGAS simulations but additionally include simulated Gaia observations of quasars from ground-based surveys. All observations are simulated with basic angle variations. To obtain a full astrometric solution for the quasars in TGAS we explore the use of prior information for their proper motions. Results. It is possible to determine the parallax zero-point for the quasars with a few {\mu}as uncertainty, and it agrees to a similar precision with the zero-point for the Tycho-2 stars. The proposed strategy is robust even for quasars exhibiting significant fictitious proper motion due to a variable source structure, or when the quasar subset is contaminated with stars misidentified as quasars. Conclusions. Using prior information about quasar proper motions we could provide an independent verification of the parallax zero-point in early solutions based on less than one year of Gaia data.

Gaia's astrometric solution aims to determine at least five parameters for each star, together with appropriate estimates of their uncertainties and correlations. This requires at least five distinct observations per star. In the early data reductions the number of observations may be insufficient for a five-parameter solution, and even after the full mission many stars will remain under-observed, including faint stars at the detection limit and transient objects. In such cases it is reasonable to determine only the two position parameters. Their formal uncertainties would however grossly underestimate the actual errors, due to the neglected parallax and proper motion. We aim to develop a recipe to calculate sensible formal uncertainties that can be used in all cases of under-observed stars. Prior information about the typical ranges of stellar parallaxes and proper motions is incorporated in the astrometric solution by means of Bayes' rule. Numerical simulations based on the Gaia Universe Model Snapshot (GUMS) are used to investigate how the prior influences the actual errors and formal uncertainties when different amounts of Gaia observations are available. We develop a criterion for the optimum choice of priors, apply it to a wide range of cases, and derive a global approximation of the optimum prior as a function of magnitude and galactic coordinates. The feasibility of the Bayesian approach is demonstrated through global astrometric solutions of simulated Gaia observations. With an appropriate prior it is possible to derive sensible positions with realistic error estimates for any number of available observations. Even though this recipe works also for well-observed stars it should not be used where a good five-parameter astrometric solution can be obtained without a prior. Parallaxes and proper motions from a solution using priors are always biased and should not be used.

The ESA Gaia mission uses two telescopes to create the most ambitious survey of the Galaxy. The angle between them must be known with exquisite precision and accuracy. An interferometer: the Basic Angle Monitoring system measures its variations. High quality data have been retrieved and analysed for more than a year. A summary of the in-orbit performance and some early results are presented

Context. The first release of astrometric data from Gaia will contain the mean stellar positions and magnitudes from the first year of observations, and proper motions from the combination of Gaia data with Hipparcos prior information (HTPM). Aims. We study the potential of using the positions from the Tycho-2 Catalogue as additional information for a joint solution with early Gaia data. We call this the Tycho-Gaia astrometric solution (TGAS). Methods. We adapt Gaia's Astrometric Global Iterative Solution (AGIS) to incorporate Tycho information, and use simulated Gaia observations to demonstrate the feasibility of TGAS and to estimate its performance. Results. Using six to twelve months of Gaia data, TGAS could deliver positions, parallaxes and annual proper motions for the 2.5 million Tycho-2 stars, with sub-milliarcsecond accuracy. TGAS overcomes some of the limitations of the HTPM project and allows its execution half a year earlier. Furthermore, if the parallaxes from Hipparcos are not incorporated in the solution, they can be used as a consistency check of the TGAS/HTPM solution.

We provide a revised assessment of the number of exoplanets that should be discovered by Gaia astrometry, extending previous studies to a broader range of spectral types, distances, and magnitudes. Our assessment is based on a large representative sample of host stars from the TRILEGAL Galaxy population synthesis model, recent estimates of the exoplanet frequency distributions as a function of stellar type, and detailed simulation of the Gaia observations using the updated instrument performance and scanning law. We use two approaches to estimate detectable planetary systems: one based on the S/N of the astrometric signature per field crossing, easily reproducible and allowing comparisons with previous estimates, and a new and more robust metric based on orbit fitting to the simulated satellite data. With some plausible assumptions on planet occurrences, we find that some 21,000 (+/-6000) high-mass (1-15M_J) long-period planets should be discovered out to distances of ~500pc for the nominal 5-yr mission (including at least 1000-1500 around M dwarfs out to 100pc), rising to some 70,000 (+/-20,000) for a 10-yr mission. We indicate some of the expected features of this exoplanet population, amongst them ~25-50 intermediate-period (P~2-3yr) transiting systems.

High-precision astrometric and radial-velocity observations require accurate modelling of stellar motions in order to extrapolate measurements over long time intervals, and to detect deviations from uniform motion caused for example by unseen companions. We aim to explore the simplest possible kinematic model of stellar motions, namely that of uniform rectilinear motion relative to the Solar System Barycentre, in terms of observable quantities including error propagation. The apparent path equation for uniform rectilinear motion is solved analytically in a classical (special-relativistic) framework, leading to rigorous expressions which relate the (apparent) astrometric parameters and radial velocity to the (true) kinematic parameters of the star in the barycentric reference system. We present rigorous and explicit formulae for the transformation of stellar positions, parallaxes, proper motions, and radial velocities from one epoch to another, assuming uniform rectilinear motion and taking into account light-time effects. The Jacobian matrix of the transformation is also given, allowing accurate and reversible propagation of errors over arbitrary time intervals. The light-time effects are generally very small but exceeds 0.1 mas or 0.1 m/s over 100 yr for at least 33 stars in the Hipparcos Catalogue. For high-velocity stars within a few tens of pc from the Sun light-time effects are generally more important than the effects of the curvature of their orbits in the Galactic potential.

The first release of astrometric data from Gaia is expected in 2016. It will contain the mean stellar positions and magnitudes from the first year of observations. For more than 100 000 stars in common with the Hipparcos Catalogue it will be possible to compute very accurate proper motions due to the time difference of about 24 years between the two missions. This Hundred Thousand Proper Motions (HTPM) project will be part of the first release. Our aim is to investigate how early Gaia data can be optimally combined with information from the Hipparcos Catalogue in order to provide the most accurate and reliable results for HTPM. The Astrometric Global Iterative Solution (AGIS) was developed to compute the astrometric core solution based on the Gaia observations and will be used for all releases of astrometric data from Gaia. We adapt AGIS to process Hipparcos data in addition to Gaia observations, and use simulations to verify and study the joint solution method. For the HTPM stars we predict proper motion accuracies between 14 and 134 muas/yr, depending on stellar magnitude and amount of Gaia data available. Perspective effects will be important for a significant number of HTPM stars, and in order to treat these effects accurately we introduce a scaled model of kinematics. We define a goodness-of-fit statistic which is sensitive to deviations from uniform space motion, caused for example by binaries with periods of 10-50 years. HTPM will significantly improve the proper motions of the Hipparcos Catalogue well before highly accurate Gaia- only results become available. Also, HTPM will allow us to detect long period binary and exoplanetary candidates which would be impossible to detect from Gaia data alone. The full sensitivity will not be reached with the first Gaia release but with subsequent data releases. Therefore HTPM should be repeated when more Gaia data become available.

Starting in 2013, Gaia will deliver highly accurate astrometric data, which eventually will supersede most other stellar catalogues in accuracy and completeness. It is, however, lim- ited to observations from magnitude 6 to 20 and will therefore not include the brightest stars. Nano-JASMINE, an ultrasmall Japanese astrometry satellite, will observe these bright stars, but with much lower accuracy. Hence, the Hipparcos catalogue from 1997 will likely remain the main source of accurate distances to bright nearby stars. We are investigating how this might be improved by optimally combining data from all three missions in a joint astrometric solu- tion. This would take advantage of the unique features of each mission: the historic bright-star measurements of Hipparcos, the updated bright-star observations of Nano-JASMINE, and the very accurate reference frame of Gaia. The long temporal baseline between the missions pro- vides additional benefits for the determination of proper motions and binary detection, which indirectly improve the parallax determination further. We present a quantitative analysis of the expected gains based on simulated data for all three missions.

The Gaia payload ensures maximum passive stability using a single material, SiC, for most of its elements. Dedicated metrology instruments are, however, required to carry out two functions: monitoring the basic angle and refocusing the telescope. Two interferometers fed by the same laser are used to measure the basic angle changes at the level of $\mu$as (prad, micropixel), which is the highest level ever achieved in space. Two Shack-Hartmann wavefront sensors, combined with an ad-hoc analysis of the scientific data are used to define and reach the overall best-focus. In this contribution, the systems, data analysis, procedures and performance achieved during commissioning are presented

While the precise relationship between the Milky Way disk and the symmetry planes of the dark matter halo remains somewhat uncertain, a time-varying disk orientation with respect to an inertial reference frame seems probable. Hierarchical structure formation models predict that the dark matter halo is triaxial and tumbles with a characteristic rate of ~2 rad/Hubble time (~30 muas/yr). These models also predict a time-dependent accretion of gas, such that the angular momentum vector of the disk should be misaligned with that of the halo. These effects, as well as tidal effects of the LMC, will result in the rotation of the angular momentum vector of the disk population with respect to the quasar reference frame. We assess the accuracy with which the positions and proper motions from Gaia can be referred to a kinematically non-rotating system, and show that the spin vector of the transformation from any rigid self-consistent catalog frame to the quasi-inertial system defined by quasars should be defined to better than 1 muas/yr. Determination of this inertial frame by Gaia will reveal any signature of the disk orientation varying with time, improve models of the potential and dynamics of the Milky Way, test theories of gravity, and provide new insights into the orbital evolution of the Sagittarius dwarf galaxy and the Magellanic Clouds.

A number of large spectroscopic surveys of stars in the Milky Way are under way or are being planned. In this context it is important to discuss the extent to which elemental abundances can be used as discriminators between different (known and unknown) stellar populations in the Milky Way. We aim to establish the requirements in terms of precision in elemental abundances, as derived from spectroscopic surveys of the Milky Way's stellar populations, in order to detect interesting substructures in elemental abundance space. We present a simple relation between the minimum number of stars needed to detect a given substructure and the precision of the measurements. The results are in agreement with recent small- and large-scale studies, with high and low precision, respectively. Large-number statistics cannot fully compensate for low precision in the abundance measurements and each survey should carefully evaluate what the main science drivers are for the survey and ensure that the chosen observational strategy will result in the precision necessary to answer the questions posed.

Context. Kinematical data such as the mean velocities and velocity dispersions of stellar samples are useful tools to study galactic structure and evolution. However, observational data are often incomplete (e.g., lacking the radial component of the motion) and may have significant observational errors. For example, the majority of faint stars observed with Gaia will not have their radial velocities measured. Aims. Our aim is to formulate and test a new maximum likelihood approach to estimating the kinematical parameters for a local stellar sample when only the transverse velocities are known (from parallaxes and proper motions). Methods. Numerical simulations using synthetically generated data as well as real data (based on the Geneva-Copenhagen survey) are used to investigate the statistical properties (bias, precision) of the method, and to compare its performance with the much simpler "projection method" described by Dehnen & Binney (1998). Results. The maximum likelihood method gives more correct estimates of the dispersion when observational errors are important, and guarantees a positive-definite dispersion matrix, which is not always obtained with the projection method. Possible extensions and improvements of the method are discussed.

The Hipparcos mission (1989-1993) resulted in the first space-based stellar catalogue including measurements of positions, parallaxes and annual proper motions accurate to about one milli-arcsecond. More space astrometry missions will follow in the near future. The ultra-small Japanese mission Nano-JASMINE (launch in late 2013) will determine positions and annual proper motions with some milli-arcsecond accuracy. In mid 2013 the next-generation ESA mission Gaia will deliver some tens of micro-arcsecond accurate astrometric parameters. Until the final Gaia catalogue is published in early 2020 the best way of improving proper motion values is the combination of positions from different missions separated by long time intervals. Rather than comparing positions from separately reduced catalogues, we propose an optimal method to combine the information from the different data sets by making a joint astrometric solution. This allows to obtain good results even when each data set alone is insufficient for an accurate reduction. We demonstrate our method by combining Hipparcos and simulated Nano-JASMINE data in a joint solution. We show a significant improvement over the conventional catalogue combination.

The ESA space astrometry mission Gaia, planned to be launched in 2013, has been designed to make angular measurements on a global scale with micro-arcsecond accuracy. A key component of the data processing for Gaia is the astrometric core solution, which must implement an efficient and accurate numerical algorithm to solve the resulting, extremely large least-squares problem. The Astrometric Global Iterative Solution (AGIS) is a framework that allows to implement a range of different iterative solution schemes suitable for a scanning astrometric satellite. In order to find a computationally efficient and numerically accurate iteration scheme for the astrometric solution, compatible with the AGIS framework, we study an adaptation of the classical conjugate gradient (CG) algorithm, and compare it to the so-called simple iteration (SI) scheme that was previously known to converge for this problem, although very slowly. The different schemes are implemented within a software test bed for AGIS known as AGISLab, which allows to define, simulate and study scaled astrometric core solutions. After successful testing in AGISLab, the CG scheme has been implemented also in AGIS. The two algorithms CG and SI eventually converge to identical solutions, to within the numerical noise (of the order of 0.00001 micro-arcsec). These solutions are independent of the starting values (initial star catalogue), and we conclude that they are equivalent to a rigorous least-squares estimation of the astrometric parameters. The CG scheme converges up to a factor four faster than SI in the tested cases, and in particular spatially correlated truncation errors are much more efficiently damped out with the CG scheme.

The Gaia satellite will observe about one billion stars and other point-like sources. The astrometric core solution will determine the astrometric parameters (position, parallax, and proper motion) for a subset of these sources, using a global solution approach which must also include a large number of parameters for the satellite attitude and optical instrument. The accurate and efficient implementation of this solution is an extremely demanding task, but crucial for the outcome of the mission. We provide a comprehensive overview of the mathematical and physical models applicable to this solution, as well as its numerical and algorithmic framework. The astrometric core solution is a simultaneous least-squares estimation of about half a billion parameters, including the astrometric parameters for some 100 million well-behaved so-called primary sources. The global nature of the solution requires an iterative approach, which can be broken down into a small number of distinct processing blocks (source, attitude, calibration and global updating) and auxiliary processes (including the frame rotator and selection of primary sources). We describe each of these processes in some detail, formulate the underlying models, from which the observation equations are derived, and outline the adopted numerical solution methods with due consideration of robustness and the structure of the resulting system of equations. Appendices provide brief introductions to some important mathematical tools (quaternions and B-splines for the attitude representation, and a modified Cholesky algorithm for positive semidefinite problems) and discuss some complications expected in the real mission data.

The Gaia mission has been designed to perform absolute astrometric measurements with unprecedented accuracy; the end-of-mission parallax standard error is required to be 30 micro-arcseconds for a G2V type star of magnitude 15. These requirements set a stringent constraint on the accuracy of the estimation of the location of the stellar image on the CCD for each observation: e.g., 0.3 milli-arseconds (mas) or 0.005 pixels for the same V=15 G2V star. However the Gaia CCDs will suffer from charge transfer inefficiency (CTI) caused by radiation damage that will degrade the stellar image quality and may degrade the astrometric performance of Gaia if not properly addressed. For the first time at this level of detail, the potential impact of radiation damage on the performance of Gaia is investigated. In this first paper we focus on the evaluation of the CTI impact on the image location accuracy. We show that CTI decreases the stellar image signal-to-noise ratio and irreversibly degrades the image location estimation precision. As a consequence the location estimation standard errors increase by up to 6% for a radiation damage level equivalent to the end-of-mission. In addition the CTI-induced image distortion introduces a systematic bias in the image location estimation (up to 0.05 pixels or 3 mas in the Gaia operating conditions). We present a novel approach to CTI mitigation that enables, without correction of the raw data, the unbiased estimation of the image location and flux from damaged observations. Its implementation reduces the maximum measured location bias for the faintest magnitude to 0.005 pixels (~4e-4 pixels at magnitude 15). In a second paper we will investigate how the CTI effects affect the final astrometric accuracy of Gaia by propagating residual errors through the astrometric solution.

This paper provides a description of the Java software framework which has been constructed to run the Astrometric Global Iterative Solution for the Gaia mission. This is the mathematical framework to provide the rigid reference frame for Gaia observations from the Gaia data itself. This process makes Gaia a self calibrated, and input catalogue independent, mission. The framework is highly distributed typically running on a cluster of machines with a database back end. All code is written in the Java language. We describe the overall architecture and some of the details of the implementation.

Context. The Gaia satellite will measure highly accurate absolute parallaxes of hundreds of millions of stars by comparing the parallactic displacements in the two fields of view of the optical instrument. The requirements on the stability of the 'basic angle' between the two fields are correspondingly strict, and possible variations (on the microarcsec level) are therefore monitored by an on-board metrology system. Nevertheless, since even very small periodic variations of the basic angle might cause a global offset of the measured parallaxes, it is important to find independent verification methods. Aims. We investigate the potential use of Galactic Cepheids as standard candles for verifying the Gaia parallax zero point. Methods. We simulate the complete population of Galactic Cepheids and their observations by Gaia. Using the simulated data, simultaneous fits are made of the parameters of the period-luminosity relation and a global parallax zero point. Results. The total number of Galactic Cepheids is estimated at about 20 000, of which nearly half could be observed by Gaia. In the most favourable circumstances, including negligible intrinsic scatter and extinction errors, the determined parallax zero point has an uncertainty of 0.2 microarcsec. With more realistic assumptions the uncertainty is several times larger, and the result is very sensitive to errors in the applied extinction corrections. Conclusions. The use of Galactic Cepheids alone will not be sufficient to determine a possible parallax zero-point error to the full potential systematic accuracy of Gaia. The global verification of Gaia parallaxes will most likely depend on a combination of many different methods, including this one.

Current optical space telescopes rely upon silicon Charge Coupled Devices (CCDs) to detect and image the incoming photons. The performance of a CCD detector depends on its ability to transfer electrons through the silicon efficiently, so that the signal from every pixel may be read out through a single amplifier. This process of electron transfer is highly susceptible to the effects of solar proton damage (or non-ionizing radiation damage). This is because charged particles passing through the CCD displace silicon atoms, introducing energy levels into the semi-conductor bandgap which act as localized electron traps. The reduction in Charge Transfer Efficiency (CTE) leads to signal loss and image smearing. The European Space Agency's astrometric Gaia mission will make extensive use of CCDs to create the most complete and accurate stereoscopic map to date of the Milky Way. In the context of the Gaia mission CTE is referred to with the complementary quantity Charge Transfer Inefficiency (CTI = 1-CTE). CTI is an extremely important issue that threatens Gaia's performances. We present here a detailed Monte Carlo model which has been developed to simulate the operation of a damaged CCD at the pixel electrode level. This model implements a new approach to both the charge density distribution within a pixel and the charge capture and release probabilities, which allows the reproduction of CTI effects on a variety of measurements for a large signal level range in particular for signals of the order of a few electrons. A running version of the model as well as a brief documentation and a few examples are readily available at http://www.strw.leidenuniv.nl/~prodhomme/cemga.php as part of the CEMGA java package (CTI Effects Models for Gaia).

We investigate whether it is possible that viable microbes could have been transported to Earth from the planets in extra-solar systems by means of natural vehicles such as ejecta expelled by comet or asteroid impacts on such planets. The probabilities of close encounters with other solar systems are taken into account as well as the limitations of bacterial survival times inside ejecta in space, caused by radiation and DNA decay. The conclusion is that no potentially DNA/RNA life-carrying ejecta from another solar system in the general Galactic star field landed on Earth before life already existed on Earth, not even if microbial survival time in space is as long as tens of millions of years. However, if the Sun formed initially as a part of a star cluster, as is commonly assumed, we cannot rule out the possibility of transfer of life from one of the sister systems to us. Likewise, there is a possibility that some extra-solar planets carry life that originated in our solar system. It will be of great interest to identify the members of the Sun's birth cluster of stars and study them for evidence for planets and life on the planets. The former step may be accomplished by the GAIA mission, the latter step by the SIM and DARWIN missions. Therefore it may not be too long until we have experimental knowledge on the question whether the natural transfer of life from one solar system to another has actually taken place.

We investigate the astrometric effects of stellar surface structures as a practical limitation to ultra-high-precision astrometry, e.g. in the context of exoplanet searches, and to quantify the expected effects in different regions of the HR-diagram. Stellar surface structures are likely to produce fluctuations in the integrated flux and radial velocity of the star, as well as a variation of the observed photocentre, i.e. astrometric jitter, and closure phase. We use theoretical considerations supported by Monte Carlo simulations to derive statistical relations between the corresponding astrometric, photometric, and radial-velocity effects. For most stellar types the astrometric jitter due to stellar surface structures is expected to be of order 10 micro-AU or greater. This is more than the astrometric displacement typically caused by an Earth-size exoplanet in the habitable zone, which is about 1-4 micro-AU for long-lived main-sequence stars. Only for stars with extremely low photometric variability (<0.5 mmag) and low magnetic activity, comparable to that of the Sun, will the astrometric jitter be of order 1 micro-AU, suffcient to allow the astrometric detection of an Earth-sized planet in the habitable zone. While stellar surface structure may thus seriously impair the astrometric detection of small exoplanets, it has in general negligible impact on the detection of large (Jupiter-size) planets and on the determination of stellar parallax and proper motion. From the starspot model we also conclude that the commonly used spot filling factor is not the most relevant parameter for quantifying the spottiness in terms of the resulting astrometric, photometric and radial-velocity variations.

Gaia is ESA's ambitious space astrometry mission the main objective of which is to astrometrically and spectro-photometrically map 1000 Million celestial objects (mostly in our galaxy) with unprecedented accuracy. The announcement of opportunity for the data processing will be issued by ESA late in 2006. The Gaia Data Processing and Analysis Consortium (DPAC) has been formed recently and is preparing an answer. The satellite will downlink close to 100 TB of raw telemetry data over 5 years. To achieve its required accuracy of a few 10s of Microarcsecond astrometry, a highly involved processing of this data is required. In addition to the main astrometric instrument Gaia will host a Radial Velocity instrument, two low-resolution dispersers for multi-color photometry and two Star Mappers. Gaia is a flying Giga Pixel camera. The various instruments each require relatively complex processing while at the same time being interdependent. We describe the overall composition of the DPAC and the envisaged overall architecture of the Gaia data processing system. We shall delve further into the core processing - one of the nine, so-called, coordination units comprising the Gaia processing system.

(Abridged) We investigate the effects of metallicity on the broad-band photometric colors of late-type giants, and make a comparison of synthetic colors with observed photometric properties of late-type giants over a wide range of effective temperatures (T_eff=3500-4800 K) and gravities (log g=0.0-2.5), at [M/H]=-1.0 and -2.0. The influence of metallicity on the synthetic photometric colors is generally small at effective temperatures above \~3800 K, but the effects grow larger at lower T_eff, due to the changing efficiency of molecule formation which reduces molecular opacities at lower [M/H]. To make a detailed comparison of the synthetic and observed photometric colors of late type giants in the T_eff--color and color--color planes, we derive a set of new T_eff--log g--color relations based on synthetic photometric colors, at [M/H]=-0.5, -1.0, -1.5, and -2.0. While differences between the new T_eff--color relations and those available from the literature are typically well within ~100 K, effective temperatures predicted by the scales based on synthetic colors tend to be slightly higher than those resulting from the T_eff--color relations based on observations, with the offsets up to ~100 K. This is clearly seen both at [M/H]=-1.0 and -2.0, especially in the T_eff--(B-V) and T_eff--(V-K) planes. The consistency between T_eff--log g--color scales based on synthetic colors calculated with different stellar atmosphere codes is very good, with typical differences being well within \Delta T_eff~70 K at [M/H]=-1.0 and \Delta T_eff~40 K at [M/H]=-2.0.

To assess the current status in the theoretical modeling of the spectral properties of late-type giants, we provide a comparison of synthetic photometric colors of late-type giants (calculated with PHOENIX, MARCS and ATLAS model atmospheres) with observations, at [M/H]=0.0 and -2.0. Overall, there is a good agreement between observed and synthetic colors, and synthetic colors and published Teff-color relations, both at [M/H]=0.0 and -2.0. Deviations from the observed trends in Teff-color planes are generally within \pm 150K (or less) in the effective temperature range Teff=3500-4800K. Synthetic colors calculated with different stellar atmosphere models typically agree to ~100K, within a large range of effective temperatures and gravities. Some discrepancies are seen in the Teff-(B-V) plane below Teff~3800K at [M/H]=0.0, due to difficulties in reproducing the 'turn-off' to the bluer colors which is seen in the observed data at Teff~3600K. Note that at [M/H]=-2.0 effective temperatures given by the scale of Alonso et al. (1999) are generally lower than those resulting from other Teff-color relations based both on observed and synthetic colors.

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the proven principles of ESA's Hipparcos mission to create an all-sky survey of about one billion stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. Through its massive measurement of stellar distances, motions and multi-colour photometry it will provide fundamental data necessary for unravelling the structure, formation and evolution of the Galaxy. This paper presents the design and performance of the broad- and medium-band set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands (extending from the ultraviolet to the far-red), the criteria, and the methodology on which this choice has been based are discussed in detail. We analyse the photometric capabilities for characterizing the luminosity, temperature, gravity and chemical composition of stars. We also discuss the automatic determination of these physical parameters for the large number of observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. Finally, the capability of the photometric system to deal with the main Gaia science case is outlined.

We present new synthetic broad-band photometric colors for late-type giants based on synthetic spectra calculated with the PHOENIX model atmosphere code. The grid covers effective temperatures T_eff=3000...5000K, gravities log g=-0.5...+3.5, and metallicities [M/H]=+0.5...-4.0. We show that individual broad-band photometric colors are strongly affected by model parameters such as molecular opacities, gravity, microturbulent velocity, and stellar mass. Our exploratory 3D modeling of a prototypical late-type giant shows that convection has a noticeable effect on the photometric colors too, as it alters significantly both the vertical and horizontal thermal structures in the outer atmosphere. The differences between colors calculated with full 3D hydrodynamical and 1D model atmospheres are significant (e.g., \Delta(V-K)~0.2 mag), translating into offsets in effective temperature of up to ~70K. For a sample of 74 late-type giants in the Solar neighborhood, with interferometric effective temperatures and broad-band photometry available in the literature, we compare observed colors with a new PHOENIX grid of synthetic photometric colors, as well as with photometric colors calculated with the MARCS and ATLAS model atmosphere codes. (abridged)

Late type giants (RGB/AGB stars) will be important tracers of the Galactic morphology and evolution in the framework of Gaia, as they are intrinsically bright and thus can probe distant stellar populations or those obscured by interstellar extinction. A realistic representation of their atmospheres and spectra with stellar atmosphere models is thus of crucial importance, both for the design and optimization of Gaia instruments, as well as the interpretation of provided astrophysical data. Our analysis of synthetic photometric colors of late-type giants based on PHOENIX, MARCS and ATLAS model atmospheres indicates a general agreement between the current theoretical predictions and observations in the framework of stationary 1-D model atmospheres. Presently available models allow temperature determinations of RGB/AGB stars to an accuracy of ~\pm100 K. In an exploratory study we try to quantify possible residual systematic effects due to the approximations made in 1-D models using full 3-D hydrodynamical models. We find that differences in broad-band photometric colors calculated with 1-D and 3-D models are significant, translating to the offsets in effective temperature of up to ~70 K. Clearly, full 3-D hydrodynamical models will help to alleviate such ambiguities in current theoretical modeling. Additionally, they will allow to study new phenomena, to open qualitatively new windows for stellar astrophysics in the Gaia-era.

Gaia will play an important role in providing information about star formation histories, merging events, intergalactic streams etc., for nearby galaxies of the Local Group. One of the most crucial contributions will be proper motions, especially for stars in the outermost parts of the galaxies, obtainable for stellar populations to ~150 kpc with RGB stars. Together with radial velocities for the brightest giants <80 kpc, this will provide membership information for individual stars and global kinematical picture of the most nearby galaxies, including the Magellanic Clouds (MCs). Gaia will also provide photometric metallicities (\sigma([M/H])<0.3) for individual giants and/or supergiants in dwarf galaxies to ~200 kpc. MSTO ages will be possible for the youngest stellar populations in the most nearby galaxies (e.g., MCs), whereas stars on RGB/AGB may provide age estimates for populations to ~150 kpc. Gaia will allow to study the outermost parts of the galaxies, which (because of their large spatial extent) are difficult to assess from the ground. Apart from allowing to clarify the structure and evolution of the dwarf galaxies, this will also make it possible to investigate galactic tidal debris, thus providing additional details for the global picture of formation and evolution of the Milky Way Galaxy.

Spectral-line displacements away from the wavelengths naively expected from the Doppler shift caused by stellar radial motion may originate as convective shifts (correlated velocity and brightness patterns in the photosphere), as gravitational redshifts, or perhaps be induced by wave motions. Absolute lineshifts, in the past studied only for the Sun, are now accessible also for other stars thanks to astrometric determination of stellar radial motion, and spectrometers with accurate wavelength calibration. Comparisons between spectroscopic apparent radial velocities and astrometrically determined radial motions reveal greater spectral blueshifts in F-type stars than in the Sun (as theoretically expected from their more vigorous convection), further increasing in A-type stars (possibly due to atmospheric shockwaves). An important near-future development to enable a further analysis of stellar surface structure will be the study of wavelength variations across spatially resolved stellar disks, e.g., the center-to-limb wavelength changes along a stellar diameter, and their spatially resolved time variability.

Spectroscopic radial velocities for several nearby open clusters suggest that spectra of (especially earlier-type) rapidly rotating stars are systematically blueshifted by 3 km/s or more, relative to the spectra of slowly rotating ones. Comparisons with astrometrically determined radial motions in the Hyades suggests this to be an absolute blueshift, relative to wavelengths naively expected from stellar radial motion and gravitational redshift. Analogous trends are seen also in most other clusters studied (Pleiades, Coma Berenices, Praesepe, Alpha Persei, IC 2391, NGC 6475, IC 4665, NGC 1976 and NGC 2516). Possible mechanisms are discussed, including photospheric convection, stellar pulsation, meridional circulation, and shock-wave propagation, as well as effects caused by template mismatch in determining wavelength displacements. For early-type stars, a plausible mechanism is shock-wave propagation upward through the photospheric line-forming regions. Such wavelength shifts thus permit studies of certain types of stellar atmospheric dynamics and - irrespective of their cause - may influence deduced open-cluster membership (when selected from common velocity) and deduced cluster dynamics (some types of stars might show fortuitous velocity patterns).

We discuss the tracing of star formation histories with ESA's space astrometry mission GAIA, emphasizing the advantages of AGB stars for this purpose. GAIA's microarcsecond-level astrometry, multi-band photometry and spectroscopy will provide individual distances, motions, effective temperatures, gravities and metallicities for vast numbers of AGB stars in the Galaxy and beyond. Reliable ages of AGB stars can be determined to distances of \~200 kpc in a wide range of ages and metallicities, allowing star formation histories to be studied in a diversity of astrophysical environments.

This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were adopted during the 24th General Assembly in Manchester, 2000 and provides details and explanations for these Resolutions. It is explained why they present significant progress over the corresponding IAU 1991 Resolutions and why they are necessary in the light of present accuracies in astrometry, celestial mechanics and metrology. In fact most of these Resolutions are consistent with astronomical models and software already in use. The metric tensors and gravitational potentials of both the Barycentric Celestial Reference System and Geocentric Celestial Reference System are defined and discussed. The necessity and relevance of the two celestial reference systems are explained. The transformations of coordinates and gravitational potentials are discussed. Potential coefficients parameterizing the post-Newtonian gravitational potentials are expounded. Simplified versions of the time transformations suitable for modern clock accuracies are elucidated. Various approximations used in the Resolutions are explicated and justified. Some models (e.g. for higher spin moments) that serve the purpose for estimating orders of magnitude have actually never been published before.

Accuracy levels of metres per second require the fundamental concept of 'radial velocity' for stars and other distant objects to be examined, both as a physical velocity, and as measured by spectroscopic and astrometric techniques. Aiming at definitions that are unambiguous at accuracy levels of 1 m/s, we analyse different concepts of radial velocity and their interrelations. We distinguish between the geometrically defined 'kinematic radial velocity' and 'astrometric radial velocity' on one hand, and the spectroscopically defined 'barycentric radial-velocity measure' on the other. The last two concepts are defined by recent resolutions adopted by the International Astronomical Union (IAU), the motives and consequences of which are explained in this paper.

We argue that tracing star formation histories with ESA's space mission GAIA using main sequence turn-off (MSTO) point dating will mainly be effective in cases of mild interstellar extinction (E(B-V)<0.5). For higher reddenings the MSTO approach will be severely limited both in terms of traceable ages (t<0.5 Gyr at 8 kpc; E(B-V)=1.0) and/or distances (d<2 kpc for t=15 Gyr; E(B-V)=1.0), since the MSTO will be located at magnitudes too faint for GAIA. AGB stars may alternatively provide precise population ages with GAIA for a wide range of ages and metallicities, with traceable distances of up to 250 kpc at E(B-V)=0 (d=15 kpc if E(B-V)=2.0). It is essential however that effective temperatures precise to 0.01 dex, metallicities to 0.2 dex, and E(B-V) to 0.03 are derived for individual stars, in order to obtain their ages precise to 0.2 dex. This task is quite challenging for GAIA photometry and spectroscopy, though preliminary tests show that comparable precisions may be achieved with GAIA medium band photometry.

Radial motions of stars in nearby moving clusters are determined from accurate proper motions and trigonometric parallaxes, without any use of spectroscopy. Assuming that cluster members share the same velocity vector (apart from a random dispersion), we apply a maximum-likelihood method on astrometric data from Hipparcos to compute radial and space velocities (and their dispersions) in the Ursa Major, Hyades, Coma Berenices, Pleiades, and Praesepe clusters, and for the alpha Persei, Scorpius-Centaurus, and `HIP 98321' associations. The radial motion of the Hyades cluster is determined to within 0.47 km/s (standard error), and that of its individual stars to within 0.6 km/s. For other clusters, Hipparcos data yield astrometric radial velocities with typical accuracies of a few km/s. A comparison of these astrometric values with spectroscopic radial velocities in the literature shows a good general agreement and, in the case of the best-determined Hyades cluster, also permits searches for subtle astrophysical differences, such as evidence for enhanced convective blueshifts of F-dwarf spectra, and decreased gravitational redshifts in giants. Similar comparisons for the Scorpius OB2 complex indicate some expansion of its associations, albeit slower than expected from their ages. As a by-product from the radial-velocity solutions, kinematically improved parallaxes for individual stars are obtained, enabling Hertzsprung-Russell diagrams with unprecedented accuracy in luminosity. For the Hyades (parallax accuracy 0.3 mas), its main sequence resembles a thin line, possibly with wiggles in it.

The GAIA astrometric mission has recently been approved as one of the next two `cornerstones' of ESA's science programme, with a launch date target of not later than mid-2012. GAIA will provide positional and radial velocity measurements with the accuracies needed to produce a stereoscopic and kinematic census of about one billion stars throughout our Galaxy (and into the Local Group), amounting to about 1 per cent of the Galactic stellar population. GAIA's main scientific goal is to clarify the origin and history of our Galaxy, from a quantitative census of the stellar populations. It will advance questions such as when the stars in our Galaxy formed, when and how it was assembled, and its distribution of dark matter. The survey aims for completeness to V=20 mag, with accuracies of about 10 microarcsec at 15 mag. Combined with astrophysical information for each star, provided by on-board multi-colour photometry and (limited) spectroscopy, these data will have the precision necessary to quantify the early formation, and subsequent dynamical, chemical and star formation evolution of our Galaxy. Additional products include detection and orbital classification of tens of thousands of extra-Solar planetary systems, and a comprehensive survey of some 10^5-10^6 minor bodies in our Solar System, through galaxies in the nearby Universe, to some 500,000 distant quasars. It will provide a number of stringent new tests of general relativity and cosmology. The complete satellite system was evaluated as part of a detailed technology study, including a detailed payload design, corresponding accuracy assesments, and results from a prototype data reduction development.

The Hipparcos Transit Data are a collection of partially reduced, fully calibrated observations of (mostly) double and multiple stars obtained with the ESA Hipparcos astrometry satellite. The data are publicly available, as part of the CD-ROM set distributed with the Hipparcos and Tycho Catalogues (ESA SP--1200, 1997), for about a third of the Hipparcos Catalogue entries including all confirmed or suspected non-single stars. The Transit Data consist of signal modulation parameters derived from the individual transits of the targets across the Hipparcos focal grid. The Transit Data permit re-reduction of the satellite data for individual objects, using arbitrarily complex object models in which time-variable photometric as well as geometric characteristics may be taken into account. We describe the structure and contents of the Transit Data files and give examples of how the data can be used. Some of the applications use standard astronomical software: Difmap or AIPS for aperture synthesis imaging, and GaussFit for detailed model fitting. Fortran code converting the data into formats suitable for these application programs has been made public in order to encourage and facilitate the use of Hipparcos Transit Data.

High-accuracy astrometry permits the determination of not only stellar tangential motion, but also the component along the line-of-sight. Such non-spectroscopic (i.e. astrometric) radial velocities are independent of stellar atmospheric dynamics, spectral complexity and variability, as well as of gravitational redshift. Three methods are analysed: (1) changing annual parallax, (2) changing proper motion and (3) changing angular extent of a moving group of stars. All three have significant potential in planned astrometric projects. Current accuracies are still inadequate for the first method, while the second is marginally feasible and is here applied to 16 stars. The third method reaches high accuracy (<1 km/s) already with present data, although for some clusters an accuracy limit is set by uncertainties in the cluster expansion rate.

GAIA is a short-listed candidate for the ESA Cornerstone mission C5, meeting the ESA Survey Committee requirement for an observatory mission, dedicated to astrometry, providing 10 micro-arcsecond accuracy at 15th magnitude. The GAIA mission concept follows the dramatic success of the ESA HIPPARCOS mission, utilising a continuously scanning spacecraft, accurately measuring 1-D coordinates along great circles, in (at least) two simultaneous fields of view, separated by a known angle. These 1-D relative coordinates are later converted to the five astrometric parameters of position and motions in a global analysis. GAIA will provide precise astrometry and multi-colour photometry for all the one billion stars, quasars, and compact galaxies to I=20 on the sky. GAIA will additionally provide the sixth phase-space parameter, radial velocity, from a slitless spectroscopic survey of most stars brighter than about magnitude 17. The technical challenges are considerable, but achievable. The scientific returns are spectacular, with greatest impact in the study of stellar populations and dynamical structure of the galaxies of our Local Group, and in providing the first complete census of the stars and massive planets in the Solar neighbourhood. GAIA will revolutionise our knowledge of the origin and evolution of our Milky Way Galaxy, and of the distribution of planetary systems around other stars.

This contribution is intended as a `rough guide' to the Hipparcos Catalogue for the non-expert user. Some general aspects of the use of astrometric data are discussed as well as Hipparcos-specific applications. We discuss when and at what level one may expect systematic errors to occur in the Hipparcos Catalogue. Next we discuss the question of the interpretation of the measured parallaxes in terms of distances and luminosities of stars. What are the biases one should be aware of and how can these be corrected? When using the astrometric data to study the statistics of stars one should take the full covariance matrix of the errors on the astrometric parameters into account. We explain how to do this and discuss the specific case of a moving cluster. Finally, we address the question of the correlation of astrometric parameters over a given region of the sky. At present the Hipparcos Catalogue contains no identified systematic errors.