In this proceeding, we present a short review of the fascinating nebulosities of symbiotic binary R Aquarii. The R Aquarii system, comprising the central binary and surrounding nebular material, has been the subject of near-continuous study since its discovery, with a few hundred papers listed in ADS. As such, it is impossible to provide here the comprehensive review that R Aquarii deserves, instead we chose to focus on the nebulosities -- covering both our own research and other relevant results from the literature.

Yellow hypergiants are evolved massive stars that were suggested to be in post-red supergiant stage. Post-red supergiants that evolve back to the blue, hot side of the Hertzsprung-Russell diagram can intersect a temperature domain in which their atmospheres become unstable against pulsations (the Yellow Void or Yellow Wall), and the stars can experience outbursts with short, but violent mass eruptions. The yellow hypergiant Rho Cas is famous for its historical and recent outbursts, during which the star develops a cool, optically thick wind with a very brief but high mass-loss rate, causing a sudden drop in the light curve. Here we report on a new outburst of Rho Cas which occurred in 2013, accompanied by a temperature decrease of ~3000 K and a brightness drop of 0.6 mag. During the outburst TiO bands appear, together with many low excitation metallic atmospheric lines characteristic for a later spectral type. With this new outburst, it appears that the time interval between individual events decreases, which might indicate that Rho Cas is preparing for a major eruption that could help the star to pass through the Yellow Void. We also analysed the emission features that appear during phases of maximum brightness and find that they vary synchronous with the emission in the prominent [CaII] lines. We conclude that the occasionally detected emission in the spectra of Rho Cas, as well as certain asymmetries seen in the absorption lines of low to medium-excitation potential, are circumstellar in nature, and we discuss the possible origin of this material.

R Aquarii is a symbiotic binary surrounded by a large and complex nebula with a prominent curved jet. It is one of the closest known symbiotic systems, and therefore offers a unique opportunity to study the central regions of these systems and the formation and evolution of astrophysical jets. We studied the evolution of the central jet and outer nebula of R Aqr taking advantage of a long term monitoring campaign of optical imaging, as well as of high-resolution integral field spectroscopy. Narrow-band images acquired over a period of more than 21 years are compared in order to study the expansion and evolution of all components of the R Aqr nebula. The magnification method is used to derive the kinematic ages of the features that appear to expand radially. Integral field spectroscopy of the OIII 5007A emission is used to study the velocity structure of the central regions of the jet. New extended features, further out than the previously known hourglass nebula, are detected. The kinematic distance to R Aqr is calculated to be 178 pc using the expansion of the large hourglass nebula. This nebula of R Aqr is found to be roughly 650 years old, while the inner regions have ages ranging from 125 to 290 years. The outer nebula is found to be well described by a ballistic expansion, while for most components of the jet strong deviations from such behaviour are found. We find that the Northern jet is mostly red-shifted while its Southern part is blue-shifted, apparently at odds with findings from previous studies but almost certainly a consequence of the complex nature of the jet and variations in ionisation and illumination between observations.

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.

B[e] supergiants and yellow hypergiants share a number of common properties regarding their circumstellar environments. Using the forbidden [O I] and [Ca II] lines as disk tracers, we suggest the presence of a Keplerian disk or ring around the yellow hypergiant V509 Cas and confirm the pole-on inner disk around V1302 Aql. These findings indicate a change in mass-loss behavior from spherical in cooler yellow hypergiants to axisymmetric in the hotter ones during the passage through the Yellow Void. The accumulation of material in the equatorial plane reminds of the disks of B[e] supergiants, supporting the suggestion that yellow hypergiants might appear as B[e] supergiants after they reach the blue edge of the yellow instability domain.

Spectroscopic monitoring of the yellow hypergiant $\rho$ Cas revealed a new outburst in 2013, which is obvious from the development of TiO bands in the spectra. Also many atmospheric lines characteristic for a later spectral type appear. This spectroscopic outburst is in agreement with the photometric light curve, which displays a drop by about 0.6 mag during the same period.

