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⁻¹.     

The H i opacity of the intergalactic medium at redshifts 1.6 < z < 3.2

We use high-quality echelle spectra of 24 quasi-stellar objects to provide a calibrated measurement of the total amount of Lyα forest absorption (DA) over the redshift range  2.2 < z < 3.2  . Our measurement of DA excludes absorption from metal lines or the Lyα lines of Lyman-limit systems and damped Lyα systems. We use artificial spectra with realistic flux calibration errors to show that we are able to place continuum levels that are accurate to better than 1 per cent. When we combine our results with our previous results between  1.6 < z < 2.2  , we find that the redshift evolution of DA is well described over  1.6 < z < 3.2  as   A (1 + z )^γ  , where   A = 0.0062  and  γ= 2.75  . We detect no significant deviations from a smooth power-law evolution over the redshift range studied. We find less H  i absorption than expected at   z = 3  , implying that the ultraviolet background is  ∼40  per cent higher than expected. Our data appears to be consistent with an H  i ionization rate of  Γ∼ 1.4 × 10⁻¹² s⁻¹  .     

The relation between Lyman α absorbers and gas-rich galaxies in the local Universe

We use high-resolution hydrodynamical simulations to investigate the spatial correlation between weak  ( N _{H  i} < 10¹⁵ cm⁻²)  Lyα absorbers and gas-rich galaxies in the local Universe. We confirm that Lyα absorbers are preferentially expected near gas-rich galaxies and that the degree of correlation increases with the column density of the absorber. The real-space galaxy auto-correlation is stronger than the cross-correlation (correlation lengths   r _0,gg= 3.1 ± 0.1 Mpc  h ⁻¹  and   r _0,ag= 1.4 ± 0.1 Mpc  h ⁻¹  , respectively), in contrast with the recent results of Ryan-Weber, and the auto-correlation of absorbers is very weak. These results are robust to the presence of strong galactic winds in the hydrodynamical simulations. In redshift space, a further mismatch arises since at small separations the distortion pattern of the simulated galaxy–absorber cross-correlation function is different from the one measured by Ryan-Weber. However, when sampling the intergalactic medium along a limited number of lines-of-sight, as in the real data, uncertainties in the cross-correlation estimates are large enough to account for these discrepancies. Our analysis suggests that the statistical significance of difference between the cross-correlation and auto-correlation signal in current data sets is ∼1σ only.     

Semi-analytic simulations of galactic winds: volume filling factor, ejection of metals and parameter study

We present a semi-analytic treatment of galactic winds within high-resolution, large-scale cosmological N -body simulations of a Λ cold dark matter (ΛCDM) universe. The evolution of winds is investigated by following the expansion of supernova-driven superbubbles around the several hundred thousand galaxies that form in an approximately spherical region of space with diameter 52  h ⁻¹ Mpc and mean density close to the mean density of the universe. We focus our attention on the impact of winds on the diffuse intergalactic medium. Initial conditions for mass loss at the base of winds are taken from Shu, Mo & Mao. Results are presented for the volume filling factor and the mass fraction of the intergalactic medium (IGM) affected by winds, and their dependence on the model parameters is carefully investigated. The mass-loading efficiency of bubbles is a key factor to determine the evolution of winds and their global impact on the IGM: the higher the mass loading, the later the IGM is enriched with metals. Galaxies with 10⁹ < M _★ < 10¹⁰ M_⊙ are responsible for most of the metals ejected into the IGM at   z = 3  , while galaxies with   M _★ < 10⁹ M_⊙   give a non-negligible contribution only at higher redshifts, when larger galaxies have not yet assembled. We find a higher mean IGM metallicity than Lyα forest observations suggest, and we argue that the discrepancy may be explained by the high temperatures of a large fraction of the metals in winds, which may not leave detectable imprints in absorption in the Lyα forest.     

ZEN2: a narrow J-band search for z∼ 9 Lyα emitting galaxies directed towards three lensing clusters

