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 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability 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.

The ESA Gaia mission provides a unique time-domain survey for more than one billion sources brighter than G=20.7 mag. Gaia offers the unprecedented opportunity to study variability phenomena in the Universe thanks to multi-epoch G-magnitude photometry in addition to astrometry, blue and red spectro-photometry, and spectroscopy. Within the Gaia Consortium, Coordination Unit 7 has the responsibility to detect variable objects, classify them, derive characteristic parameters for specific variability classes, and provide global descriptions of variable phenomena. We describe the variability processing and analysis that we plan to apply to the successive data releases, and we present its application to the G-band photometry results of the first 14 months of Gaia operations that comprises 28 days of Ecliptic Pole Scanning Law and 13 months of Nominal Scanning Law. Out of the 694 million, all-sky, sources that have calibrated G-band photometry in this first stage of the mission, about 2.3 million sources that have at least 20 observations are located within 38 degrees from the South Ecliptic Pole. We detect about 14% of them as variable candidates, among which the automated classification identified 9347 Cepheid and RR Lyrae candidates. Additional visual inspections and selection criteria led to the publication of 3194 Cepheid and RR Lyrae stars, described in Clementini et al. (2016). Under the restrictive conditions for DR1, the completenesses of Cepheids and RR Lyrae stars are estimated at 67% and 58%, respectively, numbers that will significantly increase with subsequent Gaia data releases. Data processing within the Gaia Consortium is iterative, the quality of the data and the results being improved at each iteration. The results presented in this article show a glimpse of the exceptional harvest that is to be expected from the Gaia mission for variability phenomena. [abridged]

The nova outburst of V407 Cyg in 2010 Mar. 10 was the first observed for this star but its close resemblance to the well known symbiotic-like recurrent nova RS Oph suggests that it is also a member of this rare type of Galactic novae. The nova was the first detected at $\gamma$-ray energies and is the first known nova explosion for this system. The extensive multiwavelength coverage of this outburst makes it an ideal comparison with the few other outbursts known for similar systems. We extend our previous analysis of the Mira and the expanding shock from the explosion to detail the time development of the photoionized Mira wind, circumstellar medium, and shocked circumstellar environment to derive their physical parameters and how they relate to large scale structure of the environment, extending the previous coverage to more than 500 days after outburst. Absorption lines of Fe-peak ions formed in the Mira wind were visible as P Cyg profiles at low velocity before Day 69, around the time of the X-ray peak and we identified many absorption transitions without accompanying emission for metal lines. The H Balmer lines showed strong P Cyg absorption troughs that weakened during the 2010 observing period, through Day 128. We distinguish the components from the shock, the photoionized environment, and the chromosphere and inner Mira wind using spectra taken more than one year after outburst. The multiple shells and radiative excitation phenomenology are similar to those recently cited for GRBs and SNIa (severely truncated)

H{\alpha} emission V/R variations caused by a discontinous mass transfer in interacting binaries with a rapidly rotating accreting star are modelled qualititatively for the first time. The program ZEUS-MP was used for a non-linear 3-D hydrodynamical modeling of a development of a blob of gaseous material injected into an orbit around a star. It resulted in the formation of an elongated disk with a slow prograde revolution. The LTE radiative transfer program SHELLSPEC was used to calculate the H{\alpha} profiles originating in the disk for several phases of its revolution. The profiles have the form of a double emission and exhibit V/R and radial velocity variations. However, these variations should be a temporal phenomenon since imposing a viscosity in given model would lead to a circularization of the disk and fading-out of given variations.

