• 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.
• ### Cosmic-rays, gas, and dust in nearby anti-centre clouds : III -- Dust extinction, emission, and grain properties(1803.05686)

March 15, 2018 astro-ph.GA, astro-ph.HE
We have explored the capabilities of dust extinction and $\gamma$ rays to probe the properties of the interstellar medium in the nearby anti-centre region. We have jointly modelled the $\gamma$-ray intensity and the stellar reddening, E(B-V) as a combination of H$_{\rm I}$-bright, CO-bright, and ionised gas components. The complementary information from dust reddening and $\gamma$ rays is used to reveal the dark gas not seen, or poorly traced, by H$_{\rm I}$, free-free, and $^{12}$CO emissions. We compare the total gas column densities, $N_{\rm{H}}$, derived from the $\gamma$ rays and stellar reddening with those inferred from a similar analysis (Remy et al. 2017) of $\gamma$ rays and of the optical depth of the thermal dust emission, $\tau_{353}$, at 353 GHz. We can therefore compare environmental variations in specific dust reddening, E(B-V)/$N_{\rm H}$, and in dust emission opacity (dust optical depth per gas nucleon), $\tau_{353}/N_{\rm{H}}$. Over the whole anti-centre region, we find an average E(B-V)/$N_{\rm H}$ ratio of $(2.02\pm0.48)\times$ $10^{-22}$~mag~cm$^2$, with maximum local variations of about $\pm30\%$ at variance with the two to six fold coincident increase seen in emission opacity as the gas column density increases. In the diffuse medium, the small variations in specific reddening, E(B-V)/$N_{\rm H}$ implies a rather uniform dust-to-gas mass ratio in the diffuse parts of the anti-centre clouds. The small amplitude of the E(B-V)/$N_{\rm H}$ variations with increasing $N_{\rm{H}}$ column density confirms that the large opacity $\tau_{353}/N_{\rm{H}}$ rise seen toward dense CO clouds is primarily due to changes in dust emissivity. The environmental changes are qualitatively compatible with model predictions based on mantle accretion on the grains and the formation of grain aggregates.
• ### Cosmic-rays, gas, and dust in nearby anticentre clouds : II -- Interstellar phase transitions and the Dark Neutral Medium(1711.05506)

Nov. 15, 2017 astro-ph.HE
H$_{\rm I}$ 21-cm and $^{12}$CO 2.6-mm line emissions respectively trace the atomic and molecular gas phases, but they miss most of the opaque H$_{\rm I}$ and diffuse H$_2$ present in the Dark Neutral Medium (DNM) at the transition between the H$_{\rm I}$-bright and CO-bright regions. Jointly probing H$_{\rm I}$, CO, and DNM gas, we aim to constrain the threshold of the H$_{\rm I}$-H$_2$ transition. We also aim to measure gas mass fractions in the different phases and to test their relation to cloud properties. We have used dust optical depth measurements at 353 GHz, $\gamma$-ray maps at GeV energies, and H$_{\rm I}$ and CO line data to trace the gas column densities and map the DNM in nearby clouds toward the Galactic anticentre and Chamaeleon regions. The onset of detectable CO intensities varies by only a factor of 4 from cloud to cloud, between 0.6$\times 10^{21}$ cm$^{-2}$ and 2.5$\times 10^{21}$ cm$^{-2}$ in total gas column density. We observe larger H$_2$ column densities than linearly inferred from the CO intensities at $A_{\rm{V}}> 3$ magnitudes because of the large CO optical thickness: the additional H$_2$ mass in this regime represents on average 20$\%$ of the CO-inferred molecular mass. In the DNM envelopes, we find that the fraction of diffuse CO-dark H$_2$ in the molecular column densities decreases with increasing $A_{\rm{V}}$ in a cloud. For a half molecular DNM, the fraction decreases from more than 80$\%$ at 0.4 magnitude to less than 20$\%$ beyond 2 magnitudes. In mass, the DNM fraction varies with the cloud properties. The mass present in the DNM envelopes appears to scale with the molecular mass seen in CO as $M_{\rm H}^{\rm DNM} = (62 \pm 7) {M_{\rm H_2}^{\rm CO}}^{\,0.51 \pm 0.02}$ across two decades in mass. The phase transitions in these clouds show both common trends and environmental differences.
• 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.
• ### Cosmic rays, gas and dust in nearby anticentre clouds : I -- CO-to-H2 conversion factors and dust opacities(1703.05237)