Yellow Hypergiants (YHGs) and B[e] supergiants (B[e]SGs), though in different phases in their evolution, display many features in common. This is partly due to the fact that both types of objects undergo strong, often asymmetric mass loss, and the ejected material accumulates in shells, rings, or disk-like structures, giving rise to emission from warm molecules and dust. We performed an optical spectroscopic survey of northern Galactic emission-line stars aimed at identifying tracers for the structure and kinematics of circumstellar environments. We identified two sets of lines, [O I] and [Ca II], which originate from the discs of B[e]SGs. The same set of lines is observed in V1302 Aql and V509 Cas, which are both hot YHGs. While V1302 Aql is known to have a disc-like structure, the kinematical broadening of the lines in V509 Cas suggest a Keplerian disk or ring around this star alike their hotter B[e]SG counterparts.

The stellar evolution theory of massive stars remains uncalibrated with high-precision photometric observational data mainly due to a small number of luminous stars that are monitored from space. Automated all-sky surveys have revealed numerous variable stars but most of the luminous stars are often overexposed. Targeted campaigns can improve the time base of photometric data for those objects. The aim of this investigation is to study the variability of luminous stars at different timescales in young open clusters and OB associations. We monitored 22 open clusters and associations from 2011 to 2013 using a 0.25-m telescope. Variable stars were detected by comparing the overall light-curve scatter with measurement uncertainties. Variability was analysed by the light curve feature extraction tool FATS. Periods of pulsating stars were determined using the discrete Fourier transform code SigSpec. We then classified the variable stars based on their pulsation periods and available spectral information. We obtained light curves for more than 20000 sources of which 354 were found to be variable. Amongst them we find 80 eclipsing binaries, 31 $\alpha$ Cyg, 13 $\beta$ Cep, 62 Be, 16 slowly pulsating B, 7 Cepheid, 1 $\gamma$ Doradus, 3 Wolf-Rayet and 63 late-type variable stars. Up to 55% of these stars are potential new discoveries as they are not present in the Variable Star Index (VSX) database. We find the cluster membership fraction for variable stars to be 13% with an upper limit of 35%.

Context. Mass loss is an important property in evolution models of massive stars. As up to 90% of the massive stars have a visual or spectroscopic companion and many of them exhibit mass exchange, mass-loss rates can be acquired through the period study of massive binaries. Aims. Using our own photometric observations as well as archival data, we look for variations in orbital periods of seven massive eclipsing binary systems in the Cygnus OB2 association and estimate their mass-loss rates and stellar parameters. Methods. We use a Bayesian parameter estimation method to simultaneously fit the period and period change to all available data and a stellar modelling tool to model the binary parameters from photometric and radial-velocity data. Results. Four out of the seven selected binaries show non-zero period change values at two-sigma confidence level. We also report for the first time the eclipsing nature of a star MT059.

Due to their large expansion speed, the apparent growth of nearby nova remnants such as GK Per can be easily resolved from ground-based optical imagery on a timescale of months. If the expansion in the plane of the sky is coupled with the Doppler shift velocities, an almost complete dynamical picture is drawn. We will discuss our latest results of such a study on GK Per.

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=15-17) 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.

We present a kinematical study of the optical ejecta of GK Per. It is based on proper motions measurements of 282 knots from ~20 images spanning 25 years. Doppler-shifts are also computed for 217 knots. The combination of proper motions and radial velocities allows a unique 3-D view of the ejecta to be obtained. The main results are: (1) the outflow is a thick shell in which knots expand with a significant range of velocities, mostly between 600 and 1000 km/s; (2) kinematical ages indicate that knots have suffered only a modest deceleration since their ejection a century ago; (3) no evidence for anisotropy in the expansion rate is found; (4) velocity vectors are generally aligned along the radial direction but a symmetric pattern of non-radial velocities is also observed at specific directions; (5) the total Halpha+[NII] flux has been linearly decreasing at a rate of 2.6 % per year in the last decade. The Eastern nebular side is fading at a slower rate than the Western one. Some of the knots displayed a rapid change of brightness during the 2004-2011 period. Over a longer timescale, a progressive circularization and homogenization of the nebula is taking place; (6) a kinematic distance of 400+-30 pc is determined. These results raise some problems to the previous interpretations of the evolution of GK Per. In particular, the idea of a strong interaction of the outflow with the surrounding medium in the Southwest quadrant is not supported by our data.

