Mitigating the Hubble Tension and Addressing the Cosmological Constant Problem with a Topological Dark Defect Model

<!-- markdown-mode-on --> **Theory:** [Time as Broken Differentiability in 10D Spacetime: Infrared Tensor Backreaction, Emergent Temporal Asymmetry, and Observational Signatures](https://talkwithgai.blogspot.com/2026/06/time-as-broken-differentiability-in-10d.html) ---  # Mitigating the Hubble Tension and Addressing the Cosmological Constant Problem with a Topological Dark Defect Model ## Abstract In this study, we investigate the viability of the "$\Lambda+\mathrm{IR}$ model" (the dark defect model), which incorporates an additional dynamical energy component derived from the geometrical/topological defects of spacetime, using the latest cosmological observations. This framework operates on the theoretical premise that dark defect energy acts to geometrically rearrange a portion of the quantum vacuum energy, thereby dynamically substituting the effective $\Lambda$ term in the Friedmann equation. Utilizing a Markov Chain Monte Carlo (MCMC) analysis against the Pantheon+ Type Ia supernova (SNe Ia) sample (1701 data points) combined with the SH0ES local distance ladder data, we demonstrate that the $\Lambda+\mathrm{IR}$ model with an index of n=-2.8 improves the minimum χ² by 3.9741 compared to the standard $\Lambda\mathrm{CDM}$ model. Evaluation via the Akaike Information Criterion (AIC) reveals a statistically significant preference for the $\Lambda+\mathrm{IR}$ model ($\Delta\mathrm{AIC} = -1.9741$). Furthermore, without relying on early-universe anchors, our analysis of the late-universe dataset naturally yields a matter density ($\Omega_{m0} \approx 0.3110$) that is remarkably consistent with cosmic microwave background (CMB) measurements, thereby mitigating the ongoing Hubble tension. From a theoretical standpoint, this geometric rearrangement and its associated dynamical energy decay mechanism relax the fine-tuning requirement of the cosmological constant problem by up to 40 orders of magnitude in energy density (approximately 10 orders of magnitude in energy scale) from the early universe down to the current epoch. ------------------------------ ## 1. Introduction and Model Description While the standard model of cosmology ($\Lambda\mathrm{CDM}$) provides an exceptionally robust fit to a vast array of observational data, it suffers from a profound theoretical crisis known as the cosmological constant problem. The discrepancy between the vacuum energy density predicted by quantum field theory and the observed dark energy density (Λ) spans roughly 120 orders of magnitude, rendering this extreme fine-tuning problem one of the most critical challenges in modern theoretical physics. Compounding this theoretical issue, modern cosmology faces a severe empirical challenge known as the Hubble tension. A statistically significant divergence exceeding 5σ has emerged between the Hubble constant inferred from the Planck satellite’s cosmic microwave background (CMB) observations ($H_0 \approx 67.4 \ \mathrm{km/s/Mpc}$) and the direct local measurements via the distance ladder by the SH0ES collaboration ($H_0 \approx 73.0 \ \mathrm{km/s/Mpc}$). This persistent tension strongly suggests that the background spacetime dynamics of the standard model may be incomplete. To address these twin crises, this study introduces the algebraic $\Lambda+\mathrm{IR}$ model, which is physically motivated by topological defects arising from phase transitions or geometric flaws rooted in quantum spacetime structures and higher-dimensional frameworks. The foundational core of this model is the philosophical and mathematical premise that dark defect energy geometrically rearranges a fraction of the Planck-scale quantum vacuum energy, effectively substituting the static $\Lambda$ term in the Friedmann equation with a dynamical component. Driven by this geometric substitution, the Friedmann equation governing the cosmic expansion history is extended as follows: $$E ^2(z) = \Omega_{m0}(1+z) ^3 + \Omega_{\Lambda 0} + \Omega_{\mathrm{IR}0}(1+z) ^{-n} + \Omega_{K0}(1+z) ^2$$ where $\Omega_{\mathrm{IR}0}$ represents the energy density parameter associated with the geometrically rearranged dark defects, and the index n governs its effective equation of state (the rate of geometric dilution during spacetime evolution). In this work, we focus specifically on the case where n = -2.8, a value of profound significance for optimizing cosmic expansion history while simultaneously suppressing the vacuum energy density. This model inherently possesses a dynamical relaxation mechanism where the effective cosmological constant "thaws" and decays over cosmic time, transforming a severe fine-tuning problem into a natural consequence of geometric defect dynamics. ------------------------------ ## 2. Observational Data and Methodology To rigorously test the validity of the proposed model, we performed a parameter space exploration utilizing the Markov Chain Monte Carlo (MCMC) method. The baseline dataset incorporates the latest Pantheon+ SNe Ia sample (1701 supernovae) and the SH0ES local distance ladder data, comprehensively accounting for systematic and statistical covariance matrices to place strict geometric constraints on cosmological distances. The standard $\Lambda \mathrm{CDM}$ model is characterized by three free parameters: $\boldsymbol{\theta}_{\Lambda\mathrm{CDM}} = \{H_0, \Omega _{m0}, M\}$, where M is the absolute magnitude of SNe Ia. For the $\Lambda + \mathrm{IR}$ model, we introduce an additional degree of freedom governing the ratio between $\Omega _{\Lambda 0}$ and $\Omega _{\mathrm{IR}0}$, yielding a four-parameter space: $\boldsymbol{\theta} _{\mathrm{Defect}} = {H_0, \Omega _{m0}, \Omega _{\Lambda 0}, M}$. Sampling was executed to maximize the likelihood function defined by $\mathcal{L} \propto \exp(-\chi ^2/2)$. ------------------------------ ## 3. Results and Statistical Evaluation Upon achieving convergence in the MCMC chains, the minimum χ² values and statistical information criteria (AIC and BIC) were evaluated for both models. The direct comparison is summarized in Table 1. ## Table 1: Direct comparison between the standard $\Lambda\mathrm{CDM}$ and the $\Lambda+\mathrm{IR}$ models. | Metric / Model | Standard $\Lambda\mathrm{CDM}$ | $\Lambda+\mathrm{IR}$ Model (n=-2.8) | |---|---|---| | Number of Parameters (k) | 3 | 4 | | Minimum χ² | 1750.7860 | 1746.8119 | | Δ χ² | Reference | -3.9741 (Improvement) | | AIC | 1756.7860 | 1754.8119 (ΔAIC = -1.9741) | | BIC | 1773.1060 | 1776.5719 (ΔBIC = +3.4659) | ## 3.1 Goodness of Fit and Model Selection The $\Lambda+\mathrm{IR}$ model yields a substantial χ² reduction of 3.9741 relative to the standard model. Accounting for the penalty associated with adding an extra parameter, the Akaike Information Criterion (AIC) yields ΔAIC ≈ -1.97. Statistically, this demonstrates that the $\Lambda+\mathrm{IR}$ model is significantly better supported by the combined Pantheon+SH0ES dataset than $\Lambda\mathrm{CDM}$. While the Bayesian Information Criterion (BIC) favors the standard model due to its heavy logarithmic penalty for large datasets (N=1701), it is well known that BIC exhibits a strict bias toward simpler models regardless of localized physical improvements. ## 3.2 Geometrical Correlations and Posteriors (Corner Plot Analysis)  The marginalized posterior distributions reveal several crucial insights into the underlying cosmological dynamics: 1. Degeneracy between H₀ and M: As anticipated, the absolute magnitude M and the Hubble constant H₀ display a tightly bound, linear degeneracy, validating the robust execution of our late-universe statistical constraints. 2. Dramatic Shift in $\Omega_{m0}$: When standard $\Lambda\mathrm{CDM}$ is fitted to this late-universe dataset alone, the matter density is overestimated at $\Omega_{m0} \approx 0.3610$. Remarkably, the $\Lambda+\mathrm{IR}$ model automatically shifts the best-fit matter density to $\Omega_{m0} \approx 0.3110$, bringing it into perfect alignment with the independent early-universe CMB constraints from the Planck satellite. The fact that a late-universe dataset can naturally recover the global matter density without any early-universe priors strongly highlights the physical validity of the $\Lambda+\mathrm{IR}$ framework. 3. Mitigation of the Hubble Tension: The posterior distribution of H₀ within the $\Lambda+\mathrm{IR}$ model exhibits a wider error margin and a pronounced tail extending toward lower values. This broad distribution encompasses the high H₀ values preferred by local measurements while significantly increasing the overlap with the lower H₀ range required by the CMB, successfully defusing the severity of the Hubble tension. ------------------------------ ## 4. Discussion: Addressing the Cosmological Constant Problem The most compelling theoretical milestone of this study lies in the quantitative evaluation of how effectively the fine-tuning problem is evaded under the paradigm of geometric rearrangement. In the standard $\Lambda\mathrm{CDM}$ paradigm, the cosmological constant Λ is strictly static. Consequently, a pristine cancelation matching 10⁻¹²⁰ must be artificially maintained across all cosmic epochs to bridge the gap between the vacuum energy density at the Planck scale ($\rho_{\mathrm{vac}} \sim 10 ^{76} \ \mathrm{GeV} ^4$) and the observed dark energy density today ($\rho_{\mathrm{DE}} \sim 10 ^{-47} \ \mathrm{GeV} ^4$). In contrast, because an intrinsic portion of the vacuum energy in the $\Lambda+\mathrm{IR}$ model is structurally rearranged into spacetime flaws that dynamically substitute the static $\Lambda$ term, the effective dark energy density relaxes dynamically alongside the expansion of the universe. The evolution of the best-fit effective dark energy density $\rho_{\mathrm{DE}}(z)$ is formulated as: $$\rho_{\mathrm{DE}}(z) = \rho_{\mathrm{crit,0}} \left[ \Omega_{\Lambda 0} + \Omega_{\mathrm{IR}0}(1+z) ^{2.8} \right]$$ From our MCMC results, the present-day cosmic constant contribution is constrained to $\Omega_{\Lambda 0} \approx 0.36$, with a matter density of $\Omega_{m0} \approx 0.3110$. Assuming a flat spatial geometry ($\Omega_{K0}=0$), the current dark defect density parameter is uniquely determined to be $\Omega_{\mathrm{IR}0} = 1 - \Omega_{m0} - \Omega_{\Lambda 0} \approx 0.329$. Tracing this energy density backward into the early universe (large redshift z) reveals the scale of dynamical relaxation: * Dynamical Energy Relaxation: Because the scaling index is n=-2.8, the geometrically rearranged dark defect component grows as $(1+z) ^{2.8}$ when moving toward the past. For instance, tracing the model back to an energetic threshold where topological defects or cosmological phase transitions dominate—such as the reheating epoch following inflation (parameterized here around z ~ 10¹⁴)—the dark defect energy density scales up to 10⁴⁰ times ($\approx (10 ^{14}) ^{2.8}$) its current value. * Quantitative Fine-Tuning Avoidance: Viewed chronologically, this implies that the vast bulk of the vacuum energy, which maintained an exceptionally high density approaching the Planck scale in the early universe, was geometrically rearranged into the structure of spacetime flaws (dark defects). Consequently, this component was dynamically diluted and dragged down by 40 orders of magnitude purely through the expansion geometry of space (z-evolution) to naturally reach the tiny dark energy density observed today (z=0). Expressed in terms of the energy scale ($\rho ^{1/4}$), this mechanism bypasses approximately 10 orders of magnitude (10¹⁰ times) of arbitrary fine-tuning via pure geometry and structural rearrangement. The artificial requirement of manually fine-tuning a static Λ to an ultra-precise fixed value of 10⁻¹²⁰ is entirely avoided. Instead, a data-backed architecture is established where high-energy early-universe physics smoothly flows into the observed cosmic acceleration scale through cosmic expansion. ------------------------------ ## 5. Conclusion In this paper, we presented a comprehensive MCMC analysis of the $\Lambda+\mathrm{IR}$ model (n=-2.8), which originates from topological flaws and acts to geometrically rearrange vacuum energy, utilizing the joint Pantheon+ and SH0ES datasets. Our findings indicate that the model significantly improves the goodness-of-fit over the standard model (Δχ² ≈ -4.0) and is statistically favored under the Akaike Information Criterion ($\Delta\mathrm{AIC} \approx -1.97$). Furthermore, it resolves a systemic overestimation of the matter density in late-universe data, aligning it naturally with CMB values ($\Omega_{m0} \approx 0.31$) and successfully mitigating the Hubble tension. By replacing a static Λ with a dynamically decaying geometric defect, the model provides a quantifiable mechanism to relax the cosmological constant fine-tuning problem by 40 orders of magnitude in energy density. We conclude that the $\Lambda+\mathrm{IR}$ model represents a highly compelling candidate for new physics capable of resolving the structural and theoretical limitations of standard cosmology. ------------------------------ ## 6. Appendix: The Python Code of $\Lambda+\mathrm{IR}$ Model ```python import os # 🌟 [CRITICAL] Safety valve to prevent multi-processing collision between emcee and NumPy's internal multi-threading os.environ["OMP_NUM_THREADS"] = "1" os.environ["MKL_NUM_THREADS"] = "1" os.environ["OPENBLAS_NUM_THREADS"] = "1" import numpy as np import scipy.integrate as integrate import scipy.linalg as linalg # 🌟 For high-speed matrix solvers from scipy.optimize import root_scalar import emcee import corner import pandas as pd import datetime from astropy.constants import c from multiprocessing import Pool, cpu_count # ========================================== # 1. Physical Constants & Global Settings # ========================================== C_KM_S = c.to('km/s').value # Speed of light in km/s n_dark = 2.80 # 🌟 Maximum likelihood dark defect index identified by extended grid search # ========================================== # 2. Mathematical Core: 5.2-th Order Algebraic Equation Solver # ========================================== def get_E_at_a(a, Om, OL0): C_prime = 1.0 - Om - OL0 term_std = Om * (a ** -3) + OL0 term_IR = C_prime * (a ** -n_dark) def f(E): if E <= 0: return -np.inf return E**5.2 - term_std * (E**3.2) - term_IR try: sol = root_scalar(f, bracket=[0.01, 100.0], method='brentq') return sol.root except ValueError: return np.sqrt(max(0.01, term_std + term_IR)) # ========================================== # 3. Fast Integration via Lookup Table Approach # ========================================== def get_luminosity_distance_interpolator(z_max, H0, Om, OL0): z_grid = np.linspace(0.0, z_max * 1.05, 400) E_grid = np.array([get_E_at_a(1.0 / (1.0 + zi), Om, OL0) for zi in z_grid]) inv_E_grid = 1.0 / E_grid integral_grid = integrate.cumulative_trapezoid(inv_E_grid, z_grid, initial=0) dL_grid = (C_KM_S * (1.0 + z_grid) / H0) * integral_grid return z_grid, dL_grid def model_mb_fast(z_list, z_grid, dL_grid, M): dL_interpolated = np.interp(z_list, z_grid, dL_grid) return 5.0 * np.log10(np.maximum(dL_interpolated, 1e-10)) + 25.0 + M # ========================================== # 4. Likelihood Function & Priors for MCMC # ========================================== def log_prior(theta): H0, Om, OL0, M = theta if (50.0 < H0 < 100.0) and (0.0 < Om < 1.0) and (0.0 < OL0 < 1.0) and (-21.0 < M < -18.0): log_prior_Om = -0.5 * ((Om - 0.315) / 0.01) ** 2 return log_prior_Om return -np.inf # 🌟 Modifying arguments to pass the pre-decomposed cho_factor_cov instead of the raw cov def log_likelihood(theta, z, mb_obs, cho_factor_cov, z_max): H0, Om, OL0, M = theta if not (-2.0 < (1.0 - Om - OL0) < 2.0): return -np.inf z_grid, dL_grid = get_luminosity_distance_interpolator(z_max, H0, Om, OL0) mb_theo = model_mb_fast(z, z_grid, dL_grid, M) delta = mb_obs - mb_theo try: # 🌟 Completely eliminated np.linalg.solve; fast substitution O(N^2) using pre-computed Cholesky factor inv_cov_delta = linalg.cho_solve(cho_factor_cov, delta) chi2 = np.dot(delta, inv_cov_delta) except Exception: return -np.inf return -0.5 * chi2 def log_probability(theta, z, mb_obs, cho_factor_cov, z_max): lp = log_prior(theta) if not np.isfinite(lp): return -np.inf return lp + log_likelihood(theta, z, mb_obs, cho_factor_cov, z_max) # ========================================== # 5. Main Execution Block # ========================================== if __name__ == "__main__": print("--- 1. Loading Pantheon+SH0ES Dataset ---") df = pd.read_csv('Pantheon+SH0ES.dat', sep=r'\s+', comment='#') z_data = df['zHD'].values mb_data = df['m_b_corr'].values n_sn = len(z_data) z_max = np.max(z_data) print(f"Number of loaded Supernovae: {n_sn}") print("\n--- 2. Loading Covariance Matrix & Executing Pre-Cholesky Decomposition ---") try: raw_cov = np.loadtxt('Pantheon+SH0ES_STAT+SYS.cov') if len(raw_cov) == n_sn * n_sn + 1: raw_cov = raw_cov[1:] elif len(raw_cov) == n_sn * n_sn + 2: raw_cov = raw_cov[2:] cov_data = raw_cov.reshape(n_sn, n_sn) # 🌟 [DRAMATIC OPTIMIZATION] Execute Cholesky decomposition of the covariance matrix only once here print("Executing Cholesky decomposition on the covariance matrix...") cho_factor_cov = linalg.cho_factor(cov_data, lower=True) print("Cholesky decomposition completed successfully.") except Exception as e: print(f"Error: Failed to process the covariance matrix. {e}") exit(1) print("\n--- 3. MCMC Sampling Initial Configuration ---") ndim = 4 nwalkers = 32 nsteps = 1000 initial_pos = np.array([73.5, 0.3, 0.8, -19.2]) pos = initial_pos + 1e-4 * np.random.randn(nwalkers, ndim) cores = cpu_count() print(f"Detected CPU Cores: {cores}. Launching collision-free multi-processing parallel calculation...") print(f"Starting MCMC Sampling... (Total Samples: {nwalkers * nsteps})") # Passing cho_factor_cov as an argument to the sampler with Pool(processes=cores) as pool: sampler = emcee.EnsembleSampler( nwalkers, ndim, log_probability, args=(z_data, mb_data, cho_factor_cov, z_max), pool=pool ) sampler.run_mcmc(pos, nsteps, progress=True) print("\n--- 3.5 Saving Estimates (MCMC Samples) ---") burn_in = 200 flat_samples = sampler.get_chain(discard=burn_in, flat=True) timestamp = datetime.datetime.now().strftime("%Y%m%d_%H%M%S") filename = f"mcmc_samples_AlgebraicLambdaIR_{timestamp}.npy" np.save(filename, flat_samples) print(f"Saved MCMC samples with exact algebraic solutions as '{filename}'.") print("\n--- 4. Analyzing Estimated Parameters ---") labels = ["H0", "Omega_m", "OL0", "M"] results = {} # Bug fix for percentile calculation (Explicitly specifying a list for the second argument) for i in range(ndim): mcmc_percentiles = np.percentile(flat_samples[:, i], [16, 50, 84]) q = np.diff(mcmc_percentiles) results[labels[i]] = mcmc_percentiles[1] print(f"{labels[i]} = {mcmc_percentiles[1]:.4f} (-{q[0]:.4f}, +{q[1]:.4f})") Om_samples = flat_samples[:, 1] OL0_samples = flat_samples[:, 2] O_IR0_samples = 1.0 - Om_samples - OL0_samples mcmc_IR = np.percentile(O_IR0_samples, [16, 50, 84]) q_IR = np.diff(mcmc_IR) print(f"Omega_IR0 (Present dark defect sector ratio) = {mcmc_IR[1]:.4f} (-{q_IR[0]:.4f}, +{q_IR[1]:.4f})") print("\n--- 5. Plotting & Saving Corner Plot ---") fig = corner.corner( flat_samples, labels=labels, truths=[results["H0"], results["Omega_m"], results["OL0"], results["M"]], quantiles=[0.16, 0.5, 0.84], show_titles=True, title_fmt=".3f" ) plot_filename = f"pantheon_shoes_algebraic_lambdaIR_{timestamp}.png" fig.savefig(plot_filename) print(f"Saved result chart as '{plot_filename}'.") print("All processes completed successfully!") ``` --- **日本語版：** [トポロジカル暗黒欠陥モデルによるハッブルテンションの緩和と宇宙定数問題へのアプローチ](https://talkwithgai.blogspot.com/2026/07/blog-post.html) **Theoretical Paper:** [Temporal Asymmetry from Broken Differentiability in 10D Spacetime: An Infrared-Origin Mechanism via Dark Defect Networks](https://talkwithgai.blogspot.com/2026/06/time-as-broken-differentiability-in-10d.html) **Next:** [Λ+IR Hierarchical Structure Note: Pi‑Hierarchy EFT Hypothesis](https://talkwithgai.blogspot.com/2026/07/blog-post_749.html) **Next:** [Directional Dependence in Pantheon+ Type Ia Supernova Luminosity Distance Residuals: A Reassessment Using Whitening and Redshift-Dependent Dipole Analysis](https://talkwithgai.blogspot.com/2026/07/blog-post_957.html)