We assess the ability to image the photosphere of red supergiants and, in particular Betelgeuse, through the modelling of the observed linear polarization in atomic spectral lines. We also aim to analyse the resulting images over time, to measure the size and dynamics of the convective structures in these stars. Rayleigh scattering polarizes the continuum and spectral lines depolarize it.This depolarization is seen as a linear polarization signal parallel to the radial direction on the stellar disk. Integrated over the disk, it would result in a null signal, except if brightness asymmetries/inhomogeneities are present. This is the basic concept behind our imaging technique. The several tests and comparisons performed prove that our technique reliably retrieves the salient brightness structures in the photosphere of Betelgeuse, and should be relevant to other red supergiants. We demonstrate that these structures are convective cells, with a characteristic size of more than 60% of the stellar radius. We also derive the characteristic upflow and downflow speeds, 22 and 10 km/s respectively. We find weak magnetic fields concentrated in the downflow lanes in between granules. We follow those convective structures in time. Changes happen on timescales of one week, but individual structures can be tracked over four years of observations.

Betelgeuse is an M-type supergiant that presents a circularly polarized (Stokes V) signal in its line profiles, interpreted in terms of a surface magnetic field. The weak circular polarization signal has been monitored over 7.5 years in order to follow its evolution on different timescales, and eventually to determine its physical origin. Linear polarization measurements have also been obtained regularly in the last few years. We used both the ESPaDOnS and Narval spectropolarimeters to obtain high signal-to-noise ratio (S/N) spectra, which were processed by means of the least-squares deconvolution (LSD) method. In order to ensure the reality of the very weak circular polarization, special care has been taken to limit instrumental effects. In addition, several tests were performed on the Stokes V signal to establish its stellar and Zeeman origin. We confirm the magnetic nature of the circular polarization, pointing to a surface magnetic field of the order of 1G. The Stokes V profiles present variations over different timescales, the most prominent one being close to the long secondary period (LSP; around 2000d for Betelgeuse) often invoked in red evolved stars. This long period is also dominant for all the other Stokes parameters. The circular polarization is tentatively modeled by means of magnetic field concentrations mimicking spots, showing in particular that the velocity associated with each "spot" also follows the long timescale, and that this signal is nearly always slightly redshifted. From the coupled variations of both linear and circular polarization signatures in amplitude, velocity and timescale, we favour giant convection cells as the main engine at the origin of polarization signatures and variations in all the Stokes parameters. This strengthens support for the hypothesis that large convective cells are at the origin of the LSP.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

RSG stars are very massive cool evolved stars. Recently, a weak magnetic field was measured at the surface of $\alpha$ Ori and this is so far the only M-type supergiant for which a direct detection of a surface magnetic field has been reported. By extending the search for surface magnetic field in a sample of late-type supergiants, we want to determine whether the surface magnetic field detected on $\alpha$ Ori is a common feature among the M-type supergiants. With the spectropolarimeter Narval at TBL we undertook a search for surface magnetic fields in a sample of cool supergiant stars, and we analysed circular polarisation spectra using the least-squares deconvolution technique. We detect weak Zeeman signatures of stellar origin in the targets CE Tau, $\alpha^1$ Her and $\mu$ Cep. For the latter star, we also show that cross-talk from the strong linear polarisation signals detected on this star must be taken into account. For CE Tau and $\mu$ Cep, the longitudinal component of the detected surface fields is at the Gauss-level, such as in $\alpha$~Ori. We measured a longitudinal field almost an order of magnitude stronger for $\alpha^1$ Her. We also report variability of the longitudinal magnetic field of CE Tau and $\alpha^1$ Her, with changes in good agreement with the typical atmospheric dynamics time-scales. We also report a non-detection of magnetic field at the surface of the yellow supergiant star $\rho$ Cas. The two RSG stars of our sample, CE Tau and $\mu$ Cep, display magnetic fields very similar to that of $\alpha$ Ori. The non-detection of a magnetic field on the post-RSG star $\rho$ Cas suggests that the magnetic field disappears, or at least becomes undetectable with present methods, at later evolutionary stages. Our analysis of $\alpha^1$ Her supports the proposed reclassification of the star as an AGB star.

