
We highlight the power of the Gaia DR2 in studying many fine structures of
the HertzsprungRussell 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, postAGB stars, planetary nebulae, and white dwarfs are also
analysed, as well as lowmass 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
HRDbased studies.

Context. Stellar physical properties of star clusters are poorly known. Aims.
Our goals are to perform a spectrophotometric study of the B star population in
open clusters to derive accurate stellar parameters, search for the presence of
circumstellar envelopes, and discuss the characteristics of these stars.
Methods. The BCD spectrophotometric system is a powerful tool to obtain
fundamental parameters and infer the main properties of open clusters: distance
modulus, color excess, and age. We inspected the Balmer discontinuity to seek
circumstellar disks and identify Bestar candidates. High resolution spectra in
the H$\alpha$ region are used to confirm the Be nature. Results. We provide
Teff , log g, Mv , Mbol and spectral types for a sample of 68 stars in the
field of the open clusters NGC 6087, NGC 6250, NGC 6383, and NGC 6530, as well
as the cluster distances, ages and reddening. Then, based on a sample of 230 B
stars in the direction of the 11 open clusters studied along this series of
three papers, we report 6 new Be stars, 4 blue straggler candidates, and 15
Btype stars (called Bdd) with a double Balmer discontinuity. We also find that
the majority of the Be stars are dwarfs and present a maximum at the spectral
type B2B4 in young and intermediateage open clusters. Another maximum of Be
stars is observed at the spectral type B6B8 in open clusters older than 40
Myr, where the population of Bdd stars also becomes relevant. Conclusions. Our
results support previous statements that the Be phenomenon is present along the
whole main sequence band and occurs in very different evolutionary states. We
find clear evidence of an increase of stars with circumstellar envelopes with
cluster age. The Be phenomenon reaches its maximum in clusters of intermediate
age and the number of B stars with circumstellar envelopes (Be+Bdd stars) is
also high for the older clusters.

