Overtones of silicate and aluminate minerals and the 5–8 μm ice bands of deeply embedded objects

The distinct patterns, relatively low intensities and peak positions of overtone-combination bands of silicates and oxides suggest that the 5–8 μm spectral region can provide clues for the dust composition when near optically thick conditions exist for the 10-μm silicate feature. We present 1000–2500 cm⁻¹ room-temperature laboratory spectra obtained from powders of silicate, aluminate and nitride minerals and silicate glasses. The spectra exhibit overtone absorption bands with mass absorption coefficients ∼100 times weaker than the fundamentals. These data are compared with the 5–8 μm spectra of deeply embedded young stellar objects observed with the Short Wavelength Spectrometer on the Infrared Space Observatory . Fits of the laboratory data to the observations, after subtraction of the 6.0-μm H₂O ice feature and the 6.0-μm feature identified with organic refractory material, indicate that crystalline melilite (a silicate) or metamict hibonite (a radiation-damaged crystalline aluminate) may be responsible for much of the 6.9-μm absorption feature in the observations, with melilite providing the best match. A weaker 6.2-μm absorption in the young stellar object spectra is well matched by the spectra of hydrous crystalline amphibole silicates (actinolite and tremolite). Relative abundances of Si–O in room-temperature amphiboles to low-temperature H₂O ice are in the range 0.46–3.9 and in melilite are in the range 2.5–8.6. No astronomical feature was matched by the overtones of amorphous silicates because these bands are too broad and peak at the wrong wavelength. Hence, this analysis is consistent with the 10-μm features of these objects being due to a mixture of crystalline and amorphous silicates, rather than only amorphous silicates.     

A study of the 2.5–25 μm spectrum of H2O ice

We present new absorbance spectra of the 3-, 6- and 12-μm features of amorphous and crystalline H₂O ice obtained between 10 and 140 K. Three sets of measurements have been made. In series I, the ice film was initially deposited on to a CsI substrate at 10 K and successive spectra were then obtained at intermediate temperatures as the ice was warmed up to 140 K. The second set, series II, comprises spectra for ice films deposited and measured at temperatures between 10 and 140 K. In the third set of measurements, series III, spectra were obtained for an ice film deposited at 140 K and then at intermediate temperatures as the film was cooled down to 10 K. The series I and II results show that the ice undergoes an amorphous-to-crystalline phase transition in the 110–120 K range. The 3- and 12-μm bands have similar trends in full width at half-maximum (FWHM) and opposite peak wavelength shifts. The temperature behaviour of the 6-μm band is different, as no clear phase transition temperature can be discerned from its FWHM and peak wavelength position data. In the series III spectra the peak wavelength positions and FWHM of the three bands remain relatively constant, thus demonstrating the stability of the crystalline phase against thermal cycling. A comparison between the laboratory results and astronomical spectra suggests that the identification of the librational band of H₂O ice in OH 231.8 + 4.2 may be incorrect.     

A shallow though extensive H2 2.122-μm imaging survey of Taurus–Auriga–Perseus – I. NGC 1333, L1455, L1448 and B1

We discuss wide-field near-infrared (near-IR) imaging of the NGC 1333, L1448, L1455 and B1 star-forming regions in Perseus. The observations have been extracted from a much larger narrow-band imaging survey of the Taurus–Auriga–Perseus complex. These H₂ 2.122-μm observations are complemented by broad-band K imaging, mid-IR imaging and photometry from the Spitzer Space Telescope , and published submillimetre CO   J = 3–2  maps of high-velocity molecular outflows. We detect and label 85 H₂ features and associate these with 26 molecular outflows. Three are parsec-scale flows, with a mean flow lobe length exceeding 11.5 arcmin. 37 (44 per cent) of the detected H₂ features are associated with a known Herbig–Haro object, while 72 (46 per cent) of catalogued HH objects are detected in H₂ emission. Embedded Spitzer sources are identified for all but two of the 26 molecular outflows. These candidate outflow sources all have high near-to-mid-IR spectral indices (mean value of  α∼ 1.4  ) as well as red IRAC 3.6–4.5 μm and IRAC/MIPS 4.5–24.0 μm colours: 80 per cent have [3.6]–[4.5] > 1.0 and [4.5]–[24] > 1.5. These criteria – high α and red [4.5]–[24] and [3.6]–[4.5] colours – are powerful discriminants when searching for molecular outflow sources. However, we find no correlation between α and flow length or opening angle, and the outflows appear randomly orientated in each region. The more massive clouds are associated with a greater number of outflows, which suggests that the star formation efficiency is roughly the same in each region.     

