The cooling function of HD molecule revisited

We report new calculations of the cooling rate of primordial gas by the HD molecule, taking into account its ro-vibrational structure. The HD cooling function is calculated including radiative and collisional transitions for   J ≤ 8  rotational levels, and for the vibrational levels v = 0, 1, 2 and 3. The ro-vibrational level population is calculated from the balance equation assuming steady state. The cooling function is evaluated in the ranges of the kinetic temperatures, T _k, from 10² to  2 × 10⁴ K  and the number densities, n _H, from 1 to  10⁸ cm⁻³  . We find that the inclusion of collisional ro-vibrational transitions increases significantly the HD cooling efficiency, in particular for high densities and temperatures. For   n _H≳ 10⁵  and   T _k∼ 10⁴ K  the cooling function becomes more than an order of magnitude higher than previously reported. We give also the HD cooling rate in the presence of the cosmic microwave radiation field for radiation temperatures of 30, 85 and 276 K (redshifts of 10, 30 and 100). The tabulated cooling functions are available at http://www.cifus.uson.mx/Personal_Pages/anton/DATA/HD_cooling/HD_cool.html . We discuss the relevance to explore the effects of including our results into models and simulations of galaxy formation, especially in the regime when gas cools down from temperatures above ∼3000 K.     

ISO observations of molecular hydrogen in the DR 21 bipolar outflow

We demonstrate that a wide range of molecular hydrogen excitation can be observed in protostellar outflows at wavelengths in excess of 5 μm. Cold H₂ in DR 21 is detected through the pure rotational transitions in the ground vibrational level (0–0). Hot H₂ is detected in pure rotational transitions within higher vibrational levels (1–1, 1–2, etc.). Although this emission is relatively weak, we have detected two 1–1 lines in the DR 21 outflow with the ISO SWS instrument. We thus investigate molecular excitation over energy levels corresponding to the temperature range 1015–15 722 K, without the uncertainty introduced by differential extinction when employing near-infrared data.
This gas is thermally excited. We uncover a rather low H₂ excitation in the DR 21 West Peak. The line emission cannot be produced from single C-shocks or J-shocks; a range of shock strengths is required. This suggests that bow shocks and/or bow-generated supersonic turbulence is responsible. We are able to distinguish this shock-excited gas from the fluoresced gas detected in the K band, providing support for the dual-excitation model of Fernandes, Brand & Burton.     

The thermal balance of the first structures in the primordial gas

The contraction of matter in the primordial medium, to form the first gravitationally bound structures, was mediated by radiative cooling of the gas by H₂ and HD. We have computed the initial phases of free-fall collapse, incorporating the results of quantum mechanical calculations of rate coefficients for collisional excitation of H₂ and HD by the principal perturbers, H, He, H₂ and H⁺. The structure of shock waves produced when the collapse speed exceeds the local sound speed is determined. In the post-shock gas, radiative cooling by H₂ exceeds that by HD, but by a factor of only 4. The intensities of the strongest emission lines of H₂– rotational transitions within the vibrational ground state – are calculated. Even with coarse spectral and angular resolution, these transitions might be observable as inhomogeneities in the cosmic background radiation.     

The cooling of astrophysical media by HD

The results of recent quantum mechanical calculations of cross-sections for rotational transitions within the vibrational ground state of HD are used to evaluate the rate of radiative energy loss from gas containing HD, in addition to H, He and H₂. The cooling function for HD (i.e. the rate of cooling per HD molecule) is evaluated in steady state on a grid of values of the relevant parameters of the gas, namely the gas density and temperature, the atomic to molecular hydrogen abundance ratio and the ortho:para-H₂ density ratio. The corresponding cooling function for H₂, previously computed by Le Bourlot et al., is slightly revised to take account of transitions induced by collisions with ground-state ortho-H₂ ( J =1). The cooling functions and the data required for their calculation are available from http://ccp7.dur.ac.uk/. We then make a study of the rate of cooling of the primordial gas through collisions with H₂ and HD molecules. In this case, radiative transitions induced by the cosmic background radiation field and, in the case of H₂, collisional transitions induced by H⁺ ions should additionally be included.     

The rotational excitation of H2 by H2

We present rate coefficients for rotational transitions induced in collisions between H₂ molecules. Rotational levels J  ≤ 8 and kinetic temperatures T  ≤ 1000 K are considered. The interaction potential computed by Schwenke has been used, together with the quantal coupled channels method of calculating the cross-sections. Comparison is made with the more recent of previous results.     

Ultra-high-resolution observations of circumstellar K i and C2 around the post-AGB star HD 56126

