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.     

Observable effects of dust formation in dynamic atmospheres of M-type Mira variables

The formation of dust with temperature-dependent non-grey opacity is considered in a series of self-consistent model atmospheres at different phases of an O-rich Mira variable of mass  1.2 M_⊙  . Photometric and interferometric properties of these models are predicted under different physical assumptions regarding the dust formation. The iron content of the initial silicate that forms and the availability of grain nuclei are found to be critical parameters that affect the observable properties. For certain plausible parameter values, dust would form at 2–3 times the average continuum photospheric radius. This work provides a consistent physical explanation for the larger apparent size of Mira variables at wavelengths shorter than 1 μm than that predicted by dust-free fundamental-mode pulsation models.     

Determining temperature amplitude as a function of depth in the atmospheres of rapidly oscillating Ap stars

We present numerical models based on realistic treatment of the intensity spectrum (from model atmospheres), and demonstrate that they are consistent with Kurtz and Medupe's recent formula in showing that limb darkening is too small an effect to explain the observed sharp decline of pulsation light amplitude with wavelength in rapidly oscillating Ap stars. Kurtz and Medupe's formula is shown to be a special form of Watson's earlier general formula for non-radial light variations of a star pulsating in any mode ( l m ). Using a technique suggested by Kurtz and Medupe we derive temperature semi-amplitude as a function of depth in the atmospheres of α Cir and HR 3831, assuming that we can neglect non-adiabatic effects.     

Near Infrared Observations of Sakurai's Object during 1998–99

We present new near-infrared observations of Sakurai's Object obtainedduring 1998–99 when this final helium shell flash object was in the dustcondensation phase. The infrared colours have reddened compared to earlierepochs, indicating increased dust condensation. The infrared spectrareveal all the features of a carbon star superposed on a dust continuum.     

SCUBA photometry of candidate Vega-like sources

New SCUBA measurements at millimetre wavelengths are presented for a sample of Vega-like stars. Six stars were detected, while sensitive upper limits were obtained for a further 11 sources. Most of the sample selected from a recent catalogue of Vega-like stars have infrared excesses similar to those of the prototype Vega-like stars α Lyr and α PsA. Their IR–submm spectral indices are steep, indicating that the submm emission from the discs is dominated by grains which are smaller than the wavelength of observation and that only small grains exist in those dusty discs. HD 98800 has an IR–submillimetre spectral index of less than two, which suggests that grains have grown to more than 0.3 mm in size. Hipparcos parallax data for HD 42137 and HD 123160 suggest that these two stars are giants rather than dwarfs, similar to the situation previously found for HD 233517. Dust masses, or upper limits, were derived for the sample; these indicate that most of the sources do not have as much dust as Herbig Ae/Be or T Tauri stars, but are likely to have dust masses comparable to those of the prototype Vega-like stars.     

Mid-infrared imaging of WL 16: the spatial distribution of the hydrocarbon emission features

We present arcsec-resolution images at 8.2, 10.0 and 11.3 μm of the unusual young object WL 16 in Ophiuchus, which has an extended envelope of fluorescing hydrocarbon molecules. To the limit of achieved sensitivity, the faint 10.0-μm continuum has a surface-brightness distribution that is not distinguishable from those at 8.2 and 11.3 μm, where the luminosity is known to be dominated by the polycyclic aromatic hydrocarbon (PAH) emission features. We conclude that the 10-μm continuum either arises from non-equilibrium heating of small dust grains that are well mixed with the hydrocarbons or is quasi-continuous emission from the PAH particles themselves, rather than thermal equilibrium emission from macroscopic dust grains, and that there is no significant silicate absorption variation across the source. The extended hydrocarbon emission may trace a flattened, equatorial distribution of circumstellar material or arise in bipolar lobes. The former case is slightly favoured, based on currently available data, and would imply that WL 16 is a relatively evolved Herbig Ae star, the equatorial plane of which has been almost cleared of normal dust, leaving only fluorescing hydrocarbons and larger coagulated particles as a possibly transient fossil of the original circumstellar disc.     

The accurate determination of stellar angular diameters using broad-band stellar interferometry

The angular diameter of a star can be estimated from interferometric observations by fitting the data with the visibility function for a uniformly illuminated disc and then using published correction factors to convert the uniform-disc angular diameter to the limb-darkened angular diameter. The correction factors are strictly valid only for monochromatic light. We investigate the effect of using a broad bandwidth, and present a simple method for calculating broad-band correction factors from the monochromatic factors.
The technique of fitting the data with a uniform-disc visibility function is only useful for stars with compact atmospheres and 'typical' limb-darkening profiles. It should not be applied to stars with extended atmospheres or that show extreme limb darkening. These stars have visibility functions that are qualitatively different from a uniform-disc visibility function, so they can be distinguished observationally from compact-atmosphere 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.     

