The stakes come from a subtle but consequential feature of the DESI data. When the latest baryon-acoustic-oscillation measurements are combined with cosmic microwave background constraints and Type Ia supernova distances, the data can be summarized by a phenomenological form in which the dark energy equation of state, written as w(a) = w0 + wa(1 − a), prefers w0 slightly greater than −1 and wa negative. That combination implies that the effective equation of state crosses the so-called phantom divide at w = −1 at some point in the recent expansion history. For theorists, this crossing is notoriously difficult. A canonical scalar field with a standard kinetic term can never push w below −1, while a field with a wrong-sign kinetic term can do so but at the cost of ghost instabilities, quantum states of negative norm that signal a breakdown of the theory. Stable phantom-crossing models therefore demand extra structure: additional degrees of freedom, effective-field-theory operators, modified gravity, or scalar-tensor dynamics.

Li and Liu’s proposal builds on the running-curvaton framework, in which a light scalar field called the curvaton contributes to the primordial curvature perturbation that seeded all structure in the universe. In the original, minimally coupled version of the model, the same field was expected to survive to the present day and act as dark energy, but its late-time behavior resembled thawing quintessence, a form of dynamical dark energy that cannot naturally reproduce the DESI-motivated phantom-crossing region. The authors’ key move is to add a Jordan-frame non-minimal coupling, a term of the form ξχ²R that ties the curvaton field directly to the Ricci scalar, a measure of spacetime curvature. This coupling is not exotic by the standards of modern field theory; scalar fields coupled to curvature have been studied for decades in inflation and scalar-tensor gravity, and such theories belong to the well-studied Horndeski class. The novelty lies in embedding the coupling within the running-curvaton scenario in a way that preserves the model’s early-universe successes.

That preservation is the first technical hurdle the paper clears. During slow-roll inflation, the Ricci scalar is approximately 12H², where H is the Hubble expansion rate, so the non-minimal coupling shifts the effective mass of the curvaton fluctuation from its original value to m_eff² = 6H²(g0 + 2ξ), where g0 parametrizes the inflaton-curvaton interaction. Crucially, the observable spectral tilt of the curvaton perturbations depends only on the combination g0 + 2ξ. The authors show that this shift can be absorbed, at leading order, by retuning the mass parameter: one simply defines an observed coupling g0_obs = g0 + 2ξ and adjusts g0 accordingly. The predicted spectral index n_χ − 1 ≈ 4g0_obs then matches the minimally coupled model exactly. Similarly, the local non-Gaussianity parameter f_NL, which measures departures from Gaussian statistics in the primordial perturbations, is controlled mainly by the curvaton’s energy fraction at decay, with f_NL ≈ 5/(4r_dec) in the sudden-decay approximation. Because spacetime is nearly flat during radiation domination, R ≈ 0, the curvature-induced mass correction is suppressed precisely when it matters most for non-Gaussianity, and the Planck satellite’s observational bounds remain satisfied. Tensor perturbations fare equally well: since the curvaton is a subdominant spectator during inflation, the correction to the primordial gravitational-wave amplitude is of order ξχ²/M_P², negligible for sub-Planckian field values.

The second hurdle is the late-time dynamics, and this is where the non-minimal coupling earns its keep. After the inflaton decays, the surviving part of the curvaton potential, a phenomenological form V(χ) = V1(1 − V0 e^(−λχ/M_P)), drives the field’s dark-energy-like evolution. When the modified Friedmann equations are recast in the standard Einstein form, the effective dark-energy pressure acquires geometric terms involving the coupling, such as −ξχ²(2Ḣ + 3H²) and −2ξ(χχ̈ + χ̇²). The sum of the effective dark-energy density and pressure, which determines whether w is above or below −1, becomes ρ+p = (1 − 2ξ)χ̇² + 2ξ(Hχχ̇ − χχ̈ − χ²Ḣ). The second group of terms can become sufficiently negative during late-time evolution to push the effective equation of state across the phantom divide, all without introducing a phantom kinetic term. The theory retains a canonical kinetic structure, and because the Horndeski function G4 depends on the field but not on its kinetic variable, the speed of gravitational waves remains exactly luminal, c_T² = 1, consistent with the landmark constraint from the neutron-star merger GW170817.

