We aim to explore the variable photometric and stellar properties of four yellow hypergiants (YHGs), HR8752, HR 5171A, $\rho$ Cas, and HD 179821, and their pulsations of hundreds of days, and long-term variations (LTVs) of years. We tackled multi-colour and visual photometric data sets, looked for photometric indications betraying eruptions or enhanced mass-loss episodes, calculated stellar properties mainly using a published temperature calibration, and investigated the nature of LTVs and their influence on quasi-periods and stellar properties. The $BV$ photometry revealed a high-opacity layer in the atmospheres. When the temperature rises the mass loss increases as well, consequently, as the density of the high-opacity layer. As a result, the absorption in $B$ and $V$ grow. The absorption in $B$, presumably of the order of one to a few 0\fm1, is always higher than in $V$. This difference renders redder and variable $(B-V)$ colour indexes, but the absorption law is unknown. This property of YHGs is unpredictable and explains why spectroscopic temperatures are always higher than photometric ones. We propose shorter distances for $\rho$ Cas and HR 5171A than the accepted ones. Therefore, a correction to decrease the blue luminescence of HR 5171A by polycyclic aromatic hydrocarbon (PAH) molecules is necessary, and HR 5171A would no longer be a member of the cluster Gum48d. HR 5171A is only subject to one source of light variation, not by two as the literature suggests. Eruptive episodes of YHGs prefer relatively cool circumstances when a red evolutionary loop (RL) has shifted the star to the red on the HR diagram. After the eruption, a blue loop evolution (BL) is triggered lasting one to a few decades. The reddening episode of HR 5171A between 1960 and 1974 was most likely due to a red loop evolution, and the reddening after the 1975 eruption was likely due to a shell ejection.

Fundamental atomic transition parameters, such as oscillator strengths and wavelengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work under way to produce these parameters for many ions, uncertainties in these parameters remain large and can limit the accuracy of chemical abundance determinations. The Belgian repository of fundamental atomic data and stellar spectra (BRASS) aims to provide a large systematic and homogeneous quality assessment of the atomic data available for quantitative spectroscopy. BRASS shall compare synthetic spectra against extremely high quality observed spectra, at a resolution of ~85000 and signal-noise ratios of ~1000, for around 20 bright BAFGK spectral type stars, in order to evaluate the atomic data available for over a thousand potentially useful spectral lines. A large-scale homogeneous selection of atomic lines is performed by synthesising theoretical spectra of literature atomic lines, for FGK-type stars including the Sun, resulting in a selection of 1091 theoretically deep and unblended lines, in the wavelength range 4200-6800~\AA, which may be suitable for quality assessment. Astrophysical log(gf) values are determined for the 1091 transitions using two commonly employed methods. The agreement of these log(gf) values are used to select well-behaving lines for quality assessment. 845 atomic lines were found to be suitable for quality assessment, of which 408 were found to be robust against systematic differences between analysis methods. Around 54% of the quality-assessed lines were found to have at least one literature log(gf) value in agreement with our derived values, though the remaining values can disagree by as much as 0.5 dex. Only 38% of FeI lines were found to have sufficiently accurate log(gf) values, increasing to ~70-75% for the remaining Fe-group lines.

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.

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.

Fundamental atomic parameters, such as oscillator strengths, play a key role in modelling and understanding the chemical composition of stars in the universe. Despite the significant work underway to produce these parameters for many astrophysically important ions, uncertainties in these parameters remain large and can propagate throughout the entire field of astronomy. The Belgian repository of fundamental atomic data and stellar spectra (BRASS) aims to provide the largest systematic and homogeneous quality assessment of atomic data to date in terms of wavelength, atomic and stellar parameter coverage. To prepare for it, we first compiled multiple literature occurrences of many individual atomic transitions, from several atomic databases of astrophysical interest, and assessed their agreement. Several atomic repositories were searched and their data retrieved and formatted in a consistent manner. Data entries from all repositories were cross-matched against our initial BRASS atomic line list to find multiple occurrences of the same transition. Where possible we used a non-parametric cross-match depending only on electronic configurations and total angular momentum values. We also checked for duplicate entries of the same physical transition, within each retrieved repository, using the non-parametric cross-match. We report the cross-matched transitions for each repository and compare their fundamental atomic parameters. We find differences in log(gf) values of up to 2 dex or more. We also find and report that ~2% of our line list and Vienna Atomic Line Database retrievals are composed of duplicate transitions. Finally we provide a number of examples of atomic spectral lines with different log(gf) values, and discuss the impact of these uncertain log(gf) values on quantitative spectroscopy. All cross-matched atomic data and duplicate transitions are available to download at brass.sdf.org.

The long-term X-ray lightcurves of classical supergiant X-ray binaries and supergiant fast X-ray transients show relatively similar super-orbital modulations, which are still lacking a sound interpretation. We propose that these modulations are related to the presence of corotating interaction regions (CIRs) known to thread the winds of OB supergiants. To test this hypothesis, we couple the outcomes of 3-D hydrodynamic models for the formation of CIRs in stellar winds with a simplified recipe for the accretion onto a neutron star. The results show that the synthetic X-ray light curves are indeed modulated by the presence of the CIRs. The exact period and amplitude of these modulations depend on a number of parameters governing the hydrodynamic wind models and on the binary orbital configuration. To compare our model predictions with the observations, we apply the 3-D wind structure previously shown to well explain the appearance of discrete absorption components in the UV time series of a prototypical B0.5I-type supergiant. Using the orbital parameters of IGRJ16493-4348 which has the same B0.5I donor spectral type, the period and modulations in the simulated X-ray light-curve are similar to the observed ones, thus providing support to our scenario. We propose, that the presence of CIRs in donor star winds should be considered in future theoretical and simulation efforts of wind-fed X-ray binaries.

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.

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 medium-resolution 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 ground-based 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 mid-resolution spectra obtained with the GIRAFFE and UVES spectrographs.

Context. Luminous blue variables (LBVs) are rare massive stars with very high luminosity. They are characterized by strong photo-metric and spectroscopic variability related to transient eruptions. The mechanisms at the origin of these eruptions is not well known. In addition, their formation is still problematic and the presence of a companion could help to explain how they form. Aims. This article presents a study of seven LBVs (about 20% of the known Galactic population), some Wolf-Rayet stars, and massive binaries. We probe the environments that surround these massive stars with near-, mid-, and far-infrared images, investigating potential nebula/shells and the companion stars. Methods. To investigate large spatial scales, we used seeing-limited and near diffraction-limited adaptive optics images to obtain a differential diagnostic on the presence of circumstellar matter and to determine their extent. From those images, we also looked for the presence of binary companions on a wide orbit. Once a companion was detected, its gravitational binding to the central star was tested. Tests include the chance projection probability, the proper motion estimates with multi-epoch observations, flux ratio, and star separations. Results. We find that two out of seven of LBVs may have a wide orbit companion. Most of the LBVs display a large circumstellar envelope or several shells. In particular, HD168625, known for its rings, possesses several shells with possibly a large cold shell at the edge of which the rings are formed. For the first time, we have directly imaged the companion of LBV stars.

