{"id":"52954bdb-e47c-49e2-af14-029fc7c86840","arxiv_id":"2501.14909","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"A cuscuton-like kinetic term is combined with quintessence fields in a first-order formalism, producing analytical dark-energy histories, but the advertised phantom phases rely on an inconsistent sign choice and the AIC comparison does not actually favor the model.","lead":"This paper adds a cuscuton term, a nonstandard kinetic contribution, to the standard dark energy scalar field and builds a first-order framework that yields analytical cosmic expansion histories. It also fits one model to supernova, BAO, and cosmic chronometer data and compares it with ΛCDM using the Akaike criterion, making it relevant to the question of whether dark energy is exactly constant.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The hyperbolic-potential solutions rely on an unacknowledged sign flip: Eq. (34) has φdot<0 while the first-order framework assumes sgn(φdot)>0, so the phantom-phase scenarios do not solve the stated field equations.","rationale":"The paper's central new results are the first-order framework and analytical solutions that interpolate between matter/radiation phases and ω=−1, including phantom crossings for hyperbolic potentials. The sign assumption is not cosmetic: the cuscuton term contains |φdot|, so its equation-of-motion contribution is ±α depending on the sign of φdot. The paper explicitly assumes sgn(φdot)>0, yet the hyperbolic solution (34) has φdot<0 on t>0, making Eqs. (6b), (21), (35), and the phantom-phase interpretation inconsistent with the stated Lagrangian. This is a genuine internal inconsistency, not merely a deviation from consensus, and it invalidates the main advertised phenomenological novelty. The reader's weakest_assumption identifies exactly this issue, and my independent reading of Sec. III.B confirms it. A secondary but independent problem is the AIC interpretation in Sec. V: the reported positive ΔAIC values, with ΛCDM as reference, indicate that ΛCDM has lower AIC, so the abstract's claim of 'strong support' for the cuscuton model is contradicted by the paper's own numbers. Both issues are load-bearing, but the sign error is more fundamental because it undermines the theoretical derivation of the phantom phase. The exponential single-field model, by contrast, satisfies the sign assumption and remains internally consistent, so the framework has a valid core; however, the hyperbolic models and the phantom-phase conclusions drawn from them are not solutions of the equations as written. Since the reader's REJECT verdict is supported by this same concern, no adjustment to the verdict is needed.","tokens_in":18893,"tokens_out":7963,"duration_ms":74887,"concrete_test":"Symbolically substitute the hyperbolic solution φ(t)=2/B arctanh(e^{−AB²t}), its derivative, and H(t) from Eq. (36) into the original EOM (5) with sgn(φdot) evaluated as −1 (the actual sign of the solution). Compute the residual; if it does not vanish identically, the solution does not satisfy the model. Equivalently, verify that the relation φdot = −AB sinh(Bφ) is incompatible with the correctly signed first-order equation φdot = α − Wφ derived from Eq. (5) under sgn(φdot)=−1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After Eq. (7) the paper fixes sgn(φdot)>0. This assumption enters the first-order equation (21), φdot = -(Wφ + α), and the pressure (6b). For the hyperbolic model W = A cosh(Bφ) − αφ, Eq. (21) gives φdot = −AB sinh(Bφ), which is negative for the stated solution (34) on t>0. Thus the solution contradicts the assumption under which the framework was derived. Repeating the derivation with the correct sign, sgn(φdot) = −1, changes the first-order relation to φdot = α − Wφ; the advertised solution satisfies neither that relation nor the original second-order EOM (5) once sgn(φdot) is evaluated properly. Consequently the phantom phase shown in Fig. 2, and the analogous two-field hyperbolic cases in Secs. IVB and IVC, are not solutions of the model defined by Eqs. (4)–(6). The exponential single-field model respects the sign assumption, so the framework itself is not uniformly invalid; but the paper's headline phantom-phase results rest on this sign error. Independently, the AIC comparison in Sec. V misreads positive ΔAIC relative to ΛCDM as support for the cuscuton model: positive values mean the cuscuton model has larger AIC, so the abstract's 'strong support' claim is also not supported by the quoted numbers.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies a canonical scalar field plus a cuscuton-like term in a flat FLRW background. The authors introduce a first-order formalism H=W(φ), φ̇=-(W_φ+α) under the assumption sgn(φ̇)>0, and use it to construct exponential and hyperbolic potentials for single- and two-field models. They report analytical solutions, discuss transitions to ω=-1, and identify phantom phases for hyperbolic potentials. The paper then adopts an H(z) parametrization inspired by the exponential single-field solution and fits it to cosmic chronometers, Pantheon+, BAO, and geometric CMB data, using AIC to compare with ΛCDM.","tokens_in":19191,"tokens_out":9025,"duration_ms":79080,"significance":"The exponential single-field construction is self-consistent: the solution (27) has φ̇>0 as assumed, and Eq. (72) provides a simple two-parameter H(z) extension of ΛCDM that can be constrained by data. If the phantom-phase claims were valid, the paper would be a useful contribution. However, the hyperbolic-potential solutions violate the sign assumption at the root of the first-order reduction, and the AIC paragraph in Sec. V misstates what positive ΔAIC values mean. These are load-bearing problems for the main theoretical and statistical conclusions.","major_comments":[{"comment":"The first-order reduction is derived under sgn(φ̇)>0 (stated after Eq. (7)), giving Eq. (21) as φ̇=-(W_φ+α). For W=A cosh(Bφ)-αφ this yields φ̇=-AB sinh(Bφ), and the advertised solution (34) has φ̇(t)<0 for all t>0. Hence Eq. (34) does not solve Eq. (5) when the absolute value in the cuscuton term is evaluated with the correct sign, and it does not satisfy Eq. (21) under the assumption used to derive the framework. Repeating the derivation for sgn(φ̇)=-1 gives φ̇=α-W_φ=2α-AB sinh(Bφ), which is not Eq. (34). The phantom-phase curve in Fig. 2 therefore is not a solution of the model defined by Eqs. (4)-(6).","section":"§III.B, Eq. (34)"},{"comment":"In Sec. V the paper reports ΔAIC=0.4, 0.9, 1.6 for the cuscuton-like model with ΛCDM as the reference. Under the stated definition ΔAIC=AIC_i-AIC_min (with ΛCDM as the best model), these positive values mean the cuscuton-like model has a larger AIC than ΛCDM; they do not indicate 'strong support' for the cuscuton-like model. At most ΔAIC≤2 shows that the two models are comparable by the usual rule of thumb. The abstract's claim of strong support for the model is therefore not supported by the quoted values; the authors should either recompute the comparison with the cuscuton model as the reference and report the actual AIC values, or substantially soften the claim.","section":"§V, AIC paragraph after Table I"},{"comment":"The same sign inconsistency affects the two-field hyperbolic models. The φ solution used there is again (53a) or (63a), which has φ̇<0 for t>0, while the equations of motion (38a) assume sgn(φ̇)>0. Consequently the e.o.s. results in Figs. 5 and 8, including the phantom phase in Fig. 5, should not be presented as consequences of the cuscuton Lagrangian (37) until the sign branch is handled correctly. The paper needs either to solve the models with the proper sgn(φ̇) branch or to restrict the claim to models that satisfy the stated assumption.","section":"§IV.B and §IV.C, Eqs. (53)-(55) and (63)-(64)"}],"minor_comments":[{"comment":"The text repeatedly has 'FLR W' where 'FLRW' is intended; this typo should be corrected throughout.","section":"Abstract and §II"},{"comment":"There is a typo 'custucon-like' in the paragraph on future work; it should be 'cuscuton-like'.","section":"Final remarks"},{"comment":"The column header '2√8 α∗' is not defined in the table or caption; the authors should state the derived quantity being reported (for example, in terms of α/B and H0).","section":"Table I"},{"comment":"The denominator contains (1-e^{αBt}); consistency with (59a) and the surrounding algebra suggests this should be (1-e^{-αBt}). Please check the sign convention.","section":"Eq. (60)"}],"recommendation":"reject","confidential_remarks":"The sign error in the hyperbolic models is the central technical problem, and the AIC interpretation is also incorrect. I would be willing to reconsider a revised version that either proves the phantom behavior with the correct sgn(φ̇) branch or removes the hyperbolic claims and corrects the statistical interpretation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Rough read: the paper is a mixed bag. The genuinely new thing is the first-order Hamilton-Jacobi-style construction for a scalar with a cuscuton kinetic term, H=W(φ), φdot=-(Wφ+α), with the potential reconstructed from W. That is a clean extension of Bazeia et al., and it works for the exponential single-field case: the analytic solution, the H(z) parametrization in Eqs. (31)/(72), and the MCMC constraints are internally consistent. Credit where due: the stability discussion via Horndeski sound speeds is competently done, and the exponential model gives a useful analytic dark-energy parametrization that reduces to ΛCDM at z=0 while allowing w(z) to evolve.\n\nThe soft spots are real. The hyperbolic-potential section violates the paper's own sgn(φdot)>0 assumption. Eq. (34) has φdot<0 for t>0; Eq. (33) gives φdot=-AB sinh(Bφ), which is negative. The first-order formalism was derived assuming positive sign, so the phantom phase in Fig. 2 and the analogous two-field cases are not solutions of the Lagrangian defined in Eq. (4) once the absolute value is evaluated correctly. This is not parameter tuning; it is a sign inconsistency inside the model. The exponential model is unaffected, but the advertised phantom crossing is not established.\n\nSecond, the AIC section misreads the numbers. They take ΛCDM as reference and quote ΔAIC=0.4, 0.9, 1.6. Since ΔAIC is defined as AIC_model minus AIC_min, positive values mean the cuscuton model is slightly worse than the reference, not strongly supported. The abstract's 'strong support' claim is contradicted by their own table. This is fixable in wording but it is embarrassing.\n\nThird, the two-field sections inherit the same sign problem whenever the hyperbolic field is used. The exponential+exponential case is fine; the hyperbolic cases are not.\n\nBottom line: the exponential single-field model and the first-order tool are worth having, and the observational constraints on that model are reasonable. The hyperbolic/phantom claims should not be published as they stand. A referee should ask for the sign issue to be fixed or the phantom sections removed, and for the AIC interpretation to be corrected. I would send this to review rather than desk-reject: the framework is novel enough and the exponential part is sound, but revision needs to be substantial.","headline":"A genuinely useful first-order framework for cuscuton dark energy, undermined by a sign inconsistency in the hyperbolic phantom solutions and an AIC misreading that overstates the data support.","tokens_in":19757,"tokens_out":8687,"would_cite":false,"duration_ms":69917,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k","95.36.+x"],"model":"deepseek-v4-flash","headline":"Adding a cuscuton-like term to quintessence dark energy makes the equations first-order and solvable, lets the equation of state cross below −1 without ghosts, and yields AIC support over ΛCDM.","keywords":["cuscuton","dark energy","first-order formalism","equation of state","phantom phase","FLRW cosmology","Monte Carlo Markov Chain","Akaike Information Criterion"],"falsifier":"Take the hyperbolic solution $\\phi(t)=\\frac{2}{B}\\operatorname{arctanh}(e^{-AB^{2}t})$, note that $\\dot{\\phi}<0$, and insert it into the original equation of motion with the cuscuton term evaluated as $|\\dot{\\phi}|$ rather than $+\\dot{\\phi}$; the equation is not satisfied. A direct numerical integration of Eq. (5) with the absolute value for the hyperbolic potential, or an independent check of whether $\\omega<-1$ persists with the correct sign, would settle the central phantom-phase claim.","tokens_in":18678,"feed_emoji":"🌌","tokens_out":6823,"duration_ms":55024,"temperature":0.7,"pith_summary":"This paper argues that adding a cuscuton-like term, $\\alpha\\sqrt{|\\partial_\\mu\\phi\\,\\partial^\\mu\\phi|}$, to the standard kinetic term of a quintessence scalar field changes the pressure without changing the energy density, and that this minimal modification is enough to let dark energy evolve from a matter- or radiation-dominated phase to a cosmological-constant-like phase. Under the assumption that the field velocity is positive, the equations of motion reduce to a first-order system, $H=W(\\phi)$ and $\\dot{\\phi}=-(W_\\phi+\\alpha)$, for which the authors obtain analytic solutions for exponential and hyperbolic potentials. For the hyperbolic potential the solutions cross the phantom divide, $\\omega<-1$, and the paper shows that the model has positive scalar and tensor sound speeds, so no ghosts or gradient instabilities appear. The paper also fits the single-field exponential version to geometrical cosmological data and, using the Akaike Information Criterion, reports strong support for the cuscuton-like model over $\\Lambda$CDM when supernova, BAO, and CMB data are included.","feed_headline":"Cuscuton term steers dark energy through a phantom phase","feed_subtitle":"A tiny added term makes quintessence analytically solvable and can cross ω = −1 without ghosts.","key_machinery":"The load-bearing object is the first-order pair $H=W(\\phi)$ and $\\dot{\\phi}=-(W_\\phi+\\alpha)$, together with the potential $V=\\tfrac{3}{2}W^{2}-\\tfrac{1}{2}(W_\\phi+\\alpha)^{2}$ derived from the Friedmann equation. It converts the second-order equation of motion into a quadrature, so each choice of $W$ gives an analytical cosmic history; the cuscuton constant $\\alpha$ enters additively in the field velocity and in the pressure but not in the energy density, which is what allows $\\omega$ to drop below $-1$. The stability argument is carried by the identification of the Lagrangian with a subclass of Horndeski theory, which yields $c_s^{2}=1+\\alpha/\\sqrt{2X}>0$ and $c_{\\rm GW}^{2}=1$ for one field, and $v^{2}=0$ for the second field in the two-field case, ruling out ghosts and gradient instabilities in the models considered.","core_discovery":"The central claim is that the cuscuton-like addition to the scalar-field dark energy Lagrangian supports a first-order framework analogous to the Hamilton-Jacobi formalism: taking the Hubble parameter as $H=W(\\phi)$ makes the field equation $\\dot{\\phi}=-(W_\\phi+\\alpha)$, with the potential fixed by $V=\\tfrac{3}{2}W^{2}-\\tfrac{1}{2}(W_\\phi+\\alpha)^{2}$. This turns the search for cosmic histories into the choice of $W$, and the paper constructs models where the equation-of-state parameter $\\omega$ evolves from the radiation or matter value toward $-1$. In the single-field hyperbolic model, and in two-field models with a hyperbolic component, the evolution crosses into the phantom regime $\\omega<-1$ before settling at the cosmological constant; the paper verifies that the scalar and tensor propagation speeds remain positive and that the non-dynamical cuscuton behavior is recovered in the two-field case, so the phantom phase is not accompanied by ghosts. Using the single-field exponential model parametrized by $\\rho_{\\rm CL}$ in terms of redshift, the paper constrains the extra parameter $B$ and the derived cuscuton parameter $\\alpha$ with a background-only Monte Carlo Markov Chain analysis of geometrical probes, and finds AIC differences of $0.4$–$1.6$ relative to $\\Lambda$CDM for datasets that include supernovae, BAO, and CMB, which it interprets as strong support.","pith_inferences":["The paper fixes $\\operatorname{sgn}(\\dot{\\phi})>0$ when writing the cuscuton term, but its hyperbolic-potential solution $\\phi(t)=\\frac{2}{B}\\operatorname{arctanh}(e^{-AB^{2}t})$ has $\\dot{\\phi}<0$; evaluating the absolute value with the correct sign would reverse the cuscuton contribution in Eq. (5), so the phantom-phase solutions appear to be solutions of a different equation than the one stated","The AIC comparison uses only background evolution, as the paper itself acknowledges; a full CMB likelihood and perturbation evolution, including the matter sound speed and growth of structure, could change the statistical verdict.","A natural testable extension is to compute the growth rate $f\\sigma_8$ in the cuscuton-like model: because the second-field sound speed is $v^{2}=0$, the model predicts a specific non-standard clustering signature that geometric-only data cannot see.","If the sign issue is repaired (for example, by allowing $\\operatorname{sgn}(\\dot{\\phi})$ to become time-dependent or by choosing a formulation that makes the sign consistent), the first-order construction itself is a useful template for generating analytic late-time cosmologies with non-canonical kinetic terms."],"forward_implications":["In the single-field exponential model, $\\omega$ becomes time-dependent and asymptotes to $-1$, so a scalar field that would have constant equation of state in standard dynamics can describe the transition from deceleration to acceleration.","In the hyperbolic models, the cuscuton term produces a phantom phase ($\\omega<-1$) before the cosmological constant phase, and the stability analysis implies this phase has no ghost or gradient instabilities.","In two-field models, the cuscuton parameter $\\alpha$ controls the duration of transitions between phases in the exponential case, and the amplitude of fluctuations around the initial phase plus the presence of a phantom phase in the hyperbolic case.","The observational fit yields values of the matter density and $H_0$ compatible with $\\Lambda$CDM within $1\\sigma$, while $\\omega$ excludes $-1$ at $1\\sigma$ when only cosmic chronometers are used; extended datasets restore $\\omega\\approx-1$.","AIC comparison gives $\\Delta\\rm AIC=0.4$–$1.6$ for datasets including supernovae, BAO, and CMB, which the paper reads as strong support for the cuscuton-like model over $\\Lambda$CDM, although with chronometers alone the model is moderately disfavored ($\\Delta\\rm AIC=4.1$)."],"supporting_citations":[{"why":"Introduces the cuscuton as a dark energy model and supplies the non-dynamical term whose modification the paper studies.","marker":"[24, 25]"},{"why":"Provides the Horndeski-theory stability conditions and sound-speed formulas used to show the absence of ghosts.","marker":"[44]"},{"why":"Establishes the first-order formalism for dark energy that the paper adapts with the cuscuton term.","marker":"[45]"},{"why":"Gives the bent-brane first-order construction with a hyperbolic potential that motivates the hyperbolic model.","marker":"[58]"},{"why":"Supplies the fast background-only Monte Carlo Markov Chain sampler and likelihoods used for the parameter constraints.","marker":"[66, 67]"},{"why":"Provides the cosmic chronometer Hubble data used in the observational fits.","marker":"[69]"},{"why":"Supplies the type Ia supernova compilation used in the distance constraints.","marker":"[70]"},{"why":"Provide the baryon acoustic oscillation measurements used in the combined datasets.","marker":"[71-76]"},{"why":"Supplies the geometric CMB distance-ratio likelihood included in the full dataset.","marker":"[77]"},{"why":"Defines the Akaike Information Criterion used for the model comparison with $\\Lambda$CDM.","marker":"[78]"}],"fun_headline_variants":["Cuscuton term lets dark energy cross phantom divide cleanly","Analytical dark energy with cuscuton: phantom phase without ghosts","Cuscuton dark energy: data favor it over ΛCDM","New dark energy model crosses ω=-1, no ghosts, fits data","Cuscuton term opens analytical window on phantom crossing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the unvarying sign choice $\\operatorname{sgn}(\\dot{\\phi})>0$ made when the cuscuton term is written as $+\\alpha\\dot{\\phi}$, because the hyperbolic-potential solutions that produce the phantom phase have $\\dot{\\phi}<0$ and therefore do not solve the original absolute-value equation of motion unless that sign choice is silently abandoned.","fun_headline_variants_meta":{"raw":{"variants":["Cuscuton term lets dark energy cross phantom divide cleanly","Analytical dark energy with cuscuton: phantom phase without ghosts","Cuscuton dark energy: data favor it over ΛCDM","New dark energy model crosses ω=-1, no ghosts, fits data","Cuscuton term opens analytical window on phantom crossing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000307,"raw_usage":{"total_tokens":1784,"prompt_tokens":996,"completion_tokens":788,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":696}},"tokens_in":612,"tokens_out":788,"duration_ms":6845,"temperature":1.0,"reasoning_tokens":696,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T14:49:09.811119+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the hyperbolic solution $\\phi(t)=\\frac{2}{B}\\operatorname{arctanh}(e^{-AB^{2}t})$, note that $\\dot{\\phi}<0$, and insert it into the original equation of motion with the cuscuton term evaluated as $|\\dot{\\phi}|$ rather than $+\\dot{\\phi}$; the equation is not satisfied. A direct numerical integration of Eq. (5) with the absolute value for the hyperbolic potential, or an independent check of whether $\\omega<-1$ persists with the correct sign, would settle the central phantom-phase claim.","supporting_citations":[{"cited_title":"Zhong, Fei-Yu Li, Xu-Dong Liu, K-field kinks in two- dimensional dilaton gravity , Phys","cited_arxiv_id":null,"evidence_quote":"Provides the Horndeski-theory stability conditions and sound-speed formulas used to show the absence of ghosts."}],"review_version":1}