{"id":"ccc453cf-1696-48dd-b881-d9124f6b4d26","arxiv_id":"2506.16004","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"A fractional holographic dark energy model in an anisotropic Bianchi I background reduces to the flat wCDM expansion when fit to current data.","lead":"This paper fits a fractional holographic dark energy model to Hubble, BAO, and supernova data and finds parameter values close to the standard cosmological model. The Bianchi anisotropy drops out of the fitted equations, so the practical result is a standard wCDM fit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fitted H(z) is standard wCDM with no n dependence; the LRS Bianchi pressure equations are not satisfied, so the data never test the proposed anisotropic FHDE model.","rationale":"The reader's weakest assumption identifies the same load-bearing concern. Re-deriving Eqs. (13)-(25) confirms that Eq. (24) follows from Eq. (13) plus the two continuity equations, without using the pressure equations (14) and (15). Inserting Eq. (24) into either pressure equation reduces, for n>2, to omega_d Omega_d = 1, which is incompatible with the fitted values (omega_d ~ -0.78, Omega_d ~ 0.73, product ~ -0.57) and would imply q = 2. The spatial part of the Einstein equations is therefore not satisfied by the model used in the fit, so the Bianchi structure is effectively dropped. I also confirmed that Eq. (29) is normalized to wCDM with n cancelling; the observational constraints cannot distinguish the claimed anisotropic model from a flat isotropic wCDM fit. The reader's additional point that alpha = 1.01 lies outside the 1-sigma range is less secure, since 1.01 sits at the upper 1-sigma edge of the reported errors in all rows of Table 1, so I do not rest the rejection on that point. The constant-versus-varying omega_d issue is real but secondary to the field-equation inconsistency. These findings support the REJECT verdict, so no adjustment is needed.","tokens_in":15068,"tokens_out":13595,"duration_ms":136381,"concrete_test":"Take the reported best-fit values, e.g. alpha = 0.87 and Omega_d0 = 0.73, and substitute Eq. (24), equivalently Hdot/H^2 = -(3/2)(1+omega_d Omega_d), into both pressure equations (14) and (15). For n>2, any solution of the Bianchi field equations must satisfy omega_d Omega_d = 1, hence Hdot/H^2 = -3 and q = 2. Evaluate Eq. (25) at the same parameters, giving Hdot/H^2 = -3alpha(1+omega_d)/(3alpha-2) ~ -0.54. If the two values differ, the fitted H(z) is not a solution of the Bianchi Einstein equations. A decisive numerical version is to refit the data using the full system (13), (16), (17), (20) together with (14) and (15), treating the fractional-HDE parameter c^2 and n as free; if the resulting constraints differ materially from Table 1 or no acceptable solution exists for n>2, the central viability claim is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central viability claim is carried by Eq. (29), H(z)=H0[Omega_m0(1+z)^3+(1-Omega_m0)(1+z)^{3(1+omega_d)}]^{1/2}. This is the standard flat wCDM Hubble law: the Bianchi index n cancels and the fractional-HDE density (20) is replaced by a constant-EoS power law (27). The MCMC constraints in Table 1 therefore do not test the LRS Bianchi geometry or the fractional density; any isotropic wCDM fit would give the same H(z). More seriously, the assumed Bianchi model is not a solution of its own field equations. The derivation obtains Eq. (24) from Eq. (13) and the conservation equations, but never enforces the pressure equations (14) and (15). Substituting Eq. (24) into either (14) or (15) gives omega_d Omega_d = 1 for n>2, so Eq. (24) forces Hdot/H^2 = -3 and q = 2. Equation (25) with the fitted values (alpha ~ 0.87, Omega_d0 ~ 0.73) gives Hdot/H^2 ~ -0.54. The fitted expansion history cannot satisfy the Bianchi Einstein equations. The fit additionally treats omega_d as constant in Eq. (27), while the model's own Eq. (26) is a function of Omega_d(z). Thus the constrained object is a wCDM proxy, not the fractional HDE model, and the Lambda-CDM-like plots at alpha = 1.01 do not establish viability of the proposed anisotropic fractional model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a locally rotationally symmetric (LRS) Bianchi type-I cosmological model filled with fractional holographic dark energy, derives a Hubble parameter H(z), and fits it to Hubble, BAO, and Pantheon+SH0ES data using MCMC. The authors report constraints on H0, the fractional parameter α, and Ω_m0, and use these to study the density parameters, deceleration parameter, equation of state, and statefinder diagnostics. The central claim is that the model is observationally viable and closely mimics ΛCDM, with a deceleration transition at z_t ≈ 0.55 and a statefinder pair (r,s) near (0.74,0.07). My assessment is that the central derivation is inconsistent with the paper's own field equations and that the fitted Hubble law is a standard wCDM proxy rather than the proposed Bianchi fractional-HDE model.","tokens_in":15445,"tokens_out":10509,"duration_ms":106489,"significance":"If the central claim were sound, the paper would be a useful test of whether a fractional-calculus holographic dark energy density can be embedded in an anisotropic Bianchi spacetime while remaining consistent with late-time cosmological data. The use of multiple public data sets and a standard MCMC pipeline is appropriate, and the reported H0 and Ω_m0 values are plausible. However, the load-bearing derivation is not a solution of the Einstein equations stated in Section 3, and the fitted H(z) is exactly the isotropic flat wCDM expression with no dependence on the Bianchi index n or on the fractional density evolution. The paper therefore does not establish the observational viability of the proposed model; it establishes only that a wCDM-like expansion history fits the data.","major_comments":[{"comment":"The full LRS Bianchi field equations are not satisfied by the fitted expansion history. Eq. (24) follows from the 00 equation (13) together with the continuity equations, but the pressure equations (14) and (15) impose an additional constraint. Substituting Eq. (24) into either (14) or (15) gives ω_d Ω_d = 1 for n > 2, which forces Ȟ/H² = −3 and hence q = 2. The fitted values (ω_d ≈ −0.78, Ω_d ≈ 0.73) give Ȟ/H² ≈ −0.54, so the H(z) used in the likelihood is not a solution of the Einstein equations presented in Section 3. Equivalently, subtracting (14) from (15) yields Ȟ/H² = −3n/(n+2), which for n > 2 gives q = −1 + 3n/(n+2) > 0.8, showing that an accelerating solution of the full system does not exist under the ansatz A1 ∝ A2^n with n > 2.","section":"§3, Eqs. (13)–(15) and (24)"},{"comment":"The fitted Hubble law is a wCDM proxy, not the fractional-HDE model. Eq. (26) gives ω_d as a function of Ω_d (and therefore of z), and it is equivalent to the conservation equation (16) combined with ρ_d = 3c² H^{(3α−2)/α}. Eq. (27), by contrast, integrates (16) with ω_d held constant, which is an additional assumption not implied by the model. The text then takes a single value of ω_d from Eq. (26) and inserts it into Eq. (29). Consequently α enters the likelihood only through a constant equation-of-state parameter, the fractional density (20) is never evolved, and the Bianchi parameter n cancels from Eq. (29). The MCMC constraints in Table 1 therefore characterize an isotropic flat wCDM model, not the LRS Bianchi fractional-HDE model claimed in the abstract.","section":"§3, Eqs. (26)–(29)"},{"comment":"All physical predictions are made with α = 1.01, whereas the MCMC constraints in Table 1 give α ≈ 0.86–0.885 for every data combination. No derivation of the value 1.01 is provided, and it is not the best-fit value. The claimed transition redshift z_t = 0.55, present values q ≈ −0.35, ω_d ≈ −0.77, and statefinder (r,s) = (0.74,0.07) are therefore not the predictions of the fitted model. This disconnect undermines the central viability claim in the abstract, which is based on the α = 1.01 curves rather than on the observationally constrained parameters.","section":"§5, Figs. 7–11; Table 1"}],"minor_comments":[{"comment":"Equation (29) is printed without the exponent 1/2, making it dimensionally inconsistent as a formula for H(z); it should read H(z) = H0 [Ω_m0(1+z)^3 + (1−Ω_m0)(1+z)^{3(1+ω_d)}]^{1/2}.","section":"§3, Eq. (29)"},{"comment":"The text reports α = 0.885^{+0.121}_{-0.149} for the H(z) dataset, but Table 1 lists 0.885^{+0.141}_{-0.149}; these should be reconciled.","section":"§4.6, Table 1"},{"comment":"The text states that α = 1.01 and Ω_d0 = 0.73 are 'based on observational analysis', but neither value appears in Table 1, and the uncertainty on Ω_d0 is not propagated into the plotted curves for q, ω_d, and the statefinder.","section":"§5.1 and §5.2"},{"comment":"'No data are associated with the manuscript' is inaccurate, since the analysis uses the 46-point Hubble sample, 15 BAO points, and the Pantheon+SH0ES catalog; the specific data versions or repository links should be cited.","section":"Data Availability"}],"recommendation":"reject","confidential_remarks":"The central calculation error is load-bearing: the full Einstein equations force ω_d Ω_d = 1 under the stated ansatz, which is incompatible with the fitted parameters, and the fitted H(z) is simply the flat wCDM expression. These are not presentation issues; they invalidate the main claims. The manuscript would need to be substantially reworked, either by solving the full Bianchi system consistently or by reframing the contribution as a wCDM parametrization without the anisotropic and fractional-model claims, before it could be considered for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the fractional-HDE density in Eq. (5) is a genuinely new construct, but the paper's observational analysis never tests it. Eq. (29) is exactly flat wCDM: the Bianchi exponent n cancels and the fractional density is replaced by a constant-EoS power law. So the MCMC constraints on H0, alpha, Omega_m0 are constraints on wCDM, not on the proposed anisotropic model.\n\nWhat the paper does well: deriving an HDE density from fractional black-hole entropy is a novel idea worth a footnote, and the MCMC pipeline (emcee, 46 H(z), 15 BAO, Pantheon+Shoes) is competently executed. The reported H0 around 67.7 matches Planck, but that consistency is a property of the wCDM proxy, not of the Bianchi geometry.\n\nThe soft spots are load-bearing. The assumed A1 ∝ A2^n relation with a single isotropic pressure cannot solve the field equations. Equating the pressure equations (14) and (15) forces omega_d * Omega_d = 1 for n > 2; the fitted values (omega_d ≈ -0.8, Omega_d ≈ 0.73) give about -0.58, so the expansion history used in the fits is not a solution of the model's own Einstein equations. The paper derives Eq. (24) from the 00 equation and conservation, but never enforces (14) and (15). Separately, the fit treats omega_d as constant while Eq. (26) gives a redshift-dependent omega_d(z); these are not the same object. And the figures use alpha = 1.01, which is higher than the MCMC best fits (0.86–0.89) for most datasets; describing it as 'based on observational analysis' is misleading. The statefinder pair (0.74, 0.07) is computed from the same fitted parameters, so it is a consistency check of the fitting formula, not an independent test.\n\nVerdict: the paper reduces to a standard wCDM fit with a new label. It does not establish the observational viability of anisotropic fractional HDE. The internal inconsistency with the pressure equations is a fatal flaw that a revision cannot fix short of abandoning the Bianchi setup. I would not send it to peer review. The audience is limited to people cataloguing HDE constructions, who might cite Eq. (5) as an alternative entropy-based density.","headline":"A competent wCDM fit wearing an anisotropic fractional-HDE costume: the Bianchi geometry drops out of the observable and the model contradicts its own pressure equations.","tokens_in":15978,"tokens_out":2555,"would_cite":false,"duration_ms":26532,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k","95.36.+x"],"model":"deepseek-v4-flash","headline":"The paper claims a fractional holographic dark-energy model with one free parameter reproduces ΛCDM's late-time expansion when fitted to Hubble, BAO, and supernova data.","keywords":["fractional holographic dark energy","LRS Bianchi type-I model","dark energy","observational constraints","Markov Chain Monte Carlo","statefinder diagnostics","deceleration parameter","ΛCDM comparison"],"falsifier":"Evaluate the two pressure equations of the LRS Bianchi field equations, Eqs. (14) and (15), at the paper's fitted parameters $\\alpha \\approx 0.87$, $\\Omega_{d0} \\approx 0.73$, $\\omega_d \\approx -0.77$: consistency of the metric with a single isotropic pressure forces $\\omega_d \\Omega_d = 1$, whereas the fitted values give $\\omega_d \\Omega_d \\approx -0.56$. That direct substitution settles whether the claimed Bianchi solution exists at all, independently of the goodness of the data fits.","tokens_in":14851,"feed_emoji":"🌌","tokens_out":21813,"duration_ms":205375,"temperature":0.7,"pith_summary":"This paper argues that a new 'fractional holographic dark energy' — an energy density obtained from the holographic principle together with fractional-calculus corrections to black-hole entropy — can drive the universe's late-time acceleration as effectively as the cosmological constant. The authors place this density in a locally rotationally symmetric Bianchi type-I universe, derive a parameterized Hubble law, and fit its three parameters ($H_0$, the fractional index $\\alpha$, and the matter density $\\Omega_{m0}$) to Hubble, BAO, and Pantheon+SH0ES data using Markov Chain Monte Carlo. Across four dataset combinations the fits return $H_0 \\approx 67.7$–$67.8$ km/s/Mpc, $\\Omega_{m0} \\approx 0.26$–$0.27$, and $\\alpha \\approx 0.86$–$0.89$. Taking $\\alpha = 1.01$, the model produces a deceleration-parameter transition at $z_t = 0.55$, an equation-of-state parameter that approaches $-1$, and statefinder values $(r, s) = (0.74, 0.07)$ heading toward the $\\Lambda$CDM fixed point, which the authors read as demonstrating the observational viability of fractional holographic dark energy.","feed_headline":"Fractional dark-energy model fits cosmic data like ΛCDM","feed_subtitle":"MCMC fits to Hubble, BAO and supernova data put the transition to acceleration at z=0.55, matching the standard model.","key_machinery":"The load-bearing object is the fractional holographic dark-energy density $\\rho_d = 3c^2 L^{-(3\\alpha-2)/\\alpha}$, formed by inserting the fractional entropy law $S_h \\propto A^{(\\alpha+2)/2\\alpha}$ into the holographic inequality and taking the Hubble horizon $L = H^{-1}$ as the infrared cutoff; the fractional order $\\alpha \\in (1, 2]$ is the single new parameter, with $\\alpha = 2$ recovering ordinary holographic dark energy. Combined with the LRS Bianchi type-I metric $ds^2 = -dt^2 + A_1^2 dx^2 + A_2^2 (dy^2 + dz^2)$ and the assumption $A_1 \\propto A_2^n$ with $n > 2$, this density yields the field equations, the dark-energy equation of state $\\omega_d = -1 + (3\\alpha - 2)(1 - \\Omega_d)/(2\\alpha - (3\\alpha - 2)\\Omega_d)$, and, treating dark matter as pressureless and $\\omega_d$ as constant, the fitted Hubble law $H(z) = H_0[\\Omega_{m0}(1+z)^3 + (1-\\Omega_{m0})(1+z)^{3(1+\\omega_d)}]$. The deceleration parameter $q$ and the statefinder pair $(r, s)$ derived from the same equations supply the diagnostics used to compare the model with $\\Lambda$CDM.","core_discovery":"The central claim is that fractional holographic dark energy, with density $\\rho_d = 3c^2 H^{(3\\alpha-2)/\\alpha}$ under the Hubble-horizon cutoff, is an observationally viable alternative to the cosmological constant. The fractional index $\\alpha$ enters through the entropy law $S_h \\propto A^{(\\alpha+2)/2\\alpha}$ obtained from the fractional Wheeler-DeWitt equation, and $\\rho_d$ reduces to the standard holographic dark energy density when $\\alpha = 2$. Fitted to the combined data, the model gives a present dark-energy density parameter $\\Omega_{d0} \\approx 0.73$, an equation-of-state parameter $\\omega_d \\approx -0.77$ that drifts toward $-1$ at late times, and a deceleration parameter that crosses from deceleration to acceleration at $z_t \\approx 0.55$, with the statefinder pair $(r, s) = (0.74, 0.07)$ approaching the $\\Lambda$CDM point $(1, 0)$. The authors conclude that the LRS Bianchi type-I fractional holographic dark energy model is observationally viable, that its fitted $H_0$ agrees with early-universe CMB measurements, and that it closely reproduces the $\\Lambda$CDM expansion history.","pith_inferences":["Because the fitted Hubble law is algebraically a flat wCDM model with a constant dark-energy equation of state, current distance data cannot distinguish fractional holographic dark energy from a simpler two-parameter dark-energy model; matching $\\Lambda$CDM is a consistency check rather than a detection of fractional physics.","The paper's dynamical plots use $\\alpha = 1.01$, yet the fitted values cluster near $\\alpha \\approx 0.87$, outside the model's stated range $1 < \\alpha \\le 2$; recomputing the deceleration and statefinder trajectories at the central fitted value is a direct check of whether the $\\Lambda$CDM-like conclusions survive.","The paper treats $\\omega_d$ as constant in deriving $H(z)$ even though its own Eq. (26) defines a redshift-dependent $\\omega_d$; redoing the fit with an evolving $\\omega_d(z)$ is a natural extension that could shift the best-fit parameters and the transition redshift.","The Bianchi geometry carries a shear that the fitting procedure ignores, so the same model could be tested against cosmic-microwave-background isotropy bounds on anisotropy, an independent check the paper does not pursue."],"forward_implications":["If the model is correct, the fitted $H_0 \\approx 67.7$–$67.8$ km/s/Mpc matches early-universe CMB estimates, meaning this dark-energy model does not by itself resolve the Hubble tension.","The deceleration-to-acceleration transition at $z_t \\approx 0.55$ falls inside the observationally favored range $0.5 \\le z_t \\le 0.8$, consistent with supernova and BAO evidence.","The dark-energy equation of state that starts near zero at high redshift, reaches about $-0.77$ today, and approaches $-1$ at late times reproduces the quintessence-to-$\\Lambda$CDM trajectory without introducing a cosmological constant.","The consistent parameter estimates across all four dataset combinations indicate that the $\\Lambda$CDM-like behavior is not an artifact of any single probe."],"supporting_citations":[{"why":"Formulates the holographic principle that lets the dark-energy density be written with a single length scale.","marker":"[17]"},{"why":"Supplies the holographic inequality that the fractional dark-energy density is built from.","marker":"[19]"},{"why":"Derives the fractional black-hole entropy whose power-law form generates the fractional exponent in the new dark-energy density.","marker":"[30]"},{"why":"Defines the standard holographic dark energy framework that the fractional model extends and recovers at $\\alpha = 2$.","marker":"[18]"},{"why":"Introduces the statefinder pair $(r, s)$ used to compare the model's trajectory with the $\\Lambda$CDM fixed point.","marker":"[42]"},{"why":"Provides the 46 cosmic-chronometer $H(z)$ measurements that anchor the Hubble likelihood.","marker":"[33]"},{"why":"Defines the comoving distance and dilation-scale quantities used in the BAO likelihood.","marker":"[34]"},{"why":"Supplies the Pantheon+SH0ES supernova sample of 1701 distance moduli used in the fits.","marker":"[36]"},{"why":"Gives the early-universe CMB values of $H_0$ and $\\Omega_{m0}$ against which the fitted parameters are judged consistent.","marker":"[5]"}],"fun_headline_variants":["Fractional dark energy passes cosmic tests, mimics ΛCDM","Fractional entropy law leads to viable dark energy model","Cosmic data favors fractional dark energy, close to ΛCDM","Fractional dark energy model mirrors ΛCDM expansion","New dark energy model fits cosmic data, resembles ΛCDM"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumption that the LRS Bianchi type-I field equations are a consistent basis for the fitted expansion history, but equating the two pressure equations forces $\\omega_d \\Omega_d = 1$ while the fitted values ($\\omega_d \\approx -0.77$, $\\Omega_d \\approx 0.73$) give $\\omega_d \\Omega_d \\approx -0.56$, so the Bianchi structure drops out of the dynamics.","fun_headline_variants_meta":{"raw":{"variants":["Fractional dark energy passes cosmic tests, mimics ΛCDM","Fractional entropy law leads to viable dark energy model","Cosmic data favors fractional dark energy, close to ΛCDM","Fractional dark energy model mirrors ΛCDM expansion","New dark energy model fits cosmic data, resembles ΛCDM"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001085,"raw_usage":{"total_tokens":4690,"prompt_tokens":1254,"completion_tokens":3436,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":870,"completion_tokens_details":{"reasoning_tokens":3353}},"tokens_in":870,"tokens_out":3436,"duration_ms":30088,"temperature":1.0,"reasoning_tokens":3353,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:46:31.442769+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the two pressure equations of the LRS Bianchi field equations, Eqs. (14) and (15), at the paper's fitted parameters $\\alpha \\approx 0.87$, $\\Omega_{d0} \\approx 0.73$, $\\omega_d \\approx -0.77$: consistency of the metric with a single isotropic pressure forces $\\omega_d \\Omega_d = 1$, whereas the fitted values give $\\omega_d \\Omega_d \\approx -0.56$. That direct substitution settles whether the claimed Bianchi solution exists at all, independently of the goodness of the data fits.","supporting_citations":[{"cited_title":"Effective field t heory, black holes, and the cosmological constant","cited_arxiv_id":null,"evidence_quote":"Supplies the holographic inequality that the fractional dark-energy density is built from."},{"cited_title":"Prospecting black hole ther- modynamics with fractional quantum mechanics","cited_arxiv_id":null,"evidence_quote":"Derives the fractional black-hole entropy whose power-law form generates the fractional exponent in the new dark-energy density."},{"cited_title":"Wang, Yi Wang and Miao Li, Holographic Dark Energy, Phys","cited_arxiv_id":null,"evidence_quote":"Defines the standard holographic dark energy framework that the fractional model extends and recovers at $\\alpha = 2$."},{"cited_title":"Statefinder—A new geometrical diagnost ic of dark energy","cited_arxiv_id":null,"evidence_quote":"Introduces the statefinder pair $(r, s)$ used to compare the model's trajectory with the $\\Lambda$CDM fixed point."},{"cited_title":"Cosmic chronometers: constraining the equat ion of state of dark energy. I: H(z) measurements","cited_arxiv_id":null,"evidence_quote":"Provides the 46 cosmic-chronometer $H(z)$ measurements that anchor the Hubble likelihood."},{"cited_title":"Baryon acous tic oscillations in the Sloan Digital Sky Survey Data Release 7 galaxy sample","cited_arxiv_id":null,"evidence_quote":"Defines the comoving distance and dilation-scale quantities used in the BAO likelihood."},{"cited_title":"Scolnic et al., Astrophys","cited_arxiv_id":null,"evidence_quote":"Supplies the Pantheon+SH0ES supernova sample of 1701 distance moduli used in the fits."},{"cited_title":"Aghanim, Y","cited_arxiv_id":null,"evidence_quote":"Gives the early-universe CMB values of $H_0$ and $\\Omega_{m0}$ against which the fitted parameters are judged consistent."}],"review_version":1}