Dust and the spectral energy distribution of the OH/IR star OH 127.8+0.0: Evidence for circumstellar metallic iron

We present a fit to the spectral energy distribution of OH 127.8+0.0, a typical asymptotic giant branch star with an optically thick circumstellar dust shell. The fit to the dust spectrum is achieved using non-spherical grains consisting of metallic iron, amorphous and crystalline silicates and water ice. Previous similar attempts have not resulted in a satisfactory fit to the observed spectral energy distributions, mainly because of an apparent lack of opacity in the 3--8 micron region of the spectrum. Non-spherical metallic iron grains provide an identification for the missing source of opacity in the near-infrared. Using the derived dust composition, we have calculated spectra for a range of mass-loss rates in order to perform a consistency check by comparison with other evolved stars. The L-[12 micron] colours of these models correctly predict the mass-loss rate of a sample of AGB stars, strengthening our conclusion that the metallic iron grains dominate the near-infrared flux. We discuss a formation mechanism for non-spherical metallic iron grains.Comment: 10 pages, 6 figures, accepted for publication by A&

Does the momentum flux generated by gravitational contraction drive AGB mass-loss?

Gravitational contraction always generates a radially directed momentum flux. A particularly simple example occurs in the electron-degenerate cores of AGB stars, which contract steadily under the addition of helium ashes from shell hydrogen burning. The resulting momentum flux is quantified here. And since the cores of AGB stars lack efficient momentum cancellation mechanisms, they can maintain equilibrium by exporting their excess momentum flux to the stellar envelope, which disposes of much of it in a low velocity wind. Gravitational contraction easily accounts for the momentum flux in the solar wind, as well as the flux required to lift mass into the dust formation zone of every AGB star, whereon radiation pressure continues its ejection as a low velocity wind. This mechanism explains the dependence of the AGB mass-loss rate on core mass; its generalization to objects with angular momentum and/or strong magnetic fields suggests a novel explanation of why most planetary nebulae and proto planetary nebulae exhibit axial symmetry. Quasistatic contraction is inherently biased to the generation of the maximum possible momentum flux. Its formalism is therefore readily adapted to providing an upper limit to the momentum flux needed to sustain mass loss when this begins from a semi-continuous rather than impulsive process.Comment: 35 pages, including 1 fig and 2 tables, to appear in Astrophysical Journal -- ps documen

Near and mid-IR sub-arcsecond structure of the dusty symbiotic star R Aqr

The results of a high-resolution interferometric campaign targeting the symbiotic long-period variable (LPV) R~Aqr are reported. With both near-infrared measurements on baselines out to 10m and mid-infrared data extending to 32m, we have been able to measure the characteristic sizes of regions from the photosphere of the LPV and its extended molecular atmosphere, out to the cooler circumstellar dust shell. The near-infrared data were taken using aperture masking interferometry on the Keck-I telescope and show R~Aqr to be partially resolved for wavelengths out to 2.2 microns but with a marked enlargement, possibly due to molecular opacity, at 3.1 microns. Mid-infrared interferometric measurements were obtained with the U.C. Berkeley Infrared Spatial Interferometer (ISI) operating at 11.15 microns from 1992 to 1999. Although this dataset is somewhat heterogeneous with incomplete coverage of the Fourier plane and sampling of the pulsation cycle, clear changes in the mid-infrared brightness distribution were observed, both as a function of position angle on the sky and as a function of pulsation phase. Spherically symmetric radiative transfer calculations of uniform-outflow dust shell models produce brightness distributions and spectra which partially explain the data, however limitations to this approximation are noted. Evidence for significant deviation from circular symmetry was found in the mid-infrared and more tentatively at 3.08 microns in the near-infrared, however no clear detection of binarity or of non-LPV elements in the symbiotic system is reported.Comment: Accepted to Astrophysical Journal. To appear in volume 534. 14 pages; 3 postscript figure

The last gasps of VY CMa: Aperture synthesis and adaptive optics imagery