We present the results of a continuing survey to detect Lyα emitting galaxies at redshifts   z ∼ 9  : the ' z equals nine' (ZEN) survey. We have obtained deep VLT Infrared Spectrometer and Array Camera observations in the narrow J -band filter NB119 directed towards three massive lensing clusters: Abell clusters 1689, 1835 and 114. The foreground clusters provide a magnified view of the distant Universe and permit a sensitive test for the presence of very high redshift galaxies. We search for   z ∼ 9 Lyα  emitting galaxies displaying a significant narrow-band excess relative to accompanying J -band observations that remain undetected in Hubble Space Telescope ( HST )/Advanced Camera for Surveys (ACS) optical images of each field. No sources consistent with this criterion are detected above the unlensed 90 per cent point-source flux limit of the narrow-band image,   F _NB= 3.7 × 10⁻¹⁸ erg s⁻¹ cm⁻²  . To date, the total coverage of the ZEN survey has sampled a volume at   z ∼ 9  of approximately 1700 comoving Mpc³ to a Lyα emission luminosity of  10⁴³ erg s⁻¹  . We conclude by considering the prospects for detecting   z ∼ 9 Lyα  emitting galaxies in light of both observed galaxy properties at   z < 7  and simulated populations at   z > 7  .     

The evolution of the Lyα forest effective optical depth following He ii reionization

Three independent observational studies have now detected a narrow  (Δ z ≃ 0.5)  dip centred at   z = 3.2  in the otherwise smooth redshift evolution of the Lyα forest effective optical depth. This feature has previously been interpreted as an indirect signature of rapid photoheating in the intergalactic medium (IGM) during the epoch of He  ii reionization. We examine this interpretation using a semi-analytic model of inhomogeneous He  ii reionization and high-resolution hydrodynamical simulations of the Lyα forest. We instead find that a rapid  (Δ z ≃ 0.2)  boost to the IGM temperature  (Δ T ≃ 10⁴ K)  beginning at   z = 3.4  produces a well understood and generic evolution in the Lyα effective optical depth, where a sudden reduction in the opacity is followed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low-density IGM. This behaviour is inconsistent with the narrow feature in the observational data. If photoheating during He  ii reionization is instead extended over several redshift units, as recent theoretical studies suggest, then the Lyα opacity will evolve smoothly with redshift. We conclude that the sharp dip observed in the Lyα forest effective optical depth is instead most likely due to a narrow peak in the hydrogen photoionization rate around   z = 3.2  , and suggest that it may arise from the modulation of either reprocessed radiation during He  ii reionization, or the opacity of Lyman limit systems.     

Lyα absorption in the nearby Universe: the sightline to Q1821+643

We present the first results of a survey designed to understand the origin of Lyα-forest absorption systems at low redshift. Using the WYFFOS and HYDRA multifibre spectrographs on the William Herschel and Wisconsin Indiana Yale NOAO (WIYN) telescopes, we have identified 51 galaxies brighter than b _j := :18.5 within 30 arcmin of the sightline of the QSO 1821+643. We find three galaxies within 500 h ⁻¹ kpc of the QSO sightline; the nearest galaxy is 104 h ⁻¹ kpc away from the line of sight, and is at the same redshift as a strong ( W _r :=:0.63Å) Lyα absorption line. The remaining two galaxies have no corresponding absorption to extremely low equivalent width limits (< 0.05 Å). Beyond 500 h ⁻¹ kpc, Lyα absorption lines are found at redshifts similar to those of several galaxies, but we show that these coincidences are likely to be accidental.   Half of the Lyα systems for which we could have found at least an L ^* galaxy have no galaxies at the redshifts of the absorbers. For the majority of the remainder, we show that any apparent association with galaxies is probably coincidental. These Lyα systems are characterized by their weak equivalent widths ( W _r :<:0.2Å), and we conclude that this population of absorbers is uncorrelated, or at best weakly correlated, with galaxies.     

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.     

Glimpsing through the high-redshift neutral hydrogen fog