We have used the Ultra-High-Resolution Facility (UHRF) at the AAT, operating at a resolution of 0.35 km s⁻¹ (FWHM), to observe K  i and C₂ absorption lines arising in the circumstellar environment of the post-AGB star HD 56126. We find three narrow circumstellar absorption components in K  i , two of which are also present in C₂. We attribute this velocity structure to discrete shells resulting from multiple mass-loss events from the star. The very high spectral resolution has enabled us to resolve the intrinsic linewidths of these narrow lines for the first time, and we obtain velocity dispersions ( b -values) of 0.2–0.3 km s⁻¹ for the K  i components, and 0.54±0.03 km s⁻¹ for the strongest (and best defined) C₂ component. These correspond to rigorous kinetic temperature upper limits of 211 K for K  i and 420 K for C₂, although the b -value ratio implies that these two species do not co-exist spatially. The observed degree of rotational excitation of C₂ implies low kinetic temperatures ( T _k≈10 K) and high densities ( n ≈10⁶ to 10⁷ cm⁻³) within the shell responsible for the main C₂ component. Given this low temperature, the line profiles then imply either mildly supersonic turbulence or an unresolved velocity gradient through the shell.     

The distribution of the warm and dense molecular gas around Cepheus A HW 2

We present VLA observations of the ( J , K )=(1, 1), (2, 2), (3, 3) and (4, 4) inversion transitions of NH₃ toward the HW 2 object in Cepheus A, with 1-arcsec angular resolution. Emission is detected in the main hyperfine line of the first three transitions. The NH₃(2, 2) emission shows a non-uniform 'ring' structure, which is more extended (3 arcsec) and intense than the emission seen in the (1, 1) and (3, 3) lines. A rotational temperature of ∼ 30–50 K and a lower limit to the mass of ∼ 1 ( X _NH3/10⁻⁸)⁻¹ M_⊙ are derived for the ring structure. The spatio-kinematical distribution of the NH₃ emission does not seem to be consistent with a simple circumstellar disc around the HW 2 thermal biconical radio jet. We suggest that it represents the remnant of the parental core from which both the inner 300-au (0.4 arcsec) disc, traced by the water maser spots previously found in the region, and the central object have formed. The complex velocity field of this core is probably produced from bound motions (similar to those of the inner disc) and from interaction with outflowing material.     

A submillimetre survey of the kinematics of the Perseus molecular cloud – I. Data

We present submillimetre observations of the   J = 3 → 2  rotational transition of ¹²CO, ¹³CO and C¹⁸O across over 600 arcmin² of the Perseus molecular cloud, undertaken with the Heterodyne Array Receiver Programme (HARP), a new array spectrograph on the James Clerk Maxwell Telescope. The data encompass four regions of the cloud, containing the largest clusters of dust continuum condensations: NGC 1333, IC348, L1448 and L1455. A new procedure to remove striping artefacts from the raw HARP data is introduced. We compare the maps to those of the dust continuum emission mapped with the Submillimetre Common-User Bolometer Array (SCUBA; Hatchell et al.) and the positions of starless and protostellar cores (Hatchell et al.). No straightforward correlation is found between the masses of each region derived from the HARP CO and SCUBA data, underlining the care that must be exercised when comparing masses of the same object derived from different tracers. From the ¹³CO/C¹⁸O line ratio the relative abundance of the two species  ([¹³CO]/[C¹⁸O]∼ 7)  and their opacities (typically τ is 0.02–0.22 and 0.15–1.52 for the C¹⁸O and ¹³CO gas, respectively) are calculated. C¹⁸O is optically thin nearly everywhere, increasing in opacity towards star-forming cores but not beyond  τ₁₈∼ 0.9  . Assuming the ¹²CO gas is optically thick, we compute its excitation temperature, T _ex (around 8–30 K), which has little correlation with estimates of the dust temperature.     

Fe xiii emission lines in active region spectra obtained with the Solar Extreme-Ultraviolet Research Telescope and Spectrograph

Recent fully relativistic calculations of radiative rates and electron impact excitation cross-sections for Fe  xiii are used to generate emission-line ratios involving 3s²3p²–3s3p³ and 3s²3p²–3s²3p3d transitions in the 170–225 and 235–450 Å wavelength ranges covered by the Solar Extreme-Ultraviolet Research Telescope and Spectrograph (SERTS). A comparison of these line ratios with SERTS active region observations from rocket flights in 1989 and 1995 reveals generally very good agreement between theory and experiment. Several new Fe  xiii emission features are identified, at wavelengths of 203.79, 259.94, 288.56 and 290.81 Å. However, major discrepancies between theory and observation remain for several Fe  xiii transitions, as previously found by Landi and others, which cannot be explained by blending. Errors in the adopted atomic data appear to be the most likely explanation, in particular for transitions which have 3s²3p3d ¹D₂ as their upper level. The most useful Fe  xiii electron-density diagnostics in the SERTS spectral regions are assessed, in terms of the line pairs involved being (i) apparently free of atomic physics problems and blends, (ii) close in wavelength to reduce the effects of possible errors in the instrumental intensity calibration, and (iii) very sensitive to changes in N _e over the range  10⁸–10¹¹ cm⁻³  . It is concluded that the ratios which best satisfy these conditions are 200.03/202.04 and 203.17/202.04 for the 170–225 Å wavelength region, and 348.18/320.80, 348.18/368.16, 359.64/348.18 and 359.83/368.16 for 235–450 Å.     

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.     

The continuing saga of Sakurai's object (V4334 Sgr): dust production and helium line emission

