The Evolution QSOs — The Optical Luminosity Function

Using the data of the 2 Degree Field (2dF) QSO Redshift Survey (2QZ) and the associated 6 Degree Field (6dF) QSO Redshift Survey (6QZ) of the Anglo-Australian Telescope, a study of the Optical Luminosity Function of Quasi Stellar Objects (QSOs) has been made to understand the evolutionary scenario of QSOs. Different models for the QSO evolution are studied. The two-power law model of the optical luminosity function with second order polynomial evolution is found to fit best the observed QSO optical luminosity function. We have also determined an improved evolutionary model which fits better than the second order polynomial evolution model. The best fit parameters for the observed optical luminosity function have been determined using the Levenberg-Marquardt algorithm of non-linear least square fit for a flat universe i.e., ΩΛ + Ωm = 1 and Ho = 70 km·s·Mpc. The observed slope of the log N − m curve i.e., 1.10 ± 0.01 reveals that there were more QSOs at larger distances (or look back times) than there are locally which, in turn, indicates that the QSOs are evolving. The observed value of <V/Vmax> is found to be greater than 0.5 for different values of the cosmological constant which also hints strong evolution of the QSOs.


Introduction
The Quasi Stellar Objects (QSOs) constitute an evolving population.Their cosmological evolution was established by Schmidt (1968) [1] by applying the "V/V max " test.If the evolution of an ensemble of the QSOs can be satisfactorily described, some clues as to the evolution of the individual sources may be apparent and these may lead to a better understanding of the physics of the QSO phenomenon [2].Therefore, it is expected that the QSO evolution would throw light on the energy generation mechanism in the Active Galactic Nuclei (AGNs).Such an evolutionary scenario would also reveal information on the relationship between the QSOs and the normal galaxies.The QSOs represent an important phase of galactic evolution.Some papers e.g.[3][4][5][6] proposed joint evolutionary models of galaxies and QSOs in which mergers between gas rich galaxies drive nuclear inflow of gas fueling the growth of supermassive black holes until feedback (e.g.galactic wind) energy from black hole growth expels the surrounding gas rendering the QSO as a briefly visible bright optical source.Lapi et al. (2008) [7] proposed a model in which the QSO activity marked a transition between an earlier phase of violent and heavily dust shrouded starburst activity which promotes rapid black hole growth, and a later phase of almost passive evolution.
The luminosity function provides a basis for determining the amount of matter in different forms in the universe [8].The optical luminosity function of the QSOs and its evolution with redshift provides fundamental information on the demographics of the overall population of the AGNs.It is a useful instrument to quantify the activity statistics of the AGNs [9].The optical luminosity function of the QSOs and the Seyfert nuclei join smoothly.The shape, normalization and evolution of the luminosity function are among the most basic descriptions of the QSO population and hence of the AGNs in general [10].The optical luminosity function of QSOs as a function of redshift provides essential constraints on how the population characteristics of the QSOs have changed with time (Koo & Kron, 1988) [2] and also provides constraints on the physical models of the QSOs [11].Further, it may provide information on the models of structure formation in the early universe [12].Croom et al. 2009 [13] observed flattening of the lu-minosity function towards faint absolute magnitudes in their study of the QSO luminosity function of the 2dF-SDSS LRG and QSO (2SLAQ) survey.Clear evidence of downsizing was observed i.e. the number density of faint QSOs peaks at lower redshift than do the bright QSOs.The bright end of the QSO luminosity function tells us about the intrinsic properties of the QSO population during the time when the black holes are increasing in mass rapidly whereas the faint end of the QSO luminosity function tells us about the length of time QSOs spend at relatively low accretion rates [14].Fontanot et al. 2007 [14] computed the QSO contribution to the UV background at high redshift using the method proposed by Barger et al. 2003 [15] and confirmed that the QSO contribution to the UV background was insufficient for ionizing the intergalactic medium at these redshifts.
A brief description of the 2QZ and the 6QZ data is given in Section 2. Simple tests for the evolution of the QSOs are described in Section 3 and the study of the optical luminosity function in Section 4.