We present new observations of the red supergiant VY CMa at 1.25 micron, 1.65 micron, 2.26 micron, 3.08 micron and 4.8 micron. Two complementary observational techniques were utilized: non-redundant aperture masking on the 10-m Keck-I telescope yielding images of the innermost regions at unprecedented resolution, and adaptive optics imaging on the ESO 3.6-m telescope at La Silla attaining extremely high (~10^5) peak-to-noise dynamic range over a wide field. For the first time the inner dust shell has been resolved in the near-infrared to reveal a one-sided extension of circumstellar emission within 0.1" (~15 R_star) of the star. The line-of-sight optical depths of the circumstellar dust shell at 1.65 micron, 2.26 micron, and 3.08 micron have been estimated to be 1.86 +/- 0.42, 0.85 +/- 0.20, and 0.44 +/- 0.11. These new results allow the bolometric luminosity of VY~CMa to be estimated independent of the dust shell geometry, yielding L_star ~ 2x10^5 L_sun. A variety of dust condensations, including a large scattering plume and a bow-shaped dust feature, were observed in the faint, extended nebula up to 4" from the central source. While the origin of the nebulous plume remains uncertain, a geometrical model is developed assuming the plume is produced by radially-driven dust grains forming at a rotating flow insertion point with a rotational period between 1200-4200 years, which is perhaps the stellar rotational period or the orbital period of an unseen companion.Comment: 25 pages total with 1 table and 5 figures. Accepted by Astrophysical Journal (to appear in February 1999

On the difference between type E and A OH/IR stars

The observed SEDs of a sample of 60 OH/IR stars are fitted using a radiative transfer model of a dusty envelope. Among the whole sample, 21 stars have reliable phase-lag distances while the others have less accurate distances. L*-P,Mlr-P and Mlr-L* relations have been plotted for these stars. It is found that type E (with emission feature at 10um and type A (with absorption feature at 10um) OH/IR stars have different L*-P and Mlr-L* relations while both of them follow a single Mlr-P relation. The type E stars are proven to be located in the area without large scale dense interstellar medium while the type A stars are located probably in dense interstellar medium. It is argued here that this may indicate the two types of OH/IR stars have different chemical composition or zero age main sequence mass and so evolve in different ways. This conclusion has reinforced the argument by Chen et al.(2001) who reached a similar conclusion from the galactic distribution of about 1000 OH/IR stars with the IRAS low-resolution spectra (LRS).Comment: 6 pages, 9 figures, 2 table

Evolution of the dust mass loss with luminosity along the giant branch of the globular cluster 47 Tuc

The paper investigates the properties of the dust mass loss in stars populating the giant branch of the globular cluster 47 Tuc, by combining ISOCAM and DENIS data. Raster maps of 5 fields covering areas ranging from 4 x 4 to 15 x 15 arcmin2 at different distances from the center of the cluster have been obtained with ISOCAM at 11.5 mum (LW10 filter). The covered fields include most of the red variables known in this cluster. A detection threshold of about 0.2 mJy is achieved, allowing to detect giant stars at 11.5 mum all the way down to the horizontal branch. No dust-enshrouded asymptotic giant branch stars have been found in the observed fields, contrary to the situation encountered in LMC/SMC globular clusters with larger turnoff masses. The color index [12]-[2] (based on the ISO 11.5 mum flux and on the DENIS Ks magnitude) is used as a diagnostic of dust emission (and hence dust mass loss). Its evolution with luminosity along the giant branch reveals that dust mass loss is only present in V3 (the only cluster Mira variable observed in the present study) and in V18, a star presenting intermittent variability. This conclusion confirms the importance of stellar pulsations in the dust formation and ensuing mass loss.Comment: 16 pages, accepted in Astronomy & Astrophysic

Pareto Autonomous Local Search

This paper presents a study for the dynamic selection of operators in a local search process. The main purpose is to propose a generic autonomous local search method which manages operator selection from a set of available operators, built on neighborhood relations and neighbor selection functions, using the concept of Pareto dominance with respect to quality and diversity. The latter is measured using two different metrics. This control method is implemented using the Comet language in order to be easily introduced in various constraint local search algorithms. Focusing on permutation-based problems, experimental results are provided for the QAP and ATSP to assess the method’s effectiveness

The Physical Properties of the Red Supergiant WOH G64: The Largest Star Known?

WOH G64 is an unusual red supergiant (RSG) in the Large Magellanic Cloud (LMC), with a number of properties that set it apart from the rest of the LMC RSG population, including a thick circumstellar dust torus, an unusually late spectral type, maser activity, and nebular emission lines. Its reported physical properties are also extreme, including the largest radius for any star known and an effective temperature that is much cooler than other RSGs in the LMC, both of which are at variance with stellar evolutionary theory. We fit moderate-resolution optical spectrophotometry of WOH G64 with the MARCS stellar atmosphere models, determining an effective temperature of 3400 +/- 25 K. We obtain a similar result from the star's broadband V - K colors. With this effective temperature, and taking into account the flux contribution from the aysmmetric circumstellar dust envelope, we calculate log(L/L_sun) = 5.45 +/- 0.05 for WOH G64, quite similar to the luminosity reported by Ohnaka and collaborators based on their radiative transfer modeling of the star's dust torus. We determine a radius of R/R_sun = 1540, bringing the size of WOH G64 and its position on the H-R diagram into agreement with the largest known Galactic RSGs, although it is still extreme for the LMC. In addition, we use the Ca II triplet absorption feature to determine a radial velocity of 294 +/- 2 km/s for the star; this is the same radial velocity as the rotating gas in the LMC's disk, which confirms its membership in the LMC and precludes it from being an unusual Galactic halo giant. Finally, we describe the star's unusual nebula emission spectrum; the gas is nitrogen-rich and shock-heated, and displays a radial velocity that is significantly more positive than the star itself by 50 km/s.Comment: 25 pages, 5 figures; accepted for publication in The Astronomical Journa

The Infrared Continuum Spectrum of VY CMa

We combine spectra of VY CMa obtained with the short- and long-wavelength spectrometers, SWS and LWS, on the Infrared Space Observatory to provide a first detailed continuum spectrum of this highly luminous star. The circumstellar dust cloud through which the star is observed is partially self-absorbing, which makes for complex computational modeling. We review previous work and comment on the range of uncertainties about the physical traits and mineralogical composition of the modeled disk. We show that these uncertainties significantly affect the modeling of the outflow and the estimated mass loss. In particular, we demonstrate that a variety of quite diverse models can produce good fits to the observed spectrum. If the outflow is steady, and the radiative repulsion on the dust cloud dominates the star's gravitational attraction, we show that the total dust mass-loss rate is $\sim 4\times 10^{-6}M_{\odot}$ yr$^{-1}$, assuming that the star is at a distance of 1.5 kpc. Several indications, however, suggest that the outflow from the star may be spasmodic. We discuss this and other problems facing the construction of a physically coherent model of the dust cloud and a realistic mass-loss analysis