Lyman α radiative transfer in a multiphase medium

Hydrogen Lyman α (Lyα) is our primary emission-line window into high-redshift galaxies. Despite an extensive literature, Lyα radiative transfer in the most realistic case of a dusty, multiphase medium has received surprisingly little detailed theoretical attention. We investigate Lyα resonant scattering through an ensemble of dusty, moving, optically thick gas clumps. We treat each clump as a scattering particle and use Monte Carlo simulations of surface scattering to quantify continuum and Lyα surface scattering angles, absorption probabilities, and frequency redistribution, as a function of the gas dust content. This atomistic approach speeds up the simulations by many orders of magnitude, making possible calculations which are otherwise intractable. Our fitting formulae can be readily adapted for fast radiative transfer in numerical simulations. With these surface scattering results, we develop an analytic framework for estimating escape fractions and line widths as a function of gas geometry, motion, and dust content. Our simple analytic model shows good agreement with full Monte Carlo simulations. We show that the key geometric parameter is the average number of surface scatters for escape in the absence of absorption,     , and we provide fitting formulae for several geometries of astrophysical interest. We consider the following two interesting applications. (i) Equivalent widths ( EWs ). Lyα can preferentially escape from a dusty multiphase interstellar medium if most of the dust lies in cold neutral clouds, which Lyα photons cannot penetrate. This might explain the anomalously high EWs sometimes seen in high-redshift/submillimetre sources. (ii) Multiphase galactic outflows . We show the characteristic profile is asymmetric with a broad red tail, and relate the profile features to the outflow speed and gas geometry. Many future applications are envisaged.     

The expected detection of dust emission from high-redshift Lyman α galaxies

Recent results have shown that a substantial fraction of high-redshift Lyman α (Lyα) galaxies contain considerable amounts of dust. This implies that Lyα galaxies are not primordial, as has been thought in the past. However, this dust has not been directly detected in emission; rather it has been inferred based on extinction estimates from rest-frame ultraviolet (UV) and optical observations. This can be tricky, as both dust and old stars redden galactic spectra at the wavelengths used to infer dust. Measuring dust emission directly from these galaxies is thus a more accurate way to estimate the total dust mass, giving us real physical information on the stellar populations and interstellar medium enrichment. New generation instruments, such as the Atacama Large Millimeter Array and Sub-Millimeter Array, should be able to detect dust emission from some of these galaxies in the submillimeter. Using measurements of the UV spectral slopes, we derive far-infrared flux predictions for of a sample of  23 z ≥ 4  Lyα galaxies. We find that in only a few hours, we can detect dust emission from 39 ± 22 per cent of our Lyα galaxies. Comparing these results to those found from a sample of 21 Lyman break galaxies (LBGs), we find that LBGs are on average 60 per cent more likely to be detected than Lyα galaxies, implying that they are more dusty, and thus indicating an evolutionary difference between these objects. These observations will provide better constraints on dust in these galaxies than those derived from their UV and optical fluxes alone. Undeniable proof of dust in these galaxies could explain the larger than expected Lyα equivalent widths seen in many Lyα galaxies today.     

K-band polarimetry of seven high-redshift radio galaxies

We present the results of K -band imaging polarimetry of seven 3CR radio galaxies with 0.7 <  z  < 1.2. We find strong evidence for polarization in three sources: 3C 22, 3C 41 and 3C 114. Of these, 3C 41 shows strong evidence of having a quasar core the infrared light of which is scattered by dust. We also find some evidence for polarization in 3C 54 and in 3C 356. The two point-like sources (3C 22 and 3C 41) and the barely-elongated 3C 54 appear to have of order 10 per cent of their K -band flux contributed by scattered light from the active nucleus. We conclude that scattered nuclear light can form a significant component of the near-infrared light emitted by high-redshift radio galaxies, and discuss models in which the scattering particles are electrons and dust grains.     

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.     

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  .     

Probing the neutral fraction of the IGM with GRBs during the epoch of reionization

