Modified Chaplygin Gas in Cosmology

• Modified Chaplygin Gas (MCG) is an exotic fluid defined by a unique equation of state that interpolates between pressureless matter and a negative-pressure accelerating phase.
• It unifies dark matter and dark energy by recovering ΛCDM and GCG limits through specific parameter choices and influences the universe’s expansion history.
• MCG’s versatility is showcased in its scalar-field and perturbation analyses, though its observational viability depends on careful parameter tuning and model splitting.

Modified Chaplygin Gas (MCG) is an exotic fluid model defined by an equation of state that augments the Chaplygin inverse-density pressure term with a linear barotropic contribution. In its most common form, $p=B\rho-A/\rho^\alpha$, it is studied as a unified dark-sector description that can behave approximately like pressureless matter at early times and like a negative-pressure component at late times, thereby interpolating between matter-dominated and accelerated-expansion regimes (Benaoum, 2012). Across the literature, MCG appears both as a background cosmological model and as an effective fluid in perturbation theory, scalar-field reconstructions, modified-gravity settings, anisotropic cosmologies, viscous models, and even black-hole environments (Xu et al., 2012).

1. Equation of state, notation, and model identity

The defining MCG equation of state is convention-dependent but structurally stable across the literature. The common form is

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$

with model parameters $A$, $B$, and $\alpha$ (Xu et al., 2012). Equivalent notations also appear, such as

$p=\gamma\rho-\frac{B}{\rho^\alpha},$

or

$p=W\rho-\frac{A}{\rho^\alpha},$

with $B>0$ or $A>0$, and typically $0\le \alpha \le 1$ in many analyses (Mazumder et al., 2011).

This equation generalizes related Chaplygin models through standard limiting cases. Setting $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$0 in the $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$1 notation recovers the generalized Chaplygin gas (GCG), while $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$2 in the GCG case reproduces the original Chaplygin gas (Fabris et al., 2010). In several treatments, $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$3CDM is recovered when $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$4 and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$5, and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$6 yields a perfect fluid with constant equation-of-state parameter $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$7 (Xu et al., 2012).

A central feature of the model is its intended interpolation between regimes. At high density, the linear barotropic term dominates, so the fluid can mimic ordinary matter or radiation-like behavior; at low density, the inverse-power term becomes dominant and drives negative pressure (Paul et al., 2014). The model is therefore frequently described as a unified dark matter–dark energy fluid, or a quartessence-like model, rather than as two separately conserved sectors (Fabris et al., 2010).

Form	Parameters	Limit
$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$8	$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$9	Standard MCG form (Xu et al., 2012)
$A$0	$A$1	Same structure in alternate notation (Mazumder et al., 2011)
$A$2	$A$3	GCG when linear term vanishes (Fabris et al., 2010)

The notational variability is not merely cosmetic. In different papers, the symbol $A$4 may denote either the barotropic coefficient or the inverse-density coefficient. This suggests that comparisons across analyses require checking definitions before interpreting quoted parameter values.

2. Background cosmological dynamics

In FRW cosmology, inserting the MCG equation of state into the continuity equation yields the standard density evolution law

$A$5

with

$A$6

or equivalent expressions in redshift $A$7 (Xu et al., 2012). Positivity of the density is commonly enforced through $A$8 (Paul et al., 2013).

The corresponding equation-of-state parameter is

$A$9

which makes the interpolation explicit (Xu et al., 2012). At early times, $B$0 in the parameter regions favored by some cosmological analyses, while at late times the fluid drives acceleration (Xu et al., 2012). Other treatments emphasize a broader asymptotic statement, namely that $B$1 ranges between $B$2 and the barotropic coefficient $B$3 or $B$4, depending on notation (Benaoum, 2012).

In flat FRW cosmology with baryons and radiation, the Hubble function is written as

$B$5

with the actual analysis in one major constraint study restricted to a spatially flat universe (Xu et al., 2012). Related FRW treatments show that MCG can interpolate between radiation-like, dust-like, and $B$6CDM-like behavior, depending on parameter choices and epoch (Mazumder et al., 2011).

Several papers also analyze the model in dynamical-systems form. In one FRW study, the evolution equations are reduced to a two-dimensional autonomous system on the $B$7 phase plane, with a first integral

$B$8

