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.

The light-trajectory in the gravitational field of N extended bodies in arbitrary motion is determined in the first post-Newtonian approximation. According to the theory of reference systems, the gravitational fields of these massive bodies are expressed in terms of their intrinsic multipoles, allowing for arbitrary shape and inner structure of these bodies. The results of this investigation aim towards a consistent general-relativistic theory of light propagation in the Solar system for high-precision astrometry at sub-micro-arcsecond level of accuracy.

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.

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.

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.

Two photon correlations are studied for a binary star system. It is investigated how the differential Hanbury Brown and Twiss (HBT) approach can be used in order to determine orbital parameters of a binary star.

High precision astrometry, space missions and certain tests of General Relativity, require the knowledge of the metric tensor of the solar system, or more generally, of a gravitational system of N extended bodies. Presently, the metric of arbitrarily shaped, rotating, oscillating and arbitrarily moving N bodies of finite extension is only known for the case of slowly moving bodies in the post-Newtonian approximation, while the post-Minkowskian metric for arbitrarily moving celestial objects is known only for pointlike bodies with mass-monopoles and spin-dipoles. As one more step towards the aim of a global metric for a system of N arbitrarily shaped and arbitrarily moving massive bodies in post-Minkowskian approximation, two central issues are on the scope of our investigation: (i) We first consider one extended body with full multipole structure in uniform motion in some suitably chosen global reference system. For this problem a co-moving inertial system of coordinates can be introduced where the metric, outside the body, admits an expansion in terms of Damour-Iyer moments. A Poincare transformation then yields the corresponding metric tensor in the global system in post-Minkowskian approximation. (ii) It will be argued why the global metric, exact to post-Minkowskian order, can be obtained by means of an instantaneous Poincare transformation for the case of pointlike mass-monopoles and spin-dipoles in arbitrary motion.

The light deflection of one component of a binary system due to the gravitational field of the other component is investigated. While this relativistic effect has not been observed thus far, the question arises that whether this effect becomes detectable in view of todays high-precision astrometry which soon will reach the microarcsecond level of accuracy. The effect is studied and its observability is investigated. It turns out, that in total there are about 1000 binaries having orbital parameters such that the light deflection amounts to be at least 1 microarcsecond. Two stringent criteria for the orbital parameters are presented, by means of which one can easily determine the maximal value of light deflection effect for a given binary system. It is found, that for relevant binaries their orbital parameters must take rather extreme values in order to have a light deflection of the order of a few microarcseconds. Only in a very few and rather extreme binary systems the light deflection effect might be detectable by todays astrometry, but their existence is highly improbable. Thus, the detection of this subtle effect of relativity still remains a challenge for future astrometric missions.