We show that near-infrared observations of the red side of the Lyα line from a single gamma-ray burst (GRB) afterglow cannot be used to constrain the global neutral fraction of the intergalactic medium (IGM),     , at the GRB's redshift to better than     . Some GRB sightlines will encounter more neutral hydrogen than others at fixed     owing to the patchiness of reionization. GRBs during the epoch of reionization will often bear no discernible signature of a neutral IGM in their afterglow spectra. We discuss the constraints on     from the   z = 6.3  burst, GRB050904, and quantify the probability of detecting a neutral IGM using future spectroscopic observations of high-redshift, near-infrared GRB afterglows. Assuming an observation with signal-to-noise ratio similar to the Subaru FOCAS spectrum of GRB050904 and that the column density distribution of damped Lyα absorbers is the same as measured at lower redshifts, a GRB from an epoch when     can be used to detect a partly neutral IGM at 97 per cent confidence level ≈10 per cent of the time (and, for an observation with three times the sensitivity, ≈30 per cent of the time).     

Understanding the redshift evolution of the luminosity functions of Lyman α emitters

We present a semi-analytical model of star formation which explains simultaneously the observed ultraviolet (UV) luminosity function (LF) of high-redshift Lyman break galaxies (LBGs) and LFs of Lyman α emitters. We consider both models that use the Press–Schechter (PS) and Sheth–Tormen (ST) halo mass functions to calculate the abundances of dark matter haloes. The Lyman α LFs at   z ≲ 4  are well reproduced with only ≲10 per cent of the LBGs emitting Lyman α lines with rest equivalent width greater than the limiting equivalent width of the narrow band surveys. However, the observed LF at   z > 5  can be reproduced only when we assume that nearly all LBGs are Lyman α emitters. Thus, it appears that  4 < z < 5  marks the epoch when a clear change occurs in the physical properties of the high-redshift galaxies. As Lyman α escape depends on dust and gas kinematics of the interstellar medium (ISM), this could mean that on an average the ISM at   z > 5  could be less dusty, more clumpy and having more complex velocity field. All of these will enable easier escape of the Lyman α photons. At   z > 5  , the observed Lyman α LF are well reproduced with the evolution in the halo mass function along with very minor evolution in the physical properties of high-redshift galaxies. In particular, up to   z = 6.5  , we do not see the effect of evolving intergalactic medium opacity on the Lyman α escape from these galaxies.     

Modeling the Lyα radiation of high-redshift galaxies

The Lyα line emission of high-redshift galaxies depends on the density and temperature distribution of the gas, the kinematics and the dust content. We use a finite element method to model theLyα radiation of different 3D configurations considering complete frequency redistribution and the influence of velocity fields. Our results show that the central absorption feature of the double-peaked Lyα line profile observed in many radio galaxies with z=2-4 is probably the consequence of frequency redistribution rather than foreground absorption. The blue peak of the profile is enhanced for models with in fall motion and the red peak for models with outflow motion. In particular, we attempt to model the extendedLyα emission of high-redshift radio galaxies, where we consider results of corresponding hydrodynamical simulations to select possible model configurations. We find that Lyα photons scattered outside a jet-influenced low-density region are able to produce an extended Lyα halo. This revised version was published online in August 2006 with corrections to the Cover Date.     

An 80-kpc Lyα halo around a high-redshift type-2 quasi-stellar object

We announce the discovery of an extended emission-line region associated with a high-redshift type-2 quasi-stellar object (QSO). The halo, which was discovered in our new wide-field narrow-band survey, resides at   z = 2.85  in the Spitzer First Look Survey region and is extended over ∼80 kpc. Deep very long baseline interferometry (VLBI) observations imply that approximately 50 per cent of the radio emission is extended on scales >200 pc. The inferred active galactic nuclei (AGN) luminosity is sufficient to ionize the extended halo, and the optical emission is consistent with being triggered coevally with the radio source. The Lyα halo is as luminous as those found around high-redshift radio galaxies; however, the active nucleus is several orders of magnitude less luminous at radio wavelengths than those Fanarof–Riley type II (FRIIs) more commonly associated with extended emission-line regions. AMS05 appears to be a high-redshift analogue to the radio-quiet quasar E1821+643 which is core dominated, but which also exhibits extended Fanarof–Riley type I (FRI)-like structure and contains an optically powerful AGN. We also find evidence for more quiescent kinematics in the Lyα emission line in the outer regions of the halo, reminiscent of the haloes around the more powerful FRIIs. The optical to mid-infrared spectral energy distribution is well described by a combination of an obscured QSO  ( L _bol∼ 3.4 ± 0.2 × 10¹³ L_⊙)  and a 1.4 Gyr old simple stellar population with mass  ∼3.9 ± 0.3 × 10¹¹ M_⊙  .     