In the solar limb, linear polarisation is due to anisotropy of the radiation field induced by limb darkening. It is maximal when it is seen parallel to the limb and it vanishes when it is integrated over the spherically-symmetric solar disk. Therefore for distant stars, that present spherical symmetry, linear polarisation signatures are very difficult to observe. However strong linear polarisation features have been reported in the prototypical red supergiant star alpha Ori (Betelgeuse). With an analytical model we propose to explain them.

Betelgeuse is an M supergiant that harbors spots and giant granules at its surface and presents linear polarization of its continuum. We have previously discovered linear polarization signatures associated with individual lines in the spectra of cool and evolved stars. Here, we investigate whether a similar linearly polarized spectrum exists for Betelgeuse. We used the spectropolarimeter Narval, combining multiple polarimetric sequences to obtain high signal-to-noise ratio spectra of individual lines, as well as the least-squares deconvolution (LSD) approach. We have discovered the existence of a linearly polarized spectrum for Betelgeuse, detecting a rather strong signal (at a few times 10$^{-4}$ of the continuum intensity level), both in individual lines and in the LSD profiles. Studying its properties and the signal observed for the resonant \ion{Na}{i}\,D lines, we conclude that we are mainly observing depolarization of the continuum by the absorption lines. The linear polarization of the Betelgeuse continuum is due to the anisotropy of the radiation field induced by brightness spots at the surface and Rayleigh scattering in the atmosphere. We have developed a geometrical model to interpret the observed polarization, from which we infer the presence of two brightness spots and their positions on the surface of Betelgeuse. We show that applying the model to each velocity bin along the Stokes Q and U profiles allows the derivation of a map of the bright spots. We use the Narval linear polarization observations of Betelgeuse obtained over a period of 1.4 years to study the evolution of the spots. Our study of the linearly polarized spectrum of Betelgeuse provides a novel method for studying the evolution of brightness spots at its surface and complements quasi-simultaneous observations obtained with PIONIER at the VLTI.

This article presents an on-line tool (rttools.irap.omp.eu) and its accompanying software ressources for the numerical solution of basic radiation transfer out of local thermodynamic equilibrium (LTE). State-of-the-art stationary iterative methods such as Accelerated $\Lambda$-Iteration and Gauss-Seidel schemes, using a short characteristics-based formal solver are used. We also comment on typical numerical experiments associated to the basic non-LTE radiation problem. These ressources are intended for the largest use and benefit, in support to more classical radiation transfer lectures usually given at the Master level.

Context: The solution of the nonlocal thermodynamical equilibrium (non-LTE) radiative transfer equation usually relies on stationary iterative methods, which may falsely converge in some cases. Furthermore, these methods are often unable to handle large-scale systems, such as molecular spectra emerging from, for example, cool stellar atmospheres. Aims: Our objective is to develop a new method, which aims to circumvent these problems, using nonstationary numerical techniques and taking advantage of parallel computers. Methods: The technique we develop may be seen as a generalization of the coupled escape probability method. It solves the statistical equilibrium equations in all layers of a discretized model simultaneously. The numerical scheme adopted is based on the generalized minimum residual method. Result:. The code has already been applied to the special case of the water spectrum in a red supergiant stellar atmosphere. This demonstrates the fast convergence of this method, and opens the way to a wide variety of astrophysical problems.