To test whether these theoretical trajectories are actually compatible with observations, the authors performed a Markov-chain Monte Carlo likelihood analysis using numerical solutions of the model’s homogeneous background equations. The observational input combined three pillars: a reduced CMB prior on the acoustic angular scale and the physical baryon and total matter densities, the official DESI DR2 baryon-acoustic-oscillation data vector with its full 13-by-13 covariance matrix, and the Pantheon+ supernova compilation with 1580 objects, deliberately excluding the SH0ES Cepheid calibration to avoid entangling the analysis with the Hubble-tension debate. The sampled parameters were the dimensionless Hubble parameter h, the baryon and cold dark matter densities, and four microscopic quantities of the curvaton sector: the coupling ξ, the potential slope λ, the potential normalization V0, and the initial field value χ_i. The late-time normalization V1 was fixed for each sample by a shooting condition requiring the normalized expansion rate to equal unity today.

The results are encouraging but carefully circumscribed. The combined data yield H0 = 67.88 km/s/Mpc with an uncertainty of roughly half a percent, a matter density Ωm = 0.3072, and dark energy parameters w0 = −0.922 and wa = −0.205, placing the posterior squarely in the DESI-motivated phantom-crossing region. The best-fit point improves the chi-square relative to a flat ΛCDM reference fit by Δχ² = −6.56, but because the curvaton model carries four additional microscopic parameters, the Akaike information criterion actually favors ΛCDM by ΔAIC = +1.44. The authors are explicit about this: the model is compatible with the distance data and can modestly improve the best fit, but the analysis does not constitute a decisive statistical preference over the standard cosmological constant. The posterior bands show the median equation of state starting above −1 today and evolving toward a mildly phantom regime at intermediate redshift, while the expansion-rate history H(z)/H0 remains tightly constrained, demonstrating that the phantom-crossing mechanism does not spoil the distance information.

Equally instructive is what the data cannot yet pin down. The derived background quantities, including H0, Ωm, the sound horizon r_d, w0 and wa, are compressed sharply by the likelihood, but the microscopic parameters ξ, λ, V0 and χ_i remain broad, skewed and long-tailed in the marginalized posterior. A quantitative compression diagnostic makes the point vivid: the 68 percent posterior width of ξ spans about 61 percent of its allowed prior range, whereas the corresponding ratios for h and the cold dark matter density are only about 2 to 3 percent. A supplementary run widening the prior on ξ from 5 to 8 broadens the high-ξ tail further without the chain running away, confirming that the current bound is prior-conditional. This degeneracy is not a failure of the mechanism but a genuine limitation of background-level data, which constrain the expansion history far more tightly than the underlying microphysics that produces it.

The authors are equally candid about the analysis’s boundaries. The CMB information was compressed to three early-background numbers, so the calculation does not include the non-minimal coupling’s effects on the full temperature and polarization spectra, CMB lensing, the late integrated Sachs-Wolfe effect, or scalar perturbation evolution. A complete test will require implementing the model in a Boltzmann code, the standard machinery used to compute CMB anisotropies from first principles. Local gravity constraints also demand a separate post-Newtonian or environmental analysis: in high-density regions the scalar mass receives a contribution proportional to the local matter density, suggesting possible Yukawa suppression of fifth-force effects analogous to chameleon or symmetron screening, but no definitive Solar-System viability claim is made. Within those limits, the study delivers a clean proof of concept: the same field that may have shaped the seeds of galaxies in the infant universe could, through its intimate coupling to spacetime curvature, be steering the accelerated expansion we observe today, and it can do so while crossing the phantom divide that has long frustrated simpler models of dark energy.

Subject of Research: A non-minimally coupled running curvaton model linking primordial perturbations to phantom-crossing dynamical dark energy constrained by DESI, CMB and supernova data

