A simple model for the evolution of supermassive black holes and the quasar population

An empirically motivated model is presented for accretion-dominated growth of supermassive black holes (SMBH) in galaxies, and the implications are studied for the evolution of the quasar population in the Universe. We investigate the core aspects of the quasar population, including space density evolution, evolution of the characteristic luminosity, plausible minimum masses of quasars, the mass function of SMBH and their formation epoch distribution. Our model suggests that the characteristic luminosity in the quasar luminosity function arises primarily as a consequence of a characteristic mass scale above which there is a systematic separation between the black hole and the halo merging rates. At lower mass scales, black hole merging closely tracks the merging of dark haloes. When combined with a declining efficiency of black hole formation with redshift, the model can reproduce the quasar luminosity function over a wide range of redshifts. The observed space density evolution of quasars is well described by formation rates of SMBH above  ∼10⁸  M_⊙  . The inferred mass density of SMBH agrees with that found independently from estimates of the SMBH mass function derived empirically from the quasar luminosity function.     

Evidence for evolution in the clustering of quasars

We discuss the changes in the amplitude of the correlation function of quasar clustering with the redshift. Using the new Hewitt-Burbidge QSO catalogue, we find from a statistical analysis of the quasar pairs that the positive results of quasar clustering are mainly contributed by low-redshift quasars. There is an obvious change in the quasar distribution at z = 2. There is no clustering for z > 2. The difference in clustering is not due to any selection effect on the brightness. The evidence is that quasar clustering evolved with cosmological time.     

The evolution of star formation in quasar host galaxies

We have used far-infrared data from IRAS , Infrared Space Observatory ( ISO ), Spitzer Wide-Area Infrared Extragalactic (SWIRE), Submillimetre Common User Bolometer Array (SCUBA) and Max-Planck Millimetre Bolometer (MAMBO) to constrain statistically the mean far-infrared luminosities of quasars. Our quasar compilation at redshifts  0 < z < 6.5  and I -band luminosities  −20 < I _AB < −32  is the first to distinguish evolution from quasar luminosity dependence in such a study. We carefully cross-calibrate IRAS against Spitzer and ISO , finding evidence that IRAS 100-μm fluxes at <1 Jy are overestimated by ∼30 per cent. We find evidence for a correlation between star formation in quasar hosts and the quasar optical luminosities, varying as star formation rate (SFR)  ∝ L ^0.44±0.07_opt  at any fixed redshift below   z = 2  . We also find evidence for evolution of the mean SFR in quasar host galaxies, scaling as  (1 + z )^1.6±0.3  at   z < 2  for any fixed quasar I -band absolute magnitude fainter than −28. We find no evidence for any correlation between SFR and black hole mass at  0.5 < z < 4  . Our data are consistent with feedback from black hole accretion regulating stellar mass assembly at all redshifts.     

The low-redshift evolution of the Lyman-β forest

The low-redshift evolution of the intergalactic medium is investigated using hydrodynamic cosmological simulations. The assumed cosmological model is a critical density cold dark matter universe. The imposed uniform background of ionizing radiation has the amplitude, shape and redshift evolution as computed from the observed quasar luminosity function by Haardt &38; Madau. We have analysed simulated Lyman-α spectra using Voigt-profile fitting, mimicking the procedure with which quasar spectra are analysed. Our simulations reproduce the observed evolution of the number of Lyman-α absorption lines over the whole observed interval of z  = 0.5 to 4. In particular, our simulations show that the decrease in the rate of evolution of Lyman-α absorption lines at z  ≤ 2, as observed by the Hubble Space Telescope , can be explained by the steep decline in the photoionizing background resulting from the rapid decline in quasar numbers at low redshift.     

The time evolution of cosmological redshift as a test of dark energy

