100-yr mass-loss modulations on the asymptotic giant branch

We analyse the differences in infrared circumstellar dust emission between oxygen-rich Mira and non-Mira stars, and find that they are statistically significant. In particular, we find that these stars segregate in the K–[12] versus [12]–[25] colour–colour diagram, and have distinct properties of the IRAS LRS spectra, including the peak position of the silicate emission feature. We show that the infrared emission from the majority of non-Mira stars cannot be explained within the context of standard steady-state outflow models.
The models can be altered to fit the data for non-Mira stars by postulating non-standard optical properties for silicate grains, or by assuming that the dust temperature at the inner envelope radius is significantly lower (300–400 K) than typical silicate grain condensation temperatures (800–1000 K) . We argue that the latter is more probable and provide detailed model fits to the IRAS LRS spectra for 342 stars. These fits imply that two-thirds of non-Mira stars and one-third of Mira stars do not have hot dust (>500 K) in their envelopes.
The absence of hot dust can be interpreted as a recent (∼100 yr) decrease in the mass-loss rate. The distribution of best-fitting model parameters agrees with this interpretation and strongly suggests that the mass loss resumes on similar time-scales. Such a possibility appears to be supported by a number of spatially resolved observations (e.g. recent Hubble Space Telescope images of the multiple shells in the Egg Nebula) and is consistent with new dynamical models for mass loss on the asymptotic giant branch.     

The influence of atmospheric dust on limb darkening of M-type Mira variables

We explore the occurrence of dust in M-type Mira atmospheres and its effect on limb darkening under schematic assumptions about dust temperatures and dust particle properties. Dust particles that are thermodynamically coupled to the surrounding gas may form and may affect limb darkening, though only by very little in infrared continuum bandpasses. Dust particles that assume the equilibrium temperature given by the mean intensity of the radiation field only form under rare circumstances. Unexpectedly large or wavelength-dependent infrared continuum radii observed by interferometry are unlikely to be caused by atmospheric dust, except possibly near 1 μm; however, radius measurements may be significantly affected by molecular band contamination.     

Scattered light models of the dust shell around HD 179821

The dust shell around the evolved star HD 179821 has been detected in scattered light in near-IR imaging polarimetry observations. Here, we subtract the contribution of the unpolarized stellar light to obtain an intrinsic linear polarization of between 30 and 40 per cent in the shell that seems to increase with radial offset from the star. The J - and K -band data are modelled using a scattering code to determine the shell parameters and dust properties. We find that the observations are well described by a spherically-symmetric distribution of dust with an r ⁻² density law, indicating that when mass-loss was occurring, the mass-loss rate was constant. The models predict that the detached nature of a spherically-symmetric, optically-thin dust shell, with a distinct inner boundary, will only be apparent in polarized flux. This is in accordance with the observations of this and other optically-thin circumstellar shells, such as IRAS 17436+5003. By fitting the shell brightness we derive an optical depth to the star that is consistent with V -band observations and that, assuming a distance of 6 kpc, gives an inner-shell radius of     , a dust number density of     at r _in and a dust mass of     . We have explored axisymmetric shell models but conclude that any deviations from spherical symmetry in the shell must be slight, with an equator-to-pole density contrast of less than 2:1. We have not been able to fit simultaneously the high linear polarizations and the small     colour excess of the shell and we attribute this to the unusual scattering properties of the dust. We suggest that the dust grains around HD 179821 either are highly elongated or consist of aggregates of smaller particles.     

Cyclic variations in the angular diameter of χ Cygni