March 15, 2017 astro-ph.HE
We aim to explore the capabilities of dust emission and rays for probing the properties of the interstellar medium in the nearby anti-centre region, using gamma-ray observations with the Fermi Large Area Telescope (LAT), and the thermal dust optical depth inferred from Planck and IRAS observations. In particular, we aim at quantifying potential variations in cosmic-ray density and dust properties per gas nucleon across the different gas phases and different clouds, and at measuring the CO-to-H2 conversion factor, X$_{CO}$ , in different environments. We have separated six nearby anti-centre clouds that are coherent in velocities and distances, from the Galactic-disc background in HI 21-cm and $^{12}$CO 2.6-mm line emission. We have jointly modelled the gamma-ray intensity recorded between 0.4 and 100 GeV, and the dust optical depth at 353 GHz as a combination of HI-bright, CO-bright, and ionised gas components. The complementary information from dust emission and gamma rays was used to reveal the gas not seen, or poorly traced, by HI , free-free, and $^{12}$CO emissions, namely (i) the opaque HI and diffuse H$_2$ present in the Dark Neutral Medium at the atomic-molecular transition, and (ii) the dense H$_2$ to be added where $^{12}$CO lines saturate. The measured interstellar gamma-ray spectra support a uniform penetration of the cosmic rays with energies above a few GeV through the clouds. We find a gradual increase in grain opacity as the gas becomes more dense. The increase reaches a factor of four to six in the cold molecular regions that are well shielded from stellar radiation. Consequently, the X$_{CO}$ factor derived from dust is systematically larger by 30% to 130% than the gamma-ray estimate. We also evaluate the average gamma-ray X$_{CO}$ factorfor each cloud, and find that X$_{CO}$ tends to decrease from diffuse to more compact molecular clouds, as expected from theory.
• 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.
• ### Galactic cold cores VI. Dust opacity spectral index(1509.08023)

Sept. 26, 2015 astro-ph.GA
The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the whole sample and inside individual fields. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of beta=1.84. We observe a clear T-beta anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and beta(FIR) can rise to ~2.2 in individual clumps. The Planck 217 GHz band shows a systematic excess that is consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 um, the median spectral index values are beta(FIR) ~ 1.91 and beta(mm) ~ 1.66. The spectral index changes as a function of column density and wavelength. Beta variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater.
• ### Galactic cold cores V. Dust opacity(1501.07092)

Jan. 28, 2015 astro-ph.GA
The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observations of the reddening of background stars are converted into near-infrared optical depth. We studied the systematic errors affecting these parameters and used modelling to correct for the expected biases. The ratio of 250um and J band opacities is correlated with the cloud location and star formation activity. We find a median ratio of tau(250um)/tau(J)= (1.6+-0.2)*10^-3, which is more than three times the mean value in diffuse medium. Assuming a spectral index beta=1.8 instead of beta=2.0, the value would be lower by ~30%. No significant systematic variation is detected with Galactocentric distance or with Galactic height. The tau(250um)/tau(J) maps reveal six fields with clear increase of submillimetre opacity of up to tau(250um)/tau(J) ~ 4*10^-3. These are all nearby fields with spatially resolved clumps of high column density. We interpret the increase in the far-infrared opacity as a sign of grain growth in the densest and coldest regions of interstellar clouds.
• ### Galactic cold cores II. Herschel study of the extended dust emission around the first Planck detections(1101.3003)

Jan. 15, 2011 astro-ph.GA
Within the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold interstellar clouds detected with the Planck satellite. The three fields observed as part of the Herschel science demonstration phase (SDP) provided the first glimpse into the nature of these sources. We examine the properties of the dust emission within the fields. We determine the dust sub-millimetre opacity, look for signs of spatial variations in the dust spectral index, and estimate how the apparent variations of the parameters could be affected by different sources of uncertainty. We use the Herschel observations where the zero point of the surface brightness scale is set with the help of the Planck satellite data. We derive the colour temperature and column density maps of the regions and determine the dust opacity by a comparison with extinction measurements. By simultaneously fitting the colour temperature and the dust spectral index values we look for spatial variations in the apparent dust properties. With a simple radiative transfer model we estimate to what extent these can be explained by line-of-sight temperature variations, without changes in the dust grain properties. The analysis of the dust emission reveals cold and dense clouds that coincide with the Planck sources and confirm those detections. The derived dust opacity varies in the range kappa(250um) ~ 0.05-0.2 cm^2/g, higher values being observed preferentially in regions of high column density. The average dust spectral index beta is ~ 1.9-2.2. There are indications that beta increases towards the coldest regions. The spectral index decreases strongly near internal heating sources but, according to radiative transfer models, this can be explained by the line-of-sight temperature variations without a change in the dust properties.
• ### Distribution and characteristics of Infrared Dark Clouds using genetic forward modelling(0908.3851)

