{"id":"cd19c241-57c4-46c0-b7d4-d79449208f76","arxiv_id":"2412.05577","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"Derives a non-flat interacting Barrow holographic dark energy model and fits it to CC and Pantheon data, but the statistical claim of preference is contradicted by the paper's own AIC/BIC/DIC table.","lead":"This paper extends Barrow holographic dark energy to curved universes with dark matter-dark energy interactions, then fits it to supernova and cosmic clock data. It concludes the data prefer a slightly curved, interacting universe, but the paper's own model-comparison table actually favors the standard ΛCDM model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Model-selection claim rests on a sign misreading of Table 3; the paper's own AIC/BIC/DIC values favor ΛCDM for CC and mostly for CC+Pantheon, so 'non-flat interacting preferred' is unsupported.","rationale":"The reader's rejection is correct, and I reach the same verdict, though through the model-comparison table rather than the future-event-horizon concern identified as the reader's weakest assumption. The central claim is a statistical preference claim: 'a non-flat interacting scenario is preferred by observational data'. Preference claims of this kind are decided by the information criteria reported in Table 3. The paper's own table, read with the standard sign convention where lower AIC/BIC/DIC means the better model, does not support the abstract. For CC, the model has consistently higher AIC/BIC/DIC than ΛCDM, meaning ΛCDM fits better, despite the text claiming the opposite. For Pantheon, the model has slightly lower AIC/BIC/DIC, meaning the model fits slightly better, despite the text claiming ΛCDM is preferred. For CC+Pantheon, ΛCDM generally remains favored or comparable. The additional argument from 'nonzero' best-fit values of Γ and Ω_k0 would require confidence intervals or a posterior probability that excludes zero; neither is reported. The horizon-integral convergence issue is a legitimate secondary concern, but it is not needed to reject the headline claim, because the statistical comparison is already internally inconsistent. I therefore keep the reader's REJECT verdict unchanged, while noting that the load-bearing failure is the misinterpretation of the model-selection results in Section 6.1 and the Conclusions.","tokens_in":26274,"tokens_out":4759,"duration_ms":45572,"concrete_test":"Reproduce Table 3 from the quoted χ²_min values and sample sizes using ΔIC = IC_model − IC_ΛCDM (not the absolute value), and explicitly mark which model has the lower IC in each row; if the corrected signs show ΛCDM is preferred or statistically indistinguishable in the CC and CC+Pantheon rows, then the abstract's 'non-flat interacting scenario is preferred' is contradicted by the paper's own model-selection values.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's claim that a non-flat interacting scenario is preferred by observational data is not supported by the paper's own Section 6 model comparison. Section 6.1 defines ΔAIC = |AIC_model − AIC_ΛCDM| and states that ΔIC ≤ 2 means the model is 'strongly favoured', but Table 3 lists signed differences and the text interprets the signs in the wrong direction. For CC, Case I with k=+1 has AIC_model = 51.74 versus AIC_ΛCDM = 50.10, i.e. ΔAIC = +1.64; the model has the larger AIC, so ΛCDM is preferred, yet the text says the proposed model is 'strongly favored'. For Pantheon, Case I with k=+1 has AIC_model = 38.38 versus AIC_ΛCDM = 39.77, i.e. ΔAIC = −1.39; the proposed model has the lower AIC, yet the text says 'negative indicates that ΛCDM model is preferred'. Under the corrected sign convention, the CC and CC+Pantheon rows of Table 3 mostly leave ΛCDM with the lower or comparable AIC/BIC/DIC. Moreover, Tables 1 and 2 report best-fit Γ and Ω_k0 as 'nonzero' without confidence intervals, so no statistical significance for the nonzero values is established. The observational preference claim is therefore contradicted by the paper's own numbers rather than supported by them.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends Barrow holographic dark energy (BHDE) to an interacting, non-flat FRW background. For closed and open spatial geometries and three phenomenological interaction terms Q = -ΓHρ_DE, Q = -ΓHρ_m, and Q = -ΓH(ρ_m+ρ_DE), the authors derive coupled evolution equations for Ω_DE and Ω_m and analytic expressions for the dark-energy equation of state w_DE, using the future event horizon as the holographic cutoff. They then fit the parameters Γ, Ω_k0, and Δ to cosmic chronometer (CC) and Pantheon data via MCMC and compare the model with ΛCDM using AIC, BIC, and DIC. The abstract claims that nonzero interaction strength and curvature contribution imply that a non-flat interacting scenario is preferred by observational data, while Section 7 states statistical compatibility with ΛCDM for CC and mild tension for CC+Pantheon.","tokens_in":26621,"tokens_out":8059,"duration_ms":76005,"significance":"The analytic part of the paper is a useful extension of earlier BHDE studies: it generalizes the flat interacting BHDE model of Sheykhi et al. and the non-flat non-interacting model of Adhikary et al., and the reported limits Γ=0 and Δ=0 reduce to known equations. However, the observational significance is not as claimed. Read with the correct sign convention, the paper's own Table 3 shows that ΛCDM has lower or comparable AIC/BIC/DIC for most CC and CC+Pantheon comparisons, so the central preference claim is unsupported. The 'nonzero' best-fit values in Tables 1 and 2 are reported without confidence intervals, and several contours in Figures 5-7 include zero. The paper therefore does not currently establish its headline result; its value is primarily in the derivation and, at best, in a set of largely null observational constraints.","major_comments":[{"comment":"The model-selection interpretation is based on a sign error. The paper defines ΔAIC = |AIC_model − AIC_ΛCDM| and states that ΔIC ≤ 2 means the model is strongly favored, but Table 3 lists signed differences and the text reads the signs in the wrong direction. For CC, Case I, k=+1, AIC_model = 51.74 and AIC_ΛCDM = 50.10, so the model is worse by 1.64, yet the text says the proposed model is strongly favored. For Pantheon, Case I, k=+1, AIC_model = 38.38 and AIC_ΛCDM = 39.77, so the model is better by 1.39, yet the text says negative Δ indicates ΛCDM is preferred. Correcting the signs shows that for the CC and CC+Pantheon rows ΛCDM is generally preferred or indistinguishable, directly contradicting the abstract's claim that a non-flat interacting scenario is preferred by the data.","section":"Section 6.1, Table 3"},{"comment":"The conclusion that Γ and Ω_k0 are 'nonzero' is not statistically established. Tables 1 and 2 list only best-fit values without confidence intervals, and the contours in Figures 5-7 include Γ = 0 for several cases (for example, CC Case I in Figure 5 has a Γ axis extending from -0.2 to 0.2). Additionally, the horizontal axes of Figure 5 are labeled with Ω_k0 ranging from about 0.0290 to 0.0310 for Cases I and II, whereas Table 1 lists best-fit values around 0.0098; this apparent inconsistency makes the displayed constraints unreliable as reported.","section":"Section 6, Tables 1-2 and Figures 5-7"},{"comment":"The numerical solution of the evolution equations is not sufficiently documented. Equations (3.16)-(3.17) and (4.14)-(4.15) are integro-differential because y is defined by the improper future-event-horizon integrals in equations (3.3) and (4.3). The text states only that the equations are 'solved numerically imposing initial conditions,' without specifying the integrator, step size, convergence criteria, or how the improper integrals are evaluated. Since every fitted parameter and all AIC/BIC/DIC values depend on these solutions, the lack of numerical details is a load-bearing reproducibility gap.","section":"Sections 5.1 and 6"},{"comment":"The model-selection comparison is ambiguous about which parameter sets are being compared. Table 1 fixes Δ = 0.1 and fits (Γ, Ω_k0); Table 2 fixes Ω_k0 and fits (Γ, Δ); the text says the comparison uses the Γ-Δ case. The number of effective parameters used in the AIC/BIC/DIC calculation for the proposed model and for the ΛCDM baseline is not stated, and the absolute-value definition of ΔAIC is inconsistent with the signed entries in Table 3. This makes it difficult to verify the reported information criteria.","section":"Section 6.1, Tables 1-3"}],"minor_comments":[{"comment":"The abstract's statement that a non-flat interacting scenario is preferred should be removed or rewritten to reflect the corrected reading of Table 3; Section 7's compatibility and tension statements should also be rephrased accordingly.","section":"Abstract and Section 7"},{"comment":"The text says that for Δ = 0.1 the phantom divide is crossed at z ~ 0.5 and that for Δ = 0.2 it is crossed 'much earlier (at around z ~ 0.5)'; this is internally contradictory and should be replaced with the actual crossing redshifts read from the figures.","section":"Section 5, Figure 1 discussion"},{"comment":"The x-axis values in Figure 5 appear to be inconsistent with the best-fit values in Table 1; either the axes are mislabeled or the table is incorrect, and this needs to be fixed.","section":"Figure 5"},{"comment":"The MCMC analysis is described only by a reference; the priors, proposal scheme, burn-in length, and convergence diagnostics should be specified for reproducibility.","section":"Section 6"},{"comment":"There are several typographical and stylistic issues, including 'negative negative curvature' in the paragraph preceding Figure 11, 'T able 1' in the caption, and incomplete reference entries for [77] and [79]; these should be corrected.","section":"Various"}],"recommendation":"reject","confidential_remarks":"The paper's advertised observational conclusion is reversed by its own model-selection table, and the absence of confidence intervals on the 'nonzero' parameters removes the remaining statistical support for the headline claim. A revised version could potentially be reframed as a derivation-and-constraints paper, but that would require substantial rewriting and re-analysis, and the editor may wish to consider whether such a version fits the journal's interests."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe useful part of this paper is the first half. Adhikary and Das construct the coupled evolution equations for Ω_DE and Ω_m for Barrow holographic dark energy in closed and open FRW backgrounds, with a future event horizon cutoff and three phenomenological interaction terms. The algebra is internally consistent, the Γ=0 and Δ=0 limits reduce to the known non-flat HDE equations, and the combination of Barrow entropy + curvature + interaction is genuinely new relative to Sheykhi et al. (flat) and to the authors' own prior non-flat, non-interacting paper. That is real, reproducible derivational work and I would not want to lose it.\n\nThe soft spot is the observational section, and it is not minor. The abstract says the data prefer a non-flat interacting scenario. Table 3 says the opposite, once you read the signs correctly. The authors define ΔAIC, ΔBIC, ΔDIC as absolute differences in Section 6.1, but then interpret signed values as if they were absolute. For CC data, the model has larger AIC than ΛCDM in every case (e.g., 51.74 vs 50.10 for Case I, k=+1), so ΛCDM is preferred, yet the text calls the model 'strongly favoured.' For Pantheon, the model has lower AIC (38.38 vs 39.77) and the text says 'negative indicates that ΛCDM model is preferred' — the exact opposite of what lower AIC means. Once the signs are corrected, the model only beats ΛCDM on Pantheon alone, which is too weak to justify two extra parameters. The CC+Pantheon comparison mostly favors ΛCDM.\n\nThe 'nonzero' Γ and Ω_k0 claim is also unsupported, because the paper reports best-fit points with no confidence intervals. A point estimate different from zero is not evidence.\n\nThe conclusions section is more careful — it says 'statistically compatible' for CC and 'mild tension' for CC+Pantheon — but the abstract and the model-comparison text contradict the paper's own numbers. That is a load-bearing flaw, not a cosmetic one. The fix is straightforward: correct the sign convention, report uncertainties on Γ and Ω_k0, and rewrite the abstract to match what Table 3 actually shows.\n\nWho reads this? Someone working on holographic dark energy phenomenology who wants the coupled equations in curved space. The derivation is worth a referee's time. I would not desk-reject it, but I would send it back for major revision with a clear request to redo the statistical interpretation. My own verdict: the paper as written does not support its headline claim, but the core derivation is sound.","headline":"The derivation of interacting Barrow HDE in curved space is solid, but the abstract's observational-preference claim contradicts the paper's own Table 3 and the sign handling in the model comparison is backwards.","tokens_in":27178,"tokens_out":3793,"would_cite":false,"duration_ms":34379,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.36.+x","98.80.-k"],"model":"deepseek-v4-flash","headline":"Barrow holographic dark energy with both dark-sector interaction and spatial curvature fits the CC+Pantheon data, with $\\Gamma$ and $\\Omega_{k0}$ best-fit to nonzero values.","keywords":["Barrow holographic dark energy","interacting dark sectors","non-flat universe","future event horizon","equation of state","cosmic chronometers","Pantheon supernovae","phantom crossing"],"falsifier":"A decisive test is to repeat the MCMC fit with a different, equally motivated holographic cutoff—such as the Hubble horizon or a generalized entropy-based cutoff—and check whether the best-fit $\\Gamma$ and $\\Omega_{k0}$ remain nonzero; if they return to zero, the claimed preference for a non-flat interacting scenario is an artifact of the horizon choice. A second check is to compare the model's predicted CMB acoustic-peak positions with Planck data, since a nonzero $\\Omega_{k0}$ of order 0.01 would shift those peaks.","tokens_in":26019,"feed_emoji":"🔭","tokens_out":11487,"duration_ms":97879,"temperature":0.7,"pith_summary":"This paper tries to establish that Barrow holographic dark energy—dark energy whose density follows from a quantum-gravitational modification of the entropy-area relation—is preferred by cosmic chronometer and Pantheon data when the dark sectors interact and the spatial geometry is non-flat. It derives coupled evolution equations for the dark-energy and matter density parameters in closed and open FRW universes, for three interaction terms $Q = -\\Gamma H\\rho_{\\mathrm{DE}}$, $Q = -\\Gamma H\\rho_m$, and $Q = -\\Gamma H(\\rho_m + \\rho_{\\mathrm{DE}})$, along with analytic expressions for the equation-of-state parameter $w_{\\mathrm{DE}}$. The reported best-fit values of the interaction strength $\\Gamma$ and curvature density $\\Omega_{k0}$ are small but nonzero, which the authors read as evidence for a non-flat interacting scenario. The paper also reports that the model is statistically compatible with ΛCDM on chronometer data and mildly tensioned on the combined CC+Pantheon data.","feed_headline":"Cosmic data favor curved Barrow dark energy with interaction","feed_subtitle":"Best fits put both curvature and the dark-sector coupling away from zero, challenging flat ΛCDM.","key_machinery":"Barrow entropy $S_B \\propto (A/A_0)^{1+\\Delta/2}$, with $0 \\le \\Delta \\le 1$, changes the holographic dark-energy density to $\\rho_{\\mathrm{DE}} = C L^{\\Delta-2}$, where $L$ is the horizon length. For curved space the paper takes the future event horizon with $L = a\\sin y$ (closed) and $L = a\\sinh y$ (open), $y = R_h/a = \\int_x^\\infty dx/(aH)$; this integral ties the local density evolution to the whole future history of the scale factor. The coupled equations for $\\Omega_{\\mathrm{DE}}$ and $\\Omega_m$ and the closed-form $w_{\\mathrm{DE}}$ for each interaction term are derived from the continuity equations $\\dot{\\rho}_m + 3H\\rho_m = Q$ and $\\dot{\\rho}_{\\mathrm{DE}} + 3H(1+w_{\\mathrm{DE}})\\rho_{\\mathrm{DE}} = -Q$, with the three phenomenological choices of $Q$. The model's observational signature is carried by three parameters: the Barrow exponent $\\Delta$, the interaction strength $\\Gamma$, and the curvature density $\\Omega_{k0}$.","core_discovery":"On the paper's own terms, the central discovery is that the future-event-horizon Barrow holographic dark energy model, extended to non-flat geometry and supplemented with a linear dark-sector interaction, remains observationally viable and can cross the phantom divide at late times. For each of the three interaction choices, the evolution of $\\Omega_{\\mathrm{DE}}$ and $\\Omega_m$ is governed by a coupled pair of first-order equations whose curvature dependence enters through $y = \\int_x^\\infty dx/(aH)$ via $\\sin y$ (closed) or $\\sinh y$ (open) in the horizon length $L = a\\sin y$ or $L = a\\sinh y$. The MCMC fits to the CC and Pantheon data give best-fit $\\Gamma$ and $\\Omega_{k0}$ that are nonzero, and the authors interpret this as indicating that a non-flat interacting scenario is preferred. Their AIC/BIC/DIC comparison shows the model is statistically compatible with ΛCDM for the CC data, while the combined dataset produces mild tension for one of the three interaction forms.","pith_inferences":["The 'preferred by data' statement is based on nonzero best-fit values rather than a decisive model-selection statistic; the same paper's AIC/BIC/DIC table shows ΛCDM preferred for Pantheon alone, so the evidential weight is dataset-dependent.","Because $\\Omega_{m0}$ and $\\Omega_{\\mathrm{DE}0}$ are fixed rather than marginalized, the inferred nonzero $\\Gamma$ and $\\Omega_{k0}$ are conditional on those priors; leaving matter density free could widen the contours and shift the best fit.","The curvature preference is of order $|\\Omega_{k0}| \\sim 0.01$, which future measurements of the CMB acoustic scale and BAO angular diameter distance should detect or exclude independently of this model.","The fact that all three interaction forms are simple $Q \\propto H\\rho$ choices means the nonzero $\\Gamma$ could be a proxy for more complicated dark-energy physics; a field-theoretically derived $Q$ would clarify whether the interaction is physical."],"forward_implications":["The dark-energy equation of state $w_{\\mathrm{DE}}$ can cross the phantom divide for nonzero $\\Delta$, so an interaction can generate phantom-like behavior even though each component begins in the quintessence regime.","For fixed $\\Delta$, increasing $\\Gamma$ shifts the phantom crossing redshift and pushes $w_{\\mathrm{DE}}$ deeper into the phantom region in both closed and open geometries.","Closed models enter the phantom regime earlier than open models at equal $\\Delta$ and $\\Gamma$, giving curvature a potentially observable effect on the timing of acceleration.","If the best-fit nonzero $\\Gamma$ and $\\Omega_{k0}$ are real, the universe's energy flow runs from matter to dark energy and the spatial geometry is slightly closed or open, both of which are absent in flat ΛCDM.","Relative to ΛCDM, the model is statistically compatible for CC data but mildly tensioned for CC+Pantheon in the interaction form proportional to total density, so the current H(z)+SNe data cannot decisively choose between them."],"supporting_citations":[{"why":"Defines the Barrow entropy-area relation that motivates the modified dark-energy density used throughout the paper.","marker":"[31]"},{"why":"Establishes the Barrow holographic dark energy model and the form $\\rho_{\\mathrm{DE}} = C L^{\\Delta-2}$.","marker":"[32]"},{"why":"Supplies the holographic dark energy framework with the future event horizon as the standard cutoff.","marker":"[19]"},{"why":"Provides the earlier non-flat Barrow holographic dark energy evolution equations that this paper extends to the interacting case.","marker":"[45]"},{"why":"Reports an earlier interacting Barrow holographic dark energy study with future event horizon but in a spatially flat background, the gap this paper fills.","marker":"[55]"},{"why":"Cites evidence that a nonzero curvature universe is favored, motivating the non-flat analysis.","marker":"[56]"},{"why":"Gives the curved-space future event horizon relations $L = a\\sin y$ and $L = a\\sinh y$ used in the derivation.","marker":"[66]"},{"why":"Provides the 57-point cosmic chronometer H(z) dataset used for the CC likelihood.","marker":"[77]"},{"why":"Provides the 1048-point Pantheon supernova dataset used for the SNIa likelihood.","marker":"[78]"}],"fun_headline_variants":["Non-flat Barrow dark energy fits cosmic data with interaction","Curved universe plus coupling: Barrow DE favored by data","Barrow DE: curvature and interaction nonzero in best fits","Interacting Barrow DE crosses phantom in non-flat universe"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the future event horizon is the correct holographic cutoff, in the curved-space form $L = a\\sin y$ (closed) or $L = a\\sinh y$ (open) with $y = \\int_x^\\infty dx/(aH)$, and that the resulting coupled integro-differential system has a unique numerical solution over the fitted redshift range; if this horizon choice is wrong or the integrals do not converge for the fitted parameters, the derived $w_{\\mathrm{DE}}$ and constraints change.","fun_headline_variants_meta":{"raw":{"variants":["Non-flat Barrow dark energy fits cosmic data with interaction","Curved universe plus coupling: Barrow DE favored by data","Barrow DE: curvature and interaction nonzero in best fits","Interacting Barrow DE crosses phantom in non-flat universe"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000613,"raw_usage":{"total_tokens":2827,"prompt_tokens":899,"completion_tokens":1928,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":1858}},"tokens_in":515,"tokens_out":1928,"duration_ms":15417,"temperature":1.0,"reasoning_tokens":1858,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:35:02.241646+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is to repeat the MCMC fit with a different, equally motivated holographic cutoff—such as the Hubble horizon or a generalized entropy-based cutoff—and check whether the best-fit $\\Gamma$ and $\\Omega_{k0}$ remain nonzero; if they return to zero, the claimed preference for a non-flat interacting scenario is an artifact of the horizon choice. A second check is to compare the model's predicted CMB acoustic-peak positions with Planck data, since a nonzero $\\Omega_{k0}$ of order 0.01 would shift those peaks.","supporting_citations":[{"cited_title":"Barrow, The area of a rough black hole , Physics Letters B 808 (2020) 135643","cited_arxiv_id":null,"evidence_quote":"Defines the Barrow entropy-area relation that motivates the modified dark-energy density used throughout the paper."},{"cited_title":"Saridakis, Barrow holographic dark energy , Phys","cited_arxiv_id":null,"evidence_quote":"Establishes the Barrow holographic dark energy model and the form $\\rho_{\\mathrm{DE}} = C L^{\\Delta-2}$."},{"cited_title":"Adhikary, S","cited_arxiv_id":null,"evidence_quote":"Provides the earlier non-flat Barrow holographic dark energy evolution equations that this paper extends to the interacting case."},{"cited_title":"Huang and M","cited_arxiv_id":null,"evidence_quote":"Gives the curved-space future event horizon relations $L = a\\sin y$ and $L = a\\sinh y$ used in the derivation."},{"cited_title":"Mhamdi et al., Observational constraints on the growth index parameters in f (q) gravity, 2024","cited_arxiv_id":null,"evidence_quote":"Provides the 57-point cosmic chronometer H(z) dataset used for the CC likelihood."},{"cited_title":"Scolnic, D","cited_arxiv_id":null,"evidence_quote":"Provides the 1048-point Pantheon supernova dataset used for the SNIa likelihood."}],"review_version":1}