Probing the surfaces of interstellar dust grains: the adsorption of CO at bare grain surfaces

A solid-state feature was detected at around 2175 cm⁻¹ towards 30 embedded young stellar objects in spectra obtained using the Infrared Spectrometer and Array Camera at the European Southern Observatory Very Large Telescope. We present results from laboratory studies of CO adsorbed at the surface of zeolite wafers, where absorption bands were detected at 2177 and 2168 cm⁻¹ (corresponding to CO chemisorbed at the zeolite surface) and 2130 cm⁻¹ (corresponding to CO physisorbed at the zeolite surface), providing an excellent match to the observational data. We propose that the main carrier of the 2175-band is CO chemisorbed at bare surfaces of dust grains in the interstellar medium. This result provides the first direct evidence that gas–surface interactions do not have to result in the formation of ice mantles on interstellar dust. The strength of the 2175-band is estimated to be  ∼4 × 10⁻¹⁹ cm  molecule⁻¹. The abundance of CO adsorbed at bare grain surfaces ranges from 0.06 to 0.16 relative to H₂O ice, which is, at most, half of the abundance (relative to H₂O ice) of CO residing in H₂O-dominated ice environments. These findings imply that interstellar grains have a large (catalytically active) surface area, providing a refuge for interstellar species. Consequently, the potential exists for heterogeneous chemistry to occur involving CO molecules in unique surface chemistry pathways not currently considered in gas grain models of the interstellar medium.     

Accurate positions of H2O masers in the Large Magellanic Cloud

Positions with subarcsecond accuracy have been measured for seven 22-GHz H₂O masers associated with H  ii regions in the Large Magellanic Cloud (LMC); two of the masers are new detections. Initial position measurements were obtained with the 70-m antenna of the Canberra NASA Deep Space Network during a period of more than two years in which the antenna was used to monitor the maser emission. The positions were further improved using 22-GHz observations involving three antennas of the Australia Telescope Compact Array.
The positions have been compared with those of 1.6-GHz continuum emission and other LMC masers (of OH and CH₃OH). The H₂O maser positions range from within 1 arcsec (270 mpc) of the centre of a compact H  ii component to beyond the boundary of significant continuum emission. Three of the four masers located near continuum peaks are close to OH masers. In two cases the positional agreement is better than 0.2 arcsec (53 mpc); in the third case the agreement is worse (0.9 arcsec) but the positions of the individual H₂O features appear to be spread over more than 1 arcsec. The velocities of the OH masers are within the spread of the H₂O velocities. The three H₂O masers offset from continuum centres are located  3–7 arcsec  from optical or infrared phenomena probably associated with very early stages of star formation; no other molecular masers are known in these directions.     

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.     

Interaction between the north-eastern boundary of Sgr A East and giant molecular clouds

We have detected the   v = 1 → 0 S(1) (λ= 2.1218 μm)  and   v = 2 → 1 S(1) (λ= 2.2477 μm)  lines of H₂ in the Galactic Centre, in a  90 × 27 arcsec²  region between the north-eastern boundary of the non-thermal source Sgr A East, and the giant molecular cloud (GMC)  M−0.02 − 0.07  . The detected  H₂ v = 1 → 0  S(1) emission has an intensity of  1.6–21 × 10⁻¹⁸ W m⁻² arcsec⁻²  and is present over most of the region. Along with the high intensity, the large linewidths  (FWHM = 40–70 km s⁻¹)  and the  H₂ v = 2 → 1 S(1)  to   v = 1 → 0 S(1)  line ratios (0.3–0.5) can be best explained by a combination of C-type shocks and fluorescence. The detection of shocked H₂ is clear evidence that Sgr A East is driving material into the surrounding adjacent cool molecular gas. The H₂ emission lines have two velocity components at ∼+50 and  ∼0 km s⁻¹  , which are also present in the NH₃(3, 3) emission mapped by McGary, Coil & Ho. This two-velocity structure can be explained if Sgr A East is driving C-type shocks into both the  GMC M−0.02 − 0.07  and the northern ridge of McGary et al.     

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.     

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.     

Atmospheric transmission at submillimetre wavelengths from Mauna Kea

