How this instrument works
This instrument sums starlight the way Heinrich Olbers did in 1823: divide space into thin concentric shells centered on the observer, work out how much flux each shell contributes, and add the shells together out to some outer radius. A single star at distance r delivers flux L ⁄ (4πr²) at the observer, the inverse-square law. A shell at that radius, thickness dr, contains n · 4πr² · dr stars if the density n is uniform. Multiply the two and the 4πr² cancels exactly: each shell contributes n · L · dr, independent of how far away it sits.
That cancellation is the entire point. Integrate a constant shell contribution from 0 out to radius R and the result is F = n · L · R, flux growing in direct proportion to how far out the sum runs, not converging to some finite ceiling the way a single light source's brightness does. Push R toward the edge of the observable universe and the number stops looking like a modest astronomical quantity and starts looking absurd, which is exactly the point Olbers was making: a universe that is infinite, static, and eternally starlit would flood every line of sight with a star's surface, and the sky would glow rather than stay dark.
The formula's honesty has a boundary. It treats stars as point sources that never block one another, ignores that real starlight is absorbed by interstellar dust and redshifted by cosmic expansion, and assumes the universe has existed forever at constant density. None of that holds. The universe is about 13.8 billion years old, so light from anything farther than roughly that many light-years in travel time has not reached us yet, and space itself has been expanding the whole time, stretching distant starlight toward wavelengths the eye cannot use. Those two facts, not an error in the arithmetic here, are what keep the actual night sky dark.
- Enter the Average stellar density in stars per cubic light-year — how crowded the modeled region of space is.
- Set the Average star luminosity in watts; the Sun's nominal value is about 3.828×10²⁶ W for Sun-like stars.
- Set the Observable radius in light-years — how far the shell-by-shell sum runs before it stops.
- Read the Estimated total starlight flux in W ⁄ ly², the summed brightness of every star out to that radius.
Worked example — summing starlight out to 13.8 billion light-years
Set the Average stellar density to 0.14 stars ⁄ ly³, the Average star luminosity to the Sun's nominal 3.828×10²⁶ W, and the Observable radius to 13,800,000,000 light-years, a light-travel distance matching the universe's age. The formula multiplies straight through: F = 0.14 × 3.828×10²⁶ × 1.38×10¹⁰. The first two factors give 5.3592×10²⁵, and multiplying that by the radius gives 7.395696×10³⁵ W ⁄ ly², the Estimated total starlight flux the instrument reports.
That number has no everyday counterpart, and it is not meant to. It is what a naive, infinite, static, uniformly starlit universe would deliver to every square light-year of sky, and it is enormous precisely because the assumptions behind it are wrong. The real observable universe is finite in age and expanding, so most of that summed light either has not arrived yet or has redshifted out of the visible band, which is the actual, physical resolution of the paradox this calculator is built to demonstrate.
Questions
What is Olbers' Paradox?
If the universe were infinite, static, and uniformly filled with stars forever, every line of sight would eventually end on a star's surface, so the whole night sky should glow like a stellar surface rather than stay dark. Heinrich Olbers formalized the puzzle in 1823, though earlier astronomers including Johannes Kepler had already noticed the tension between an endless starry universe and a dark sky. The mismatch between that prediction and plain observation is the paradox.
Why does the formula multiply density, luminosity, and radius instead of dividing by distance squared?
Because the inverse-square dimming of one star exactly cancels the growth in shell volume as distance increases: a shell at radius r holds n · 4πr² · dr stars, each individually dimmed by 1 ⁄ (4πr²), so the 4πr² terms cancel and every shell contributes the same n · L · dr regardless of distance. Summing equal contributions from 0 to R gives a total proportional to R itself, which is why the formula is a plain product rather than a ratio.
Why doesn't the real night sky match this calculator's enormous number?
Because the real universe breaks both assumptions this formula makes. It has a finite age of about 13.8 billion years, so light from stars farther than that light-travel distance has not had time to reach us, capping how many shells actually contribute. It is also expanding, which redshifts light from distant, receding stars out of the visible spectrum before it arrives. This calculator deliberately omits both effects so the paradox shows up plainly in the raw arithmetic.
Does the stellar density input have to be a realistic cosmic average?
No. The field accepts whatever density you want to test, and the instrument's job is to show how sharply the total scales with that choice, not to certify one number as correct. Real stellar density varies enormously by location, from under 10⁻⁴ stars ⁄ ly³ in the thin space between galaxies to far higher figures inside a crowded galactic disk; try a few values and watch how linearly the output responds.
What happens to the total if I use a smaller radius, like the size of the Milky Way?
The total falls in direct proportion, because F = n · L · R is linear in R. Swap the observable-universe radius of 13.8 billion light-years for the Milky Way's roughly 100,000-light-year diameter, keeping density and luminosity fixed, and the flux drops by that same factor of about 138,000 — a useful way to feel how much of the golden result comes from pushing the outer radius out to cosmological distances.
Is the total starlight flux this calculator reports the actual brightness of the sky?
No. It is the sum a flat, static, infinitely old universe would produce, not a measurement or a physically achievable prediction. Actual sky brightness at night is many orders of magnitude lower, dominated by the finite number of stars whose light has had time to reach us plus faint contributions from zodiacal light and airglow — quantities this simplified shell integral does not attempt to model.