We re-examined photometry (VBLUW, UBV, uvby) of the yellow hypergiant HR 5171A made a few decades ago. In that study no proper explanation could be given for the enigmatic brightness excesses in the L band (VBLUW system, lambda_eff=3838 A). In the present paper, we suggest that this might have been caused by blue luminescence (BL), an emission feature of neutral polycyclic aromatic hydrocarbon molecules (PAHs), discovered in 2004. It is a fact that the highest emission peaks of the BL lie in the L band. Our goals were to investigate other possible causes, and to derive the fluxes of the emission. We used two-colour diagrams based on atmosphere models, spectral energy distributions, and different extinctions and extinction laws, depending on the location of the supposed BL source: either in Gum48d on the background or in the envelope of HR 5171A. False L-excess sources, such as a hot companion, a nearby star, or some instrumental effect, could be excluded. Also, emission features from a hot chromosphere are not plausible. The fluxes of the L excess, recorded in the data sets of 1971, 1973, and 1977 varied (all in units of 10^(-10) W m^(-2) micron^(-1)) between 1.4 to 21, depending on the location of the source. A flux near the low side of this range is preferred. Small brightness excesses in uv (uvby system) were present in 1979, but its connection with BL is doubtful. For the L fluxes we consider the lowest values as more realistic. The uncertainties are 20-30 %. Similar to other yellow hypergiants, HR 5171A showed powerful brightness outbursts, particularly in the 1970s. A release of stored H-ionization energy by atmospheric instabilities could create BL emitted by neutral PAHs.

We investigate the physical properties of large-scale wind structures around massive hot stars with radiatively-driven winds. We observe Discrete Absorption Components (DACs) in optical He I P Cygni lines of the LBV binary MWC 314 (Porb=60.8 d). The DACs are observed during orbital phases when the primary is in front of the secondary star. They appear at wind velocities between -100 km/s and -600 km/s in the P Cyg profiles of He I lam5875, lam6678, and lam4471, signaling high-temperature expanding wind regions of enhanced density and variable outflow velocity. The DACs can result from wave propagation linked to the orbital motion near the low-velocity wind base. The He I lines indicate DAC formation close to the primary's surface in high-temperature wind regions in front of its orbit, or in dynamical wind regions confined between the binary stars. We observed the DACs with Mercator-HERMES on 5 Sep 2009, 5 May 2012, and 6 May 2014 when the primary is in front of the secondary star. XMM-Newton observations of 6 May 2014 significantly detected MWC 314 in X-rays at an average rate of ~0.015 cts/s.

This work describes a procedure to derive several fundamental parameters such as the effective temperature, surface gravity, equatorial rotational velocity and microturbulent velocity. In this work, we have written a numerical procedure in Python which finds the best fit between a grid of synthetic spectra and the observed spectra by minimizing a standard chi-square. LTE model atmospheres were calculated using the ATLAS9 code and were used as inputs to the spectrum synthesis code SYNSPEC48 in order to compute a large grid of synthetic Balmer line profiles. This new procedure has been applied to a large number of new observations (GIRAFFE spectra) of B and A stars members of the young open cluster NGC3293. These observations are part of the GAIA ESO Survey. Takeda's procedure was also used to derive rotational velocities and microturbulent velocities. The results have been compared to previous determinations by other authors and are found to agree with them. As a first result, we concluded that using this procedure, an accuracy of +-200 K could be achieved in effective temperature and +-0.2 dex in surface gravities.

Context : We still do not know which mechanisms are responsible for the transport of angular momentum inside stars. The recent detection of mixed modes that contain the signature of rotation in the spectra of Kepler subgiants and red giants gives us the opportunity to make progress on this issue. Aims: Our aim is to probe the radial dependance of the rotation profiles for a sample of Kepler targets. For this purpose, subgiants and early red giants are particularly interesting targets because their rotational splittings are more sensitive to the rotation outside the deeper core than is the case for their more evolved counterparts. Methods: We first extract the rotational splittings and frequencies of the modes for six young Kepler red giants. We then perform a seismic modeling of these stars using the evolutionary codes CESAM2k and ASTEC. By using the observed splittings and the rotational kernels of the optimal models, we perform inversions of the internal rotation profiles of the six stars. Results: We obtain estimates of the mean rotation rate in the core and in the convective envelope of these stars. We show that the rotation contrast between the core and the envelope increases during the subgiant branch. Our results also suggest that the core of subgiants spins up with time, contrary to the RGB stars whose core has been shown to spin down. For two of the stars, we show that a discontinuous rotation profile with a deep discontinuity reproduces the observed splittings significantly better than a smooth rotation profile. Interestingly, the depths that are found most probable for the discontinuities roughly coincide with the location of the H-burning shell, which separates the layers that contract from those that expand. These results will bring observational constraints to the scenarios of angular momentum transport in stars.

We initiated long-term optical interferometry monitoring of the diameters of unstable yellow hypergiants (YHG) with the goal of detecting both the long-term evolution of their radius and shorter term formation related to large mass-loss events. We observed HR5171 A with AMBER/VLTI. We also examined archival photometric data in the visual and near-IR spanning more than 60 years, as well as sparse spectroscopic data. HR5171A exhibits a complex appearance. Our AMBER data reveal a surprisingly large star for a YHG R*=1315+/-260Rsun\ (~6.1AU) at the distance of 3.6+/-0.5kpc. The source is surrounded by an extended nebulosity, and these data also show a large level of asymmetry in the brightness distribution of the system, which we attribute to a newly discovered companion star located in front of the primary star. The companion's signature is also detected in the visual photometry, which indicates an orbital period of Porb=1304+/-6d. Modeling the light curve with the NIGHTFALL program provides clear evidence that the system is a contact or possibly over-contact eclipsing binary. A total current system mass of 39^{+40}_{-22} solar mass and a high mass ratio q>10 is inferred for the system. The low-mass companion of HR5171 A is very close to the primary star that is embedded within its dense wind. Tight constraints on the inclination and vsini of the primary are lacking, which prevents us from determining its influence precisely on the mass-loss phenomenon, but the system is probably experiencing a wind Roche-Lobe overflow. Depending on the amount of angular momentum that can be transferred to the stellar envelope, HR5171 A may become a fast-rotating B[e]/Luminous Blue Variable (LBV)/Wolf-Rayet star. In any case, HR5171 A highlights the possible importance of binaries for interpreting the unstable YHGs and for massive star evolution in general.