The submillimetre atmospheric transmission spectrum above Mauna Kea has been measured at a resolution of 0.005 cm⁻¹ (150 MHz) with a Fourier transform spectrometer at the James Clerk Maxwell Telescope, using the Sun as a source. Column abundances of O₂, H₂O and O₃ determined from these spectra are found to be in excellent agreement with independent measurements. The derived column abundances have been used as inputs to the atmospheric spectral modelling program fascod . The synthetic transmission spectrum is found to be in excellent agreement with the measured spectrum, and provides a template for submillimetre observations from the JCMT.     

The Gunn–Peterson effect and the Lyman alpha forest

We show that spatial correlations in a stochastic large-scale velocity field in an otherwise smooth intergalactic medium (homogeneous comoving density) superposed on the general Hubble flow may cause a 'line-like' structure in QSO spectra similar to the population of unsaturated Lyα forest lines which usually are attributed to individual clouds with 10¹¹ ≲ N _{H i}  5 × 10¹³ cm⁻². Therefore there is no clear observational distinction between a diffuse intergalactic medium and discrete intergalactic clouds. It follows that the H  i density in the diffuse intergalactic medium might be substantially underestimated if it is determined from the observed intensity distribution near the apparent continuum in high-resolution spectra of QSOs. Our tentative estimate implies a diffuse neutral hydrogen opacity τ_GP ∼ 0.3 at z  ∼ 3 and a current baryon density Ω_IGM ≃ 0.08, assuming a Hubble constant H ₀ = 70 km s⁻¹ Mpc⁻¹.     

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  .     

Molecular hydrogen emission in Cygnus A

We present J , H and K -band spectroscopy of Cygnus A, spanning 1.0–2.4 μm in the rest-frame and hence several rovibrational H₂, H recombination and [Fe  ii ] emission lines. The lines are spatially extended by up to 6 kpc from the nucleus, but their distinct kinematics indicate that the three groups (H, H₂ and [Fe  ii ]) are not wholly produced in the same gas. The broadest line, [Fe  ii ] λ 1.644, exhibits a non-Gaussian profile with a broad base (FWHM≃1040 km s⁻¹), perhaps because of the interaction with the radio source. Extinctions to the line-emitting regions substantially exceed earlier measurements based on optical H recombination lines.
Hard X-rays from the quasar nucleus are likely to dominate the excitation of the H₂ emission. The results of Maloney, Hollenbach & Tielens are thus used to infer the total mass of gas in H₂ v=1–0 S(1)-emitting clouds as a function of radius, for gas densities of 10³ and 10⁵ cm⁻³, and stopping column densities N _H=10²²–10²⁴ cm⁻². Assuming azimuthal symmetry, at least 2.3×10⁸ M_⊙ of such material is present within 5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno et al. actually arises in a circumnuclear torus, the implied gas mass rises to ∼10¹⁰ M_⊙. The latter plausibly accounts for 10⁹ yr of mass deposition from the cluster cooling flow, for which within this radius.     

On the possibility of observing H2 emission from primordial molecular cloud kernels

We study the prospects for observing H₂ emission during the assembly of primordial molecular cloud kernels. The primordial molecular cloud cores, which resemble those at the present epoch, can emerge around  1+ z ∼20  according to recent numerical simulations. The kernels form inside the cores, and the first stars will appear inside the kernels. A kernel typically contracts to form one of the first generation stars with an accretion rate that is as large as ∼0.01 M_⊙ yr⁻¹. This occurs owing to the primordial abundances, which result in a kernel temperature of order 1000 K, and the collapsing kernel emits H₂ line radiation at a rate ∼10³⁵ erg s⁻¹. Predominantly   J =5-3   ( v =0)  rotational emission of H₂ is expected. At redshift  1+ z ∼20  , the expected flux is ∼0.01 μJy for a single kernel. While an individual object is not observable by any facilities available in the near future, the expected assembly of primordial star clusters on subgalactic scales can result in fluxes at the sub-mJy level. This is marginally observable with ASTRO-F and ALMA. We also examine the rotational   J =2-0   ( v =0)  and vibrational   δv =1  emission lines. The former may possibly be detectable with ALMA.     

Application of a clumpy core model to NGC 2071