Center-to-limb variation of the continuum intensity and linear polarization of stars with transiting exoplanets

The limb darkening and center-to-limb variation of the continuum polarization is calculated for a grid of one-dimensional stellar model atmospheres and for a wavelength range between 300 and 950 nm. Model parameters match those of the transiting stars taken from the NASA exoplanet archive. The limb darkening of the continuum radiation for these stars is shown to decrease with the rise in their effective temperature. For the λ = 370 nm wavelength, which corresponds to the maximum of the Johnson–Cousins UX filter, the limb darkening values of the planet transiting stars lie in a range between 0.03 and 0.3. The continuum linear polarization depends not only on the effective temperature of the star but also on its gravity and metallicity. Its value decreases for increasing values of these parameters. In the UX band, the maximum linear polarization of stars with transiting planets amounts to 4%, while the minimum value is approximately 0.3%. The continuum limb darkening and the linear polarization decrease rapidly with wavelength. At the R band maximum (λ = 700 nm), the linear polarization close to the limb is in fact two orders of magnitude smaller than in the UX band. The center- to-limb variation of the continuum intensity and the linear polarization of the stars with transiting planets can be approximated, respectively, by polynomials of the fourth and the sixth degree. The coefficients of the polynomials, as well as the IDL procedures for reading them, are available in electronic form. It is shown that there are two classes of stars with high linear polarization at the limb. The first one consists of cold dwarfs. Their typical representatives are HATS-6, Kepler-45, as well as all the stars with similar parameters. The second class of stars includes hotter giants and subgiants. Among them we have CoRoT-28, Kepler-91, and the group of stars with effective temperatures and gravities of approximately 5000 K and 3.5, respectively.     

Too large and overlooked? Extended free–free emission towards massive star formation regions

We present Australia Telescope Compact Array observations towards six massive star formation regions, which, from their strong 24 GHz continuum emission but no compact 8 GHz continuum emission, appeared good candidates for hypercompact H  ii regions. However, the properties of the ionized gas derived from the 19 to 93 GHz continuum emission and  H70α+ H57α  radio recombination line data show the majority of these sources are, in fact, regions of spatially extended, optically thin free–free emission. These extended sources were missed in the previous 8 GHz observations due to a combination of spatial filtering, poor surface brightness sensitivity and primary beam attenuation.
We consider the implications that a significant number of these extended H  ii regions may have been missed by previous surveys of massive star formation regions. If the original sample of 21 sources is representative of the population as a whole, the fact that six contain previously undetected extended free–free emission suggests a large number of regions have been mis-classified. Rather than being very young objects prior to UCH  ii region formation, they are, in fact, associated with extended H  ii regions and thus significantly older. In addition, inadvertently ignoring a potentially substantial flux contribution (up to ∼0.5 Jy) from free–free emission has implications for dust masses derived from sub-mm flux densities. The large spatial scales probed by single-dish telescopes, which do not suffer from spatial filtering, are particularly susceptible and dust masses may be overestimated by up to a factor of ∼2.     

Dust coagulation in star formation with different metallicities

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution and possibly affects the thermal evolution of star-forming clouds. We here investigate dust coagulation in collapsing pre-stellar cores with different metallicities by considering the thermal motions of grains. We show that coagulation does occur even at low metallicity  ∼10⁻⁶ Z_⊙  . However, we also find (i) that the H₂ formation rate on dust grains is reduced only after the majority of H₂ is formed and (ii) that the dust opacity is modified only after the core becomes optically thick. Therefore, we conclude that the effects of dust coagulation can safely be neglected in discussing the temperature evolution of the pre-stellar cores for any metallicity as long as the grain motions are thermal.     

Hubble Space Telescope observations of SV Cam – II. First derivative light-curve modelling using phoenix and atlas model atmospheres

The variation of the specific intensity across the stellar disc is an essential input parameter in surface brightness reconstruction techniques such as Doppler imaging, where the relative intensity contributions of different surface elements are important in detecting star-spots. We use phoenix and atlas model atmospheres to model light curves derived from high precision (signal-to-noise ratio ≃ 5000) Hubble Space Telescope ( HST ) data of the eclipsing binary SV Cam (F9V+K4V), where the variation of specific intensity across the stellar disc will determine the contact points of the binary system light curve. For the first time, we use χ² comparison fits to the first derivative profiles to determine the best-fitting model atmosphere. We show the wavelength dependence of the limb darkening and that the first derivative profile is sensitive to the limb-darkening profile very close to the limb of the primary star. It is concluded that there is only a marginal difference (<1σ) between the χ² comparison fits of the two model atmospheres to the HST light curve at all wavelengths. The usefulness of the second derivative of the light curve for measuring the sharpness of the primary's limb is investigated, but we find that the data are too noisy to permit a quantitative analysis.     