We carried out an extensive observational study of the Slowly Pulsating B (SPB) star, HD 25558. The ~2000 spectra obtained at different observatories, the ground-based and MOST satellite light curves revealed that this object is a double-lined spectroscopic binary with an orbital period of about 9 years. The observations do not allow the inference of an orbital solution. We determined the physical parameters of the components, and found that both lie within the SPB instability strip. Accordingly, both show line-profile variations due to stellar pulsations. Eleven independent frequencies were identified in the data. All the frequencies were attributed to one of the two components based on Pixel-by-pixel variability analysis of the line profiles. Spectroscopic and photometric mode identification was also performed for the frequencies of both stars. These results suggest that the inclination and rotation of the two components are rather different. The primary is a slow rotator with ~6 d period, seen at ~60 deg inclination, while the secondary rotates fast with ~1.2 d period, and is seen at ~20 inclination. Spectropolarimetric measurements revealed that the secondary component has a magnetic field with at least a few hundred Gauss strength, while no magnetic field can be detected in the primary.

The unparalleled photometric data obtained by NASA's Kepler space telescope led to an improved understanding of red giant stars and binary stars. Seismology allows us to constrain the properties of red giants. In addition to eclipsing binaries, eccentric non-eclipsing binaries, exhibiting ellipsoidal modulations, have been detected with Kepler. We aim to study the properties of eccentric binary systems containing a red giant star and derive the parameters of the primary giant component. We apply asteroseismic techniques to determine masses and radii of the primary component of each system. For a selected target, light and radial velocity curve modelling techniques are applied to extract the parameters of the system. The effects of stellar on the binary system are studied. The paper presents the asteroseismic analysis of 18 pulsating red giants in eccentric binary systems, for which masses and radii were constrained. The orbital periods of these systems range from 20 to 440days. From radial velocity measurements we find eccentricities between e=0.2 to 0.76. As a case study we present a detailed analysis of KIC5006817. From seismology we constrain the rotational period of the envelope to be at least 165 d, roughly twice the orbital period. The stellar core rotates 13 times faster than the surface. From the spectrum and radial velocities we expect that the Doppler beaming signal should have a maximum amplitude of 300ppm in the light curve. Through binary modelling, we determine the mass of the secondary component to be 0.29$\pm$0.03\,$M_\odot$. For KIC5006817 we exclude pseudo-synchronous rotation of the red giant with the orbit. The comparison of the results from seismology and modelling of the light curve shows a possible alignment of the rotational and orbital axis at the 2$\sigma$ level. Red giant eccentric systems could be progenitors of cataclysmic variables and hot subdwarf B stars.

A programme to observe all A dwarfs in open clusters brighter than V=6.5 mag of various ages and in the field was initiated several years ago. In this work we present the current status of microturbulent velocity for A and F dwarfs. We have performed high resolution high signal-to-noise spectroscopy of stars well distributed in mass along the Main Sequence. Microturbulent velocities are derived iteratively by fitting grids of synthetic spectra calculated in LTE to observed spectra of 61 A field stars, 55 A and 58 F in open clusters (Pleiades, Coma Berenices, Hyades and the Ursa Major moving group). We compared our results to recent works and found a good agreement with their analytical formulation for the standard microturbulence. Our results show a broad maximum for microturbulent velocities in the range A5V to about A9V and a decrease (to $\sim 1$ km/s) for cooler and hotter stars as indicated in Smalley (2004).We also present a comparison to preliminary science results of Lobel et al. (2013) for the Gaia-ESO Public Spectroscopic Survey.

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 spectroscopic analysis of MWC 314, a luminous blue variable (LBV) candidate with an extended bipolar nebula. The detailed spectroscopic variability is investigated to determine if MWC 314 is a massive binary system with a supersonically accelerating wind or a low-mass B[e] star. We compare the spectrum and spectral energy distribution to other LBVs (such as P Cyg) and find very similar physical wind properties, indicating strong kinship. We combined long-term high-resolution optical spectroscopic monitoring and V-band photometric observations to determine the orbital elements and stellar parameters and to investigate the spectral variability with the orbital phases. We developed an advanced model of the large-scale wind-velocity and wind-density structure with 3-D radiative transfer calculations that fit the orbitally modulated P Cyg profile of He I lam5876, showing outflow velocities above 1000 km/s. We find that MWC 314 is a massive semi-detached binary system of ~1.22 AU, observed at an inclination angle of i=72.8 deg. with an orbital period of 60.8 d and e=0.23. The primary star is a low-vsini LBV candidate of m1=39.6 Msun and R1=86.8 Rsun. The detailed radiative transfer fits show that the geometry of wind density is asymmetric around the primary star with increased wind density by a factor of 3.3, which leads the orbit of the primary. The variable orientation causes the orbital modulation that is observed in absorption portions of P Cyg wind lines. Wind accretion in the system produces a circumbinary disc. MWC 314 is in a crucial evolutionary phase of close binary systems, when the massive primary star has its H envelope being stripped and is losing mass to a circumbinary disc. MWC 314 is a key system for studying the evolutionary consequences of these effects.

Spectroscopic monitoring with Mercator-HERMES over the past two years reveals that MWC 314 is a massive binary system composed of an early B-type primary LBV star and a less-luminous supergiant companion. We determine an orbital period Porb of 60.85 d from optical S II and Ne I absorption lines observed in this single-lined spectroscopic binary. We find an orbital eccentricity of e=0.26, and a large amplitude of the radial velocity curve of 80.6 km/s. The ASAS V light-curve during our spectroscopic monitoring reveals two brightness minima (\Delta V~0.1 mag.) over the orbital period due to partial eclipses at an orbital inclination angle of ~70 degrees. We find a clear correlation between the orbital phases and the detailed shapes of optical and near-IR P Cygni-type line profiles of He I, Si II, and double- or triple-peaked stationary cores of prominent Fe II emission lines. A preliminary 3-D radiative transfer model computed with Wind3D shows that the periodic P Cygni line profile variability results from an asymmetric common-envelope wind with enhanced density (or line opacity) in the vicinity of the LBV primary. The variable orientation of the inner LBV wind region due to the orbital motion produces variable P Cygni line profiles (with wind velocities of ~200 km/s) between orbital phases \phi = 0.65 to 0.85, while weak inverse P Cygni profiles are observed half an orbital period later around \phi = 0.15 to 0.35. We do not observe optical or near-IR He II, C III, and Si III lines, signaling that the LBV's spectral type is later than B0. Detailed modeling of the asymmetrical wind properties of massive binary MWC 314 provides important new physical information about the most luminous hot (binary) stars such as Eta Carinae.