Measurements of the 44-μm band of H2O ice deposited on amorphous carbon and amorphous silicate substrates

We present 20–110 µm absorbance spectra of H₂O ice, deposited on amorphous carbon and silicate substrates, obtained over the 10–140 K temperature range. The measurements have been carried out in a manner that simulates the deposition, warming and cooling of H₂O ice mantles on interstellar and circumstellar grains. For H₂O ice films deposited on these substrates we find (i) similar 44-µm-band peak wavelength temperature dependences, (ii) no bandshape differences in the respective spectra, and (iii) a structural phase transition occurring between 120 and 130 K. In comparison with published data obtained using a polyethylene substrate, the 52-µm feature (the longitudinal optical mode) observed in our spectra is less prominent. This suggests the presence of material-dependent substrate effects that can alter the appearance of the H₂O far-infrared spectrum. The crystallization temperature of H₂O ice films deposited on our amorphous silicate substrate is significantly different from that reported by Moore et al. (1994) , who found crystallization temperatures down to < 20 K for ice also deposited on an amorphous silicate substrate. This is attributed to differences in the surface structures of the respective substrates. This may indicate that, at least in the context of laboratory measurements, substrate material composition is not as significant as substrate surface structure.     

Correlated spectral variability in brown dwarfs

Models of brown dwarf atmospheres suggest they exhibit complex physical behaviour. Observations have shown that they are indeed dynamic, displaying small photometric variations over time-scales of hours. Here, I report results of infrared (0.95–1.64 μm) spectrophotometric monitoring of four field L and T dwarfs spanning time-scales of 0.1–5.5 h, the goal being to learn more about the physical nature of this variability. Spectra are analysed differentially with respect to a simultaneously observed reference source in order to remove Earth-atmospheric variations. The variability amplitude detected is typically 2–10 per cent, depending on the source and wavelength. I analyse the data for correlated variations between spectral indices. This approach is more robust than single band or  χ²  analyses, because it does not assume an amplitude for the (often uncertain) noise level (although the significance test still assumes a shape for the noise power spectrum). Three of the four targets show significant evidence for correlated variability. Some of this can be associated with specific features including Fe, FeH, VO and K  i , and there is good evidence for intrinsic variability in H₂O and possibly also CH₄. Yet some of this variability covers a broader spectral range which would be consistent with dust opacity variations. The underlying common cause is plausibly localized temperature or composition fluctuations caused by convection. Looking at the high signal-to-noise ratio stacked spectra, we see many previously identified spectral features of L and T dwarfs, such as K  i , Na  i , FeH, H₂O and CH₄. In particular, we may have detected methane absorption at 1.3–1.4 μm in the L5 dwarf SDSS 0539−0059.     

Spitzer observations of acetylene bands in carbon-rich asymptotic giant branch stars in the Large Magellanic Cloud

We investigate the molecular bands in carbon-rich asymptotic giant branch (AGB) stars in the Large Magellanic Cloud (LMC), using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope ( SST ) over the 5–38 μm range. All 26 low-resolution spectra show acetylene (C₂H₂) bands at 7 and 14 μm. The hydrogen cyanide (HCN) bands at these wavelengths are very weak or absent. This is consistent with low nitrogen abundances in the LMC. The observed 14 μm C₂H₂  band is reasonably reproduced by an excitation temperature of 500 K. There is no clear dilution of the 14 μm C₂H₂  band by circumstellar dust emission. This 14-μm band originates from molecular gas in the circumstellar envelope in these high mass-loss rate stars, in agreement with previous findings for Galactic stars. The C₂H₂ column density, derived from the 13.7 μm band, shows a gas mass-loss rate in the range 3 × 10⁻⁶ to 5 × 10⁻⁵ M_⊙ yr⁻¹. This is comparable with the total mass-loss rate of these stars estimated from the spectral energy distribution. Additionally, we compare the line strengths of the 13.7 μm C₂H₂  band of our LMC sample with those of a Galactic sample. Despite the low metallicity of the LMC, there is no clear difference in the C₂H₂  abundance among LMC and Galactic stars. This reflects the effect of the third dredge-up bringing self-produced carbon to the surface, leading to high carbon-to-oxygen ratio at low metallicity.     

Hyperluminous infrared galaxies