The Data
The 2-degree Field (2dF) QSO Redshift Survey (2QZ) and the associated 6-degree Field (6dF) QSO Redshift Survey (6QZ) catalogues are used for the study of QSOs and their evolution.These surveys provide more than 23,500 QSOs in a single homogeneous survey covering almost five magnitudes in the b J band (16 < b J < 20.85).
The QSO candidates for the 2QZ and 6QZ surveys are selected based on the broadband ub J r colours for the automated plate measurement (APM) of UK Schmidt Telescope (UKST) photographic plates.Details are given in [16].The survey region comprises of 30 UKST fields arranged in two 75˚ × 5˚ declination strips, one passing across the South Galactic Cap (SGP) centred on δ = -30˚ (SGP) strip and the other across the Northern Galactic Cap (NGP) centred on δ = 0˚ (referred to as the equatorial strip, also known as the NGP strip) The SGP strip extends from α = 21 h 40 to α = 3 h 15 and the equatorial strip from α = 9 h 50 to α = 14 h 50 (in B1950 coordinate system).The total survey area is 721.6 deg 2 , when allowance is made for regions of sky exercised around bright stars.
The QSO candidates were selected based on fulfilling at least one of the following colour criteria: u -b J ≤ -0.36; u -b J ≤ 0.12 -0.8(b J -r); b J -r < 0.05.The u -b J limit was tightened to u -b J ≤ -0.50 for the 6QZ sample (b J ≤ 18.25) so as to reduce the high fraction of contamination by galactic stars (for details see [17]).

The Log N -m Test
We can test whether the space density of a population of a particular class of objects is constant or not by plotting the logarithm of the cumulative distribution, N as a function of the apparent magnitude, m and measuring the slope.A slope dlogN(m)/dm > 0.6 indicates that the space density increases with distance [18,19].The plot of log N versus m plot for the 2QZ and 6QZ QSOs is shown in Figure 1.The apparent magnitude, m is taken in the b J band.The slope of the most luminous QSOs is observed as 1.10 ± 0.01 which indicates that there were more QSOs at larger distances ( or look back times) than there are locally, i.e. there were more QSOs in the past, indicating a strong evolution of the QSOs.

The Luminosity-Volume Test (V/V max Test)
The luminosity-volume test was devised by M. Schmidt [1] to test for the evolution of a population of objects.The test, also known as the V/V max test was first used by Schmidt to study the space distribution of a complete sample of radio quasars from the 3CR catalogue.A uniform distribution of objects is characterized by the <V/V max > value of 0.5 [18,19].For the full sample of 2QZ and 6QZ QSOs the <V/V max > is 0.716 ± 0.002 for Ω Λ = 0.7 universe.For the QSOs with 0.3 ≤ z ≤ 2.3 the <V/V max > value is 0.671 ± 0.002 for Ω Λ = 0.7.For Ω Λ = 0.3, 0.5 and 0.9 the values of <V/V max > are 0.648 ± 0.002, 0.674 ± 0.002 and 0.712 ± 0.002 respectively.These values of <V/V max > signify a strong evolution of the QSOs.For Ω Λ = 0.7 universe, the binned luminosity function is shown in Figure 2. We have used the Hubble constant, H o = 70 km•s -1 •Mpc -1 .The downturn at the faint end of the luminosity function is due to incompleteness of the corresponding bins [21].These incomplete bins are not included in our study of the QSO luminosity function.Similarly, for Ω Λ = 0.0, 0.3, 0.5 and 0.9 the incomplete bins have been ignored.We have assumed a flat Universe, i.e.Ω Λ + Ω m = 1.The routine mrqmin given in [22] based on the Levenberg-Marquardt method and the associated routines are used to fit the observed optical luminosity function to various models.Since the bright end of the luminosity function is steeper than the faint end despite the absence of a sharp break, a two-power-law function in the luminosity,