The variation of the expansion rate of the Universe with time produces an evolution in the cosmological redshift of distant sources (e.g. quasar Lyman α absorption lines) that might be directly observed by future ultrastable, high-resolution spectrographs (such as the COsmic Dynamics Experiment) coupled to extremely large telescopes (such as the European Southern Observatory's Extremely Large Telescope). This would open a new window to explore the physical mechanism responsible for the current acceleration of the Universe. We investigate the evolution of cosmological redshift from a variety of dark energy models, and compare it with simulated data. We perform a Fisher matrix analysis and discuss the prospects for constraining the parameters of these models and for discriminating among competing candidates. We find that, because of parameter degeneracies, and the inherent technical difficulties involved in this kind of observations, the uncertainties on parameter reconstruction can be rather large unless strong external priors are assumed. However, the method could be a valuable complementary cosmological tool, and give important insights on the dynamics of dark energy, not obtainable using other probes.     

Periodic redshift distribution of quasars associated with low — redshift galaxies

We applied power spectral analysis to the redshift distribution of the quasars associated with low-redshift galaxies. Periodicity is confirmed at a level of confidence >99%, with a period of Δln(1+z) = 0.206. This result is exactly the same as previously found for different quasar samples. The periodicity is not caused by any selection effect in the emission lines. The two phenomena, the periodicity in the redhsift distribution and the association of quasars with low-redhsift galaxies can both be rather well explained by a multiply-connected model of the universe.     

On the time evolution of quasar luminosity and redshift distribution

The possibility to predict the redshift distribution on the basis of certain assumptions about the time-dependence of quasar luminosity is discussed. It is found that for a certain general class of models that resembles the spinar model for quasars only specific combinations of parameters give rise to acceptable results. Further consequences are discussed.     

Low correlation between redshift and intrinsic brightness for quasars and galaxies in a de sitter universe

Since the discovery of quasars, it has been recognized that these objects must have either an extraordinary intrinsic brightness or a nonlinear redshift. The most widely accepted current belief incorporates a linear (Hubble) redshift- distance relation and time evolution within a Big Bang model. We reconsider the possibility of a nonlinear (de Sitter) redshift-distance relation and find quasar intrinsic brightness to be not at all extraordinary. Given a de Sitter law, intrinsic brightness is found to be independent of redshift over five orders of magnitude.     

Quasar evolution

I present a critical review of recent progress in the field of quasar evolution, concentrating on the methodical issues of constructing and analysing observed quasar samples. In particular, I attempt to trace the origins and shortcomings of the popular concept of ‘pure luminosity evolution’ that has dominated the discussion over the last decade. I discuss results from new surveys showing that this notion should be abandoned, and that quasar evolution is more complex than previously thought.     

Using galaxy redshift surveys to detect gravitationally lensed quasars

Gravitationally lensed quasars can be discovered as a by-product of galaxy redshift surveys. Lenses discovered spectroscopically in this way should require less observational effort per event than those found in dedicated lens surveys. Further, the lens galaxies should be relatively nearby, facilitating a number of detailed observations that are impossible for the more common high-redshift lenses. This is epitomized by the wide range of results that have been obtained from Q 2237+0305, which was discovered as part of the Center for Astrophysics redshift survey, and remains the only quasar lens discovered in this way. The likelihood of this survey yielding a lens is calculated to be ∼0.03, which is an order of magnitude larger than previous estimates due to two effects. First, the quasar images themselves increase the observed flux of the lens, so that lens galaxies up to a magnitude fainter than the nominal survey limit must be included in the calculation. Secondly, it is possible for lensed quasars with extremely faint deflectors to enter the survey due to the extended morphology of the multiple images. Extrapolating these results to future surveys, the 2 degree Field galaxy redshift survey should contain between 10 and 50 lenses and the Sloan Digital Sky Survey should yield between 50 and 300 lenses, depending on the cosmological model and the observing conditions.     

Evolution of the u-band luminosity function from redshift 1.2 to 0

We produce and analyse u -band (  λ≈ 355  nm) luminosity functions (LFs) for the red and blue populations of galaxies using data from the Sloan Digital Sky Survey (SDSS) u -band Galaxy Survey ( u GS) and Deep Evolutionary Exploratory Probe 2 (DEEP2) survey. From a spectroscopic sample of 41 575 SDSS u GS galaxies and 24 561 DEEP2 galaxies, we produce colour magnitude diagrams and make use of the colour bimodality of galaxies to separate red and blue populations. LFs for eight redshift slices in the range  0.01 < z < 1.2  are determined using the  1/ V _max  method and fitted with Schechter functions showing that there is significant evolution in   M *  , with a brightening of 1.4 mag for the combined population. The integration of the Schechter functions yields the evolution in the u -band luminosity density (LD) out to   z ∼ 1  . By parametrizing the evolution as  ρ∝ (1 + z )^β  , we find that  β= 1.36 ± 0.2  for the combined populations and  β= 2.09 ± 0.2  for the blue population. By removing the contribution of the old stellar population to the u -band LD and correcting for dust attenuation, we estimate the evolution in the star formation rate (SFR) of the Universe to be  β_SFR= 2.5 ± 0.3  . Discrepancies between our result and higher evolution rates measured using the infrared and far-UV can be reconciled by considering possibilities such as an underestimated dust correction at high redshifts or evolution in the stellar initial mass function.     

Probing the evolution of early-type galaxies using multicolour number counts and redshift distributions

We investigate pure luminosity evolution models for early-type (elliptical and S0) galaxies (i.e. no number density change or morphological transition), and examine whether these models are consistent with observed number counts in the B , I and K bands, and redshift distributions of two samples of faint galaxies selected in the I and K bands. The models are characterized by the star formation time-scale τ _SF and the time t _gw when the galactic wind starts to blow, in addition to several other conventional parameters. We find that the single-burst model ( τ _SF=0.1 Gyr and t _gw=0.353 Gyr), which is known to reproduce the photometric properties of early-type galaxies in clusters, is inconsistent with the redshift distributions of early-type galaxies in the field environment, owing to overpredictions of the number of galaxies at z ≳1.4 even with strong extinction which is at work until t _gw. In order for dust extinction to be more effective, we treat τ _SF and t _gw as free parameters, and find that models with τ _SF≳0.5 Gyr and t _gw>1.0 Gyr can be made consistent with both the observed redshift distributions and the number counts, if we introduce strong extinction [ E ( B − V )≥1 as a peak value]. These results suggest that early-type galaxies in the field environment do not have the same evolutionary history as described by the single-burst model.     

Chemo-spectrophotometric evolution of spiral galaxies – V. Properties of galactic discs at high redshift

We explore the implications for the high-redshift universe of 'state-of-the-art' models for the chemical and spectrophotometric evolution of spiral galaxies. The models are based on simple 'scaling relations' for discs, obtained in the framework of cold dark matter models for galaxy formation, and were 'calibrated' so as to reproduce the properties of the Milky Way and of nearby discs (at redshift z ∼0) . In this paper, we compare the predictions of our 'hybrid' approach to galaxy evolution to observations at moderate and high redshift. We find that the models are in fairly good agreement with observations up to z ∼1 , while some problems appear at higher redshift (provided there is no selection bias in the data); these discrepancies may suggest that galaxy mergers (not considered in this work) played a non-negligible role at z >1 . We also predict the existence of a 'universal' correlation between abundance gradients and disc scalelengths, independent of redshift.     

The UVX quasar optical luminosity function and its evolution

The recently finished Edinburgh UVX quasar survey at B <18 is used together with other complete samples to estimate the shape and evolution of the optical luminosity function in the redshift range 0.3< z <2.2. There is a significantly higher space density of quasars at high luminosity and low redshift than previously found in the PG sample of Schmidt & Green, with the result that the shape of the luminosity function at low redshifts ( z <1) is seen to be consistent with a single power law. At higher redshifts the slope of the power law at high luminosities appears to steepen significantly. There does not appear to be any consistent break feature which could be used as a tracer of luminosity evolution in the population.     

The detection and X-ray evolution of galaxy groups at high redshift

We describe some of the first X-ray detections of groups of galaxies at high redshifts  ( z ∼0.4)  , based on the UK deep X-ray survey of M^cHardy et al. Combined with other deep ROSAT X-ray surveys with nearly complete optical identifications, we investigate the X-ray evolution of these systems. We find no evidence for evolution of the X-ray luminosity function up to   z =0.5  at the low luminosities of groups of galaxies and poor clusters  ( L _X≳10^42.5 erg s^-1)  , although the small sample size precludes very accurate measurements. This result confirms and extends to lower luminosities current results based on surveys at brighter X-ray fluxes. The evolution of the X-ray luminosity function of these low-luminosity systems is more sensitive to the thermal history of the intragroup medium (IGM) than to cosmological parameters. Energy injection into the IGM (from, for example, supernovae or active galactic nuclei winds) is required to explain the X-ray properties of nearby groups. The observed lack of evolution suggests that the energy injection occurred at redshifts   z >0.5  .     

Quasar correlation function and redshift space distortion from SDSS DR5 data

For z = 0.8–2.2 redshift interval, quasar pair correlation function parameters and β redshift space distortion parameter (connected to large-scale potential flows) values are estimated. We base them on the Main QSO Sample from SDSS Data Release 5. Standard correlation function form ξ(r) = (r ₀/r)^γ is used for comoving distances r = 2–50 Mpc between quasars. We fix the parameters of the cosmological model: Ω_Λ = 1 − Ω _M = 0.726 and H ₀ = 70.5 km/(s Mpc). We come to the best-fit parameter values of γ = 1.77 ± 0.20, r ₀ = 5.52 ± 0.95 Mpc/h for r in the range 2–30 Mpc, γ = 1.91 ± 0.11, r ₀ = 5.82 ± 0.61 Mpc for r in the range 2–50 Mpc. The mean β value is β = 0.43 ± 0.22.     

Cosmological evolution and hierarchical galaxy formation

We calculate the rate at which dark matter haloes merge to form higher mass systems. Two complementary derivations using Press–Schechter theory are given, both of which result in the same equation for the formation rate. First, a derivation using the properties of the Brownian random walks within the framework of Press–Schechter theory is presented. We then use Bayes' theorem to obtain the same result from the standard Press–Schechter mass function. The rate obtained is shown to be in good agreement with results from Monte Carlo and N -body simulations. We illustrate the usefulness of this formula by calculating the expected cosmological evolution in the rate of star formation that is due to short-lived, merger-induced starbursts. The calculated evolution is well-matched to the observed evolution in ultraviolet luminosity density, in contrast to the lower rates of evolution that are derived from semi-analytic models that do not include a dominant contribution from starbursts. Hence we suggest that the bulk of the observed ultraviolet starlight at z >1 arises from merger-induced starbursts. Finally, we show that a simple merging-halo model can also account for the bulk of the observed evolution in the comoving quasar space density.     

The chemical evolution of the Universe – I. High column density absorbers

We construct a simple, robust model of the chemical evolution of galaxies from high to low redshift, and apply it to published observations of damped Lyman α quasar absorption line systems (DLAs). The elementary model assumes quiescent star formation and isolated galaxies (no interactions, mergers or gas flows). We consider the influence of dust and chemical gradients in the galaxies, and hence explore the selection effects in quasar surveys. We fit individual DLA systems to predict some observable properties of the absorbing galaxies, and also indicate the expected redshift behaviour of chemical element ratios involving nucleosynthetic time delays.
Despite its simplicity, our 'monolithic collapse' model gives a good account of the distribution and evolution of the metallicity and column density of DLAs, and of the evolution of the global star formation rate and gas density below redshifts z ∼3. However, from the comparison of DLA observations with our model, it is clear that star formation rates at higher redshifts ( z >3) are enhanced. Galaxy interactions and mergers, and gas flows very probably play a major role.     

Chemical evolution of the high redshift galaxies

We have tried to determine the rate of chemical evolution of high redshift galaxies from the observed redshift distribution of the heavy element absorption systems in the spectra of QSOs, taking into account the evolution in the intensity of the metagalactic UV ionizing radiation background, the radius and/or the co-moving number density of, and the fraction of mass in the form of gas in, the absorbers. The data for both the Lyman limit systems and the C IV systems have been fitted simultaneously. It seems that the abundance of carbon has possibly increased by about a factor of 5 to 20 from the cosmic time corresponding to the redshift ≃ 4 to 2. The data also suggest that either the radius or the co-moving number density of the galaxies increased with redshift up to z = 2.0 and decreased slowly thereafter. The total mass of the halo gas was higher in the past, almost equal to the entire mass of the galaxy at z = 4. The hydrogen column density distribution for Lyman limit systems predicted by the model is in agreement with the observed distribution.     

Seasonal dependences of geomagnetic variations in the polar region in connection with large-amplitude annual Z-variation at the geomagnetic pole

A noticeable depression of the vertical component Z of the geomagnetic field is observed in the polar cap in summer. From the statistical analysis of the equivalent overhead current patterns for daily geomagnetic variations in the summer and winter polar regions for three different conditions of IMF (interplanetary magnetic field), it was concluded that the annual variation of geomagnetic Z in the vicinity of the geomagnetic pole is attributed to the relative spatial shift of the twin-vortex current patterns over the polar cap from summer to winter. In winterthe clockwise current vortex in the dawn sector extends almost over the entire polar cap (except for the periods when the B_z component of IMF has a large positive value), and this will result in the positive variation of the Z-value at the geomagnetic pole. In summer the counter-clockwise current vortex in the dusk sector always extends over the whole polar cap even when B_z of IMF is positiveso that the variation of Z becomes negative. The persistent existence of current vortex in the dawn sector is important for the further study of magnetospheric convection when B_z is positive.