Aug. 26, 2009 astro-ph.GA
Infrared Dark Clouds (IRDCs) are dark clouds seen in silhouette in mid-infrared surveys. They are thought to be the birthplace of massive stars, yet remarkably little information exists on the properties of the population as a whole (e.g. mass spectrum, spatial distribution). Genetic forward modelling is used along with the Two Micron All Sky Survey and the Besancon Galactic model to deduce the three dimensional distribution of interstellar extinction towards previously identified IRDC candidates. This derived dust distribution can then be used to determine the distance and mass of IRDCs, independently of kinematic models of the Milky Way. Along a line of sight that crosses an IRDC, the extinction is seen to rise sharply at the distance of the cloud. Assuming a dust to gas ratio, the total mass of the cloud can be estimated. The method has been successfully applied to 1259 IRDCs, including over 1000 for which no distance or mass estimate currently exists. The IRDCs are seen to lie preferentially along the spiral arms and in the molecular ring of the Milky Way, reinforcing the idea that they are the birthplace of massive stars. Also, their mass spectrum is seen to follow a power law with an index of -1.75 +/- 0.06, steeper than giant molecular clouds in the inner Galaxy, but comparable to clumps in GMCs. This slope suggests that the IRDCs detected using the present method are not gravitationally bound, but are rather the result of density fluctuations induced by turbulence.
• ### The large scale dust lanes of the Galactic bar(0711.2471)

Nov. 15, 2007 astro-ph
(abridged) By comparing the distribution of dust and gas in the central regions of the Galaxy, we aim to obtain new insights into the properties of the offset dust lanes leading the bar's major axis in the Milky Way. On the one hand, the molecular emission of the dust lanes is extracted from the observed CO l-b-V distribution according to the interpretation of a dynamical model. On the other hand, a three dimensional extinction map of the Galactic central region constructed from near-infrared observations is used as a tracer of the dust itself and clearly reveals dust lanes in its face-on projection. Comparison of the position of both independent detections of the dust lanes is performed in the (l, b) plane. These two completely independent methods are used to provide a coherent picture of the dust lanes in the Milky Way bar. In both the gas and dust distributions, the dust lanes are found to be out of the Galactic plane, appearing at negative latitudes for l > 0 deg and at positive latitudes for l < 0 deg. However, even though there is substantial overlap between the two components, they are offset from one another with the dust appearing to lie closer to the b = 0 deg plane. Two scenarios are proposed to explain the observed offset. The first involves grain destruction by the bar shock and reformation downstream. Due to the decrease in velocity caused by the shock, this occurs at lower z. The second assumes that the gas and dust remain on a common tilted plane, but that the molecular gas decouples from the Milky Way's magnetic field, itself strong enough to resist the shear of the bar's shock. The diffuse gas and dust remain coupled to the field and are carried further downstream. This second scenario has recently been suggested in order to explain observations of the barred galaxy NGC 1097.
• ### Modelling the Galactic Interstellar Extinction Distribution in Three Dimensions(astro-ph/0604427)

April 20, 2006 astro-ph
The Two Micron All Sky Survey, along with the Stellar Population Synthesis Model of the Galaxy, developed in Besancon, is used to calculate the extinction distribution along different lines of sight. By combining many lines of sight, the large scale distribution of interstellar material can be deduced. The Galaxy model is used to provide the intrinsic colour of stars and their probable distances, so that the near infrared colour excess, and hence the extinction, may be calculated and its distance evaluated. Such a technique is dependent on the model used, however we are able to show that moderate changes in the model parameters result in insignificant changes in the predicted extinction. This technique has now been applied to over 64000 lines of sight, each separated by 15 arcmin, in the inner Galaxy (|l|<=100 deg, |b|<=10 deg). Using our extinction map, we have derived the main characteristics of the large scale structure of the dust distribution: scale height and warp of the ISM disc as well as the angle of the dust in the Galactic Bar. This resulting extinction map will be useful for studies of the inner Galaxy and its stellar populations.