Modified Generalized Chaplygin Gas (MGCG)

• Modified Generalized Chaplygin Gas (MGCG) is a cosmological fluid model defined by p = Aρ – B/ρ^α that transitions from dust-like to dark-energy regimes.
• It unifies dark matter and dark energy by employing thermodynamic treatments, scalar-field reconstructions, and perturbative analyses within Friedmann cosmology.
• Observational constraints and varying parameter choices in MGCG offer practical insights into cosmic acceleration and the evolution of large-scale structure.

The Modified Generalized Chaplygin Gas (MGCG) is a phenomenological cosmological fluid, widely written in the form $p=A\rho-B\rho^{-\alpha}$, that has been studied as a dark-energy component and as a unified dark-matter–dark-energy sector in Friedmann cosmology. Across the literature, closely related symbol conventions also appear, including $p=B\rho-A/\rho^\alpha$ and $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$; a distinct “revisited” variant further replaces the constant Chaplygin term by a redshift-dependent function $\vartheta(z)$ (Ebadi et al., 2015, Benaoum, 2012, Bouhmadi-López et al., 2015, Deng, 2011). In all of these formulations, the central idea is an interpolation between an early dust-like or barotropic regime and a late negative-pressure regime, together with a substantial literature on horizon thermodynamics, scalar-field realizations, perturbative viability, and observational constraints.

1. Definitions, notation, and limiting cases

In the notation used in the thermodynamic treatment of MGCG, the equation of state is

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$

with $B>0$, where $\alpha$ controls the interpolation between dust-like and cosmological-constant-like behavior (Ebadi et al., 2015). The same functional structure is written in other papers as

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

or, for the late-time perturbation analysis,

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$

so the literature is not uniform in its symbol assignments (Fabris et al., 2010, Bouhmadi-López et al., 2015). A further unified-dark-sector construction uses

$$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$

with $p=B\rho-A/\rho^\alpha$0 chosen so that the fluid behaves as pressureless matter at high redshift and as quiessence dark energy at low redshift (Deng, 2011).

Within the standard $p=B\rho-A/\rho^\alpha$1 parameterization, the most common special cases are the following.

Parameter choice	Recovered model	Equation of state
$p=B\rho-A/\rho^\alpha$2	Generalized Chaplygin Gas	$p=B\rho-A/\rho^\alpha$3
$p=B\rho-A/\rho^\alpha$4	Original Chaplygin Gas	$p=B\rho-A/\rho^\alpha$5
$p=B\rho-A/\rho^\alpha$6	Perfect fluid	$p=B\rho-A/\rho^\alpha$7

The corresponding equation-of-state parameter is

$p=B\rho-A/\rho^\alpha$8

so the barotropic term fixes the early-time limit while the Chaplygin term dominates at low density (Jamil et al., 2011). In this form, phantom crossing $p=B\rho-A/\rho^\alpha$9 occurs when

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$0

a condition that is explicit in the f-essence treatment of the same equation of state (Jamil et al., 2011).

2. FRW evolution and effective cosmic phases

For a spatially flat FRW universe without interaction, the continuity equation

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$1

integrates exactly for the MGCG equation of state. In the $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$2 convention one obtains

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$3

and hence

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$4

with $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$5 an integration constant (Ebadi et al., 2015). Equivalent expressions appear in the FRW and scalar-field treatments of the model, sometimes after introducing the dimensionless combination

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$6

depending on notation (Paul et al., 2014, Bhardwaj et al., 25 Sep 2025).

This solution makes the interpolation mechanism explicit. In the f-essence analysis, $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$7 at early times and $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$8 at late times, so the model moves from a barotropic or approximately dust-like phase to an asymptotic cosmological-constant regime (Jamil et al., 2011). In the mGCG late-time perturbation study, writing $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)A/\varepsilon_{\rm Ch}^\alpha$$9, one finds that $\vartheta(z)$0 for $\vartheta(z)$1, again mimicking a cosmological constant, जबकि for $\vartheta(z)$2 the fluid approaches a barotropic phase with $\vartheta(z)$3 (Bouhmadi-López et al., 2015). In the FRW scalar-field summary, the deceleration parameter tends from $\vartheta(z)$4 at early times to $\vartheta(z)$5 at late times (Benaoum, 2012).

A related late-time-acceleration construction starts not from the usual equation of state but from the ansatz

$\vartheta(z)$6

which yields radiation at early times, dust at intermediate times, an MGCG-like regime $\vartheta(z)$7, and finally a de Sitter phase (Abdussattar et al., 2016). This suggests that MGCG behavior can emerge as an effective regime even when the fundamental starting point is a more general density ansatz.

3. Thermodynamics on the apparent horizon