We report the direct detection of cyclic diameter variations in the Mira variable χ Cygni. Interferometric observations made between 1997 July and 1999 September, using the Cambridge Optical Aperture Synthesis Telescope (COAST) and the William Herschel Telescope (WHT), indicate periodic changes in the apparent angular diameter at a wavelength of 905 nm, with amplitude 45 per cent of the smallest value. The star appears largest at minimum light. Measurements made at a wavelength of 1.3 μm over the same period suggest much smaller size changes. This behaviour is consistent with a model in which most of the apparent diameter variation at 905 nm is caused by a large increase in the opacity of the outer atmospheric layers (which is mostly owing to titanium oxide) near minimum light, rather than by physical motions of the photosphere. The 1.3-μm waveband is relatively uncontaminated by TiO, and so much smaller size changes would be expected in this band. The latest non-linear pulsational models predict maximum physical size close to maximum light, and increases in opacity near minimum light that are too small to reproduce the diameter variation seen at 905 nm. This suggests either that the phase-dependence of the model pulsation is incorrect, or that the opacities in the models are underestimated. Future interferometric monitoring in uncontaminated near-infrared wavebands should resolve this question.     

Large-amplitude periodic variations in the angular diameter of R Leonis

We report the first direct detection of long-term periodic diameter variations in a Mira variable. Angular diameter measurements of the 313-d period variable R Leonis at 833 nm and 940 nm obtained between 1996 February and 1997 June using the Cambridge Optical Aperture Synthesis Telescope (COAST) and the William Herschel Telescope (WHT) show a cyclic modulation of the apparent stellar diameter by approximately 35 per cent. The agreement between these new data and archival measurements from 1992 January suggests coherence in the modulation over a 5-yr period. Our data are consistent with recent models which suggest that, in photometric bands with only weak to moderate molecular contamination, periodic variations in stellar diameter of order 50 per cent can be maintained. The measurements indicate that the apparent stellar diameter was largest at visual phase 0.5 and that any phase shifts between the visual light curve and those at 833 and 940 nm were at most 0.05. The large offset (∼ 0.25) between the phase of the observed diameter maximum and that predicted for the photospheric continuum diameter variations suggests that our observations are more sensitive to the changing temperature structure of the outer atmosphere than to the deeper continuum-forming layers.     

Mode switching in the nearby Mira-like variable R Doradus

We discuss visual observations spanning nearly 70 years of the nearby semiregular variable R Doradus. Using wavelet analysis, we show that the star switches back and forth between two pulsation modes having periods of 332 d and about 175 d, the latter with much smaller amplitude. Comparison with model calculations suggests that the two modes are the first and third radial overtone, with the physical diameter of the star making fundamental-mode pulsation unlikely. The mode changes occur on a time-scale of about 1000 d, which is too rapid to be related to a change in the overall thermal structure of the star and may instead be related to weak chaos.   The Hipparcos distance to R Dor is 62.4 ± 2.8 pc which, taken with its dominant 332-d period, places it exactly on the period–luminosity (P–L) relation of Miras in the Large Magellanic Cloud. Our results imply first-overtone pulsation for all Miras which fall on the P–L relation. We argue that semiregular variables with long periods may largely be a subset of Miras and should be included in studies of Mira behaviour. The semiregulars may contain the immediate evolutionary Mira progenitors, or stars may alternate between periods of semiregular and Mira behaviour.     

The variability of Magellanic cluster infrared stars

We report on the light variations of the infrared stars that were discovered recently in the Magellanic clusters NGC 419, 1783 and 1978. Their periods, of 528, 458 and 491 days, are among the longest known for carbon-rich Mira variables in the Clouds. All three IR stars were found to lie on the extension of the period– M _bol relation derived from the shorter-period C-rich Miras while they were 0.45–0.70 mag fainter than the extension of the period– M _K relation. Their main sequence masses were determined by isochrone fitting to be 1.5–1.6 M_⊙, consistent with the prediction of the evolutionary models of Vassiliadis & Wood.     

Wind-swept clouds and possible triggered star formation associated with the supernova remnant G357.7+0.3