39 galaxies are now known, from follow-up of faint IRAS sources and from submillimetre observations of high-redshift AGN, with far-infrared luminosities >10¹³ L_⊙. 13 of these, which have been found in 60- or 850-μm surveys, form an important unbiased subsample. 12 have been found by comparison of 60-μm surveys with quasar or radio galaxy catalogues, or from infrared surveys with colour selection biased towards AGN, while a further 14 have been found through submillimetre observations of known high-redshift AGN. In this paper I argue, on the basis of detailed modelling of the spectral energy distributions of hyperluminous galaxies with accurate radiative transfer models, and from evidence of high gas mass in several cases, that the bulk of the emission from these galaxies at rest frame wavelengths ≥50 μm is caused by star formation. Even after correction for the effects of lensing, hyperluminous galaxies with emission peaking at rest frame wavelengths ≥50 μm are therefore undergoing star formation at rates >10³ M_⊙ yr⁻¹ and are strong candidates for being primeval galaxies, in the process of a major episode of star formation.     

IRAS 22036+5306: an Al2O3 oxide-dominated post-AGB star

Radiative transfer modelling of the Infrared Space Observatory ( ISO ) spectrum of IRAS 22036+5306 has shown that its unusual 11-μm band can be suitably modelled with an alumina-olivine mixture substantially dominated by the former. The results of this work add further credence to recent findings that significant amounts of Al  ₂O₃  dust grains are present in the dust shells of stars near or beyond the tip of the asymptotic giant branch. Indeed, in the case of IRAS 22036+5306, Al  ₂O₃  dominates the dust composition to the extent that it shifts the 9.8-μm band due to amorphous silicates to 11 μm. IRAS 22036+5306 may be an unusual case in that the inner dust torus is maintained at a sufficiently high temperature for Al  ₂O₃  condensation, but not silicate.     

Testing the circumstellar disc hypothesis: a search for H2 outflow signatures from massive young stellar objects with linearly distributed methanol masers

The results of a survey searching for outflows using near-infrared imaging are presented. Targets were chosen from a compiled list of massive young stellar objects associated with methanol masers in linear distributions. Presently, it is a widely held belief that these methanol masers are found in (and delineate) circumstellar accretion discs around massive stars. If this scenario is correct, one way to test the disc hypothesis is to search for outflows perpendicular to the methanol maser distributions. The main objective of the survey was to obtain wide-field near-infrared images of the sites of linearly distributed methanol masers using a narrow-band 2.12-μm filter. This filter is centred on the  H₂ v = 1–0 S(1)  line; a shock diagnostic that has been shown to successfully trace CO outflows from young stellar objects. 28 sources in total were imaged of which 18 sources display H₂ emission. Of these, only two sources showed emission found to be dominantly perpendicular to the methanol maser distribution. Surprisingly, the H₂ emission in these fields is not distributed randomly, but instead the majority of sources are found to have H₂ emission dominantly parallel to their distribution of methanol masers. These results seriously question the hypothesis that methanol masers exist in circumstellar discs. The possibility that linearly distributed methanol masers are instead directly associated with outflows is discussed.     

2–16 μm spectroscopy of micron-sized enstatite (Mg,Fe)2Si2O6 silicates from primitive chondritic meteorites

We present mid-infrared spectra from individual enstatite silicate grains separated from primitive type 3 chondritic meteorites. The 2–16 μm transmission spectra were taken with microspectroscopic Fourier-transform infrared (FT-IR) techniques as part of a project to produce a data base of infrared spectra from minerals of primitive meteorites for comparison with astronomical spectra. In general, the wavelength of enstatite bands increases with the proportion of Fe. However, the wavelengths of the strong En₁₀₀ bands at 10.67 and 11.67 decrease with increasing Fe content. The 11.67-μm band exhibits the largest compositional wavelength shift (twice as large as any other). Our fits of the linear dependence of the pyroxene peaks indicate that crystalline silicate peaks in the 10-μm spectra of Herbig AeBe stars, HD 179218 and 104237, are matched by pyroxenes of En₉₀₋₉₂ and En₇₈₋₈₀, respectively. If these simplistic comparisons with the astronomical grains are correct, then the enstatite pyroxenes seen in these environments are more Fe-rich than are the forsterite (Fo₁₀₀) grains identified in the far-infrared which are found to be Mg end-member grains. This differs from the general composition of type 3 chondritic meteoritic grains in which the pyroxenes are more Mg-rich than are the olivines from the same meteorite.     

Sakurai's Object: characterizing the near-infrared CO ejecta between 2003 and 2007