IRAS 12311−3509: a carbon star with SiC2 emission

The optical spectrum of the carbon star IRAS 12311−3509 is dominated by the Merrill–Sanford emission bands of SiC₂, by absorption and emission in the Swan system of C₂, and by resonance emission lines of neutral metals. The infrared energy distribution is flat from 1 to 60 μm. These observations are interpreted as arising from a star with a cool dusty disc which is edge-on to the observer and obscures direct starlight. The infrared continuum is caused predominantly by absorption of stellar light by dust in the disc and re-emission at longer wavelengths. The optical stellar spectrum is seen by reflection off dusty material which lies out of the plane of the disc, and the molecular and atomic emission arises in the same geometry through resonance fluorescence. The object has similarities to the J-silicate stars, but may have a carbon-rich rather than oxygen-rich disc. A full spectroscopic assignment and discussion of the SiC₂ bands and their intensities are given. Modelling of the rotational contours of the     band yields a rotational temperature of 250 K, indicating very cool gas.     

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.     

2002 Kuiper prize lecture: Dust Astronomy

Dust particles, like photons, carry information from remote sites in space and time. From knowledge of the dust particles' birthplace and their bulk properties, we can learn about the remote environment out of which the particles were formed. This approach is called “Dust Astronomy” which is carried out by means of a dust telescope on a Dust Observatory in space. Targets for a dust telescope are the local interstellar medium and nearby star forming regions, as well as comets and asteroids. Dust from interstellar and interplanetary sources is distinguished by accurately sensing their trajectories. Trajectory sensors may use the electric charge signals that are induced when charged grains fly through the detector. Modern in-situ dust impact detectors are capable of providing mass, speed, physical and chemical information of dust grains in space. A Dust Observatory mission is feasible with state-of-the-art technology. It will (1) provide the distinction between interstellar dust and interplanetary dust of cometary and asteroidal origin, (2) determine the elemental composition of impacting dust particles, and (3) monitor the fluxes of various dust components as a function of direction and particle masses.     

Limb-darkening corrections for interferometric uniform disc stellar angular diameters

Stellar angular diameters determined interferometrically are generally established by fitting the observed visibility data with a curve appropriate for a uniformly illuminated disc. The resulting uniform-disc diameters must be corrected for the effects of limb darkening in order to determine the true angular diameters of the stars. An extensive grid of limb-darkening corrections, based directly on the centre-to-limb intensity variations for Kurucz model stellar atmospheres, has been computed without the intermediate step of a parametrized representation of the centre-to-limb variation. The limitations of this method of correction are discussed.     

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.     

Hα emission objects in Circinus

We present a new survey for Hα emission objects in the Circinus cloud complex and introduce an efficient photometric method for detecting Hα emission via observations in a narrow‐band filter. The observed flux is compared to a blackbody fit of the continuum. Our search strategy reveals 20 stars with strong Hα emission (EW > 10 Å), eight of them being new detections. All Hα stars display infrared excess corroborating their youth. On the other hand, the region contains a number of infrared excess objects that do not show Hα emission. Our results support the picture that accretion – as witnessed by Hα emission – is a highly variable phenomenon. Therefore, Hα surveys can only trace the temporarily active objects. In contrast, infrared excess is a more robust tracer that reveals most of the population of young stellar objects in a star forming region. Our analysis of the general stellar content of the region yields a reliable distance of 450 pc for the Circinus cloud. Moreover, we find a ratio of total‐to‐selective extinction of R^V = 2.8 suggesting that smaller‐than‐normal dust grains may be present. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Infrared spectroscopy of Nova Cassiopeiae 1993 – IV. A closer look at the dust

Nova Cassiopeiae 1993 (V705 Cas) was an archetypical dust-forming nova. It displayed a deep minimum in the visual light curve, and spectroscopic evidence for carbon, hydrocarbon and silicate dust. We report the results of fitting the infrared (IR) spectral energy distribution (SED) with the dusty code, which we use to determine the properties and geometry of the emitting dust. The emission is well described as originating in a thin shell whose dust has a carbon:silicate ratio of 2:1 by number (  ∼1.26:1  by mass) and a relatively flat size distribution. The 9.7- and 18-μm silicate features are consistent with freshly condensed dust and, while the lower limit to the grain size distribution is not well constrained, the largest grains have dimensions  ∼0.06 μm  ; unless the grains in V705 Cas were anomalously small, the sizes of grains produced in nova eruptions may previously have been overestimated in novae with optically thick dust shells. Laboratory work by Grishko & Duley may provide clues to the apparently unique nature of nova unidentified infrared (UIR) features.     

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.