{"id":"65027555-32bc-4490-b956-8702a307066e","arxiv_id":"2508.06586","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":11,"one_line_summary":"Bouncing universe solutions in f(Q,C) modified gravity on a Bianchi I background are claimed, but the central field equations fail the GR limit and the stability analysis is contradicted by the paper's own figures.","lead":"A cosmology model that has the universe contract and then expand again, built by changing Einstein's equations and picking a special geometry. The math looks internally inconsistent: the main equations give the wrong answer in the simplest limit, and the claimed stability check shows huge perturbations, not small ones.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Field equations are derived from isotropic Q,C,G for an anisotropic Bianchi I metric; they reduce to GR only at Ω=1, so the bounce solutions do not solve the stated f(Q,C) theory.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing defect: the paper substitutes isotropic forms for Q, C, and the Einstein tensor while analyzing an anisotropic Bianchi I spacetime. My independent check of the Ricci scalar confirms the isotropic expressions only hold for Ω=1. This is not a matter of external consensus or model choice; it is an internal inconsistency: the equations being solved do not follow from the stated action and metric. Since every result—bounce existence, NEC violation, EoS, stability—flows from Eqs. (17)–(18), the central claim is unsupported. The stability analysis additionally shows perturbation amplitudes of order 10^5–10^10, contradicting the assumed smallness, but the field-equation inconsistency is more fundamental. I find no reason to alter the reader's REJECT verdict; no independent support (e.g., reproducible code or numerical simulation) is provided to override this. Thus the verdict remains unchanged, with high confidence in rejection.","tokens_in":22746,"tokens_out":11264,"duration_ms":111019,"concrete_test":"Independently derive Q, C, and G_ab for the metric ds^2=-dt^2+a^2dx^2+b^2(dy^2+dz^2) with a=b^Ω, using either the coincident-gauge connection (all Γ=0) or the paper's stated 'vanishing affine connection' convention. Compare with Eqs. (13)–(16). Then insert f(Q,C)=Q into the full Bianchi I field equations and check whether Eqs. (17)–(18) are recovered. The test fails if any difference appears for Ω≠1; if the comparison is done numerically, use e.g. Ω=4.5 and a generic b(t)=e^{t^2/2} to see nonzero residuals.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that f(Q,C) gravity admits non-singular bouncing solutions in Bianchi I. The load-bearing step is Eqs. (13)–(16), where the authors replace the anisotropic Bianchi I geometry with isotropic expressions: G_ab = -h_ab(3H^2+2Hdot)+3H^2 u_a u_b, Q=6H^2, C=6(H^2+Hdot), with H=(Ω+2)bdot/(3b). For the metric (12), a=b^Ω, the actual Levi-Civita Ricci scalar is R=2(Ω+2)Bdot + (2Ω^2+4Ω+6)B^2 with B=bdot/b, whereas the isotropic expression 12H^2+6Hdot = 2(Ω+2)Bdot + (4/3)(Ω+2)^2 B^2. These match only at Ω=1. Likewise, the coincident-gauge non-metricity scalar Q for this diagonal anisotropic metric is not 6H^2 unless Ω=1. The Einstein tensor is anisotropic (R_xx/a^2 ≠ R_yy/b^2 generally), so Eq. (13) cannot hold. Consequently Eqs. (17)–(18) are not the field equations of f(Q,C) gravity for this spacetime. The inconsistency is fatal to the central claim: in the GR limit f(Q,C)=Q, Eqs. (17)–(18) give ρ=-(5+4Ω)H^2 and P=(7-2Ω)H^2-2Hdot, which do not reduce to the Bianchi I Friedmann equations and even produce negative density for the quoted Ω values. The bounce solutions, energy condition violations, and stability analysis are all built on these incorrect equations, so they do not support the claimed feasibility of bouncing solutions in this theory.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims to construct feasible non-singular bouncing cosmological solutions in f(Q,C) gravity for a Bianchi type-I spacetime with a perfect fluid. It imposes the constraint a(t)=b(t)^Ω, adopts an isotropic-like set of geometric quantities (G_ab, Q, C) in Eqs. (13)-(16), and derives field equations (17)-(18). Three forms of f(Q,C) are then studied: μ1Q+μ2C^2, μ3Q^{λ+1}+μ2C^2, and Q+μ4/Q+μ2C^2. The paper presents scale-factor, Hubble, deceleration, density, pressure, equation-of-state, energy-condition, Hubble-radius, and redshift plots, and it argues from a stability analysis that the bounce is stable. The central claims are that the null energy condition is violated, that ω≈-1.03, and that the models are stable under linear perturbations.","tokens_in":23369,"tokens_out":6009,"duration_ms":70195,"significance":"If correct, the paper would provide an interesting example of bouncing cosmology in a modified symmetric-teleparallel gravity and would extend the literature on f(Q,C) models to anisotropic geometries. The topic is relevant, and the authors correctly emphasize that a viable bounce requires NEC violation. However, the central derivation is inconsistent: the field equations used do not reduce to general relativity in the limit f(Q,C)=Q, and the stability analysis does not actually analyze perturbations. The reported physical conclusions are therefore not established.","major_comments":[{"comment":"The field equations used in the paper are not the field equations of f(Q,C) gravity for the stated Bianchi I geometry. Equation (13) replaces the Einstein tensor of the anisotropic metric (12) by the isotropic expression -h_ab(3H^2+2Hdot)+3H^2 u_a u_b. For a=b^Ω with Ω≠1, the Einstein tensor is anisotropic, and the non-metricity scalar Q in the coincident gauge is not generally 6H^2. More decisively, the GR limit f(Q,C)=Q in Eqs. (17)-(18) yields ρ=-(5+4Ω)H^2 and P=(7-2Ω)H^2-2Hdot. Even at Ω=1 this gives ρ=-9H^2 and P=5H^2-2Hdot instead of the standard Friedmann results ρ=3H^2 and P=-(3H^2+2Hdot). Thus the equations being solved are not those of the stated theory, and all subsequent model results, energy-condition violations, and stability statements rest on this incorrect starting point.","section":"§2, Eqs. (13)-(18)"},{"comment":"The stability analysis is not a valid linear perturbation analysis. Equation (47) assumes δΩ is a small perturbation, but the subsequent calculation substitutes the background solution into the background conservation equation (48) and solves for δΩ(t). This is an algebraic rearrangement of the background equations, not an evolution equation for perturbations. Moreover, Fig. 14 shows δΩ values of order 10^5-10^10, which are not small and violate the linearization assumption. The plotted curves increase with time, contradicting the text's statement that the perturbations 'slowly decrease and approach to zero.' The stability conclusion is therefore unsupported.","section":"§5, Eqs. (47)-(51) and Fig. 14"},{"comment":"The bounce is imposed rather than derived from the theory. The parametric scale factor (28) is chosen with specific values of θ1, θ2, θ3, n so that H crosses zero, and the subsequent NEC violation and ω≈-1.03 are consequences of that choice. Similarly, Sec. 4.8 fixes q=-0.831 from observation to set b(t)=t^{1/(1+q)} and then computes ρ and P as functions of z. These procedures are reconstructions from imposed kinematics or observational input, not predictions of the f(Q,C) dynamics. The paper would need to show that the field equations select such a bounce, or that the results are robust across a genuinely wide parameter range, before claiming that the framework 'admits' these solutions.","section":"§4.1 and §4.8, Eqs. (28) and (37)"}],"minor_comments":[{"comment":"The text states that ω=-1.03±0.03 'aligns with the results reported by the Planck collaboration' and then says this 'meets the conditions for the quintessence regime (−1<ω<−1/3).' Since −1.03 is in the phantom regime (ω<−1), this classification is internally inconsistent. The abstract and final remarks also refer to both quintessence and phantom, which should be reconciled.","section":"§4.5, Eqs. (31)-(33)"},{"comment":"The claim that increasing Ω leads to a later bounce time is confounded: in Tables 1 and 2 both Ω and θ2 (or θ3) change simultaneously. A controlled comparison varying only Ω would be needed to support that statement.","section":"§4.2, Tables 1-2"},{"comment":"The manuscript contains passages of the form 'In the revised manuscript, we have included ...', which appear to be remnants of a response to a previous referee report. These should be removed or rewritten for a self-contained submission.","section":"§2 and §6"},{"comment":"Many plotted quantities (e.g., ρ and P in Figs. 4-6 and 12-13) are displayed as large numbers without units or with unspecified units. Dimensional analysis and a statement of the chosen parameter values and their units would greatly improve reproducibility.","section":"General"}],"recommendation":"reject","confidential_remarks":"The GR-limit check of Eqs. (17)-(18) is unambiguous and sufficient for rejection: the central equations are not the field equations of the stated theory. Even if that were repaired, the stability analysis in Sec. 5 would need to be replaced by a genuine perturbation study. The manuscript also contains self-referential 'in the revised manuscript' passages that suggest it may be an incomplete resubmission; the editor may wish to verify the provenance of the text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the central claim doesn't hold. The authors set up a Bianchi I metric with a=b^Ω, then use the isotropic expressions Q=6H^2, C=6(H^2+Ḣ), and G_ab=-h_ab(3H^2+2Ḣ)+3H^2 u_a u_b, with H=(Ω+2)ḃ/(3b). Those expressions are correct for FLRW, not for this metric unless Ω=1. The actual Einstein tensor for a=b^Ω is anisotropic, and the curvature scalars differ. So Eqs. (17)-(18) are not the field equations of f(Q,C) gravity for the stated spacetime. The simplest check: set f(Q,C)=Q and the equations give ρ=-(4Ω+5)H^2, P=(7-2Ω)H^2-2Ḣ, which don't reduce to GR even in the isotropic limit—and for the quoted Ω values ρ is negative. That's fatal to the paper's stated goal.\n\nThere is some value here. The paper surveys three f(Q,C) models, uses a parametric scale factor from Singh et al., and is thorough about plotting energy conditions, EoS, Hubble radius, and redshift behavior. The algebra is extensive and internally consistent with the (incorrect) equations. The literature review is broad, though heavily self-referential in places.\n\nOther soft spots, in order of severity: the stability analysis is not a stability analysis in any meaningful sense. The perturbation δΩ is defined as small, but the plotted values are 10^5 to 10^10, and Figure 14 shows them growing or staying large while the text claims they \"slowly decrease and approach to zero.\" That is an internal contradiction. The redshift section uses the observed q=-0.831 to fix the power-law scale factor and then \"predicts\" ρ and P from it—circular. There are also leftover revision notes (\"In the revised manuscript, we have included...\") and garbled references, e.g. \"Comgamman. Math. Phys.\"\n\nI don't think this deserves referee time as it stands. The error in the field equations is load-bearing and would require redoing the whole analysis. A desk reject with a clear explanation is appropriate. That said, the paper could serve a useful pedagogical purpose: it is a textbook example of what goes wrong when you apply isotropic formulas to an anisotropic background in modified gravity.\n\nBottom line: don't cite it, don't build on it, and don't send it to referees unless the authors first fix the field equations and redo the analysis. If it comes back with corrected equations, I'd be willing to look again.","headline":"The paper's central bounce results are not supported: the field equations are the isotropic FLRW ones applied to an anisotropic Bianchi I metric, and they fail the GR limit for Ω≠1.","tokens_in":23797,"tokens_out":4250,"would_cite":false,"duration_ms":43678,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.Kd","04.20.Dw","04.40.Dg"],"model":"deepseek-v4-flash","headline":"An f(Q,C) modified-gravity theory is shown to admit non-singular bouncing solutions in a Bianchi type-I universe with a perfect fluid, with null energy condition violation and stable linear perturbations.","keywords":["f(Q,C) gravity","non-metricity","bouncing cosmology","Bianchi type-I","null energy condition","stability analysis","dark energy","cosmological singularity"],"falsifier":"Compute the exact Einstein tensor components of metric (12) with $a=b^\\Omega$ and $\\Omega\\neq1$, then insert the three $f(\\mathcal{Q},\\mathcal{C})$ models into the full field equations (11). If the resulting dynamics no longer give an $H=0$ crossing with positive density and negative pressure, the claimed bounce is an artifact of the isotropic reduction used in (17)-(18).","tokens_in":22621,"feed_emoji":"🌌","tokens_out":10339,"duration_ms":100825,"temperature":0.7,"pith_summary":"The paper tries to establish that $f(\\mathcal{Q},\\mathcal{C})$ gravity—a modified theory built from the non-metricity scalar and its boundary term—can drive a non-singular bouncing cosmology in a Bianchi type-I spacetime filled with a perfect fluid. Using three functional forms of $f$, it shows the universe can pass from contraction through $H=0$ into expansion, with positive energy density, negative pressure, and violation of the null energy condition at the bounce. It reports an equation-of-state parameter near $-1.03$, consistent with an accelerating phase, and claims that linear perturbations of the Hubble parameter decay, so the bounce is stable. If correct, this gives a singularity-free alternative to the big-bang model and a geometric source for dark-energy-like acceleration without a cosmological constant.","feed_headline":"Bounce replaces big bang in non-metric gravity","feed_subtitle":"Three f(Q,C) models cross from contraction to expansion, violate the null energy condition, and stay stable.","key_machinery":"The load-bearing machinery is the reduction of the $f(\\mathcal{Q},\\mathcal{C})$ field equations on a Bianchi type-I metric with the constraint $a=b^\\Omega$ to the two ordinary differential equations (17)-(18), using the identifications $\\mathcal{Q}=6H^2$, $\\mathcal{C}=6(H^2+\\dot H)$, and the average Hubble parameter $H=(\\Omega+2)\\dot b/(3b)$. Substituting the parametric scale factor $b=\\beta e^{\\vartheta_3 t^{n+1}/(n+1)+\\vartheta_2 t^n/n+\\vartheta_1 t}$ turns those equations into explicit density and pressure expressions for each model; the bounce is the $H=0$ crossing, and the perturbation ansatz $H_{\\rm pert}(t)=H(t)(\\delta_\\Omega(t)+1)$ together with the conservation equation produces the","core_discovery":"The central claim is that $f(\\mathcal{Q},\\mathcal{C})$ gravity—where $\\mathcal{Q}$ is the non-metricity scalar and $\\mathcal{C}=\\bar{R}-\\mathcal{Q}$ is the boundary term—admits viable non-singular bounces. For three model choices, the scale factor reaches a minimum, the Hubble parameter runs from negative through zero to positive, the energy density stays positive while pressure is negative, and the null energy condition is violated near the bounce. The equation-of-state parameter is about $-1.03$, and a linear perturbation of the Hubble parameter decays with time. The paper interprets these features as evidence that this modified-gravity framework can resolve the initial-singularity problem","pith_inferences":["Editorial extension: because $\\mathcal{Q}=6H^2$ and $\\mathcal{C}=6(H^2+\\dot H)$ are isotropic reductions, the same analysis should be repeated with the exact anisotropic Einstein tensor of metric (12); for $\\Omega\\neq1$ the two sets of equations are not guaranteed to agree, so the quantitative profiles may change.","Editorial extension: the parametric scale factor can be used in reverse—prescribing any desired bounce profile $b(t)$ determines an $f(\\mathcal{Q},\\mathcal{C})$ that realizes it—making the framework a reconstruction tool for bounce models.","Editorial extension: the stability test perturbs only the Hubble parameter; a fuller test would perturb the metric components and matter fields in the full Bianchi type-I equations, which could expose additional modes near the bounce."],"forward_implications":["If the results are correct, the big-bang singularity can be replaced by a smooth bounce in $f(\\mathcal{Q},\\mathcal{C})$ gravity, giving a non-singular early universe.","All three models violate the null energy condition near the bounce, so the breach appears to be a feature of the framework rather than of one special choice of parameters.","The equation-of-state parameter near $-1.03$ can reproduce the observed late-time acceleration, so the theory can mimic dark energy without a cosmological constant.","The decaying perturbation $\\delta_\\Omega$ indicates the bounce is robust to small fluctuations of the Hubble rate, allowing the model to survive into later cosmic epochs.","The bounce occurs in a Bianchi type-I background with small anisotropy, connecting early anisotropic conditions to the nearly isotropic late universe."],"supporting_citations":[{"why":"Introduces $f(\\mathcal{Q},\\mathcal{C})$ gravity and supplies the field equations that the paper solves.","marker":"[39]"},{"why":"Supplies the parametric scale factor used to construct the bouncing profiles.","marker":"[80]"},{"why":"Provides the two model forms $f=\\mu_1\\mathcal{Q}+\\mu_2\\mathcal{C}^2$ and $f=\\mu_3\\mathcal{Q}^{\\lambda+1}+\\mu_2\\mathcal{C}^2$.","marker":"[92]"},{"why":"Supplies the inverse-$\\mathcal{Q}$ model $f=\\mathcal{Q}+\\mu_4/\\mathcal{Q}+\\mu_2\\mathcal{C}^2$.","marker":"[93]"},{"why":"Justifies the constant shear-to-expansion ratio that underlies the constraint $a=b^\\Omega$.","marker":"[77]"},{"why":"Gives the reference equation-of-state values the paper compares its result near $-1.03$ against.","marker":"[86]"}],"fun_headline_variants":["Bouncing universe in f(Q,C) gravity skips singularity","Non-metric gravity bounces: no big bang, stable perturbation","Bianchi-I bounce: null energy violation, stable evolution","Scale factor minimum confirms bounce in f(Q,C) gravity","Modified gravity offers bounce without initial singularity"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The results rest on treating the anisotropic Bianchi type-I metric (12) with $a=b^\\Omega$ as though it had the isotropic forms $\\mathcal{Q}=6H^2$, $\\mathcal{C}=6(H^2+\\dot H)$, and $G_{ab}=-h_{ab}(3H^2+2\\dot H)+3H^2u_au_b$ with $H=(\\Omega+2)\\dot b/(3b)$; for $\\Omega\\neq1$ those reductions are not automatically the Einstein tensor of metric (12).","fun_headline_variants_meta":{"raw":{"variants":["Bouncing universe in f(Q,C) gravity skips singularity","Non-metric gravity bounces: no big bang, stable perturbation","Bianchi-I bounce: null energy violation, stable evolution","Scale factor minimum confirms bounce in f(Q,C) gravity","Modified gravity offers bounce without initial singularity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000886,"raw_usage":{"total_tokens":3641,"prompt_tokens":702,"completion_tokens":2939,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":446,"completion_tokens_details":{"reasoning_tokens":2868}},"tokens_in":446,"tokens_out":2939,"duration_ms":20135,"temperature":1.0,"reasoning_tokens":2868,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T23:00:34.827273+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact Einstein tensor components of metric (12) with $a=b^\\Omega$ and $\\Omega\\neq1$, then insert the three $f(\\mathcal{Q},\\mathcal{C})$ models into the full field equations (11). If the resulting dynamics no longer give an $H=0$ crossing with positive density and negative pressure, the claimed bounce is an artifact of the isotropic reduction used in (17)-(18).","supporting_citations":[{"cited_title":"and Saridakis, E.N.: J","cited_arxiv_id":null,"evidence_quote":"Introduces $f(\\mathcal{Q},\\mathcal{C})$ gravity and supplies the field equations that the paper solves."},{"cited_title":"et al.: J","cited_arxiv_id":null,"evidence_quote":"Supplies the parametric scale factor used to construct the bouncing profiles."},{"cited_title":"et al.: Eur","cited_arxiv_id":null,"evidence_quote":"Provides the two model forms $f=\\mu_1\\mathcal{Q}+\\mu_2\\mathcal{C}^2$ and $f=\\mu_3\\mathcal{Q}^{\\lambda+1}+\\mu_2\\mathcal{C}^2$."},{"cited_title":"et al.: Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the inverse-$\\mathcal{Q}$ model $f=\\mathcal{Q}+\\mu_4/\\mathcal{Q}+\\mu_2\\mathcal{C}^2$."},{"cited_title":"and Sachs, R.K.: J","cited_arxiv_id":null,"evidence_quote":"Justifies the constant shear-to-expansion ratio that underlies the constraint $a=b^\\Omega$."},{"cited_title":"and Vagenas, E.C.: Phys","cited_arxiv_id":null,"evidence_quote":"Gives the reference equation-of-state values the paper compares its result near $-1.03$ against."}],"review_version":1}