On 2010 Mar 10, V407 Cyg was discovered in outburst, eventually reaching V< 8 and detected by Fermi. Using medium and high resolution ground-based optical spectra, visual and Swift UV photometry, and Swift X-ray spectrophotometry, we describe the behavior of the high-velocity profile evolution for this nova during its first three months. The peak of the X-ray emission occurred at about day 40 with a broad maximum and decline after day 50. The main changes in the optical spectrum began at around that time. The He II 4686A line first appeared between days 7 and 14 and initially displayed a broad, symmetric profile that is characteristic of all species before day 60. Low-excitation lines remained comparatively narrow, with v(rad,max) of order 200-400 km/s. They were systematically more symmetric than lines such as [Ca V], [Fe VII], [Fe X], and He II, all of which showed a sequence of profile changes going from symmetric to a blue wing similar to that of the low ionization species but with a red wing extended to as high as 600 km/s . The Na I D doublet developed a broad component with similar velocity width to the other low-ionization species. The O VI Raman features were not detected. We interpret these variations as aspherical expansion of the ejecta within the Mira wind. The blue side is from the shock penetrating into the wind while the red wing is from the low-density periphery. The maximum radial velocities obey power laws, v(rad,max) t^{-n} with n ~ 1/3 for red wing and ~0.8 for the blue. (truncated)

AG Dra is one of a small group of low metallicity S-type symbiotic binaries with K-type giants that undergoes occasional short-term outbursts of unknown origin. Our aim is to study the behavior of the white dwarf during an outburst using the optical Raman lines and other emission features in the red giant wind. The goal is to determine changes in the envelope and the wind of the gainer in this system during a major outburst event and to study the coupling between the UV and optical during a major outburst. Using medium and high resolution groundbased optical spectra and comparisons with archival $FUSE$ and $HST/STIS$ spectra, we study the evolution of the Raman O VI features and the Balmer, He I, and He II lines during the outburst from 2006 Sept. through 2007 May and include more recent observations (2009) to study the subsequent evolution of the source. The O VI Raman features disappeared completely at the peak of the major outburst and the subsequent variation differs substantially from that reported during the previous decade. The He I and He II lines, and the Balmer lines, vary in phase with the Raman features but there is a double-valuedness to the He I 6678, 7065 relative to the O VI Raman 6825\AA\ variations in the period between 2006-2008 that has not been previously reported. The variations in the Raman feature ratio through the outburst interval are consistent with the disappearance of the O VI FUV resonance wind lines from the white dwarf and of the surrounding O$^{+5}$ ionized region within the red giant wind provoked by the expansion and cooling of the white dwarf photosphere.

The eclipsing and double-lined spectroscopic binary V380 Cyg is an extremely important probe of stellar evolution: its primary component is a high-mass star at the brink of leaving the main sequence whereas the secondary star is still in the early part of its main sequence lifetime. We present extensive high-resolution echelle and grating spectroscopy from Ondrejov, Calar Alto, Victoria and La Palma. We apply spectral disentangling to unveil the individual spectra of the two stars and obtain new spectroscopic elements. The secondary star contributes only about 6% of the total light, which remains the main limitation to measuring the system's characteristics. We determine improved physical properties, finding masses 13.1 +/- 0.3 and 7.8 +/- 0.1 M_sun, radii 16.2 +/- 0.3 and 4.06 +/- 0.08 R_sun, and effective temperatures 21750 +/- 280 and 21600 +/- 550 K, for the primary and secondary components respectively. We perform a detailed abundance analysis by fitting non-LTE theoretical line profiles to the disentangled spectrum of the evolved primary star, and reveal an elemental abundance pattern reminiscent of a typical nearby B star. Contrary to the predictions of recent theoretical evolution models with rotational mixing, no trace of abundance modifications due to the CNO cycle are detected. No match can be found between the predictions of these models and the properties of the primary star: a mass discrepancy of 1.5 M_sun exists and remains unexplained.