Ultraviolet and visible spectra of the symbiotic nova RR Telescopii are used to derive reference wavelengths for many forbidden and intercombination transitions of ions +1 to +6 of elements C, N, O, Ne, Na, Mg, Al, Si, P, S, Cl, Ar, K and Ca. The wavelengths are then used to determine new energy values for the levels within the ions' ground configurations or first excited configuration. The spectra were measured by the Space Telescope Imaging Spectrograph of the Hubble Space Telescope and the UltraViolet Echelle Spectrograph of the European Southern Observatory in 2000 and 1999, respectively, and cover 1140 to 6915 A. Particular care was taken to assess the accuracy of the wavelength scale between the two instruments. An investigation of the profiles of the emission lines reveals that the nebula consists of at least two plasma components at different velocities. The components have different densities, and a simple model of the lines' emissions demonstrates that most of the lines principally arise from the high density component. Only these lines were used for the wavelength study.

We develop 3-D models of the structured winds of massive hot stars with the Wind3D radiative transfer (RT) code. We investigate the physical properties of large-scale structures observed in the wind of the B-type supergiant HD 64760 with detailed line profile fits to Discrete Absorption Components (DACs) and rotational modulations observed with IUE in Si IV {\lambda}1395. We develop parameterized input models Wind3D with large-scale equatorial wind density- and velocity-structures, or so-called `Co-rotating Interaction Regions' (CIRs) and `Rotational Modulation Regions' (RMRs). The parameterized models offer important advantages for high-performance RT calculations over ab-initio hydrodynamic input models. The acceleration of the input model calculations permits us to simulate and investigate a wide variety of physical conditions in the extended winds of massive hot stars. The new modeling method is very flexible for constraining the dynamic and geometric wind properties of RMRs in HD 64760. We compute that the modulations are produced by a regular pattern of radial density enhancements that protrude almost linearly into the equatorial wind. We find that the modulations are caused by narrow `spoke-like' wind regions. We present a hydrodynamic model showing that the linearly shaped radial wind pattern can be caused by mechanical wave action at the base of the stellar wind from the blue supergiant.

We present results of a long-term spectroscopic monitoring program (since mid 2009) of Luminous Blue Variables with the new HERMES echelle spectrograph on the 1.2 m Mercator telescope at La Palma (Spain). We investigate high-resolution (R=80,000) optical spectra of two LBVs, P Cyg and HD 168607, the LBV candidates MWC 930 and HD 168625, and the LBV binary MWC 314. In P Cyg we observe flux changes in the violet wings of the Balmer H{\alpha}, H{\beta}, and He I lines between May and Sep 2009. The changes around 200 km/s to 300 km/s are caused by variable opacity at the base of the supersonic wind from the blue supergiant. We observe in MWC 314 broad double-peaked metal emission lines with invariable radial velocities over time. On the other hand, we measure in the photospheric S II {\lambda}5647 absorption line, with lower excitation energy of ~14 eV, an increase of the heliocentric radial velocity centroid from 37 km/s to 70 km/s between 5 and 10 Sep 2009 (and 43 km/s on 6 Apr 2010). The increase of radial velocity of ~33 km/s in only 5 days can confirm the binary nature of this LBV close to the Eddington luminosity limit. A comparison with VLT-UVES and Keck-Hires spectra observed over the past 13 years reveals strong flux variability in the violet wing of the H{\alpha} emission line of HD 168625, and in the absorption portion of the H{\beta} line of HD 168607. In HD 168625 we observe H{\alpha} wind absorption at velocities exceeding 200 km/s which develops between Apr and June 2010.

We develop a new method to determine oscillator strength values of atomic absorption lines with state-of-the-art detailed spectral synthesis calculations of the optical spectrum of the Sun and of standard spectral reference stars. We update the log(gf)-values of 911 neutral lines observed in the KPNO-FTS flux spectrum of the Sun and high-resolution echelle spectra (R=80,000) of Procyon (F5 IV-V) and {\epsilon} Eri (K2 V) observed with large signal-to-noise (S/N) ratios of ~2,000 using the new Mercator-Hermes spectrograph at La Palma Observatory (Spain). We find for 483 Fe I, 85 Ni I, and 51 Si I absorption lines in the sample a systematic over-estimation of the literature log(gf)-values with central line depths below 15 %. We employ a curve-of-growth analysis technique to test the accuracy of the new oscillator strength values and compare calculated equivalent line widths to the Moore, Minnaert, & Houtgast atlas of the Sun. The online SpectroWeb database at http://spectra.freeshell.org interactively displays the observed and synthetic spectra and provides the new log(gf)-values together with important atomic line data. The graphical database is under development for stellar reference spectra of every spectral sub-class observed with large spectral resolution and S/N ratios.

First results of near-IR adaptive optics (AO)-assisted imaging, interferometry, and spectroscopy of this Luminous Blue Variable (LBV) are presented. They suggest that the Pistol Star is at least double. If the association is physical, it would reinforce questions concerning the importance of multiplicity for the formation and evolution of extremely massive stars.

We present multi-instruments and multi-wavelengths observations of the famous LBV star Pistol Star. These observations are part of a larger program about early O stars at different metallicities. The Pistol star has been claimed as the most massive star known, with 250 solar masses. We present the preliminary results based on X-Shooter spectra, as well as the observations performed with the VLTI-AMBER and the VLT-NACO adaptive optics. The X-shooter spectrograph allows to obtain simultaneously a spectrum from the UV to the K-band with a resolving power of $\sim$15000. The preliminary results obtained indicate that Pistol Star has similar properties of Eta Car, including shells of matter, but also the binarity. Other objects of the program, here briefly presented, were selected for their particular nature: early O stars with mass discrepancies between stellar evolution models and observations, discrepancies with the wind momentum luminosity relation.

We present modeling research work of the winds and circumstellar environments of prototypical hot and cool massive stars using advanced radiative transfer (RT) calculations. This research aims at unraveling the detailed physics of various mass-loss mechanisms of luminous stars in the upper H-R diagram. Very recent 3-D RT calculations, combined with hydrodynamic simulations, show that radiatively-driven winds of OB supergiants are structured due to large-scale density- and velocity-fields caused by rotating bright spots. The mass-loss rates computed from matching DACs in HD 64760 (B Ib) do not reveal appreciable changes from the rates of smooth wind models. Intermediate yellow supergiants (such as Rho Cas, F-G Ia0), on the other hand, show prominent spectroscopic signatures of strongly increased mass-loss rates during episodic outbursts. Long-term spectroscopic monitoring of hypergiants near the Yellow Evolutionary Void reveals that their mass-loss rates and wind-structure are dominated by photospheric eruptions and large-amplitude pulsations that impart mechanical momentum to the circumstellar environment by propagating shock waves. In massive red supergiants, however, clear evidence for mechanical wave propagation from the sub-photospheric convection zones is lacking. Recent spatially resolved HST-STIS observations inside Betelgeuse's (M Iab) very extended chromosphere and dust envelope show evidence of warm chromospheric gas far beyond the dust condensation radius of RT models. Models for these long-term spectroscopic observations demonstrate that the chromospheric pulsations are not spherically symmetric. The STIS observations point to the importance of mechanical wave propagation for heating and sustaining chromospheric conditions in the extended winds of red supergiants.