The Gaia satellite will be launched at the end of 2011. It will observe at least 1 billion stars, and among them several million emission line stars and hot stars. Gaia will provide parallaxes for each star and spectra for stars till V magnitude equal to 17. After a general description of Gaia, we present the codes and methods, which are currently developed by our team. They will provide automatically the astrophysical parameters and spectral classification for the hot and emission line stars in the Milky Way and other close Local Group galaxies such as the Magellanic Clouds.

We provide a quantitative analysis of time-variable phenomena in the photospheric, near-star, and outflow regions of the late-B supergiant (SG) HD 199478. The analysis is based primarily on optical spectroscopic datasets secured between 1999 and 2000 from the Bulgarian NAO, Tartu, and Ritter Observatories. The temporal behaviour of HD 199478 is characterised by three key empirical properties: (i) systematic central velocity shifts in the photospheric absorption lines, including C II and He I, over a characteristic time-scale of abou 20 days; (ii) extremely strong, variable H alpha emission with no clear modulation signal, and (iii) the occurrence in 2000 of a (rare) high-velocity absorption (HVA) event in H alpha, which evolved over about 60 days, showing the clear signature of mass infall and outflows. In these properties HD 199478 resembles few other late-B SGs with peculiar emission and HVAs in H alpha (HD 91619, HD 34085, HD 96919). Non-LTE line synthesis modelling is conducted using FASTWIND for these late-B SGs to constrain and compare their fundamental parameters within the context of extreme behaviour in the H alpha lines. Our analysis indicate that at the cooler temperature edge of B SGs, there are objects whose wind properties, as traced by H alpha, are inconsistent with the predictions of the smooth, spherically symmetric wind approximation. This discordance is still not fully understood and may highlight the role of a non-spherical, disk-like, geometry, which may result from magnetically-driven equatorial compression of the gas. Ordered dipole magnetic fields may also lead to confined plasma held above the stellar surface, which ultimately gives rise to transient HVA events.

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 (O6.5f?pe-O8fp). The Balmer and HeI lines show strongly variable emission which is highly reproducible on a well-determined 538-d period. Metal lines and HeII absorptions (including many selective emission lines but excluding He II 4686A emission) are essentially constant in line strength, but are variable in velocity, establishing a double-lined binary orbit with P(orb) = 1542d, e=0.45. We conduct a model-atmosphere analysis of the primary, and find that the system is consistent with a O8: giant with a B1: main-sequence secondary. Since the periodic 538-d changes are unrelated to orbital motion, rotational modulation of a magnetically constrained plasma is strongly favoured as the most likely underlying `clock'. An upper limit on the equatorial rotation is consistent with this hypothesis, but is too weak to provide a strong constraint.

The European Gaia astrometry mission is due for launch in 2011. Gaia will rely on the proven principles of ESA's Hipparcos mission to create an all-sky survey of about one billion stars throughout our Galaxy and beyond, by observing all objects down to 20th magnitude. Through its massive measurement of stellar distances, motions and multi-colour photometry it will provide fundamental data necessary for unravelling the structure, formation and evolution of the Galaxy. This paper presents the design and performance of the broad- and medium-band set of photometric filters adopted as the baseline for Gaia. The nineteen selected passbands (extending from the ultraviolet to the far-red), the criteria, and the methodology on which this choice has been based are discussed in detail. We analyse the photometric capabilities for characterizing the luminosity, temperature, gravity and chemical composition of stars. We also discuss the automatic determination of these physical parameters for the large number of observations involved, for objects located throughout the entire Hertzsprung-Russell diagram. Finally, the capability of the photometric system to deal with the main Gaia science case is outlined.