We analyse the transmitted flux in a sample of 17 QSOs spectra at 5.74 ≤ z _em≤ 6.42 to obtain tighter constraints on the volume-averaged neutral hydrogen fraction, x _{H  i} , at z ≈ 6. We study separately the narrow transmission windows (peaks) and the wide dark portions (gaps) in the observed absorption spectra. By comparing the statistics of these spectral features with a semi-analytical model of the Lyα forest, we conclude that x _{H  i} evolves smoothly from 10^−4.4 at   z = 5.3  to 10^−4.2 at   z = 5.6  , with a robust upper limit x _{H  i} < 0.36 at   z = 6.3  . The frequency and physical sizes of the peaks imply an origin in cosmic underdense regions and/or in H  ii regions around faint quasars or galaxies. In one case (the intervening H  ii region of the faint quasar RD J1148+5253 at   z = 5.70  along the line of sight of SDSS J1148+5251 at   z = 6.42  ) the increase of the peak spectral density is explained by the first-ever detected transverse proximity effect in the H  i Lyα forest; this indicates that at least some peaks result from a locally enhanced radiation field. We then obtain a strong lower limit on the foreground QSO lifetime of t _Q > 11 Myr. The observed widths of the peaks are found to be systematically larger than the simulated ones. Reasons for such discrepancy might reside either in the photoionization equilibrium assumption or in radiative transfer effects.     

The matter power spectrum from the Lyα forest: an optical depth estimate

We measure the matter power spectrum from 31 Lyα spectra spanning the redshift range of 1.6–3.6. The optical depth, τ, for Lyα absorption of the intergalactic medium is obtained from the flux using the inversion method of Nusser & Haehnelt. The optical depth is converted to density by using a simple power-law relation,  τ∝ (1 +δ)^α  . The non-linear 1D power spectrum of the gas density is then inferred with a method that makes simultaneous use of the one- and two-point statistics of the flux and compared against theoretical models with a likelihood analysis. A cold dark matter model with standard cosmological parameters fits the data well. The power-spectrum amplitude is measured to be (assuming a flat Universe),  σ₈= (0.92 ± 0.09) × (Ω_m/0.3)^−0.3  , with α varying in the range of 1.56–1.8 with redshift. Enforcing the same cosmological parameters in all four redshift bins, the likelihood analysis suggests some evolution in the temperature–density relation and the thermal smoothing length of the gas. The inferred evolution is consistent with that expected if reionization of He  ii occurred at   z ∼ 3.2  . A joint analysis with the Wilkinson Microwave Anisotropy Probe results together with a prior on the Hubble constant as suggested by the Hubble Space Telescope key project data, yields values of Ω_m and σ₈ that are consistent with the cosmological concordance model. We also perform a further inversion to obtain the linear 3D power spectrum of the matter density fluctuations.     

The deuterium abundance in the z= 0.7 absorber towards QSO PG1718+4807

We report on a further analysis of the ratio of deuterium to hydrogen (D/H) using Hubble Space Telescope ( HST ) spectra of the   z = 0.701  Lyman limit system towards the quasi-stellar object (QSO) PG1718+481. Initial analyses of this absorber found it gave a high D/H value, 1.8–  3.1 × 10⁻⁴  , inconsistent with several higher redshift measurements. It is thus important to critically examine this measurement. By analysing the velocity widths of the D  i , H  i and metal lines present in this system, Kirkman et al. report that the additional absorption in the blue wing of the Lyα line cannot be D  i , with a confidence level of 98 per cent. Here we present a more detailed analysis, taking into account possible wavelength shifts between the three sets of HST spectra used in the analysis. We find that the constraints on this system are not as strong as those claimed by Kirkman et al. The discrepancy between the parameters of the blue wing absorption and the parameters expected for D  i is marginally worse than 1σ.
Tytler et al. commented on the first analysis of Webb et al., reporting the presence of a contaminating lower redshift Lyman limit system, with  log[ N (H  i )] = 16.7  at   z = 0.602  , which biases the N (H  i ) estimate for the main system. Here we show that this absorber actually has  log[ N (H  i )] < 14.6  and does not impact on the estimate of N (H  i ) in the system of interest at   z = 0.701  .
The purpose of the present paper is to highlight important aspects of the analysis which were not explored in previous studies, and hence to help refine the methods used in future analyses of D/H in quasar spectra.     