We present an examination of high resolution, ultraviolet spectroscopy from Hubble Space Telescope of the photospheric spectrum of the O-supergiant in the massive X-ray binary HD 226868 = Cyg X-1. We analyzed this and ground-based optical spectra to determine the effective temperature and gravity of the O9.7 Iab supergiant. Using non-local thermodynamic equilibrium (non-LTE), line blanketed, plane parallel models from the TLUSTY grid, we obtain T_eff = 28.0 +/- 2.5kK and log g > 3.00 +/- 0.25, both lower than in previous studies. The optical spectrum is best fit with models that have enriched He and N abundances. We fit the model spectral energy distribution for this temperature and gravity to the UV, optical, and IR fluxes to determine the angular size of and extinction towards the binary. The angular size then yields relations for the stellar radius and luminosity as a function of distance. By assuming that the supergiant rotates synchronously with the orbit, we can use the radius - distance relation to find mass estimates for both the supergiant and black hole as a function of the distance and the ratio of stellar to Roche radius. Fits of the orbital light curve yield an additional constraint that limits the solutions in the mass plane. Our results indicate masses of 23^{+8}_{-6} M_sun for the supergiant and 11^{+5}_{-3} M_sun for the black hole.

We present results from Hubble Space Telescope UV spectroscopy of the massive X-ray binary system, HD226868 = Cyg X-1. The spectra were obtained at both orbital conjunction phases in two separate runs in 2002 and 2003 when the system was in the X-ray high/soft state. The stellar wind lines suffer large reductions in strength when the black hole is in the foreground due to the X-ray ionization of the wind ions. A comparison of HST and archival IUE spectra shows that similar photoionization effects occur in both the X-ray states. We constructed model UV wind line profiles assuming that X-ray ionization occurs everywhere in the wind except the zone where the supergiant blocks the X-ray flux. The good match between the observed and model profiles indicates that the wind ionization extends to near to the hemisphere of the supergiant facing the X-ray source. The H-alpha emission strength is generally lower in the high/soft state compared to the low/hard state, but the He II 4686 emission is relatively constant between states. The results suggest that mass transfer in Cyg X-1 is dominated by a focused wind flow that peaks along the axis joining the stars and that the stellar wind contribution is shut down by X-ray photoionization effects. The strong stellar wind from the shadowed side of the supergiant will stall when Coriolis deflection brings the gas into the region of X-ray illumination. This stalled gas component may be overtaken by the orbital motion of the black hole and act to inhibit accretion from the focused wind. The variations in the strength of the shadow wind component may then lead to accretion rate changes that ultimately determine the X-ray state.

We have developed a systematic procedure to study the disks in Algol-type binaries using spectroscopic analysis, synthetic spectra, and tomography. We analyzed 119 H-alpha spectra of TT Hya, an Algol-type eclipsing interacting binary, collected from 1985-2001. The new radial velocities enabled us to derive reliable orbital elements, including a small non-zero eccentricity, and to improve the accuracy of the absolute dimensions of the system. High resolution IUE spectra were also analyzed to study the formation of the ultraviolet lines and continuum. Synthetic spectra of the iron curtain using our new shellspec program enabled us to derive a characteristic disk temperature of 7000K. We have demonstrated that the UV emission lines seen during total primary eclipse cannot originate from the accretion disk, but most likely arise from a hotter disk-stream interaction region. The synthetic spectra of the stars, disk, and stream allowed us to derive a lower limit to the mass transfer rate of 2e-10 solar masses per year. Doppler tomography of the observed H-alpha profiles revealed a distinct accretion disk. The difference spectra produced by subtracting the synthetic spectra of the stars resulted in an image of the disk, which virtually disappeared once the composite synthetic spectra of the stars and disk were used to calculate the difference spectra. An intensity enhancement of the resulting tomogram revealed images of the gas stream and an emission arc. We successfully modeled the gas stream using shellspec and associated the emission arc with an asymmetry in the accretion disk.

