{"id":"3e48d626-eb87-471d-bcb8-9b023b33bcf0","arxiv_id":"2504.21504","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A T^2-dependent coupling between the derivative of T_mu nu T^mu nu and the baryon current can generate the observed baryon asymmetry in GR and f(R,T^2) gravity, at the cost of fitted parameters.","lead":"This paper adds a new interaction between matter and gravity, built from the square of the energy-momentum tensor, to the standard gravitational baryogenesis mechanism. The authors show that with suitably chosen cutoff and decoupling parameters, this term can reproduce the observed matter-antimatter asymmetry while respecting Big Bang Nucleosynthesis bounds.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The proposed dimension-12 operator is introduced without a UV derivation and without a specified B- or B-L-violating sector; the quoted successes amount to fitting M_* and T_D to the observed asymmetry, so the central mechanism is conditional rather than explanatory.","rationale":"The paper's central calculation is a plausible phenomenological exercise, and I verified that the GR branch (Eq. 70) reproduces the quoted asymmetry with the stated parameters. The most load-bearing weakness is the one the reader identified: the new interaction term is an arbitrary EFT operator with a free cutoff M_*, and the required B or B-L-violating process is explicitly deferred. This is not merely a lack of UV completion; it means the mechanism cannot yet produce a net baryon number, because the partial(T^2)J coupling only biases an unspecified B-violating sector. The fact that n_b/s scales as T_D^9/M_*^8 in the GR case confirms that two free parameters are being tuned to one observed number, so the agreement is not a predictive success. Two further internal problems reinforce the conditional verdict: the n=1/2 successful benchmark uses T_D=2e16 GeV, which violates the paper's own inflationary-scale bound M_I about 1.6e16 GeV stated in Sec. 2, and the n=1 benchmark eta=989.7*M_Pl^4 is dimensionally inconsistent with Eq. (78) given the definition eta'=eta/M_Pl^6. These are fixable issues, but they show that the showcased scenarios require additional scrutiny. I found no reason to suspect bad faith; the paper is explicit about its limitations. The reader's CONDITIONAL verdict remains appropriate, and a concrete B-violating completion would settle whether the mechanism is physically realizable.","tokens_in":19199,"tokens_out":27956,"duration_ms":276602,"concrete_test":"Pick a concrete B-L-violating process, e.g., lepton-number-violating scattering mediated by a heavy Majorana neutrino of mass M_N, and impose the decoupling condition Gamma_B(T_D)=H(T_D) at the GR benchmark T_D=1.48e15 GeV and M_*=2e16 GeV. Check whether the implied M_N satisfies T_D<M_N<M_* and is compatible with neutrino-mass and proton-decay bounds. If no physically allowed M_N exists, the showcased parameter point cannot be realized by a complete mechanism; if one does, the operator's missing B-violating sector is at least completable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 introduces L_int = epsilon/M_*^8 * partial_mu(T^2) J^mu_B (Eq. 2) with the statement that it 'may arise' in braneworld, string-inspired, or quantum-gravity settings, but no explicit UV construction is provided. The operator is dimension 12, so its coefficient 1/M_*^8 is a free parameter, and the paper admits in Sec. 5 that the B or B-L-violating interactions needed to convert a CPT-biased equilibrium into a net baryon number are left to future work. In gravitational baryogenesis the partial(T^2) coupling does not violate B by itself; it only biases processes that do. Consequently, the computed n_b/s (e.g., Eqs. 70, 71, 78) is not a prediction: for fixed M_*, n_b/s scales as a high power of T_D (T_D^9 in GR), so the observed value 8.8e-11 can be reproduced by a continuum of (M_*, T_D) choices. This is a fitting exercise, not a derivation. The conditionality is sharpened by internal tensions: the n=1/2 success with M_*=1e17 GeV uses T_D=2e16 GeV, exceeding the paper's own upper bound M_I about 1.6e16 GeV from Sec. 2; and the n=1 benchmark eta=989.7*M_Pl^4 is dimensionally incompatible with Eq. (78), where eta must be dimensionless by the definition eta'=eta/M_Pl^6. These issues do not make the phenomenology wrong, but they mean the central claim is not yet supported as an explanation of the baryon asymmetry.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new effective interaction term L_int = (ε/M_*^8) ∂_μ(T^2) J^μ_B for gravitational baryogenesis, where T^2 = T_{μν}T^{μν}, and studies its consequences in General Relativity and in f(R,T^2) = R + η'(T^2)^n gravity for n = 1/2 and n = 1. It derives the cosmological background for these models, computes the resulting baryon-to-entropy ratio n_b/s for the new ∂(T^2) coupling and for the standard ∂R coupling, and checks Big Bang Nucleosynthesis constraints. The paper reports that suitable choices of the cutoff M_*, decoupling temperature T_D, and model parameter η reproduce the observed n_b/s ≈ 8.8×10^{-11} in GR, in the n = 1/2 model, and in the n = 1 model.","tokens_in":19647,"tokens_out":6520,"duration_ms":67081,"significance":"The paper is transparent in its presentation of the f(R,T^2) cosmological equations and includes an explicit BBN constraint analysis, which is a useful addition to the energy-momentum-squared gravity literature. The central proposal, however, is a conditional feasibility argument rather than a derivation: the operator (2) is assumed without a UV construction, the B/B−L violating processes are not specified, and the reported agreement with the observed asymmetry is obtained by scanning the free parameters M_*, T_D, and η. Since n_b/s scales as a high power of T_D (T_D^9 in the GR and n=1/2 cases, T_D^11 in the n=1 case), the quoted 'successes' are fits rather than independent predictions. The paper could be publishable as a parameter-space scan showing that such an operator can in principle bias baryogenesis to the observed level, but the present framing overstates the explanatory power of the mechanism.","major_comments":[{"comment":"The central claim rests on an unproven effective operator. The paper states that Eq. (2) 'may arise' in braneworld, string-inspired, or quantum-gravity settings, but no derivation from any of those frameworks is provided, and §5 explicitly leaves the B or B−L violating interaction to future work. Since the ∂_μ(T^2) coupling alone does not violate baryon number, the computed n_b/s in Eqs. (70), (71), and (78) is not a prediction: for fixed M_* and η, the asymmetry scales as T_D^9 or T_D^11, so the observed value 8.8×10^{-11} can be reproduced by a continuum of (M_*, T_D) choices. The authors should either reframe the results as a feasibility scan with explicit statements about the free parameters, or provide a concrete UV or particle-physics mechanism that fixes T_D and M_*.","section":"§2, Eq. (2); §5"},{"comment":"One of the two 'successful' benchmarks in the n = 1/2 case uses T_D = 2×10^{16} GeV with M_* = 10^{17} GeV. This decoupling temperature exceeds the paper's own upper bound M_I ≈ 1.6×10^{16} GeV, stated in §2 as the maximum energy scale based on tensor-mode constraints. The authors should either provide a benchmark that satisfies T_D ≤ M_I, or explicitly justify why the decoupling temperature can exceed the inflationary upper bound. As written, the quoted success lies outside the admitted domain of validity.","section":"§4.2(ii), Eqs. (70)-(71)"},{"comment":"The quoted benchmark η = 989.7 × M_Pl^4 is dimensionally incompatible with the definition of η in the Introduction, where η' = η M_Pl^{−6} for n = 1 implies that η is dimensionless. Equation (78) contains η^{1/2} and is dimensionally consistent only if η is dimensionless; inserting a dimensionful η = 989.7 M_Pl^4 makes the expression dimensionally inconsistent. In addition, the statement in the text that η' = η/M_Pl^4 for n = 1 contradicts the Introduction's η' = η M_Pl^{−6}. These inconsistencies undermine the n = 1 success claim and must be corrected.","section":"§4.2(iii), Eq. (78)"}],"minor_comments":[{"comment":"The paper contains several typographical errors, including 'cuto ff' in the abstract and 'is is' at the start of §4; these should be corrected in a careful proofreading pass.","section":"Abstract, §4"},{"comment":"The notation η is used both for the f(R,T^2) coupling parameter and for the baryon-to-entropy ratio η_s in the introduction; the authors should consistently distinguish the two to avoid confusion.","section":"§1, §4.2(iii)"},{"comment":"The derivation of ˙T^2 = −η'^{-1} t^{-3} from ρ(t) ≃ sqrt(3/(8η')) t^{-1} is not shown in detail, and the numerical coefficient in Eq. (78) is not transparent; a step-by-step derivation would help the reader verify the result.","section":"§4.2(iii), Eqs. (74)-(78)"},{"comment":"The term proportional to ∂^2 L_m/∂g^{μν}∂g^{αβ} is set to zero 'to mitigate divergences'; the authors should acknowledge the loss of generality that this choice entails, since the BBN constraints derived later depend on this assumption.","section":"§3.1, Eq. (14)"},{"comment":"The BBN bounds are quoted both for η and for η'; the authors should explicitly state the units of η' in each case (e.g., GeV^{−2} for n = 1/2) to avoid dimensional ambiguity.","section":"§4.1(i), Eq. (52)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript shows competence in the f(R,T^2) formalism and the BBN analysis, but the central mechanism is not yet established. The referee should require the authors to reframe the work as a feasibility study with clearly labeled free parameters, to remove the internal inconsistencies (T_D vs. M_I and the dimensionful η), and to soften the claim that the asymmetry has been 'solved'. If these revisions are made, the paper could be a useful contribution to gravitational baryogenesis phenomenology."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this paper has one genuinely new trick and a pile of honest but conditional phenomenology. The new trick is the interaction ∂_μ(T^2)J^μ_B. Because T^2 = ρ^2(1+3w^2), the coupling does not vanish in the radiation epoch (w=1/3), unlike the standard ∂_μ R J^μ term with its (1-3w) factor. That is a real observation, and I hadn't seen it in the baryogenesis literature.\n\nCredit where due: the authors derive explicit asymmetry formulas and BBN bounds for f(R,T^2)=R+η'(T^2)^n with n=1/2 and n=1, and they are transparent about the structure of the calculation. They also state plainly that the actual B or B-L violating interactions are left to future work, which is the right thing to say given the mechanism only biases pre-existing violation.\n\nNow the soft spots. The operator is dimension 12, introduced ad hoc with a 'may arise' gesture toward braneworld/string/quantum-gravity effective actions and no explicit construction. M_* and T_D are free parameters, and because n_b/s scales like T_D^9 in the GR case, the reported successes are fits, not predictions. The paper sometimes overstates this: 'the model predicts successful baryogenesis' is really 'for these chosen numbers the formula returns the observed value.'\n\nInternal consistency has two warts. One success case uses T_D=2×10^16 GeV, above the paper's own quoted inflationary bound M_I ≈1.6×10^16 GeV. And in the n=1 model, the benchmark η=989.7 M_Pl^4 appears dimensionally incompatible with the definition η'=η/M_Pl^6 and with Eq. (78), which requires η dimensionless. This is fixable but should have been caught.\n\nBottom line: the central mechanism is conditional rather than explanatory. That is not a reason to desk-reject; the (1+3w^2) observation and the explicit formulas are useful for the modified-gravity baryogenesis community. Send it to peer review, with the expectation that the referee will ask the authors to fix the dimensional slips, restate the results as parameter fits, and either produce a UV motivation or clearly frame the work as a phenomenological template.","headline":"A genuinely new (1+3w^2) trick for gravitational baryogenesis, wrapped in a parameter-fitting exercise with internal consistency slips; worth refereeing, but not yet an explanation.","tokens_in":20194,"tokens_out":6545,"would_cite":false,"duration_ms":60954,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.Kd","98.80.Cq"],"model":"deepseek-v4-flash","headline":"A new 'matter-squared' coupling to the baryon current can generate the observed baryon asymmetry, the paper argues.","keywords":["baryogenesis","gravitational baryogenesis","f(R,T^2) gravity","energy-momentum squared gravity","matter-antimatter asymmetry","Big Bang nucleosynthesis","early Universe cosmology","modified gravity"],"falsifier":"Derive the coefficient of $\\partial_\\mu(\\mathcal{T}^2)J^\\mu_B$ in a concrete braneworld, string, or quantum-gravity effective action: if it is absent or corresponds to $M_*$ far above $10^{17}$ GeV, the predicted asymmetry cannot reach $8.8\\times10^{-11}$. A second, independent check is to measure the expansion history near $T\\sim10^{15}$ GeV through primordial gravitational waves; a background inconsistent with the modified Friedmann equations used here would close the viable window.","tokens_in":18988,"feed_emoji":"🌌","tokens_out":11276,"duration_ms":102955,"temperature":0.7,"pith_summary":"The paper aims to establish that higher-order matter contributions in gravity can produce the observed excess of matter over antimatter in the Universe. It proposes the interaction $\\mathcal{L}_{\\mathrm{int}} = (\\epsilon/M_*^8)\\,\\partial_\\mu(\\mathcal{T}^2)J^\\mu_B$, in which the derivative of the squared energy-momentum tensor couples to the baryon current, and it computes the resulting baryon-to-entropy ratio in General Relativity and in $f(R,\\mathcal{T}^2)$ gravity with $n=1/2$ and $n=1$. The authors find that this term alone can reproduce $n_b/s \\simeq 8.8\\times10^{-11}$ in all these settings for suitable cutoff scales $M_*$ and decoupling temperatures $T_D$, and they check each scenario against Big Bang Nucleosynthesis bounds. The broader point is that a small departure from General Relativity, or even no departure at all, may suffice for baryogenesis.","feed_headline":"'Matter-squared' term reproduces the observed baryon asymmetry","feed_subtitle":"The asymmetry would arise from gravity's matter-squared coupling, not from new physics at the weak scale.","key_machinery":"The load-bearing object is the dimension-12 effective operator $\\mathcal{L}_{\\mathrm{int}} = (\\epsilon/M_*^8)\\,\\partial_\\mu(\\mathcal{T}^2)J^\\mu_B$, built from the scalar $\\mathcal{T}^2 \\equiv T_{\\mu\\nu}T^{\\mu\\nu} = \\rho^2(1+3w^2)$ for a perfect fluid. It is the $\\mathcal{T}^2$ analogue of the standard gravitational baryogenesis coupling $\\partial_\\mu R J^\\mu$, but with two advantages used throughout the paper: it does not vanish when $w=1/3$, and its time derivative scales as a high power of temperature, so the asymmetry inherits a strong $T_D^9/M_*^8$ dependence. The machinery also includes the modified Friedmann, acceleration, and continuity equations of $f(R,\\mathcal{T}^2)$ gravity with $f = R + \\eta'(\\mathcal{T}^2)^n$, whose functions $F_{\\mathrm{Frd}}(n,w)$ and $F_{\\mathrm{Acc}}(n,w)$ control how $\\mathcal{T}^2$ alters the expansion rate and hence the decoupling condition.","core_discovery":"The central claim is that the operator $\\mathcal{L}_{\\mathrm{int}} = (\\epsilon/M_*^8)\\,\\partial_\\mu(\\mathcal{T}^2)J^\\mu_B$ generates a net baryon asymmetry through the time derivative of $\\mathcal{T}^2$ evaluated at decoupling, giving $n_b/s \\simeq -(15g_b\\epsilon/4\\pi^2 g_*)\\,(\\dot{\\mathcal{T}}^2/M_*^8 T)\\big|_{T_D}$. For a perfect fluid $\\mathcal{T}^2 = \\rho^2(1+3w^2)$, so the coupling survives in the radiation epoch ($w=1/3$), where the original Ricci-scalar coupling $\\partial_\\mu R J^\\mu$ vanishes, and it is amplified by the quadratic density dependence. In pure General Relativity the paper obtains $n_b/s\\simeq8.79\\times10^{-11}$ for $M_*\\simeq2\\times10^{16}\\,\\mathrm{GeV}$ and $T_D\\simeq1.48\\times10^{15}\\,\\mathrm{GeV}$, and it reports successful baryogenesis in the $f(R,\\mathcal{T}^2)$ models with $n=1/2$ and $n=1$ as well. Big Bang Nucleosynthesis bounds are used to restrict the model parameter $\\eta$, leaving narrower but non-empty viable windows.","pith_inferences":["If operator (2) is real, the same construction could be tried with the lepton or $B-L$ current, and the required scales $M_*\\sim10^{16}$ GeV point toward GUT-scale or Planck-scale physics; identifying the explicit B/L-violating process is the natural next step the paper leaves open.","The strong scaling $n_b/s\\propto T_D^9/M_*^8$ makes the mechanism falsifiable by independent early-Universe probes, such as primordial gravitational-wave backgrounds, that constrain either $M_*$ or the expansion history at temperatures near $10^{15}$ GeV.","The constraint that forces $\\rho\\propto a^{-4}$ relates $\\eta$ to the equation of state; since it vanishes for $w=1/3$, a radiation-only universe would reduce to GR, suggesting that stiff-matter or other non-radiation phases, which the paper notes can enhance gravitational-wave signals, may be the natural arena for this baryogenesis channel.","A direct derivation of (2) in a particular braneworld or string-inspired model would turn the parametric success into a prediction, because the effective cutoff $M_*$ would no longer be free."],"forward_implications":["The $\\mathcal{T}^2$ coupling remains active during the radiation epoch, so gravitational baryogenesis no longer suffers the vanishing $(1-3w)$ factor that suppresses the original Ricci-scalar mechanism.","A successful asymmetry can be obtained in pure General Relativity with $M_*\\simeq2\\times10^{16}$ GeV and $T_D\\simeq1.48\\times10^{15}$ GeV, so the mechanism does not require modified gravity.","In the $n=1/2$ $f(R,\\mathcal{T}^2)$ model, Big Bang Nucleosynthesis-compatible values of $\\eta$ still yield the observed asymmetry when the cutoff is pushed to $M_*\\simeq10^{17}$ GeV and $T_D\\simeq2\\times10^{16}$ GeV, so only modest departures from GR are needed.","In the $n=1$ model, the original $\\dot{R}$ coupling generates no asymmetry because $R=0$ in the radiation epoch, but the $\\mathcal{T}^2$ coupling works and is the most sensitive to the model parameter $\\eta$.","The dimension-12 interaction is suppressed by $(T_{\\mathrm{BBN}}/M_*)^8$, so it is harmless at nucleosynthesis even for cutoff scales as low as $10^7$ GeV."],"supporting_citations":[{"why":"supplies the original gravitational baryogenesis mechanism and the $\\partial_\\mu R J^\\mu$ interaction that the new term generalizes.","marker":"[13]"},{"why":"introduced $f(R,T^2)$ gravity and the $n=1/2$ model adopted here.","marker":"[43]"},{"why":"provides the modified Friedmann, acceleration, and continuity equations used for the cosmological background.","marker":"[44]"},{"why":"provides the energy-momentum-squared gravity model for the $n=1$ case.","marker":"[45]"},{"why":"imposes the ghost-free stability bound $2>\\eta>0$ that restricts the parameter range.","marker":"[50]"},{"why":"fixes the effective relativistic degrees of freedom $g_*$ and the radiation energy density normalization used in the asymmetry formulas.","marker":"[61]"},{"why":"supplies the Big Bang Nucleosynthesis freeze-out method linking $\\delta H$ to $\\delta T_{\\mathrm{freeze}}$.","marker":"[80]"},{"why":"supplies the primordial helium abundance constraint $\\delta Y_p<10^{-4}$ underlying the BBN limits.","marker":"[83]"}],"fun_headline_variants":["Matter-squared term in gravity reproduces baryon excess","Gravity's T^2 coupling yields observed matter asymmetry","Baryogenesis from higher-order matter in gravity","Matter-squared gravity explains matter-antimatter asymmetry","New gravity term sets baryon asymmetry without weak-scale physics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mechanism rests on the assumption that the interaction (2) exists in the early-Universe effective theory with a cutoff $M_*$ in the range needed; the paper motivates it by analogy with braneworld, string-inspired, and quantum-gravity effective actions but does not derive it from any of them, and the actual $B$ or $B-L$ violating process is left unspecified.","fun_headline_variants_meta":{"raw":{"variants":["Matter-squared term in gravity reproduces baryon excess","Gravity's T^2 coupling yields observed matter asymmetry","Baryogenesis from higher-order matter in gravity","Matter-squared gravity explains matter-antimatter asymmetry","New gravity term sets baryon asymmetry without weak-scale physics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000299,"raw_usage":{"total_tokens":1769,"prompt_tokens":1024,"completion_tokens":745,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":640,"completion_tokens_details":{"reasoning_tokens":667}},"tokens_in":640,"tokens_out":745,"duration_ms":8291,"temperature":1.0,"reasoning_tokens":667,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:01:24.380742+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Derive the coefficient of $\\partial_\\mu(\\mathcal{T}^2)J^\\mu_B$ in a concrete braneworld, string, or quantum-gravity effective action: if it is absent or corresponds to $M_*$ far above $10^{17}$ GeV, the predicted asymmetry cannot reach $8.8\\times10^{-11}$. A second, independent check is to measure the expansion history near $T\\sim10^{15}$ GeV through primordial gravitational waves; a background inconsistent with the modified Friedmann equations used here would close the viable window.","supporting_citations":[],"review_version":1}