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.     

The polarization of scattered Lyα radiation around high-redshift galaxies

The high-redshift Universe contains luminous Lyα emitting sources such as galaxies and quasars. The emitted Lyα radiation is often scattered by surrounding neutral hydrogen atoms. We show that the scattered Lyα radiation obtains a high level of polarization for a wide range of likely environments of high-redshift galaxies. For example, the backscattered Lyα flux observed from galaxies surrounded by a superwind-driven outflow may reach a fractional polarization as high as ∼40 per cent. Equal levels of polarization may be observed from neutral collapsing protogalaxies. Resonant scattering in the diffuse intergalactic medium typically results in a lower polarization amplitude (≲7 per cent), which depends on the flux of the ionizing background. Spectral polarimetry can differentiate between Lyα scattering off infalling gas and outflowing gas; for an outflow, the polarization should increase towards longer wavelengths while for infall the opposite is true. Our numerical results suggest that Lyα polarimetry is feasible with existing instruments, and may provide a new diagnostic of the distribution and kinematics of neutral hydrogen around high-redshift galaxies. Moreover, polarimetry may help suppress infrared lines originating in the Earth's atmosphere, and thus improve the sensitivity of ground-based observations to high-redshift Lyα emitting galaxies outside the currently available redshift windows.     

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.     

Cosmic dynamics in the era of Extremely Large Telescopes

The redshifts of all cosmologically distant sources are expected to experience a small, systematic drift as a function of time due to the evolution of the Universe's expansion rate. A measurement of this effect would represent a direct and entirely model-independent determination of the expansion history of the Universe over a redshift range that is inaccessible to other methods. Here we investigate the impact of the next generation of Extremely Large Telescopes on the feasibility of detecting and characterizing the cosmological redshift drift. We consider the Lyα forest in the redshift range  2 < z < 5  and other absorption lines in the spectra of high-redshift QSOs as the most suitable targets for a redshift drift experiment. Assuming photon-noise-limited observations and using extensive Monte Carlo simulations we determine the accuracy to which the redshift drift can be measured from the Lyα forest as a function of signal-to-noise ratio and redshift. Based on this relation and using the brightness and redshift distributions of known QSOs we find that a 42-m telescope is capable of unambiguously detecting the redshift drift over a period of ∼20 yr using 4000 h of observing time. Such an experiment would provide independent evidence for the existence of dark energy without assuming spatial flatness, using any other cosmological constraints or making any other astrophysical assumption.     

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.     

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).     

Radiatively heated, protoplanetary discs with dead zones – I. Dust settling and thermal structure of discs around M stars

The irradiation of protoplanetary discs by central stars is the main heating mechanism for discs, resulting in their flared geometric structure. In a series of papers, we investigate the deep links between two-dimensional self-consistent disc structure and planetary migration in irradiated discs, focusing particularly on those around M stars. In this first paper, we analyse the thermal structure of discs that are irradiated by an M star by solving the radiative transfer equation by means of a Monte Carlo code. Our simulations of irradiated hydrostatic discs are realistic and self-consistent in that they include dust settling with multiple grain sizes  ( N = 15)  , the gravitational force of an embedded planet on the disc and the presence of a dead zone (a region with very low levels of turbulence) within it. We show that dust settling drives the temperature of the mid-plane from an   r ^−3/5  distribution (well mixed dust models) towards an   r ^−3/4  . The dead zone, meanwhile, leaves a dusty wall at its outer edge because dust settling in this region is enhanced compared to the active turbulent disc at larger disc radii. The disc heating produced by this irradiated wall provides a positive gradient region of the temperature in the dead zone in front of the wall. This is crucially important for slowing planetary migration because Lindblad torques are inversely proportional to the disc temperature. Furthermore, we show that low turbulence of the dead zone is self-consistently induced by dust settling, resulting in the Kelvin–Helmholtz instability (KHI). We show that the strength of turbulence arising from the KHI in the dead zone is  α= 10⁻⁵  .     

Radiative transfer in disc galaxies – I. A comparison of four methods to solve the transfer equation in plane-parallel geometry

Accurate photometric and kinematic modelling of disc galaxies requires the inclusion of radiative transfer models. Because of the complexity of the radiative transfer equation (RTE), sophisticated techniques are required. Various techniques have been employed for the attenuation in disc galaxies, but a quantitative comparison of them is difficult, because of the differing assumptions, approximations and accuracy requirements that are adopted in the literature. In this paper, we present an unbiased comparison of four methods to solve the RTE, in terms of accuracy, efficiency and flexibility. We apply them all to one problem that can serve as a first approximation of large portions of disc galaxies: a one-dimensional plane-parallel geometry, with both absorption and multiple scattering taken into account, with arbitrary vertical distributions of stars and dust and an arbitrary angular redistribution of the scattering. We find that the spherical harmonics method is by far the most efficient way to solve the RTE, whereas both Monte Carlo simulations and the iteration method, which are straightforward to extend to more complex geometries, have a cost that is about 170 times larger.     

Damped Lyα systems at high redshift and models of protogalactic discs

