• We compare predictions of large-scale cosmological hydrodynamical simulations for neutral hydrogen absorption signatures in the vicinity of 1e11 - 1e12.5 MSun haloes with observational measurements. Two different hydrodynamical techniques and a variety of prescriptions for gas removal in high density regions are examined. Star formation and wind feedback play only secondary roles in the HI absorption signatures outside the virial radius, but play important roles within. Accordingly, we identify three distinct gaseous regions around a halo: the virialized region, the mesogalactic medium outside the virial radius arising from the extended haloes of galaxies out to about two turnaround radii, and the intergalactic medium beyond. Predictions for the amount of absorption from the mesogalactic and intergalactic media are robust across different methodologies, and the predictions agree with the amount of absorption observed around star-forming galaxies and QSO host galaxies. Recovering the measured amount of absorption within the virialized region, however, requires either a higher dynamic range in the simulations, additional physics, or both.
  • The proximity profile in the spectra of z~3 quasars, where fluxes extend blueward of the He II Lya wavelength 304 (1+z) A, is one of the most important spectral features in the study of the intergalactic medium. Based on the HST spectra of 24 He II quasars, we find that the majority of them display a proximity profile, corresponding to an ionization radius as large as 20 Mpc in the source's rest frame. In comparison with those in the H i spectra of the quasars at z~6, the He II proximity effect is more prominent and is observed over a considerably longer period of reionization. The He II proximity zone sizes decrease at higher redshifts, particularly at z > 3.3. This trend is similar to that for H I, signaling an onset of He II reionization at z~4. For quasar SDSS1253+6817 (z=3.48), the He II absorption trough displays a gradual decline and serves a good case for modeling the He II reionization. To model such a broad profile requires a quasar radiation field whose distribution between 4 and 1 Rydberg is considerably harder than normally assumed. The UV continuum of this quasar is indeed exceptionally steep, and the He II ionization level in the quasar vicinity is higher than the average level in the intergalactic medium. These results are evidence that a very hard EUV continuum from this quasar produces a large ionized zone around it. Distinct exceptions are the two brightest He II quasars at z~2.8, for which no significant proximity profile is present, possibly implying that they are young.
  • We examine the impact of helium reionization on the structure of the intergalactic medium (IGM). We model the reionization using a radiative transfer (RT) code coupled to the combined gravity hydrodynamics code Enzo. Neutral hydrogen and helium are initially ionized by a starburst spectrum, which is allowed to gradually evolve into a power law spectrum over the redshift interval 3 < z < 4. The temperature-density relation of the gas is found to fan out and flatten following HeII reionization, with an inversion for overdensities above 5. Peculiar velocities of up to 10 km/s are induced by the increased pressure, with the gas density field distorted over large coherent regions by 10-20%, and the dark matter by levels of 1%. The photoionization-induced flows may thus distort the matter power spectrum at comoving wavenumbers k > 0.5 h/Mpc by a few percent by z = 2. Absorption spectra for HI and HeII are drawn from the simulations, and absorption lines are fit to the spectra. A median Doppler parameter of 35 km/s is obtained for the HI absorption systems at z = 3. Dividing into subsamples optically thick and optically thin at line centre reveals that the optically thick systems undergo only mild evolution while the optically thin systems evolve rapidly following HeII reionization. A comparison between HeII and HI absorption features shows a broad distribution in the HeII and HI column density ratio, peaking near the measured value and only slightly narrower than measured. A comparison with approximate simulation methods shows moderately good agreement in the absorption line properties, but not to the precision to which they may be measured.
  • We perform population synthesis modelling of a magnitude-limited sample of 4391 Luminous Red Galaxies selected from the Sloan Digital Sky Survey Data Release 4 (SDSS DR4). We fit measured spectral indices using a large library of high resolution spectra, covering a wide range of metallicities and assuming an exponentially decaying star-formation rate punctuated by bursts, to obtain median-likelihood estimates for the light-weighted age, metallicity, stellar mass and extinction for the galaxies. The ages lie in the range 4-10 Gyr, peaking near 6 Gyr, with metallicities in the range -0.4<[Z/H]<0.4, peaking at [Z/H] ~ 0.2. Only a few per cent of the spectra are better fit allowing for a burst in addition to continuous star-formation. The total stellar masses of all the galaxies are confined to a very narrow range. Our results broadly agree with those of previous groups using an independent population synthesis code. We find, however, that our choice in priors results in ages 1-2 Gyr smaller, decreasing the peak star formation epoch from z=2.3 to z=1.3. We develop a metal evolution model incorporating stochastic star-formation quenching motivated by recent attempts to solve the `anti-hierarchical' formation problem of ellipticals. Two scenarios emerge, a closed box with an effective stellar yield of 0.26, and an accreting box with an effective stellar yield of 0.10. Both scenarios require an IMF weighted towards massive stars and characteristic star-formation quenching times of about 100 Myr, the expected lifetime of luminous QSOs. The models predict an anti-correlation between the age and mean metallicity similar to that observed.
  • We present a spectroscopic survey of almost 15,000 candidate intermediate-redshift Luminous Red Galaxies (LRGs) brighter than i=19.8, observed with 2dF on the Anglo-Australian Telescope. The targets were selected photometrically from the Sloan Digital Sky Survey (SDSS) and lie along two narrow equatorial strips covering 180 sq deg. Reliable redshifts were obtained for 92% of the targets and the selection is very efficient: over 90% have redshifts between 0.45 and 0.8. More than 80% of the ~11,000 red galaxies have pure absorption-line spectra consistent with a passively-evolving old stellar population. The redshift, photometric and spatial distributions of the LRGs are described. The 2SLAQ data will be released publicly from mid-2006, providing a powerful resource for observational cosmology and the study of galaxy evolution.
  • We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.
  • We have used the 2dF instrument on the AAT to obtain redshifts of a sample of z<3, 18.0<g<21.85 quasars selected from SDSS imaging. These data are part of a larger joint programme: the 2dF-SDSS LRG and QSO Survey (2SLAQ). We describe the quasar selection algorithm and present the resulting luminosity function of 5645 quasars in 105.7 deg^2. The bright end number counts and luminosity function agree well with determinations from the 2dF QSO Redshift Survey (2QZ) data to g\sim20.2. However, at the faint end the 2SLAQ number counts and luminosity function are steeper than the final 2QZ results from Croom et al. (2004), but are consistent with the preliminary 2QZ results from Boyle et al. (2000). Using the functional form adopted for the 2QZ analysis, we find a faint end slope of beta=-1.78+/-0.03 if we allow all of the parameters to vary and beta=-1.45+/-0.03 if we allow only the faint end slope and normalization to vary. Our maximum likelihood fit to the data yields 32% more quasars than the final 2QZ parameterization, but is not inconsistent with other g>21 deep surveys. The 2SLAQ data exhibit no well defined ``break'' but do clearly flatten with increasing magnitude. The shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that derived from type I quasars found in hard X-ray surveys. [Abridged]
  • We present the large-scale correlation function measured from a spectroscopic sample of 46,748 luminous red galaxies from the Sloan Digital Sky Survey. The survey region covers 0.72 h^{-3} Gpc^3 over 3816 square degrees and 0.16<z<0.47, making it the best sample yet for the study of large-scale structure. We find a well-detected peak in the correlation function at 100h^{-1} Mpc separation that is an excellent match to the predicted shape and location of the imprint of the recombination-epoch acoustic oscillations on the low-redshift clustering of matter. This detection demonstrates the linear growth of structure by gravitational instability between z=1000 and the present and confirms a firm prediction of the standard cosmological theory. The acoustic peak provides a standard ruler by which we can measure the ratio of the distances to z=0.35 and z=1089 to 4% fractional accuracy and the absolute distance to z=0.35 to 5% accuracy. From the overall shape of the correlation function, we measure the matter density Omega_mh^2 to 8% and find agreement with the value from cosmic microwave background (CMB) anisotropies. Independent of the constraints provided by the CMB acoustic scale, we find Omega_m = 0.273 +- 0.025 + 0.123 (1+w_0) + 0.137 Omega_K. Including the CMB acoustic scale, we find that the spatial curvature is Omega_K=-0.010+-0.009 if the dark energy is a cosmological constant. More generally, our results provide a measurement of cosmological distance, and hence an argument for dark energy, based on a geometric method with the same simple physics as the microwave background anisotropies. The standard cosmological model convincingly passes these new and robust tests of its fundamental properties.
  • Constraints on the ionization structure of the Intergalactic Medium are derived as directly imposed by observations in conjunction with the results of numerical simulations for structure formation. Under the assumption that the population of sources dominating the UV ionizing background at z<6 is the same population which reionized the IGM, it is shown that consistency with measurements of the mean LyA transmitted flux at high redshifts requires the epoch of hydrogen reionization of the IGM to have occurred at a redshift z < 11, independent of the space density of the sources. The upper limit on the reionization redshift depends only on the shape of the UV spectra of the sources. Consistency with constraints on the reionization epoch from WMAP requires the sources of photoionization to have had hard spectra, such as QSOs or Population III stars. The only way to escape these conclusions is either: 1.the sources which dominate the photoionization background at z<6 are the remnants of a population that was a much more prodigious source of ionizing photons at earlier times, or 2.the sources responsible for the photoionization of the IGM are an unknown population that contributes negligibly to the ionization of the IGM at z<6. The contribution of QSOs to the UV background is also estimated based on recent surveys.
  • We model how repeated supernova explosions in high-redshift dwarf starburst galaxies drive superbubbles and winds out of the galaxies. We compute the efficiencies of metal and mass ejection and energy transport from the galactic potentials, including the effect of cosmological infall of external gas. The starburst bubbles quickly blow out of small, high-redshift, galactic disks, but must compete with the ram pressure of the infalling gas to escape into intergalactic space. We show that the assumed efficiency of the star formation rate dominates the bubble evolution and the metal, mass, and energy feedback efficiencies. With star formation efficiency f*=0.01, the ram pressure of infall can confine the bubbles around high-redshift dwarf galaxies with circular velocities v_c>52 km/s. We can expect high metal and mass ejection efficiencies, and moderate energy transport efficiencies in halos with v_c~30-50 km/s and f*~0.01 as well as in halos with v_c~100 km/s and f*>>0.01. Such haloes collapse successively from 1-2 sigma peaks in LambdaCDM Gaussian density perturbations as time progresses. These dwarf galaxies can probably enrich low and high-density regions of intergalactic space with metals to 10^-3-10^-2 Zsun as they collapse at z~8 and z<5 respectively. They also may be able to provide adequate turbulent energy to prevent the collapse of other nearby halos, as well as to significantly broaden Lyman-alpha absorption lines to v_rms~20-40 km/s. We compute the timescales for the next starbursts if gas freely falls back after a starburst, and find that, for star formation efficiencies as low as f*<0.01, the next starburst should occur in less than half the Hubble time at the collapse redshift. This suggests that episodic star formation may be ubiquitous in dwarf galaxies.
  • Numerical hydrodynamical simulations are used to predict the expected absorption properties of the Lyman-Alpha forest for a variety of Cold Dark Matter dominated cosmological scenarios: CHDM, OCDM, LCDM, SCDM, and tCDM. Synthetic spectra are constructed duplicating the resolution, signal-to-noise ratio, and wavelength coverage of several published high resolution spectra, and their statistical properties compared on the basis of the flux distribution of the spectra, the distribution of coefficients in a wavelet decomposition of the spectra, and the distributions of absorption line profile parameters. Agreement between the measured and predicted cumulative distributions is found at the few to several percent level. The best-fitting models to the flux distribution correspond to normalizations on the scale of the cosmological Jeans length of 1.3 < sigma_J < 1.7 at z=3. No single model provides an acceptable match to all the distributions. Larger median Doppler parameters are found in the measured spectra than predicted by all but the lowest normalization models (CHDM and tCDM), which provide poor fits to the flux distributions. The discrepancy in Doppler parameters is particularly large for absorption systems optically thin at the Lyman-Alpha line-centre. This may indicate the need to introduce additional energy injection throughout the IGM, as may be provided by late HeII reionization (z_HeII~3.5) or supernovae-driven winds from young galaxies. The models require a hydrogen ionization rate at redshifts 1.7<z<3.5 within a factor of 2 of that predicted from QSOs alone as the sources of the UV photoionization background, although with a slower rate of decline with redshift at z>3.5 than predicted from current QSO counts.
  • Numerical hydrodynamical simulations have proven a successful means of reproducing many of the statistical properties of the Lyman-Alpha forest as measured in high redshift quasar spectra. Pseudo-hydrodynamical methods based only on simulating the dark matter component have been claimed to yield a comparable level of success. We investigate the degree to which two pseudo-methods, with and without allowing for a pseudo-gas pressure, are able to match the predictions of fully hydrodynamical plus dark matter simulations. We also address the requirements for convergence to the statistics of the spectra and the inferred properties of the Lyman-Alpha forest as a function of resolution and box size. Generally we find it is possible to reach agreement with full hydrodynamic simulations at the 10% level in the cumulative distributions of the flux and absorption line parameter statistics for readily achievable particle and grid numbers, but difficult to do much better.
  • It is shown how wavelets may be used to analyse the absorption properties of the Lya forest. The Discrete Wavelet Transform of a QSO spectrum is used to decompose the light fluctuations that comprise the forest into orthogonal wavelets. It is demonstrated that most of the signal is carried by the moderate to lower frequency wavelets in high resolution spectra, and that a statistically acceptable description of even high signal-to-noise spectra is provided by only a fraction (10-30%) of the wavelets. The distributions of the wavelet coefficients provide a statistical basis for discriminating between different models of the Lya forest. The method is illustrated using the measured spectrum of Q1937-1009. The procedure described is readily automated and may be used to process both measured spectra and the large number of spectra generated by numerical simulations, permitting a fair comparison between the two.
  • We find statistically significant correlations in the cosmological matter power spectrum over the full range of observable scales. While the correlations between individual modes are weak, the band-averaged power spectrum shows strong non-trivial correlations. The correlations are significant when the modes in either one or both bands are in the non-linear regime, and approach 100% for pairs of bands in which all the modes are non-linear. The correlations are weaker, but not absent, when computed in redshift space. Since estimates of the power spectrum from galaxy surveys require band-averaging, the correlations must be taken into account when comparing a measured power spectrum with theoretical models.
  • The comparison of quasar absorption spectra with numerically simulated spectra from hierarchical cosmological models of structure formation promises to be a valuable tool to discriminate among these models. We present simulation results for the column density, Doppler b parameter, and optical depth probability distributions for five popular cosmological models.
  • We investigate the consequences of a non-negligible baryon fraction for models of structure formation in Cold Dark Matter dominated cosmologies, emphasizing in particular the existence of oscillations in the present-day matter power spectrum. These oscillations are the remnants of acoustic oscillations in the photon-baryon fluid before last scattering. For acceptable values of the cosmological and baryon densities, the oscillations modulate the power by up to 10%, with a `period' in spatial wavenumber which is close to Delta k approximately 0.05/ Mpc. We study the effects of nonlinear evolution on these features, and show that they are erased for k > 0.2 h/ Mpc. At larger scales, the features evolve as expected from second-order perturbation theory: the visibility of the oscillations is affected only weakly by nonlinear evolution. No realistic CDM parameter combination is able to account for the claimed feature near k = 0.1 h/ Mpc in the APM power spectrum, or the excess power at 100 Mpc/h wavelengths quoted by several recent surveys. Thus baryonic oscillations are not predicted to dominate existing measurements of clustering. We examine several effects which may mask the features which are predicted, and conclude that future galaxy surveys may be able to detect the oscillatory features in the power spectrum provided baryons comprise more than 15% of the total density, but that it will be a technically challenging achievement.