The multi-transitional observations of CS molecules towards the NGC 2071 core have been re-analysed by using a tri-dimensional Monte Carlo radiative transfer code. Better agreement with the observations is made by an introduction of clumpiness to this model than by smoothly varying density to the 1D microturbulent one. The best-fitting model shows that, when a unique density is assumed for clumps, the volume filling factor of the clumps varies as r ⁻² with an average of ∼5 per cent over the entire core, and that the H₂ number density and the CS abundance of the clump relative to H₂ are ∼ 2 × 10⁶ cm⁻³ and ∼ 6 × 10⁻¹⁰, respectively. The radial density gradient ∝ r ⁻² obtained from our clumpy core model is steeper than that (∝ r ^−1.3) obtained from the microturbulent model. Since all clumps are subject to random bulk motions in this 3D clumpy macroturbulent model, synthesized line profiles do not show self-absorption dips even for opaque transitions and the resulting linewidth is in good accordance with the observations.     

Hydrated sulphuric acid in dense molecular clouds

We consider sulphur depletion in dense molecular clouds, and suggest hydrated sulphuric acid, H₂SO₄ ·  n H₂O, as a component of interstellar dust in icy mantles. We discuss the formation of hydrated sulphuric acid in collapsing clouds and its instability in heated regions in terms of the existing hot core models and observations. We also show that some features of the infrared spectrum of hydrated sulphuric acid have correspondence in the observed spectra of young stellar objects.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

Local subgiants and time-scales of disc formation

This paper concentrates on the relationship between the rate of gas emission from the nucleus of Comet 9P/Tempel 1, the fraction f of the nucleus that is active, and the crater damage inflicted by the recent 2005 July 4 Deep Impact space mission. The cometary nucleus has a surface area of about  1.7 × 10⁸ m²  and a mean radius of about 3700 m. Before the impact it is estimated that only a fraction f = 0.0056 of the nucleus surface was actively producing gas and dust. The active area was about  9.4 × 10⁵ m²  . Absolute magnitudes obtained at recent perihelion passages of this comet indicate that variations in the 0.0074 > f > 0.0039 range can occur from apparition to apparition. Because of the low size of the original active area, the production of a new impact crater in the diameter range 40 to 300 m would lead to a long-term change in the cometary visual magnitude in the range 0.0018 to 0.098 respectively. This is below the limit of detectability. It has been suggested that the cometary dust is in the form of 'talcum powder' not 'beach sand'. We suggest that the dust ejected from the impact site has been broken up by the energetic impact process and thus has a different size distribution from dust locked in the snowy matrix of the nucleus and normally lifted off the nucleus by gentle sublimation processes.     

A high-accuracy computed water line list

A computed list of  H₂¹⁶O  infrared transition frequencies and intensities is presented. The list, BT2, was produced using a discrete variable representation two-step approach for solving the rotation–vibration nuclear motions. It is the most complete water line list in existence, comprising over 500 million transitions (65 per cent more than any other list) and it is also the most accurate (over 90 per cent of all known experimental energy levels are within 0.3 cm⁻¹ of the BT2 values). Its accuracy has been confirmed by extensive testing against astronomical and laboratory data.
The line list has been used to identify individual water lines in a variety of objects including comets, sunspots, a brown dwarf and the nova-like object V838 Mon. Comparison of the observed intensities with those generated by BT2 enables water abundances and temperatures to be derived for these objects. The line list can also be used to provide an opacity for models of the atmospheres of M dwarf stars and assign previously unknown water lines in laboratory spectra.     

Shocked molecular gas towards the supernova remnant G359.1–0.5 and the Snake

We have found a bar of shocked molecular hydrogen (H₂) towards the OH(1720 MHz) maser located at the projected intersection of supernova remnant (SNR)  G359.1–0.5  and the non-thermal radio filament known as the Snake. The H₂ bar is well aligned with the SNR shell and almost perpendicular to the Snake. The OH(1720 MHz) maser is located inside the sharp western edge of the H₂ emission, which is consistent with the scenario in which the SNR drives a shock into a molecular cloud at that location. The spectral line profiles of ¹²CO, HCO⁺ and CS towards the maser show broad-line absorption, which is absent in the ¹³CO spectra and most probably originates from the pre-shock gas. A density gradient is present across the region and is consistent with the passage of the SNR shock, while the H₂ filament is located at the boundary between the pre-shock and post-shock regions.     