A possible EGRET source with a two-sided radio jet structure on a parsec scale – 3EG J0808+5114

Recent observations have revealed that damped Lyα clouds (DLAs) host star formation activity. In order to examine if such star formation activity can be triggered by ionization fronts, we perform high-resolution hydrodynamics and radiative transfer simulations of the effect of radiative feedback from propagating ionization fronts on high-density clumps. We examine two sources of ultraviolet (UV) radiation field to which high-redshift ( z ∼ 3) galaxies could be exposed: one corresponding to the UV radiation originating from stars within the DLA, itself, and the other corresponding to the UV background radiation. We find that, for larger clouds, the propagating I-fronts created by local stellar sources can trigger cooling instability and collapse of significant part, up to 85 per cent, of the cloud, creating conditions for star formation in a time-scale of a few Myr. The passage of the I-front also triggers collapse of smaller clumps (with radii below ∼4 pc), but in these cases the resulting cold and dense gas does not reach conditions conducive to star formation. Assuming that 85 per cent of the gas initially in the clump is converted into stars, we obtain a star formation rate of  ∼0.25 M_⊙ yr⁻¹ kpc⁻²  . This is somewhat higher than the value derived from recent observations. On the other hand, the background UV radiation which has harder spectrum fails to trigger cooling and collapse. Instead, the hard photons which have long mean free-path heat the dense clumps, which as a result expand and essentially dissolve in the ambient medium. Therefore, the star formation activity in DLAs is strongly regulated by the radiative feedback, both from the external UV background and internal stellar sources and we predict quiescent evolution of DLAs (not starburst-like evolution).     

Fluorescent Lyman α emission from gas near a QSO at redshift 4.28

We use integral field spectroscopy with the Gemini North Telescope to detect probable fluorescent Lyman α (Lyα) emission from gas lying close to the luminous QSO PSS 2155+1358 at redshift 4.28. The emission is most likely coming not from primordial gas, but from a multiphase, chemically enriched cloud of gas lying about 50 kpc from the QSO. It appears to be associated with a highly ionized associated absorber seen in the QSO spectrum. With the exception of this gas cloud, the environment of the QSO is remarkably free of neutral hydrogen. We also marginally detect Lyα emission from a foreground subdamped Lyα absorption-line system.     

The obscured BLR in the radio galaxy 3C 234

We present spectropolarimetric observations of the radio galaxy 3C 234 at optical and, for the first time, near-IR wavelengths. In agreement with previous observations we detect broad Hα in total and polarized flux, consistent with the scenario in which the central active nucleus is hidden from view and observed via scattered radiation.
We model the flux and polarization properties at optical and near-IR wavelengths, which result in a point-source, cone-based scattering geometry in the optical, and an extended source scattering in the near-IR, with a dichroic view to the emission regions, which becomes important only at the longer wavelengths. From this model we calculate an intrinsic Hα luminosity of 4.9×10⁴⁴ erg s⁻¹, and an extinction to the near-IR emission region of 60 mag for A_V . The scattered radiation also undergoes extinction, this time by 1 mag for A_V , assuming a source function of a typical stellar-subtracted Seyfert 1 galaxy. Our modelling does not require the presence of a second featureless continuum (FC2).
Additionally, we require that the scatterers are in bulk radial outflow at approximately 600 km s⁻¹ to explain an observed increase in the intrinsic polarization of the scattered broad Hα line in the red wing, and blueshifted narrow-line components in polarized flux.     

The profiles of polarized broad Hα lines in radio galaxies

