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.     

How big were the first cosmological objects?

We calculate the cooling times at constant density for haloes with virial temperatures from 100 K to  1×10⁵ K  that originate from a 3 σ fluctuation of a CDM power spectrum in three different cosmologies. Our intention is to determine the first objects that can cool to low temperatures, but not to follow their dynamical evolution. We identify two generations of haloes: those with low virial temperatures,   T _vir≲9000 K  that remain largely neutral, and those with larger virial temperatures that become ionized. The lower temperature, lower mass haloes are the first to cool to 75 per cent of their virial temperature. The precise temperature and mass of the first objects are dependent upon the molecular hydrogen (H₂) cooling function and the cosmological model. The higher mass haloes collapse later but, in this paradigm, cool much more efficiently once they have done so, first via electronic transitions and then via molecular cooling: in fact, a greater residual ionization once the haloes cool below 9000 K results in an enhanced H₂ production and hence a higher cooling rate at low temperatures than for the lower mass haloes, so that within our constant-density model it is the former that are the first to cool to really low temperatures. We discuss the possible significance of this result in the context of CDM models in which the shallow slope of the initial fluctuation spectrum on small scales leads to a wide range of halo masses (of differing overdensities) collapsing over a small redshift interval. This 'crosstalk' is sufficiently important that both high- and low-mass haloes collapse during the lifetimes of the massive stars which may be formed at these epochs. Further investigation is thus required to determine which generation of haloes plays the dominant role in early structure formation.     

Cosmic rays and the primordial gas

One of the most-outstanding problems in the gravitational collapse scenario of early structure formation is the cooling of primordial gas to allow for small-mass objects to form. As the neutral primordial gas is a poor radiator at temperatures   T ≤ 10⁴ K  , molecular hydrogen is needed for further cooling down to temperatures   T ∼ 100 K  . The formation of molecular hydrogen is catalyzed by the presence of free electrons, which could be provided by the ionization due to an early population of cosmic rays (CRs). In order to investigate this possibility, we developed a code to study the effects of ionizing CRs on the thermal and chemical evolution of primordial gas. We found that CRs can provide enough free electrons needed for the formation of molecular hydrogen, and therefore can increase the cooling ability of such primordial gas under following conditions. A dissociating photon flux with   F < 10⁻¹⁸ erg cm⁻² Hz⁻¹ s⁻¹  , initial temperature of the gas  ∼10³ K  , total gas number densities   n ≥ 1 cm⁻³  , and cosmic-ray sources with     .     

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.     

Density gradients in Galactic ultra-compact H ii regions

We have used recent radiative transfer solutions for cavity-centred shells to investigate the prevalence of density gradients in Galactic ultra-compact H  ii regions. We find that an analysis of 5 and 1.4 GHz data, taken from the recent compilation of Giveon et al., implies that ∼76 per cent of sources may have appreciable density gradients. It would also seem that the central cavities of these sources must be relatively small, with radii no greater than ∼20 per cent of those of the outer shells. The remainder of these sources are presumably homogenous, have much larger cavities, or possess reverse density gradients (densities which increase with increasing radius). A good fraction of the H  ii regions also appear to have high brightness temperatures, implying mean electron temperatures  〈 T _e〉  of the order of  ≈1.3 × 10⁴ K  . This value is higher than has been determined for other such sources.     

Rosseland and Planck mean opacities for primordial matter

We present newly calculated low-temperature opacities for gas with a primordial chemical composition. In contrast to earlier calculations, which took a pure metal-free hydrogen/helium mixture, we take into account the small fractions of deuterium and lithium as resulting from standard big bang nucleosynthesis. Our opacity tables cover the density range  −16 < log ρ[g cm⁻³] < −2  and the temperature range of  1.8 < log  T [K] < 4.6  , while previous tables have usually been restricted to   T > 10³ K  . We find that, while the presence of deuterium does not significantly alter the opacity values, the presence of lithium gives rise to major modifications of the opacities, at some points increasing it by approximately two orders of magnitude relative to pure hydrogen/helium opacities.     

Spin temperatures and covering factors for H i 21-cm absorption in damped Lyman α systems

