In the potential theory, isochrony was introduced by Michel H\'enon in 1959 to characterize astrophysical observations of some globular clusters. Today, Michel Henon's isochrone potential is mainly used for his integrable property in numerical simulations, but is generally not really known. In a recent paper [29], we have presented new fundamental and theoretical results about isochrony that have particular importance in self-gravitating dynamics and which are detailed in this paper. In particular, new characterization of the isochrone state has been proposed which are investigated in order to analyze the product of the fast relaxation of a self-gravitating system. The general paradigm consists in considering that this product is a lowered isothermal sphere (King Model). By a detailed numerical study we show that this paradigm fails when the isochrone model succeeds in reproducing the quasi-equilibrium state obtained just after fast relaxation.

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 Gaia Data Release 2 provides a parallax of 0.734+/-0.073 mas for SDSS J102915+172927, currently the most metal-poor known object. This parallax implies that it is dwarf star, ruling out the scenario that it is a subgiant. The subgiant scenario had as a corollary that the star had been formed in a medium highly enriched in C, thus making line cooling efficient during the collapse, that was also highly enriched in Fe by Type Ia SNe. This scenario can also now be ruled out for this star, reinforcing the need of dust cooling and fragmentation to explain its formation.

For Gaia DR2 (GDR2), 280 million spectra, collected by the RVS instrument on-board Gaia, were processed and median radial velocities were derived for 9.8 million sources brighter than Grvs = 12 mag. This paper describes the validation and properties of the median radial velocities published in GDR2. Quality tests and filters are applied to select, from the 9.8 million radial velocities, those with the quality to be published in GDR2. The accuracy of the selected sample is assessed with respect to ground-based catalogues. Its precision is estimated using both ground-based catalogues and the distribution of the Gaia radial velocity uncertainties. GDR2 contains median radial velocities for 7 224 631 stars, with Teff in the range [3550, 6900] K, which passed succesfully the quality tests. The published median radial velocities provide a full sky-coverage and have a completness with respect to the astrometric data of 77.2\% (for $G \leq 12.5$ mag). The median radial velocity residuals with respect to the ground-based surveys vary from one catalogue to another, but do not exceed a few 100s m/s. In addition, the Gaia radial velocities show a positive trend as a function of magnitude, which starts around Grvs $\sim 9$ mag and reaches about $+500$ m/s at Grvs $= 11.75$ mag. The overall precision, estimated from the median of the Gaia radial velocity uncertainties, is 1.05 km/s. The radial velocity precision is function of many parameters, in particular the magnitude and effective temperature. For bright stars, Grvs in [4, 8] mag, the precision is in the range 200-350 m/s, which is about 3 to 5 times more precise than the pre-launch specification of 1 km/s. At the faint end, Grvs = 11.75 mag, the precisions for Teff = 5000 K and 6500 K are respectively 1.4 km/s and 3.7 km/s.

The Gaia Data Release 2 contains the 1st release of radial velocities complementing the kinematic data of a sample of about 7 million relatively bright, late-type stars. Aims: This paper provides a detailed description of the Gaia spectroscopic data processing pipeline, and of the approach adopted to derive the radial velocities presented in DR2. Methods: The pipeline must perform four main tasks: (i) clean and reduce the spectra observed with the Radial Velocity Spectrometer (RVS); (ii) calibrate the RVS instrument, including wavelength, straylight, line-spread function, bias non-uniformity, and photometric zeropoint; (iii) extract the radial velocities; and (iv) verify the accuracy and precision of the results. The radial velocity of a star is obtained through a fit of the RVS spectrum relative to an appropriate synthetic template spectrum. An additional task of the spectroscopic pipeline was to provide 1st-order estimates of the stellar atmospheric parameters required to select such template spectra. We describe the pipeline features and present the detailed calibration algorithms and software solutions we used to produce the radial velocities published in DR2. Results: The spectroscopic processing pipeline produced median radial velocities for Gaia stars with narrow-band near-IR magnitude Grvs < 12 (i.e. brighter than V~13). Stars identified as double-lined spectroscopic binaries were removed from the pipeline, while variable stars, single-lined, and non-detected double-lined spectroscopic binaries were treated as single stars. The scatter in radial velocity among different observations of a same star, also published in DR2, provides information about radial velocity variability. For the hottest (Teff > 7000 K) and coolest (Teff < 3500 K) stars, the accuracy and precision of the stellar parameter estimates are not sufficient to allow selection of appropriate templates. [Abridged]

Aims. The Radial Velocity Spectrometer (RVS) on board the ESA satellite mission Gaia has no calibration device. Therefore, the radial velocity zero point needs to be calibrated with stars that are proved to be stable at a level of 300 m/s during the Gaia observations. Methods. We compiled a dataset of ~71000 radial velocity measurements from five high-resolution spectrographs. A catalogue of 4813 stars was built by combining these individual measurements. The zero point was established using asteroids. Results. The resulting catalogue has seven observations per star on average on a typical time baseline of six years, with a median standard deviation of 15 m/s. A subset of the most stable stars fulfilling the RVS requirements was used to establish the radial velocity zero point provided in Gaia Data Release 2. The stars that were not used for calibration are used to validate the RVS data.

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.

We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized H\'enon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual and quantitative comparisons are performed. In particular we introduce new statistics, correlators and entropic estimators, based on the likelihood of whether N-body simulations actually trace randomly the Vlasov phase-space density. When taking into account the limits of both the N-body and the Vlasov codes, namely collective effects due to the particle shot noise in the first case and diffusion and possible nonlinear instabilities due to finite resolution of the phase-space grid in the second case, we find a spectacular agreement between both methods, even in regions of phase-space where nontrivial physical instabilities develop. However, in the colder case, R=0.1, it was not possible to prove actual numerical convergence of the N-body results after a number of dynamical times, even with N=10$^8$ particles.