We present spectropolarimetric observations of seven broad-line radio galaxies. We find significantly polarized broad Hα emission in four objects including two, Arp 102B and 3C 390.3, which have double-peaked line profiles. In these objects the prominent redshifted and blueshifted peaks of the broad Hα line have no clear counterparts in polarized flux. This conflicts with theoretical predictions for a relativistic line-emitting disc with an electron scattering atmosphere, one of the leading models advanced to account for the double-peaked lines. The shapes and widths of the polarized line profiles can be explained if, as expected in unified schemes, the scattering occurs near the poles of an obscuring torus. However, the observed polarization position angles favour geometries in which the scattering plane is perpendicular to the radio jet. A configuration in which Hα photons emitted by a biconical flow are scattered off the inner wall of the torus has this property, and would also produce a single-peaked scattered line profile. With the exception of 3C 227, the sample as a whole conforms to the general trend in powerful radio galaxies for the optical polarization to be aligned with the radio source axis, favouring toroidal rather than polar scattering.     

Contribution of galaxies to the background hydrogen-ionizing flux

We estimate the evolution of the contribution of galaxies to the cosmic background flux at 912 Å by means of a semi-analytic model of galaxy formation and evolution. Such modelling has been quite successful in reproducing the optical properties of galaxies. We assume that high-redshift damped Lyman α systems are the progenitors of present-day galaxies, and we design a series of models that are consistent with the evolution of cosmic comoving emissivities in the available near-infrared, optical, ultraviolet and far-infrared bands along with the evolution of the neutral hydrogen content and average metallicity of damped Lyman α systems. We use these models to compute the galactic contribution to the Lyman-limit emissivity and background flux for 0 ≃  z  ≤ 4. We take into account the absorption of Lyman-limit photons by H  I and dust in the interstellar medium of the galaxies. We find that the background Lyman-limit flux due to galaxies might dominate (or be comparable to) the contribution from quasars at almost all redshifts if the absorption by H  I in the interstellar medium is neglected. Such H  I absorption would result in a severe diminishing of this flux — by almost three orders of magnitude at high redshifts and by one to two orders at z  ≃ 0. Though the resulting galaxy flux is completely negligible at high redshifts, it is comparable to the quasar flux at z  ≃ 0.     

Damped Lyman-α Systems

Damped Lyman-α systems observed in the spectra of high-redshift quasars are considered as the progenitors of present-day galaxies.Indeed, the large neutral hydrogen column densities observed and the presence of metals imply that the gas is somehow closely associated with regions of star formation. The nature of the absorbing object is unclear however.The discussion of whether high-redshift damped Lyman-α systems are produced throughlarge, fast-rotating, thick disks or through building-blocks of galaxies is important since it is related to how present-day galaxies form, either through initial formation of large disks and subsequent accretion of gas or as a result of merging of pregalactic clumps.We argue here that damped Lyman-α systems are progenitors of any kind of galaxy (dwarfs and large spirals); they trace peaks in the spatial density distribution of the gas. This revised version was published online in July 2006 with corrections to the Cover Date.     

Physical effects on the Lyα forest flux power spectrum: damping wings, ionizing radiation fluctuations and galactic winds

We explore several physical effects on the power spectrum of the Lyα forest transmitted flux. The effects we investigate here are not usually part of hydrodynamic simulations and so need to be estimated separately. The most important effect is that of high column density absorbers with damping wings, which add power on large scales. We compute their effect using the observational constraints on their abundance as a function of column density. Ignoring their effect leads to an underestimation of the slope of the linear theory power spectrum. The second effect we investigate is that of fluctuations in the ionizing radiation field. For this purpose we use a very large high-resolution N -body simulation, which allows us to simulate both the fluctuations in the ionizing radiation and the small-scale Lyα forest within the same simulation. We find an enhancement of power on large scales for quasars and a suppression for galaxies. The strength of the effect rapidly increases with increasing redshift, allowing it to be uniquely identified in cases where it is significant. We develop templates that can be used to search for this effect as a function of quasar lifetime, quasar luminosity function and attenuation length. Finally, we explore the effects of galactic winds using hydrodynamic simulations. We find the wind effects on the Lyα forest power spectrum to be degenerate with parameters related to the temperature of the gas that are already marginalized over in cosmological fits. While more work is needed to conclusively exclude all possible systematic errors, our results suggest that, in the context of data analysis procedures, where parameters of the Lyα forest model are properly marginalized over, the flux power spectrum is a reliable tracer of cosmological information.     

The probability distribution function of the Lyman α transmitted flux from a sample of Sloan Digital Sky Survey quasars