Article Title: Non-minimally coupled running curvaton for DESI-motivated dynamical dark energy

Article References: Li, B., & Liu, L.-H. (2026). Non-minimally coupled running curvaton for DESI-motivated dynamical dark energy. The European Physical Journal C, 86(10), Article 1156. https://doi.org/10.1140/epjc/s10052-026-16396-3

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16396-3

Keywords: dark energy, curvaton, non-minimal coupling, DESI, phantom divide, cosmology, scalar-tensor gravity, baryon acoustic oscillations, Pantheon+, cosmic microwave background, Horndeski theory, MCMC

The study, led by Fan Yang of Hunan University of Science and Technology together with colleagues in China, examines interacting dark energy, a class of models in which the dark sector is not passive but conversational. Instead of evolving independently, dark matter and dark energy exchange energy through a non-gravitational coupling, a handoff that could naturally produce the effective dynamical dark energy behavior hinted at by recent surveys. The researchers considered two versions of this idea: an interacting variant of wCDM, where dark energy’s pressure-to-density ratio, its equation of state w, is a constant, and an interacting version of the Chevallier-Polarski-Linder model, in which w evolves over cosmic time according to two parameters, w0 and wa.

The technical engine of the analysis is the growth of matter density perturbations, the subtle wrinkles in the early universe that gravity amplified into galaxies and galaxy clusters. Cosmologists summarize this growth with the growth rate f, defined as the logarithmic derivative of the matter density contrast with respect to the scale factor. A widely used shorthand links f to the matter density fraction via a power law, f approximately equal to Omega_m raised to the growth index gamma. In standard Lambda-CDM governed by General Relativity, gamma sits at roughly 0.545. The beauty of this number is its diagnostic power: modified gravity theories predict different values, with the DGP braneworld model yielding about 0.6875 and viable f(R) gravity models clustering near 0.40 to 0.43.

Here is where the physics becomes genuinely subtle. The team derived a second-order approximation for the growth index in the interacting wCDM scenario and uncovered a clean theoretical result: the dark-sector coupling alpha induces a nearly universal correction to the growth index of about 1.1 times alpha, in both the wCDM and CPL frameworks. That means a positive coupling, in which energy flows from dark energy into dark matter, pushes the growth index upward toward modified gravity territory, while a negative coupling drags it down toward the f(R) range. Theoretical work dating back to studies of interacting quintessence has shown that such a coupling can even be engineered to reproduce both the expansion history and the growth history of DGP gravity exactly, rendering the two fundamentally different explanations observationally indistinguishable from expansion and growth data alone. If the coupling is free to vary, the growth index loses its diagnostic punch.

To determine whether observations actually permit this mimicry, the researchers confronted both models with an arsenal of modern cosmological data. The background expansion probes included the Pantheon+ compilation of 1701 Type Ia supernova light curves spanning redshifts from 0.001 to 2.26, thirty-two Hubble parameter measurements from cosmic chronometers, compressed distance priors from the Planck 2018 cosmic microwave background release, and two alternative baryon acoustic oscillation datasets, one drawn from all four generations of the Sloan Digital Sky Survey and one from the second data release of the Dark Energy Spectroscopic Instrument. Crucially, they folded in redshift-space distortion measurements of the quantity f sigma 8, which combines the growth rate with the amplitude of matter fluctuations and is immune to the galaxy bias that plagues raw growth observations.

On the theoretical side, the team was careful to handle a known pitfall. In the interacting case, the simple constant-gamma approximation breaks down: at low redshifts the discrepancy between the parameterization and the true growth rate reaches roughly 12 percent for a coupling of 0.02, comparable to the typical 15 percent observational error on f sigma 8. Their improved second-order formula for gamma(z), which explicitly includes the coupling alpha, cuts that maximum error to about 5 percent, a roughly 60 percent improvement. They validated this analytical shortcut against exact numerical integration of the modified growth equation within Markov Chain Monte Carlo analyses, finding robust agreement in both central values and uncertainties, and gaining substantial computational efficiency in the process. The analysis also explicitly included the radiation component in the background equations to guarantee accuracy at high redshifts, which matters for the CMB distance priors.