The structures of the outer atmospheres of red giants are very complex. The notion of large optically thick molecular spheres around the stars (MOLspheres) has been invoked in order to explain e.g. spectro-interferometric observations. However, high-resolution spectra in the mid-IR do not easily fit into this picture. They rule out any large sphere of water vapour in LTE surrounding red giants. Our aim here is to investigate high-resolution, mid-infrared spectra for a range of red giants, from early-K to mid M. We have recorded 12 microns spectra of 10 well-studied bright red giants, with TEXES on the IRTF. We find that all giants in our study cooler than 4300 K, spanning a range of effective temperatures, show water absorption lines stronger than expected. The strengths of the lines vary smoothly with spectral type. We identify several spectral features in the wavelength region that undoubtedly are formed in the photosphere. From a study of water-line ratios of the stars, we find that the excitation temperatures, in the line-forming regions, are several hundred Kelvin lower than expected from a classical photospheric model. This could either be due to an actually lower temperature structure in the outer regions of the photospheres caused by, for example, extra cooling, or due to non-LTE level populations, affecting the source function and line opacities. We have demonstrated that these diagnostically interesting water lines are a general feature of red giants across spectral types, and we argue for a general explanation of their formation rather than explanations requiring specific properties. Since the water lines are neither weak (filled in by emission) nor appear in emission, as predicted by LTE MOLsphere models in their simplest forms, the evidence for the existence of such large optically-thick, molecular spheres enshrouding the stars is weakened. (abbreviated)

The physical structures of the outer atmospheres of red giants are not known. They are certainly complex and a range of recent observations are showing that we need to embrace to non-classical atmosphere models to interpret these regions. This region's properties is of importance, not the least, for the understanding of the mass-loss mechanism for these stars, which is not still understood. Here, we present observational constraints of the outer regions of red giants, based on mid-IR, high spectral resolution spectra. We also discuss possible non-LTE effects and highlight a new non-LTE code that will be used to analyse the spectra of these atmospheric layers. We conclude by mentioning our new SOFIA/EXES observations of red giants at 6 microns, where the vibration-rotation lines of water vapour can be detected and spectrally resolved for the first time.

The growth of luminous structures and the building blocks of life in the Universe began as primordial gas was processed in stars and mixed at galactic scales. The mechanisms responsible for this development are not well understood and have changed over the intervening 13 billion years. To follow the evolution of matter over cosmic time, it is necessary to study the strongest (resonance) transitions of the most abundant species in the Universe. Most of them are in the ultraviolet (UV; 950A-3000A) spectral range that is unobservable from the ground. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the Universe. Habitable planets grow in protostellar discs under ultraviolet irradiation, a by-product of the star-disk interaction that drives the physical and chemical evolution of discs and young planetary systems. The electronic transitions of the most abundant molecules are pumped by the UV field, providing unique diagnostics of the planet-forming environment that cannot be accessed from the ground. Earth's atmosphere is in constant interaction with the interplanetary medium and the solar UV radiation field. A 50-100 times improvement in sensitivity would enable the observation of the key atmospheric ingredients of Earth-like exoplanets (carbon, oxygen, ozone), provide crucial input for models of biologically active worlds outside the solar system, and provide the phenomenological baseline to understand the Earth atmosphere in context. In this white paper, we outline the key science that such a facility would make possible and outline the instrumentation to be implemented.

The interpretation of water lines in red supergiant stellar atmospheres has been much debated over the past decade. The introduction of the so-called MOLspheres to account for near-infrared "extra" absorption has been controversial. We propose that non-LTE effects should be taken into account before considering any extra-photospheric contribution. After a brief introduction on the radiative transfer treatment and the inadequacy of classical treatments in the case of large-scale systems such as molecules, we present a new code, based on preconditioned Krylov subspace methods. Preliminary results suggest that NLTE effects lead to deeper water bands, as well as extra cooling.