We present a measurement of the probability distribution function (PDF) of the transmitted flux in the Lyman α (Lyα) forest from a sample of 3492 quasars included in the Sloan Digital Sky Survey data release 3 (SDSS DR3). Our intention is to investigate the sensitivity of the Lyα flux PDF as measured from low-resolution and low signal-to-noise ratio data to a number of systematic errors such as uncertainties in the mean flux, continuum and noise estimate. The quasar continuum is described by the superposition of a power law and emission lines. We perform a power-law continuum fitting on a spectrum-by-spectrum basis, and obtain an average continuum slope of  α_ν= 0.59 ± 0.36  in the redshift range  2.5 < z < 3.5  . We take into account that the variation in the continuum indices increases the mean flux by 3 and 7 per cent at   z = 3  and 2.4, respectively, as compared to the values inferred with a single (mean) continuum slope. We compare our measurements to the PDF obtained with mock lognormal spectra, whose statistical properties have been constrained to match the observed Lyα flux PDF and power spectrum of high-resolution data. Using our power-law continuum fitting and the SDSS pipeline noise estimate yields a poor agreement between the observed and mock PDFs. Allowing for a break in the continuum slope and, more importantly, for residual scatter in the continuum level substantially improves the agreement. A decrease of ∼10–15 per cent in the mean quasar continuum with a typical rms variance at the 20 per cent level can account for the data, provided that the noise excess correction is no larger than ≲10 per cent.     

A method to measure the mass of damped Lyα absorber host galaxies using fluctuations in 21-cm emission

Observations of damped Lyα absorbers (DLAs) indicate that the fraction of hydrogen in its neutral form (H  i ) is significant by mass at all redshifts. This gas represents the reservoir of material that is available for star formation at late times. As a result, observational identification of the systems in which this neutral hydrogen resides is an important missing ingredient in models of galaxy formation. Precise identification of DLA host mass via traditional clustering studies is not practical owing to the small numbers of known systems being spread across sparsely distributed sightlines. However, following the completion of re-ionization, 21-cm surface brightness fluctuations will be dominated by neutral hydrogen in DLAs. No individual DLAs could be detected in 21-cm emission. Rather, observations of these fluctuations will measure the combined clustering signal from all DLAs within a large volume. We show that measurement of the spherically averaged power spectrum of 21-cm intensity fluctuations due to DLAs could be used to measure the galaxy bias for DLA host galaxies when combined with an independent measurement of the cosmological H  i mass density from quasar absorption studies. Utilizing this technique, the low-frequency arrays now under construction could measure the characteristic DLA host mass with a statistical precision as low as 0.3 dex at z ≳ 4. In addition, high signal-to-noise ratio observations of the peculiar-motion-induced anisotropy of the power spectrum would facilitate measurement of both the DLA host mass and the cosmic H  i density directly from 21-cm fluctuations. By exploiting this anisotropy, a second generation of low-frequency arrays with an order of magnitude increase in collecting area could measure the values of cosmic H  i density and DLA host mass, with uncertainties of a few per cent and a few tens of per cent, respectively.     

21-cm fluctuations from inhomogeneous X-ray heating before reionization

Many models of early structure formation predict a period of heating immediately preceding reionization, when X-rays raise the gas temperature above that of the cosmic microwave background. These X-rays are often assumed to heat the intergalactic medium (IGM) uniformly, but in reality will heat the gas more strongly closer to the sources. We develop a framework for calculating fluctuations in the 21-cm brightness temperature that originate from this spatial variation in the heating rate. High-redshift sources are highly clustered, leading to significant gas temperature fluctuations (with fractional variations ∼40 per cent, peaking on   k ∼ 0.1 Mpc⁻¹  scales). This induces a distinctive peak-trough structure in the angle-averaged 21-cm power spectrum, which may be accessible to the proposed Square Kilometre Array. This signal reaches the ∼10 mK level, and is stronger than that induced by Lyα flux fluctuations. As well as probing the thermal evolution of the IGM before reionization, this 21-cm signal contains information about the spectra of the first X-ray sources. Finally, we consider disentangling temperature, density and Lyα flux fluctuations as functions of redshift.     

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.