Luminosity function (astronomy)

The main sequence luminosity function maps the distribution of main sequence stars according to their luminosity. It is used to compare star formation and death rates rates, and evolutionary models, with observations. Main sequence luminosity functions vary depending on their host galaxy and on selection criteria for the stars, for example in the Solar neighbourhood or the Small Magellanic Cloud.[2]

The white dwarf luminosity function (WDLF) gives the number of white dwarf stars with a given luminosity. As this is determined by the rates at which these stars form and cool, it is of interest for the information it gives about the physics of white dwarf cooling and the age and history of the Galaxy.[3][4]

The Schechter luminosity function[5] provides a parametric description of the space density of galaxies as a function of their luminosity. The form of the function is

${\displaystyle n(L)\ \mathrm {d} L=\phi ^{*}\left({\frac {L}{L^{*}}}\right)^{\alpha }\mathrm {e} ^{-L/L^{*}}{\frac {\mathrm {d} L}{L^{*}}},}$

where ${\displaystyle L}$ is galaxy luminosity, and ${\displaystyle L^{*}}$ is a characteristic galaxy luminosity where the power-law form of the function cuts off. The parameter ${\displaystyle \,\!\phi ^{*}}$ has units of number density and provides the normalization.

Equivalently, this equation can be expressed in terms of log-quantities with

${\displaystyle n(L)\ \mathrm {d} \left(\log _{10}L\right)=\ln(10)\phi ^{*}\left({\frac {L}{L^{*}}}\right)^{\alpha +1}\mathrm {e} ^{-L/L^{*}}\mathrm {d} \left(\log _{10}L\right).}$

The galaxy luminosity function may have different parameters for different populations and environments; it is not a universal function. One measurement from field galaxies is ${\displaystyle \alpha =-1.25,\ \phi ^{*}=1.2\times 10^{-2}\ h^{3}\ \mathrm {Mpc} ^{-3}}$.[6]

It is often more convenient to rewrite the Schechter function in terms of magnitudes, rather than luminosities. In this case, the Schechter function becomes:

${\displaystyle n(M)\ \mathrm {d} M=(0.4\ \ln 10)\ \phi ^{*}\ [10^{0.4(M^{*}-M)}]^{\alpha +1}\exp[-10^{0.4(M^{*}-M)}]\ \mathrm {d} M.}$

Note that because the magnitude system is logarithmic, the power law has logarithmic slope ${\displaystyle \alpha +1}$. This is why a Schechter function with ${\displaystyle \alpha =-1}$ is said to be flat.

Integrals of the Schechter function can be expressed via the incomplete gamma function

${\displaystyle \int _{a}^{b}x^{\alpha }e^{-x}\mathrm {d} x=\Gamma (\alpha +1,a)-\Gamma (\alpha +1,b)}$