interpretable as motion in a one-dimensional potential (Mazumder et al., 2011). In that analysis, the physically relevant finite critical point is a saddle, while critical points at infinity appear as nodes in the accelerating regime (Mazumder et al., 2011).

A distinct late-time analytical treatment uses a first-order approximation to obtain

$B$9

with explicit expressions for $\alpha$0 and $\alpha$1, and derives a flip time from deceleration to acceleration (Panigrahi et al., 6 Dec 2025). In that treatment, $\alpha$2 and $\alpha$3 are identified as viable ranges, and a non-singular emergent-universe scenario appears when the integration constant is negative (Panigrahi et al., 6 Dec 2025).

3. Perturbations, sound speed, and structure growth

The adiabatic sound speed is a central diagnostic of MCG viability. In unified-fluid analyses it is written as

$\alpha$4

or equivalently

$\alpha$5

depending on notation (Xu et al., 2012). Stability requires $\alpha$6, and some analyses also note that causality ideally favors $\alpha$7 (Xu et al., 2012).

One major MCMC study treats the MCG explicitly as a single unified fluid without splitting it into dark matter and dark energy, evolves perturbations in synchronous gauge with $\alpha$8 and adiabatic initial conditions, and imposes $\alpha$9 numerically as a hard filter during sampling (Xu et al., 2012). That analysis finds the best-fit model very close to $p=\gamma\rho-\frac{B}{\rho^\alpha},$0CDM, with small positive sound speed and early-time cold-dark-matter-like behavior (Xu et al., 2012).

Growth-based studies use the linear perturbation equation

$p=\gamma\rho-\frac{B}{\rho^\alpha},$1

together with the logarithmic growth rate

$p=\gamma\rho-\frac{B}{\rho^\alpha},$2

and the Wang–Steinhardt ansatz

$p=\gamma\rho-\frac{B}{\rho^\alpha},$3

to connect growth observables to the MCG parameters (Paul et al., 2013). In this framework, the zeroth-order growth index for MCG is

$p=\gamma\rho-\frac{B}{\rho^\alpha},$4

with a first-order correction also given in the literature (Paul et al., 2014).

The empirical picture is mixed. One perturbative and matter-power-spectrum study concludes that the hydrodynamical MCG is not a successful cosmic medium unless $p=\gamma\rho-\frac{B}{\rho^\alpha},$5, in which case it reduces to the GCG limit (Fabris et al., 2010). That paper reports acceptable agreement at about $p=\gamma\rho-\frac{B}{\rho^\alpha},$6 only for

$p=\gamma\rho-\frac{B}{\rho^\alpha},$7

typically with either $p=\gamma\rho-\frac{B}{\rho^\alpha},$8 or $p=\gamma\rho-\frac{B}{\rho^\alpha},$9, and interprets this as extreme fine-tuning (Fabris et al., 2010).

By contrast, growth-plus-background analyses report viable best fits with small $p=W\rho-\frac{A}{\rho^\alpha},$0 and $p=W\rho-\frac{A}{\rho^\alpha},$1. One such study finds, for growth + $p=W\rho-\frac{A}{\rho^\alpha},$2 + OHD,

$p=W\rho-\frac{A}{\rho^\alpha},$3

and reports present-day

$p=W\rho-\frac{A}{\rho^\alpha},$4

with $p=W\rho-\frac{A}{\rho^\alpha},$5 small and positive, roughly in the range $p=W\rho-\frac{A}{\rho^\alpha},$6 to $p=W\rho-\frac{A}{\rho^\alpha},$7 (Paul et al., 2013). A broader analysis combining background and growth tests finds

$p=W\rho-\frac{A}{\rho^\alpha},$8

with

$p=W\rho-\frac{A}{\rho^\alpha},$9

for the best-fit MCG model (Paul et al., 2014).

The resulting tension is methodological as much as phenomenological. Hydrodynamical perturbation analyses strongly constrain the barotropic correction, whereas combined background-growth fits can still admit small but nonzero values. This suggests that MCG viability depends sensitively on whether the model is treated as a strict adiabatic fluid, a more general effective component, or an alternative field realization.

4. Observational constraints and statistical status

Observational constraints on MCG span supernovae, BAO, CMB, Hubble data, growth measurements, and more specialized datasets. A widely cited unified-fluid MCMC analysis using WMAP 7-year full CMB temperature and polarization spectra, BAO, and Union2 supernovae with 557 SNe Ia obtains

$B>0$0

$B>0$1