On the polarization of Hα scattered by neutral hydrogen in active galactic nuclei

Raman scattering by atomic hydrogen converts the UV continuum around Lyβ into optical continuum around Hα, and the basic atomic physics has been discussed in several works on symbiotic stars. We propose that the same process may operate in active galactic nuclei (AGN) and calculate the linear polarization of the broad emission lines Raman-scattered by high-column neutral hydrogen component. The conversion efficiency of the Raman scattering process is discussed and the expected scattered flux is computed using the spectral energy distribution of an AGN given by a typical power law. The high-column H  i component in AGN is suggested by many observations, encompassing the radio through UV and X-ray ranges.   When neutral hydrogen component with a column density ∼10²² cm⁻² is present around the active nucleus, it is found that the scattered Hα is characterized by a very broad width ∼20 000 km s⁻¹ and that the strength of the polarized flux is comparable to that of the electron-scattered flux expected from a conventional unified model of narrow-line AGN. The width of the scattered flux is mainly determined by the column density of the neutral scatterers where the total scattering optical depth becomes of order unity. The asymmetry in the Raman scattering cross-section around Lyβ introduces red asymmetric polarized profiles around Hα. The effects of the blended Lyβ and O  vi 1034 doublet are also investigated.   We briefly discuss the spectropolarimetric observations performed on the Seyfert galaxy IRAS 110548-1131 and the narrow line radio galaxy Cyg A. Several predictions regarding the scattering by the high-column neutral hydrogen component in AGN are discussed.     

Seismological studies of ZZ Ceti stars – II. Application to the ZZ Ceti class

We used the detected pulsation modes and adiabatic pulsation models to do seismology of the class of ZZ Ceti stars and measure the H layer mass for 83 stars. We found the surface hydrogen layer to be within the range  10^−9.5≥ M _H/ M _*≥ 10⁻⁴  , with an average of   M _H/ M _*= 10^−6.3  , which is thinner than the predicted value of   M _H/ M _*= 10⁻⁴  , indicating that the stars lose more mass during their evolution than previously expected. These results are preliminary and do not include the possible effects of realistic C/O profiles on the fits.     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

Ly+ and ultraviolet emission from high-redshift gamma-ray burst hosts: to what extent do gamma-ray bursts trace star formation?

We report the result of a search for Lyα emission from the host galaxies of the gamma-ray bursts  (GRBs) 030226 ( z = 1.986), 021004 ( z = 2.335)  and  020124 ( z = 3.198)  . We find that the host galaxy of GRB 021004 is an extended (around 8 kpc) strong Lyα emitter with a rest-frame equivalent width (EW) of 68⁺¹²₋₁₁Å, and a star formation rate of  10.6 ± 2.0 M_⊙ yr⁻¹  . We do not detect the hosts of GRB 030226 and GRB 020124, but the upper limits on their Lyα fluxes do not rule out large rest-frame EWs. In the fields of GRB 021004 and GRB 030226 we find seven and five other galaxies, respectively, with excess emission in the narrow-band filter. These galaxies are candidate Lyα-emitting galaxies in the environment of the host galaxies. We have also compiled a list of all   z ≳ 2  GRB hosts, and demonstrate that a scenario where they trace star formation in an unbiased way is compatible with current observational constraints. Fitting the   z = 3  luminosity function (LF) under this assumption results in a characteristic luminosity of   R *= 24.6  and a faint-end slope of  α=−1.55  , consistent with the LF measured for Lyman-break galaxies.     

Revision of VLT/UVES constraints on a varying fine-structure constant

