Dark Models in Cosmology

• Dark models are theoretical frameworks that extend the standard cosmological model by incorporating interactions and modifications to explain dark matter and dark energy.
• They employ interacting dark energy/dark matter models, modified gravity, and particle physics constructions to address cosmic acceleration, the coincidence problem, and the Hubble tension.
• Observational constraints from SNIa, BAO, and Hubble datasets reveal varying coupling parameters, offering insights into late-time universe behavior and structure formation.

Dark Models

A “dark model” is a theoretical or phenomenological framework posited to explain the nature, composition, and dynamics of the universe’s dark sector—namely dark matter (DM) and dark energy (DE)—by extending the Standard Model of particle physics and/or classical general relativity. The class of dark models is broad, encompassing interacting dark energy/dark matter models, modified gravity, non-minimal couplings, and particle physics constructions for dark matter, each designed to address observational puzzles such as the cosmic acceleration, the coincidence problem, or the microphysical nature of DM.

1. Interacting Dark Energy and Dark Matter Models

Interacting dark sector models generalize the standard cosmological scenario by allowing non-gravitational energy exchange between DE and DM. The fundamental premise is that the total energy-momentum tensor remains covariantly conserved, but individual DM and DE stress-energy tensors may exchange energy via a current $Q^\mu$: $$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$ In homogeneous, isotropic cosmology, this reduces to coupled continuity equations: $$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$ where $Q$ encodes the interaction rate.

Phenomenological choices for $Q$ include:

• $Q = 3\beta H\rho_m$ (“β-model”)
• $Q = 3\eta H\rho_{de}$ (“η-model”)
• $Q = \gamma H(\rho_m + \rho_{de})$

The sign of $\beta,\eta$ sets the direction of energy flow: positive values correspond to DE$\to$DM transfer, while negative values correspond to DM$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$0DE (Rugg et al., 2024, Zimdahl, 2012).

Dynamical systems formulations extend these models, including non-minimal couplings (e.g., $$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$1 for a scalar DE (Ashmita et al., 2024)), leading to sectors with rich critical point structure and diverse late-time behaviors.

2. Observational Constraints and Phenomenology

Analyses employing Type Ia supernovae (SNIa) luminosity distances, BAO, and Hubble expansion datasets constrain interacting dark models. MCMC fits with flat priors on coupling parameters ($$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$2 or $$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$3) and cosmological parameters yield the following:

Model/Data	$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$4	Coupling	$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$5 [km/s/Mpc]
$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$6CDM / Pantheon	0.280(10)	--	71.85(22)
$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$7CDM / OHD	0.250(20)	--	70.79(123)
β-model / Pantheon	0.171(54)	$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$8	72.16(28)
β-model / OHD	0.555(128)	$$\nabla_\mu T^{\mu\nu}_m = Q^\nu, \quad \nabla_\mu T^{\mu\nu}_{de} = -Q^\nu$$9	65.77(262)
η-model / Pantheon	0.186(55)	$$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$0	72.37(36)
η-model / OHD	0.445(115)	$$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$1	65.81(322)

Parentheses indicate $$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$2 uncertainties. Pantheon data prefer negative couplings (DM$$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$3DE), low $$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$4–$$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$5, and high $$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$6 km/s/Mpc. OHD/BAO data favor positive couplings (DE$$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$7DM), high $$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$8–$$\dot\rho_m + 3H\rho_m = Q \qquad \dot\rho_{de} + 3H(\rho_{de}+p_{de}) = -Q$$9, and low $Q$0 km/s/Mpc. The tension between datasets limits the statistical significance of $Q$1 at $Q$2. Combined likelihoods mildly prefer DE$Q$3DM (Rugg et al., 2024).

Observational consequences of $Q$4 include measurable differences in the expansion history $Q$5, the growth rate of cosmological structures, and the effective equation of state parameter $Q$6, which can vary significantly from the native $Q$7 of the uncoupled DE model (Avelino et al., 2012, Zimdahl, 2012). Large-scale surveys (Euclid, DESI, LSST, SKA, CMB-S4) are projected to push coupling constraints to sub-percent levels (Zimdahl, 2012).

3. Addressing the Coincidence and $Q$8 Problems

A principle motivation for dark models with couplings is the coincidence problem: the near-equality of DM and DE energy densities today, which in $Q$9CDM is an unexplained temporal coincidence. In the β- or η-models, the evolution of the ratio $Q$0 can be substantially slowed. For $Q$1, $Q$2 evolves slowly, maintaining $Q$3 values over longer epochs—alleviating fine-tuning (Rugg et al., 2024).

Moreover, the impact of couplings on late-time expansion allows dark models to partially relax the Hubble tension by yielding a dynamically evolving, rather than fixed, DE density. For specified $Q$4, the background evolution can mimic that of phantom ($Q$5) models, even when $Q$6, thereby producing late-time $Q$7 values closer to local vs. early-universe inferences (Rugg et al., 2024, Avelino et al., 2012).

4. Theoretical Frameworks and Model Construction

Model-building approaches for dark models span:

