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.

We study the formation of voids in a modified gravity model in which gravity is generically stronger or weaker on large scales. We show that void abundances provide complementary information to halo abundances: if normalized to the CMB, models with weaker large-scale gravity have smaller large scale power, fewer massive halos and fewer large voids, although the scalings are not completely degenerate with $\sigma_8$. Our results suggest that, in addition to their abundances, halo and void density profiles may also provide interesting constraints on such models: stronger large scale gravity produces more concentrated halos, and thinner void walls. This potentially affects the scaling relations commonly assumed to translate cluster observables to halo masses, potentially making these too, useful probes of gravity.

Modifications to the gravitational potential affect the nonlinear gravitational evolution of large scale structures in the Universe. To illustrate some generic features of such changes, we study the evolution of spherically symmetric perturbations when the modification is of Yukawa type; this is non-trivial, because we should not and do not assume that Birkhoff's theorem applies. We then show how to estimate the abundance of virialized objects in such models. Comparison with numerical simulations shows reasonable agreement: When normalized to have the same fluctuations at early times, weaker large scale gravity produces fewer massive halos. However, the opposite can be true for models that are normalized to have the same linear theory power spectrum today, so the abundance of rich clusters potentially places interesting constraints on such models. Our analysis also indicates that the formation histories and abundances of sufficiently low mass objects are unchanged from standard gravity. This explains why simulations have found that the nonlinear power-spectrum at large k is unaffected by such modifications to the gravitational potential. In addition, the most massive objects in CMB-normalized models with weaker gravity are expected to be similar to the high-redshift progenitors of the most massive objects in models with stronger gravity. Thus, the difference between the cluster and field galaxy populations is expected to be larger in models with stronger large-scale gravity.

In studies of the environmental dependence of structure formation, the large scale environment is often thought of as providing an effective background cosmology: e.g. the formation of structure in voids is expected to be just like that in a less dense universe with appropriately modified Hubble and cosmological constants. However, in the excursion set description of structure formation which is commonly used to model this effect, no explicit mention is made of the effective cosmology. Rather, this approach uses the spherical evolution model to compute an effective linear theory growth factor, which is then used to predict the growth and evolution of nonlinear structures. We show that these approaches are, in fact, equivalent: a consequence of Birkhoff's theorem. We speculate that this equivalence will not survive in models where the gravitational force law is modified from an inverse square, potentially making the environmental dependence of clustering a good test of such models.

Halo model interpretations of the luminosity dependence of galaxy clustering assume that there is a central galaxy in every sufficiently massive halo, and that this central galaxy is very different from all the others in the halo. The halo model decomposition makes the remarkable prediction that the mean luminosity of the non-central galaxies in a halo should be almost independent of halo mass: the predicted increase is about 20% while the halo mass increases by a factor of more than 20. In contrast, the luminosity of the central object is predicted to increase approximately linearly with halo mass at low to intermediate masses, and logarithmically at high masses. We show that this weak, almost non-existent mass-dependence of the satellites is in excellent agreement with the satellite population in group catalogs constructed by two different collaborations. This is remarkable, because the halo model prediction was made without ever identifying groups and clusters. The halo model also predicts that the number of satellites in a halo is drawn from a Poisson distribution with mean which depends on halo mass. This, combined with the weak dependence of satellite luminosity on halo mass, suggests that the Scott effect, such that the luminosities of very bright galaxies are merely the statistically extreme values of a general luminosity distribution, may better apply to the most luminous satellite galaxy in a halo than to BCGs. If galaxies are identified with halo substructure at the present time, then central galaxies should be about 4 times more massive than satellite galaxies of the same luminosity, whereas the differences between the stellar M/L ratios should be smaller. Therefore, a comparison of the weak lensing signal from central and satellite galaxies should provide useful constraints. [abridged]