The QSO Optical Luminosity Function
In terms of magnitude, The evolution of the luminosity function is given by the redshift dependence of the characteristic or break luminosity, J J or break magnitude, . Different models of the luminosity evolution are studied to determine the best fit parameters for different cosmologies Ω Λ = 0.0, 0.3, 0.5, 0.7 and 0.9.These models include 1)  , the power law evolution with redshift, z; 2) , or equivalently, , the second-order polynomial evolution model; 3) , the exponential evolution with look-back time, τ.These three models have already been used in several well-known papers e.g.[9,17,23,24].
To find an alternative model the exponent "2" of z in the second order polynomial evolution model is replaced by a parameter which is varied to give a value that minimizes the χ 2 .It is observed that the minimum of χ 2 is obtained when the parameter is 1.47.Thus, another luminosity evolution model given by is determined.The luminosity functions and the predictions of the best fitting luminosity evolution model and the best fit second order polynomial evolution model are shown in Figures 3 and 4 respectively.The best-fitting parameters, the χ 2 values, the degrees of freedom ν, P(χ 2 ) and the cosmological constant Ω Λ for different evolution models are given in Table 1.  .47  ( The χ 2 values for the power law evolution with redshift are exceedingly high and the P(χ 2 ) values are very small.And for the exponential evolution model convergence of the χ 2 is not achieved.However, for the sake of completeness the estimated parameters are listed.
From Table 1, it is observed that luminosity evolution of the QSO optical luminosity function is best described by the new model with P(χ 2 ) And the second order polynomial evolution model is also acceptable with P(χ 2 ) = 7.56 × 10 -5 .The value of k 1 is larger than the values reported in the earlier papers e.g.[9,13,17,23,24].The best fit bright end slope of the luminosity function is steeper as compared to the earlier results.The inclusion of QSOs with lower spectroscopic completeness, a large fraction of which goes to the fainter magnitudes might have contributed to the steepening of the luminosity function.However, the slope of the faint end of the luminosity function agrees well with the earlier results.The values of the constant k for the model of exponential evolution with look-back time are found to be exceedingly high compared to the values reported in the above papers.The failure of the convergence of the χ 2 to less than or equal to 0.01, which The number of the observed, N obs and predicted, N pred QSOs for the best-fitting luminosity evolution model k 1 z + k 2 z 1.47 are given in Table 2.

Conclusion
A study of the evolutionary scenario of the 2QZ and the associated 6QZ QSOs is carried out.A modest attempt is made to parameterize the QSO evolution using the optical luminosity function.An improved evolutionary model in the framework of the two-power-law model of the optical luminosity function which fits the observed QSO optical luminosity function, is determined.The best fit parameters for the QSO optical luminosity function are also determined.A complete understanding of the optical luminosity function would reveal information on the AGN phenomena, the relation between the normal galaxies and the AGNs and also about the structure formation in the early universe.Of the five values of the cosmological constant, Ω Λ , the χ 2 value is found to be minimum for Ω Λ = 0.7 which agrees with the most widely accepted value.This shows that the study of the optical luminosity function might provide a tool for estimation of the value of the cosmological constant.It is hoped that information contained in the QSO population in the form of the optical luminosity function would reveal the nature of the QSO energy source and help to constrain the values of some of the cosmological parameters.

Acknowledgements
The 2dF QSO Redshift Survey (2QZ) was compiled by the 2QZ survey team from observations made with the 2degree Field on the Anglo-Australian Telescope.The au-

AFigure 1 .
Figure 1.Log N -m plot for the 2QZ and 6QZ QSOs.Here N is the cumulative number distribution of the QSOs and m, the apparent magnitude taken in the b J band.The line shown is best fit line for the luminous QSOs.

Figure 2 .
Figure 2. Binned estimate of the optical luminosity function,   J b M  as a function of the absolute magnitude,

Figure 3 .
Figure 3. Binned estimate of the optical luminosity function,   J b M  as a function of the absolute magnitude, .5 for a flat universe i.e.Ω Λ + Ω m = 1 and H o = 70 km•s -1 •Mpc -1 .The cosmological constant, Ω Λ , the evolution models, the χ 2 , degrees of freedom ν and P(χ 2 ) are also listed.
J b