We report further UKIRT spectroscopic observations of Sakurai's object (V4334 Sgr) made in 1999 April/May in the 1–4.75 μm range, and find that the emission is dominated by amorphous carbon at T _d~600 K. The estimated maximum grain size is 0.6 μm, and the mass lower limit is 1.7±0.2×10⁻⁸ M_⊙ to 8.9±0.6×10⁻⁷ M_⊙ for distances of 1.1–8 kpc. For 3.8 kpc the mass is 2.0±0.1×10⁻⁷ M_⊙.
We also report strong He  i emission at 1.083 μm, in contrast to the strong absorption in this line in 1998. We conclude that the excitation is collisional, and is probably caused by a wind, consistent with the P Cygni profile observed by Eyres et al. in 1998.     

The minimum mass ratio of W UMa-type binary systems

When the total angular momentum of a binary system   J _tot= J _orb+ J _spin  is at a certain critical (minimum) value, a tidal instability occurs which eventually forces the stars to merge into a single, rapidly rotating object. The instability occurs when   J _orb= 3 J _spin  , which in the case of contact binaries corresponds to a minimum mass ratio   q _min≈  0.071–0.078. The minimum mass ratio is obtained under the assumption that stellar radii are fixed and independent. This is not the case with contact binaries where, according to the Roche model, we have   R ₂= R ₂( R ₁, a , q )  . By finding a new criterion for contact binaries, which arises from  d J _tot= 0  , and assuming   k ²₁≠ k ²₂  for the component's dimensionless gyration radii, a theoretical lower limit   q _min= 0.094–0.109  for overcontact degree   f = 0–1  is obtained.     

Cosmic ray generation by quasar remnants: constraints and implications

We report the first detection of CO in the bulge of M31. The ¹²CO (1–0) and (2–1) lines are both detected in the dust complex D395A/393/384, at 1.3 arcmin (∼0.35 kpc) from the centre. From these data and from visual extinction data, we derive a CO luminosity to reddening ratio (and a CO luminosity to H₂ column density ratio) quite similar to that observed in the local Galactic clouds. The (2–1) to (1–0) line intensity ratio points to a CO rotational temperature and a gas kinetic temperature of >10 K. The molecular mass of the complex, inside a 25-arcsec (100 pc) region, is 1.5×10⁴ M_⊙.     

Electron-impact rotational excitation of linear molecular ions

Molecular R -matrix calculations are performed to give rotational excitation rates for electron collisions with linear molecular ions. Results are presented for CO⁺, HCO⁺, NO⁺ and H₂⁺ up to electron temperatures of 10 000 K. De-excitation rates and critical electron densities are also given. It is shown that the widely used Coulomb–Born approximation is valid for Δ j =1 transitions when the molecular ion has a dipole greater than about 2D, but otherwise is not reliable for studying electron-impact rotational excitation. In particular, transitions with Δ j >1 are found to have appreciable rates and are found to be entirely dominated by short-range effects.     

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.     

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.     

Near-infrared properties of i-drop galaxies in the Hubble Ultra Deep Field

We analyse near-infrared Hubble Space Telescope ( HST )/Near-Infrared Camera and Multi-Object Spectrometer F 110 W ( J ) and F 160 W ( H ) band photometry of a sample of 27 i '-drop candidate   z ≃ 6  galaxies in the central region of the HST /Advanced Camera for Surveys Ultra Deep Field . The infrared colours of the 20 objects not affected by near neighbours are consistent with a high-redshift interpretation. This suggests that the low-redshift contamination of this i '-drop sample is smaller than that observed at brighter magnitudes, where values of 10–40 per cent have been reported. The J – H colours are consistent with a slope flat in   f_ν ( f_λ ∝λ⁻²)  , as would be expected for an unreddened starburst. However, there is evidence for a marginally bluer spectral slope  ( f_λ ∝λ^−2.2)  , which is perhaps indicative of an extremely young starburst (∼10 Myr old) or a top heavy initial mass function and little dust. The low levels of contamination, median photometric redshift of   z ∼ 6.0  and blue spectral slope, inferred using the near-infrared data, support the validity of the assumptions in our earlier work in estimating the star formation rates, and that the majority of the i -drop candidates galaxies lie at   z ∼ 6  .     

Extreme ultraviolet transitions of Fe xxi in solar, stellar and laboratory spectra

Recent R -matrix calculations of electron impact excitation rates for transitions among the 2s²2p², 2s2p³ and 2p⁴ levels of Fe  xxi are used to derive theoretical electron density ( N _e) sensitive emission-line ratios involving 2s²2p²–2s2p³ transitions in the ∼98–146Å wavelength range. A comparison of these with observations from the PLT tokamak plasma, for which the electron density has been independently determined, reveals generally very good agreement between theory and experiment, and in some instances removes discrepancies found previously. The observed Fe  xxi ratios for a solar flare, obtained with the OSO–5 satellite, imply electron densities which are consistent, with discrepancies that do not exceed 0.2 dex. In addition, the derived values of N _e are similar to those estimated for the high-temperature regions of other solar flares. The good agreement between theory and observation, in particular for the tokamak spectra, provides experimental support for the accuracy of the present line-ratio calculations, and hence for the atomic data on which they are based.     

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.