pith. sign in

arxiv: 2605.18240 · v1 · pith:YUAIVSN7new · submitted 2026-05-18 · 🌌 astro-ph.CO

Quadratic Dark Energy Phase-Space Dynamics and Analysis

Sahar Mohammadi , Ebrahim Yusofi , Kosar Asadi This is my paper

Pith reviewed 2026-05-20 00:35 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords quadratic dark energyphase space analysisphantom attractorsdynamical systemsequation of statecosmological dynamicslate-time accelerationDESI observations
0
0 comments X

The pith

Negative values of the dynamical parameter steer quadratic dark energy to stable phantom attractors approached asymptotically without crossing.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies a minimal extension of dark energy in which the pressure contains a term quadratic in the energy density. This introduces a dynamical parameter η that controls transitions between cosmic regimes. Phase-space analysis of the resulting autonomous system shows that negative η produces stable attractors with phantom-like behavior while positive η produces unstable repellers. The effective equation of state approaches minus one from both the quintessence and phantom sides but never actually crosses it. If this holds, the model can generate stronger late-time acceleration and higher Hubble rates that observations could test.

Core claim

In the quadratic dark energy model the pressure includes a nonlinear term proportional to the square of the energy density. Dynamical systems analysis identifies critical points whose stability depends on the sign of η, with negative values corresponding to stable phantom sinks and positive values to unstable sources. Trajectories in phase space approach the phantom divide at w_eff equals minus one from both quintessence and phantom regimes without crossing, producing enhanced Hubble expansion and more pronounced late-time acceleration.

What carries the argument

The autonomous dynamical system obtained from the Friedmann equations with the added quadratic pressure term, where η acts as the parameter that sets the stability of the late-time attractors.

If this is right

  • Stable phantom attractors produce higher expansion rates at late times than a cosmological constant.
  • The non-crossing approach to w_eff equals minus one supplies an alternative to models that require phantom crossing to fit evolving dark energy data.
  • Positive η values make the corresponding regimes transient and repelling.
  • The enhanced acceleration from negative-η attractors can be directly compared with DESI constraints on late-time expansion.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Quadratic pressure extensions could be examined in other dark energy models to achieve similar non-crossing phantom-like behavior.
  • High-precision Hubble parameter measurements at low redshifts would provide a direct test of the predicted enhancement from stable attractors.
  • Adding early-universe degrees of freedom might alter the late-time attractor structure in ways not explored here.

Load-bearing premise

The pressure contains a nonlinear term exactly proportional to the square of the energy density and the resulting autonomous system fully captures the late-time cosmological evolution.

What would settle it

Future measurements that detect an actual crossing of the phantom divide in the effective equation of state at low redshift would rule out the claimed asymptotic non-crossing behavior.

Figures

Figures reproduced from arXiv: 2605.18240 by Ebrahim Yusofi, Kosar Asadi, Sahar Mohammadi.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
read the original abstract

We present a comprehensive phase-space analysis of a quadratic dark energy model where the pressure includes a nonlinear term proportional to the square of the energy density. This minimal extension beyond the $\Lambda$CDM framework introduces a dynamical parameter $\eta(z)$ that governs transitions between different cosmological regimes. Through dynamical systems theory, we identify critical points and their stability properties, revealing that negative $\eta$ values drive the system toward stable phantom attractors (sinks), while positive values correspond to unstable repellers (sources). The model exhibits a distinctive asymptotic approach to the phantom divide ($w_{\rm eff}=-1$) from both quintessence and phantom sides without actual crossing, providing a non-crossing alternative to the phantom-crossing behavior preferred by recent DESI DR2 constraints. Our analysis shows that stable phantom attractors produce enhanced Hubble expansion rates and more pronounced late-time acceleration, features that can be compared with recent DESI observations suggesting evolving dark energy.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a phase-space analysis of a quadratic dark energy model in which the pressure includes a nonlinear term proportional to the square of the energy density. It introduces a redshift-dependent dynamical parameter η(z) that governs transitions between regimes and applies dynamical systems methods to locate critical points. The central claims are that negative values of η drive trajectories to stable phantom attractors (sinks), positive values produce unstable repellers (sources), and the effective equation of state approaches w_eff = −1 asymptotically from both quintessence and phantom sides without crossing, yielding enhanced late-time acceleration that can be compared with DESI DR2 data.

