Hubble-tension toolkit: quadrature, σ-tension, distance modulus, Hubble's law
DefinitionHubbleTension_Defscosmologyhubble-constant
Definitions shared by the mission:
- , the combination in quadrature of two independent uncertainties (e.g. stat and sys);
- , the discrepancy between independent measurements and in standard deviations;
- for a distance in parsecs;
- parsecs;
- , Hubble's law with in km/s and given in parsecs (converted to Mpc).
Definition code
import Mathlib
/-!
# Definitions for the Hubble-tension mission
Source: W. L. Freedman and B. F. Madore, *Progress in Direct Measurements of the Hubble
Constant*, arXiv:2309.05618v2 (2023).
-/
namespace HubbleTension
/-- Combination in quadrature of two independent uncertainties (e.g. a statistical and a
systematic error quoted as `x ± s (stat) ± t (sys)`): `√(s² + t²)`. -/
noncomputable def quadSum (s t : ℝ) : ℝ := Real.sqrt (s ^ 2 + t ^ 2)
/-- The tension, in units of standard deviations ("σ"), between two independent measurements
`a ± σa` and `b ± σb`: `|a - b| / √(σa² + σb²)`. -/
noncomputable def tension (a σa b σb : ℝ) : ℝ := |a - b| / quadSum σa σb
/-- The distance modulus `μ = 5 log₁₀ (d / 10 pc)` of an object at distance `d` parsecs
(meaningful for `d > 0`). -/
noncomputable def distanceModulus (d : ℝ) : ℝ := 5 * Real.logb 10 (d / 10)
/-- The distance in parsecs corresponding to a distance modulus `μ`: `d = 10^(μ/5 + 1)`. -/
noncomputable def distanceFromModulus (μ : ℝ) : ℝ := (10 : ℝ) ^ (μ / 5 + 1)
/-- Hubble's law `v = H₀ d`: the Hubble constant (in km s⁻¹ Mpc⁻¹) inferred from a recession
velocity `v` (km s⁻¹) and a distance `d` in parsecs, i.e. `H₀ = v / (d / 10⁶)`. -/
noncomputable def hubbleConstant (v d : ℝ) : ℝ := v / (d / 1e6)
end HubbleTension
Source
W. L. Freedman and B. F. Madore, Progress in Direct Measurements of the Hubble Constant, JCAP (2023), arXiv:2309.05618v2, https://arxiv.org/abs/2309.05618, §§1, 4.1, 5, 10 (standard conventions of the field)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.