We investigate the practice of assigning high spin temperatures to damped Lyman α absorption systems (DLAs) not detected in H  i 21-cm absorption. In particular, Kanekar & Chengalur have attributed the mix of 21-cm detections and non-detections in low-redshift  ( z _abs≤ 2.04) DLAs  to a mix of spin temperatures, while the non-detections at high redshift were attributed to high spin temperatures. Below   z _abs= 0.9  , where some of the DLA host galaxy morphologies are known, we find that 21-cm absorption is normally detected towards large radio sources when the absorber is known to be associated with a large intermediate (spiral) galaxy. Furthermore, at these redshifts, only one of the six 21-cm non-detections has an optical identification and these DLAs tend to lie along the sight-lines to the largest background radio continuum sources. For these and many of the high-redshift DLAs occulting large radio continua, we therefore expect covering factors of less than the assumed/estimated value of unity. This would have the effect of introducing a range of spin temperatures considerably narrower than the current range of  Δ T _s≳ 9000 K  , while still supporting the hypothesis that the high-redshift DLA sample comprises a larger proportion of compact galaxies than the low-redshift sample.     

Rovibrational excitation of HD in collisions with atomic and molecular hydrogen

We have computed cross-sections and rate coefficients for rovibrational transitions in HD, induced by collisions with atomic and molecular hydrogen. We employed fully quantum-mechanical methods and the potential of Boothroyd et al. for H–HD, and that of Schwenke for H₂–HD. The rate coefficients for vibrational relaxation v =1→0 of HD are compared with the corresponding values for H₂. The influence of vibrationally excited channels on the rate coefficients for rotational transitions within the v =0 vibrational ground state of HD is shown to be small at T =500 K, where T is the kinetic temperature. The rate coefficients, for 100 T 2000 K, are available from http://ccp7.dur.ac.uk/.     

The temperature of the warm neutral medium in the Milky Way

We report high-spectral-resolution Australia Telescope Compact Array (ATCA) H  i 21-cm observations resulting in the detection of the warm neutral medium (WNM) of the Galaxy in absorption against two extragalactic radio sources, PKS 1814−637 and PKS 0407−658. The two lines of sight were selected on the basis of the simplicity of their absorption profiles and the strength of the background sources; the high velocity resolution of the spectra then enabled us to estimate the kinetic temperatures of the absorbing gas by fitting multiple Gaussians to the absorption profiles. Four separate WNM components were detected towards the two sources, with peak optical depths  τ_max= (1.0 ± 0.08) × 10⁻², (1.4 ± 0.2) × 10⁻³, (2.2 ± 0.5) × 10⁻³  and  (3.4 ± 0.5) × 10⁻³  and kinetic temperatures   T _k= 3127 ± 300, 3694 ± 1595, 3500 ± 1354  and  2165 ± 608 K  , respectively. All four components were thus found to have temperatures in the thermally unstable range  500 < T _k < 5000 K  ; this suggests that thermal equilibrium has not been reached throughout the WNM.     

Magnetic hydrogen atmosphere models and the neutron star RX J1856.5–3754

RX J1856.5−3754 is one of the brightest nearby isolated neutron stars (INSs), and considerable observational resources have been devoted to it. However, current models are unable to satisfactorily explain the data. We show that our latest models of a thin, magnetic, partially ionized hydrogen atmosphere on top of a condensed surface can fit the entire spectrum, from X-rays to optical, of RX J1856.5−3754, within the uncertainties. In our simplest model, the best-fitting parameters are an interstellar column density   N _H≈ 1 × 10²⁰ cm⁻²  and an emitting area with   R ^∞≈ 17 km  (assuming a distance to RX J1856.5−3754 of 140 pc), temperature   T ^∞≈ 4.3 × 10⁵ K  , gravitational redshift   z _g ∼ 0.22  , atmospheric hydrogen column   y _H≈ 1 g cm⁻²  , and magnetic field   B ≈ (3–4) × 10¹² G  ; the values for the temperature and magnetic field indicate an effective average over the surface. We also calculate a more realistic model, which accounts for magnetic field and temperature variations over the NS surface as well as general relativistic effects, to determine pulsations; we find that there exist viewing geometries that produce pulsations near the currently observed limits. The origin of the thin atmospheres required to fit the data is an important question, and we briefly discuss mechanisms for producing these atmospheres. Our model thus represents the most self-consistent picture to date for explaining all the observations of RX J1856.5−3754.     