We present evidence for interaction between the supernova remnant (SNR) G357.7+0.3 and nearby molecular clouds, leading to the formation of wind-swept structures and bright emission rims. These features are not observed at visual wavelengths, but are clearly visible in mid-infrared mapping undertaken using the Spitzer Space Telescope . Analysis of one of these clouds, the bright cometary structure G357.46+0.60, suggests that it contains strong polycyclic aromatic hydrocarbon emission features in the 5.8 and 8.0 μm photometric bands, and that these are highly variable over relatively small spatial scales. The source is also associated with strong variations in electron density; a far-infrared continuum peak associated with dust temperatures of ∼30 K; and has previously been observed in the 1720 MHz maser transition of OH, known to be associated with SNR shock excitation of interstellar clouds. This source also appears to contain a young stellar object (YSO) within the bright rim structure, with a steeply rising spectrum between 1.25 and 24 μm. If the formation of this star has been triggered recently by the SNR, then YSO modelling suggests a stellar mass  ∼5–10 M_⊙  , and luminosity   L _YSO∼10²–2 × 10³ L_⊙  .
Finally, it is noted that a further, conical emission region appears to be associated with the Mira V1139 Sco, and it is suggested that this may represent the case of a Mira outflow interacting with a SNR. If this is the case, however, then the distance to the SNR must be ∼half of that determined from CS   J = 2–1  and 3–2 line radial velocities.     

Quantifying dust and the ultraviolet radiation density in the local Universe

A sample of local galaxies for which far-infrared and ultraviolet fluxes are available is used to estimate the characteristic dust extinction in galaxies and to test whether standard dust properties are plausible. Assuming galaxies can be characterized by a single dust optical depth (certainly not valid for galaxies with a dominant starburst component), the infrared excess and ultraviolet colours of local galaxies are found to be consistent with normal Milky Way dust, with a mean value for   E ( B − V )  of 0.16. A significant fraction of the dust heating is caused by older, lower-mass stars, and this fraction increases towards earlier galaxy types.
Analysis of  ( F _FIR/ F _UV)  versus ultraviolet colour diagrams for starburst galaxies in terms of a simple screen dust model does not support a Calzetti (1997) rather than a Milky Way extinction law, though the absence of the expected strong 2200- Å feature in several galaxies with IUE spectra does show that more detailed radiative transfer models are needed – probably with nonspherical geometry.
A simple treatment in which the  100/60-μm  flux ratio is used to subtract the optically thick starburst contribution to the far-infrared radiation results in lower extinction estimates for the optically thin cirrus component, with a mean   E ( B − V )  of 0.10.
The ultraviolet luminosity density, corrected for dust extinction, is derived and a value for the local mean star formation rate inferred. This is consistent with previous estimates from ultraviolet surveys and from  Hα  surveys.     

Continuum emission around AGB stars at 1.2 mm

It is generally acknowledged that the mass-loss of asymptotic giant branch (AGB) stars undergoes variations on different time-scales. We constructed models for the dust envelopes for a sample of AGB stars to assess whether mass-loss variations influence the spectral energy distribution. To constrain the variability, extra observations at millimetre wavelengths (1.2 mm) were acquired. From the analysis of the dust models, two indications for the presence of mass-loss variations can be found, being (1) a dust temperature at the inner boundary of the dust envelope that is far below the dust condensation temperature and (2) an altered density distribution with respect to  ρ( r ) ∝ r ⁻²  resulting from a constant mass-loss rate. For five out of the 18 studied sources a two-component model of the envelope is required, consisting of an inner region with a constant mass-loss rate and an outer region with a less steep density distribution. For one source an outer region with a steeper density distribution was found. Moreover, in a search for time variability in our data set at 1.2 mm, we found that WX Psc shows a large relative time variation of 34 per cent which might partially be caused by variable molecular line emission.     

SiO maser emission in Miras

We describe a combined dynamic atmosphere and maser propagation model of SiO maser emission in Mira variables. This model rectifies many of the defects of an earlier model of this type, particularly in relation to the infrared (IR) radiation field generated by dust and various wavelength-dependent, optically thick layers. Modelled masers form in rings with radii consistent with those found in very long baseline interferometry (VLBI) observations and with earlier models. This agreement requires the adoption of a radio photosphere of radius approximately twice that of the stellar photosphere, in agreement with observations. A radio photosphere of this size renders invisible certain maser sites with high amplification at low radii, and conceals high-velocity shocks, which are absent in radio continuum observations. The SiO masers are brightest at an optical phase of 0.1–0.25, which is consistent with observed phase lags. Dust can have both mild and profound effects on the maser emission. Maser rings, a shock and the optically thick layer in the SiO pumping band at 8.13 μm appear to be closely associated in three out of four phase samples.     