We present detailed semi-empiric models for rotational modulations observed in ultraviolet wind lines of B0.5 supergiant HD 64760. We model the Rotational Modulation Regions (RMRs) with advanced 3-D radiative transfer calculations in the stellar wind and quantitatively fit the time-evolution of the Si IV 1395 resonance line. We find that the RMRs are due to linearly-shaped narrow sector-like density enhancements in the equatorial wind. Unlike the Co-rotating Interaction Regions (CIRs) which produce Discrete Absorption Components in the line, the RMRs do not spread out with larger distance above the stellar surface. The detailed best fit shows that the RMRs of HD 64760 have maximum density enhancements of ~17% above the surrounding smooth wind density, about twice smaller than the hydrodynamic models of CIRs that warp around the star. The semi-empiric 3-D transfer modelling reveals that the narrow spoke-like RMRs must co-exist with broader and curved large-scale CIR wind density structures in the equatorial plane of this fast rotating Ib-supergiant.

We model the detailed time-evolution of Discrete Absorption Components (DACs) observed in P Cygni profiles of the Si IV lam1400 resonance doublet lines of the fast-rotating supergiant HD 64760 (B0.5 Ib). We adopt the common assumption that the DACs are caused by Co-rotating Interaction Regions (CIRs) in the stellar wind. We perform 3D radiative transfer calculations with hydrodynamic models of the stellar wind that incorporate these large-scale density- and velocity-structures. We develop the 3D transfer code Wind3D to investigate the physical properties of CIRs with detailed fits to the DAC shape and morphology. The CIRs are caused by irregularities on the stellar surface that change the radiative force in the stellar wind. In our hydrodynamic model we approximate these irregularities by circular symmetric spots on the stellar surface. We use the Zeus3D code to model the stellar wind and the CIRs, limited to the equatorial plane. We constrain the properties of large-scale wind structures with detailed fits to DACs observed in HD 64760. A model with two spots of unequal brightness and size on opposite sides of the equator, with opening angles of 20 +/- 5 degr and 30 +/- 5 degr diameter, and that are 20 +/- 5 % and 8 +/- 5 % brighter than the stellar surface, respectively, provides the best fit to the observed DACs. The recurrence time of the DACs compared to the estimated rotational period corresponds to spot velocities that are 5 times slower than the rotational velocity. The mass-loss rate of the structured wind model for HD 64760 does not exceed the rate of the spherically symmetric smooth wind model by more than 1 %. The fact that DACs are observed in a large number of hot stars constrains the clumping that can be present in their winds, as substantial amounts of clumping would tend to destroy the CIRs.

We update the online SpectroWeb database of spectral standard reference stars with 1178 oscillator strength values of atomic absorption lines observed in the optical spectrum of the Sun and Procyon (Alpha CMi A). The updated line oscillator strengths are measured with best fits to the disk-integrated KPNO-FTS spectrum of the Sun observed between 4000 A and 6800 A using state-of-the-art detailed spectral synthesis calculations. A subset of 660 line oscillator strengths is validated with synthetic spectrum calculations of Procyon observed with ESO-UVES between 4700 A and 6800 A. The new log(gf)-values in SpectroWeb are improved over the values offered in the online Vienna Atomic Line Database (VALD). We find for neutral iron-group elements, such as Fe I, Ni I, Cr I, and Ti I, a statistically significant over-estimation of the VALD log(gf)-values for weak absorption lines with normalized central line depths below 15 %. For abundant lighter elements (e.g. Mg I and Ca I) this trend is statistically not significantly detectable, with the exception of Si I for which the log(gf)-values of 60 weak and medium-strong lines are substantially decreased to best fit the observed spectra. The newly measured log(gf)-values are available in the SpectroWeb database at http://spectra.freeshell.org which interactively displays the observed and computed stellar spectra, together with corresponding atomic line data.

SpectroWeb is an online maintained interactive graphical database of digital spectral atlases of spectral standard stars at http://spectra.freeshell.org . It is an efficient and user-friendly research tool for accurate analyses of stellar spectra observed with large spectral resolution, including the solar spectrum. The web-interface displays observed and theoretical stellar spectra, and comprehensively provides detailed atomic and molecular line information via user interaction. It fully integrates interactive spectrum visualization tools for the analysis, management, and maintenance of large volumes of spectral line-identification, -transition, and -property data. SpectroWeb 1.0 currently offers optical (3300-6800 A) flux normalized high-resolution spectra of Betelgeuse (M2), Arcturus (K1), The Sun (G2), Beta Aqr (G0), Procyon (F5), and Canopus (F0). The provided line identifications are based on state of the art spectrum synthesis calculations. The graphical database is under permanent development as an online repository of identified (absorption) lines in spectral standard reference stars, covering a broad range of stellar spectral types. Its object-oriented (Java) implementation offers future expansion capabilities to link and read stellar spectral atlases from various public internet sites.

We discuss the development of the new radiative transfer code Wind3D. It solves the non-LTE radiative transport problem in moving stellar atmosphere models in three geometric dimensions. The code accepts arbitrary 3D velocity fields in Cartesian geometry without assumptions of axial symmetry. Wind3D is currently implemented as a fully parallelized (exact) accelerated lambda iteration scheme with a two level atom formulation. The numerical transfer scheme is efficient and very accurate to trace small variations of local velocity gradients on line opacity in strongly scattering dominated extended stellar winds. We investigate the detailed formation of P Cygni line profiles observed in ultraviolet spectra of massive stars. We compute the detailed shape of these resonance lines to model local enhancements of line opacity that can for instance be caused by clumping in supersonically expanding winds. Wind3D will be applied to hydrodynamic models to investigate physical properties of discrete absorption line components.

We apply the 3-dimensional radiative transport code Wind3D to 3D hydrodynamic models of Corotating Interaction Regions to fit the detailed variability of Discrete Absorption Components observed in Si IV UV resonance lines of HD 64760 (B0.5 Ib). We discuss important effects of the hydrodynamic input parameters on these large-scale equatorial wind structures that determine the detailed morphology of the DACs computed with 3D transfer. The best fit model reveals that the CIR in HD 64760 is produced by a source at the base of the wind that lags behind the stellar surface rotation. The non-corotating coherent wind structure is an extended density wave produced by a local increase of only 0.6 % in the smooth symmetric wind mass-loss rate.

