
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.

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.

The paper deals with a special kind of problems that appear in solutions of
Einstein's field equations for extended bodies: many structuredependent terms
appear in intermediate calculations that cancel exactly in virtue of the local
equations of motion or can be eliminated by appropriate gauge transformations.
For a single body at rest these problems are well understood for both the
postNewtonian and the postMinkowskian cases. However, the situation is still
unclear for approximations of higher orders. This paper discusses this problem
for a "body" of spherical symmetry to postlinear order. We explicitly
demonstrate how the usual Schwarzschild field can be derived directly from the
field equations in the postlinear approximation in the harmonic gauge and for
an arbitrary spherically symmetric matter distribution. Both external and
internal solutions are considered. The case of static incompressible fluid is
then compared to the wellknown results from the literature. The results of
this paper can be applied to generalize the wellknown postNewtonian and
postMinkowskian multipole expansions of the metric in the postlinear
approximation.

Comparison of stellar catalogues with position and proper motion components
using a decomposition on a set of orthogonal vector spherical harmonics. We
show the theoretical and practical advantages of this technique as a result of
invariance properties and the independence of the decomposition from a prior
model. We describe the mathematical principles used to perform the spectral
decomposition, evaluate the level of significance of the multipolar components
and examine the transformation properties under space rotation. The principles
are illustrated with a characterisation of the systematic effects in the FK5
catalogue compared to Hipparcos and with an application to the extraction of
the rotation and dipole acceleration in the astrometric solution of QSOs
expected from Gaia.

ASTROD I is a planned interplanetary space mission with multiple goals. The
primary aims are: to test General Relativity with an improvement in sensitivity
of over 3 orders of magnitude, improving our understanding of gravity and
aiding the development of a new quantum gravity theory; to measure key solar
system parameters with increased accuracy, advancing solar physics and our
knowledge of the solar system and to measure the time rate of change of the
gravitational constant with an order of magnitude improvement and the anomalous
Pioneer acceleration, thereby probing dark matter and dark energy
gravitationally. It is an international project, with major contributions from
Europe and China and is envisaged as the first in a series of ASTROD missions.
ASTROD I will consist of one spacecraft carrying a telescope, four lasers, two
event timers and a clock. Twoway, twowavelength laser pulse ranging will be
used between the spacecraft in a solar orbit and deep space laser stations on
Earth, to achieve the ASTROD I goals. A second mission, ASTROD II is envisaged
as a threespacecraft mission which would test General Relativity to one part
per billion, enable detection of solar gmodes, measure the solar
LenseThirring effect to 10 parts per million, and probe gravitational waves at
frequencies below the LISA bandwidth. In the third phase (ASTROD III or
SuperASTROD), larger orbits could be implemented to map the outer solar system
and to probe primordial gravitationalwaves at frequencies below the ASTROD II
bandwidth.

The first order post Newtonian scheme in multiple systems presented by
DamourSoffelXu is extended to the second order one for light propagation
without changing the advantage of the scheme on the linear partial differential
equations of potential and vector potential. The spatial components of the
metric tensor are extended to the second order level both in the global
coordinates ($q_{ij}/c^4$ term) and in a local coordinates ($Q_{ab}/c^4$ term).
The equations of $q_{ij}$ (or $Q_{ab}$) are deduced from Einstein field
equations. The linear relationship between $q_{ij}$ and $Q_{ab}$ are presented
also. The 2PN equations of light ray based on the extended scheme are deduced
by means of the iterative method. We also use parametrized second post
Newtonian metric tensor to substitute into the null geodetic equations to
obtain the parametrized second order equations of light ray which might be
useful in the observation and measurement in the future space missions.

In this paper, the dynamical equations and junction conditions at the
interface between adjacent layers of different elastic properties for an
elastic deformable astronomical body in the first postNewtonian approximation
of Einstein theory of gravity are discussed in both rotating Cartesian
coordinates and rotating spherical coordinates. The unperturbed rotating body
(the ground state) is described as uniformly rotating, stationary and
axisymmetric configuration in an asymptotically flat spacetime manifold.
Deviations from the equilibrium configuration are described by means of a
displacement field. In terms of the formalism of relativistic celestial
mechanics developed by Damour, Soffel and Xu, and the framework established by
Carter and Quintana the post Newtonian equations of the displacement field and
the symmetric tracefree shear tensor are obtained. Corresponding
postNewtonian junction conditions at interfaces also the outer surface
boundary conditions are presented. The PN junction condition is an extension of
Wahr's one which is a Newtonian junction conditions without rotating.