Relativistic ionized accretion disc models of MCG–6-30-15

We present BeppoSAX observations of Nova Velorum 1999 (V382 Vel), carried out in a broad X-ray band covering 0.1–300 keV only 15 d after the discovery and again after 6 months. The nova was detected at day 15 with the BeppoSAX instruments which measured a flux F _x≃1.8×10⁻¹¹ erg cm⁻² s⁻¹ in the 0.1–10 keV range and a 2 σ upper limit F _x<6.7×10⁻¹² erg cm⁻² s⁻¹ in the 15–60 keV range. We attribute the emission to shocked nebular ejecta at a plasma temperature kT ≃6 keV . At six months no bright component emerged in the 15–60 keV range, but a bright central supersoft X-ray source appeared. The hot nebular component previously detected had cooled to a plasma temperature kT <1 keV . There was strong intrinsic absorption of the ejecta in the first observation and not in the second, because the column density of neutral hydrogen decreased from N (H)≃1.7×10²³ to N (H)≃10²¹ cm⁻² (close to the interstellar value). The unabsorbed X-ray flux also decreased from F _x=4.3×10⁻¹¹ to F _x≃10⁻¹² erg cm⁻² s⁻¹ .     

Theoretical calculations of the H i, He i and He ii free–bound continuous emission spectra

We present coefficients for the calculation of the continuous emission spectra of H  i , He  i and He  ii due to electron–ion recombination. Coefficients are given for photon energies from the first ionization threshold for each ion to the   n = 20  threshold of hydrogen  (36.5 μm)  , and for temperatures  100 ≤  T ≤ 10⁵ K  . The emission coefficients for He  i are derived from accurate ab initio photoionization data. The coefficients are scaled in such a way that they may be interpolated by a simple scheme with uncertainties less than 1 per cent in the whole temperature and wavelength domain. The data are suitable for incorporation into photoionization/plasma codes and should aid with the interpretation of spectra from the very cold ionized gas phase inferred to exist in a number of gaseous clouds.     

New limits on the possible variation of physical constants

Recent detections of high-redshift absorption by both atomic hydrogen and molecular gas in the radio spectra of quasars have provided a powerful tool for measuring possible temporal and spatial variations of physical 'constants' in the universe.
We compare the frequency of high-redshift hydrogen 21-cm absorption with that of associated molecular absorption in two quasars to place new (1σ) upper limits on any variation in y≡g_pα² (where α is the fine-structure constant, and g_p is the proton g -factor of
   
at redshifts z = 0.25 and 0.68. These quasars are separated by a comoving distance of 3000 Mpc ( H ₀= 75 km s⁻¹ Mpc⁻¹ and q ₀). We also derive limits on the time rates of change
   
    between the present epoch and z = 0.68. These limits are more than an order of magnitude smaller than previous results derived from high-redshift measurements.     

Supergiant temperatures and linear radii from near-infrared interferometry

We present angular diameters for 42 Luminosity Class (LC) I stars and 32 LC II stars that have been interferometrically determined with the Palomar Testbed Interferometer. Derived values of radius and effective temperature are established for these objects, and an empirical calibration of these parameters for supergiants will be presented as a function of spectral type and colours. For the effective temperature versus  ( V − K )₀  colour, we find an empirical calibration with a median deviation of  Δ T = 70 K  in the range of  0.7 < ( V − K )₀ < 5.1  for LC I stars; for LC II, the median deviation is  Δ T = 120 K  from  0.4 < ( V − K )₀ < 4.3  . Effective temperature as a function of spectral type is also calibrated from these data, but shows significantly more scatter than the T _EFF versus  ( V − K )₀  relationship. No deviation of T _EFF versus spectral type is seen for these high-luminosity objects relative to LC II giants. Directly determined diameters range up to  400 R_⊙  , though are limited by poor distance determinations, which dominate the error estimates. These temperature and radii measures reflect a direct calibration of these parameters for supergiants from empirical means.     

The first galaxies: assembly, cooling and the onset of turbulence

We investigate the properties of the first galaxies at   z ≳ 10  with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds  ≃10⁴ K  , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a  ≃5 × 10⁷ M_⊙  system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures  ≳10⁴ K  catalyses the formation of  H₂  and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.     

Boundary layer on the surface of a neutron star

In an attempt to model the accretion on to a neutron star in low-mass X-ray binaries, we present 2D hydrodynamical models of the gas flow in close vicinity of the stellar surface. First, we consider a gas pressure-dominated case, assuming that the star is non-rotating. For the stellar mass we take   M _star= 1.4 × 10⁻² M_⊙  and for the gas temperature   T = 5 × 10⁶ K  . Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of T ≃ 10⁸ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.     

Quasi-stars: accreting black holes inside massive envelopes

We study the structure and evolution of 'quasi-stars', accreting black holes embedded within massive hydrostatic gaseous envelopes. These configurations may model the early growth of supermassive black hole seeds. The accretion rate on to the black hole adjusts so that the luminosity carried by the convective envelope equals the Eddington limit for the total mass,   M _*+ M _BH≈ M _*  . This greatly exceeds the Eddington limit for the black hole mass alone, leading to rapid growth of the black hole. We use analytic models and numerical stellar structure calculations to study the structure and evolution of quasi-stars. We show that the photospheric temperature of the envelope scales as   T _ph∝ M ^−2/5_BH M ^7/20_*  , and decreases with time while the black hole mass increases. Once   T _ph < 10⁴ K  , the photospheric opacity drops precipitously and T _ph hits a limiting value, analogous to the Hayashi track for red giants and protostars, below which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) opacities, this limiting temperature is approximately 4000 K. After a quasi-star reaches this limiting temperature, it is rapidly dispersed by radiation pressure. We find that black hole seeds with masses between 10³ and  10⁴ M_⊙  could form via this mechanism in less than a few Myr.     

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 brightness temperature problem in extreme intra-day variable quasars: a model for PKS 0405-385