The results are a triumph of precision. Using the DESI DR2 baryon acoustic oscillation data combined with the other background probes and redshift-space distortions, the interacting wCDM model yields a coupling constant of alpha = 0.0230 plus or minus 0.0076, consistent with zero at approximately the three sigma confidence level, and an equation of state w = -0.948 plus or minus 0.025, consistent with the cosmological constant value of -1 at two sigma. The dynamical CPL version tells a similar story: alpha = 0.0173 plus or minus 0.0095, with w0 = -0.888 plus or minus 0.059 and wa = -0.28 plus or minus 0.25, all consistent with their Lambda-CDM values within two sigma. In short, there is no definitive statistical evidence that dark matter and dark energy are talking to each other, though the DESI-based constraint on the coupling is roughly 85 percent tighter than comparable earlier analyses.

The implications for the degeneracy with modified gravity are decisive. Translating the observational bounds on alpha into the growth index, the DESI-based wCDM analysis restricts gamma to the interval from about 0.5469 to 0.5974 at three sigma, while the CPL model yields roughly 0.53 to 0.59. Taken together, both models confine the growth index to a common window of approximately 0.53 to 0.60. That window is far too narrow to reach either the DGP prediction near 0.69 or the f(R) predictions near 0.40 to 0.43. The region of parameter space where interacting dark energy could masquerade as modified gravity is strongly disfavored by current data, which means the growth index regains its power as a diagnostic despite the theoretical degeneracy.

The analysis also revealed telling correlations among the dark sector parameters. In the interacting wCDM model, the coupling alpha and the equation of state w are positively correlated: a stronger energy transfer to dark matter prefers a less negative equation of state to preserve the observed expansion history. In the CPL model the pattern grows more intricate, with a negative correlation between w0 and alpha but a positive one between wa and alpha. Interestingly, allowing the equation of state to vary in time partially absorbs the signatures that would otherwise be attributed to the interaction, shifting the best-fit coupling closer to zero and slightly widening its uncertainty. This warns that fixing any of these parameters to their Lambda-CDM values could bias the others and mask the true underlying mechanism, so future analyses must constrain the dark sector jointly.

Looking ahead, the study sharpens the interpretive framework for the coming decade of structure surveys. The Dark Energy Spectroscopic Instrument’s dedicated redshift-space distortion analysis, the Euclid space telescope, and the Square Kilometre Array will measure the growth of cosmic structure with unprecedented precision. The message from this work is clear: if those experiments detect a statistically significant deviation of the growth index from the Lambda-CDM value of about 0.55, the culprit will most likely be a genuine modification of gravity on cosmological scales rather than a clandestine exchange between dark matter and dark energy, whose permitted shifts now appear too small to account for such a signal. Meanwhile, with recent DESI and Dark Energy Survey analyses pointing toward mildly evolving dark energy, the question of what actually drives cosmic acceleration remains wide open, and the universe’s dark components will keep their secrets a little longer.

Subject of Research: Observational constraints on interacting dark energy models and the growth index of matter density perturbations

Article Title: Impact of interacting dark energy on the growth of matter density perturbations: observational constraints from DESI and multi-probe data

Article References: Yang, F., Zhai, R., Fu, X., Xu, B., Liu, K., Ding, C., & Huang, Y. (2026). Impact of interacting dark energy on the growth of matter density perturbations: observational constraints from DESI and multi-probe data. The European Physical Journal C, 86(9), Article 1115. https://doi.org/10.1140/epjc/s10052-026-16366-9

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16366-9

Keywords: dark energy, dark matter, DESI, growth index, modified gravity, baryon acoustic oscillations, redshift-space distortions, Pantheon+, cosmic microwave background, Lambda-CDM, CPL parameterization, cosmology