The luminosity function is a distribution function, specifically the distribution of luminosities of objects in a sample. Luminosity itself is difficult to measure, since the total or bolometric luminosity requires both accurate distance measurements and integration of each object's spectrum over all frequencies. Generally we measure specific luminosity, L, over a given band, or range of frequency, so the units of L are erg sec-1 Hz-1. Then the luminosity function is n(L), where n(L) dL is the number of galaxies with luminosity in the range L to L + dL; n(L) has units Mpc-3 (watt Hz-1)-1. The integral of n(L) dL over all luminosities is just the density of galaxies, n.
Note that if n(L) is described by a power law, for L to be finite requires that the power law index be less than -1 at high luminosities, and greater than -1 at low luminosities, where L is the total luminosity of all galaxies in a unit volume, i.e.,
A simple combination of power laws is often a useful approximate form to
assume for a luminosity function, i.e.,
What we know about luminosity functions generally comes from the
distribution of fluxes in a survey which is complete in some way. The
best is a volume-limited sample (meaning every galaxy in the volume
has been measured). This is rarely available, since all surveys have
some minimum detectable flux, Smin, which translates into a cutoff
distance rcut = sqrt(L / 4
Smin) which is a function of luminosity. So the sample
size is a function of luminosity. This bias can be corrected if we
know that the distribution of galaxies is homogeneous, i.e., that the
total density
is independent of position. In general n is a strong function of
position; it varies by several orders of magnitude between rich
clusters and voids; indeed, the study of this function is the main
subject of this course.
If the distribution of objects is homogeneous, so that n is a
constant independent of position, which presumably olds when we
average over very large scales, then we can easily evaluate some
simple integrals of the luminosity function which apply to a flux
limited sample. The number of objects brighter than the minimum flux
Smin is just given by
where we have simply interchanged the integration over volume and over
luminosity. (Note: the distance r, which on small scales is simply
cz / H0, generalizes in a Friedman universe with = 0, = 1 to
see Condon (1984a).
In the simplest case of ``standard candles'' the
luminosity function is
so that
In magnitude notation N>(Smin) Smin-3/2
becomes N<(m)
m0.6. The
``differential source count`` function, n(Smin) =
-dN> (Smin) / dSmin depends
on Smin to the minus 5/2 power, and the total flux from all
sources is
proportional to S
n(S) dS which diverges. This is a statement of
Olber's paradox, which can only be resolved in a Euclidian universe if
the luminosity function evolves, i.e., changes with time. The only
tracers for which we have complete flux limited samples reaching to
large redshifts (z > 1) are QSO's and radio sources, both of which
apparently have luminosity functions which evolve strongly on
cosmological time scales.
In a flux limited sample the contribution to n(s) of sources of
different luminosities is most easily seen in
von Hoerner's (1973)
``visibility function'':
whose dimensions are watt3/2 Hz-3/2
pc-3, which is usually converted to
Jy3/2. A plot of log
(L) vs. log L immediately shows what range of
luminosities contribute most to a flux limited sample, since
(see Condon 1984a,
b).
A useful integral of the luminosity function gives the median
distance to objects in a flux limited sample, r1/2, given by
where again we can interchange integration to get
where L1/2 = 4 Smin r1/22. This can easily be
evaluated for n(L) having the
simple form of equation 1. For example, if n2 = 0 and 1 = 0 [i.e., a
step function n(L)], we find simply that r1/2 is 0.64
sqrt(L0 / 4
Smin). Using
(L) equation 6 becomes simply
where x = log L1/2, which is the obvious median value of (L) when plotted vs. log L.