I re-examine the brightness temperature problem in PKS 0405-385, which is an extreme intra-day variable radio quasar with an inferred brightness temperature of  ∼5 × 10¹⁴ K  at 5 GHz, well above the Compton catastrophe limit of  ∼10¹¹ K  that is reached when the synchrotron photon energy density exceeds the energy density of the magnetic field. If one takes into account the uncertainty in the distance to the ionized clouds responsible for interstellar scintillation causing rapid intra-day variability in PKS 0405-385, it is possible that the brightness temperature could be as low as  ∼10¹³ K  at 5 GHz, or even lower. The radio spectrum can be fitted by optically thin emission from mono-energetic electrons, or an electron spectrum with a low-energy cut-off such that the critical frequency of the lowest energy electrons is above the radio frequencies of interest. If one observes optically thin emission along a long narrow emission region, the average energy density in the emission region can be many orders of magnitude lower than calculated from the observed intensity if one assumed a spherical emission region. I discuss the physical conditions in the emission region and find that the Compton catastrophe can then be avoided using a reasonable Doppler factor. I also show that MeV to 100-GeV gamma-ray emission at observable flux levels should be expected from extreme intra-day variable sources such as PKS 0405-385.     

The far-infrared signature of dust in high-latitude regions

We present ISOPHOT observations of eight interstellar regions in the 60–200 μm wavelength range. The regions belong to mostly quiescent high-latitude clouds and have optical extinction peaks from   A_V ∼1–6 mag  . From the 150- and 200-μm emission, we derived colour temperatures for the classical big grain component which show a clear trend of decreasing temperature with increasing 200-μm emission. The 200-μm emission per unit   A_V   , however, does not drop at lower temperatures. This fact can be interpreted in terms of an increased far-infrared (FIR) emissivity of the big grains. We developed a two-component model including warm dust with the temperature of the diffuse interstellar medium (ISM) of   T = 17.5 K  , and cold dust with   T = 13.5 K  and FIR emissivity increased by a factor of >4. A mixture of the two components can reproduce the observed colour variations and the ratios   I ₂₀₀/ A_V   and  τ₂₀₀/ A_V   . The relative abundance of small grains with respect to the big grains shows significant variations from region to region at low column densities. However, in lines of sight of higher column density, our data indicate the disappearance of small grains, perhaps a signature of adsorption/coagulation of dust. The larger size and porous structure could also explain the increased FIR emissivity. Our results from eight independent regions suggest that these grains might be ubiquitous in the galactic ISM.