$B>0$2

with

$B>0$3

slightly better than the corresponding $B>0$4CDM fit to the same data, $B>0$5 (Xu et al., 2012). That study interprets the result in terms of dark degeneracy: MCG remains viable, but current data do not strongly distinguish it from $B>0$6CDM (Xu et al., 2012).

An earlier combined fit using 182 Gold SNe Ia, WMAP 3-year, and the SDSS baryon acoustic peak reports best-fit parameters

$B>0$7

with

$B>0$8

and finds that the best-fit effective equation of state crosses $B>0$9 at about $A>0$0, with present-day value

$A>0$1

and $A>0$2 range

$A>0$3

(Lu et al., 2010). The same paper states that MCG has the smallest $A>0$4 among eight models considered, while the Akaike Information Criterion places it in the same support class as the most favored model rather than uniquely selecting it (Lu et al., 2010).

Later background-plus-growth fits generally push the model toward much smaller $A>0$5 and $A>0$6, thereby making it more $A>0$7CDM-like (Paul et al., 2014). This convergence toward small deviations from $A>0$8CDM is a recurring result across independent analyses (Xu et al., 2012).

More recent work continues this pattern while extending the framework. A Pantheon+ plus Cepheid-calibrated analysis of $A>0$9CDM, GCG, MCG, and an altered Chaplygin gas reports for MCG

$0\le \alpha \le 1$0

$0\le \alpha \le 1$1

$0\le \alpha \le 1$2

$0\le \alpha \le 1$3

together with

$0\le \alpha \le 1$4

in that paper’s notation (Yi-Syuan et al., 18 Sep 2025). The same study emphasizes broader posteriors for MCG than for $0\le \alpha \le 1$5CDM or the altered Chaplygin model because of parameter degeneracies (Yi-Syuan et al., 18 Sep 2025).

The statistical status of MCG is therefore not uniform. Some datasets and model-selection criteria place it close to, or slightly ahead of, $0\le \alpha \le 1$6CDM in raw fit quality (Lu et al., 2010), while others emphasize that its best-fit region lies very near the $0\le \alpha \le 1$7CDM limit and is weakened by degeneracy or perturbative tension (Xu et al., 2012).

5. Field-theoretic realizations and theoretical reformulations

MCG has repeatedly been re-expressed in terms of more fundamental degrees of freedom. A review treatment reconstructs the model as a homogeneous minimally coupled scalar field with

$0\le \alpha \le 1$8

and derives a self-interacting potential that reproduces the MCG background evolution exactly (Benaoum, 2012). The same work also studies a tachyonic realization and the mapping between scalar-field and tachyonic-field descriptions (Benaoum, 2012).

A more recent scalar-field construction formulates both GCG and MCG using a canonical Lagrangian

$0\le \alpha \le 1$9

with analytic expressions for $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$00, $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$01, and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$02 (Yi-Syuan et al., 18 Sep 2025). In that framework the MCG equation of state is written as

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$03

and the density evolution is derived as

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$04

(Yi-Syuan et al., 18 Sep 2025). This suggests that the unified-fluid picture can be embedded in canonical scalar dynamics rather than treated only as a phenomenological fluid.

A distinct theoretical route derives MCG from geometrothermodynamics (GTD). In that approach, a fundamental entropy

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$05

induces the Chaplygin-type equation of state, and for $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$06 it reduces to the standard MCG form

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$07

(Benaoum et al., 2019). The GTD scalar curvature is generically nonzero and is interpreted as internal thermodynamic interaction (Benaoum et al., 2019).

MCG has also been embedded in fermionic effective-fluid models. In f-essence cosmology, the pressure is identified with the fermionic Lagrangian $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$08, the energy density with $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$09, and the MCG relation is imposed directly, leading to explicit solvable branches and phantom-crossing behavior for suitable parameter choices (Jamil et al., 2011).

These constructions do not remove the fluid description from the theory; rather, they provide alternative microscopic or effective realizations of the same background equation of state. A plausible implication is that some perturbative objections directed at the purely hydrodynamical MCG need not transfer identically to scalar or spinor realizations.

6. Extensions, modified settings, and nonstandard applications