We employ observationally determined intrinsic velocity widths and column densities of damped Lyman alpha (Lyα) systems at high redshift to investigate the distribution of baryons in protogalaxies within the context of a standard cold dark matter (CDM) model. We proceed under the assumption that damped Lyα systems represent a population of cold, rotationally supported, protogalactic discs, and that the abundance of dark matter haloes is well approximated by a CDM model with critical density and vanishing cosmological constant. Using conditional cross-sections to observe a damped system with a given velocity width and column density, we compare observationally inferred velocity width and column density distributions to the corresponding theoretically determined distributions for a variety of disc parameters and CDM normalizations. In general, we find that the observations cannot be reproduced by the models for most disc parameters and CDM normalizations. Whereas the column density distribution favours small discs with large neutral gas fraction, the velocity width distribution favours large and thick discs with small neutral gas fraction. The possible resolutions of this problem in the context of this CDM model may be (1) an increased contribution of rapidly rotating discs within massive dark matter haloes to damped Lyα absorption, or (2) the abandoning of simple disc models within this CDM model for damped Lyα systems at high redshift. Here the first possibility may be achieved by supposing that damped Lyα system formation occurs only in haloes with fairly large circular velocities, and the second possibility may result from a large contribution of mergers and double discs to damped Lyα absorption at high redshift.     

The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).     

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.     

Energy transfer by the scattering of resonant photons

A formal derivation is presented of the energy transfer rate between radiation and matter due to the scattering of an isotropic distribution of resonance line photons. The derivation is developed in the context of the two-level atom in the absence of collisions and radiative transitions to and from the continuum, but includes the full angle-averaged redistribution function for photon scattering. The result is compared with previous derivations, all of which have been based on a Fokker–Planck approximation (FPA) to the radiative transfer equation. A new FPA, including a Kramers–Moyal extension to higher (post-diffusive) orders, is derived to solve the radiative transfer equation, and time-dependent numerical solutions are found. The relaxation of the colour temperature to the matter temperature is computed as the radiation field approaches statistical equilibrium through scattering. The results are discussed in the context of the Wouthuysen–Field mechanism for coupling the 21-cm spin temperature of neutral hydrogen to the kinetic temperature of the gas through Lyα scattering. The evolution of the heating rate is also computed, and shown to diminish as the gas approaches statistical equilibrium.     

A reliable cluster detection technique using photometric redshifts: introducing the 2TecX algorithm

We present a new cluster detection algorithm designed for finding high-redshift clusters using optical/infrared imaging data. The algorithm has two main characteristics. First, it utilizes each galaxy's full redshift probability function, instead of an estimate of the photometric redshift based on the peak of the probability function and an associated Gaussian error. Second, it identifies cluster candidates through cross-checking the results of two substantially different selection techniques (the name 2TecX representing the cross-check of the two techniques). These are adaptations of the Voronoi Tesselations and Friends-Of-Friends methods. Monte Carlo simulations of mock catalogues show that cross-checking the cluster candidates found by the two techniques significantly reduces the detection of spurious sources. Furthermore, we examine the selection effects and relative strengths and weaknesses of either method. The simulations also allow us to fine-tune the algorithm's parameters, and define completeness and mass limit as a function of redshift. We demonstrate that the algorithm isolates high-redshift clusters at a high level of efficiency and low contamination.     

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.     

Detecting non-sinusoidal periodicities in observational data using multiharmonic periodograms

We address the problem of assessing the statistical significance of candidate periodicities found using the so-called 'multiharmonic' periodogram, which is being used for detection of non-sinusoidal signals and is based on the least-squares fitting of truncated Fourier series. The recent author's investigation made for the Lomb–Scargle periodogram is extended to the more general multiharmonic periodogram. As a result, closed and efficient analytic approximations to the false alarm probability, associated with multiharmonic periodogram peaks, are obtained. The resulting analytic approximations are tested under various conditions using Monte Carlo simulations. The simulations showed a nice precision and robustness of these approximations.     

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_⊙  .     

The non-linear probability distribution function in models with local primordial non-Gaussianity

We use the spherical evolution approximation to investigate non-linear evolution from the non-Gaussian initial conditions characteristic of the local f _nl model. We provide an analytic formula for the non-linearly evolved probability distribution function (PDF) of the dark matter which shows that the underdense tail of the non-linear PDF in the f _nl model should differ significantly from that for Gaussian initial conditions. Measurements of the underdense tail in numerical simulations may be affected by discreteness effects, and we use a Poisson counting model to describe this effect. Once this has been accounted, our model is in good quantitative agreement with the simulations. In principle, our calculation is an important first step in programs which seek to reconstruct the shape of the initial PDF from observations of large-scale structures in the Lyα forest and the galaxy distribution at later times.     

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 GLARE Survey – II. Faint z≈ 6 Lyα line emitters in the HUDF

The galaxy population at   z ≈ 6  has been the subject of intense study in recent years, culminating in the Hubble Ultra Deep Field (HUDF) – the deepest imaging survey yet. A large number of high-redshift galaxy candidates have been identified within the HUDF, but until now analysis of their properties has been hampered by the difficulty of obtaining spectroscopic redshifts for these faint galaxies. Our 'Gemini Lyman-Alpha at Reionization Era' (GLARE) project has been designed to undertake spectroscopic follow-up of faint  ( z ' < 28.5)  i '-drop galaxies at   z ≈ 6  in the HUDF. In a previous paper we presented preliminary results from the first 7.5 h of data from GLARE. In this paper we detail the complete survey. We have now obtained 36 h of spectroscopy on a single GMOS slitmask from Gemini-South, with a spectral resolution of  λ/Δλ_FWHM≈ 1000  . We identify five strong Lyα emitters at   z > 5.5  , and a further nine possible line emitters with detections at lower significance. We also place tight constraints on the equivalent width of Lyα emission for a further ten i '-drop galaxies and examine the equivalent width distribution of this faint spectroscopic sample of   z ≈ 6  galaxies. We find that the fraction of galaxies with little or no emission is similar to that at   z ≈ 3  , but that the   z ≈ 6  population has a tail of sources with high rest-frame equivalent widths. Possible explanations for this effect include a tendency towards stronger line emission in faint sources, which may arise from extreme youth or low metallicity in the Lyman-break population at high redshift, or possibly a top-heavy initial mass function.