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.

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.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

One of the most accurate models currently used to represent the gravity field of irregular bodies is the polyhedral approach. In this model, the mass of the body is assumed to be homogeneous, which may not be true for a real object. The main goal of the present paper is to study the dynamical effects induced by three different internal structures (uniform, three- and four-layers) of asteroid (21) Lutetia, an object that recent results from space probe suggest being at least partially differentiated. The Mascon gravity approach used in the present work, consists of dividing each tetrahedron into eight parts to calculate the gravitational field around the asteroid. The zero-velocity curves show that the greatest displacement of the equilibrium points occurs in the position of the E4 point for the four-layers structure and the smallest one occurs in the position of the E3 point for the three-layers structure. Moreover, stability against impact shows that the planar limit gets slightly closer to the body with the four-layered structure. We then investigated the stability of orbital motion in the equatorial plane of (21) Lutetia and propose numerical stability criteria to map the region of stable motions. Layered structures could stabilize orbits that were unstable in the homogeneous model.

The aim of this study is first to determine the gravity field of the comet 67P Churyumov-Gerasimenko and second to derive the solar component of the precession rate and nutation coefficients of the spin axis of the comet nucleus, i.e. without the direct, usually larger, effect of outgassing. The gravity field, and related moments of inertia, are obtained from two polyhedra, that are provided by the OSIRIS and NAV-CAM experiments on Rosetta, and are based on the assumption of uniform density for the comet nucleus. We also calculate the forced precession rate as well as the nutation coefficients on the basis of Kinoshita's theory of rotation of the rigid Earth and adapted it to be able to indirectly include the effect of outgassing on the rotational parameters. The 2nd degree denormalized Stokes coefficients of comet 67P C-G turn out to be (bracketed numbers refer to second shape model) C20 = -6.74 [-7.93] * 10^-2, C22 = 2.60 [2.71] * 10^-2 consistent with normalized principal moments of inertia A/MR^2 = 0.13 [0.11], B/MR^2 = 0.23 [0.22], with polar moment c = C/MR^2 = 0.25, depending on the choice of the polyhedron model. The obliquity between the rotation axis and the mean orbit normal is 52 degree, and the precession rate only due to solar torques becomes 20-30 "/y. Oscillations in longitude caused by the gravitational pull of the Sun turn out to be of the order of 1', oscillations in obliquity can be estimated to be of the order of 0.5'.

The goal of this presentation is to report the latest progress in creation of the next generation of VLBI-based International Celestial Reference Frame, ICRF3. Two main directions of ICRF3 development are improvement of the S/X-band frame and extension of the ICRF to higher frequencies. Another important task of this work is the preparation for comparison of ICRF3 with the new generation optical frame GCRF expected by the end of the decade as a result of the Gaia mission.

% context heading (optional) {Thanks to the Venus Express Mission, new data on the properties of Venus could be obtained in particular concerning its rotation.} % aims heading (mandatory) {In view of these upcoming results, the purpose of this paper is to determine and compare the major physical processes influencing the rotation of Venus, and more particularly the angular rotation rate.} % methods heading (mandatory) {Applying models already used for the Earth, the effect of the triaxiality of a rigid Venus on its period of rotation are computed. Then the variations of Venus rotation caused by the elasticity, the atmosphere and the core of the planet are evaluated.} % results heading (mandatory) {Although the largest irregularities of the rotation rate of the Earth at short time scales are caused by its atmosphere and elastic deformations, we show that the Venus ones are dominated by the tidal torque exerted by the Sun on its solid body. Indeed, as Venus has a slow rotation, these effects have a large amplitude of 2 minutes of time (mn). These variations of the rotation rate are larger than the one induced by atmospheric wind variations that can reach 25-50 seconds of time (s), depending on the simulation used. The variations due to the core effects which vary with its size between 3 and 20s are smaller. Compared to these effects, the influence of the elastic deformation cause by the zonal tidal potential is negligible.} % conclusions heading (optional), leave it empty if necessary {As the variations of the rotation of Venus reported here are of the order 3mn peak to peak, they should influence past, present and future observations providing further constraints on the planet internal structure and atmosphere.}