We critically review the current null results on a varying fine-structure constant, α, derived from Very Large Telescope (VLT)/Ultraviolet and Visual Echelle Spectrograph (UVES) quasar absorption spectra, focusing primarily on the many-multiplet analysis of 23 absorbers from which Chand et al. reported a weighted mean relative variation of  Δα/α= (−0.06 ± 0.06) × 10⁻⁵  . Our analysis of the same reduced data , using the same fits to the absorption profiles , yields very different individual  Δα/α  values with uncertainties typically larger by a factor of ∼3. We attribute the discrepancies to flawed parameter estimation techniques in the original analysis and demonstrate that the original  Δα/α  values were strongly biased towards zero. Were those flaws not present, the input data and spectra should have given a weighted mean of  Δα/α= (−0.44 ± 0.16) × 10⁻⁵  . Although this new value does reflect the input spectra and fits (unchanged from the original work – only our analysis is different), we do not claim that it supports previous Keck/High Resolution Echelle Spectrograph (HIRES) evidence for a varying α: there remains significant scatter in the individual  Δα/α  values which may stem from the overly simplistic profile fits in the original work. Allowing for such additional, unknown random errors by increasing the uncertainties on  Δα/α  to match the scatter provides a more conservative weighted mean,  Δα/α= (−0.64 ± 0.36) × 10⁻⁵  . We highlight similar problems in other current UVES constraints on varying α and argue that comparison with previous Keck/HIRES results is premature.     

The local star formation rate and radio luminosity density

We present a new determination of the local volume-averaged star formation rate from the 1.4-GHz luminosity function of star forming galaxies. Our sample, taken from the   B ≤12  Revised Shapley–Ames catalogue (231 normal spiral galaxies over an effective area of 7.1 sr) has ≃100 per cent complete radio detections and is insensitive to dust obscuration and cirrus contamination. After removal of known active galaxies, the best-fitting Schechter function has a faint-end slope of  −1.27±0.07  in agreement with the local H α luminosity function, characteristic luminosity   L ∗=(2.6±0.7)×10²² W Hz⁻¹  and density   φ ∗=(4.8±1.1)×10⁻⁴ Mpc⁻³.  The inferred local radio luminosity density of  (1.73±0.37±0.03)×10¹⁹ W Hz⁻¹ Mpc⁻³  (Poisson noise, large-scale structure fluctuations) implies a volume-averaged star formation rate ∼2 times larger than the Gallego et al. H α estimate, i.e.   ρ _1.4 GHz=(2.10±0.45±0.04)×10⁻² M_⊙ yr⁻¹ Mpc⁻³  for a Salpeter initial mass function from  0.1–125 M_⊙  and Hubble constant of 50 km s⁻¹ Mpc⁻¹. We demonstrate that the Balmer decrement is a highly unreliable extinction estimator, and argue that optical–ultraviolet (UV) star formation rates (SFRs) are easily underestimated, particularly at high redshift.     

An excess of damped Lyman α galaxies near quasi-stellar objects

We present a sample of 33 damped Lyman α systems (DLAs) discovered in the Sloan Digital Sky Survey (SDSS) whose absorption redshifts ( z _abs) are within 6000 km s⁻¹ of the quasi-stellar object's (QSO) systemic redshift ( z _sys). Our sample is based on  731 2.5 < z _sys < 4.5  non-broad absorption line (non-BAL) QSOs from Data Release 3 (DR3) of the SDSS. We estimate that our search is ≈100 per cent complete for absorbers with N (H  i )  ≥ 2 × 10²⁰ cm⁻²  . The derived number density of DLAs per unit redshift, n ( z ), within  Δ v < 6000 km s⁻¹  is higher (3.5σ significance) by almost a factor of 2 than that of intervening absorbers observed in the SDSS DR3, i.e. there is evidence for an overdensity of galaxies near the QSOs. This provides a physical motivation for excluding DLAs at small velocity separations in surveys of intervening 'field' DLAs. In addition, we find that the overdensity of proximate DLAs is independent of the radio-loudness of the QSO, consistent with the environments of radio-loud and radio-quiet QSOs being similar.     

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.     

AMI limits on 15-GHz excess emission in northern H ii regions

We present observations between 14.2 and 17.9 GHz of 16 Galactic H  ii regions made with the Arcminute Microkelvin Imager. In conjunction with data from the literature at lower radio frequencies we investigate the possibility of a spinning dust component in the spectra of these objects. We conclude that there is no significant evidence for spinning dust towards these sources and measure an average spectral index of  α= 0.15 ± 0.07  (where   S ∝ν^−α  ) between 1.4 and 17.9 GHz for the sample.