A major line of work interprets MGCG thermodynamically on the apparent horizon of FRW spacetime. For $\vartheta(z)$8, the apparent-horizon radius is $\vartheta(z)$9, and in equilibrium the fluid temperature equals the horizon temperature,

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$0

Using Gibbs’ law,

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$1

the noninteracting entropy variation becomes

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$2

for the standard MGCG equation of state (Ebadi et al., 2015).

When MGCG exchanges energy with pressureless dark matter through

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$3

the continuity equations are modified and the entropy balance acquires an additive interaction correction,

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$4

In the same framework, first-order thermal fluctuations generate the logarithmic correction

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$5

where $$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$6 is the heat capacity at constant volume (Ebadi et al., 2015). The paper explicitly ties the strength of the mutual interaction $$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$7 to the size of this fluctuation correction.

The same thermodynamic analysis also addresses the coincidence problem. If dark energy is required to decay into dark matter, then $$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$8, so

$$p_D=A\,\rho_D-\frac{B}{\rho_D^\alpha}, \qquad 0\le \alpha\le 1,\qquad A\neq -1,$$9

Writing $B>0$0, this becomes

$B>0$1

The resulting bound constrains admissible interaction couplings if one attempts to use the interaction to alleviate coincidence (Ebadi et al., 2015).

A different thermodynamic route derives the modified Chaplygin gas from geometrothermodynamics (GTD). In that construction the entropy is taken as

$B>0$2

and for $B>0$3 one recovers the standard Modified Chaplygin Gas form $B>0$4 (Benaoum et al., 2019). The equilibrium-manifold scalar curvature is nonzero in general, and in GTD this is interpreted as a signal of internal thermodynamic interaction.

4. Scalar-field, tachyonic, and f-essence representations

The MGCG equation of state admits several field-theoretic realizations. In the homogeneous canonical-scalar description,

$B>0$5

and matching to MGCG gives

$B>0$6

A closed-form self-interacting potential can be written in terms of $B>0$7, making explicit the interpolation between an early-time regime and an asymptotic constant-potential regime (Benaoum, 2012).

A more systematic scalar-field reconstruction writes the equation of state as

$B>0$8

with the identifications $B>0$9, $\alpha$0, and $\alpha$1. The corresponding potential is

$\alpha$2

where

$\alpha$3

Using bounded variables

$\alpha$4

the exact modified Chaplygin gas perfect-fluid solution appears as the straight invariant line

$\alpha$5

in the phase plane (Uggla, 2013).

The dynamical-systems analysis yields a sharp statement about exactness versus approximation. No other solutions stay close to the perfect-fluid orbit during their entire temporal evolution, but there exists an open subset of solutions that stay arbitrarily close during an intermediate time interval, and into the future when the scalar-field potential has a global minimum (Uggla, 2013). In the same study, the future asymptotics depend on the sign of

$\alpha$6

For $\alpha$7, the de Sitter point $\alpha$8 is a unique global attractor; for $\alpha$9, the symmetric de Sitter points $p=B\rho-\frac{A}{\rho^\alpha}$0 are the sinks (Uggla, 2013).

Beyond canonical scalars, the same equation of state has been embedded in f-essence cosmology through the action

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

with

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

and explicit positive- and negative-pressure solutions for suitable free parameters (Jamil et al., 2011). A tachyonic-field mapping has also been developed, yielding an exact hypergeometric relation between the tachyon $p=B\rho-\frac{A}{\rho^\alpha}$3 and the scalar-field variable, together with a quadratic slow-roll tachyon potential $p=B\rho-\frac{A}{\rho^\alpha}$4 (Benaoum, 2012).

5. Perturbations, clustering, and the viability debate

The perturbative status of MGCG is one of the most disputed aspects of the model. In a hydrodynamical unified-dark-sector treatment with equation of state $p=B\rho-\frac{A}{\rho^\alpha}$5, synchronous-gauge perturbations and the 2dFGRS matter power spectrum lead to a very strong suppression of the barotropic term: the best fit drives $p=B\rho-\frac{A}{\rho^\alpha}$6, and at $p=B\rho-\frac{A}{\rho^\alpha}$7

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

The same analysis argues that even tiny nonzero $p=B\rho-\frac{A}{\rho^\alpha}$9 raises the adiabatic sound speed enough to generate acoustic oscillations or excessive suppression of sub-horizon density fluctuations, so the model is “not a successful candidate for the cosmic medium unless $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$0” (Fabris et al., 2010).

A different late-Universe treatment uses the mechanical approach inside the cell of uniformity and studies mGCG fluctuations on top of CDM, radiation, and discrete structures. In that framework, consistency of the Poisson-type equation at large scale factor restricts

$$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$1