Mass loss of red supergiants (RSG) is important for the evolution of massive stars, but is not fully explained. Several empirical prescriptions have been proposed, trying to express the mass-loss rate (Mdot) as a function of fundamental stellar parameters (mass, luminosity, effective temperature). Our goal is to test whether the de Jager et al. (1988) prescription, used in some stellar evolution models, is still valid in view of more recent mass-loss determinations. By considering 40 Galactic RSGs presenting an infrared excess and an IRAS 60-mu flux larger than 2 Jy, and assuming a gas-to-dust mass ratio of 200, it is found that the de Jager rate agrees within a factor 4 with most Mdot estimates based on the 60-mu signal. It is also in agreement with 6 of the only 8 Galactic RSGs for which Mdot can be measured more directly through observations of the circumstellar gas. The two objects that do not follow the de Jager prescription (by an order of magnitude) are mu Cep and NML Cyg. We have also considered the RSGs of the Magellanic Clouds. Thanks to the works of Groenewegen et al. (2009) and Bonanos et al. (2010), we find that the RSGs of the SMC have Mdots consistent with the de Jager rate scaled by (Z/Zsun)**(alpha), where Z is the metallicity and alpha is 0.7. The situation is less clear for the LMC RSGs. In particular, for luminosties larger than 1.6E+05 Lsun, one finds numerous RSGs (except WOH-G64) having Mdot significantly smaller than the de Jager rate, and indicating that Mdot would no longer increase with L. Before this odd situation is confirmed through further analysis of LMC RSGs, we suggest to keep the de Jager prescription unchanged at solar metallicity in the stellar evolutionary models and to apply a (Z/Zsun)**0.7 dependence.

Synthetic spectra are needed to determine fundamental stellar and wind parameters of all types of stars. They are also used for the construction of theoretical spectral libraries helpful for stellar population synthesis. Therefore, a database of theoretical spectra is required to allow rapid and quantitative comparisons to spectroscopic data. We provide such a database offering an unprecedented coverage of the entire Hertzsprung-Russell diagram. We present the POLLUX database of synthetic stellar spectra. For objects with Teff < 6 000 K, MARCS atmosphere models are computed and the program TURBOSPECTRUM provides the synthetic spectra. ATLAS12 models are computed for stars with 7 000 K <Teff <15 000 K. SYNSPEC gives the corresponding spectra. Finally, the code CMFGEN provides atmosphere models for the hottest stars (Teff > 25 000 K). Their spectra are computed with CMF_FLUX. Both high resolution (R>150 000) optical spectra in the range 3 000 to 12 000 A and spectral energy distributions extending from the UV to near--IR ranges are presented. These spectra cover the HR diagram at solar metallicity. We propose a wide variety of synthetic spectra for various types of stars in a format that is compliant with the Virtual Observatory standards. A user--friendly web interface allows an easy selection of spectra and data retrieval. Upcoming developments will include an extension to a large range of metallicities and to the near--IR high resolution spectra, as well as a better coverage of the HR diagram, with the inclusion of models for Wolf-Rayet stars and large datasets for cool stars. The POLLUX database is accessible at http://pollux.graal.univ-montp2.fr/ and through the Virtual Observatory.

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with CO5BOLD to various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection pattern on this star. Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive interferometric observables using the post-processing radiative transfer code OPTIM3D. The synthetic visibility curves and closure phases are compared to observations. Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (30 mas) convective cell. In this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null of the visibility curve (i.e. the apparent stellar radius).

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.

We report the detection, from observations using the James Clerk Maxwell Telescope, of CO J $=$ 3$\to$ 2 transition lines in six carbon stars, selected as members of the Galactic Halo and having similar infrared colors. Just one Halo star had been detected in CO before this work. Infrared observations show that these stars are red (J-K $>$3), due to the presence of large dusty circumstellar envelopes. Radiative transfer models indicates that these stars are losing mass with rather large dust mass-loss rates in the range 1--3.3 $\times$$10^{-8}$M$_{\odot}$yr$^{-1}$, similar to what can be observed in the Galactic disc. We show that two of these stars are effectively in the Halo, one is likely linked to the stream of the Sagittarius Dwarf Spheroidal galaxy (Sgr dSph), and the other three stars certainly belong to the thick disc. The wind expansion velocities of the observed stars are low compared to carbon stars in the thin disc and are lower for the stars in the Halo and the Sgr dSph stream than in the thick disc. We discuss the possibility that the low expansion velocities result from the low metallicity of the Halo carbon stars. This implies that metal-poor carbon stars lose mass at a rate similar to metal-rich carbon stars, but with lower expansion velocities, as predicted by recent theoretical models. This result implies that the current estimates of mass-loss rates from carbon stars in Local Group galaxies will have to be reconsidered.