• Phenomenological continuity equations: Effective ansätze for $Q$8 without explicit microphysics, used in current data-driven analyses and dynamical system studies (Rugg et al., 2024, Ashmita et al., 2024).
• Field-theoretic realizations: Scalar or tachyonic field Lagrangians with explicit DE–DM coupling (e.g., Yukawa interactions $Q$9 for a DE scalar $Q = 3\beta H\rho_m$0 and DM fermion $Q = 3\beta H\rho_m$1) (Micheletti, 2010, Pourtsidou et al., 2013, Zimdahl, 2012).
• Modified gravity: Generalizations such as $Q = 3\beta H\rho_m$2 or $Q = 3\beta H\rho_m$3 gravity incorporate effective dark-sector interactions at the metric/action level, yielding “effective DE” (Astashenok, 2013, Yoo et al., 2012).
• Holographic models: Enforcement of holographic DE constraints on the potential energy, combined with dark-sector coupling in Lagrangian field theory (Micheletti, 2010).

Entries in this taxonomy (see table below) capture the key mathematical structure:

Model Type	Interaction Term $Q = 3\beta H\rho_m$4	Reference
β-model	$Q = 3\beta H\rho_m$5	(Rugg et al., 2024, Zimdahl, 2012)
η-model	$Q = 3\beta H\rho_m$6	(Rugg et al., 2024, Zimdahl, 2012)
Two-param ansatz	$Q = 3\beta H\rho_m$7	(Avelino et al., 2012)
Scalar field coupling	$Q = 3\beta H\rho_m$8	(Ashmita et al., 2024, Pourtsidou et al., 2013)
Holographic Yukawa	$Q = 3\beta H\rho_m$9	(Micheletti, 2010)
Decaying vacuum	$Q = 3\eta H\rho_{de}$0	(Zimdahl, 2012)

These frameworks can support both analytic solution for background evolution and full Boltzmann integration for CMB and LSS predictions.

5. Implications for Structure Formation and High-Precision Probes

Dark models directly affect cosmic expansion and structure-growth histories. Key observable signatures include:

• Modifications to the matter/DE density evolution and $Q = 3\eta H\rho_{de}$1, affecting SN distances, BAO positions, and SNIa statefinder diagnostics (e.g., jerk parameter) (Zimdahl, 2012).
• Altered growth rates for linear perturbations, potentially detectable in redshift-space distortions and weak lensing, especially for $Q = 3\eta H\rho_{de}$2 (Zimdahl, 2012, Yoo et al., 2012).
• Changes to the CMB angular power spectrum, notably the integrated Sachs–Wolfe effect and clustering on large scales if DE perturbations or non-adiabatic pressure contributions are enhanced by the interaction (Zimdahl, 2012).
• Non-trivial degeneracies between interacting DE and phantom models: IDE models with $Q = 3\eta H\rho_{de}$3 can exactly mimic the background and CMB signatures of a $Q = 3\eta H\rho_{de}$4 scenario if CMB-inferred $Q = 3\eta H\rho_{de}$5 is shifted appropriately (Avelino et al., 2012).
• Implications for neutrino phenomenology: dark models coupling DE to neutrinos can alter redshift-dependent neutrino oscillation probabilities, potentially distinguishable in next-generation high-z neutrino telescopes (Khalifeh et al., 2021).

6. Extensions Beyond Interacting Dark Sector Models

Beyond classical coupled dark sector models, the “dark model” umbrella encompasses:

• Modified gravity theories ($Q = 3\eta H\rho_{de}$6, $Q = 3\eta H\rho_{de}$7): provide effective dark components through generalized Einstein–Hilbert actions, yielding rich cosmological dynamics and often evading standard no-go theorems for acceleration (Astashenok, 2013, Yoo et al., 2012).
• Non-minimal dark sector constructions: Composite dark matter, accidental symmetry models, secluded dark sectors with their own gauge dynamics, and Stueckelberg or hidden photon extensions (Fortes et al., 2017, Palmisano et al., 2024). These model the particle physics of DM, including late-time decay, self-interactions, and baryonic portal phenomenology.
• Cosmological inhomogeneity (e.g., Lemaitre–Tolman–Bondi): attempts to mimic cosmic acceleration without DE by positing large-scale inhomogeneities (Yoo et al., 2012).
• Machine learning frameworks: Model selection employing VAE–GAN hybrids for reconstructing SN distance moduli and discriminating dark energy models from data (Li et al., 2019).

7. Theoretical and Experimental Status

Current cosmological observations provide strong constraints on allowed interactions in dark models, typically requiring coupling strengths $Q = 3\eta H\rho_{de}$8 for $Q = 3\eta H\rho_{de}$9–type scenarios (Rugg et al., 2024, Zimdahl, 2012). Standard $Q = \gamma H(\rho_m + \rho_{de})$0CDM remains observationally favored, but interacting models, particularly those aimed at addressing the coincidence or $Q = \gamma H(\rho_m + \rho_{de})$1 tension, remain viable in limited parameter ranges. Future experiments, especially those with sensitivity to detailed structure growth, expansion history, or novel dark sector signatures (e.g., fifth force, neutrino oscillations), will further constrain or discover departures from $Q = \gamma H(\rho_m + \rho_{de})$2CDM.

Open challenges in dark model research include unambiguous observational discrimination of small coupling strengths, robust UV completions embedding the phenomenological interactions, and joint parameter estimation across extended cosmological datasets (Rugg et al., 2024, Zimdahl, 2012, Yoo et al., 2012).