
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
preprocessed 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 fiveparameter 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 nonrotating 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.

We highlight the power of the Gaia DR2 in studying many fine structures of
the HertzsprungRussell 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, postAGB stars, planetary nebulae, and white dwarfs are also
analysed, as well as lowmass 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
HRDbased studies.

The Einstein Equivalence Principle (EEP) is one of the foundations of the
theory of General Relativity and several alternative theories of gravitation
predict violations of the EEP. Experimental constraints on this fundamental
principle of nature are therefore of paramount importance. The EEP can be split
in three subprinciples: the Universality of Free Fall (UFF), the Local Lorentz
Invariance (LLI) and the Local Position Invariance (LPI). In this paper we
propose to use stable clocks in eccentric orbits to perform a test of the
gravitational redshift, a consequence of the LPI. The best test to date was
performed with the Gravity Probe A (GPA) experiment in 1976 with an
uncertainty of $1.4\times10^{4}$. Our proposal considers the opportunity of
using Galileo satellites 5 and 6 to improve on the GPA test uncertainty. We
show that considering realistic noise and systematic effects, and thanks to a
highly eccentric orbit, it is possible to improve on the GPA limit to an
uncertainty around $(34)\times 10^{5}$ after one year of integration of
Galileo 5 and 6 data.

Determining range, Doppler and astrometric observables is of crucial interest
for modelling and analyzing space observations. We recall how these observables
can be computed when the travel time of a light ray is known as a function of
the positions of the emitter and the receiver for a given instant of reception
(or emission). For a long time, such a functioncalled a reception (or
emission) time transfer functionhas been almost exclusively calculated by
integrating the null geodesic equations describing the light rays. However,
other methods avoiding such an integration have been considerably developped in
the last twelve years. We give a survey of the analytical results obtained with
these new methods up to the third order in the gravitational constant $G$ for a
mass monopole. We briefly discuss the case of quasiconjunctions, where
higherorder enhanced terms must be taken into account for correctly
calculating the effects. We summarize the results obtained at the first order
in $G$ when the multipole structure and the motion of an axisymmetric body is
taken into account. We present some applications to ongoing or future missions
like Gaia and Juno. We give a short review of the recent works devoted to the
numerical estimates of the time transfer functions and their derivatives.

Given the extreme accuracy of modern space science, a precise relativistic
modeling of observations is required. We use the Time Transfer Functions
formalism to study light propagation in the field of uniformly moving
axisymmetric bodies, which extends the field of application of previous works.
We first present a spacetime metric adapted to describe the geometry of an
ensemble of uniformly moving bodies. Then, we show that the expression of the
Time Transfer Functions in the field of a uniformly moving body can be easily
derived from its wellknown expression in a stationary field by using a change
of variables. We also give a general expression of the Time Transfer Function
in the case of an ensemble of arbitrarily moving point masses. This result is
given in the form of an integral easily computable numerically. We also provide
the derivatives of the Time Transfer Function in this case, which are mandatory
to compute Doppler and astrometric observables. We particularize our results in
the case of moving axisymmetric bodies. Finally, we apply our results to study
the different relativistic contributions to the range and Doppler tracking for
the JUNO mission in the Jovian system.

Given the extreme accuracy of modern space science, a precise relativistic
modeling of observations is required. In particular, it is important to
describe properly light propagation through the Solar System. For two decades,
several modeling efforts based on the solution of the null geodesic equations
have been proposed but they are mainly valid only for the first order
PostNewtonian approximation. However, with the increasing precision of ongoing
space missions as Gaia, GAME, BepiColombo, JUNO or JUICE, we know that some
corrections up to the second order have to be taken into account for future
experiments. We present a procedure to compute the relativistic coordinate time
delay, Doppler and astrometric observables avoiding the integration of the null
geodesic equation. This is possible using the Time Transfer Function formalism,
a powerful tool providing key quantities such as the time of flight of a light
signal between two pointevents and the tangent vector to its nullgeodesic.
Indeed we show how to compute the Time Transfer Functions and their derivatives
(and thus range, Doppler and astrometric observables) up to the second
postMinkowskian order. We express these quantities as quadratures of some
functions that depend only on the metric and its derivatives evaluated along a
Minkowskian straight line. This method is particularly well adapted for
numerical estimations. As an illustration, we provide explicit expressions in
static and spherically symmetric spacetime up to second postMinkowskian
order. Then we give the order of magnitude of these corrections for the
range/Doppler on the BepiColombo mission and for astrometry in a GAMElike
observation.

A lot of fundamental tests of gravitational theories rely on highly precise
measurements of the travel time and/or the frequency shift of electromagnetic
signals propagating through the gravitational field of the Solar System. In
practically all of the previous studies, the explicit expressions of such
travel times and frequency shifts as predicted by various metric theories of
gravity are derived from an integration of the null geodesic differential
equations. However, the solution of the geodesic equations requires heavy
calculations when one has to take into account the presence of mass multipoles
in the gravitational field or the tidal effects due to the planetary motions,
and the calculations become quite complicated in the postpostMinkowskian
approximation. This difficult task can be avoided using the time transfer
function's formalism. We present here our last advances in the formulation of
the oneway frequency shift using this formalism up to the
postpostMinkowskian approximation.

In this paper we review the current predictions of numerical simulations for
the origin and observability of the warm hot intergalactic medium (WHIM), the
diffuse gas that contains up to 50 per cent of the baryons at z~0. During
structure formation, gravitational accretion shocks emerging from collapsing
regions gradually heat the intergalactic medium (IGM) to temperatures in the
range T~10^510^7 K. The WHIM is predicted to radiate most of its energy in the
ultraviolet (UV) and Xray bands and to contribute a significant fraction of
the soft Xray background emission. While O VI and C IV absorption systems
arising in the cooler fraction of the WHIM with T~10^510^5.5 K are seen in
FUSE and HST observations, models agree that current Xray telescopes such as
Chandra and XMMNewton do not have enough sensitivity to detect the hotter
WHIM. However, future missions such as ConstellationX and XEUS might be able
to detect both emission lines and absorption systems from highly ionised atoms
such as O VII, O VIII and Fe XVII.

In this paper we consider Hopf's Lemma and the Strong Maximum Principle for
supersolutions to a class of non elliptic equations. In particular we prove a
sufficient condition for the validity of Hopf's Lemma and of the Strong Maximum
Principle and we give a condition which is at once necessary for the validity
of Hopf's Lemma and sufficient for the validity of the Strong Maximum
Principle.