Beyond standard FRW cosmology, MCG has been studied in several extended settings. In Horava–Lifshitz gravity, observational fits using $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$10, BAO, and the CMB shift parameter show that the effective dark-radiation sector strongly affects the allowed range of the matter-like parameter $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$11; the qualitative result is that “greater dark radiation less the matter contribution in MCG” (Paul et al., 2012). In RS II brane cosmology, Stern $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$12, BAO, and CMB data are used to constrain $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$13 and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$14 for fixed $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$15 and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$16, and the resulting model is described as perfectly consistent with the Union2 sample (Ranjit et al., 2013).

Bulk viscosity and matter creation have also been combined with MCG. In a $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$17-dimensional FRW model with effective pressure

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$18

bulk viscosity is reported to damp structure-growth oscillations and to keep the deceleration parameter negative across the plotted range in the viscous case (Dhankar et al., 25 Apr 2025). In a flat FLRW late-time model with

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$19

the addition of the R22 prior shifts the best-fit Hubble constant upward and improves competitiveness with $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$20CDM in AIC terms (Bhardwaj et al., 25 Sep 2025).

Anisotropic realizations have likewise been explored. In Kantowski–Sachs spacetime, a massless nonlinear spinor field can be chosen so that its effective stress-energy reproduces the MCG equation of state

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$21

and observational constraints yield

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$22

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$23

for the full data combination considered there (Goray et al., 5 May 2026). The present-day values $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$24 and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$25 in that study indicate late-time acceleration together with effective isotropization (Goray et al., 5 May 2026).

MCG has also appeared outside the usual dark-energy context. A magnetogenesis experiment uses a modified Chaplygin gas as one component of a two-fluid plasma system, with equation of state

$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$26

and reports that mixtures containing MCG produce much stronger magnetic fields than a pure plasma fluid, with the strongest fields arising in a $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$27 MCG / $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$28 plasma mixture (Adams et al., 2014). The proposed mechanism is that the Chaplygin component changes the plasma response to gravity, thereby enhancing the Biermann battery source term (Adams et al., 2014).

In black-hole physics, an MCG-like dark fluid has been used to modify the geometry of static charged AdS black holes, affecting geodesics, shadow radius, Hawking temperature, greybody bounds, and quasinormal modes, with constraints drawn from EHT observations of M87* and Sgr A* (Zare et al., 2024). This use is not a cosmological background model in the usual sense, but it extends the Chaplygin-fluid idea into astrophysical spacetime phenomenology.

7. Viability, controversies, and current interpretation

The main conceptual divide in the MCG literature concerns whether the model should be treated as a single unified fluid or decomposed into effective dark-matter and dark-energy pieces. One major analysis explicitly rejects the split, arguing that such decompositions are not unique and can alter perturbation evolution; instead, the entire MCG is treated as one unified dark fluid with its own perturbations (Xu et al., 2012). Other works do perform a split and construct an effective $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$29 for the dark-energy part (Lu et al., 2010).

A second controversy concerns perturbative viability. Background fits alone often make MCG appear observationally acceptable, and in some cases slightly favored in $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$30 relative to $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$31CDM (Lu et al., 2010). However, matter-power-spectrum analyses can be much more restrictive and in one prominent hydrodynamical treatment rule out the modified term unless the model effectively collapses to GCG, $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$32 (Fabris et al., 2010). This is not a contradiction in the narrow sense: the datasets, perturbative assumptions, and model realizations differ. It does indicate that the status of MCG depends strongly on how its perturbations are modeled.

A third recurring theme is parameter degeneracy. Multiple observational analyses find best-fit values of $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$33 and $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$34 very close to zero, implying a phenomenology almost indistinguishable from $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$35CDM at both the background and perturbative levels (Xu et al., 2012). More flexible formulations can fit data well but typically produce broader posterior distributions than simpler models (Yi-Syuan et al., 18 Sep 2025).

The present research picture is therefore internally differentiated rather than uniform. MCG remains a mathematically versatile unified-fluid framework with exact background solutions, scalar and spinor realizations, and successful fits in several observational settings (Benaoum, 2012). At the same time, its hydrodynamical perturbation sector is strongly constrained, and some analyses conclude that only the near-$p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$36CDM or GCG-like corner of parameter space survives precision structure data (Fabris et al., 2010). This suggests that MCG is best understood not as a settled replacement for $p_{MCG}=B\rho_{MCG}-\frac{A}{\rho_{MCG}^{\alpha}},$37CDM, but as a family of closely related unified dark-sector models whose empirical viability depends on the adopted realization, perturbative prescription, and dataset combination.