Significance. If the derivations are internally consistent, the work supplies a minimal, non-crossing dynamical alternative to phantom-crossing models that may align with recent indications of evolving dark energy. The explicit linkage between the sign of η and attractor stability, together with the predicted enhancement of the Hubble rate at late times, offers concrete, observationally testable features.

major comments (2)
  1. [§2] §2 (Model setup): The pressure is taken to contain a term exactly proportional to ρ_de², after which the system is declared autonomous. However, η is introduced as η(z) with explicit redshift dependence. No auxiliary equation for dη/dN is supplied to close the system; without it the vector field acquires explicit time dependence, the critical points are not stationary, and the eigenvalue signs cannot be assigned solely by the instantaneous value of η.
  2. [§4] §4 (Linear stability analysis): The classification that negative η produces sinks while positive η produces sources treats η as a fixed bifurcation parameter. If η(z) varies along trajectories, a given orbit may cross regions of changing stability; the abstract’s stability conclusions therefore require an explicit demonstration that the system remains autonomous or that the auxiliary dynamics for η preserve the reported eigenvalue structure.
minor comments (2)
  1. [Abstract] The abstract states stability conclusions and the non-crossing property but supplies neither the phase-space variables (x, y, …) nor the explicit autonomous equations; adding one or two key equations would improve readability.
  2. [Figures] Figure captions should explicitly state the sign of η used for each phase portrait and indicate whether the flow is shown for constant or varying η.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which help clarify the presentation of the dynamical systems analysis. We address each major comment below and indicate the revisions planned for the next version.

read point-by-point responses

  1. Referee: [§2] §2 (Model setup): The pressure is taken to contain a term exactly proportional to ρ_de², after which the system is declared autonomous. However, η is introduced as η(z) with explicit redshift dependence. No auxiliary equation for dη/dN is supplied to close the system; without it the vector field acquires explicit time dependence, the critical points are not stationary, and the eigenvalue signs cannot be assigned solely by the instantaneous value of η.

    Authors: We appreciate this observation. The parameter η(z) is introduced phenomenologically to govern transitions between regimes. For the phase-space analysis in the manuscript, we treat η as a fixed bifurcation parameter within each regime to locate and classify the critical points. This is a standard approach when exploring the qualitative behavior for different signs of η. We acknowledge that, strictly, the absence of an auxiliary equation for η renders the system non-autonomous. In the revised manuscript we will clarify this point explicitly in §2, state that the reported critical points and stability apply for constant η, and add a short discussion of how slow variation of η would affect trajectories. We will also outline a minimal auxiliary equation (e.g., a simple linear evolution for η with N) that could close the system in an extended autonomous formulation, without altering the core stability results for fixed η. revision: yes

  2. Referee: [§4] §4 (Linear stability analysis): The classification that negative η produces sinks while positive η produces sources treats η as a fixed bifurcation parameter. If η(z) varies along trajectories, a given orbit may cross regions of changing stability; the abstract’s stability conclusions therefore require an explicit demonstration that the system remains autonomous or that the auxiliary dynamics for η preserve the reported eigenvalue structure.

    Authors: We agree that continuous variation of η could in principle allow orbits to traverse regions of differing stability. Our analysis and numerical examples focus on the instantaneous eigenvalue structure for a given sign of η and show rapid convergence to the attractor once dark energy dominates. In the revised version we will expand §4 with a discussion of the relevant time-scale separation (late-time acceleration occurs on a timescale shorter than significant η evolution in the model) and add supporting numerical evidence that trajectories reach the reported attractor before η changes appreciably. This preserves the abstract’s conclusions for the cosmologically relevant epochs while acknowledging the limitation for fully time-dependent η. revision: yes

Circularity Check

0 steps flagged

Phase-space stability analysis is self-contained and independent of inputs

