--- title: ZTF Dipole Investigation of Cosmic Expansion url: https://www.emergentmind.com/topics/ztf-dipole-investigation type: topic --- # ZTF Dipole Investigation of Cosmic Expansion The ZTF Dipole Investigation encompasses the methodological and statistical framework for testing cosmic expansion isotropy through Type Ia supernovae (SNe Ia) observed by the Zwicky Transient Facility (ZTF). It is motivated by growing interest in potential large-scale anisotropies that would violate the cosmological principle, with the specific focus on directional variations in the Hubble constant, $H_0$. Using ZTF’s wide sky coverage and realistic SN simulations, the approach develops, injects, and detects a velocity dipole in the expansion rate, quantifying both amplitude and sky direction, while rigorously controlling for systematics linked to cadence, calibration, and peculiar velocities [2601.11139]. ## 1. Theoretical Framework for Hubble Flow Anisotropy In this framework, the local Hubble constant is modeled as a directional quantity: $$ H_0'(n̂) = H_0 + ΔH_0 (n̂·p̂) $$ where $H_0$ is the isotropic (monopole) value, $ΔH_0$ the dipole amplitude, $n̂$ the supernova’s sky direction, and $p̂$ the dipole’s maximum. In equatorial coordinates $(α, δ)$, $$ \cos Δθ = \sin δ \sin δ_0 + \cos δ \cos δ_0 \cos(α - α_0) $$ so that $$ H_0'(α, δ) = H_0 + ΔH_0 \cos Δθ. $$ The luminosity distance in flat $\Lambda$CDM incorporating this anisotropy is: $$ d_{L}'(z; α, δ) = \frac{c(1 + z)}{H_0'(α, δ)} \int_0^z \frac{dζ}{\sqrt{(1+ζ)^3Ω_{m} + (1−Ω_{m})}} $$ with the corresponding distance modulus: $$ μ_{\mathrm{cosmo}}' = 5 \log_{10}[d_L'/(1\,\mathrm{Mpc})] + 25. $$ This dipole parameterization provides a rigorous test of the cosmological principle at scales probed by ZTF SNe Ia. ## 2. Dipole Injection in ZTF Supernova Simulations Mock ZTF catalogs are synthesized using the skysurvey package and SALT2.4 light-curve model, emulating the cadence, depth, and spectroscopic selection of ZTF’s DR2.5. Two main simulation sets are created: "unclustered" (SNe Ia placed isotropically, no large-scale structure) and "clustered" (positions and peculiar velocities from the Uchuu N-body halo catalog, with 27 sub-boxes to z_max ≈ 0.11 and random Gaussian velocities with σ=250 km/s beyond). Three dipole injection schemes are compared: - **d_L–method**: Recalculates $d_L'$ using $H_0'$ in the integral (computationally exact, but impractical for multiple directions). - **z–method**: Replaces $z$ with $z' = z (H_0'/H_0)$ from cz=H₀d (simple but produces a ≥6% bias even at $z \approx 0.06$). - **m_B–method**: Shifts the simulated apparent magnitude $m_B$ post hoc: $$ m_B' = m_B + 5 \log_{10}[H_0/H_0'(α, δ)] $$ This is algebraically equivalent to shifting $μ$ and remains unbiased at <0.5% while being highly efficient. Subsequent analyses adopt the m_B–method with a reference injected $ΔH_0 = 3$ km s⁻¹ Mpc⁻¹ and various dipole orientations. ## 3. Dipole Recovery: Fitting Procedures and Likelihood Construction Recovered distances are standardized via the Tripp relation: $$ μ_{\mathrm{obs}} = m_B - M_B + α x_1 - βc $$ where $(x_1, c, m_B)$ are SALT2 parameters, and $(M_B, α, β, σ_{\mathrm{int}})$ are nuisance parameters. The dipole recovery likelihood is: $$ \mathcal{L}(Θ) = -\frac{1}{2} \sum_{i=1}^N \left[\frac{μ_{i,\mathrm{obs}}(Θ) - μ_{i,\mathrm{cosmo}}(Θ)}{σ_i}\right]^2 - \frac{1}{2} \sum_{i} \ln(2\pi σ_i^2) $$ with $Θ = \{M_B, α, β, σ_{\mathrm{int}}, ΔH_0, α_0, δ_0\}$. To resolve the well-known $M_B \| H_0$ degeneracy, the minimization is performed in four structured steps: 1. Fit $(M_B, α, β, σ_{\mathrm{int}})$ with no dipole. 2. Coarse scan over 12 Healpix directions, fitting only $ΔH_0$ to localize $p̂$. 3. Full parameter fit $(M_B, α, β, ΔH_0, α_0, δ_0)$ from previous results. 4. Final refit of $σ_{\mathrm{int}}$ to mitigate joint-fit bias. This approach has minimal dependence on the assumed $H_0$ and robustly converges to the injected $ΔH_0$ and dipole sky position. ## 4. Error Modeling: Statistical and Systematic Budget Statistical uncertainty is quantified by the median absolute deviation (MAD) from multiple independent realizations. The chief systematic arises from residuals in $ΔH_0$ recovery dependent on the dipole direction—attributed primarily to imprints from observing cadence. Errors are combined in quadrature: $$ σ_{\mathrm{tot}}^2 = σ_{\mathrm{stat}}^2 + σ_{\mathrm{sys}}^2(ΔH_0) $$ where $σ_{\mathrm{sys}}(ΔH_0)$ is empirically modeled across varying $ΔH_0$. For the clustered ZTF-like configuration with $ΔH_0 = 3$ km s⁻¹ Mpc⁻¹, the achieved precisions are: - $ΔH_0 = 3.05 \pm 0.33$ km s⁻¹ Mpc⁻¹ - Directional precisions: $\pm 3.4^\circ$ (right ascension), $\pm 6.1^\circ$ (declination) With no injected dipole, $σ_{\mathrm{stat}} \approx 0.16$ km s⁻¹ Mpc⁻¹ and $σ_{\mathrm{sys}} \approx 0.44$ km s⁻¹ Mpc⁻¹, giving a combined $σ_{\mathrm{tot}} \approx 0.47$ km s⁻¹ Mpc⁻¹, ensuring no false >1$σ$ signals in null cases. ## 5. Robustness Tests and Sensitivity Analysis Multiple robustness tests validate the stability of the pipeline: - **Monopole independence**: $H_0$ values of 67, 70, 73 km s⁻¹ Mpc⁻¹ produce identical $ΔH_0$ detection and dipole direction. - **Amplitude recovery**: $ΔH_0$ of 1, 2, 3 km s⁻¹ Mpc⁻¹ are unbiasedly recovered; sensitivity, $ΔH_0/σ_{\mathrm{tot}}$, is $\sim$3$σ$ even for $ΔH_0=1$. - **Sky coverage**: Restriction to ZTF’s footprint, MW masking, completeness modeling, and real observing logs do not bias $ΔH_0$ or direction; only cadence imprints induce mild latitude-dependent scatter well within $σ_{\mathrm{tot}}$. - **Large-scale structure**: N-body derived peculiar velocities add scatter but do not bias dipole recovery. - **Volume-limited samples**: A $z<0.06$ subsample ($\sim$1,000 SNe) yields the same dipole detection but with 25–30% larger errors, as expected from reduced statistics. This suite of tests affirms internal consistency, lack of bias from prior parameter choices, and control of major known systematics. ## 6. Prospects and Recommendations for ZTF Dipole Searches For imminent ZTF analyses, the recommended pipeline begins with a bias-controlled, volume-limited SN Ia sample ($0.02