The post-main sequence evolution of stars of intermediate or large masses is notoriously complex. In the recent past, a number of workshops and meetings have focused on either the Asymptotic Giant Branch of intermediate mass stars, or the evolution of massive stars. But how well defined is the boundary between these categories of objects defined? How would an observer proceed to classify stars into one or the other category? How do objects near the boundary evolve, die, and contribute to the chemical evolution of their environment? During this 3-day international workshop, 26 high quality presentations were given by specialists in the relevant fields of astrophysics, and stimulating discussions followed. It is technically impossible to provide an exhaustive census of the results and ideas that emerged. In this brief article, we choose to point to key elements of the workshop, some of which are now the topic of new collaborations and will lead to publications elsewhere. For the sake of brevity, we deliberately cite only the contributors to the workshop and no external references. Many bibliographic references can be found in the original presentations, which can be retrieved through: http://astro.u-strasbg.fr/observatoire/obs/stars2009/stars2009.html The programme workshop, which includes the titles of the individual contributions, is provided as an appendix.

Context. It has been suggested that convection in Red Supergiant (RSG) stars gives rise to large-scale granules causing observable surface inhomogeneities. This convection is also extremely vigorous, and suspected to be one of the causes of mass-loss in RSGs. It must thus be understood in details. Evidence has been accumulated that there are asymmetries in the photospheres of RSGs, but detailedstudies of granulation are still lacking. Interferometric observations offer an exciting possibility to tackle this question, but they are still often interpreted using smooth symmetrical limb-darkened intensity distributions, or very simple spotted ad hoc models. Aims. We explore the impact of the granulation on visibility curves and closure phases using the radiative transfer code OPTIM3D. We simultaneously assess how 3D simulations of convection in RSG with CO5BOLD can be tested against these observations. Methods. We use 3D radiative-hydrodynamics (RHD) simulations of convection to compute intensity maps at various wavelengths and time, from which we derive interferometric visibility amplitudes and phases. We study their behaviour with time, position angle, and wavelength, and compare them to observations of the RSG alpha Ori Results. We provide average limb-darkening coefficients for RSGs. We detail the prospects for the detection and characterization of granulation (contrast, size) on RSGs. We demonstrate that our RHD simulations provide an excellent fit to existing interferometric observation of alpha Ori, contrary to limb darkened disks. This confirms the existence of large convective cells on the surface of Betelgeuse.

The GIII red giant star epsilon Oph has been found to exhibit several modes of oscillation by the MOST mission. We interpret the observed frequencies of oscillation in terms of theoretical radial p-mode frequencies of stellar models. Evolutionary models of this star, in both shell H-burning and core He-burning phases of evolution, are constructed using as constraints a combination of measurements from classical ground-based observations (for luminosity, temperature, and chemical composition) and seismic observations from MOST. Radial frequencies of models in either evolutionary phase can reproduce the observed frequency spectrum of epsilon Oph almost equally well. The best-fit models indicate a mass in the range of 1.85 +/- 0.05 Msun with radius of 10.55 +/- 0.15 Rsun. We also obtain an independent estimate of the radius of epsilon Oph using high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR instrument. The measured limb darkened disk angular diameter of epsilon Oph is 2.961 +/- 0.007 mas. Together with the Hipparcos parallax, this translates into a photospheric radius of 10.39 +/- 0.07 Rsun. The radius obtained from the asteroseismic analysis matches the interferometric value quite closely even though the radius was not constrained during the modelling.

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

We have written a 3D radiative transfer code that computes emerging spectra and intensity maps. We derive from radiative hydrodynamic (RHD) simulations of RSG stars carried out with CO5BOLD (Freytag et al. 2002) observables expected for red supergiant stars (RSG) especially for interferometric observations, with emphasis on small scale structures. We show that the convection-related surface structures are detectable in the H band with today's interferometers and that the diameter measurement should not be too dependent on the adopted model. The simulations are a great improvement over parametric models for the interpretation of interferometric observations.