The impact of meridional circulation on stellar butterfly diagrams and polar caps

Observations of rapidly rotating solar-like stars show a significant mixture of opposite-polarity magnetic fields within their polar regions. To explain these observations, models describing the surface transport of magnetic flux demand the presence of fast meridional flows. Here, we link subsurface and surface magnetic flux transport simulations to investigate (i) the impact of meridional circulations with peak velocities of  ≤125 m s⁻¹  on the latitudinal eruption pattern of magnetic flux tubes and (ii) the influence of the resulting butterfly diagrams on polar magnetic field properties. Prior to their eruption, magnetic flux tubes with low field strengths and initial cross-sections below  ∼300 km  experience an enhanced poleward deflection through meridional flows (assumed to be polewards at the top of the convection zone and equatorwards at the bottom). In particular, flux tubes which originate between low and intermediate latitudes within the convective overshoot region are strongly affected. This latitude-dependent poleward deflection of erupting magnetic flux renders the wings of stellar butterfly diagrams distinctively convex. The subsequent evolution of the surface magnetic field shows that the increased number of newly emerging bipoles at higher latitudes promotes the intermingling of opposite polarities of polar magnetic fields. The associated magnetic flux densities are about 20 per cent higher than in the case disregarding the pre-eruptive deflection, which eases the necessity for fast meridional flows predicted by previous investigations. In order to reproduce the observed polar field properties, the rate of the meridional circulation has to be of the order of 100 m s⁻¹, and the latitudinal range from which magnetic flux tubes originate at the base of the convective zone (≲50°) must be larger than in the solar case (≲35°).     

Dust and hydrogen molecules in the extremely metal-poor dwarf galaxy SBS 0335–052

During the early stages of galaxy evolution, the metallicity is generally low and nearby metal-poor star-forming galaxies may provide templates for primordial star formation. In particular, the dust content of such objects is of great importance, because early molecular formation can take place on grains. To gain insight into primeval galaxies at high redshift, we examine the dust content of the nearby extremely low-metallicity galaxy SBS  0335–052  which hosts a very young starburst (≲10⁷ yr). In young galaxies, the dust formation rate in Type II supernovae governs the amount of dust, and by incorporating recent results on dust production in Type II supernovae we model the evolution of dust content. If the star-forming region is compact (≲100 pc), as suggested by observations of SBS  0335–052  , our models consistently explain the quantity of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. We also discuss the H₂ abundance. The compactness of the region is important to H₂ formation, because the optical depth of dust for UV photons becomes large and H₂ dissociation is suppressed. We finally focus on implications for damped Ly α systems.     

Non-emission-line young stars of intermediate mass

We present optical spectra of four intermediate-mass candidate young stellar objects that have often been classified as Herbig Ae/Be stars. Typical Herbig Ae/Be emission features are not present in the spectra of these stars. Three of them, HD 36917, HD 36982 and HD 37062, are members of the young Orion nebula cluster (ONC). This association constrains their ages to be ≲1 Myr. The lack of appreciable near-infrared excess in them suggests the absence of hot dust close to the central star. However, they do possess significant amounts of cold and extended dust as revealed by the large excess emission observed at far-infrared wavelengths. The fractional infrared luminosities  ( L _ir/ L _★)  and the dust masses computed from IRAS fluxes are systematically lower than those found for Herbig Ae/Be stars but higher than those for Vega-like stars. These stars may thus represent the youngest examples of the Vega phenomenon known so far. In contrast, the other star in our sample, HD 58647, is more likely to be a classical Be star, as is evident from the low   L _ir/ L _★  , the scarcity of circumstellar dust, the low polarization, the presence of H α emission and near-infrared excess, and the far-infrared spectral energy distribution consistent with free–free emission similar to other well-known classical Be stars.     

Episodic dust formation by HD 192641 (WR 137) – II