We report on multi-epoch optical and near-infrared spectroscopy around the first overtone ro-vibrational band of CO in the pulsating yellow hypergiant Rho Cas, one of the most massive stars in the Galaxy and a candidate SN II progenitor. We argue that the double cores of the CO absorption lines, that have previously been attributed to separate circumstellar shells expelled during its recurrent outbursts, result in fact from a superposition of a wide absorption line and a narrow central emission line. The CO line doubling returns over subsequent pulsation cycles, where the superposed line emission assumes its largest intensity near phases of maximum light. We find that the morphology and behavior of the CO band closely resemble the remarkable "line-splitting phenomenon" also observed in optical low-excitation atomic lines. Based on radiative transport calculations we present a simplified model of the near-IR CO emission emerging from cooler atmospheric layers in the immediate vicinity of the photosphere. We speculate that the kinetic temperature minimum in our model results from a periodical pulsation-driven shock wave. We further discuss a number of alternative explanations for the origin of the ubiquitous emission line spectrum, possibly due to a quasi-chromosphere or a steady shock wave at the interface of a fast expanding wind and the ISM. We present a number of interesting spectroscopic similarities between Rho Cas and other types of cool variable supergiants such as the RV Tau and R CrB stars. We further propose a possibly common mechanism for the enigmatic outburst behavior of these luminous pulsating cool stars.

We present an overview of radiative transfer modeling efforts to interpret spectra of a variety of stellar objects observed with FUSE. Detailed radiative transfer modeling of high ion emission line profiles of C III and O VI observed in the far-UV spectrum provides a powerful means to probe the thermal and dynamic properties of high-temperature plasmas in the atmospheres of stars. We model asymmetric emission lines of C III lam977 observed in spectra of luminous cool stars such as Alpha Aqr, to infer the wind- and microturbulence velocity structures of the upper chromosphere. Semi-empiric radiative transfer models that include transition region temperature conditions, are further developed based on detailed fits to O VI resonance emission lines in the supergiant Alpha Aqr, the classical Cepheid variable Beta Dor, and to self-absorbed O VI emission lines in the cataclysmic variable SW UMa. We observe that the C III resonance line profile of Alpha Aqr assumes a remarkable asymmetric shape, reminiscent of P Cygni type profiles observed in hot luminous supergiants. The model calculations indicate outflow velocities above ~140 km/s at kinetic temperatures of 65 kK and higher. Based on detailed model fits to the narrow red-shifted and self-absorbed O VI emission lines of SW UMa we compute that the gas- and electron-density exceed the density conditions of the upper solar transition region by about three orders of magnitude. We discuss how detailed semi-empiric fits to emission lines observed with the high spectral resolution of FUSE can provide reliable constraints on the mass-loss or mass-accretion rates in these objects.

The yellow hypergiant Rho Cassiopeiae (F-G Ia0) has recently become very active with a tremendous outburst event in the fall of 2000. During the event the pulsating supergiant dimmed by more than a visual magnitude, while its effective temperature decreased from 7000 K to below 4000 K over about 200 d, and we directly observed the largest mass-loss rate of about 5% of the solar mass in a single stellar outburst so far. Over the past three years since the eruption we observed a very prominent inverse P Cygni profile in Balmer H alpha, signaling a strong collapse of the upper atmosphere, also observed before the 2000 event. Continuous spectroscopic monitoring reveals that the H alpha line profile has transformed into a P Cygni profile since June 2003, presently (Sept 2004) signaling supersonic expansion velocities up to ~120 km/s in the extended upper atmosphere. Based on the very recent unique spectral evolution we observed the far-UV spectrum with the FUSE satellite in July 2004. The FUSE spectrum reveals that high-temperature plasma emission lines of O VI and C III are absent in the hypergiant, also observed for the red supergiant Alpha Ori (M2 Iab). On the other hand, we observe prominent transition region emission lines in the smaller (less luminous) classical Cepheid variable Beta Dor (F-G Iab-Ia), indicating that the mean atmospheric extension and surface gravity acceleration (as compared to effective temperature and atmospheric pulsation) play a major role for the formation of high-temperature stellar atmospheric plasmas. We present an overview of the recent spectral variability phases of Rho Cas with enhanced mass-loss from this enigmatic cool star.

The Hubble Space Telescope observed red supergiant Betelgeuse with the Space Telescope Imaging Spectrograph to investigate the outer atmosphere from spatially resolved spectra. We present a new set of seven high-resolution near-UV spectra observed with HST-STIS in fall 2002 and spring 2003, by scanning at chromospheric intensity peak-up position and six off-limb target positions up to three arcseconds away from the star. A small aperture is used to study and determine the thermal conditions and flow dynamics in its upper chromosphere and inner circumstellar dust envelope. We provide the first evidence for the presence of warm chromospheric plasma at least 3 arcseconds away from Betelgeuse at ~120 R* (1 R* ~= 700 Rsun) based on detailed spectroscopic observations of the Mg II h & k emission lines. Many other weak chromospheric emission lines as Fe II lam2716, C II lam2327, Al II] lam2669, and Fe I lam2823 are detected out to at least 1 arcsecond in the spatially resolved STIS observations. The profiles of the Mg II h & k and the Si I resonance emission lines reveal a strong increase of asymmetry by scanning off-limb, signaling the outward acceleration of wind expansion in Betelgeuse's upper chromosphere beyond 200 mas (~8 R*). We discuss detailed radiative transfer models that fit the STIS observations showing that the local kinetic gas temperature in the upper chromosphere exceeds 2600 K. Our radiation transport models for the IR silicate dust emission at 9.8 um in the upper chromosphere show however that the ambient gas temperature remains below 600 K to sustain the presence of dust grains. Hence, the STIS spectra of Betelgeuse's upper chromosphere directly demonstrate that warm chromospheric plasma must co-exist with cool dusty plasma in its outer atmosphere.

FUSE ultraviolet spectra of 8 giant and supergiant stars reveal that high temperature (3 X 10^5 K) atmospheres are common in luminous cool stars and extend across the color-magnitude diagram from Alpha Car (F0 II) to the cool giant Alpha Tau (K5 III). Emission present in these spectra includes chromospheric H-Lyman Beta, Fe II, C I, and transition region lines of C III, O VI, Si III, Si IV. Emission lines of Fe XVIII and Fe XIX signaling temperatures of ~10^7 K and coronal material are found in the most active stars, Beta Cet and 31 Com. A short-term flux variation, perhaps a flare, was detected in Beta Cet during our observation. Stellar surface fluxes of the emission of C III and O VI are correlated and decrease rapidly towards the cooler stars, reminiscent of the decay of magnetically-heated atmospheres. Profiles of the C III (977A) lines suggest that mass outflow is underway at T~80,000 K, and the winds are warm. Indications of outflow at higher temperatures (3 X 10^5K) are revealed by O VI asymmetries and the line widths themselves. High temperature species are absent in the M-supergiant Alpha Ori. Narrow fluorescent lines of Fe II appear in the spectra of many giants and supergiants, apparently pumped by H Lyman Alpha, and formed in extended atmospheres. Instrumental characteristics that affect cool star spectra are discussed.