The excitation and kinematical properties of H2 and [Fe ii] in the HH 46/47 bipolar outflow

Long-slit spectra of the molecular outflow Herbig–Haro (HH) 46/47 have been taken in the J and K near-infrared bands. The observed H₂ line emission confirms the existence of a bright and extended redshifted counter-jet outflow south-west of HH 46. In contrast with the optical appearance of this object, we show that this outflow seems to be composed of two different emission regions characterized by distinct heliocentric velocities. This implies an acceleration of the counter-jet.
The observed [Fe  ii ] emission suggests an average extinction of 7–9 visual magnitudes for the region associated with the counter-jet.
Through position–velocity diagrams, we show the existence of different morphologies for the H₂ and [Fe  ii ] emission regions in the northern part of the HH 46/47 outflow. We have detected for the first time high-velocity (−250 km s⁻¹) [Fe  ii ] emission in the region bridging HH 46 to HH 47A. The two strong peaks detected can be identified with the optical positions B8 and HH 47B.
The H₂ excitation diagrams for the counter-jet shock suggest an excitation temperature for the gas of T _ex≈2600 K . The lack of emission from the higher energy H₂ lines, such as the 4–3 S(3) transition, suggests a thermal excitation scenario for the origin of the observed emission. Comparison of the H₂ line ratios with various shock models yielded useful constraints about the geometry and type of these shocks. Planar shocks can be ruled out whereas curved or bow shocks (both J- and C-type) can be parametrized to fit our data.     

A spectroscopic study of IRAS F10214 + 4724

The z  = 2.286  IRAS galaxy F10214 + 4724 remains one of the most luminous galaxies in the Universe, despite its gravitational lens magnification. We present optical and near-infrared spectra of F10214 + 4724, with clear evidence for three distinct components: lines of width ∼ 1000 km s⁻¹ from a Seyfert 2 nucleus; ≲ 200 km s⁻¹ lines which are likely to be associated with star formation; and a broad (∼ 4000 km s⁻¹) C  III ] 1909-Å emission line which is blueshifted by ∼ 1000 km s⁻¹ with respect to the Seyfert 2 lines. Our study of the Seyfert 2 component leads to several new results. (i) From the double-peaked structure in the Lyα line, and the lack of Lyβ, we argue that the Lyα photons have emerged through a neutral column of N _H ∼ 2.5 × 10²⁵ m⁻², possibly located within the AGN narrow-line region, as proposed for several high-redshift radio galaxies. (ii) The resonant O  VI 1032, 1036-Å doublet (previously identified as Lyβ) is in an optically thick (1:1) ratio. At face value this implies an extreme density ( n _e ∼ 10¹⁷ m⁻³) more typical of broad-line region clouds. However, we attribute this instead to the damping wings of Lyβ from the resonant absorption. (iii) A tentative detection of He  II 1086 suggests little extinction in the rest frame ultraviolet.     

Exciting molecular hydrogen in the central galaxies of cooling flows

The origin of rovibrational H₂ emission in the central galaxies of cooling flow clusters is poorly understood. Here we address this issue using data from our near-infrared spectroscopic survey of 32 of the most line-luminous such systems, presented in the companion paper by Edge et al.
We consider excitation by X-rays from the surrounding intracluster medium (ICM), ultra-violet (UV) radiation from young stars, and shocks. The   v = 1–0  K -band lines with upper levels within  10⁴ K  of the ground state appear to be mostly thermalized (implying gas densities  ≳10⁵ cm⁻³  ), with the excitation temperature typically exceeding 2000 K, as found earlier by Jaffe, Bremer & van der Werf. Together with the lack of strong   v = 2–0  lines in the H -band, this rules out UV radiative fluorescence.
Using the cloudy photoionization code, we deduce that the H₂ lines can originate in a population of dense clouds, exposed to the same hot  ( T ∼ 50 000 K)  stellar continuum as the lower density gas which produces the bulk of the forbidden optical line emission in the Hα-luminous systems. This dense gas may be in the form of self-gravitating clouds deposited directly by the cooling flow, or may instead be produced in the high-pressure zones behind strong shocks. Furthermore, the shocked gas is likely to be gravitationally unstable, so collisions between the larger clouds may lead to the formation of globular clusters.