We present observations of Sakurai's Object obtained at 1–5 μm between 2003 and 2007. By fitting a radiative transfer model to an echelle spectrum of CO fundamental absorption features around  4.7 μm  , we determine the excitation conditions in the line-forming region. We find  ¹²C/¹³C = 3.5^+2.0_−1.5  , consistent with CO originating in ejecta processed by the very late thermal pulse, rather than in the pre-existing planetary nebula. We demonstrate the existence of  2.2 × 10⁻⁶≤ M _CO≤ 2.7 × 10⁻⁶ M_⊙  of CO ejecta outside the dust, forming a high-velocity wind of  500 ± 80 km s⁻¹  . We find evidence for significant weakening of the CO band and cooling of the dust around the central star between 2003 and 2005. The gas and dust temperatures are implausibly high for stellar radiation to be the sole contributor.     

Ice in the Southern Coalsack

To better understand the conditions under which ice mantles form on grains in molecular clouds, three globules in the Southern Coalsack have been searched for the presence of H₂O ice. Given the total lack of star formation in the Coalsack, it is an ideal site for studying unprocessed icy molecular mantles. In our sample of eight field stars lying behind the Coalsack we detect strong H₂O ice absorption in the lines of sight to two stars and possible weak absorption in four others. We estimate H₂O ice column densities or upper limits for these lines of sight. Compared to dark clouds such as Taurus, the Coalsack H₂O ice column densities are lower than expected given the quiescent nature of the Coalsack region. It is possible that the chemical evolution of the Coalsack may simply be at too early a stage for significant ice mantles to appear on the grains, except perhaps in the densest parts of some of the globules. Alternatively, the presence or absence of ice absorption may be related to the distribution of dust along each line of sight, specifically, the relative contributions of dense globules and a more extended diffuse component. For example, our observations are consistent with an ice threshold extinction similar to that observed in the Taurus dark cloud if extinction amounting to   A _V∼5  towards Globules 2 and 3 arises in the extended component. Globule 1 appears to have no extended component.     

Limits on the 2.2-μm contrast ratio of the close-orbiting planet HD 189733b

We obtained 238 spectra of the close-orbiting extrasolar giant planet HD 189733b with resolution   R ∼ 15 000  during one night of observations with the Near-Infrared High-Resolution Spectrograph (NIRSPEC), at the Keck II Telescope. We have searched for planetary absorption signatures in the  2.0–2.4 μm  region where H₂O and CO are expected to be the dominant atmospheric opacities. We employ a phase-dependent orbital model and tomographic techniques to search for the planetary absorption signatures in the combined stellar and planetary spectra. Because potential absorption signatures are hidden in the noise of each single exposure, we use a model list of lines to apply a spectral deconvolution. The resulting mean profile possesses a signal-to-noise ratio (S/N) that is 20 times greater than that found in individual lines. Our spectral time series thus yields spectral signatures with a mean S/N = 2720. We are unable to detect a planetary signature at a contrast ratio of  log₁₀( F _p/ F _*) =−3.40  , with 63.8 per cent confidence. Our findings are not consistent with model predictions which nevertheless give a good fit to mid-infrared observations of HD 189733b. The 1σ result is a factor of 1.7 times less than the predicted 2.185-μm planet/star flux ratio of  log₁₀( F _p/ F _*) ∼−3.16  .     

Submillimetre continuum images of the NGC 2024 star-forming ridge

We present 450- and 800-μm images, made with the James Clerk Maxwell Telescope, of the NGC 2024 molecular ridge. The seven previously known compact cores, FIR1–7, have been detected, and FIR5 has been resolved into a compact object and an associated extended source to the east. The estimated masses of the dense cores vary between 1.6 and 5.1 M⊙ per 14-arcsec beam, assuming a dust temperature of 30 K and a dust opacity of κ_800 μm = 0.002 m² kg⁻¹. A spectral index map made from the 450- and 800-μm images shows spatial variations, with the spectral index, α ( F _ν ∝ ν^α), being systematically lower towards the dense cores. We interpret this as evidence for a lower value of the frequency dependence of the dust opacity, β, towards the denser cores relative to the surrounding molecular material. This may indicate that grain growth is occurring in the cores, prior to planetesimal formation. By comparing the high-resolution 450-μm image with interferometer maps of the integrated CS(2–1) emission, the previously reported discrepancy between dust continuum emission and molecular line emission is found to be very localized. Depletion and temperature variations are discussed as possible explanations.     

The initial conditions of isolated star formation – VII. Spitzer mapping of pre-stellar cores

