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.

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 function--called a reception (or emission) time transfer function--has 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 quasi-conjunctions, where higher-order 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 on-going 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.

A new derivation of the propagation direction of light is given for a 3-parameter family of static, spherically symmetric space-times within the post-post-Newtonian framework. The emitter and the observer are both located at a finite distance. The case of a ray emitted at infinity is also treated.

A new method based on Synge's world function is developed for determining within the WKB approximation the gravitationally induced quantum phase shift of a particle propagating in a stationary spacetime. This method avoids any calculation of geodesics. A detailed treatment is given for relativistic particles within the weak field, linear approximation of any metric theory. The method is applied to the calculation of the oscillation terms governing the interference of neutrinos considered as a superposition of two eigenstates having different masses. It is shown that the neutrino oscillations are not sensitive to the gravitomagnetic components of the metric as long as the spin contributions can be ignored. Explicit calculations are performed when the source of the field is a spherical, homogeneous body. A comparison is made with previous results obtained in Schwarzschild spacetime.

Future space astrometry missions are planned to measure positions and/or parallaxes of celestial objects with an accuracy of the order of the microarcsecond. At such a level of accuracy, it will be indispensable to take into account the influence of the mass multipole structure of the giant planets on the bending of light rays. Within the parametrized post-Newtonian formalism, we present an algorithmic procedure enabling to determine explicitly this influence on a light ray connecting two points located at a finite distance. Then we specialize our formulae in the cases where 1) the light source is located at space infinity, 2) both the light source and the observer are located at space infinity. We examine in detail the cases where the unperturbed ray is in the equatorial plane or in a meridian plane.

Almost all of the studies devoted to relativistic astrometry are based on the integration of the null geodesic differential equations. However, the gravitational deflection of light rays can be calculated by a different method, based on the determination of the bifunction giving half the squared geodesic distance between two arbitrary points-events, the so-called Synge's world function. We give a brief review of the main results obtained by this method.

In suitably chosen domains of space-time, the world function may be a powerful tool for modelling the deflection of light and the time/frequency transfer. In this paper we work out a recursive procedure for expanding the world function into a perturbative series of ascending powers of the Newtonian gravitational constant G. We show rigorously that each perturbation term is given by a line integral taken along the unperturbed geodesic between two points. Once the world function is known, it becomes possible to determine the time transfer functions giving the propagation time of a photon between its emission and its reception. We establish that the direction of a light ray as measured in the 3-space relative to a given observer can be derived from these time transfer functions, even if the metric is not stationary. We show how to deduce these functions up to any given order in G from the perturbative expansion of the world function. To illustrate the method, we carry out the calculation of the world function and of the time transfer function outside a static, spherically symmetric body up to the order G^2, the metric containing three arbitrary parameters.

Within the weak-field, post-Newtonian approximation of the metric theories of gravity, we determine the one-way time transfer up to the order 1/c^4, the unperturbed term being of order 1/c, and the frequency shift up to the order 1/c^4. We adapt the method of the world-function developed by Synge to the Nordtvedt-Will PPN formalism. We get an integral expression for the world-function up to the order 1/c^3 and we apply this result to the field of an isolated, axisymmetric rotating body. We give a new procedure enabling to calculate the influence of the mass and spin multipole moments of the body on the time transfer and the frequency shift up to the order 1/c^4. We obtain explicit formulas for the contributions of the mass, of the quadrupole moment and of the intrinsic angular momentum. In the case where the only PPN parameters different from zero are beta and gamma, we deduce from these results the complete expression of the frequency shift up to the order 1/c^4. We briefly discuss the influence of the quadrupole moment and of the rotation of the Earth on the frequency shifts in the ACES mission.

A deep HST image of the MS 2137-2353 core has revealed detailed morphological structures in two arc systems, which are modelled and well reproduced after a complete analysis of the lensing properties of the dark matter component. Latter could have a simple elliptical mass distribution with ellipticity and angular orientation similar to those of the visible and X-ray light, which suggests that the MS 2137-2353 is a relaxed cluster at z=0.313. The predicted density profile ($\rho$ $\sim$ $r^{-1.56\pm0.1}$ with $r_{c}$ $\le$ 22.5$h_{50}^{-1}$ kpc) within 150$h_{50}^{-1}$ kpc implies increasing M/L ratio with the radius, and could be in agreement with predictions from standard CDM simulations. At least two faint sources (unlensed magnitude, R=23.9 and 26, respectively) are aligned with the cluster core and are responsible of the arc systems. They have been reconstructed with details as small as 0".02 (or 160$h_{50}^{-1}$ pc in the source assumed at z= 1), one could be a nearly edge-on barred spiral galaxy, and the other has a more complex morphology, which could be related to a close interacting pair and/or to dust. They show strong signs of star formation indicated by compact HII regions well off their center. The observation of giant luminous arcs by HST could even resolve the size of giant HII regions at z$\sim$ 1.