Thirteen transitions within the ground 3d^2 configuration of Fe VII are identified in ultraviolet and optical spectra of the symbiotic star RR Telescopii obtained with the STIS instrument of the Hubble Space Telescope. The line fluxes are compared with theoretical data computed with the recent atomic data of K.A. Berrington et al., and high resolution optical spectra from VLT/UVES are used to identify blends. Seven branching ratios are measured, with three in good agreement with theory and one affected by blending. The lambda5277/lambda4943 branching ratio is discrepant by > 3 sigma, indicating errors in the atomic data for the lambda5277 line. A least-squares minimization scheme is used to simultaneously derive the temperature, T, and density, N_e, of the RR Tel nebula, and the interstellar extinction, E(B-V), towards RR Tel from the complete set of emission lines. The derived values are: log T/K = 4.50 +/- 0.23, log N_e/cm^-3=7.25 +/- 0.05, and E(B-V)<0.27. The extinction is not well-constrained by the Fe VII lines, but is consistent with the more accurate value E(B-V)=0.109^{+0.052}_{-0.059} derived here from the Ne V lambda2974/lambda1574 ratio in the STIS spectrum. Large differences between the K.A. Berrington et al. electron excitation data and the earlier F.P. Keenan & P.H. Norrington data-set are demonstrated, and the latter is shown to give worse agreement with observations.

We present spectroscopic observations of the short-period cataclysmic variable SW Ursa Majoris, obtained by the Far Ultraviolet Spectroscopic Explorer (FUSE) satellite while the system was in quiescence. The data include the resonance lines of O VI at 1031.91 and 1037.61 A. These lines are present in emission, and they exhibit both narrow (~ 150 km/s) and broad (~ 2000 km/s) components. The narrow O VI emission lines exhibit unusual double-peaked and redshifted profiles. We attribute the source of this emission to a cooling flow onto the surface of the white dwarf primary. The broad O VI emission most likely originates in a thin, photoionized surface layer on the accretion disk. We searched for emission from H_2 at 1050 and 1100 A, motivated by the expectation that the bulk of the quiescent accretion disk is in the form of cool, molecular gas. If H_2 is present, then our limits on the fluxes of the H_2 lines are consistent with the presence of a surface layer of atomic H that shields the interior of the disk. These results may indicate that accretion operates primarily in the surface layers of the disk in SW UMa. We also investigate the far-UV continuum of SW UMa and place an upper limit of 15,000 K on the effective temperature of the white dwarf.

We present spatially resolved spectra observed with HST-STIS of the upper chromosphere and dust envelope of Alpha Orionis (M2 Iab). In the fall of 2002 a set of five high-resolution near-UV spectra was obtained by scanning at intensity peak-up position and four off-limb target positions up to one arcsecond, using a small aperture, to investigate the thermal conditions and flow dynamics in the outer atmosphere of this important nearby cool supergiant star. Based on Mg II h & k, Fe II 2716 A, C II 2327 A, and Al II ] 2669 A emission lines we provide the first evidence for the presence of warm chromospheric plasma at least 1 arcsecond away from the star at ~40 R* (1 R*~700 Rsun). The STIS spectra reveal that Betelgeuse's upper chromosphere extends far beyond the circumstellar H alpha envelope of ~5 R*, determined from previous ground-based imaging. The flux in the broad and self-absorbed resonance lines of Mg II decreases by a factor of ~700 compared to the flux at chromospheric disk center. We observe strong asymmetry changes in the Mg II h and Si I resonance line profiles when scanning off-limb, signaling the outward acceleration of gas outflow in the upper chromosphere. From the radial intensity distributions of Fe I and Fe II emission lines we determine the radial non-LTE iron ionization balance. We compute that the local kinetic gas temperatures of the warm chromospheric gas component in the outer atmosphere exceed 2600 K, when assuming local gas densities of the cool gas component we determine from radiative transfer models that fit the 9.7 um silicate dust emission feature. The spatially resolved STIS spectra directly demonstrate that warm chromospheric plasma co-exisists with cool gas in Betelgeuse's circumstellar dust envelope.

In the summer and fall of 2000 the yellow hypergiant Rho Cassiopeiae dimmed by more than a visual magnitude, while its effective temperature decreased from \~7000 K to below 4000 K over ~200 d. We observed the highest mass-loss rate of ~5 % of the solar mass per year in a single stellar eruption so far (Lobel et al. 2003, ApJ, 583, 923). We determine from newly formed TiO bands a gas mass-loss rate of M_dot ~ 5.4 10^{-2} Msun/yr, which is of the same order of magnitude as has been proposed for the giant outbursts of the Luminous Blue Variable Eta Carinae. Over the past two years since the outburst event we observe a very prominent inverse P Cygni profile in Balmer H alpha. Strong Halpha emission has not before been observed in the hypergiant over this long period of time, signaling an unusual strong collapse of the upper hydrogen atmosphere, which we also observed in the months before the 2000 outburst. Since the brightness decrease of March 2003 we observe a transformation of the H alpha profile into a P Cygni profile, signaling a strong expansion of Rho Cas's upper atmosphere, which could possibly result in a new outburst event. High-resolution spectroscopic observations reveal that the lower photosphere rapidly expanded until May-June 2003. Very recent high-resolution observations of September 2003 show however that the photospheric absorption lines did not shift far blueward as observed in July 2000 before the strong brightness decrease of the outburst. The Fe I 5572 A is redshifted, signaling the collapse of the lower photosphere. A new strong brightness decrease by more than a magnitude in V for the fall and winter of 2003 is therefore not expected.

We present an overview of the spectral variability of the peculiar F-type hypergiant Rho Cas, obtained from our long-term monitoring campaigns over the past 8.5 years with four spectrographs in the northern hemisphere. Between 2000 June and September an exceptional variability phase occurred when the V-brightness dimmed by about a full magnitude. The star recovered from this deep minimum by 2001 April. It is the third outburst of Rho Cas on record in the last century. We observe TiO absorption bands in high-resolution near-IR spectra obtained with the Utrecht Echelle Spectrograph during the summer of 2000. TiO formation in the outer atmosphere occurred before the deep brightness minimum. Atmospheric models reveal that the effective temperature decreases by at least 3000 K, and the TiO shell is driven supersonically with Mdot ~= 5.4 10^{-2} Msun/yr. Strong episodic mass loss and TiO have also been observed during the outbursts of 1945-47 and 1985-86. A detailed analysis of the exceptional outburst spectra is provided, by comparing with high-resolution optical spectra of the early M-type supergiants Mu Cep (Ia) and Betelgeuse (Iab). The outburst spectra indicate the formation of a low-temperature, optically thick circumstellar gas shell of 3 10^{-2} Msun during ~200 d, caused by dynamic instability of the upper atmosphere of this pulsating massive supergiant near the Eddington luminosity limit. We observe that the mass-loss rate during the outburst is of the same order of magnitude as has been proposed for the outbursts of Eta Carinae. We present calculations that correctly predict the outburst time-scale, whereby the shell ejection is driven by the release of hydrogen ionization-recombination energy.