We present new infrared photometry of the WC7-type Wolf–Rayet star HD 192641 (WR 137) from 1985 to 1999. These data track the cooling of the dust cloud formed in the 1982–84 dust-formation episode from 1985 to 1991, the increase of the infrared flux from 1994.5 to a new dust-formation maximum in 1997 and its subsequent fading. From these and earlier data we derive a period of 4765±50 d (13.05±0.15 yr) for the dust-formation episodes. Between dust-emission episodes, the infrared spectral energy distribution has the form of a power law, λF _λ ∝ λ ^−1.86. The rising branch of the infrared light curve (1994–97) differs in form from that of the episodic dust-maker WR 125. Time-dependent modelling shows that this difference can be attributed to a different time dependence of dust formation in WR 137, which occurred approximately ∝ t ² until maximum, whereas that of WR 125 could be described by a step function, akin to a threshold effect. For an adopted distance of 1.6 kpc, the rate of dust formation was found to be 5.0×10⁻⁸ M_⊙ yr⁻¹ at maximum, accounting for a fraction f _C≈1.5×10⁻³ of the carbon flowing in the stellar wind. The fading branches of the light curves show evidence for secondary 'mini-eruptions' in 1987, 1988 and 1990, behaviour very different from that of the prototypical episodic dust-maker HD 193793 (WR 140), and suggesting the presence in the WR 137 stellar wind of large-scale structures that are crossed by the wind–wind collision region.     

XMM–Newton X-ray observations of the Wolf–Rayet binary system WR 147

We present results of an ≈20-ks X-ray observation of the Wolf–Rayet (WR) binary system WR 147 obtained with XMM–Newton . Previous studies have shown that this system consists of a nitrogen-type WN8 star plus an OB companion whose winds are interacting to produce a colliding wind shock. X-ray spectra from the pn and MOS detectors confirm the high extinction reported from infrared studies and reveal hot plasma including the first detection of the Fe Kα line complex at 6.67 keV. Spectral fits with a constant-temperature plane-parallel shock model give a shock temperature   kT _shock= 2.7  keV (   T _shock≈ 31  MK), close to but slightly hotter than the maximum temperature predicted for a colliding wind shock. Optically thin plasma models suggest even higher temperatures, which are not yet ruled out. The X-ray spectra are harder than can be accounted for using 2D numerical colliding wind shock models based on nominal mass-loss parameters. Possible explanations include: (i) underestimates of the terminal wind speeds or wind abundances, (ii) overly simplistic colliding wind models or (iii) the presence of other X-ray emission mechanisms besides colliding wind shocks. Further improvement of the numerical models to include potentially important physics such as non-equilibrium ionization will be needed to rigorously test the colliding wind interpretation.     

The STAGES view of red spirals and dusty red galaxies: mass-dependent quenching of star formation in cluster infall

We investigate the properties of optically passive spirals and dusty red galaxies in the A901/2 cluster complex at redshift ∼0.17 using rest-frame near-ultraviolet–optical spectral energy distributions, 24-μm infrared data and Hubble Space Telescope morphologies from the STAGES data set. The cluster sample is based on COMBO-17 redshifts with an rms precision of  σ _cz ≈ 2000 km s⁻¹  . We find that 'dusty red galaxies' and 'optically passive spirals' in A901/2 are largely the same phenomenon, and that they form stars at a substantial rate, which is only four times lower than that in blue spirals at fixed mass. This star formation is more obscured than in blue galaxies and its optical signatures are weak. They appear predominantly in the stellar mass range of  log  M _*/M_⊙=[10, 11]  where they constitute over half of the star-forming galaxies in the cluster; they are thus a vital ingredient for understanding the overall picture of star formation quenching in clusters. We find that the mean specific star formation rate (SFR) of star-forming galaxies in the cluster is clearly lower than in the field, in contrast to the specific SFR properties of blue galaxies alone, which appear similar in cluster and field. Such a rich red spiral population is best explained if quenching is a slow process and morphological transformation is delayed even more. At  log  M _*/M_⊙ < 10  , such galaxies are rare, suggesting that their quenching is fast and accompanied by morphological change. We note that edge-on spirals play a minor role; despite being dust reddened they form only a small fraction of spirals independent of environment.     

