{"id":"8d8b8057-e625-4846-b483-ae41d76da1d9","arxiv_id":"2505.04056","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":8,"one_line_summary":"In asymmetric dark matter with self-annihilation, the shear-dominated universe permits a higher self-annihilation cross section and a lower wino mass limit than standard cosmology, while the Gauss-Bonnet braneworld does the opposite.","lead":"This paper calculates how dark matter particle-antiparticle asymmetry survives self-annihilation in two non-standard early-universe cosmologies: a shear-dominated universe and a Gauss-Bonnet braneworld. It finds that faster early expansion allows larger self-annihilation cross sections and shifts the lower bound on the wino mass relative to standard cosmology.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quoted cross-section limits are not reproducible as written: Fig. 1's g*=20 is inconsistent with the xfw=5.86 check in Sec. 2.2, which only holds for g*=90; the qualitative ordering survives, but the numerical constraints need recomputation.","rationale":"The central physical claim is robust: a faster-expanding shear universe freezes out the asymmetry washout earlier, so a larger self-annihilation cross section is allowed, and the opposite holds for the Gauss-Bonnet braneworld. This follows directly from the A_s,g factors in Eq. (11). My concern is not with this ordering but with the specific quantitative limits quoted in the abstract and Eqs. (16)-(17). The g* inconsistency is a concrete, checkable flaw: the same a=4.37e-16 produces xfw=5.86 only for g*=90, whereas Fig. 1 states g*=20. This means either the figure, the analytic check, or the quoted limits is wrong. Because the paper provides no code or numeric tables, readers cannot resolve it. The secondary issue about the initial condition is inherited from the ADM self-annihilation literature and does not affect the comparative claim, but it should be stated more carefully. I do not think the abrupt-transition assumption is the most load-bearing concern: for xe=130,630 and xt=50,200, the modified phase extends past the washout freeze-out at xfw around 6, so the transition occurs after the epoch of interest. Overall, the paper's contribution is a modest, physically expected extension; it should remain conditionally accepted pending a corrected and reproducible set of numbers.","tokens_in":9933,"tokens_out":29441,"duration_ms":292890,"concrete_test":"Recompute the R(a) curves in Fig. 1 for the standard cosmology and xe=130 using Eq. (11) with both g*=20 and g*=90, and extract a at R=0.1; then substitute each a into Eq. (12) to compute xfw and compare with the stated xfw=5.86 and R_approx=0.11. If the quoted a values reproduce only with g*=90, the Fig. 1 caption is wrong; if they reproduce with g*=20, the approximate Eq. (12) is not consistent with the numerical solution and the validation claim in Sec. 2.2 fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline numbers (standard a <= 4.37e-16, shear a <= 8.84e-15 at xe=130) come from Fig. 1, whose caption fixes g*=20. Yet the internal consistency check in Sec. 2.2 quotes xfw=5.86 at a=4.37e-16 for standard cosmology. Solving Eq. (12) with g*=90 gives xfw about 5.86, but with g*=20 gives xfw about 6.6; the latter implies R_approx about 0.06, not the claimed 0.10-0.12. Thus the numerical and approximate solutions are claimed to be in strong consistency only if Fig. 1's g* is a typo and should be 90. Since the M2 limits in Eqs. (16)-(17) are derived from the same cross-section bounds, they inherit this ambiguity. A secondary concern is that Delta_-in is set to the equilibrium value Delta_-eq(1) with mu/T=1e-9; with this initialization, R about 0.1 reflects the Boltzmann suppression of the equilibrium asymmetry at freeze-out rather than the survival fraction of an independently specified asymmetry, so the absolute limits are convention-dependent. Neither issue reverses the central ordering (shear-enhanced H allows larger a, GB-weakened H requires smaller a), which follows directly from A_s>1>A_g in Eq. (11).","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the evolution of the relic density of asymmetric dark matter (ADM) when self-annihilation processes are included, in two non-standard cosmological scenarios: a shear-dominated Bianchi type I universe and a Gauss-Bonnet braneworld. Using the Hubble rates of Eqs. (2) and (3), the authors derive the asymmetry evolution equation (11) and solve it both numerically and with a sudden-freeze-out approximation. They find that, for the same final-to-initial asymmetry ratio R, the shear-dominated universe (enhanced H) permits a larger self-annihilation cross section than standard cosmology, while the Gauss-Bonnet braneworld (weakened H) requires a smaller one. They then translate these cross-section bounds into lower limits on the wino mass M2 for sneutrino and higgsino ADM, quoted in Eqs. (16) and (17). The central qualitative ordering is physically transparent and robust to the numerical details.","tokens_in":10301,"tokens_out":4824,"duration_ms":47191,"significance":"If the numerical results are correct, the paper provides concrete, model-dependent constraints: for R > 0.1, the standard-cosmology s-wave limit a ≲ 4.37e-16 is raised to a ≲ 8.84e-15 (xe=130) and a ≲ 4.20e-14 (xe=630) in the shear-dominated universe, and lowered to a ≲ 1.14e-16 (xt=50) and a ≲ 4.57e-17 (xt=200) in the Gauss-Bonnet braneworld. The corresponding M2 lower limits in Eqs. (16)-(17) extend earlier ADM analyses to these non-standard histories. The derivation of Eq. (11) from the coupled Boltzmann equations is clear, and the approximate solution of Eqs. (12)-(13) is a useful check once the g* inconsistency is resolved. The main strength is that the ordering of limits follows directly from A_s > 1 > A_g in Eq. (11), so the qualitative claim does not depend on disputed numerical details.","major_comments":[{"comment":"The numerical constraints quoted in the text are not reproducible as written because Fig. 1 fixes g* = 20 while the consistency check in Sec. 2.2 (xfw = 5.86 at a = 4.37e-16 and Rapp = 0.11) and Fig. 3 use g* = 90. With g* = 20, solving Eq. (12) at the quoted standard value a = 4.37e-16 gives xfw ≈ 6.6 rather than 5.86, and Eq. (13) then gives Rapp ≈ 0.06, not the claimed 0.10-0.12. Since the headline limits and the M2 bounds in Eqs. (16)-(17) are read from Fig. 1, this inconsistency affects all quantitative results and must be resolved by recomputing with one stated g* and correcting the figure captions and the consistency check accordingly.","section":"Sec. 2.2, Figs. 1-3"},{"comment":"The initialization Delta_-in = Delta_-eq(1) with mu/T = 1e-9 means that the reported ratio R is dominated by the Boltzmann suppression of the equilibrium asymmetry at freeze-out rather than by the survival of an independently specified input asymmetry. Consequently, the absolute limits such as a ≲ 4.37e-16 are convention-dependent; a different choice of mu/T or xin would shift them, even though the relative ordering between cosmologies would remain. The authors should either justify this initialization as physical or demonstrate how the cross-section and M2 limits change with the initial condition.","section":"Sec. 2.2, initial conditions"},{"comment":"The choice R > 0.1 is arbitrary and no sensitivity study is provided. All quoted limits in Sec. 2.2 and Eqs. (16)-(17) use this threshold, so the quantitative constraints depend on an unmotivated criterion. The authors should motivate a physically meaningful minimum R (for example, one consistent with the observed DM asymmetry) or present R-dependent results so that readers can assess how the limits vary with the assumed preserved fraction.","section":"Sec. 2.2, benchmark R > 0.1"},{"comment":"The modified Hubble rates are assumed to hold throughout the entire asymmetry-washout epoch, with an abrupt transition to the standard radiation-dominated rate at x_e or x_t. If the shear-dominated or Gauss-Bonnet phase ended before the washout process froze out, or if the transition were gradual, the cross-section and M2 limits would shift. The paper should state the assumed duration of the modified phase relative to the freeze-out point and, ideally, test sensitivity to the transition; otherwise the quoted upper limits are conditional on this unquantified assumption.","section":"Eqs. (2)-(3) and Sec. 2.2"}],"minor_comments":[{"comment":"The quantity T_e appearing in x_e = sqrt(g*/g_e^*) m/T_e is not defined precisely; it should be stated explicitly that T_e is the temperature at which the shear energy density equals the radiation energy density.","section":"Below Eq. (2)"},{"comment":"There is a typo in the sentence beginning 'When R = 0.1, we fnd...' ('fnd' should be 'find'). The claim that the constraints on a and 6b/x 'must be same' should be phrased more precisely: at fixed R, Eq. (13) fixes xfw and Eq. (12) then fixes the combination ⟨sigma v⟩, which for s-wave and p-wave gives equal effective cross sections at that xfw.","section":"Sec. 2.2, near Eq. (13)"},{"comment":"The M2 limits in Eqs. (16)-(17) are quoted with little explanation of the input spectrum; specifying M1 = M2/2 and tan^2 theta_W = 0.3 in the text is helpful, but the resulting values would benefit from a brief comparison with existing constraints in the standard cosmological case to anchor the non-standard results.","section":"Fig. 4 and Eqs. (16)-(17)"}],"recommendation":"major_revision","confidential_remarks":"The g* inconsistency (Fig. 1 at g* = 20 vs. the g* = 90 consistency check) is likely a typo, but it is load-bearing because all quoted numerical limits and M2 bounds are read from the affected figures. The qualitative conclusion is sound and follows directly from As > 1 > Ag, so I do not see a need to reject; however, the numerical constraints must be recomputed with a single, consistently used g* and the benchmark/convention issues should be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Nick,\n\nThis is a modest but useful quantitative paper. It takes the existing ADM self-annihilation washout formalism (Ellwanger–Mitropoulos, Liu–Iminniyaz) and applies it to shear-dominated and Gauss–Bonnet cosmologies. The central result—shear allows a larger self-annihilation cross section, Gauss–Bonnet requires a smaller one, with corresponding shifts in the wino mass lower bound—is physically sensible and cleanly argued. The Boltzmann equations are standard and well presented, and the ordering follows directly from A_s > 1 > A_g in Eq. (11). The citation pattern is fine; the one self-citation supplies the freeze-out approximation method, which is legitimate.\n\nThe soft spots are real but not disqualifying. Most serious: Fig. 1's caption fixes g* = 20, while the consistency check in Sec. 2.2 (xfw = 5.86, Rapp = 0.11 for the standard cosmology) only works with g* = 90. With g* = 20, Eq. (12) gives xfw about 6.6 and Rapp about 0.06, not the claimed 0.10–0.12. Since the M2 limits in Eqs. (16)–(17) are derived from the same cross-section bounds, they inherit the ambiguity. This is likely a typo—Figs. 2–4 use g* = 90—but as written the headline numbers need recomputation. Second, the initialization of Delta_-in at the equilibrium value with mu/T = 1e-9 means the R > 0.1 benchmark reflects the Boltzmann suppression of the equilibrium asymmetry at freeze-out, not the survival of an independently specified initial asymmetry. The absolute limits are convention-dependent, even though the qualitative ordering is robust. Minor: no code or numeric tables, the R > 0.1 choice is arbitrary, and the abrupt switch from modified to standard Hubble rate is imported without modeling a gradual transition.\n\nWho this is for: ADM model builders, wino DM phenomenologists, anyone working on non-standard early-universe cosmologies. It deserves a serious referee; the formalism is checkable and the qualitative result is solid. With the g* inconsistency fixed and the Delta_in convention acknowledged, I would be happy to see it published in a specialized journal.\n\nRecommendation: engage with it, but ask for the recomputation and a clarifying note on the convention before relying on the numbers.","headline":"Useful, clearly argued extension of ADM self-annihilation constraints to two non-standard cosmologies, but the headline cross-section limits are compromised by a g* inconsistency and a convention-dependent asymmetry benchmark.","tokens_in":10858,"tokens_out":3742,"would_cite":false,"duration_ms":33176,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Asymmetric dark matter's allowed self-annihilation cross section shifts upward in a faster-expanding shear-dominated early universe and downward in a slower Gauss-Bonnet braneworld, relative to standard cosmology.","keywords":["asymmetric dark matter","self-annihilation","relic density","shear-dominated universe","Gauss-Bonnet braneworld","Boltzmann equation","wino mass","freeze-out"],"falsifier":"A measurement or calculation that fixes the transition temperatures $T_e$ (shear) and $T_t$ (Gauss-Bonnet) relative to the freeze-out temperature $T_{fw}$ would settle the claim: if either transition occurred before freeze-out, the corresponding cross-section and $M_2$ bounds would revert to standard values rather than shift as reported.","tokens_in":9734,"feed_emoji":"⚛️","tokens_out":7648,"duration_ms":69552,"temperature":0.7,"pith_summary":"This paper asks how much asymmetric dark matter (ADM) may self-annihilate without erasing the particle–antiparticle asymmetry that sets its relic density, in two nonstandard early-universe cosmologies. It finds that a shear-dominated universe, whose Hubble rate is enhanced over the standard one, freezes out the asymmetry-washout process earlier and permits s-wave self-annihilation cross sections up to about twenty to a hundred times larger than in standard cosmology. A Gauss-Bonnet braneworld, with a weakened Hubble rate, delays freeze-out and forces cross sections several times smaller. These shifts translate into wino-mass lower limits for sneutrino and higgsino ADM that are lower in the shear case and higher in the Gauss-Bonnet case. The paper also gives approximate analytic formulas that reproduce the numerical freeze-out constraints.","feed_headline":"Expansion speed sets dark matter self-annihilation limits","feed_subtitle":"A faster early universe allows larger cross sections; a slower one demands smaller. Wino-mass bounds shift accordingly.","key_machinery":"The engine is the asymmetry variable $\\Delta_- = Y_\\chi - Y_{\\bar\\chi}$ and its evolution equation $d\\Delta_-/dx = -(\\lambda/(x^2 A_{s,g}))\\langle\\sigma_{\\chi\\chi} v\\rangle \\Delta_+^{eq}(\\Delta_- - \\Delta_-^{eq})$, where $A_s = \\sqrt{1+x_e^2/x^2}$ for the shear-dominated cosmology and $A_g = (x/x_t)^{2/3}$ for the Gauss-Bonnet braneworld modify the standard Boltzmann rate. $A_s > 1$ slows the washout of the asymmetry, while $A_g < 1$ speeds it up, producing the opposite shifts in the cross-section limits. The companion approximate formulas give the freeze-out point $x_{fw}$ and the surviving asymmetry ratio $R$ in closed form, and the paper shows these reproduce the numerical solution at $R=0.1$ with $R_{\\text{approx}} = 0.10$–$0.12$.","core_discovery":"Under the requirement that the ratio of final asymmetry to initial asymmetry satisfies $R > 0.1$, the paper reports that standard cosmology allows an s-wave self-annihilation cross section $a \\lesssim 4.37\\times10^{-16}$; the shear-dominated universe allows $a \\lesssim 8.84\\times10^{-15}$ at $x_e=130$ and $a \\lesssim 4.20\\times10^{-14}$ at $x_e=630$, while the Gauss-Bonnet braneworld only allows $a \\lesssim 1.14\\times10^{-16}$ at $x_t=50$ and $a \\lesssim 4.57\\times10^{-17}$ at $x_t=200$. For the same $R$, the lower limit on the wino mass $M_2$ in the s-wave case shifts from $M_2 \\gtrsim 3.5\\times10^6$ GeV in standard cosmology to $M_2 \\gtrsim 8.3\\times10^5$ GeV at $x_e=130$ and $M_2 \\gtrsim 3.8\\times10^5$ GeV at $x_e=630$ in the shear case, and to $M_2 \\gtrsim 7.2\\times10^6$ GeV at $x_t=50$ and $M_2 \\gtrsim 1.1\\times10^7$ GeV at $x_t=200$ in the Gauss-Bonnet case. The same direction of shift holds for p-wave self-annihilation, with similar numerical values.","pith_inferences":["If the shear-dominated or Gauss-Bonnet phase ended before the asymmetry freeze-out temperature, the quoted bounds would relax back toward standard values; the paper's numbers therefore presuppose that the modified era spans the entire wash-out epoch.","Because the direction of the shift is controlled by whether the modified Hubble rate is above or below the standard one, the same two equations can be applied to any cosmology with a known expansion history, giving a quick estimate of whether it loosens or tightens ADM self-annihilation bounds.","The $M_2$ bounds around $10^5$–$10^7$ GeV fall in a range where electroweak gaugino searches at colliders could eventually probe them; this is not tested in the paper, but follows from the mass scales it constrains."],"forward_implications":["For $R>0.1$, the shear-dominated universe raises the maximum s-wave self-annihilation cross section by a factor of about 20 at $x_e=130$ and about 96 at $x_e=630$ relative to standard cosmology.","The Gauss-Bonnet braneworld lowers the same limit by factors of about 3.8 at $x_t=50$ and about 9.6 at $x_t=200$.","The wino-mass lower limit for sneutrino and higgsino ADM drops by roughly an order of magnitude in the strongest shear case and rises by roughly a factor of three in the strongest Gauss-Bonnet case.","The approximate freeze-out equations give essentially the same cross-section bounds as the full numerical integration, so the constraints do not depend on solving the complete Boltzmann system.","Both velocity-independent (s-wave) and velocity-suppressed (p-wave) self-annihilation shift in the same direction, making the qualitative conclusion insensitive to the partial-wave structure."],"supporting_citations":[{"why":"Supplies the shear-dominated Hubble expansion rate $H_s = H\\sqrt{1+x_e^2/x^2}$ that drives the enhanced freeze-out.","marker":"[26, 27]"},{"why":"Supplies the Gauss-Bonnet braneworld Hubble rate $H_g = H(x/x_t)^{2/3}$ that drives the delayed freeze-out.","marker":"[31]"},{"why":"Introduces the self-annihilation washout of the ADM asymmetry in standard cosmology, the baseline being extended.","marker":"[20]"},{"why":"Provides the method for approximating the freeze-out point $x_{fw}$ of the wash-out process used here.","marker":"[21]"},{"why":"Establishes ADM relic density in the shear-dominated universe without self-annihilation, the starting point for adding the self-annihilation term.","marker":"[28]"},{"why":"Supplies the p-wave higgsino self-annihilation cross-section formula used for the $M_2$ constraint.","marker":"[33]"},{"why":"Supplies the sneutrino ADM candidate and the self-annihilation cross-section via electroweak gaugino exchange.","marker":"[15, 16, 32]"}],"fun_headline_variants":["Expansion speed dictates dark matter annihilation limits","Shear universe raises, braneworld lowers dark matter bounds","Faster expansion expands dark matter cross-section options","Cosmic expansion rate sets dark matter self-annihilation cap","Early-universe expansion alters dark matter cross-section limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The nonstandard expansion histories are assumed to last through the entire asymmetry-washout era, switching abruptly at a single transition temperature; if the shear or Gauss-Bonnet phase ended before freeze-out, or transitioned gradually, all quoted limits would shift.","fun_headline_variants_meta":{"raw":{"variants":["Expansion speed dictates dark matter annihilation limits","Shear universe raises, braneworld lowers dark matter bounds","Faster expansion expands dark matter cross-section options","Cosmic expansion rate sets dark matter self-annihilation cap","Early-universe expansion alters dark matter cross-section limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000702,"raw_usage":{"total_tokens":3219,"prompt_tokens":1048,"completion_tokens":2171,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":664,"completion_tokens_details":{"reasoning_tokens":2094}},"tokens_in":664,"tokens_out":2171,"duration_ms":15793,"temperature":1.0,"reasoning_tokens":2094,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:39:22.610172+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement or calculation that fixes the transition temperatures $T_e$ (shear) and $T_t$ (Gauss-Bonnet) relative to the freeze-out temperature $T_{fw}$ would settle the claim: if either transition occurred before freeze-out, the corresponding cross-section and $M_2$ bounds would revert to standard values rather than shift as reported.","supporting_citations":[{"cited_title":"Dark matter relic density in Gauss-Bonnet braneworld cosmology","cited_arxiv_id":"1404.4424","evidence_quote":"Supplies the Gauss-Bonnet braneworld Hubble rate $H_g = H(x/x_t)^{2/3}$ that drives the delayed freeze-out."},{"cited_title":"Upper Bounds on Asymmetric Dark Matter Self Annihilation Cross Sections","cited_arxiv_id":"1205.0673","evidence_quote":"Introduces the self-annihilation washout of the ADM asymmetry in standard cosmology, the baseline being extended."},{"cited_title":"Constraints on Asymmetric Dark Matter Self Annihilation Cross Sections in Non-standard Cosmological Scenarios","cited_arxiv_id":"2309.09155","evidence_quote":"Provides the method for approximating the freeze-out point $x_{fw}$ of the wash-out process used here."},{"cited_title":"Asymmetric Dark Matter in the Shear--dominated Universe","cited_arxiv_id":"1604.04251","evidence_quote":"Establishes ADM relic density in the shear-dominated universe without self-annihilation, the starting point for adding the self-annihilation term."}],"review_version":1}