We present the first K'-band, long-baseline interferometric observations of the northern Be stars gamma Cas, phi Per, zeta Tau, and kappa Dra. The measurements were made with multiple telescope pairs of the CHARA Array interferometer, and in every case the observations indicate that the circumstellar disks of the targets are resolved. We fit the interferometric visibilities with predictions from a simple disk model that assumes an isothermal gas in Keplerian rotation. We derive fits of the four model parameters (disk base density, radial density exponent, disk normal inclination, and position angle) for each of the targets. The resulting densities are in broad agreement with prior studies of the IR excess flux and the resulting orientations generally agree with those from interferometric H-alpha and continuum polarimetric observations. We find that the angular size of the K' disk emission is smaller than that determined for the H-alpha emission, and we argue that the difference is the result of a larger H-alpha opacity and the relatively larger neutral hydrogen fraction with increasing disk radius. All the targets are known binaries with faint companions, and we find that companions appear to influence the interferometric visibilities in the cases of phi Per and kappa Dra. We also present contemporaneous observations of the H-alpha, H-gamma, and Br-gamma emission lines. Synthetic model profiles of these lines that are based on the same disk inclination and radial density exponent as derived from the CHARA Array observations match the observed emission line strength if the disk base density is reduced by approximately 1.7 dex.

We have carried out a spectroscopic survey of several B, Be, and shell stars in optical and near-infrared regions. Line profiles of the H-alpha line and of selected Fe II and O I lines are presented.

Recent theoretical calculations of stellar evolutionary tracks for rotating high-mass stars suggests that the chemical composition of the surface layers changes even whilst the star is evolving on the Main Sequence. The abundance analysis of binary components with precisely known fundamental stellar quantities allows a powerful comparison with theory. The observed spectra of close binary stars can be separated into the individual spectra of the component stars using the method of spectral disentangling on a time-series of spectra taken over the orbital cycle. Recently, Pavlovski & Hensberge (2005, A&A, 439, 309) have shown that, even with moderately high line-broadening, metal abundances can be derived from disentangled spectra with a precision of 0.1 dex. In a continuation of this project we have undertaken a detailed abundance analysis of the components of another two high-mass binaries, V453 Cyg, and V380 Cyg. Both binaries are well-studied systems with modern solutions. The components are close to the TAMS and therefore very suitable for an observational test of early mixing in high-mass stars.

We present an analysis of new spectroscopic observations of the bright Be star kappa Dra obtained at the Ondrejov observatory during 1992 - 2003 and UBV photometric observations secured at several observatories. The fundamental stellar parameters have been obtained from a comparison with a grid of NLTE model atmospheres. The best fit was found for Teff=14000K, log g = 3.5, and v sin i = 170km/s. These values together with a Hipparcos parallax lead to a stellar mass M=4.8+-0.8 M_solar and radius R=6.4+-0.5 R_solar. Long-term variations of kappa Dra were analysed using measurements of equivalent widths, central intensities, peak intensities of emission lines and emission peak velocity differences for H-alpha, H-beta, H-gamma, H-delta, and some helium, silicon, and iron lines. It turned out that the previously reported period of 23 years in the variation of the emission strength is probably a cyclic, not a strictly periodic phenomenon. An attempt to find out a period from all available records of the H-beta emission strength led to a value of 8044+-167 days (22.0 years) but the phase plots show that each cycle has a different shape and length. For the moment, the nature and origin of the disk around kappa Dra still remains unknown. From the comparison of the electronic spectra obtained at different phases of the long-term cycle and synthetic spectra it appears that there are no detectable changes in the photospheric part of the Balmer lines related to variations in the Balmer emission strength which could be attributed to an extended photosphere corresponding to inner parts of the disk, optically thick in continuum.

An analysis of accurate radial velocities (RVs) of the Be star gamma Cas from 295 Reticon spectrograms secured between October 1993 and May 2000 allowed us to prewhiten the RVs for the long-term changes and to obtain the first orbital RV curve of this star. The orbital period is 203d59 and the orbit has an eccentricity of 0.26. The orbital motion is detectable even in the published velocities, based on photographic spectra. This implies that gamma Cas is a primary component of a spectroscopic binary. The secondary has a mass of about 1 solar mass, appropriate for a white dwarf or a neutron star, but it could also be a normal late-type dwarf. The ultimate solution of the dispute whether the observed X-ray emission is associated with the secondary or with the primary will need further dedicated studies.