to three admissible cases: $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$2 with $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$3, $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$4, or $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$5; all other ranges, including $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$6, $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$7, and $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$8, are ruled out (Bouhmadi-López et al., 2015). The physical content of these cases differs: for $$p_{\rm Ch}=\beta\,\varepsilon_{\rm Ch}-(1+\beta)\frac{A}{\varepsilon_{\rm Ch}^\alpha},$$9 the mGCG behaves as an unclustered cosmological constant, for $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$0 fluctuations may exist but the source of the gravitational potential vanishes, and for $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$1 one obtains a clustered consistent solution (Bouhmadi-López et al., 2015).

By contrast, when MGCG is treated as a dark-energy component rather than the entire dark sector, combined background and growth analyses can favor small but nonzero barotropic corrections. A representative fit gives

$$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$2

with

$$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$3

and reports positive $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$4, a present growth rate $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$5, growth index $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$6, and present deceleration $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$7 (Paul et al., 2014).

This suggests that the empirical status of MGCG depends sensitively on how the fluid is deployed. The strongest exclusions arise when a hydrodynamical MGCG is required to account simultaneously for dark matter and dark energy in the perturbation sector, whereas more permissive results appear when the model is used as a dark-energy component on top of CDM, or when late-time perturbations are analyzed in alternative approximations (Fabris et al., 2010, Paul et al., 2014, Bouhmadi-López et al., 2015).

6. Observational constraints, diagnostics, and later extensions

Several data combinations have been used to constrain MGCG and related mGCG parameterizations. The following results are representative rather than exhaustive.

Reference	Setup	Selected results
(Paul et al., 2014)	Background + growth + $$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$8	$$p=-\frac{\alpha}{1+\alpha}\rho-\frac{\vartheta(z)}{1+\alpha}\rho^{-\alpha},$$9, $p=B\rho-A/\rho^\alpha$00, $p=B\rho-A/\rho^\alpha$01
(Deng, 2011)	Unified DM–DE revisited	$p=B\rho-A/\rho^\alpha$02, $p=B\rho-A/\rho^\alpha$03
(Benaoum et al., 2019)	GTD + Union 2.1 SNe Ia	$p=B\rho-A/\rho^\alpha$04, $p=B\rho-A/\rho^\alpha$05, $p=B\rho-A/\rho^\alpha$06, $p=B\rho-A/\rho^\alpha$07
(Bhardwaj et al., 25 Sep 2025)	MGCG + matter creation + bulk viscosity	DS1: $p=B\rho-A/\rho^\alpha$08, $p=B\rho-A/\rho^\alpha$09, $p=B\rho-A/\rho^\alpha$10

In the unified revisited model, the present deceleration parameter and transition redshift are reported as

$p=B\rho-A/\rho^\alpha$11

while equality between effective dark matter and dark energy occurs at

$p=B\rho-A/\rho^\alpha$12

The same paper applies the Statefinder diagnostic $p=B\rho-A/\rho^\alpha$13, where $p=B\rho-A/\rho^\alpha$14CDM is the fixed point $p=B\rho-A/\rho^\alpha$15, quiessence corresponds to vertical lines $p=B\rho-A/\rho^\alpha$16, and MGCG trajectories start at $p=B\rho-A/\rho^\alpha$17 at early times (Deng, 2011). It also finds that density fluctuations begin to deviate from linear growth around $p=B\rho-A/\rho^\alpha$18, attributed to the onset of dark-energy dominance (Deng, 2011).

The GTD derivation of the modified Chaplygin gas produces a normalized Hubble function

$p=B\rho-A/\rho^\alpha$19

and reports a high-quality fit to the Union 2.1 supernova compilation, “virtually indistinguishable from $p=B\rho-A/\rho^\alpha$20CDM for $p=B\rho-A/\rho^\alpha$21” (Benaoum et al., 2019). In the dissipative extension with matter creation rate $p=B\rho-A/\rho^\alpha$22 and bulk viscous pressure $p=B\rho-A/\rho^\alpha$23, the total entropy rate $p=B\rho-A/\rho^\alpha$24 remains positive throughout cosmic history, while $p=B\rho-A/\rho^\alpha$25 changes sign around $p=B\rho-A/\rho^\alpha$26 (Bhardwaj et al., 25 Sep 2025).

Recent work has also transported MGCG beyond homogeneous cosmology. In a modified-gravity interpretation, a black-hole spacetime associated with MGCG is asymptotically non-flat, possesses two distinct horizons, and remains stable under scalar and electromagnetic perturbations; the quasinormal-mode spectrum depends sensitively on the MGCG parameters, suggesting a possible ringdown probe of the model (Bohra, 17 Sep 2025). Together with the thermodynamic, scalar-field, and perturbative literature, these developments place MGCG at the intersection of unified dark-sector phenomenology, horizon thermodynamics, and effective-field reconstructions rather than within a single settled cosmological paradigm.