full rationale

The paper defines a quadratic dark energy model with pressure containing a nonlinear term proportional to energy density squared, introduces the dynamical parameter η(z) to govern regime transitions, and applies standard linear stability analysis to identify critical points and their attractors in the phase space. The classification of negative η leading to stable phantom sinks and positive η to repellers follows directly from the Jacobian eigenvalues at those points under the model's autonomous equations, without any reduction of the output to a fitted input, self-definition, or self-citation chain. No load-bearing step equates a prediction to its own construction; the derivation remains independent and externally falsifiable via comparison to cosmological data such as DESI constraints.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that a quadratic term in pressure plus a single dynamical parameter η(z) is sufficient to describe the phase-space structure; no new particles or forces are postulated, but the functional form of the quadratic term and the definition of η(z) function as modeling choices.

free parameters (1)
  • η(z)
    Dynamical parameter whose sign determines attractor versus repeller behavior; its explicit redshift dependence is not derived from first principles in the abstract.
axioms (1)
  • domain assumption Standard flat FLRW metric and continuity equations for dark energy fluid
    Invoked to set up the autonomous dynamical system for phase-space analysis.

pith-pipeline@v0.9.0 · 5690 in / 1380 out tokens · 49313 ms · 2026-05-20T00:35:37.924000+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · 5 internal anchors

  1. [1]

    The quadratic model provides a natural mechanism for such evolution while avoiding theoretical pathologies associated with phantom divide crossing

    suggests a preference for evolving dark energy, making our dynamical analysis particu- larly timely. The quadratic model provides a natural mechanism for such evolution while avoiding theoretical pathologies associated with phantom divide crossing. Crucially, our model’s asymptotic approach without crossing aligns with DESI DR2 findings that, while favori...

    work page

  2. [2]

    Negative values ofηdrive the system toward stable phantom attractors (sinks), while positive values correspond to unstable repellers (sources) or saddle points in the phase space

    The quadratic parameterη(z) functions as a fundamental control parameter governing cosmic energy flow. Negative values ofηdrive the system toward stable phantom attractors (sinks), while positive values correspond to unstable repellers (sources) or saddle points in the phase space

    work page

  3. [3]

    This provides a viable 12 (a)w=−0.87,w a =−0.46,w eff <−1 (b)w=−0.8,w a =−0.1,w eff >−1 FIG

    The model exhibits asymptotic approach to the phantom divide (w eff =−1) from both quintessence and phantom sides without actual crossing. This provides a viable 12 (a)w=−0.87,w a =−0.46,w eff <−1 (b)w=−0.8,w a =−0.1,w eff >−1 FIG. 4:Effective equation of state (Eq. 5), for the quadratic dark energy model. The phantom case (left) remains beloww eff =−1 an...

    work page

  4. [4]

    Stable phantom attractors produce distinctive observational signatures, including en- hanced Hubble expansion rates at low redshifts and more pronounced late-time accel- eration, as evidenced in the deceleration parameterq(z)

    work page

  5. [5]

    The quadratic model provides a natural mechanism for this behavior within a minimal phenomenological framework

    The dynamical features align remarkably well with recent DESI constraints, which indicate a preference for evolving dark energy that approaches, but does not necessarily cross, the phantom divide. The quadratic model provides a natural mechanism for this behavior within a minimal phenomenological framework. The quadratic dark energy model thus represents ...

    work page

  6. [6]

    Measurements of Ω and Λ from 42 High-Redshift Supernovae,

    S. Perlmutter, G. Aldering, G. Goldhaber, R. A. Knop, P. Nugent, P. G. Castro, S. Deustua, S. Fabbro, A. Goobar, D. E. Groom, I. M. Hook, A. G. Kim, M. Y. Kim, J. C. Lee, N. J. Nunes, R. Pain, C. R. Pennypacker, R. Quimby, C. Lidman, R. S. Ellis, M. Irwin, R. G. McMahon, P. Ruiz-Lapuente, N. Walton, B. Schaefer, B. J. Boyle, A. V. Filippenko, T. Matheson,...

    work page 1999

  7. [7]

    Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant,

    A. G. Riess, A. V. Filippenko, P. Challis, A. Clocchiatti, A. Diercks, P. M. Garnavich, R. L. Gilliland, C. J. Hogan, S. Jha, R. P. Kirshner, B. Leibundgut, M. M. Phillips, D. Reiss, B. P. Schmidt, R. A. Schommer, R. C. Smith, J. Spyromilio, C. Stubbs, N. B. Suntzeff, and J. Tonry, “Observational Evidence from Supernovae for an Accelerating Universe and a...

    work page 1998

  8. [8]

    Dynamics of Dark Energy,

    E. J. Copeland, M. Sami, and S. Tsujikawa, “Dynamics of Dark Energy,” International Journal of Modern Physics D15, 1753 (2006)

    work page 2006

  9. [9]

    DESI DR2 Results II: Measurements of Baryon Acoustic Oscillations and Cosmological Constraints

    M. Abdul Karim et al. (DESI Collaboration), “DESI DR2 Results II: Measurements of Baryon Acoustic Oscillations and Cosmological Constraints,” arXiv:2503.14738 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  10. [10]

    Dynamical systems applied to cosmology: dark energy and modified gravity

    S. Bahamonde, C. G. B¨ ohmer, S. Carloni, E. J. Copeland, W. Fang, and N. Tamanini, “Dy- namical systems applied to cosmology: Dark energy and modified gravity,” Physics Reports 775-777, 1 (2018), arXiv:1712.03107

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  11. [11]

    Cosmological imprint of an energy compo- nent with general equation of state,

    R. R. Caldwell, R. Dave, and P. J. Steinhardt, “Cosmological imprint of an energy compo- nent with general equation of state,” Physical Review Letters80, 1582 (1998), arXiv:astro- ph/9708069. 15

    work page arXiv 1998

  12. [12]

    Graduated inflationary universes,

    J. D. Barrow, “Graduated inflationary universes,” Physics Letters B235, 40 (1990)

    work page 1990

  13. [13]

    The final state and thermodynamics of dark energy universe

    S. Nojiri and S. D. Odintsov, “The final state and thermodynamics of dark energy universe,” Physical Review D70, 103522 (2004), arXiv:hep-th/0408170

    work page internal anchor Pith review Pith/arXiv arXiv 2004

  14. [14]

    "Expansion" around the vacuum equation of state - sudden future singularities and asymptotic behavior

    H. Stefancic, “Expansion around the vacuum equation of state – Sudden future singularities and asymptotic behavior,” Physical Review D71, 084024 (2005), arXiv:astro-ph/0411630

    work page internal anchor Pith review Pith/arXiv arXiv 2005

  15. [15]

    Cosmo-dynamics and dark energy with a quadratic EoS: anisotropic models, large-scale perturbations and cosmological singularities

    K. N. Ananda and M. Bruni, “Cosmo-dynamics and dark energy with a quadratic EoS: Anisotropic models, large-scale perturbations and cosmological singularities,” Physical Re- view D74, 023524 (2006), arXiv:gr-qc/0603131

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  16. [16]

    Growth of perturbations in an expanding universe with Bose-Einstein con- densate dark matter,

    P. H. Chavanis, “Growth of perturbations in an expanding universe with Bose-Einstein con- densate dark matter,” Astronomy and Astrophysics537, A127 (2012)

    work page 2012

  17. [17]

    On a quadratic equation of state and a universe mildly bouncing above the Planck temperature,

    J. Berteaud, J. Pasquet, T. Sch¨ ucker, and A. Tilquin, “On a quadratic equation of state and a universe mildly bouncing above the Planck temperature,” JCAP10, 069 (2019), arXiv:1807.05068

    work page arXiv 2019

  18. [18]

    Stability and Thermodynamics of a Generalized Power-Law Dark Energy Model,

    S. Kazemi, M. A. Ramzanpour, E. Yusofi, and A. R. Amani, “Stability and Thermodynamics of a Generalized Power-Law Dark Energy Model,” International Journal of Geometric Methods in Modern Physics0, 2650064 (2025), https://doi.org/10.1142/S0219887826500647

    work page doi:10.1142/s0219887826500647 2025

  19. [19]

    Quadratic equation of state for cosmic acceleration: Entropy evolution and phantom crossing,

    A. Shahriar, M. Abbasiyan-Motlaq, M. Mohsenzadeh, and E. Yusofi, “Quadratic equation of state for cosmic acceleration: Entropy evolution and phantom crossing,” Physical Review D 112, 083520 (2025)

    work page 2025

  20. [20]

    Testing the dynamical stability and the validity of the generalized second law within the phantom dynamical dark energy model,

    Naseeba K. M., S. Nelleri, and N. Poonthottathil, “Testing the dynamical stability and the validity of the generalized second law within the phantom dynamical dark energy model,” Physical Review D109, 063528 (2024), arXiv:2308.03084

    work page arXiv 2024

  21. [21]

    Dynamical stability and phase space analysis of an emergent Universe with non-interacting and interacting fluids,

    B. C. Roy, A. Chanda, and B. C. Paul, “Dynamical stability and phase space analysis of an emergent Universe with non-interacting and interacting fluids,” Classical and Quantum Gravity41, 175009 (2024), arXiv:2401.00782

    work page arXiv 2024

  22. [22]

    Stability analysis of a dark energy model in Rastall gravity,

    Shaily, A. Singh, J. K. Singh, S. Hussain, and R. Myrzakulov, “Stability analysis of a dark energy model in Rastall gravity,” arXiv:2402.08709 (2024)

    work page arXiv 2024

  23. [23]

    Stability analysis of warm quintessential dark energy model,

    S. Das, S. Hussain, D. Nandi, R. O. Ramos, and R. Silva, “Stability analysis of warm quintessential dark energy model,” Physical Review D108, 083517 (2023), arXiv:2306.09369

    work page arXiv 2023

  24. [24]

    Modeling and stability analysis of dark energy ultra-compact objects in extended teleparallel gravity,

    M. R. Shahzad, L. Fakhar, K. Nabi, Z. Amjad, A. M. Mubaraki, and A. A. Yagob, “Modeling and stability analysis of dark energy ultra-compact objects in extended teleparallel gravity,” 16 European Physical Journal Plus140, 690 (2025)

    work page 2025

  25. [25]

    A Dark Matter Model with Quadratic Equation of State: Background Evolution and Structure Formation,

    K. Rezazadeh, E. Yusofi, and A. Talebian, “A Dark Matter Model with Quadratic Equation of State: Background Evolution and Structure Formation,” arXiv:2509.11138 (2025)

    work page arXiv 2025

  26. [26]

    Observational constraints on the dark energy with a quadratic equation of state,

    H. Moshafi, A. Talebian, E. Yusofi, and E. Di Valentino, “Observational constraints on the dark energy with a quadratic equation of state,” Physics of the Dark Universe45, 101524 (2024), arXiv:2403.02000

    work page arXiv 2024

  27. [27]

    Surface tension of cosmic voids as a possible source for dark energy,

    E. Yusofi, M. Khanpour, B. Khanpour, M. A. Ramzanpour, and M. Mohsenzadeh, “Surface tension of cosmic voids as a possible source for dark energy,” Monthly Notices of the Royal Astronomical Society511, L82 (2022), arXiv:1907.12418

    work page arXiv 2022

  28. [28]

    A possible role for the merger of clusters/voids in the cosmological expansion,

    S. Mohammadi, E. Yusofi, M. Mohsenzadeh, and M. K. Salem, “A possible role for the merger of clusters/voids in the cosmological expansion,” Monthly Notices of the Royal Astronomical Society525, 3274 (2023)

    work page 2023

  29. [29]

    Favale, A

    A. Favale, A. G´ omez-Valent, and M. Migliaccio, “Cosmic chronometers to calibrate the ladders and measure the curvature of the Universe. A model-independent study,” Monthly Notices of the Royal Astronomical Society523, 3406 (2023), arXiv:2301.09591

    work page arXiv 2023