The interpretation of stellar apparent fundamental parameters must be treated
consistently with the characteristics of their surface rotation law. We develop
a model to determine the distribution of the effective temperature and gravity,
which depend on the surface differential rotation law and on the stellar
external geometry. The basic assumptions are: a) the atmosphric layers are in
radiative equilibrium; b) the bolometric flux is antiparallel with the
effective gravity; c) the angular velocity in the surface is
Omega=Omega(theta,alpha,k)] where (alpha,k) are free parameters. The effective
temperature varies with colatitude `theta' and on the differentialrotation
law through the effective gravity and the gravitydarkening function (GDF).
Although the derived expressions can be treated numerically, for some low
integer values of k, analytical forms of the `integral ofcharacteristic
curves', on which the determination of the GDF relies, are obtained. The
effects of the quantities (eta,alpha,k) (eta=ratio between centrifugal and
gravitational accelerations at the equator) on the determination of the Vsini
parameter and on the `gravitydarkening exponent' are studied. Depending on the
values of (eta,alpha,k) the velocity V in the derived Vsini may strongly
deviate from the equatorial rotational velocity. It is shown that the von
Zeipel'slike gravitydarkening exponent $\beta_1$ depends on all parameters
$(\eta,\alpha,k)$ and that its value also depends on the viewingangle i.
Hence, there no unique interpretation of this exponent determined empirically
in terms of (i,alpha). The data on rotating stars should be analyzed by taking
into account the rotational effects through the GDF, by assuming k=2 as a first
approximation. Instead of the classical pair (eta,beta1), it would be more
useful to determine the quantities (eta,alpha,i) to characterize stellar
rotation.

We assume that stars may undergo surface differential rotation to study its
impact on the interpretation of $V\!\sin i$ and on the observed distribution
$\Phi(u)$ of ratios of true rotational velocities $u=V/V_\rm c$ ($V_\rm c$ is
the equatorial critical velocity). We discuss some phenomena affecting the
formation of spectral lines and their broadening, which can obliterate the
information carried by $V\!\sin i$ concerning the actual stellar rotation. We
studied the line broadening produced by several differential rotational laws,
but adopted Maunder's expression
$\Omega(\theta)=\Omega_o(1+\alpha\cos^2\theta)$ as an attempt to account for
all of these laws with the lowest possible number of free parameters. We
studied the effect of the differential rotation parameter $\alpha$ on the
measured $V\!\sin i$ parameter and on the distribution $\Phi(u)$ of ratios
$u=V/V_\rm c$. We conclude that the inferred $V\!\sin i$ is smaller than
implied by the actual equatorial linear rotation velocity $V_\rm eq$ if the
stars rotate with $\alpha<0$, but is larger if the stars have $\alpha>0$. For a
given $\alpha$ the deviations of $V\!\sin i$ are larger when $\alpha<0$. If
the studied Be stars have on average $\alpha<0$, the number of rotators with
$V_\rm eq\simeq0.9V_\rm c$ is larger than expected from the observed
distribution $\Phi(u)$; if these stars have on average $\alpha>0$, this number
is lower than expected. We discuss seven phenomena that contribute either to
narrow or broaden spectral lines, which blur the information on the rotation
carried by $V\!\sin i$ and, in particular, to decide whether the Be phenomenon
mostly rely on the critical rotation. We show that twodimensional radiation
transfer calculations are needed in rapid rotators to diagnose the stellar
rotation more reliably.

Gaia is a space mission currently measuring the five astrometric parameters
as well as spectrophotometry of at least 1 billion stars to G = 20.7 mag with
unprecedented precision. The sixth parameter in phase space (radial velocity)
is also measured thanks to mediumresolution spectroscopy being obtained for
the 150 million brightest stars. During the commissioning phase, two fields,
one around each ecliptic pole, have been repeatedly observed to assess and to
improve the overall satellite performances as well as the associated reduction
and analysis software. A groundbased photometric and spectroscopic survey was
therefore initiated in 2007, and is still running in order to gather as much
information as possible about the stars in these fields. This work is of
particular interest to the validation of the Radial Velocity Spectrometer (RVS)
outputs. The paper presents the radial velocity measurements performed for the
Southern targets in the 12  17 R magnitude range on high to midresolution
spectra obtained with the GIRAFFE and UVES spectrographs.

We present a semiautomatic procedure to obtain fundamental physical
parameters and distances of classical Be (CBe) stars, based on the
BarbierChalongeDivan (BCD) spectrophotometric system. Our aim is to apply
this procedure to a large sample of CBe stars detected by the IPHAS photometric
survey, to determine their fundamental physical parameters and to explore their
suitability as galactic structure tracers. In this paper we describe the
methodology used and the validation of the procedure by comparing our results
with those obtained from different independent astrophysical techniques for
subsamples of stars in common with other studies. We also present a test case
study of the galactic structure in the direction of the Perseus Galactic Arm,
in order to compare our results with others recently obtained with different
techniques and the same sample of stars. We did not find any significant
clustering of stars at the expected positions of the Perseus and Outer Galactic
Arms, in agreement with previous studies in the same area that we used for
verification.

The Gaia satellite will survey the entire celestial sphere down to 20th
magnitude, obtaining astrometry, photometry, and low resolution
spectrophotometry on one billion astronomical sources, plus radial velocities
for over one hundred million stars. Its main objective is to take a census of
the stellar content of our Galaxy, with the goal of revealing its formation and
evolution. Gaia's unique feature is the measurement of parallaxes and proper
motions with hitherto unparalleled accuracy for many objects. As a survey, the
physical properties of most of these objects are unknown. Here we describe the
data analysis system put together by the Gaia consortium to classify these
objects and to infer their astrophysical properties using the satellite's data.
This system covers single stars, (unresolved) binary stars, quasars, and
galaxies, all covering a wide parameter space. Multiple methods are used for
many types of stars, producing multiple results for the end user according to
different models and assumptions. Prior to its application to real Gaia data
the accuracy of these methods cannot be assessed definitively. But as an
example of the current performance, we can attain internal accuracies (RMS
residuals) on F,G,K,M dwarfs and giants at G=15 (V=1517) for a wide range of
metallicites and interstellar extinctions of around 100K in effective
temperature (Teff), 0.1mag in extinction (A0), 0.2dex in metallicity ([Fe/H]),
and 0.25dex in surface gravity (logg). The accuracy is a strong function of the
parameters themselves, varying by a factor of more than two up or down over
this parameter range. After its launch in November 2013, Gaia will nominally
observe for five years, during which the system we describe will continue to
evolve in light of experience with the real data.

In previous works of this series, we have shown that late B and early Atype
stars have genuine bimodal distributions of rotational velocities and that late
Atype stars lack slow rotators. The distributions of the surface angular
velocity ratio \Omega/\Omega_crit (\Omega_crit is the critical angular
velocity) have peculiar shapes according to spectral type groups, which can be
caused by evolutionary properties. We aim to review the properties of these
rotational velocity distributions in some detail as a function of stellar mass
and age. We have gathered v sin i for a sample of 2014 B6 to F2type stars. We
have determined the masses and ages for these objects with stellar evolution
models. The (Teff, log L/Lsun)parameters were determined from the uvby\beta
photometry and the HIPPARCOS parallaxes. The velocity distributions show two
regimes that depend on the stellar mass. Stars less massive than 2.5 Msun have
a unimodal equatorial velocity distribution and show a monotonical acceleration
with age on the main sequence (MS). Stars more massive have a bimodal
equatorial velocity distribution. Contrarily to theoretical predictions, the
equatorial velocities of stars from about 1.7 Msun to 3.2 Msun undergo a strong
acceleration in the first third of the MS evolutionary phase, while in the last
third of the MS they evolve roughly as if there were no angular momentum
redistribution in the external stellar layers. The studied stars might start in
the ZAMS not necessarily as rigid rotators, but with a total angular momentum
lower than the critical one of rigid rotators. The stars seem to evolve as
differential rotators all the way of their MS life span and the variation of
the observed rotational velocities proceeds with characteristic time scales
\delta(t)\sim 0.2 t_MS, where t_MS is the time spent by a star in the MS.

Aims. What is the origin of the largeamplitude variability in WolfRayet WN8
stars in general and WR123 in particular? A dedicated spectroscopic campaign
targets the tenhour period previously found in the highprecision photometric
data obtained by the MOST satellite. Methods. In JuneAugust 2003 we obtained a
series of high signaltonoise, midresolution spectra from several sites in
the {\lambda}{\lambda} 4000  6940 A^{\circ} domain. We also followed the star
with occasional broadband (Johnson V) photometry. The acquired spectroscopy
allowed a detailed study of spectral variability on timescales from \sim 5
minutes to months. Results. We find that all observed spectral lines of a given
chemical element tend to show similar variations and that there is a good
correlation between the lines of different elements, without any significant
time delays, save the strong absorption components of the Hei lines, which tend
to vary differently from the emission parts. We find a single sustained
periodicity, P \sim 9.8 h, which is likely related to the relatively stable
pulsations found in MOST photometry obtained one year later. In addition,
seemingly stochastic, largeamplitude variations are also seen in all spectral
lines on timescales of several hours to several days.

We present the results from the spectroscopic followup of WR140 (WC7 + O45)
during its last periastron passage in January 2009. This object is known as the
archetype of colliding wind binaries and has a relatively large period (~ 8
years) and eccentricity (~ 0.89). We provide updated values for the orbital
parameters, new estimates for the WR and O star masses and new constraints on
the massloss rates.

Since the external regions of the envelopes of rapidly rotating earlytype
stars are unstable to convection, a coupling may exist between the convection
and the internal rotation. We explore what can be learned from spectroscopic
and interferometric observations about the properties of the rotation law in
the external layers of these objects. Using simple relations between the
entropy and specific rotational quantities, some of which are found to be
efficient at accounting for the solar differential rotation in the convective
region, we derived analytical solutions that represent possible differential
rotations in the envelope of earlytype stars. A surface latitudinal
differential rotation may not only be an external imprint of the inner
rotation, but induces changes in the stellar geometry, the gravitational
darkening, the aspect of spectral line profiles, and the emitted spectral
energy distribution. By studying the equation of the surface of stars with
nonconservative rotation laws, we conclude that objects undergo geometrical
deformations that are a function of the latitudinal differential rotation able
to be scrutinized both spectroscopically and by interferometry. The combination
of Fourier analysis of spectral lines with model atmospheres provides
independent estimates of the surface latitudinal differential rotation and the
inclination angle. Models of stars at different evolutionary stages rotating
with internal conservative rotation laws were calculated to show that the Roche
approximation can be safely used to account for the gravitational potential.
The surface temperature gradient in rapid rotators induce an acceleration to
the surface angular velocity. A nonzero differential rotation parameter may
indicate that the rotation is neither rigid nor shellular underneath the
stellar surface.

A variety of physical processes leading to different types of pulsations and
chemical compositions is observed among A and Ftype stars. To investigate the
underlying mechanisms responsible for these processes in stars with similar
locations in the HR diagram, an accurate abundance determination is needed,
among others. Here, we describe a semiautomatic procedure developed to
determine chemical abundances of various elements ranging from helium to
mercury for this type of stars. We test our procedure on synthetic spectra,
demonstrating that our procedure provides abundances consistent with the input
values, even when the stellar parameters are offset by reasonable observational
errors. For a fastrotating star such as Vega, our analysis is consistent with
those carried out with other planeparallel model atmospheres. Simulations show
that the offsets from the input abundances increase for stars with low
inclination angle of about 4 degrees. For this inclination angle, we also show
that the distribution of the iron abundance found in different regions is
bimodal. Furthermore, the effect of rapid rotation can be seen in the peculiar
behaviour of the H_beta line.

Effective temperatures of earlytype supergiants are important to test
stellar atmosphere and internal structuremodels of massive and intermediate
mass objects at different evolutionary phases. However, these Teff values are
more or less discrepant depending on the method used to determine them. We aim
to obtain a new calibration of the Teff parameter for earlytype supergiants as
a function of observational quantities that are highly sensitive to the
ionization balance in the photosphere and its gas pressure, independent of the
interstellar extinction, and as much as possible modelindependent. The
observational quantities that best address our aims are the (lambda1, D)
parameters of the BCD spectrophotometric system. They describe the energy
distribution around the Balmer discontinuity, which is highly sensitive to Teff
and logg. We perform a calibration of the (lambda1, D) parameters into Teff
using effective temperatures derived with the bolometricflux method for 217
program stars, whose individual absolute relative uncertainties are on average
0.05. The final comparison of calculated with obtained Teff values in the
(lambda1,D) calibration further show that the latter have total relative
uncertainties, which on average are of the order of 0.05 for all spectral types
and luminosity classes. The effective temperatures of OB supergiants derived in
this work agree on average within some 2000 K with other determinations found
in the literature, except those issued from windfree nonLTE planeparallel
models of stellar atmospheres, which produce effective temperatures that can be
overestimated by up to more than 5000 K near Teff = 25000 K.

High resolution synthetic stellar libraries are of fundamental importance for
the preparation of the Gaia Mission. We present new sets of spectral stellar
libraries covering two spectral ranges: 300 1100 nm at 0.1 nm resolution, and
840  890 nm at 0.001 nm resolution. These libraries span a large range in
atmospheric parameters, from supermetalrich to very metalpoor (5.0
$<$[Fe/H]$<$+1.0), from cool to hot (\teff=300050000 K) stars, including
peculiar abundance variations. The spectral resolution, spectral type coverage
and number of models represent a substantial improvement over previous
libraries used in population synthesis models and in atmospheric analysis.

We carried out low resolution spectroscopic observations in the wavelength
range 34004700 A of 20 Heweak and 8 Hestrong stars to determine their
fundamental parameters by means of the DivanChalongeBarbier (BCD)
spectrophotometric system. For a few Heweak stars we also estimate the
effective temperatures and the angular diameters by integrating absolute fluxes
observed over a wide spectral range. NonLTE model calculations are carried out
to study the influence of the He/H abundance ratio on the emergent radiation of
Hestrong stars and on their Teff determination. We find that the effective
temperatures, surface gravities and bolometric absolute magnitudes of Heweak
stars estimated with the BCD system and the integrated flux method are in good
agreement between each other, and they also agree with previous determinations
based on several different methods. The mean discrepancy between the visual
absolute magnitudes derived using the Hipparcos parallaxes and the BCD values
is on average 0.3 mag for Heweak stars, while it is 0.5 mag for Hestrong
stars. For Hestrong stars, we note that the BCD calibration, based on stars in
the solar environment, leads to overestimated values of Teff. By means of model
atmosphere calculations with enhanced He/H abundance ratios we show that larger
He/H ratios produce smaller BD which naturally explains the Teff
overestimation. We take advantage of these calculations to introduce a method
to estimate the He/H abundance ratio in Hestrong stars. The BD of HD 37479
suggests that the Teff of this star remains fairly constant as the star
spectrum undergoes changes in the intensity of H and He absorption lines. Data
for the Hestrong star HD 66765 are reported for the first time.

The new mathematical framework based on the free energy of pure classical
fluids presented in [R. D. Rohrmann, Physica A 347, 221 (2005)] is extended to
multicomponent systems to determine thermodynamic and structural properties of
chemically complex fluids. Presently, the theory focuses on $D$dimensional
mixtures in the lowdensity limit (packing factor $\eta < 0.01$). The formalism
combines the freeenergy minimization technique with space partitions that
assign an available volume $v$ to each particle. $v$ is related to the
closeness of the nearest neighbor and provides an useful tool to evaluate the
perturbations experimented by particles in a fluid. The theory shows a close
relationship between statistical geometry and statistical mechanics. New,
unconventional thermodynamic variables and mathematical identities are derived
as a result of the space division. Thermodynamic potentials $\mu_{il}$,
conjugate variable of the populations $N_{il}$ of particles class $i$ with the
nearest neighbors of class $l$ are defined and their relationships with the
usual chemical potentials $\mu_i$ are established. Systems of hard spheres are
treated as illustrative examples and their thermodynamics functions are derived
analytically. The lowdensity expressions obtained agree nicely with those of
scaledparticle theory and PercusYevick approximation. Several pair
distribution functions are introduced and evaluated. Analytical expressions are
also presented for hard spheres with attractive forces due to K\^actails and
squarewell potentials. Finally, we derive general chemical equilibrium
conditions.

We show that Be stars belong to a high velocity tail of a single Btype star
rotational velocity distribution in the MS. This implies that: 1) the number
fraction N(Be)/N(Be+B) is independent of the mass; 2) Bn stars having ZAMS
rotational velocities higher than a given limit might become Be stars.

Using 2D models of rotating stars, the interferometric measurements of alpha
Eri and its fundamental parameters corrected for gravitational darkening
effects we infer that the star might have a core rotating 2.7 times faster than
the surface. We explore the consequences on spectral lines produced by surface
differential rotation combined with the effects due to a kind of internal
differential rotation with rotational energies higher than allowed for rigid
rotation which induce geometrical deformations that do not distinguish strongly
from those carried by the rigid rotation.

A sample of 97 galactic field Be stars were studied by taking into account
the effects induced by the fast rotation on their fundamental parameters. All
program stars were observed in the BCD spectrophotometric system in order to
minimize the perturbations produced by the circumstellar environment on the
spectral photospheric signatures. This is one of the first attempts at
determining stellar masses and ages by simultaneously using model atmospheres
and evolutionary tracks, both calculated for rotating objects. The stellar ages
($\tau$) normalized to the respective inferred time that each rotating star can
spend in the main sequence phase ($\tau\_{\rm MS}$) reveal a massdependent
trend. This trend shows that: a) there are Be stars spread over the whole
interval $0 \la \tau/\tau\_{\rm MS} \la 1$ of the main sequence evolutionary
phase; b) the distribution of points in the ($\tau/\tau\_{\rm
MS},M/M\_{\odot}$) diagram indicates that in massive stars ($M \ga
12M\_{\odot}$) the Be phenomenon is present at smaller $\tau/\tau\_{\rm MS}$
age ratios than for less massive stars ($M \la 12M\_{\odot}$). This
distribution can be due to: $i$) higher massloss rates in massive objets,
which can act to reduce the surface fast rotation; $ii$) circulation time
scales to transport angular momentum from the core to the surface, which are
longer the lower the stellar mass.

This work is the second part of the set of measurements of vsini for Atype
stars, begun by Royer et al. (2002). Spectra of 249 B8 to F2type stars
brighter than V=7 have been collected at Observatoire de HauteProvence (OHP).
Fourier transforms of several line profiles in the range 42004600 A are used
to derive vsini from the frequency of the first zero. Statistical analysis of
the sample indicates that measurement error mainly depends on vsini and this
relative error of the rotational velocity is found to be about 5% on average.
The systematic shift with respect to standard values from Slettebak et al.
(1975), previously found in the first paper, is here confirmed. Comparisons
with data from the literature agree with our findings: vsini values from
Slettebak et al. are underestimated and the relation between both scales
follows a linear law: vsini(new) = 1.03 vsini(old) + 7.7. Finally, these data
are combined with those from the previous paper (Royer et al. 2002), together
with the catalogue of Abt & Morrell (1995). The resulting sample includes some
2150 stars with homogenized rotational velocities.