The impact of dust on the scaling properties of galaxy clusters

We investigate the effect of dust on the scaling properties of galaxy clusters based on hydrodynamic N -body simulations of structure formation. We have simulated five dust models plus radiative cooling and adiabatic models using the same initial conditions for all runs. The numerical implementation of dust was based on the analytical computations of Montier & Giard. We set up dust simulations to cover different combinations of dust parameters that make evident the effects of size and abundance of dust grains. Comparing our radiative plus dust cooling runs with a purely radiative cooling simulation, we find that dust has an impact on cluster scaling relations. It mainly affects the normalization of the scalings (and their evolution), whereas it introduces no significant differences in their slopes. The strength of the effect critically depends on the dust abundance and grain size parameters as well as on the cluster scaling. Indeed, cooling due to dust is effective in the cluster regime and has a stronger effect on the 'baryon driven' statistical properties of clusters such as   L _X– M , Y – M , S – M   scaling relations. Major differences, relative to the radiative cooling model, are as high as 25 per cent for the   L _X– M   normalization, and about 10 per cent for the Y – M and S – M normalizations at redshift zero. On the other hand, we find that dust has almost no impact on the 'dark matter driven'   T _mw– M   scaling relation. The effects are found to be dependent in equal parts on both dust abundances and grain size distributions for the scalings investigated in this paper. Higher dust abundances and smaller grain sizes cause larger departures from the radiative cooling (i.e. with no dust) model.     

Asymptotic giant branch variables in Baade's Windows

In this work, a sample of luminous M-type giants in the Baade's Windows towards the inner Galactic bulge is investigated in the near-infrared. The ISOGAL survey at 7 and 15 μm has given information concerning the mass-loss rates of these stars and their variability characteristics have been extracted from the MACHO data base. Most are known to be semi-regular variables (SRVs). Here we discuss how their IJHK _S -region colours depend on period and on the presence or absence of mass loss, using results mainly taken from the DENIS and 2MASS surveys.
In order to compare their colours with solar neighbourhood stars, photometric colours on the DENIS, 2MASS and ESO photometric systems have been synthesized for objects in the spectrophotometric atlas of Lançon & Wood. In addition, they have been used to predict the differences in colour indices when stars with strong molecular bands are observed using different photometric systems.
The SRVs are found to inhabit the upper end of the     , K _S colour–magnitude diagram, lying just below the Miras. High mass-loss rates are associated with high luminosity. The near-infrared colours of the SRVs increase in a general way with period and are the reddest for stars with significant mass loss. The average colours of Mira variables, whose periods start at around 200 d in the bulge, are bluer than those of the semi-regulars at this period, particularly in     , thanks to the association of deep water vapour bands with large amplitude.     

Long Secondary Periods in variable red giants

We present a study of a sample of Large Magellanic Cloud red giants exhibiting Long Secondary Periods (LSPs). We use radial velocities obtained from VLT spectral observations and MACHO and OGLE light curves to examine properties of the stars and to evaluate models for the cause of LSPs. This sample is much larger than the combined previous studies of Hinkle et al. and Wood, Olivier & Kawaler.
Binary and pulsation models have enjoyed much support in recent years. Assuming stellar pulsation, we calculate from the velocity curves that the typical fractional radius change over an LSP cycle is greater than 30 per cent. This should lead to large changes in T _eff that are not observed. Also, the small light amplitude of these stars seems inconsistent with the radius amplitude. We conclude that pulsation is not a likely explanation for the LSPs. The main alternative, physical movement of the star – binary motion – also has severe problems. If the velocity variations are due to binary motion, the distribution of the angle of periastron in our large sample of stars has a probability of  1.4 × 10⁻³  that it comes from randomly aligned binary orbits. In addition, we calculate a typical companion mass of  0.09 M_⊙  . Less than 1 per cent of low-mass main-sequence stars have companions near this mass  (0.06–0.12 M_⊙)  whereas ∼25–50 per cent of low-mass red giants end up with LSPs. We are unable to find a suitable model for the LSPs and conclude by listing their known properties.