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.

We study the radial migration of stars as a function of orbital action as well as the structural properties of a large suite of N-body simulations of isolated disc galaxies. Our goal is to establish a relationship between the radial migration efficiency of stars and their vertical action. We aim to describe how that relationship depends on the relative gravitational dominance between the disc and the dark matter halo. By changing the mass ratio of our disc and dark matter halo we find a relationship between disc dominance, number and strength of spiral arms, and the ensuing radial migration as a function of the vertical action. We conclude that the importance of migration at large vertical action depends on the strength of the spiral arms and therefore the dominance of the disc. Populations with more radial action undergo less radial migration, independently of disc dominance. Our results are important for the future of analytical modelling of radial migration in galaxies and furthers the understanding of radial migration which is a key component of the restructuring of galaxies, including the Milky Way.

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 a revolutionary space mission developed by ESA and is delivering 5 parameter astrometry, photometry and radial velocities over the whole sky with astrometric accuracies down to a few tens of micro-arcseconds. 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. This problem can be overcome by switching to the Near Infra-Red (NIR) but this is not possible with silicon CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires the detectors operate in Time Delayed Integration (TDI) mode or similar.

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 introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The One Thousand and One Magellanic Fields (1001MC) survey aims to measure the kinematics and elemental abundances of many different stellar populations that sample the history of formation and interaction of the Magellanic Clouds. The survey will collect spectra of about half a million stars with $G < 19.5$ magnitudes (Vega) distributed over an area of about 1000 square degrees and will provide an invaluable dataset for a wide range of scientific applications.

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.

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.

(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).

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

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

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.

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.

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.

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.

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.

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 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 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 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.

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

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.