A generalized lens equation for weak gravitational fields in Schwarzschild metric and valid for finite distances of source and observer from the light deflecting body is suggested. The magnitude of neglected terms in the generalized lens equation is estimated to be smaller than or equal to 15 Pi/4 (m/d')^2, where m is the Schwarzschild radius of massive body and d' is Chandrasekhar's impact parameter. The main applications of this generalized lens equation are extreme astrometrical configurations, where 'Standard post-Newtonian approach' as well as 'Classical lens equation' cannot be applied. It is shown that in the appropriate limits the proposed lens equation yields the known post-Newtonian terms, 'enhanced' post-post-Newtonian terms and the Classical lens equation, thus provides a link between these both essential approaches for determining the light deflection.

Although the formulas for the light deflection due to quadrupole gravitational field of deflecting bodies are well known, the formulas are rather complicated, so that massive computations of quadrupole light deflection (e.g., in the framework of astrometric survey missions like Gaia) are time-consuming. Considering an observer situated within a few million kilometers from the Earth (clearly the most practical case), we derive the simplest possible form of the relevant formulas still having numerical accuracy of 1 micro-arcsecond. This form leads to simple upper estimates for the quadrupole light deflection in various cases allowing one to relate the magnitude of the actual quadrupole deflection with the corresponding monopole deflection due to the same body. These upper estimates can be used to decide if, for a given configuration, the actual quadrupole deflection should be computed for a given accuracy goal.

Efficient computation of the quadrupole light deflection for quasars/quasars and solar system objects within the framework of the baseline Gaia relativity model (GREM) is discussed. Two refinements have been achieved with the goal to improve the performance of the model: First, the quadrupole deflection formulas for both cases are simplified as much as possible considering the Gaia nominal orbit (only approximate minimal distances between Gaia and the giant planets were used here), physical parameters of the giant planets and the envisaged accuracy of 1 microarcsecond for individual systematic effects. The recommended formulas are given by Eq.(40) for stars/quasars and by Eq.(81) for solar system objects. Second, simple expressions for the upper estimate of the quadrupole light deflection have been found allowing, with a few additional arithmetical operations, to judge a priori if the quadrupole light deflection should be computed or not for a given source and for a given requested accuracy. The recommended criteria are given by Eq.(45) for stars/quasars and by Eq.(92) for solar system objects. Additionally, the quadrupole Shapiro effect for solar system objects is reconsidered. A strict upper bound of quadrupole Shapiro effect for solar system objects is given in Eq.(109).

This paper summarizes the analysis of the consequences of the violation of the Local Lorentz Invariance (LLI) on astrometric observations. We demonstrate that from the point of view of the LLI astrometric observations represent an experiment of Michelson-Morley type. The future high-accuracy astrometric projects (e.g., Gaia) will be used to test the LLI.

Numerical integration of the differential equations of light propagation in the Schwarzschild metric shows that in some situations relevant for practical observations the well-known post-Newtonian solution for light propagation has an error up to 16 microarcsecond. The aim of this work is to demonstrate this fact, identify the reason for this error and to derive an analytical formula accurate at the level of 1 microarcsecond as needed for high-accuracy astrometric projects (e.g., Gaia). An analytical post-post-Newtonian solution for the light propagation for both Cauchy and boundary problems is given for the Schwarzschild metric augmented by the PPN and post-linear parameters $\beta$, $\gamma$ and $\epsilon$. Using analytical upper estimates of each term we investigate which post-post-Newtonian terms may play a role for an observer in the solar system at the level of 1 microarcsecond and conclude that only one post-post-Newtonian term remains important for this numerical accuracy. In this way, an analytical solution for the boundary problem for light propagation is derived. That solution contains terms of both post-Newtonian and post-post-Newtonian order, but is valid for the given numerical level of 1 microarcsecond. The derived analytical solution has been verified using the results of a high-accuracy numerical integration of differential equations of light propagation and found to be correct at the level well below 1 microarcsecond for arbitrary observer situated within the solar system. Furthermore, the origin of the post-post-Newtonian terms relevant for the microarcsecond accuracy is elucidated. We demonstrate that these terms result from an inadequate choice of the impact parameter in the standard post-Newtonian formulas.

In [2] we have seen that among the post-post-Newtonian terms in the light deflection formulas there are "enhanced" ones that may become much larger than the other post-post-Newtonian terms. It is demonstrated here that these "enhanced" terms result from an inadequate choice of the impact parameter. Introducing another impact parameter, that can be considered as coordinate-independent, we demonstrate that all "enhanced" terms disappear from the formulas.

Numerical integration of the differential equations of light propagation in the Schwarzschild metric shows that in some extreme situations relevant for practical observations (e.g. for Gaia) the well-known standard post-Newtonian formula for the boundary problem has an error up to 16 \muas. The aim of this note is to identify the reason for this error and to derive an extended formula accurate at the level of 1 \muas as needed e.g. for Gaia. The analytical parametrized post-post-Newtonian solution for light propagation derived by \citet{report1} gives the solution for the boundary problem with all analytical terms of order $\OO4$ taken into account. Giving an analytical upper estimates of each term we investigate which post-post-Newtonian terms may play a role for an observer in the solar system at the level of 1 \muas. We conclude that only one post-post-Newtonian term remains important for this numerical accuracy and derive a simplified analytical solution for the boundary problem for light propagation containing all the terms that are indeed relevant at the level of 1 \muas. The derived analytical solution has been verified using the results of a high-accuracy numerical integration of differential equations of light propagation and found to be correct at the level well below 1 \muas for arbitrary observer situated within the solar system.

A rigorous analytical solution of light propagation in Schwarzschild metric in post-post Newtonian approximation has been presented in \cite{report1}. In \cite{report2} it has been asserted that the sum of all those terms which are of order ${{\cal O} (\frac{m^2}{d^2})}$ and ${{\cal O}(\frac{m^2}{d_\sigma^2})}$ is not greater than $15/4 \pi \frac{m^2}{d^2}}$ and $15/4 \pi \frac{m^2}{d_\sigma^2}}$, respectively. All these terms can be neglected on microarcsecond level of accuracy, leading to considerably simplified analytical transformations of light propagation. In this report, we give formal mathematical proofs for the inequalities used in the appendices of \cite{report2}.

An analytical solution for light propagation in the post-post-Newtonian approximation is given for the Schwarzschild metric in harmonic gauge augmented by PPN and post-linear parameters $\beta$, $\gamma$ and $\epsilon$. The solutions of both Cauchy and boundary problem are given. The Cauchy problem is posed using the initial position of the photon $\ve{x}_0 = \ve{x}(t_0)$ and its propagation direction \ve{\sigma} at minus infinity: $\ve{\sigma} = {1\over c} \lim\limits_{t \to -\infty}\dot{\ve{x}}(t)$. An analytical expression for the total light deflection is given. The solutions for $t - t_0$ and $\ve{\sigma}$ are given in terms of boundary conditions $\ve{x}_0 = \ve{x} (t_0)$ and $\ve{x} = \ve{x}(t)$.