Luminosity Function Calculator

Calculator Inputs

Form

Choose how the Schechter function is evaluated.

φ* (Mpc⁻³)

Normalization of the luminosity function.

α (faint-end slope)

Typical ranges: −0.5 to −1.6.

L* (same units as L)

Characteristic luminosity scale.

Input Luminosity L

Value where φ(L) is evaluated.

Lmin for n(>Lmin)

Lower cutoff used in the integral.

M* (absolute magnitude)

Characteristic magnitude scale.

Input Magnitude M

Value where φ(M) is evaluated.

Magnitude limit Mlim (brighter than)

Computes n(M<Mlim), i.e., brighter objects.

Compute L/L☉ from magnitude?

Uses L/L☉ = 10^{0.4(M☉ − M)}.

M☉ (solar absolute magnitude)

Default 4.83 (bolometric). Choose your band if needed.

Tip: You can enter both sets of fields. The selected form controls what is computed. For stable integrations, keep Lmin or Mlim meaningful for your dataset.

Formula Used

Schechter luminosity form

φ(L) = (φ*/L*) (L/L*)^α exp(−L/L*)

Here φ* sets normalization, L* is the characteristic luminosity, and α controls the faint-end slope.

Schechter magnitude form

φ(M) = 0.4 ln(10) φ* 10^{0.4(α+1)(M*−M)} exp(−10^{0.4(M*−M)})

Magnitudes map to luminosity ratio x = 10^{0.4(M*−M)} = L/L*.

Cumulative number density

n(>Lmin) = φ* Γ(α+1, Lmin/L*)

Γ(a, x) is the upper incomplete gamma function. The calculator shifts parameters safely to handle common α values.

How to Use This Calculator

Select Luminosity form φ(L) or Magnitude form φ(M).
Enter φ* and α. These control normalization and slope.
For luminosity form, set L*, the target L, and Lmin.
For magnitude form, set M*, the target M, and Mlim.
Press Calculate. Results appear above the form.
Use Download CSV or Download PDF to export your results.

Example Data Table

Form	φ*	α	L* or M*	Input	Typical Output
Luminosity	0.004	−1.1	L* = 1×10¹⁰	L = 5×10⁹	φ(L) ~ few ×10⁻¹³ (per unit L)
Magnitude	0.004	−1.1	M* = −20.5	M = −19	φ(M) ~ 10⁻³ to 10⁻² (mag⁻¹)
Magnitude	0.004	−1.1	M* = −20.5	Mlim = −18	n(M<Mlim) shown as Mpc⁻³

These are illustrative. Real surveys use band-specific M☉ and carefully chosen limits.

Professional Article

1) Why luminosity functions matter

A luminosity function describes how many astrophysical sources exist per unit luminosity or per unit magnitude. In galaxy surveys, it links observed counts to underlying formation physics and selection effects. A well-fit function supports comparisons across redshift, environment, or wavelength band.

2) Schechter parameters and interpretation

The Schechter form is widely used because it captures a power-law rise at faint luminosities and an exponential cutoff at the bright end. The faint-end slope α shapes low-luminosity behavior, φ* sets the overall density scale, and L* (or M*) marks the transition region where the exponential suppression becomes important.

3) Luminosity versus magnitude input

This calculator supports both φ(L) and φ(M). The magnitude form uses the mapping x = 10^{0.4(M*−M)} = L/L*. Small changes in magnitude imply multiplicative changes in luminosity, so consistent band definitions and k-corrections are essential when comparing studies.

4) Typical parameter ranges in practice

For many low-redshift optical samples, α commonly lies near −1.0 to −1.3, indicating a rising population of faint objects. Characteristic magnitudes often fall around M* ≈ −20 to −21 in common bands, while φ* values are frequently of order 10^{-3} to 10^{-2} Mpc^{-3}, depending on selection and cosmology choices.

5) Evaluating φ at a chosen point

When you enter a target luminosity L or magnitude M, the tool computes φ(L) or φ(M) at that point. The dimensionless ratio x = L/L* (or x = 10^{0.4(M*−M)}) helps you judge where you sit on the curve. Values near x ≈ 1 probe the transition regime, while x ≪ 1 emphasizes the faint-end slope.

6) Cumulative counts by integration

Survey questions often focus on “how many sources above a threshold.” The calculator integrates the Schechter form using the upper incomplete gamma function: n(>Lmin) = φ* Γ(α+1, Lmin/L*), and similarly for magnitude limits. This approach is standard for converting fitted parameters into number densities above a defined cut.

7) Practical data hygiene and units

Keep units consistent. If L and L* are in solar luminosities, then x is unitless and φ(L) scales as Mpc^{-3} per luminosity unit. If you work in magnitudes, use band-appropriate solar absolute magnitude M☉ when computing L/L☉. Also record the assumed h or cosmological parameters used to quote φ* and M*.

8) Reporting results and exports

Professional reporting benefits from transparent metadata: the chosen form, parameter values, thresholds, and resulting densities. Use the built-in CSV export for spreadsheets and reproducible pipelines, and the PDF export for quick sharing in notes or proposals. For publications, consider quoting uncertainties from your fit alongside these derived quantities.

FAQs

1) What does α control?

α sets the faint-end slope. More negative α means more faint objects relative to bright ones. It strongly influences integrated counts when thresholds probe far below L*.

2) Why do I need L* or M*?

L* or M* defines the turnover between the power-law and exponential regimes. It anchors the scale where bright objects become increasingly rare.

3) What is φ* physically?

φ* is a normalization with units of number density (typically Mpc⁻³). It scales the overall abundance of sources for a given fitted shape.

4) Why is the cumulative density written with Γ?

The Schechter integral reduces to an upper incomplete gamma function when expressed in x = L/L*. This provides a stable, standard way to compute n(>Lmin) from fitted parameters.

5) My φ(L) looks extremely small. Is that normal?

Yes. In luminosity form, φ(L) is “per unit luminosity,” so its numeric value can be tiny when luminosity units are large. Compare using x = L/L* or convert to binned counts.

6) Should I use bolometric M☉ = 4.83?

Use the solar absolute magnitude appropriate to your bandpass. 4.83 is common for bolometric or V-like contexts, but survey bands often use different M☉ values.

7) What if α < −1 causes divergence?

The total number density diverges only if you integrate to L→0. Real surveys apply a finite Lmin or completeness limit, which makes n(>Lmin) well-defined and physically meaningful.