Red supergiant stars represent a key phase in the evolution of massive stars. Recent radiative hydrodynamic simulations suggest that their atmospheres may be the location of large-scale convective motions. As supergiant convection is expected to generate supersonic motions and shocks, we seek constraints on these atmospheric motions and their possible relation with mass-loss rates. We present high-resolution, visible spectroscopy of a sample of red supergiants (spectral type M I) and analyse them with a tomographic technique. We observe steep velocity gradients, characterising both upward and downward supersonic motions, which are time variable on time scales of a few hundred days. These convective motions will generate turbulent pressure, which will strongly decrease the effective gravity. We suggest that this decrease, combined with radiative pressure on molecular lines, initiate the mass loss in red supergiant stars.

Deriving the physical properties of red supergiants (RSGs) depends upon accurate corrections for reddening by dust. We use our recent modeling of the optical spectra of RSGs to address this topic. We find: (1) Previous broad-band studies have underestimated the correction for extinction in the visible, and hence the luminousities, if derived from V. (2) A significant fraction of RSGs in Galactic OB associations and clusters show up to several magnitudes of excess visual extinction compared to OB stars in the same regions; we argue that this is likely due to circumstellar dust around the RSGs. (3) RSGs contribute dust grains at the rate of $3 \times 10^{-8} M_\odot$ yr$^{-1}$ kpc$^{-2}$ in the solar neighborhood, comparable to what we estimate for late-type WC Wolf-Rayet stars, $1 \times 10^{-7} M_\odot$ yr$^{-1}$ kpc$^{-2}$. In the solar neighborhood this represents only a few percent of the dust production (which is dominated by low-mass AGBs), but we note that in low-metallicity starbursts, dust production by RSGs would likely dominate over other sources.

We use moderate-resolution optical spectrophotometry and the new MARCS stellar atmosphere models to determine the effective temperatures of 74 Galactic red supergiants. From these we find a new effective temperature scale that is significantly warmer than those in the literature. We show that this temperature scale, along with the newly derived bolometric corrections, gives much better agreement between our red supergiants and stellar evolutionary tracks. This agreement provides an independent verification of our new temperature scale. The combination of effective temperature and bolometric luminosities allows us to calculate stellar radii; the coolest and most luminous stars have radii of roughly 1500 solar radii (7 AU), in excellent accordance with the largest stellar radii predicted from current evolutionary theory. We find that similar results are obtained for the effective temperatures and bolometric luminosities using only the de-reddened V-K colors, providing a powerful demonstration of the self-consistency of the MARCS models.

We report high-sensitivity millimetre observations of several molecular species (13CO, HCN, HNC, CN, HCO+ and N2H+) in a sample of compact planetary nebulae. Some species such as HCO+ and CN are particularly abundant compared to envelopes around AGB stars or even interstellar clouds. We have estimated the following average values for the column densities ratios: CN/HCN~2.6, HCO+/HCN~0.5, and HNC/HCN~0.4. Thus, the chemical composition of the molecular envelopes in these compact PNe appears somewhat intermediate between the composition of proto-PNe (such as CRL 2688 or CRL 618) and well evolved PNe (such as the Ring, M4--9, or the Helix). From observations of the CO isotopomers, we have estimated that the 12C/13C ratio is in the range 10 ~< 12C/13C ~< 40. These values are below those expected from standard asymptotic giant branch models and suggest non-standard mixing processes. The observed molecular abundances are compared to very recent modelling work, and we conclude that the observations are well explained, in general terms, by time-dependent gas-phase chemical models in which the ionization rate is enhanced by several orders of magnitude with respect to the average interstellar value. Thus, our observations confirm that the chemistry in the neutral shells of PNe is essentially governed by the high energy radiation from the hot central stars. The complexity of the chemical processes is increased by numerous factors linked to the properties of the central star and the geometry and degree of clumpiness of the envelope. Several aspects of the PN chemistry that remains to be understood are discussed within the frame of the available chemical models.