We have retrieved Spitzer archive data of pre-stellar cores taken with the Multiband Imaging Photometer for Spitzer (MIPS) at a wavelength of 160 μm. Seventeen images, containing 18 cores, were constructed. Flux densities were measured for each core, and background estimates were made. Mean off-source backgrounds were found to be 48 ± 10 MJy sr⁻¹ in Taurus and 140 ± 55 MJy sr⁻¹ in Ophiuchus. Consistency was found between the MIPS 170-μm and ISOPHOT 160-μm calibrations. Fourteen cores were detected both by MIPS and by our previous submillimetre surveys. Spectral energy distributions were made for each core, using additional 24- and 70-μm data from the Spitzer data archive, as well as previous infrared and submillimetre data. Previous temperature estimates were refined, and new temperature estimates were made where no Infrared Space Observatory ( ISO ) data exist. A temperature range of 8–18 K was found for the cores, with most lying in the range 10–13 K. We discount recent claims that a large number of pre-stellar cores may have been misclassified and in fact contain low-luminosity protostars detectable only by Spitzer . We find no new protostars in our sample other than that previously reported in L1521F. It is shown that this has a negligible effect on pre-stellar lifetime estimates.     

Infrared mapping of the dust disc around Vega

ISOPHOT has been used to perform high-resolution 60-μm scans of Vega, and these have been compared with those from γDra, to obtain a Gaussian width of 22±2 arcsec. The dust disc around Vega has been mapped, resolving it at 60 and 90 μm with ISOPHOT. At 90 μm a Gaussian width of 36±3 arsec has been derived. In addition, multi-filter photometry is presented, at 25, 60, 80, 100, 120, 150, 170 and 200 μm. The data are fitted by a modified blackbody with a temperature of 73 K [ Q (λ)∝1/λ^1.1]. The dust disc has a luminosity L _IR/ L _*∼3×10⁻⁵. Using a distance of 7.8 pc, 22 arcsec corresponds to a distance of 86 au, and 36 arcsec to a distance of 140 au.     

The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

Thermal desorption of water ice in the interstellar medium

Water (H₂O) ice is an important solid constituent of many astrophysical environments. To comprehend the role of such ices in the chemistry and evolution of dense molecular clouds and comets, it is necessary to understand the freeze-out, potential surface reactivity and desorption mechanisms of such molecular systems. Consequently, there is a real need from within the astronomical modelling community for accurate empirical molecular data pertaining to these processes. Here we give the first results of a laboratory programme to provide such data. Measurements of the thermal desorption of H₂O ice, under interstellar conditions, are presented. For ice deposited under conditions that realistically mimic those in a dense molecular cloud, the thermal desorption of thin films (≪50 molecular layers) is found to occur with zeroth-order kinetics characterized by a surface binding energy, E _des, of 5773 ± 60 K, and a pre-exponential factor, A , of 10^30 ± 2 molecules cm⁻² s⁻¹. These results imply that, in the dense interstellar medium, thermal desorption of H₂O ice will occur at significantly higher temperatures than has previously been assumed.     

Unlocking the Keyhole: H2 and PAH emission from molecular clumps in the Keyhole Nebula

To better understand the environment surrounding CO emission clumps in the Keyhole Nebula, we have made images of the region in H₂ 1–0 S(1) (2.122-μm) emission and polycyclic aromatic hydrocarbon (PAH) emission at 3.29 μm. Our results show that the H₂ and PAH emission regions are morphologically similar, existing as several clumps, all of which correspond to CO emission clumps and dark optical features. The emission confirms the existence of photodissociation regions (PDRs) on the surface of the clumps. By comparing the velocity range of the CO emission with the optical appearance of the H₂ and PAH emission, we present a model of the Keyhole Nebula whereby the most negative velocity clumps are in front of the ionization region, the clumps at intermediate velocities are in it and those which have the least negative velocities are at the far side. It may be that these clumps, which appear to have been swept up from molecular gas by the stellar winds from η  Car, are now being overrun by the ionization region and forming PDRs on their surfaces. These clumps comprise the last remnants of the ambient molecular cloud around η Car.     

Introducing a hybrid radiative transfer method for smoothed particle hydrodynamics

A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods – the polytropic cooling approximation as devised by Stamatellos et al. and flux-limited diffusion. This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of H₂ dissociation, H⁰ ionization, He⁰ and He⁺ ionization, ice evaporation, dust sublimation, molecular absorption, bound-free and free–free transitions and electron scattering. The method is tested in several scenarios, including (i) the evolution of a  0.07 M_⊙  protoplanetary disc surrounding a  0.5 M_⊙  star; (ii) the collapse of a  1 M_⊙  protostellar cloud and (iii) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.