a fundamental tool to study the evolution of active galactic nuclei ( agns ) is their luminosity function ( lf ) and its evolution with redshift .
this enables us to probe the accretion history of supermassive black holes ( smbhs ) which , in turn , set constraints on galaxy evolution , as smbhs are believed to be ubiquitous and tightly linked with all stellar bulges and spheroids today ( e.g. ferrarese & ford 2005 and references therein ) . moreover constraining the agn lf is essential to understand the agn contribution to the total emissivity background in different bands , their environmental impact , their clustering properties and the role agns had , if any , in the cosmological re - ionization of hydrogen and helium .
x - ray observations carried out with _ chandra _ , _ asca _ , _ heao-1 _ are statistically complete in x - ray to a level of @xmath11 up to @xmath12 ( e.g. la franca _ et al .
these observations have assessed that agns statistically follow a luminosity - dependent density evolution ( e.g. ueda et al . 2003 ) , in which lower luminous agns peak at @xmath13 , while higher luminous ones , usually optically bright , peak at higher redshifts , @xmath14 ( e.g. osmer 2004 ) . however
how far in time can we track agn evolution ?
the picture is still not yet complete at higher redshifts , where x - ray measurements of the agn lf are lacking .    on the optical side
the sdss team ( richard _ et al .
_ 2006 ) , and earlier works by kennefick _
( 1995 ) , combo-17 ( wolf _ et al .
_ 2003 ) have probed the bright tail of the optical agn lf up to @xmath15 down to magnitudes of @xmath16 showing that a simple power - law can easily fit the data . the goods team ( cristiani _ et al .
_ 2004 ; fontanot _ et al . _
2006a ) has been able to set constraints on the shape of the optical lf down to @xmath17 and up to @xmath18 , showing evidence for a flattening at the faint - end .
( 2004 ; f04 from here on ) have instead probed the very luminous tail of the @xmath0 optical lf up to @xmath19 .
here we use existing observations of optically selected agn and _ chandra _ deep fields to determine constraints on the redshift 6 shape of both the optical and x - ray lf .
we will in particular show that x - ray observations allow us to probe the @xmath0 lf to @xmath20 magnitudes fainter than the deepest optical observations . in section
[ sec|method ] we describe our method , present the constraints on the faint end lf in section  [ sec|results ] , their implications in sections  [ sec|implic ] and  [ sec|reionization ] , and conclude in section  [ concl ] .
we use concordance cosmological parameters @xmath21 , @xmath22 , @xmath23 , @xmath24 km s@xmath25 mpc@xmath25 , @xmath26 , and @xmath27 ( spergel _ et al .
in the following we will always refer to the total population of active galactic nuclei as agns .
we will call `` qsos '' ( or `` quasars '' ) the type i , optical agns , and type ii agns those which have a column density above @xmath28 . the latter kind of sources are often very faint or even missed in optical surveys , but still detected in x - rays .
several deep x - ray and optical surveys have now probed several portions of the agn luminosity function at high redshifts , however often finding a scarce number of faint sources . in this paper
we will show that the paucity of sources at these faint fluxes gives the most up - to - date stringent constraints on the knee and the faint - end slope of the qso lf and give hints to the overall agn distribution at @xmath0 .
similar studies , though with looser upper limits , have been carried out by sharp _
( 2004 ) , willott _
( 2005 ) , wyithe ( 2004 ) and dijkstra & wyithe ( 2006 ) .
to convert x - ray to optical rest - frame luminosities we will use the equation ( see richards _
_ 2005 ) @xmath29 between monochromatic luminosities and a constant photon index @xmath30 to compute the broad - band x - ray luminosities .
many authors ( e.g. vignali _ et al .
_ 2003 ) have found that the spectral index @xmath31 depends significantly on luminosity . in particular , the latest results by steffen _
( 2006 ) , who have collected a large number of optically and x - ray selected agns at several redshifts , claim a significant correlation of the spectral index with luminosity and a marginal one ( 1.1@xmath32 ) with reshift , i.e. @xmath33 where @xmath34 is the monochromatic luminosity at @xmath35 , in units of @xmath36 .    in the following
we will use results from several optical / nir surveys carried out on different areas . to compare the luminosities probed by different observations
, we will first need to convert the broadband z@xmath37 , where most of the observations are carried out , to the apparent magnitude at wavelength @xmath38 , i.e. ab@xmath39 , and finally convert the latter to the rest - frame luminosities m@xmath40 .
we convert z@xmath37 to ab@xmath39 magnitude following the spectrum described in cristiani & vio ( 1990 ) and updated in cristiani _ et al . _
( 2004 ) . to convert then to absolute magnitudes m@xmath40
, we note that in the ab system @xmath41 , where the monochromatic flux @xmath42 is in units of @xmath43 , and the absolute magnitude is m@xmath44 , where l@xmath45 is the luminosity at rest frame frequency @xmath46 , in units of @xmath47 .
using @xmath48 leads to ( see also stern _
2000b ) @xmath49 and ( fan _ et al .
_ 2001 ) @xmath50 where @xmath51 is the absolute b magnitude calibrated for a @xmath52 ( where we set @xmath53 ) spectrum and the factor 0.12 comes from the zero point difference between the ab and vega - based system ( fan _ et al .
_ 2001 ) . to constrain the shape of the @xmath0 agn lf , following willott _
( 2005 ) , we first define a double power - law model lf of the type @xmath54 where @xmath55 , with @xmath56 , is the average redshift of the f04 lf .
with such a parametrization we assume that the high - redshift agn lf evolves with a pure density evolution ( pde ) beyond redshift of 5 at the same rate as the lf in fan _
( 2001 ) . following the observations by f04
, @xmath57 has been calibrated to match the cumulative number density of @xmath58 for sources with @xmath59 at @xmath60 ( these are the limit and number density for @xmath61 , @xmath62 , @xmath63 and have been converted to our cosmology ) .
the bright - end slope and the break magnitude have been fixed to @xmath64 and @xmath65 respectively , as this provides a very good match to the cool et al .
2006 observations ( see below ) .
the faint - end slope @xmath66 has been left variable and constrained by matching the observed number density of sources in each survey . at @xmath18 the lf in equation  ( [ eq|modellf ] ) is in agreement with the cumulative number density determined by fan _
( 2001 ) .    thus given a model
lf we integrate it within the effective volume of the survey using the following equation ( see e.g. fan _ et al .
_ 2001 ) @xmath67 and tune the faint - end slope @xmath66 in equation  ( [ eq|modellf ] ) to match the observed number of sources @xmath68 in the survey .
@xmath69 refers to the redshift dependent absolute luminosity , corresponding to the minimum optical apparent magnitude ( or the minimum flux ) considered .
we will in general assume a constant selection function p(l , z , sed ) for each survey , as indicated by the authors or derived from their papers . in the following we will always use the form in equation  ( [ eq|modellf ] ) for @xmath70 , unless in some specified cases , in which we will adopt the fan et al .
( 2001 ) lf . the latter is a single power - law function valid in the range @xmath71 , in which we let the slope @xmath72 to be variable , while we keep the redshift evolution and normalization at @xmath73 as in fan _
( 2001 ; both @xmath73 and normalization have been converted to our cosmology , which is required at each redshift bin during the evolution ) .
a significant improvement with respect to previous works has been to fix the knee of the lf in equation  ( [ eq|modellf ] ) to the value recently inferred by cool _
these authors discovered three objects down to z@xmath37=21.93 , ( i.e @xmath74 at z@xmath75 using equation  ( [ eq|mzmbm1450 ] ) ) in the agn and galaxy evolution survey ( ages ) spectroscopic observations of the noao deep wide - field survey ( ndwfs ) .
they covered 7.7 deg@xmath76 of the ndwfs field with an average completeness level of 54% down to @xmath77 mag .
these authors claim their observations agree with the f04 lf with a slope of @xmath78 and a redshift evolution as in fan et al .
( 2001 ) and in equation  ( [ eq|modellf ] ) .
using equation  ( [ eq|integrallf ] ) , we find @xmath79 sources with z@xmath80 at @xmath81 , in very good agreement with observations and with the authors . the integration was pursued up to @xmath82 , as beyond this limit the optical selection of quasars becomes impossible because the ly@xmath83 line redshifs out of the z@xmath37-band ( see cool _ et al .
_ 2006 ) .
following a similar procedure we find significant evidence for a flattening of the qso lf below @xmath74 .
recently , mahabal _
( 2005 , m05 hereafter ) in the field around the @xmath84 quasar sdss j1148 + 5251 discovered a quasar at redshift @xmath85 with luminosity of @xmath86 , i.e. z@xmath87 .
this quasar is the faintest known optically selected quasar at @xmath88 , allowing determination of the faint end of the quasar lf to fainter luminosities .
the authors claim the source , identified in a @xmath89 deg@xmath76 survey , has a surface density consistent with the wolf _
( 2003 ) lf extrapolated through a pde within their effective volume , which follows the @xmath90-band drop - out criterion of stern _ et al . _
( 2000 ) , i.e. @xmath91 .
we have checked that adopting the latter determination for the effective volume , we recover the results of the authors only if an average @xmath92 selection detection probability is used . however using the wolf _
( 2003 ) lf , which is defined only for @xmath93 , needs an extrapolation to higher redshifts .
moreover the functional form in equation  ( [ eq|modellf ] ) has been calibrated on the data by cool _
( 2006 ) at @xmath94 .
we therefore replace equation  ( [ eq|modellf ] ) with the fan _
( 2001 ) qso lf , more appropriate for the redshift bin probed by the @xmath90-band drop - out technique .    keeping the same value for the selection function we find that the detection of one object in the effective volume probed by m05 , is consistent with a number density at @xmath86 of @xmath95\sim -6.70 $ ] at the 99% confidence level ( cl ; with @xmath96 ) , down to @xmath95\sim -6.95 $ ] at the 90% cl ( with @xmath97 ; see table  [ table|survey ]
; we plot the m05 data in figure  [ fig|lfz6 ] shifted by -0.15 magnitudes for avoiding overcrowding with willott et al.s data ) .
their data , shown as solid triangles in figure  ( [ fig|lfz6 ] ) , are consistent with extrapolating equation  ( [ eq|modellf ] ) to fainter optical magnitudes using @xmath98 and @xmath99 at the 99% and 90% cl respectively ( see table  [ table|survey ] ) .
note that these are conservative estimates .
extending the computation of the lf down to the limiting magnitude of the m05 survey ( z@xmath100 , @xmath101 ) or using a higher completeness function would yield even lower limits on @xmath66 .
( 2005 ) undertook a multicolor quasar optical imaging with the canada - france - hawaii telescope .
they found a lack of quasars down to a limiting magnitude of z@xmath10223.35 ( i.e. @xmath103 at @xmath104 ; see figure  [ fig|lfz6 ] ) .
the few detected quasar candidates had been subsequently rejected as being identified with low - mass stars in the near - infrared follow - up imaging . using equations  ( [ eq|modellf ] ) and  ( [ eq|integrallf ] )
, we find that the willott _ et al . _
( 2005 ) results are consistent with a faint - end slope of @xmath105 at 99% cl , down to @xmath106 at the 90% cl .
we use the selection function computed from simulated colors for high - redshift qsos in the redshift range @xmath107 ( see figure 2 in willott _ et al .
this is the most conservative selection function from willott _
( 2005 ) , determined for very extreme red color cuts . using a less conservative estimate ( solid line in their figure 2 ) our results drop to @xmath108 at 99% cl , down to @xmath2 at the 90% cl .
the above findings are shown in figure  [ fig|lfz6 ] .
the gray area shows the data with its @xmath109-error bars by f04 for the @xmath110 lf , while the solid dots show the extrapolation of the f04 lf to faint luminosities assuming a single slope of @xmath78 .
the presence of a significant downturn in the qso lf for magnitudes below @xmath74 is evident .
in particular the m05 observations are significantly stringent in constraining the faint - end slope .
ccccccc mahabal et al . & & 1.7 & 1.9 & 2.3 & 2.5 sq .
28 % + willott et al . & & 2.27 & 2.57 & 2.95 & 3.32 sq .
variable + barger et al . ( single source ) & & 2.32 & 2.38 & 2.48 & 160 arcmin@xmath76 & 50 % + [ table|survey ]    in figure  [ fig|lfz6 ] we also show the interesting result by stern _ et al .
_ ( 2000 ) .
these authors have discovered a quasar at @xmath111 from deep , multicolor , ground - based observations covering 74 arcmin@xmath76 with z@xmath112 . to estimate their average selection function we redo the calculations in stern _
_ ( 2000 ) and integrate down to @xmath113 ( in their cosmology ) the boyle _
( 1991 ) @xmath114 quasar lf , evolved as in schmidt _
the effective volume probed is again defined by the @xmath90-band drop - out as in m05 ( stern _ et al .
we find that an average completeness of p@xmath115 recovers their quoted @xmath1160.15 quasars in the field of view . as for m05
, we then use the more up - to - date fan _
( 2001 ) lf , more suitable for the redshift range probed by this survey , and equation  ( [ eq|integrallf ] ) . even in this case
we let the slope to be variable , @xmath117 in the specific case , while we keep the redshift evolution and normalization at @xmath73 as in fan _
( 2001 ) .
we integrate down to @xmath118 ( see table 1 in stern _
et al . _ 2000 ) , instead of using the measured z@xmath112 .
this is because the authors find their source to have unusual spectral features for which they derive a special continuum model which is significantly different from the average spectral template used in this paper .
the detection of one source in their field of view is consistent with a comoving density of @xmath119\sim -5.66\ , { \rm mpc^{-3}\ , mag^{-1}}$ ] , lower than the f04 extrapolation at the @xmath120 80% cl ( see figure  [ fig|lfz6 ] ) . the number density we find still shows evidence for a turnover in the qso lf expressed in equation  ( [ eq|modellf ] ) , with @xmath3 , though less significant than the above results , in agreement with the authors .
we now turn to the analysis of x - ray data .
( 2003 ; see also cowie _ et al .
_ 2003 ) have used deep multicolor optical data to search for @xmath81 agns in the 2 ms x - ray _ chandra _ deep field - north exposures .
they found 500 sources in their field of view out of which 423 have z@xmath121 .
color analysis was carried out on all the sources , while 249 sources have spectroscopically confirmed redshifts . of these
, only one object was found at @xmath122 and none at @xmath123 .
they found a faint agn at @xmath124 with z@xmath125 , corresponding to an absolute magnitude of @xmath126 , i.e. nearly a magnitude fainter than the m05 source , showing the power of x - ray observations .
( 2003 ) claim that none of the spectroscopically unidentified sources with z@xmath121 is red enough to have @xmath81 . with the aid of the @xmath127 method the authors estimated the lf in the comoving volume @xmath128 to be consistent with the extrapolation of the fan _
( 2001 ) lf with a slope of @xmath129 in the luminosity bin @xmath130 ( in @xmath131 ) . here
we proceed with a similar calculation but using equation  ( [ eq|integrallf ] ) and extrapolating the fan _
( 2001 ) lf at redshifts beyond @xmath18 , in the same manner we computed the number density for stern _
et al . _ ( 2000 ) .
we find that the number of observed sources down to z@xmath132 integrated between @xmath133 , with @xmath134 and a completeness of p@xmath135 ( cowie _ et al .
_ 2002 ) is @xmath136 , in very good agreement with observations .
we considered an area of 160 arcmin@xmath76 , corresponding to the solid angle where the hard x - ray flux limit has a signal - to - noise ratio of s / n=3 ( see figure 19 in alexander _
the x - ray flux considered here is @xmath137 ( alexander _ et al .
_ 2003 ; barger _
this flux limit is the minimum level of x - ray flux which corresponds to an optical magnitude brighter than z@xmath138 ( the optical detection limit in barger _
( 2003 ) ) within the whole effective volume considered .
this is required as all the 500 sources in the cdfn have optical colors .
we find a number density of @xmath139 @xmath140 @xmath141 $ ] @xmath142 corresponding to @xmath143 @xmath144 ( equation  ( [ eq|aoxvslandz ] ) ) , a factor of @xmath145 higher , but broadly consistent , with the results by barger _
( 2003 ; plotted in figure  [ fig|lfz6 ] as an open circle with its 99% cl error bars ) .
indeed deep x - ray observations effectively allow us to probe even fainter optical luminosities .
( 2002 ; brandt _ et al .
_ 2002 ) detected three @xmath0 quasars discovered by fan _
( 2001b ) with _
chandra _ with exposure times of less than 10 ksec .
the objects are detected in x - rays up to 55 kev in rest frame , demonstrating the incredible sensitivity of _ chandra _ to detect faint sources .
the faintest quasar in their sample , sdss 1030@xmath1460524 , at @xmath147 was detected in about 8 ksec with an observed flux of @xmath148 .
the _ chandra _
deep field observations have exposure time of 2 msec , so could detect down to @xmath149 times fainter sources with a flux limit of @xmath150
. this however would be the sensitivity limit at the aim point , close to the aim point sensitivity reported by alexander _
therefore we expect the _ chandra _ observations to be sensitive in detecting sources down the minimum flux considered , i.e. , we can extend the computation using equations  ( [ eq|modellf ] ) and  ( [ eq|integrallf ] ) down to the limit of @xmath151 ( i.e. down to z@xmath152 at the average redshift of @xmath153 ) . in figure
[ fig|lfz6 ] and table  [ table|survey ] we report our findings .
using a completeness of about 50% , as computed by cowie _
( 2002 ) for faint fluxes ( notice that the use of an higher completeness level would strengthen our results ) , we find that the detection of a single source in the _ chandra _ deep field in the range @xmath128 , is consistent with a faint slope of @xmath154 at the 99% cl , down to @xmath155 at the 90% cl ( shown as solid squares ) .
( 2003 ) claim that none of the unidentified sources is red enough to be at @xmath123 .
they however also point out that up to nine out of the optically faint spectroscopically unidentified sources with z@xmath156 , could still possibly lie in the range @xmath128 . therefore assuming these sources to rely in the same luminosity bin as the single @xmath124 source , they find an upper limit to their comoving lf about an order of magnitude higher than their previous estimate .
following the same approach undertaken for the case of a single source , we apply equations  ( [ eq|modellf ] ) and  ( [ eq|integrallf ] ) down to @xmath151 and using a completeness of @xmath157 , proper for faint fluxes ( cowie _ et al .
we consider all the nine sources not to be heavily obscured , i.e. they are type i agn .
we find that the comoving number density of the lf to be @xmath158\sim -4.52 \ , ( { \rm mpc^{-3}\ , mag^{-1}})$ ] , using a slope @xmath159 ( shown as an open square with its 99% cl in figure  [ fig|lfz6 ] ) . our estimate is higher but still consistent with the one by barger _
( 2003 ) . to notice also that our estimate is also quite conservative , lower limits would be inferred extending the computation to lower magnitudes , close to the survey limit ( z@xmath138 ) .
the 1-@xmath32 uncertainties in equation  ( [ eq|aoxlxlo ] ) produce uncertainties of @xmath160 magn in the final computation of the optical @xmath161 magnitude , as shown in figure  [ fig|lfz6 ] .
we also notice that our conclusions are quite insensitive to the actual fit used for the optical - to - x - ray flux ratios .
for example , similar results are also obtained using the relation between luminosity in the rest - frame @xmath90 band , @xmath162 , and the hard x - ray luminosity @xmath163 calibrated by la franca _ et al . _
( 2005 ; see their equation ( 4 ) ) for type i agns .
the vignali _
( 2003 ) fit instead of equation  ( [ eq|aoxlxlo ] ) would yield even fainter optical magnitudes for the same x - ray flux and number density , thus strengthening our results for a presence of a flattening at the faint - end of the qso lf .
we conclude this section noting that all the data are consistent with a break in the qso lf at @xmath164 .
the faint - end slope of the _ optical _ type i agn lf , is consistent with being on average @xmath165 , at the 90% cl and @xmath166 , at the 99% cl ( see figure  [ fig|lfz6 ] ) .
our overall findings are consistent with the fan _
( 2001 ) lf with @xmath167 extrapolated to @xmath104 ( shown with `` plus '' signs in figure  [ fig|lfz6 ] ) .
this clearly shows how the high - redshift lf is more consistent with a pde beyond @xmath18 .
( 2006 ) have recently found from a sample of sdss qsos , that the high - redshift qso lf bright slope flattens at @xmath168 from @xmath169 down to @xmath170 at @xmath18 .
we suggest instead that a double power - law , steeper at bright magnitudes and shallower at fainter ones , is more indicative for the qso lf at @xmath0 .
our results on the optical lf are instead in good agreement with those by fontanot _
( 2006a ) , who compiled the qso lf from a sample of goods and sdss qsos in the range @xmath171 .
they find significant evidence for a bright - end slope of @xmath172 , a break in the lf at @xmath173 and a strong flattening below this limit .      in the previous section we have considered the detected sources in the cdfn as type i agns .
the single high - redshift source detected with spectroscopic redshift of @xmath124 can in fact be safely considered as a type i agn . converting its measured hard x - ray flux @xmath174 to optical magnitudes ( using equation  ( [ eq|aoxvslandz ] ) ) we get z@xmath175 , very close to its actual magnitude of z@xmath176 .
however some of the nine optically faint spectroscopically unidentified objects found by barger _
( 2003 ) could be considered moderately obscured sources , with very faint optical fluxes ( though still detectable ) . in this section we give constraints on the faint - end shape for the agn lf , including the cdfn spectroscopically unidentified sources assuming that a certain fraction of them are type ii agns lying at @xmath0 .
indeed the @xmath177 band corresponds to the @xmath178 intrinsic band at @xmath179 , which is not much affected by absorption for column densities @xmath180 .
here we notice however that we do not expect the detection of any very obscured , compton - thick source in the cdfn , with @xmath181 . from wilman & fabian ( 1999 )
we have taken the spectrum of two typical agns with column densities @xmath182 and @xmath183 , derived from a monte - carlo realization which includes compton down - scattering which reduces the intensity even in the highest bins of energy .
we then computed the ratio between absorbed and unabsorbed luminosities ( the former is the actual _ observed _ flux ) for the various spectra in the @xmath179 redshifted , restframe energy hard band @xmath184 ( the least absorbed band of the cdfn ) .
we find that the ratio @xmath185 is @xmath186 for @xmath182 and @xmath9 for @xmath183 ( to notice for comparison that an agn with @xmath187 has a ratio of @xmath188 ) .
a typical source with intrinsic luminosity @xmath189 ( fainter luminosities would be below the x - ray flux limit ) in the observed frame will in turn produce observed fluxes @xmath190 and @xmath191 for respectively the @xmath182 and @xmath183 case .
therefore sources with @xmath192 @xmath193 and @xmath194 , would be hardly visible in the cdfn field of view ( we remind that the flux limit is @xmath195 ) . on the other hand
moderately obscured sources ( with @xmath196 ) would definitely be visible at these high redshifts , and some of them might have been already discovered in the deep _ chandra _ fields , among the optically faint barger _
( 2003 ) sample .
other possible indications of the presence of type ii sources in the high redshift universe also come from koekemoer _ et al . _
( 2004 ) . from the 2 msec hubble deep field - north and the 1 msec cdfs ,
these authors have robustly detected a sample of seven sources with extreme optical / x - ray flux ratios @xmath19710 at a flux limit above @xmath198 .
although these kind of sources are generally rare , if they are agns lying at @xmath123 such that their ly@xmath83 emission has been redshifted out of the z@xmath37 bandpass , then they would definitely be high - redshift obscured sources .
we can not directly make a prediction with their data as we are converting to optical luminosities and these are special sources , strong outliers with respect to equation  ( [ eq|aoxvslandz ] ) ( also a conversion to bolometric luminosity would be difficult in this case ) .    in this section
we still adopt the same x - ray @xmath199-correction as for unabsorbed sources with @xmath30 to convert flux to luminosities .
we still use equation  ( [ eq|aoxvslandz ] ) , which holds for type i agns , for converting x - ray fluxes to intrinsic optical luminosities for type ii agns , assuming that unification models are on average valid .
first we consider _ all _ the sources detected by barger et al .
( 2003 ) as optically faint type ii agns lying at @xmath81 .
we then integrate equation  ( [ eq|modellf ] ) in the redshift range @xmath133 above a flux limit of @xmath200 corresponding to an effective area of 40 arcmin@xmath76 ( see figure 19 , bottom panel in alexander et al .
the adopted flux limit corresponds to apparent magnitudes z@xmath201 at @xmath202 ( equation  ( [ eq|aoxvslandz ] ) ) . using @xmath203 we obtain @xmath204 sources in the field of view , corresponding to a comoving density of @xmath205\sim -3.85 \ , { \rm mpc^{-3}\ , mag^{-1}}$ ] .
the results are shown with an open square with its 99% confidence level in figure  [ fig|lfz6obscu ] with its error bar on the optical luminosity when using equation  ( [ eq|aoxvslandz ] ) .
if the nine sources are confirmed to be type ii agns at @xmath0 , this would show evidence for an extraordinary contribution of type ii agns at the faint end of the agn lf . assuming all these faint sources to lie at @xmath0
would imply a fraction obscured - to - unobscured agns of order 9 .
x - ray surveys at lower redshifts found much more moderate numbers for this ratio ( e.g. ueda _ et al .
_ 2003 , la franca _ et al .
_ 2005 , gilli _ et al .
_ 2006 ) .
in fact the above scenario is an extreme one , as most probably many of these optically faint sources probably belong to lower redshifts ( see discussions in barger _
_ 2003 and cowie _ et al .
therefore we more safely assume that up to 3 sources , out of the 9 detected , i.e. three times the surface density of the type i lf ( made up of only the single @xmath124 source ) , are actual contributors to @xmath0 agn lf with @xmath206 , following the results by la franca et al .
( 2005 ) , who find an almost flat distribution of the fraction of sources per logarithm bin of @xmath207 . in this second case
we get a comoving density of @xmath205\sim -4.3 \ , { \rm mpc^{-3}\ , mag^{-1}}$ ] , shown as a solid square in figure  [ fig|lfz6obscu ] with its 99% cl error bars .
we then compare our agn lf estimate with what would be obtained by multiplying the average qso lf ( with a faint slope of @xmath208 ; section  [ sec|typeilf ] ) by a factor @xmath209 , where @xmath90 is the ratio between type ii and type i agns , as calibrated by several groups for agns detected at lower redshifts . for computing @xmath90
we also include sources up to @xmath210 .
we use the results by ueda _
( 2003 ) , la franca _ et al . _
( 2005 ) and gilli _ et al . _
as specified above , we consider type ii agns all those sources with @xmath211 .
therefore by definition the number ratio type ii - to - type i will be given by @xmath212/n[\log ( n_h/{\rm cm^{-2}})<22]$ ] .
for ueda et al .
( 2003 ) and la franca et al .
( 2005 ) we use their quoted distributions for the fraction of objects up to @xmath213 .    gilli _ et al . _
( 2006 ) developed an updated model for the synthesis of the x - ray background , for which they adopted a different definition for the number ratio between obscured and unobscured agns , i.e. @xmath214/n[\log ( n_h/{\rm cm^{-2}})<21]$ ] . according to their model
this ratio is redshift - independent , luminosity - dependent and , translated to optical magnitudes , it reads as @xmath215 with @xmath216 , @xmath217 and @xmath218 , i.e. @xmath219 , where @xmath220 and @xmath221 are the @xmath222 $ ] @xmath223 $ ] ratio at the low- and high - luminosity regimes respectively , and @xmath218 is the characteristic luminosity dividing the two regimes ( gilli _ et al . _ 2006 ) . to include also compton - thick sources with column densities up to @xmath224
, we consider double this ratio ( see figure 7 in gilli _
_ 2006 ) .    the resulting agn lf from the gilli et al . ( 2006 ) model is shown with a dashed line in figure  [ fig|lfz6obscu ] . in the same figure we also plot the results of using @xmath90 as calibrated from ueda _ et al . _
( 2003 ; dot - dashed line ) and la franca _ et al . _
( 2005 ; long - dashed line ) .
all the models for the distribution of sources as a function of @xmath207 considered provide similar estimates for the total agn lf , and are consistent with the extrapolated data at the faint end , when using 3 sources as proxy for the @xmath0 agn population within the effective area .
we remind here that our estimate for the faint - end agn lf could easily be underestimated , as the cdfn can not detect compton - thick sources ( those with @xmath225 ) . to notice that all the models also provide a significant number of sources at the bright end of the agn lf in excess of the qso lf . due to the small volume probed by the cdfn we however can not provide any kind of constraint on the actual contribution of intrinsically very luminous and more rare obscured sources . anyway in the following we will consider as our best estimate of the total agn lf at @xmath0 , the qso lf multiplied by ( 1+@xmath226 ) , where @xmath226 is the type ii - to - type i ratio as extracted from ueda et al . (
2003 ) . in the following section we will compare theoretical models with such a determination of the agn lf
the high redshift agn lf is a very sensitive tool to probe the connections between supermassive black holes ( smbhs ) and their host dark matter ( dm ) halos . in this section
we discuss possible physical implications for smbh evolution from our estimated agn lf .
we compute the _ total _ agn lf from the formation rates of dm halos simply as ( e.g. , mahmood _ et al .
_ 2005 ) @xmath227 here @xmath228 is the formation rate of dm halos with mass @xmath229 at redshift @xmath230 , massive and virialized enough to host a galaxy / agn ; @xmath231 in this model is the average timescale during which an agn is active ( but not necessarily visible ) .
we neglect any dependence on mass for this timescale and just set it equal to @xmath232 .
@xmath233 is the smbh mass producing the luminosity @xmath234 . in the following for simplicity
we will use extended press - schechter ( eps ) theory to compute @xmath228 and we will consider only halos @xmath235 . the limiting mass of @xmath236 as host for the most massive galaxies and most luminous quasars has been constrained through clustering ( e.g. croom _ et al .
_ 2005 ) , statistical ( e.g. shankar _ et al .
_ 2006 ) , and semi - analytical ( e.g. granato _ et al .
_ 2004 ) techniques .
most probably the eps theory is a rather poor approximation for describing the statistics of high redshifts agn host halos , moreover a more physical treatment of the baryonic processes is required when comparing with observations ( e.g. lapi _ et al .
_ 2006 , fontanot _ et al . _
2006b ) . however here
we neglect such issues and focus on the most suitable relation @xmath237 , among those suggested in the literature , which better compares with observations within the very basic model expressed in equation  ( [ eq|philz ] ) .
we therefore consider different recipes for the jacobian relation in equation  ( [ eq|philz ] ) relating dm and smbh mass .
we first discuss the wyithe & loeb ( 2003 ; see also mahmood _ et al .
_ 2005 ) model which connects the smbh and dm halo masses as @xmath238 , derived from a basic agn feedback - constrained model
. we shall refer to this model as model a in the rest of the paper . using model a
, we calculate the @xmath104 lf , shown as a solid line in figure  [ fig|lfz6models1 ] .
this model is consistent with the agn lf estimated in section  [ sec|typeiilf ] .
several other groups have also formulated theoretical models for describing smbh evolution in the context of cosmological theories ( e.g. haenelt & rees 1993 , cavaliere & vittorini 2000 , vittorini _ et al .
_ 2005 , hopkins _
here we adopt a simple version of one of these models , viz .
the model by kauffmann & haenelt ( 2000 ; see also croton 2005 ) who parameterize the link between smbh mass and halo mass as : @xmath239^{-1}\ , m_{\rm cold}\ , ,      \label{eq|mbhmhalokh}\ ] ] where @xmath240 represents the merger mass ratio between the central and satellite galaxies , i.e. @xmath241 , @xmath242 is the virial velocity of the agn host halo and @xmath243 is the amount of cold baryonic matter ( which will eventually form stars ) .
the equation  ( [ eq|mbhmhalokh ] ) prescription with @xmath244 , consisting of only major mergers , is our model b.    we then consider three cases in which the cold fraction is @xmath245 , @xmath246 , and @xmath247 .
the resulting lfs are plotted in figure  [ fig|lfz6models2 ] .
the first case , in which the amount of baryons per halo is close to the cosmological value , lies well above the observed lf .
the second case fits the bright - end lf , assuming a baryonic fraction of a few percent , consistent with results by shankar _
it is interesting to note that a strong decrease of baryonic mass fraction with decreasing halo mass is needed to reproduce the agn lf , down to only 0.5% at the very faint end .
we now test a model based on ferrarese _ et al . _
( 2002 ; see also shankar _ et al .
we call this model c , and is the most empirical way to track the @xmath237 relation .
we relate the smbh mass to bulge velocity dispersion using the @xmath248 relation . the @xmath32 is then connected to the circular velocity of the host galaxy by using the empirical @xmath249 relation ( from baes _ et al .
the circular velocity is then scaled to the virial velocity assuming a navarro , frenk & white ( 1996 ; nfw ) profile and a concentration parameter @xmath250 from bullock _
finally the virial velocity is converted to halo mass using the virial theorem , @xmath251 ( we use the relations in loeb & peebles 2003 given for each @xmath252 ) , which is computed at the moment of the virialization of the host dm halo ( i.e. we set @xmath253 ) .
we also add a scatter of @xmath254 in the resulting @xmath237 relation , as empirically suggested by ferrarese ( 2002 ) and theoretically supported by lapi _
our result reads as @xmath255 which is close to the one calibrated locally through statistical arguments by shankar _
( 2006 ) and to the one obtained by granato _
( 2004 ) and lapi _ et al . _
( 2006 ) in semi - analytical models .
one uncertainty in the above computations is the radius at which one measures the circular velocity .
ferrarese ( 2002 ) uses @xmath256 ( also used in equation  [ eq|mbhmhaloc ] ) , where @xmath220 is the scale radius of the nfw dm profile . in figure  [ fig|lfz6models3 ]
we compare the predictions of model c with the calibrated high-@xmath230 agn lf .
the model matches the data when
@xmath257 and it has the advantage of being empirically - based , the only assumption here is that the @xmath248 and @xmath258 relations are almost constant in time .
our findings therefore differ from those in willott _
( 2005b ) , who attempted a similar calculation and found that the number of dm halos at @xmath0 falls short by several orders of magnitude in accounting for the number of high - redshift agns .
their conclusions were based on the use of the virial theorem calibrated at @xmath259 , which inevitably leads to significantly underestimating the bh mass ( and therefore agn luminosity ) for a given halo .
our findings are consistent with @xmath260 .
the latter relation in turn has been used by cirasuolo _
( 2005 ) to show that an impressive fit to the dispersion velocity function ( sheth et _ et al .
_ 2003 ; shankar _ et al .
_ 2004 ) is recovered when integrating the dm formation rates @xmath261 and setting @xmath262 ( see also loeb & peebles 2003 and shankar _ et al .
_ 2006 ) .
the simple approach undertaken in this section shows that a , redshift - dependent , super - linear relation of the type @xmath263 , with @xmath264 , as the ones proposed in model a and c , is best suited to represent the agn lf determined in section  [ sec|typeiilf ] .
alternative scenarios would require a very strong variation of baryonic content and/or agn activity timescale with halo mass .
in this section we compute the possible contribution of agns and stars to the reionization of the universe . in the following we will consider only type i agns as contributors of the reionization along our line of sight .
the transition from a neutral to a fully ionized inter galactic medium ( igm ) is usually statistically described by a differential equation for the time evolution of the volume filling factor of the medium @xmath265 .
the latter quantity quantifies the level of the igm _ porosity _ created by the hii , heii and heiii ionization regions around radiative sources such as qsos and galaxies .    when the clumping factor of the medium is high , @xmath266 , which is always a good approximation at high redshifts ( the simulations by gnedin & ostriker 1997 give @xmath267 ) , following madau _
( 1999 ) one can approximate @xmath268 , where @xmath269 is the volume - averaged mean recombination rate , @xmath270 is the total rate of ionizing photons and @xmath271 is the mean hydrogen density of the expanding igm .
the condition for reionization will then read @xmath272 , which translates into a total required ionizing photon rate ( madau _ et al .
_ 1999 ; fan _ et al .
_ 2001 ) : @xmath273    we compare such a rate with those from qsos and galaxies at similar redshifts .
we start computing the photoionization rate from qsos , which is given by ( e.g. madau _ et al .
_ 1999 , stiavelli _ et al .
_ 2004 ) @xmath274 where @xmath275 is the planck s constant , @xmath276 is the frequency at the lyman limit ( 912  , i.e. one rydberg ) .
the volume emissivity of the total population of sources is given as @xmath277 where @xmath278 is in @xmath279 .
the integration in equation  ( [ eq|enuagn ] ) is usually performed up to 4 rydbergs , as photons more energetic than that are preferentially absorbed by he atoms .
anyway we have checked that extending the computation to @xmath280 rydbergs for example , increases our final results by just @xmath281 .
@xmath282 in equation  ( [ eq|enu ] ) is instead the minimum luminosity probed by the survey , which we fix to @xmath283 , following the results in section  [ sec|typeilf ] .
we thus find that the agn contribution to reionization at redshift @xmath0 integrating down to @xmath283 , varies from @xmath284 to @xmath285 , using a faint - end slope @xmath66 of -2.2 , -2.8 respectively .
these results are rather independent of the details for the agn spectrum . for the computations above we have used the spectrum described in schirber & bullock ( 2003 ) , however using , for example , the three - power law spectrum adopted by madau _
( 1999 ) , reduces the agn contribution by @xmath286 .
we follow meiksin ( 2005 ) for computing the photoionization rate from galaxies @xmath287 here @xmath288 is the production rate of ionizing photons per solar mass of stars formed , @xmath289 is the star formation rate ( sfr ) in units of @xmath290 , and @xmath291 is the average fraction of ionizing photons which can freely escape the star - forming regions into the igm .
the photons per solar mass has been taken from smith _ et al . _
( 2002 ; see also meiksin 2005 ) for a salpeter initial mass function ( imf ) and a metallicity of 20% the solar .
the sfr was instead recently derived by richard _
( 2006 ) to be @xmath292 from a deep survey of lensing clusters aimed at constraining the abundance of star - forming galaxies at redshift @xmath293 .
we then take results from the recent paper by shapley _
( 2006 ) , who have constrained the escape fraction in the ionizing lyman continuum from a sample of 14 @xmath294 star - forming galaxies , finding an average @xmath295 , and we use the average value of @xmath296 . in the following
we neglect the diffuse emission from the igm itself which however could add an additional photoionization rate ( up to 50% , e.g. meiksin & madau 1993 ) over the direct one from the sources .
we find that galaxies can contribute about @xmath297 at @xmath104 to reionization .
we now take a step forward and integrate the contribution to reionization from both type of sources within the range of the _ chandra _ deep fields , i.e. @xmath298 , and compare the inferred number of ionizing photons with the number of available neutral hydrogen atoms at @xmath0 , @xmath299 ( see stiavelli _
we find that the contribution from agns to grow up to @xmath300 while the one from galaxies being about @xmath301 , keeping the sfr constant in this redshift bin , as suggested by observations .
the pde evolution of the agn lf for @xmath302 makes the contribution from agns grow more than the one from galaxies .
the photons produced by galaxies reach an average contribution to reionization of @xmath303 , implying a minor contribution from qsos , a conclusion in very good agreement with previous calculations ( e.g. meiksin 2005 , fan _ et al .
_ 2001 , barger _
et al . _ 2003 ) .
however we also notice that bouwens _
( 2005 ) from the deep nicmos fields in the hubble ultra deep fields ( hudf ) have detected 5 objects in 5.7 sq .
deg . , inferring a sfr about a factor of @xmath145 below the one found by richard _
this would drop the galaxy contribution down to @xmath304 .
similar points were also recently addressed by mannucci _
( 2006 ) who claim that the high - redshift sfr could fall short of completely reionizing the universe .
they however also point out that the bouwens _ et al . _
( 2005 ) estimate should be further increased by a factor of 2 for dust extinction and/or accounting for luminosity evolution in the estimate .
using therefore the bouwens _ et al . _
( 2005 ) result in equation  ( [ eq|enugalx ] ) would still leave room for a significant contribution from qsos .
other ways to fill the possible `` missing '' photons from galaxies could be to alter significantly the photon luminosity from massive stars , invoking extra igm diffuse light , or including other photoionizing sources ( such as population iii stars and/or miniquasars powered by intermediate - mass bhs ) .
however the contribution from the latter sources must be @xmath305 at the most .
in a recent paper , cool _ et al . _ . ( 2006 ) found three qsos at high redshifts ( @xmath306 ) down to magnitudes of @xmath307 .
the detection of these quasars is roughly consistent with the extrapolation of the bright steep slope of the quasar lf by f04 with @xmath308 . at lower magnitudes not many qsos
are observed , as recently concluded by the optical / nir surveys conducted by m05 , stern _
( 2005 ) and willott _ et al . _
( 2005 ) . in this paper
we show that such _
non_-detection of qsos implies a significant downturn in the optical qso lf at lower magnitudes with respect to those probed by cool _
moreover we show that the _ chandra _ deep field observations probe about @xmath309 magnitudes deeper than the optical ones ( down to @xmath310 ) . the inferred faint end co - moving number density is lower by at least an order of magnitude with respect to the extrapolation of the f04 lf , at a high significance level .
if a fraction of the nine _ chandra _ deep field sources are type ii agns detected at x - ray fluxes @xmath311 and with very faint optical magnitudes ( z@xmath156 ) , this would place a constraint on the shape of the type i plus type ii agns at faint magnitudes , being consistent with a slope close to the extrapolation of the f04 lf , i.e. @xmath312 .
this estimate of the agn lf is however representative of only the sources with @xmath313 , as more obscured sources would be below the cdfn flux limit .
x - ray surveys at lower redshifts ( e.g. la franca et al . 2005 ) have shown that an almost flat distribution in the number of sources is observed up to column densities of @xmath313 .
we can therefore safely assume that only 3 , out of the nine optically faint sources to be agns at @xmath0 .
this yields a total agn lf consistent with updated synthesis models for the xrbg .
we then compare this result with predictions from three simple physically and empirically motivated basic models for agn evolution in dm halos .
we generally find that a redshift - dependent relation of the kind @xmath314 , with @xmath264 , produces models more consistent with the agn lf .
such relations were also found from independent statistical studies by shankar _
( 2006 ) and rotation curve extrapolations in spirals and ellipticals ( ferrarese 2002 ) .
alternative models , such as varying the amount of cold gas fraction or agn active timescale in smaller halos , would require extreme fine tuning and a large variation in the parameters with dark matter halo mass .
here we also notice that if we reduce the agn duty - cycle just to the visible optical phase , i.e. substituting @xmath231 with @xmath315 in equation  ( [ eq|modellf ] ) , the match between models and the type i agn lf derived in section  [ sec|typeilf ] , would yield similar results on the type of @xmath237 relation .
finally we have shown that agns can account for at most @xmath316 of the total reionization of the universe at redshifts around @xmath317 .
however we also notice that a strong caveat on the actual role of galaxies in reionizing the universe , is set by the normalization of the cosmological sfr at high redshifts , which is still uncertain by factors of about @xmath318 .
we acknowledge c. s. kochanek , p. s. osmer , d. weinberg and p. martini for helpful discussions .
fs also acknowledges s. cristiani and f. fontanot for sharing data on the quasar template spectrum .
the authors are highly indebted to the referee for his / her very useful comments and suggestions .
this work was supported in part by nasa grant grt000001640 .
alexander , d. m. , et al .
2003 , aj , 126 , 539 baes , m. , buyle , p. , hau , g. k. t. , & dejonghe , h. 2003 , mnras , 341 , l44 barger , a. j. , et al .
2003 , aj , 126 , 632 boyle , b. j. , et al .
1991 , in asp conf .
21 , the space distribution of quasars , ed .
d. crampton ( san francisco : asp ) , 191 brandt , w. n. et al .
2002 , apj , 569 , l5 bullock , j. s. , kolatt , t. s. , rachel , y. s. , somerville , s. , kravtsov , a. v. , klypin , a. a. , primack , j. r. , & dekel , a. 2001 , mnras , 321 , 559 cavaliere , a. , & vittorini , v. 2000 , apj , 543 , 599 cimatti , a. , et al .
2002 , a&a , 392 , 395 cirasuolo , m. , shankar , f. , granato , g. l. , de zotti , g. , & danese , l. 2005 , apj , 629 , 816 cool , r. j. , kochanek , c. s. , eisenstein , d. j. , stern , d. , brand , k. , brown , m. j. i. , dey , a. , eisenhardt , p. r. , fan , x. , gonzalez , a. h. , green , r. f. , jannuzi , b. t. , mckenzie , e. h. , rieke , g. h. , rieke , m. , soifer , b. t. , spinrad , h. , & elston , r. j. 2006 , aj , in press , astroph/0605030 cowie , l. l. , et al .
1995 , nature , 377 , 603 cowie , l. l. , et al .
2002 , apj , 566 , l5 cowie , l. l. , et al .
2003 , apj , 584 , l57 cristiani , s. , & vio , r. 1990 , a&a , 227 , 385 cristiani , s. , et al .
2004 , apj , 600 , l119 croom , s. m. , et al .
2005 , mnras , 356 , 415 croton , d. j. 2005 , astroph/0512375 dijkstra , m. , & wyithe , j. s. b. 2006 , astroph/0606334 fan , x. , et al .
2001 , aj , 121 , 54 fan , x. , et al . 2001b , aj , 122 , 2833 fan , x. , et al .
2004 , apj , 128 , 515 ferrarese , l. 2002 , apj , 578 , 90 ferrarese , f. , & ford , h. c. 2005 , ssrv , 116 , 523 fiore , f. , et al .
2003 , a&a , 409 , 79 fontanot , f. , cristiani , s. , monaco , p. , nonino , m. , vanzella , e. , brandt , w. n. , grazian , a. , & mao , j. 2006a , a&a , accepted , astroph/0608664 fontanot , f. s. , monaco , p. , cristiani , s. , & tozzi , p. 2006b ,
mnras , accepted , astroph/0609823 gehrels , n. 1986 , apj , 303 , 336 gilli , r. , comastri , a. , & hasinger , g. 2006 , a&a , accepted , astroph/0610939 granato , g. l. , silva , l. , monaco , p. , panuzzo , p. , salucci , p. , de zotti , g. , & danese , l. 2001 , mnras , 324 , 757 granato , g. l. , de zotti , g. , silva , l. , bressan , a. , & danese , l. 2004 , apj , 600 , 580 haehnelt , m. g. , rees , m.j .
1993 , mnras , 263 , 168 hopkins , p. f. , hernquist , l. , cox , t. j. , di matteo , t. , robertson , b. , & springel , v. 2005 , apjs , in press ( astroph/0506398 ) kauffmann , g. , & haehnelt , m. 2000 , mnras , 311 , 576 kennefick , j. d. , djorgovski , s. g. , & de carvalho , r. r. 1995 , aj , 110 , 2553 la franca , f. , et al .
2005 , apj , 635 , 864 lapi , a. , shankar , f. , mao , j. , granato , g. l. , silva , l. , de zotti , g. , & danese , l. apj , in press , astroph/0603819 loeb , a. , & peebles , p. j. e. 2003 , apj , 589 , 29 mahabal , a. , stern , d. , bogosavljevic , m. , djorgovski , s.g . , &
thompson , d. 2005 , apj , 634 , l9 mahmood , a. , devriendt , j. e. g. , & silk , j. 2005 , mnras , 359 , 1363 mannucci , f. , buttery , h. , maiolino , r. , marconi , a. , & pozzetti , l. 2006 , astroph/0607143 mathur , s. , wilkes , b. j. , & ghosh , h. 2002 , apj , 570 , l5 monaco , p. , salucci , p. , & danese , l. 2000 , mnras , 311 , 279 navarro , j. f. , frenk , c. s. , & white , s. d. m. 1997 , apj , 490 , 493 ( nfw ) osmer , p. s. 2004 , in ho l.c . , ed . , carnegie observatories astrophysics series , vol .
1 , coevolution of black holes and galaxies . cambridge univ . press , cambridge , p. 324
pei , y. c. 1995 , apj , 438 , 623 richards , g. t. , et al .
2005 , mnras , 360 , 839 richards , g. t. , et al .
2006 , aj , in press ( astroph/0601434 ) shankar , f. , salucci , p. , granato , g. l. , de zotti , g. , & danese , l. 2004 , mnras , 354 , 1020 shankar , f. , lapi , a. , salucci , p. , de zotti , g. , & danese , l. 2006 , apj , 643 , 14 shapley , a. e. , steidel , c. c. , pettini , m. , adelberg , k. l. , & erb , d. k. 2006 , apj in press , astroph/0606635 sharp , r. g. , crampton , d. , hook , i. m. , & mcmahon , r. g. 2004 , mnras , 350 , 449 sheth , r. k. , & tormen , g. 2002 , mnras , 329 , 61l schirber , m. , & bullock , j. s. 2003 , apj , 584 , 110 spergel , d. n. , et al .
2003 , apjs , 148 , 175 springel , v. , et al .
2005 , nature , 435 , 629 schmidt , m , schneider , d. p. , & gunn , j. e. 1995 , aj , 110 , 68 steffen , a. t. , strateva , i. , brandt , w. n. , alexander , d. m. , koekemoer , a. m. , lehmer , b. d. , schneider , d. p , & vignali , c. 2006 , aj , 131 , 2826 stern , d. , et al .
2000 , apj , 533 , l75 stern , d. , et al . 2000b , aj , 119 , 1526 stiavelli , m. , fall , s. m. , & panagia , n. 2004 , apj , 600 , 508 tezca , m. , et al .
2004 , apj , 605 , l109 ueda , y. , akiyama , m. , ohta , k. , & miyaji , t. 2003 , apj , 598 , 886 vignali , c. , brandt , w. n. , & schneider , d.p .
2003 , aj , 125 , 433 vittorini , v. , shankar , f. , & cavaliere , a. 2005 , mnras , 363 , 1376 willott , c. j. , delfosse , x. , forveille , t. , delorme , p. , & gwyn , s. d. j. 2005 , apj , 633 , 630 willott , c. j. , et al . 2005b , apj , 626 , 657 wilman , r. j. , & fabian , a. c. 1999 , apj , 309 , 862 wolf , c. , et al .
2003 , a&a , 408 , 499 wyithe , j. s. b. 2004 , mnras , 351 , 1266 wyithe , j. s. b. , & loeb , a. 2003 , apj , 595 , 614


the diffuse spectrum of ultra - high energy ( uhe ) cosmic rays ( cr ) is expected to consist of a superposition of fluxes from many individual point - sources distributed throughout the universe .
the absence of significant event clustering across the sky sets a lower limit on the local source density or , equivalently , a lower limit on the number of sources that effectively contribute to the spectrum .
typically this number is very large and the uhe cr spectrum is expected to reflect the average contribution of these sources .
it is hence common practice for theoretical studies to approximate the distribution of cr sources via a spatially homogeneous and isotropic emission rate density that reproduces the average source spectrum .
this treatment greatly simplifies the study of uhe crs and can reproduce various spectral features from a simple power - law source spectrum .
the energy spectrum of cosmic rays follows a simple power - law over many energy decades .
small variations of the spectral index can be interpreted either as a transition between cr populations or as an imprint of cr propagation effects .
the _ ankle _  a hardening of the spectrum at @xmath1  ev  could be formed naturally by the superposition of two power - law fluxes and serves as a candidate of the transition between galactic heavy nuclei and extra - galactic cosmic ray protons  @xcite .
it has also been advocated that this feature could be well reproduced by a proton - dominated power - law spectrum , where the ankle is formed as a _ dip _ in the spectrum from the energy loss of protons via bethe - heitler pair production  @xcite . in this case
extra - galactic protons would already have started to dominate the spectrum beyond the _
2nd knee _ ,
a feature which corresponds to a slight _ softening _ of the spectrum at @xmath2 .
proton - dominance beyond the ankle is ultimately limited by the onset of photopion production on the cosmic microwave background ( cmb ) , whereas dominance of a heavy composition is restricted by nucleus photodisintegration through the giant dipole resonance ( gdr )  @xcite  the so called _ greisen - zatsepin - kuzmin _ ( gzk ) suppression at around @xmath3  ev .
indeed , a flux suppression in this energy region has been observed in the hires and auger data with a high statistical significance  @xcite .
as noted elsewhere  @xcite , secondary neutrinos and @xmath4-rays of these hadronic interactions can serve as additional discriminators between various cr models .    in order to collect the elusive events above @xmath5 ( which present an integrated flux of less than 1 event per km@xmath6 per steradian and per century ) observatories with large apertures and long exposure time
are needed .
today , the leading role in cr is played by ground based facilities that cover vast areas with particle detectors overlooked by fluorescence telescopes .
the largest is the pierre auger observatory , with a surface detector array of @xmath7 water cherenkov tanks covering @xmath8  km@xmath6 , which accumulates annually about @xmath9 of exposure  @xcite .
the more recently constructed telescope array ( ta ) covers 700  km@xmath6 with 507 scintillator detectors  @xcite , which should accumulate annually about @xmath10 of exposure .    in the near future
, the jem - euso mission will orbit the earth on board the international space station at an altitude of of about 400  km .
whilst in the `` nadir '' mode , the remote - sensing space instrument ( with @xmath11 field of view ) will monitor an area of approximately @xmath12 , recording video clips of fast uv flashes by sensing the fluorescence light produced through charged particle interactions .
this innovative pathfinder mission will observe approximately @xmath13 annually  @xcite , a factor of 10 above auger .
in this paper we elaborate on the question as to what extent the spectral information in the gzk region can be used to discriminate between different cr source composition models . due to the strength of the gzk mechanism the spectrum in this region
is dominated by ( and requires the presence of ) local sources  @xcite . in this case the flux from a few cr sources can significantly fluctuate from a homogeneous distribution that is typically assumed in cr flux predictions  @xcite . in contrast to poisson fluctuations in the gzk region
@xcite the manifestations of ensemble fluctuations persist in the limit of large event statistics .
we will quantify these stochastic fluctuations in the following utilizing an analytic solution to the flux of cr nuclei derived in refs .
@xcite .
we will start in sec .
[ sec : propagation ] with a brief review of the propagation of cr nuclei and the calculation of the mean observed fluxes . in sec .
[ sec : variation ] these results will be used to derive an analytic approximation of the flux and mean mass variations due to the distribution of sources .
we summarize our findings in sec .
[ sec : conclusion ] .
since cosmic rays are subject to deflections in galactic and inter - galactic magnetic fields the observed cr events do not point directly back to their sources .
the identification of the cr sources is hence experimentally challenging and has so far proved inconclusive .
there is a general consensus that the sources which are responsible of the uhe cr spectrum are of extra - galactic origin .
these sources are expected to follow a spatially homogeneous distribution and the mean ( ensemble - averaged ) flux of uhe crs ( of type @xmath14 ) follows a set of ( boltzmann ) continuity equations of the form : @xmath15 together with the friedman - lematre equations describing the cosmic expansion rate @xmath16 as a function of red - shift @xmath17 .
we follow the usual cosmological _ concordance model _ dominated by a cosmological constant with @xmath18 and a ( cold ) matter component , @xmath19 where , normalized to its value today of @xmath20 kms@xmath21mpc@xmath21  @xcite .
the time - dependence of the red - shift can be expressed via @xmath22 .
the first and second terms on the r.h.s .  of eq .
( [ eq : diff0 ] ) describe , respectively , red - shift and other continuous energy losses ( cel ) with rate @xmath23 .
the third and fourth terms describe more general interactions involving particle losses ( @xmath24 anything ) with total interaction rate @xmath25 , and particle generation of the form @xmath26 with differential interaction rate @xmath27 .
the last term on the r.h.s .  corresponds to the emission rate of crs of type @xmath14 per co - moving volume , depending on the emission rate @xmath28 per source and their density @xmath29 .
the two main reactions of uhe cr nuclei during their cosmic evolution are photodisintegration  @xcite and bethe - heitler pair production  @xcite with the cosmic radiation background .
in addition to the dominant contribution of the cmb we also include the infra - red / optical background from ref .
@xcite in our calculation of interaction and energy loss rates .
photodisintegration is dominated by the gdr with main branches @xmath30 and @xmath31 where @xmath32 indicates a proton or neutron  @xcite .
the gdr peak in the rest frame of the nucleus lies at at about @xmath33  mev for one - nucleon emission , corresponding to @xmath34  ev in the cosmic frame with photon energies @xmath35  mev . at energies below 10  mev there
exist typically a number of discrete excitation levels that can become significant for low mass nuclei . above 30  mev ,
where the photon wavelength becomes comparable or smaller than the size of the nucleus , the photon interacts via substructures of the nucleus . out of these
the interaction with quasi - deuterons is typically most dominant and forms a plateau of the cross section up to the photopion production threshold at @xmath36  mev .
bethe - heitler pair production can be treated as a continuous energy loss process with rate @xmath37 , where @xmath38 is the energy loss rate of protons  @xcite .
the ( differential ) photodisintegration rate @xmath39 ( @xmath40 ) is discussed in more detail in ref .
@xcite .
the evolution of the spectra proceeds very rapidly on cosmic time scales and the diffuse flux of secondary nuclei , @xmath41 , looks generally quite different from the initial source injection spectrum  @xcite .
the reaction network of nuclei depend in general on a large number of stable or long - lived isotopes .
if the life - time of an isotope is much shorter than its photodisintegration rate it can be effectively replaced by its long - lived decay products in the network  ( [ eq : diff0 ] ) .
typically , neutron - rich isotopes @xmath42-decay to a stable or long - lived nucleus with the same mass number . in most cases
there is only one stable nucleus per mass number below @xmath43fe with the exception of the pairs @xmath44cr/@xmath44fe , @xmath45ca/@xmath45ti , @xmath46ar/@xmath46ca and @xmath47s/@xmath47ar .
we follow here the approach of puget , stecker and bredekamp ( psb )  @xcite and consider only a single nucleus per mass number @xmath48 in the decay chain of primary iron @xmath43fe .
this psb - chain of nuclei linked by one - nucleon losses is discussed in more detail in ref .
@xcite .    note that the set of boltzmann eqs .
( [ eq : diff0 ] ) does not take into account the deflection of charged cr nuclei during their propagation through magnetic fields .
magnetic scattering can be viewed as a diffusion process that depends on the particle s larmor radius @xmath49 and the coherence length of the magnetic field . at energies where the diffusion length becomes larger than the distance to the nearest cr source
the cr spectrum is expected to be suppressed .
it has been speculated that for particularly strong inter - galactic magnetic fields of strength @xmath50  ng and coherence length of @xmath50 mpc ( see _ e.g. _  ref .
@xcite ) , the diffusive propagation of cr protons can start to affect the spectrum below about @xmath51  ev  @xcite . for heavier nuclei
diffusive propagation can in principle remain important up to higher energies due to the dependence @xmath52  @xcite .
the results of this paper assume that the contribution of inter - galactic or galactic magnetic fields can be neglected for the calculation of the uhe cr spectrum .
as an example we show in fig .
[ fig1 ] the solution of eq .
( [ eq : diff0 ] ) for an iron source model using a power - law spectral emission rate @xmath53 with @xmath54 and @xmath55  ev .
this model is motivated by a previous study  @xcite and reproduces the auger data above the ankle within systematic uncertainties .
the dashed black line corresponds to source contributions above @xmath56  mpc extending up to redshift @xmath57 where no red - shift evolution of the emission is assumed , _
i.e. _  @xmath58 .
the solid black line marks the local contribution up to 1  gpc calculated with the same method .
this local contribution where red - shift scaling of energies and interaction rates is suppressed can be approximated by an analytic solution which is shown as the green solid line in the plot . we will discuss this method in the next section .
gev and power index @xmath54 .
the upper solid ( green ) line shows the result of the analytic approximation ( [ eq : meann ] ) for a homogeneous distribution of sources @xmath59 with @xmath60  mpc and @xmath61  gpc . the lower solid ( black ) line is the corresponding numerical calculation based on eqs .
( [ eq : diff0 ] ) including redshift scaling of interaction rates and energies .
the relative difference between these calculations can be traced back to the onset of redshift scaling for propagation distances of the order of gpc ( see main text ) .
the dashed black line shows a model including sources up to @xmath57 .
also shown is recent data from hires  @xcite , auger  @xcite and the telescope array ( ta )  @xcite . ]
photodisintegration and photopion interactions of nuclei happen on timescales much shorter than the hubble scale .
since we are interested in the variation of the average flux and mass composition from the distribution of local cr sources we will neglect redshift scalings in the following . in this case the green s function of the boltzmann equations  ( [ eq : diff0 ] ) can be expressed in a simple analytic form , as discussed in refs .
@xcite .
it is convenient to study the flux of nuclei with mass number @xmath48 in terms of the energy per _ nucleon _ @xmath62 that is conserved under photodisintegration . introducing the cr flux @xmath63 per nucleon energy bin @xmath14 with bin - width @xmath64 ( centered at a nucleon energy @xmath65 ) and corresponding emission rates @xmath66 we find an analytic solution of the form
@xmath67 the sum on the r.h.s
.  runs over all possible production paths @xmath68 of a cr nucleus with final mass number @xmath48 and nucleon energy bin @xmath14 . each of the @xmath69 elements @xmath70 of the path
@xmath68 consists of a doublet @xmath71 denoting the mass number @xmath72 and nucleon energy bin @xmath73 of intermediate nuclei .
the first element @xmath74 of the paths corresponds to the nucleus emitted from a source at a rate @xmath75 ( possibly equal to zero depending on source composition ) and the last element is fixed at @xmath76 .
each path is weighted by a set of @xmath69 dimensionless amplitudes @xmath77 that are independent of the source distance @xmath78 . for more details
we refer to ref .
@xcite and appendix  [ sec : analytic ] .
we now want to study the statistical mean and variation of the aggregated flux of @xmath79 local cr sources denoted by @xmath80 as well as its corresponding mass composition .
herein we assume that the probability distribution function ( pdf ) for local ( @xmath81 ) sources is flat in euclidean space .
let @xmath79 sources be distributed between redshift @xmath82 and @xmath83 .
the number of sources can then be expressed via the ( local ) source density @xmath84 as @xmath85 .
the pdf of a single source is then given by @xmath86 note that the local density @xmath84 is limited by auto - correlation studies of uhe cr events ; the lower @xmath84 the larger the average emission rate of the sources and the greater the chance of local event clusters across the sky . a local source density of @xmath87  mpc@xmath88 is consistent with the absence of `` repeaters '' in cr data  @xcite .
values as low as @xmath89  mpc@xmath88 are still marginally consistent with auto - correlation studies of uhe cr nuclei  @xcite .
we will consider these two cases as fiducial values in the following .
following ref .
@xcite the ensemble - average of a quantity @xmath90 depending on the distance of the @xmath79 sources can then be expressed as @xmath91 the ensemble - average of the local flux of cr nuclei using eqs .
( [ eq : nai ] ) and ( [ eq : fullsol ] ) is then simply @xmath92 using the abbreviation ( [ eq : a ] ) we can then write this in an analytic form as @xmath93 where we define @xmath94    from the experimental point of view the interesting quantities are the ensemble - averaged total flux of nuclei @xmath95 and mean mass number @xmath96 at the highest cr energies @xmath97 . the mean total flux is simply given by eq .
( [ eq : mean ] ) as is a short - hand notation for @xmath98 with @xmath99 . ]
@xmath100 we now return to the example of an iron source model shown in fig .
the solid green line in fig .
[ fig1 ] shows the result of the analytic approximation  ( [ eq : mean ] ) .
this approximation agrees with the numerical result ( solid black line ) within a factor two or better depending on the energy .
the relative difference is expected from the onset of redshift scaling for propagation distances of the order of gpc ( @xmath101 ) . at low energies
this introduces a relative upward shift of the flux of @xmath102% for a source model with @xmath54 and @xmath58 .
this agrees well with the result of the calculations .
in addition , threshold effects that lead to breaks in the spectrum scale with redshift as @xmath103 and are shifted to lower energies by up to @xmath104% .
this effect can also be noticed by the relative position of the break in fig .
[ fig1 ] .
this example illustrates the limitations of the analytic solution for the calculation of large scale ( early - time ) contributions to the cr flux .
however , the analytic approximation provides a convenient description of the nuclei cascades that happen on small scales and depend on the local source distribution .
in particular , it enables us to study ensemble - variations of the flux . defining the residual @xmath105 , _ etc .
_ , we have as usual @xmath106 . using eqs .
( [ eq : fullsol ] ) and ( [ eq : nai ] ) we can then write the variation of the cr flux in an explicit analytic form as @xmath107 where @xmath108 is defined as the expression @xmath109 note that the last term in eq .
( [ eq : var ] ) is sometimes omitted since the number of sources @xmath79 is expected to be large , but we keep it in our calculations . based on these definitions we can express the relative variation of the total flux via the two - point density perturbations ( [ eq : var ] ) as @xmath110    with eq .
( [ eq : meann ] ) we can also _ define _ the mean mass number as @xmath111 note that eq .
( [ eq : meana ] ) is in the strict sense not the ensemble - average but serves as a first order estimator . for small fluctuations around the mean value we can approximate the relative variation of the mean mass number ( [ eq : meana ] ) via the two - point correlation function ( [ eq : var ] ) as @xmath112    in figure  [ fig2 ]
we show contour plots of the relative error of the flux ( red ; left plots ) and mass composition ( blue ; right plots ) for the case of iron sources with different model parameters indicated above the plots .
the axes show the observed cr energy vs.  the exponential cutoff energy @xmath113 . for the calculation we introduced logarithmic energy bins of @xmath114 and smoothed the result with a gaussian kernel to account for an experimental resolution of @xmath115 .
this procedure smoothes out features in the relative variance that are beyond the experimental resolution and result from the rapid nuclei transitions in the gzk region in combination with the relative cr energy shift with the mass of the daughter nuclei .
the results do not strongly depend on the spectral index as can be seen for the cases @xmath54 and @xmath116 with otherwise equal parameters that are shown in the plots of the first two rows .
this set of plots assumes a local distribution between 10  mpc and 1  gpc .
the relative errors are significantly reduced as we increase the distance to the closest source to 100  mpc as shown in the plots of the last row in fig .
however , this is only marginally reproducing the uhe cr spectrum as pointed out in ref .
note that the green lines in the plots of the top row mark the contribution for the iron source model with @xmath55  ev considered in fig .
[ fig1 ] .
one can also notice from the lower plots of fig .
[ fig2 ] that the relative error of the average mass composition is below 1% for cr energies below @xmath117  ev
. this energy marks the end of the energy region where cr event statistics allow an inference of the mass composition from cr data .
note that all the plots in fig .
[ fig2 ] show the case of a local source density of @xmath87  mpc@xmath88 and levels increase as @xmath118 .
hence , even for a density @xmath89  mpc@xmath88 , still marginally consistent with the data , the ensemble fluctuation on the average mass composition probed by present generation experiments may be safely neglected .
so far we have considered the case that @xmath119 valid at the highest cr energies where we can neglect cosmological contributions of the sources and treat the problem as effectively local . in the opposite case ,
@xmath120 , it is also possible to give an analytic expression of the statistical variation of the cr flux .
we assume that @xmath79 sources are isotropically distributed between redshift @xmath121 and @xmath122 with comoving density @xmath123 , and local density @xmath84 . the number of sources is then given by @xmath124 where the comoving volume is @xmath125 with comoving distance ( in a flat universe ) @xmath126 .
the pdf of a single source as a function of redshift is then a simple generalization of eq .
( [ eq : prob1 ] ) , @xmath127 as long as only adiabatic redshift scaling is involved the flux of a single source is given by @xmath128 .
assuming @xmath129 we then obtain a relative variation of @xmath130 ^ 2}{\left[\int{\rm d } z\,p(z)(1+z)^{-\gamma}/d^2_c(z)\right]^2}-\frac{1}{n_s}\,.\ ] ] we show values for @xmath79 and @xmath131 in table  [ tab1 ] .
obviously , the statistical variation of the average mass number is negligible in this regime .
.estimated source number and adiabatic variation defined in eqs .
( [ eq : ns ] ) and ( [ eq : adisig ] ) .
we show results for a power - law cosmological evolution as @xmath132 with @xmath133 and for the star formation rate ( sfr ) according to ref .
@xcite ( with cutoff @xmath134 ) .
in all cases we assume a local density of @xmath87  mpc@xmath88 . [
cols="^,^,^,^,^,^",options="header " , ]         the two limiting behaviors of the relative flux variation ( [ eq : locsig ] ) and ( [ eq : adisig ] ) motivates the following treatment of the overall relative flux variation . using the numerical result of the set of boltzmann eqs .
( [ eq : diff0 ] ) we can calculate the average flux for the local ( @xmath135  gpc ) , and global source distribution .
this corresponds to the solid and dashed lines , respectively , in fig .
the flux variation can then be approximated via the superposition @xmath136 where @xmath137 is calculated from the analytic approximation ( [ eq : locsig ] ) , @xmath138 the adiabatic limit ( [ eq : adisig ] ) for @xmath139  gpc ( @xmath140 ) and @xmath141 from the numerical evaluation .    in the left plot of fig .
[ fig3 ] we show the result of this procedure for the previous example of iron sources shown in fig .  [ fig1 ] . the solid line ( corresponding to the dashed line in fig .  [ fig1 ] )
shows the contribution from the cosmological distribution of iron sources .
the shaded areas show the flux within its variation based on eq .
( [ eq : combinedsig ] ) for a local source density @xmath87  mpc@xmath88 ( dark gray ) and @xmath89  mpc@xmath88 ( light gray ) .
note that the relative size of the error bands increases as @xmath118 .
the right plot of fig .
[ fig3 ] shows the corresponding result of a proton source model with similar source parameters .
the amplitude of the ensemble fluctuation is comparable to the case of iron and does not serve as a direct measure of the source composition .
however , the spectral feature of the gzk - suppression is significantly different to the case of the iron model .
in particular , the iron source model experiences a much steeper suppression at the upper end of the uhe cr spectrum .
this poses a challenge to account for the most extreme uhe events like the @xmath142  ev fly s eye event  @xcite , if sources are too distant  @xcite .
candidate nearby sources of heavy nuclei include starburst galaxies  @xcite and ultra - fast spinning newly - born pulsars  @xcite .
our method provides a tool to distinguish ensemble fluctuations from spectral features of the the average source contribution on a statistical basis .
proposed future space - based cr observatories can reach integrated exposures of @xmath143  @xcite .
this is a factor of about 100 larger than the integrated exposure reached by present ground - based air shower arrays .
the statistical error of uhe cr measurements can thus be improved by a factor 10 .
with such a resolution the observed uhe cr spectrum can show a significant deviation from the ensemble mean due to the discreteness of close - by source contributions .
for instance , the spectrum could exhibit `` spectral wiggles '' beyond the gzk suppression .
our method describes a way to quantify the bin - to - bin amplitude of these modulations for the case of a ultra - high energy cosmic ray nuclei .
in this paper we have studied the ensemble fluctuations of the mean flux and average mass number of uhe cr nuclei from the distribution of sources .
we have derived an analytic expression for the relative errors which applies to the cr data at the highest energies dominated by the local sub - gpc ( low redshift ) source distribution . for lower energies
another analytic approximation can be derived through the use of adiabatic scaling of the source contributions .
this method can be easily generalized to the case of ensemble fluctuations due to the source emission parameters such as the spectral index and the maximal energy .    as an illustration , we applied these results to a fiducial iron source model for which a homogeneous distribution of sources had previously been found to successfully reproduce the cr data within systematic uncertainties . for the case of a local source density @xmath87  mpc@xmath88 ,
the resultant ensemble fluctuations of the average mass composition on top of the ensemble mean were found to exist at a tolerable level for the analysis of present generation cr composition data below @xmath117  ev ( see right panels of fig.[fig2 ] ) .
the relative ensemble fluctuation around the mean flux increases with energy and rises above the level of 10% at about @xmath144 ( see left panels in fig .
[ fig2 ] and left plot in fig.[fig3 ] ) .
this flux variation is beyond the sensitivity of present day cr observatories as indicated by the range of the dark gray shaded bands in fig .
[ fig3 ] in comparison to the statistical errors of cr data .
the amplitude of the flux variation is similar in the case of proton models as shown in the right plot of fig .
however , the level of ensemble fluctuations for present generation spectral studies of the gzk suppression are potentially not ignorable if a smaller source density of @xmath89  mpc@xmath88 is assumed , which still remains marginally consistent with angular correlation studies assuming heavy nuclei .
this is shown as the light gray shaded bands in fig .
[ fig3 ] .
unless the actual source densities are much larger than those considered here , next generation experiments reaching accumulated exposures of @xmath143 should be sufficiently sensitive to potentially discern these fluctuations . for instance , the surface detector array of the auger observatory observed 25 events between @xmath145  ev to @xmath146  ev with an integrated exposure of @xmath147  @xcite . with an almost fifty times larger integrated exposure of future observatories the relative poisson error at this energy should drop below 3% and hence would be smaller than the ensemble fluctuation , @xmath148 , even for a local density of @xmath87  mpc@xmath88 , see fig .
[ fig2 ] .
we have seen that ensemble fluctuations of the uhe cr spectrum increase as we decrease the minimal distance @xmath82 to the extra - galactic source population , while keeping the local source density constant .
this behavior is of course expected since the relative abundance of local sources is small @xmath149 and hence more susceptible for fluctuation , but with a higher flux weight @xmath150 compared to more distant populations .
in fact , the ensemble variance formally diverges at _ all _ energies as we set @xmath151 .
this is a well - known effect in studies of galactic ensemble fluctuations  @xcite .
the minimal distance @xmath82 hence serves here as a regulator . even for a moderate minimal distance of 10mpc and a source density of @xmath152  mpc@xmath88
the ensemble fluctuation can reach large values @xmath153 in the gzk region as indicated in the left column of fig .
this is an indication that very few or just one local source can entirely dominate the spectrum at these energies .
some authors have considered even the extreme case that a local source like the radio galaxy centaurus a at a distance of 3 - 4  mpc can be responsible for the _ entire _ uhe cr spectrum above the ankle  @xcite .
we have checked that for @xmath154 and @xmath87  mpc@xmath88 the ensemble fluctuation below @xmath155  ev are less than 10% and this particular ensemble realization of extra - galactic sources to appear by chance is hence unlikely within our setup .
in addition , the missing anisotropy of uhe cr events challenges this interpretation .
possible caveats to this result could be large ensemble - fluctuations of intrinsic source properties that have been omitted in this study ( but could easily be included ) and/or strong deflections in inter - galactic magnetic fields  @xcite ( see further below ) .
all calculations in this study consider the continuous and isotropic emission of cr sources .
however , it is easy to extend the discussion for the case of episodal or anisotropic emission . if @xmath156 is the typical timescale of the source emission and @xmath157 the total experimental observation time we can define an effective local source density as @xmath158
this compensates the reduced continuous - equivalent emission rate @xmath159 of the single sources averaged over the duration of the experiment .
note that the emission rate density @xmath160 that is fixed by the observed cr spectrum is independent of this rescaling .
an analogous argument can be given for an anisotropic cr emitter with a preferred local ( but not global ) emission direction .
this applies to jet - like source emission like in blazars or gamma - ray bursts .
if the emission is concentrated in a cone @xmath161 we can again define a reduced effective local density @xmath162 that compensates the enhanced isotropic - equivalent emission rate @xmath163 .    we can hence account for time - variable or anisotropic sources by rescaling the local source density .
these two effects work in opposite directions .
time variability increases the number of sources contributing to the spectrum and lowers the amplitude of the ensemble fluctuation by @xmath164 . on the other hand
, anisotropy decreases the source number and increases the fluctuation by @xmath165 .
note , however , that autocorrelation studies of uhe cr constrain the number of sources and limit the effective source density @xmath166 . hence , the lower limit of @xmath167 used in this analysis does still apply .    in this context
it is interesting to consider the effect of magnetic fields on the propagation of uhe crs .
this has been so far neglected in this study . as mentioned before , a magnetic field ( regular or turbulent ) has no effect on a spatially homogeneous distribution of sources .
however , we break homogeneity by introducing a minimal distance @xmath82 in our calculation .
it was shown for the case of protons sources that diffusion in intergalactic magnetic fields can have a significant effect on the extra - galactic cr spectrum at lower energies .
the limited distance to the closest source introduces a low energy suppression , sometimes called `` anti - gzk - cutoff ''  @xcite . in this case
, ensemble - fluctuations at the low energies close the cr ankle can increase dramatically beyond the adiabatic value ( [ eq : adisig ] ) considered in this paper .    on top of these spectral variations
there is also an effect on the the effective source density @xmath166 . with magnetic diffusion the time - scale @xmath156 and emission cone
@xmath161 are expected to disperse .
eventually , for strong diffusion the effective local density @xmath166 is simply the true density @xmath84 .
since the diffusion coefficient does in general depend on the rigidity of the cr nuclei this effect would introduce an additional energy dependence of the variation .    in conclusion , future space - based observatories with colossal exposures , @xmath143 , will provide the required large statistics at the high - energy end of the cr spectrum , allowing identification of ensemble fluctuation from the gzk suppression features on a statistical basis . in combination with information on the arrival - direction distribution of crs and on the secondary fluxes of @xmath4-rays and neutrinos these spectral features can provide a coherent picture for an _ indirect _ determination of the uhe cr nuclear composition  @xcite and will naturally complement the current _ direct _ measurements of extensive air shower observables .
we would like to thank tom paul for valuable comments on the manuscript .
ma acknowledges support by a john bahcall fellowship for neutrino astronomy of the wisconsin icecube particle astrophysics center ( wipac ) and by the u.s
. national science foundation ( nsf ) under grants opp-0236449 and phy-0236449 .
laa is supported by the nsf under career grant phy-1053663 , the national aeronautics and space administration ( nasa ) grant 11-apra11 - 0058 , and the uwm research growth initiative .
amt acknowledges support by a schroedinger fellowship at dublin institute for advanced study .
any opinions , findings , and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of nasa or nsf .
the boltzmann equations  ( [ eq : diff0 ] ) describe the time - evolution of an isotropic flux of crs from a homogenous distribution of isotropic cr emitters . in the following we want to study the contribution of a single close - by ( @xmath81 ) source that is continuously emitting cr nuclei . as a first step we can hence neglect the redshift loss term @xmath168 in eqs .
( [ eq : diff0 ] ) and regard @xmath169 as the local density of cr particles . for the case of an inhomogeneous distribution of sources we have to add a cr convection term in eqs .
( [ eq : diff0 ] ) , that replaces the time - derivative for continuous emission . since we are not interested in the direction of the source we can imagine sitting at a center of a sphere with radius @xmath170 with surface emission rate spectrum @xmath171  @xcite .
the convection term in this spherically symmetric setup reduces then to @xmath172 .     between an initial configuration ( blue dot ) and a final configuration ( magenta dot ) including one - nucleon losses ( red solid arrows ) , two - nucleon losses ( red dotted arrows ) and continuous energy loss ( green horizontal arrows ) .
for the exact analytic solution  ( [ eq : fullsol ] ) all possible transition chains of this type are taken into account .
( figure adopted from ref .
@xcite.),height=192 ]    the secondary nuclei produced via photodisintegration carry approximately the same lorentz factor as the initial nucleus .
as already discussed in section  [ sec : variation ] it is hence convenient to express the energy of a nucleus with mass number @xmath48 as @xmath173 where @xmath174 denotes the energy _ per nucleon_. the differential interaction rate in the set of eqs .
( [ eq : diff0 ] ) corresponding to the production of a nucleus with mass number @xmath72 and energy @xmath175 from a nucleus with mass number @xmath176 and energy @xmath177 can be approximated as @xmath178 where @xmath179 is the partial width of the transition .    finally , introducing the binned cr flux @xmath63 , and corresponding emission rates , @xmath180 we can describe the point - source flux as a solution of the compact set of equations @xmath181 where we define the rates : @xmath182 note that eqs .
( [ eq : diff1 ] ) holds for nuclei heavier than beryllium .
we can easily compensate for the process @xmath183be @xmath184
@xmath185he + @xmath185he + n of the psb chain by re - defining @xmath186 for @xmath187 and @xmath188 for other nuclei .
similarly , the chains with nucleons as a final particle are re - weighted by the corresponding multiplicity in the case of intermediate nuclei lighter than @xmath183be ( see ref .
@xcite for further details ) .    in ref .
@xcite it was shown that the general solution of eq .
( [ eq : diff1 ] ) can be written in the form of eq .
( [ eq : fullsol ] ) . the first sum in eq .
( [ eq : fullsol ] ) runs over all possible production chains @xmath68 of nuclei with mass @xmath48 in the nucleon energy bin @xmath14 .
the intermediate configurations @xmath70 of this chain are parametrized by the tuple @xmath71 denoting the mass number @xmath72 and nucleon energy bin @xmath73 .
the links of the chain correspond to the partial width @xmath189 .
the dimensionless amplitudes in eq .
( [ eq : fullsol ] ) are defined as @xmath190 the total width @xmath191 is simply the total rate including photo - disintegration and cel terms . these amplitudes satisfy the identity @xmath192 independent of the path @xmath68 ( see refs .
@xcite ) .
an example of a chain is shown in fig .
the elements @xmath70 of the chain correspond to intermediate configurations of mass number and energy bin on a grid .
the first element @xmath74 ( blue dot ) corresponds to the cr particle emitted at the source ; the last element @xmath193 ( magenta dot ) is the observed ( final ) configuration @xmath194 .
the example shows one - nucleon disintegration ( solid red line ) , two - nucleon disintegration ( dashed red line ) and continuous energy loss ( green lines ) .
in general , the number of possible paths is quite large and the computation of eq .
( [ eq : fullsol ] ) numerically expensive . for the calculation in this study
we have applied two more approximations to the solution ( [ eq : fullsol ] ) to reduce the number of terms .
firstly , we reduced the cel ( last two terms on the r.h.s .
( [ eq : diff1 ] ) ) to the effective loss term @xmath195 this introduces a relative error of the order @xmath196 , which is small in the relevant energy region of @xmath197  ev . and ,
secondly , we considered only the one - nucleon loss chain of the psb approximation .
it was shown in refs .
@xcite that this is a good approximation for cr nuclei with a large mass number , which is the focus of this analysis .
j.  linsley , _ proceedings of icrc 1963 , jaipur , india _ , ( commercial printing press , ltd . ,
bombay , india , 1964 ) , pp .
c.  t.  hill and d.  n.  schramm , phys .
d * 31 * , 564 ( 1985 ) .
m. hillas , phys .
a * 24 * , 677 ( 1967 ) .
v.  berezinsky , a.  z.  gazizov and s.  i.  grigorieva , phys .
d * 74 * , 043005 ( 2006 ) [ arxiv : hep - ph/0204357 ] .
k.  greisen , phys .
rev .  lett .
* 16 * , 748 ( 1966 ) .
g.  t.  zatsepin and v.  a.  kuzmin , jetp lett .   * 4 * , 78 ( 1966 ) [ pisma zh .
* 4 * , 114 ( 1966 ) ] .
r.  abbasi _ et al . _ [ hires collaboration ] , phys .
lett .   *
100 * , 101101 ( 2008 ) [ arxiv : astro - ph/0703099 ] .
j.  abraham _ et al .
_ [ pierre auger collaboration ] , phys .  rev .
* 101 * , 061101 ( 2008 ) [ arxiv:0806.4302 [ astro - ph ] ] .
j.  abraham _ et al . _
[ pierre auger collaboration ] , phys .
b * 685 * , 239 ( 2010 ) [ arxiv:1002.1975 [ astro-ph.he ] ] .
m.  ahlers , l.  a.  anchordoqui and s.  sarkar , phys .
d * 79 * , 083009 ( 2009 ) [ arxiv:0902.3993 [ astro-ph.he ] ] .
d.  hooper , a.  m.  taylor and s.  sarkar , astropart .
phys .   * 34 * , 340 ( 2011 ) [ arxiv:1007.1306 [ astro-ph.he ] ] .
j.  abraham _ et al . _
[ pierre auger collaboration ] , nucl .
instrum .
a * 613 * , 29 ( 2010 ) [ arxiv:1111.6764 [ astro-ph.im ] ] .
t.  abu - zayyad _ et al .
_ [ telescope array collaboration ] , nucl .
instrum .
a * 689 * , 87 ( 2012 ) [ arxiv:1201.4964 [ astro-ph.im ] ] .
j.  h.  adams jr . _
_ , [ arxiv:1203.3451 [ astro-ph.im ] ]
. a.  m.  taylor , m.  ahlers and f.  a.  aharonian , phys .
d * 84 * , 105007 ( 2011 ) [ arxiv:1107.2055 [ astro-ph.he ] ] .
j.  miralda - escude and e.  waxman , astrophys .
j.   * 462 * , l59 ( 1996 ) [ arxiv : astro - ph/9601012 ] .
l.  a.  anchordoqui , m.  t.  dova , l.  n.  epele and j.  d.  swain , phys .
d * 57 * , 7103 ( 1998 ) [ arxiv : astro - ph/9708082 ] .
g.  a.  medina - tanco , astrophys .
j.   * 510 * , l91 ( 1999 ) [ arxiv : astro - ph/9810366 ] .
r.  aloisio , d.  boncioli and d.  boncioli , astropart .
phys .   * 35 * , 152 ( 2011 ) [ arxiv:1002.4134 [ astro-ph.he ] ] .
d.  de marco , p.  blasi and a.  v.  olinto , astropart .
* 20 * , 53 ( 2003 ) [ arxiv : astro - ph/0301497 ] .
d.  hooper , s.  sarkar and a.  m.  taylor , phys .
d * 77 * , 103007 ( 2008 ) [ arxiv:0802.1538 [ astro - ph ] ] .
m.  ahlers and a.  m.  taylor , phys .
d * 82 * , 123005 ( 2010 ) [ arxiv:1010.3019 [ astro-ph.he ] ] .
j.  beringer _
[ particle data group ] , phys .
d * 86 * , 010001 ( 2012 ) . f.  w.  stecker , phys .
* 180 * , 1264 ( 1969 ) .
j.  l.  puget , f.  w.  stecker and j.  h.  bredekamp , astrophys .
j.   * 205 * , 638 ( 1976 ) .
l.  n.  epele and e.  roulet , jhep * 9810 * , 009 ( 1998 ) [ arxiv : astro - ph/9808104 ]
. j.  p.  rachen , _ interaction processes and statistical properties of the propagation of cosmic - rays in photon backgrounds _ , phd thesis of the bonn university , 1996 .
e.  khan _ et al .
_ , astropart .
* 23 * , 191 ( 2005 ) [ arxiv : astro - ph/0412109 ] .
d.  allard , e.  parizot , e.  khan , s.  goriely and a.  v.  olinto , astron .
astrophys .   *
443 * , l29 ( 2005 ) [ arxiv : astro - ph/0505566 ] .
s.  goriely , s.  hilaire and a.  j.  koning , astron .
astrophys .
* 487 * , 767 ( 2008 ) [ arxiv:0806.2239 [ astro - ph ] ] .
g.  r.  blumenthal , phys .
d * 1 * , 1596 ( 1970 ) .
a.  franceschini , g.  rodighiero and m.  vaccari , astron .
astrophys .   * 487 * , 837 ( 2008 ) [ arxiv:0805.1841 [ astro - ph ] ] .
kronberg , rept .
phys .   * 57 * , 325 ( 1994 ) .
r.  aloisio and v.  s.  berezinsky , astrophys .  j.   * 625 * , 249 ( 2005 ) [ arxiv : astro - ph/0412578 ] .
d.  harari , s.  mollerach and e.  roulet , jhep * 9908 * , 022 ( 1999 ) [ arxiv : astro - ph/9906309 ] .
t.  abu - zayyad _ et al . _ [ telescope array collaboration ] , [ arxiv:1205.5067 [ astro-ph.he ] ] .
e.  waxman , k.  b.  fisher and t.  piran , astrophys .
j.   * 483 * , 1 ( 1997 ) [ arxiv : astro - ph/9604005 ] .
t.  kashti and e.  waxman , jcap * 0805 * , 006 ( 2008 ) [ arxiv:0801.4516 [ astro - ph ] ] .
h.  takami , s.  inoue and t.  yamamoto , astropart .  phys .   * 35 * , 767 ( 2012 ) [ arxiv:1202.2874 [ astro-ph.he ] ] . m.  a.  lee , astrophys .  j.  * 229 * , 424 ( 1979 ) .
a.  m.  hopkins and j.  f.  beacom , astrophys .
j.   * 651 * , 142 ( 2006 ) [ arxiv : astro - ph/0601463 ] .
h.  yuksel , m.  d.  kistler , j.  f.  beacom and a.  m.  hopkins , astrophys .
j.   * 683 * , l5 ( 2008 ) [ arxiv:0804.4008 [ astro - ph ] ] .
d.  j.  bird , astrophys .
j.   * 441 * , 144 ( 1995 ) [ arxiv : astro - ph/9410067 ] .
l.  a.  anchordoqui , g.  e.  romero and j.  a.  combi , phys .
d * 60 * , 103001 ( 1999 ) [ arxiv : astro - ph/9903145 ] .
p.  blasi , r.  i.  epstein and a.  v.  olinto , astrophys .
j.   * 533 * , l123 ( 2000 ) [ arxiv : astro - ph/9912240 ] .
k.  fang , k.  kotera and a.  v.  olinto , astrophys .
j.   * 750 * , 118 ( 2012 ) [ arxiv:1201.5197 [ astro-ph.he ] ] .
a.  santangelo and a.  petrolini , new j.  phys .   *
11 * , 065010 ( 2009 ) [ arxiv:0909.5370 [ astro-ph.he ] ] .
p.  abreu _ et al .
_ [ pierre auger collaboration ] , [ arxiv:1107.4809 [ astro-ph.he ] ] .
p.  blasi and e.  amato , jcap * 1201 * , 010 ( 2012 ) [ arxiv:1105.4521 [ astro-ph.he ] ] .
g.  bernard , t.  delahaye , p.  salati and r.  taillet , astron .
astrophys .
* 544 * , a92 ( 2012 ) [ arxiv:1204.6289 [ astro-ph.he ] ] .
p.  l.  biermann and v.  de souza , astrophys .
j.   * 746 * , 72 ( 2012 ) [ arxiv:1106.0625 [ astro-ph.he ] ] .
l.  a.  anchordoqui , h.  goldberg and t.  j.  weiler , phys .
d * 84 * , 067301 ( 2011 ) [ arxiv:1103.0536 [ astro-ph.he ] ] .
h.  yuksel , t.  stanev , m.  d.  kistler and p.  p.  kronberg , astrophys .  j.   * 758 * , 16 ( 2012 ) [ arxiv:1203.3197 [ astro-ph.he ] ] .
m.  ahlers , l.  a.  anchordoqui and a.  m.  taylor , in preparation