Although the rotation of some Saturn's satellites in spin-orbit has already been studied by several authors, this is not the case of the rotation of Phoebe, which has the particularity of being non resonant. The purpose of the paper is to determine for the first time and with precision its precession-nutation motion. We adopt an Hamiltonian formalism of the motion of rotation of rigid celestial body set up by Kinoshita (1977) based on Andoyer variables and canonical equations. First we calculate Phoebe's obliquity at J2000,0 from available astronomical data as well as the gravitational perturbation due to Saturn on Phoebe rotational motion. Then we carry out a numerical integration and we compare our results for the precession rate and the nutation coefficients with pure analytical model. Our results for Phoebe obliquity (23{\deg}95) and Phoebe precession rate (5580".65/cy) are very close to the respective values for the Earth. Moreover the amplitudes of the nutations (26" peak to peak for the nutaton in longitude and 8" for the nutation in obliquity) are of the same order as the respective amplitudes for the Earth. We give complete tables of nutation, obtained from a FFT analysis starting from the numerical signals. We show that a pure analytical model of the nutation is not accurate due to the fact that Phoebe orbital elements e, M and Ls are far from having a simple linear behaviour. The precession and nutation of Phoebe have been calculated for the first time in this paper. We should keep on the study in the future by studying the additional gravitational effects of the Sun, of the large satellites as Titan, as well as Saturn dynamical ellipticity.

% context :The precise and accurate modelling of a terrestrial planet like Venus is an exciting and challenging topic, all the more interesting since it can be compared with that of the Earth for which such a modelling has already been achieved at the milliarcsecond level % aims: We want to complete a previous study (Cottereau and Souchay, 2009), by determining at the milliarcsecond level the polhody, i.e. the torque-free motion of the axis of angular momentum of a rigid Venus in a body-fixed frame, as well as the nutation of its third axis of figure in space, which is fundamental from an observational point of view. results :In a first part we have computed the polhody, i.e. the respective free rotational motion of the axis of angular momentum of Venus with respect to a body-fixed frame. We have shown that this motion is highly elliptical, with a very long period of 525 cy to be compared with 430 d for the Earth. This is due to the very small dynamical flattening of Venus in comparison with our planet. In a second part we have computed precisely the Oppolzer terms which allow to represent the motion in space of the third Venus figure axis with respect to Venus angular momentum axis, under the influence of the solar gravitational torque. We have determined the corresponding tables of coefficients of nutation of the third figure axis both in longitude and in obliquity due to the Sun, which are of the same order of amplitude as for the Earth. We have shown that the coefficients of nutation for the third figure axis are significantly different from those of the angular momentum axis on the contrary of the Earth. Our analytical results have been validated by a numerical integration which revealed the indirect planetary effects.

Context: With the increasing knowledge of the terrestrial planets due to recent space probes it is possible to model their rotation with increasing accuracy. Despite that fact, an accurate determination of Venus precession and nutation is lacking. Aims : Although Venus rotation has been studied in several aspects, a full and precise analytical model of its precession-nutation motion remains to be constructed. We propose to determine this motion with up-to-date physical parameters of the planet Methods: We adopt a theoritical framework already used for a precise precession-nutation model of the Earth, based on a Hamiltonian formulation, canonical equations and an accurate development of the perturbing function due to the Sun. Results: After integrating the disturbing function and applying the canonical equations, we can evaluate the precession constant $\dot{\Psi}$ and the coefficients of nutation, both in longitude and in obliquity. We get $\dot{\Psi}=4474".35/Jcy \pm 66.5 $, corresponding to a precession period of $28965.10 \pm 437$ years. This result, based on recent estimations of the Venus moment of inertia is significantly different from previous estimations. The largest nutation coefficient in longitude with an argument $2L_{S}$ (where $L_{S}$ is the longitude of the Sun) has a 2"19 amplitude and a 112.35 d period. We show that the coefficients of nutation of Venus due to its triaxiality are of the same order of amplitude as these values due to its dynamical flattening, unlike of the Earth, for which they are negligible.

This contribution is a first step aiming to address a general question: what can be concluded on impact craters which exist on various planetary system objects, by combining astronomical data and known theoretical results from solid state physics. Assuming that the material of the target body is of crystaline structure,it is shown that a simple calculation gives the possibility of estimating the speed of the impactor responsible for the creation of a crater.A test value,calculated using observed data on the composition of some asteroids,gives a value of the speed in good agreement with results of celestial mechanics.

The large number and all-sky distribution of quasars from different surveys, along with their presence in large, deep astrometric catalogs,enables the building of an optical materialization of the ICRS following its defining principles. Namely: that it is kinematically non-rotating with respect to the ensemble of distant extragalactic objects; aligned with the mean equator and dynamical equinox of J2000; and realized by a list of adopted coordinates of extragalatic sources. Starting from the updated and presumably complete LQAC list of QSOs, the initial optical positions of those quasars are found in the USNO B1.0 and GSC2.3 catalogs, and from the SDSS DR5. The initial positions are next placed onto UCAC2-based reference frames, following by an alignment with the ICRF, to which were added the most precise sources from the VLBA calibrator list and the VLA calibrator list - when reliable optical counterparts exist. Finally, the LQRF axes are inspected through spherical harmonics, contemplating to define right ascension, declination and magnitude terms. The LQRF contains J2000 referred equatorial coordinates for 100,165 quasars, well represented across the sky, from -83.5 to +88.5 degrees in declination, and with 10 arcmin being the average distance between adjacent elements. The global alignment with the ICRF is 1.5 mas, and the individual position accuracies are represented by a Poisson distribution that peaks at 139 mas in right ascension and 130 mas in declination. It is complemented by redshift and photometry information from the LQAC. The LQRF is designed to be an astrometric frame, but it is also the basis for the GAIA mission initial quasars' list, and can be used as a test bench for quasars' space distribution and luminosity function studies.

The International Celestial Reference Frame (ICRF, Ma et al. 1998) is currently the best realization of a quasi-inertial reference system. It is based on more than 10 years of cumulated geodetic and astrometric VLBI observations of compact extragalactic objects at centimetric wavelengths. In the perspective of the realization of an accurate optical counterpart of the ICRF using future space astrometry missions like GAIA or SIM, this paper investigates the consistency of celestial reference frames realized through the same subset of compact extragalactic radio sources at optical wavelengths. Celestial reference frames realized in radio wavelengths with the VLBA Calibrator Survey (VCS) data and in optical wavelengths with the Sloan Digital Sky Survey (SDSS) data (DR3 quasar catalogue and DR5) are compared in terms of radio-optical distances between the common sources, global rotation of the axes and offset of the equator.