We present a high-resolution spectral analysis of Betelgeuse (M2 Iab). Between 1998 January and 1999 March four spatially resolved raster scans (http://cfa-www.harvard.edu/cfa/ep/pressrel/alobel0100.html) have been obtained with the STIS on the Hubble Space Telescope. The near-UV echelle spectra reveal double-peaked permitted emission lines of neutral and singly ionized metals, with self-absorbed line cores. We observe reversals in the intensity of both emission line components when scanning across the UV disk, for four unsaturated lines of Si I, Fe I, Al II], and Fe II. We model the Si I lam2516 resonance line with detailed non-LTE radiative transport calculations in spherical geometry, and constrain the mean velocity structure in the projected aperture area, for each scan position on the chromospheric disk. We infer the spatial velocity structure of Betelgeuse's extended chromosphere, which reveals localized upflows in the western front hemisphere in 1998 September, that expand further toward the eastern hemisphere in 1999 March. The spatial scans exhibit simultaneous up- and downflows across the lower chromosphere with mean velocities of ~2 km/s. We infer non-radial (or non-coherent) mass movements during certain phases of the stellar variability cycle from these subsonic flows. We present a discussion of constructing semi-empiric models for the chromosphere of this cool supergiant, and of its temporal variability.

Between 2000 June and September an exceptional variability phase occurred in the peculiar F-type hypergiant Rho Cas, when the V-brightness dimmed by at least a full magnitude. The star recovered from this deep minimum by 2001 April. It is the third outburst of Rho Cas on record in the last century. We observe TiO absorption bands in high-resolution optical and near-IR spectra obtained with the Utrecht Echelle Spectrograph during the summer of 2000. TiO formation in the outer atmosphere occurred before the deep brightness minimum.Atmospheric models reveal that the effective temperature decreases by at least 3000 K, and the TiO shell is driven supersonically with Mdot =~5.4 10^-2 Msun/yr. Strong episodic mass loss and TiO have also been observed during the outbursts of 1945-47 and 1985-86. An analysis of the exceptional outburst spectra of 2000-01 is provided, by comparing with high-resolution optical spectra of the early M-type supergiants Mu Cep (Ia) and Betelgeuse (Iab). During the outburst, central emission appears above the local continuum level in the split Na D lines. The outburst spectra indicate the formation of a low- temperature, optically thick circumstellar gas shell of 3 10^-2 Msun during \~200 d, caused by dynamic instability of the upper atmosphere of this pulsating massive supergiant near the Eddington luminosity limit. We present an equation that correctly predicts the outburst time-scale, whereby the shell ejection is driven by the release of hydrogen recombination energy. We observe that the mass-loss rate during the outburst is of the same order of magnitude as has been proposed for the outbursts of Eta Carinae. The research results in this paper are described in further detail in Lobel et al. (2003). A spectral movie sequence of the outburst is available at http://cfa-www.harvard.edu/~alobel/

This article is a brief review of the research by Dr. C. de Jager and co-workers over the past twenty years into the physics of hypergiant atmospheres. Various important results on the microturbulence, mass-loss, circumstellar environment, atmospheric velocity fields, and the Yellow Evolutionary Void of these enigmatic stars are summarized. Aspects of recent developments and future work are also communicated.

We have developed a new formalism to compute the thermodynamic coefficient Gamma1 in the theory of stellar and atmospheric stability. We generalize the classical derivation of the first adiabatic index, which is based on the assumption of thermal ionization and equilibrium between gas and radiation temperature, towards an expression which incorporates photo-ionization due to radiation with a temperature T_rad different from the local kinetic gas temperature.Our formalism considers the important non-LTE conditions in the extended atmospheres of supergiant stars. An application to the Kurucz grid of cool supergiant atmospheres demonstrates that models with T_rad =~ T_eff between 6500 K and 7500 K become most unstable against dynamic perturbations, according to Ledoux' stability integral <Gamma1>. This results from Gamma1 and <Gamma1> acquiring very low values, below 4/3, throughout the entire stellar atmosphere, which causes very high gas compression ratios around these effective temperatures. Based on detailed NLTE-calculations, we discuss atmospheric instability of pulsating massive yellow supergiants, like the hypergiant rho Cas (Ia+), which exist in the extension of the Cepheid instability strip, near the Eddington luminosity limit.

Four spatially resolved near-UV raster scans across the chromospheric disk of alpha Ori, obtained with the Space Telescope Imaging Spectrograph on the Hubble Space Telescope, reveal mean chromospheric infall from 1998 January to 1998 April, which reversed to upflow in deeper layers between 1998 September and 1999 March. In 1998 September we detect systematic reversals in the component maxima of four double-peaked emission lines of Si I (UV 1), Fe II (UV 36), Fe II (UV 61), and Al II] (UV 1), when scanning across the UV disk. Detailed modeling of the Si I lambda 2516 resonance line with radiative transport calculations in spherical geometry constrain the mean radial velocity structure in the projected slit area (25 by 100 mas) for different aperture positions, observed off-limb to 157.5 mas. Hence we determine with semi-empirical models that these spatial reversals of emission line components correspond with average opposite flow velocities of ~2 km/s across the chromospheric disk. We determine that the chromospheric velocity field can not be represented by a unique radial velocity structure across the stellar disk in order to match the observed peak ratios of this raster scan. These sub-sonic velocities indicate (local) non-radial movements of chromospheric fluid in confined regions during a chromospheric oscillation phase, which reverses from global contraction into expansion over this monitoring period of 15 months.

High-resolution near-ultraviolet spectra of the yellow hypergiants HR 8752 and rho Cassiopeiae indicate high effective temperatures placing both stars near the T_eff border of the ``yellow evolutionary void''. At present, the temperature of HR 8752 is higher than ever. For this star we found Teff=7900+-200 K, whereas rho Cassiopeiae has Teff=7300+-200 K. Both, HR 8752 and rho Cassiopeiae have developed strong stellar winds with Vinf ~ 120 km/s and Vinf ~ 100 km/s, respectively. For HR 8752 we estimate an upper limit for the spherically symmetric mass-loss of 6.7X10^{-6}M_solar/yr. Over the past decades two yellow hypergiants appear to have approached an evolutionary phase, which has never been observed before. We present the first spectroscopic evidence of the blueward motion of a cool super/hypergiant on the HR diagram.