
(abridged) We develop Bayesian methods and detection criteria for orbital
fitting, and revise the detectability of exoplanets in light of the inflight
properties of Gaia. Limiting ourselves to oneplanet 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 WatanabeAkaike
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 Jupitermass and Neptunemass planet,
respectively, assuming a 10 yr mission, a 4 au semimajor 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
semimajor 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 Earthlike
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 microarcsecond level. This is about 1000
times better than ESA/Gaia's accuracy for the brightest objects and represents
a factor 1030 improvement for the faintest stars (depending on the exact
observational program). In the version submitted to ESA, we proposed an optical
(3501000nm) onaxis 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.

Determination of absolute parallaxes by means of a scanning astrometric
satellite such as Hipparcos or Gaia relies on the shortterm stability of the
socalled 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.

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/Tycho2 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 nonrotating 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 infrared 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.

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
fiveparameter solution, and even after the full mission many stars will remain
underobserved, 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 underobserved 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
wellobserved stars it should not be used where a good fiveparameter
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 Tycho2
Catalogue as additional information for a joint solution with early Gaia data.
We call this the TychoGaia 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 Tycho2
stars, with submilliarcsecond 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
goodnessoffit statistic which is sensitive to deviations from uniform space
motion, caused for example by binaries with periods of 1050 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. NanoJASMINE, 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 brightstar measurements of Hipparcos,
the updated brightstar observations of NanoJASMINE, 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 (19891993) resulted in the first spacebased stellar
catalogue including measurements of positions, parallaxes and annual proper
motions accurate to about one milliarcsecond. More space astrometry missions
will follow in the near future. The ultrasmall Japanese mission NanoJASMINE
(launch in late 2013) will determine positions and annual proper motions with
some milliarcsecond accuracy. In mid 2013 the nextgeneration ESA mission Gaia
will deliver some tens of microarcsecond 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 NanoJASMINE data in a joint solution. We show a
significant improvement over the conventional catalogue combination.

(abridged) A complete census of planetary systems around a volumelimited
sample of solartype stars (FGK dwarfs) in the Solar neighborhood with uniform
sensitivity down to Earthmass 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 spaceborne
extremelyhighprecision 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 Earthmass 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 offaxis parabola singlemirror 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 Ktype 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 (lowmass
stars, young stars, etc.) discovered by Gaia, groundbased highprecision
radialvelocity 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 onboard 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
periodluminosity 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 zeropoint 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.