<<

arXiv:astro-ph/0108358v1 22 Aug 2001 rttp bet fteylo yegat (F yellow the of objects prototype pcrltp n uioiy(eJgr1998). Jager the (de of those and than relativ broader type hyper a of significantly of spectral aspect are or spectroscopic lines important extended absorption Another an of loss. mass signature a is which otmsie(20 massive most n xeddcruselrenvelopes. circumstellar q extended u causes and velocit of which large-scale range developed unstable, strongly wide variability, are tion a exhibit Their and properties. limit, luminosity Eddington the hwoeo oebodeiso opnnsi h amrH Balmer the in propert components spectral emission their broad in more differ or they but one former, show the than nous a h sa Ia usual The Ia. ∼ n18 rf onlsd ae tteUiest fUrctp Utrecht of University the at Jager de Cornelis entitled Prof. 1980 In Introduction 1. fslce eeecst hi ok hc iletuethe enthuse will which work, astrophysics. their in howe to on subject co-workers, remarkable based references chronologi many developments, or selected and recent scientific Jager of important the de most all h the in Dr. highlight which complete by be answers time to (partial) of attempting course the the of past wh in f the some compilation has during reviews which a questions and ‘festschrift’ these is research, about 1980) progress massive Jager fundamental of st (de aspect of every upp monograph limit nearly certain His upper a the than brighter determine(s) stability? mechanism(s) no what are words, there why of question 000K hyaepyial ieetfo uegat flu of supergiants from different physically are They K. 10,000 olhprins(Ia hypergiants Cool et frcn eeomnsadftr okaeas commun also are work summa future and are developments stars recent enigmatic of these pects of Void Evolutionary phys Yellow veloc the mic atmospheric the environment, into on circumstellar years results mass-loss, twenty important Various past the atmospheres. pergiant over co-workers and Jager Abstract. USA Str MA, Garden 02138 60 Cambridge Astrophysics, for Center Harvard-Smithsonian Lobel Alex Research of Decades Two h rgts Stars Brightest The + − einto osntma htte r lasmr lumi- more always are they that mean not does designation 40 hsatcei re eiwo h eerhb r .de C. Dr. by research the of review brief a is article This M ⊙ + olsaspeetykon hs tr xs near exist stars These known. presently stars cool ) rI0 r ls fsprin tr with stars supergiant of class a are Ia0) or tfloe er fsinicrsac nothe into research scientific of years followed It . 1 − Ia G ρ ed,ecsiemass-loss, excessive fields, y a n R85 r the are 8752 HR and Cas + .Te r mn the among are They ). a eal.W also We details. cal aiproi pulsa- uasi-periodic v enproposed been ave insi httheir that is giants re h ai of basis the ormed w eae.This decades. two asaso similar of stars Ia iie number limited a la atmospheric ellar t ed,adthe and fields, ity rlmt nother In limit. er cmo stellar ncommon edrwt this with reader l ag aeof rate large ely e.Teestars These ies. bihdabook a ublished ioiyclass minosity α c addressed ich roturbulence, e without ver ieprofile, line eet, ie.As- rized. T c fhy- of ics icated. eff below We summarize three fundamental questions for the study of hypergiants: i. Why do hypergiants of similar bolometric luminosity and spectral type as normal supergiants display the signatures of strong large-scale photospheric mo- tions and enhanced mass loss? Why do these properties appear to be absent for normal supergiants? ii. What is the evolutionary status of the yellow hypergiants, and how do they fit into the broader scheme of massive star evolution? Are they the likely pro- genitors of (Type II) supernovae in the not too distant future? iii. Can an understanding of their atmospheric dynamics and circumstellar envi- ronments provide clues about the origin of the luminosity boundary for evolved stars? Do these dynamics characterize the hypergiants uniquely, or do they oc- cur modified in other classes of luminous stars as well?

2.

A few years after The Brightest Stars de Jager (1984) published a theoreti- cal paper which marked the beginning of a thorough study about the stabil- ity of hypergiant atmospheres, and which has consequently been pursued by many other workers in this field. It is put forward for the first time in this paper that atmospheric dynamics is a basic ingredient for understanding a temperature-dependent maximum luminosity limit extending across the entire upper Hertzsprung-Russell (HR) diagram. The classical Eddington limit is de- fined as the line in the HR-diagram where equals the force produced by the radiation pressure due to electron scattering. Neglecting other atmospheric opacity sources, this results in a nearly horizontal line in the HR-diagram for stars with a single mass-luminosity relation. Hence, the classic limit cannot ex- 6 plain the absence of cool supergiants with L⋆∼>10 L⊙, because hot stars that are more luminous do exist. De Jager proposed instead that the local stability of these extended atmospheres is greatly altered by an outward acceleration due to the gradient of turbulent pressure which results from the local field of turbu- lent motions. The importance of the latter scales with luminosity. It causes the effective acceleration, which is due to the balance between gravity, radiation- and turbulent-pressure, to produce higher mass-loss rates for the intrinsically brighter stars. The influence of turbulence is particularly important for the cool supergiants for which the radiation pressure remains rather small. For this rea- son the theoretical luminosity limit of cool stars is occasionally termed the ‘de Jager limit’ (for turbulent pressure). Studies of strong ‘microturbulence’ observed in the spectra of supergiants continue as an important subject of hypergiant research for obvious reasons. It increases the equivalent width of Fraunhofer lines by stochastic small-scale mo- tions with wavelengths smaller than the optical scale height of the atmosphere. The line-of-sight velocity component of these movements samples only a limited part of the over which the spectral lines are formed. A pertinent problem in pinpointing the physical nature of stellar microturbulence is that it is often observed to exceed the sound speed in supergiant atmospheres. Mi- croturbulence generally increases towards the hotter and more luminous stars.

2 It is frequently larger than the isothermal velocity of sound in the Ia and Iab supergiants. In 1991 de Jager chaired an international congress in Amsterdam about instabilities in evolved super- and hypergiants. At the conference he presented the idea that stellar microturbulence is strongly related to the physics of at- mospheric shock waves. Through numerical integrations of the hydrodynamics equations an observable microturbulent velocity is computed, assuming that the dependence of shock wave energy on the wavenumber is of the Kolmogoroff-type (a power law spectrum) (de Jager 1992a). Shock trains propagating outwards are a natural consequence of acoustic energy generation in cool stars with a subphotospheric zone. De Jager et al. (1997) subsequently showed that the thin heated layers trailing these shocks constitute the main part of the observed micro-‘turbulence’. A complex of successive hot sheets, distributed stochastically over the atmosphere, yields a stochastic distribution of thermal motions that, observationally, cannot be distinguished from a random field of small-scale hydrodynamic motions (de Jager 1998). The application of this theory to the ρ Cas revealed that a field of atmospheric shock waves which is critically supersonic, or have a mean Mach number only slightly exceeding unity, reproduces the observed supersonic microturbulence velocity of 11 km s−1. Remarkably, the hydrodynamic velocity field in the shocks remains very small, of the order of only 0.5 kms−1. Nev- ertheless, the momentum provided by such weak shock wave trains, assuming spherical geometric symmetry, is sufficient to initiate the large mass-loss rates −5 −1 (typically a few times 10 M⊙yr ) observed in cool hypergiants (de Jager 1996). These rates partly result from the very large stellar surface area (the mean photospheric radius of ρ Cas ranges between 400 and 500 R⊙), and the sonic point which is situated well inside the photosphere (de Jager et al. 2001). The field of stochastic shock waves provides on average outward directed mo- mentum which diminishes the effective atmospheric acceleration and increases the atmospheric scale height.

3. Mass Loss

Evolutionary calculations show that yellow supergiants evolve along tracks of nearly constant luminosity, beneath an observational luminosity limit. This limit has been called the Humphreys-Davidson (HD) limit, when it was pointed out that the luminosity of the brightest blue stars is about a factor six higher than that of the luminous red stars. Very evolved luminous supergiants move bluewards in the HR-diagram because they shed copious amounts of mass to the circumstellar/interstellar environment in the red supergiant phase. During redward evolution the very high mass-loss rates (M˙ ) reduce the stability of the convective layer, and below a critical envelope mass it contracts into a thinner radiative envelope, which causes rapid blueward evolution. For only a limited range of initial masses stars can become red supergiants, and evolve back far bluewards. The late stages in the fast evolution of these massive stars are rel- atively short. Typical values computed for stars with masses between 20 and 40 M⊙ are around one million years, with blueward tracks which last only a few

3 times 104 years. The high M˙ for very massive stars is however thought to have a decisive influence on these late evolutionary stages. De Jager et al. (1988) presented a comprehensive table of empirical mass- loss rates for 271 objects all over the HR-diagram, together with a parameter- ization which holds for stars from early to late spectral types (Nieuwenhuijzen & de Jager 1990). M˙ primarily depends on the stellar luminosity L⋆. For hot stars there appears to be a rather linear dependence, whereas K and M stars deviate by about 3 dex from the average hot star dependence. This points to other parameters than L⋆ (and Teff ) which influence the mass-loss rate (de Jager 1993). It seems imperative to ascribe the difference to the different phys- ical mechanisms which drive the massive stellar winds in the hot and in the cool supergiants. These differences are apparent from the emission line com- ponents observed with high dispersion in ρ Cas at either or both sides of the Doppler displaced Hα absorption core (Lobel et al. 1998). These complex line profiles indicate (quasi-) cyclic variations of the atmospheric circumstances and dynamics, which causes variable mass-loss rates that can occasionally increase −4 −1 to extreme values of a few 10 M⊙ yr . The temporal variability observed in photospheric absorption lines is not correlated with the changes of Hα; the latter is at least twice longer than the former. Since Hα forms higher, over a larger fraction of the atmosphere, it indicates that the atmosphere is dynami- cally velocity stratified, and the pulsation strongly enhances the density scale height compared to hydrostatic models. Yellow hypergiants exhibit ‘eruptions’, which occur on time-scales much longer than the pulsation quasi-period (ca. half a century, say). In an outburst of 1945-46 ρ Cas suddenly dimmed and displayed TiO-bands in its spectrum, characteristic of the photospheric temperatures of M-type stars. Within a couple of years the star brightened up by nearly a (April ’47), and a mid G-type spectrum was recovered around 1950. In 1985-86 the star showed a larger-than-average amplitude in the light curve, which can be associated with shell ejection events. Pulsation variability seems to be a common phenomenon in these stars. In his review on yellow hypergiants de Jager (1998) hypothesized that a distinction should be made between the periods in which these stars can be considered as fairly quietly pulsating, and when they are in the ‘explosive’ phase. During that phase the effective atmospheric acceleration practically vanishes, resulting in strongly enhanced mass ejection. Interestingly, it is suggested that the mechanism for outburst may be related to increased pulsation velocities during a period of instability that occurs when the star, during rapid blueward evolution, increases its pulsation amplitudes. The broader view emerges here that the exceptional atmospheric dynamics manifested by the yellow hypergiants is intricately related to their late evolutionary stage.

4. Circumstellar Environment

Little is known about the circumstellar environments of yellow hypergiants. The mid-IR silicate emission feature of ρ Cas is very weak, indicating the formation of dust only at very large distance from the surface. However, IRC+10420 is of later spectral type and possesses a very extended dust shell. It is commonly as- sumed that mass-loss from cool massive stars is primarily spherically symmetric.

4 The absence of observational evidence for strong magnetic fields in these super- giants and slow rotation does not readily support wind flows into a preferential direction. It is know that for most if not all stars wind flow is not spherically symmetric. Magnetic fields and rotation have a strong tendency to concentrate the gas that is leaving the star towards the equatorial plane, thus producing a disk rather than a spherical wind structure (de Jager 1993). HST-WFPC images of IRC+10420 reveal a circumstellar environment with a bipolar symmetry, and jet-like features and arcs produced primarily in the equatorial region. A bipolar wind has also been proposed for ρ Cas (Lobel 1997), based on the phenomenon of line core splitting observed in optical absorption lines of very low excitation energy. ρ Cas’ emission line spectrum is excited in the circumstellar environment with an excitation temperature of 3050 K. This emis- sion spectrum can emerge from a fast bipolar wind which collides with denser circumstellar (or perhaps interstellar) material, expelled during its violent mass- loss history. Similar split line cores were also observed for HR 8752 in 1960-70, when it was the spectroscopic twin of ρ Cas. The line core splitting phenomenon is of a bewildering complexity, because the central reversals are permitted emis- sion lines which appear in the core of metal absorption lines of resonance and subordinate transitions. ρ Cas’ emission spectrum is not excited by a normal stellar because its spectrum lacks the classic chromospheric indi- cators, except for Hα. Long-term spectroscopic monitoring of the split profiles reveals that the central emission peak is static, while the intensity of either of the adjacent absorption components periodically intensifies and weakens. This supports the formation of the emission line spectrum at large distance from the star, where it does not share the photospheric pulsation movements below. The blue-shifted appearance of the central emission cores indicates their formation in an optically thick wind. On the other hand, we have also discovered forbidden [Ca ii] emission lines in the spectrum, which reveals a very extended tenuous gas envelope surrounding ρ Cas. Prominent [N ii] line emission is also observed in HR 8752. Surprizingly however, the same metal core splitting phenomenon is observed in FU Ori objects. These young stellar objects have bipolar outflows, produced by disk , with powerful outbursts of mass-loss. They are pre- objects with optical spectra resembling F- and G-supergiants. Inter- estingly, it has also been proposed that FUor eruptions are due to perturbations in a circumstellar disk caused by the passage of a compact companion on an eccentric orbit. There is presently no indication for a compact companion in ρ Cas, but clear evidence for an early B-type main sequence star in HR 8752 was obtained from continuous near-UV spectra observed with IUE. Whether split absorption profiles are the signature of circumstellar disks which can obstruct strong wind outflow remains currently unanswered and re- quires further research. It appears however that they are spectral indicators of complex and dynamic circumstellar environments, which characterize the yellow hypergiants.

5 5. Atmospheric Velocity Fields

In 1990 de Jager published a short paper in Solar Physics. The paper is re- markable in that it provides an explanation for the problem of the ‘granulation boundary’ in the HR-diagram. This boundary has been inferred from the shape of spectral absorption line bisectors which incline oppositely when moving from the lower to the higher effective temperature regions in the HR-diagram. For the Ib supergiants the bisector turnover occurs around spectral types G0-5, but is poorly defined for the more luminous Ia and Ia+ stars. The spectra of F-type supergiants display uncommonly broad absorption lines which are ascribed to a well developed large-scale velocity field. The large macrobroadening cannot be attributed to rapid rotation (i.e. large vsini-values) for these evolved and very extended stars. For example, the modeling of line profiles in ρ Cas yields highly supersonic macroturbulent velocities of ∼25 kms−1. De Jager (1990) proposed that for stars to the left of the boundary, in the optically thin regions of the atmosphere above the , internal gravity waves develop which involve a large part of the line forming region. To the right of the granula- tion boundary the atmospheres are in convective equilibrium, and gravity waves would only affect the smallest optical depths, which hardly influences the line profiles. In a major paper which was published the year after, de Jager et al. (1991) provided the mathematical and physical foundation for these insights. The most important wavelengths in the spectrum of gravity waves are close to the atmo- spheric cut-off wavelength. The latter results from radiative damping and the atmospheric curvature. It is pointed out that in the highest atmospheric lev- els only the very low-mode (long wavelength) gravity waves are possible. For supergiant pulsation, high l modes are not expected because the wavelength of internal gravity waves propagating in these extended atmospheres should ex- ceed this cut-off wavelength. In deeper layers the damping cut-off wavelength becomes gradually shorter and higher gravity waves become possible. The more extended the atmosphere the lower the mode-number of gravity waves that may develop. However, gravity waves cannot propagate inside convectively stable regions where their minimum period (the Brunt-V¨ais¨al¨aperiod) becomes real and non-zero. Gravity waves propagate mainly in the horizontal direction, along the stellar circumference. They do not propagate strictly vertically. The application of this theory (de Jager 1998) to the yellow hypergiants ρ Cas and HR 8752, and to the white hypergiant HD 337579, provided a number of impressive theoretical results. Based on the diagnostic diagram, which plots the periods of pressure- and gravity-waves against their respective spectra of possible wavelengths and the observed quasi-period(s), he showed that these atmospheres are dominated by systems of gravity waves with in some cases shock waves superimposed. For ρ Cas the long-period waves are due to gravity waves of the shortest allowed wavelengths and periods, with fairly high l-values. The short-period components show the presence of pressure waves with wavelengths of about half the stellar radius. l-values exceeding unity are consistent with the finding for ρ Cas (Lobel et al. 1994), from an improved Baade-Wesselink pulsation test, that photospheric radial pulsations must be discarded to combine Teff -changes with the co-eval - and light-curve, observed over a complete variability period. The quasi-periodicity observed in the light-curves 6 of cool hypergiants therefore results from a complex, but law-abiding, interplay of atmospheric (non-radial) pulsation modes due to gravity and acoustic waves.

6. The Yellow Evolutionary Void

In 1995 Nieuwenhuijzen & de Jager published a classic paper on the atmospheric accelerations and stability of supergiant atmospheres. Their research provides a well-founded theoretical explanation for the conspicuously sparse amount of white (late A to early F-type) supergiants in the HR-diagram. At the Amster- dam conference de Jager (1992b) raised the question of why most of the known yellow hypergiants appear to cluster around 6500 K ≤ Teff ≤ 7500 K with differ- ent , while the area around 10,000 K and log(L⋆/L⊙)≥5.7 is nearly devoid of stars. Due to the natural deficiency of cool supergiants it is possibly hard to judge the reality of their absence around the upturn of the HD limit, but current population studies support only a minimal number of very luminous white supergiants. This area was tentatively baptized the ‘evolutionary void’ after detailed calculations revealed that hydrostatically stable solutions cannot be computed for the atmospheres of blueward evolving yellow hypergiants at a lower temperature boundary of Teff ≃8300 K. These calculations involve a time- independent solution of the momentum equation, which considers the Newtonian gravity acceleration derived from the evolutionary mass, the gas, radiation and turbulent pressure gradients, and the momentum of the . The latter requires the observed mass-loss rate. An ‘effective acceleration’ for the atmo- sphere is obtained iteratively, but which becomes negative at the cool border of the void. The outwards directed net force causes an unstable atmosphere above 8300 K. A remarkable byproduct of this study is that as stars move along their track, for time scales longer than the dynamic time scale of the atmosphere, the atmosphere continuously adapts to the new (L⋆, Teff )-values, and remains stable without reaching the Eddington limit. Current practice of determining the stability limit by extrapolating hydrostatic models to a ‘modified’ Eddington limit (i.e. utilizing a flux-mean opacity for the atmosphere which includes line opacity) is therefore not justified. When a real star slowly changes its effective gravity acceleration to zero during evolution the atmosphere puffs up, thereby decreasing the opacity to the classic electron scattering value. Over the past few years Drs. de Jager and Nieuwenhuijzen commenced a profound investigation of several individual cool hypergiants to which their theoretical developments are being applied (Nieuwenhuijzen & de Jager 2000). Atmospheric model parameters of HD 33579 have been determined via the ‘spec- tral diagnostics method’ and reveal that it is a redward evolving (high-mass) su- pergiant inside the void. For such a star the void is not forbidden. HR 8752 and IRC+10420 have lower masses and evolve bluewards, currently approaching the void. Several independent studies have demonstrated a significant overabun- dance of in ρ Cas. It is brought to the photosphere by dredging-up, which would support its post-red supergiant phase. De Jager & Nieuwenhuijzen (1997) suggested that the recurrent eruptions in yellow hypergiants occur when these stars approach the cool boundary of the void, and ‘bounce off’ redward (for a graph see Fig. 1 of Nieuwenhuijzen & de Jager 2000). This bouncing against the void may explain why most of the cool luminous hypergiants cluster

7 near its low-temperature boundary, while the identification of hypergiants of later spectral type, possibly like VY CMa (of M-type), is seldom. An increase of the photospheric temperature of HR 8752 by 3000 K−4000 K has been ob- served based on high-resolution spectra collected over the past thirty years. The effective temperature of IRC+10420 has been as low as 5270 K and increased to ∼7000 K in 1992. The natural question arises whether (apparent) evolution- ary changes on human time-scales result from an actual active reconstruction of the stellar interior, or by some hitherto unknown envelope mechanism? Is there a global instability mechanism in these stars other than , which changes their spectral types, perhaps triggered by some internal or exter- nal cause? Is the void an eternal obstacle for their evolution and will the yellow hypergiants eventually fade as supernovae, or will they rapidly evolve across the hypergiant void into early A-type supergiants (e.g. the least luminous S Dor stars), after having lost large amounts of mass? Time will tell.

7. Future Work

As is the case with all research at the forefront of science, the present-day de- velopments are the most exciting. To avoid disclosing research results prior to the main publication (de Jager et al. 2001), we summarize some important general aspects of hypergiant atmospheric stability, and the lines along which further investigations are planned. New studies show that different areas in the upper HR-diagram can be delineated where the effective acceleration in the upper part of the photosphere is negative, or directed outwards. Stars like ρ Cas, HR 8752, and IRC+10420 are situated close to this line. Another aspect of their dynamic stability is that the sonic point (the point where the wind velocity equals the sound speed) of the stellar wind is situated inside the photosphere. Matter which is driven beyond this point escapes freely into space, causing high mass-loss rates. Parallel to these developments is the study of the behavior of the adiabatic indices in the atmospheres of luminous supergiants. An early study (Lobel et al. 1992) showed that Γ1 assumes values below 4/3 in low-gravity hydrostatic model atmospheres. This corresponds to zones of increased atmospheric destabilization due to partial hydrogen and helium ionization. The study is presently further improved by evaluating Ledoux’ stability integral <Γ1> for realistic envelope models (Lobel 2001), including the effects of photo-ionization and the dilution of radiation pressure in the extended atmospheres of these stars. Important dependencies on the atmospheric acceleration and effective temperature have been demonstrated so far. Future work will focus on the interrelations between these different stability criteria, and their various specific implications. It appears to us that consider- able progress will be accomplished by integrating the fully time-dependent hy- drodynamic equations, and incorporating the various acceleration mechanisms in the equation of motion, together with the increased atmospheric compress- ibility due to partial LTE- and NLTE-ionization. These simulations could for example model the velocity stratification observed for ρ Cas’ atmosphere. The further refinement of regions of dynamic instability in the upper portion of the

8 Figure 1. Kees de Jager inspecting the Observatori del Teide at Izana on the island of Tenerife (Canary Islands, Oct. 1999). The Gregory Coud´esolar Telescope towers in the background.

HR-diagram, and a continuous spectral monitoring program of HR 8752 (de Jager et al. 2001), for which an outburst seems imminent, is in progress. To conclude this overview, it is safe to say that our understanding of hy- pergiants has considerably changed and improved over the past two decades. At the same time, lots of new and pertinent questions have been raised about these exceptional objects. We owe much of this progress to the dedicated and continuing work of one person; Cornelis de Jager, who is a guide and has been the mentor of many researchers in and astrophysics. Acknowledgments. This article honors the eightieth birthday of Prof. C. de Jager. The author has been a Ph.D. student at the Space Research Organi- zation of the Netherlands - Utrecht and the Free University Brussels.

References de Jager, C., The Brightest Stars, Geophysics and Astrophysics Monographs, v. 19, Reidel, Dordrecht, 1980 de Jager, C., The Stability Limit of Hypergiant , Astronomy & Astrophysics, 138, 246, 1984 de Jager, C., Nieuwenhuijzen, H., and van der Hucht, K. A., Mass Loss Rates in the Hertzsprung-Russell Diagram, Astronomy & Astrophysics Suppl. Ser. , 72, 259, 1988 9 de Jager, C., An Explanation of the ’Granulation Boundary’ in the HR Diagram, Solar Physics, 126, 201, 1990 de Jager, C., de Koter, A., Carpay, J., and Nieuwenhuijzen, H., The Atmospheric Motion Field in Super- and Hypergiants, Astronomy & Astrophysics, 244, 131, 1991 de Jager, C., Spectroscopic Diagnostics of a Field of Shock Waves, In C. de Jager and H. Nieuwenhuijzen, editors, Instabilities in Evolved Super- and Hypergiants, pages 129 – 130., Proceedings of the Royal Dutch Academy of Sciences, Sect. Physics, First Series, part 36., Amsterdam, The Nether- lands, 26 February - 1 March 1991. ISBN 0-444-85736-2, North Holland, Amsterdam, 1992 de Jager, C., Dynamics and Stability of Supergiant Atmospheres, In C. de Jager and H. Nieuwenhuijzen, editors, Instabilities in Evolved Super- and Hy- pergiants, pages 98 – 103., Proceedings of the Royal Dutch Academy of Sciences, Sect. Physics, First Series, part 36. Amsterdam, The Nether- lands, 26 February - 1 March 1991. ISBN 0-444-85736-2, North Holland, Amsterdam, 1992 de Jager, C., Massive Stars: Setting the Stage, Space Science Reviews, 66, Nos. 1-4, 1993 de Jager, C., A White and a Red Star (Special Historical Review), Solar Physics, 169, 443, 1996 de Jager, C., Lobel, A., and Israelian, G., Turbulence, mass loss and Hα emis- sion by stochastic shocks in the hypergiant ρ Cassiopeiae, Astronomy & Astrophysics, 325, 714, 1997 de Jager, C., and Nieuwenhuijzen, H., An Obstacle to the Late Evolution of Massive Stars, Monthly Notices of the Royal Astronomical Society, 290, L50, 1997 de Jager, C., The Yellow Hypergiants, Astronomy & Astrophysics Rev., 8, 145, 1998 de Jager, C., Lobel, A., Nieuwenhuijzen, H., and Stothers, R., Instability Regions in the Upper HR Diagram, Monthly Notices of the Royal Astronomical Society, 2001, in press Lobel, A., Achmad, L., de Jager, C., and Nieuwenhuijzen, H., On the Instability Parameters of Stellar Atmospheres, Astronomy & Astrophysics, 264, 147, 1992 Lobel, A., de Jager, C., Nieuwenhuijzen, H., Smolinski, J., and Gesicki, K., Pulsation of the yellow hypergiant ρ Cassiopeiae in 1970, Astronomy & Astrophysics, 291, 226, 1994 Lobel, A., Pulsation and Atmospherical Instability of Luminous F- and G-type Stars, Ph.D. thesis, Free University Brussels, Belgium, 1997, ISBN 90- 423-0014-0, Shaker Publishing, Maastricht 10 Lobel, A., Israelian, G., de Jager, C., Musaev, F., Parker, J. Wm., and Mavro- giorgou, A., The spectral variability of the cool hypergiant ρ Cassiopeiae, Astronomy & Astrophysics, 330, 659, 1998 Lobel, A., On the Dynamic Stability of Cool Supergiant Atmospheres, The As- trophysical Journal, 2001, 558, in press for Sept. 10 issue Nieuwenhuijzen, H., and de Jager, C., Parameterization of stellar rates of mass loss as functions of the fundamental stellar parameters M, L, and R, Astronomy & Astrophysics, 231, 134, 1990 Nieuwenhuijzen, H., and de Jager, C., Atmospheric Accelerations and the Sta- bility of Dynamic Supergiant Atmospheres, Astronomy & Astrophysics, 302, 811, 1995 Nieuwenhuijzen, H., and de Jager, C., Checking the yellow evolutionary void. Three evolutionary critical Hypergiants: HD 33579, HR 8752 & IRC+10420, Astronomy & Astrophysics, 353, 163, 2000

11