{"id":"03e1d6ce-b796-4c2c-8934-8c58125b8384","arxiv_id":"1908.11325","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In a U(1)_{B-L} model with a Dirac singlet dark matter fermion, relic abundance constraints fix combinations of the dark matter and gauge couplings, and freeze-in scenarios leave parameter space that FASER, Belle II, SHiP, and LDMX can probe.","lead":"This paper maps the allowed mass and coupling regions for a dark matter particle and a new B-L gauge boson in a specific extension of the Standard Model. It combines thermal freeze-out and freeze-in relic calculations with direct detection, CMB, and upcoming lifetime frontier experiment sensitivities.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Freeze-in and lifetime-frontier constraints assume zero SM kinetic mixing for Z_BL; because kinetic mixing is a free parameter in U(1)_{B-L}, the claimed reach in Fig. 8 is conditional on ε being negligible.","rationale":"I agree with the reader that the kinetic-mixing assumption is the weakest point. It is explicitly a simplification rather than a symmetry, and it matters precisely in the small-g_BL regime the paper targets. The alternative candidate concerns are the numerical inconsistencies (m_zeta >= 100 vs 200 GeV in Sec. III C versus the abstract, and 1.5 vs 2.5 TeV in the B2 threshold in Sec. IV A). These are real and should be corrected, but they are boundary-level typos: the qualitative exclusion of low-mass thermal dark matter and the two-regime freeze-in structure survive either choice. The kinetic-mixing issue, by contrast, controls every observable used to map the surviving parameter space. The paper deserves credit for giving explicit cross sections and for noting that radiative mixing is smaller than the direct coupling; that statement is correct for the loop contribution. The gap is that no symmetry in the model forbids a tree-level epsilon, so the plotted experimental reach is conditional. A single numerical scan over epsilon/g_BL would settle whether the concern changes the figures. On this basis the reader's CONDITIONAL verdict is appropriate; I do not recommend ACCEPT or REJECT.","tokens_in":18634,"tokens_out":12072,"duration_ms":118535,"concrete_test":"Recompute the key bounds with an explicit kinetic-mixing term epsilon F_{mu nu} Z_BL^{mu nu}, scanning epsilon = c g_BL with c = 0, 0.1, 1, 10: (i) derive the one-loop SM contribution to epsilon and confirm it is about 2 x 10^{-3} g_BL; (ii) recompute sigma_SI from Eq. (14) with the epsilon e g_zeta amplitude included; (iii) recompute the FASER/SHiP/Belle II/LDMX event rates adding the photon-mixing production amplitude. If for any c <~ 1 a relic-density line in Fig. 8 moves out of the experimentally open region, the assumption is load-bearing; if only the loop-induced epsilon is used, it is not.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is that Z_BL has no kinetic mixing with the SM photon or Z, introduced in Sec. I ('We also ignore mixings...') and encoded in Eq. (3). The freeze-out relic line, the freeze-in relations (18)-(23), the direct-detection bound of Eq. (14), and the FASER/SHiP/Belle II/LDMX curves in Fig. 8 all use only the direct g_BL couplings to SM fermions and g_zeta to DM. In a general U(1)_{B-L} model, the operator (epsilon/2) F_{mu nu} Z_BL^{mu nu} is not forbidden and can have a tree-level coefficient independent of g_BL. The paper's justification that mixing effects 'are loop suppressed and therefore smaller' bounds only the radiative contribution, epsilon_loop ~ (g_Y g_BL)/(16 pi^2) ~ 2 x 10^{-3} g_BL, which is indeed subdominant for the plotted couplings. But if a tree-level epsilon is present at the same order as g_BL, the DM-nucleus scattering amplitude becomes (epsilon e g_zeta)/M^2, the Z_BL production in beam dumps gains a photon-mixing channel, and the XENON1T boundary at g_BL <~ 8.9 x 10^{-7} plus the experimental sensitivity lines in Fig. 8 shift, possibly closing the open window. Thus the central freeze-in claim is not robust to an unconstrained parameter of the gauge sector.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a U(1)_{B-L} extension of the Standard Model with a Dirac fermion dark matter candidate ζ carrying an arbitrary B-L charge Q. It derives relic-density constraints for two regimes: thermal freeze-out and freeze-in. For freeze-out, it finds an approximate relation g_ζ ≈ 0.016 sqrt(m_ζ[GeV]) and a thermal-equilibrium lower bound g_BL ≥ 2.7×10^{-8} sqrt(m_ζ[GeV]); combined with direct and indirect detection constraints, this restricts the allowed masses. For freeze-in, it derives coupling relations depending on whether Z_BL is in equilibrium with the SM (case A) or not (cases B1/B2), and displays the parameter space that can be probed by FASER, SHiP, Belle II, LHCb and LDMX.","tokens_in":18976,"tokens_out":7154,"duration_ms":65323,"significance":"The paper provides a useful and mostly standard mapping of the parameter space of a well-motivated B-L portal dark matter model. The Boltzmann treatment and the cross-section formulas in the appendix are explicit, and the use of external constraints (XENON1T, CMB, AMS-02, LEP) is appropriate. The freeze-in analysis, including the sequential freeze-in discussion, is a valuable addition, and the comparison of the resulting parameter space with the reach of lifetime-frontier experiments is timely. The main results are reproducible in structure, though the paper would be strengthened by resolving internal inconsistencies and by clarifying the kinetic-mixing assumption.","major_comments":[{"comment":"The abstract states that the allowed mass regions are limited to m_ζ ≳ 200 GeV and M_ZBL ≳ 10 GeV, while Sec. III C and the caption of Fig. 4 both state that the allowed green region corresponds to m_ζ ≳ 100 GeV and M_ZBL ≳ 10 GeV. This is a direct inconsistency in the headline result of the paper; the threshold for the dark matter mass must be corrected and made uniform across the abstract, the body, and the figure.","section":"Abstract vs Sec. III C and Fig. 4"},{"comment":"The paper ignores kinetic mixing between Z_BL and the SM photon/Z, stating that 'these mixing effects are loop suppressed and therefore smaller.' This statement applies only to the radiative contribution, which is of order g_Y g_BL/(16π²) ≈ 2×10^{-3} g_BL. A tree-level kinetic mixing coefficient ε is a free parameter in U(1)_{B-L} and is not loop suppressed. If ε is comparable to g_BL, the DM-nucleon scattering cross section gains a photon-mediated contribution ∝ (ε e g_ζ)²/M_ZBL^4, and Z_BL production in beam dumps and colliders gains a photon-mixing channel. Both would shift the XENON1T bound (g_BL ≲ 8.9×10^{-7} in Sec. IV B) and the sensitivity curves in Fig. 8. The freeze-in constraints in Eqs. (18)–(23) and the plotted experimental reach are therefore conditional on ε = 0; the manuscript should either justify an approximate upper bound on ε from other constraints or explicitly present the results as the ε = 0 slice of the parameter space.","section":"Sec. I and Eq. (3)"},{"comment":"The summary bullet for Case (B2) in Sec. IV A states 'for m_ζ >~ 1.5 TeV, we find g_ζ² g_BL² ≃ 8.2×10^{-24}', whereas the abstract and the introduction state that the mass-independent relation applies for m_ζ ≳ 2.5 TeV. This is an internal inconsistency in the definition of the freeze-in cases; the threshold should be stated consistently as 2.5 TeV (or, if 1.5 TeV is intended, the abstract and introduction should be revised accordingly).","section":"Sec. IV A, after Eq. (23)"}],"minor_comments":[{"comment":"The coefficient 0.82/1.2 appearing in the freeze-in relic condition is unexplained. It would be helpful to state that it arises from the numerical integration of the Boltzmann equation and to use a more transparent notation.","section":"Eq. (18)"},{"comment":"The sentence 'unless the gauge coupling is below 10^{-10} GeV' appears to have a dimensional error; the quantity should probably be the dimensionless coupling g_BL < 10^{-10}.","section":"Sec. IV C"},{"comment":"The thermal-equilibrium condition for f fbar → Z_BL γ uses n_ZBL(T) times the cross section, whereas the production rate is more properly n_f n_γ / n_ZBL times ⟨σv⟩. The estimate is acceptable because all number densities are comparable in radiation domination, but this step could be clarified.","section":"Eq. (6) and surrounding text"},{"comment":"The notation for the Z_BL mass is inconsistent: M_ZBL in the text and equations, m_ZBL in Fig. 4 and in Sec. III C. Please unify the notation.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The note added in proof mentions the similar study in arXiv:1908.09834 without any comparison. For a journal submission, the authors should discuss the relationship and differences in their results. The internal inconsistencies (200 vs 100 GeV, 1.5 vs 2.5 TeV) should be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid, useful constraints paper and the central relic-density logic holds up. It also has two real soft spots, one numerical and one conceptual, that a referee should push on before acceptance.\n\nWhat is actually new: most B-L dark matter studies either fix the dark matter charge or stop at freeze-out. Here they keep Q arbitrary and split the freeze-in regime into A, B1, B2, including the sequential freeze-in contribution with a 2.5 TeV crossover. The compact relic formulas in Eqs. (18), (19), and (23) are useful, and the target lines for FASER, FASER2, Belle II, SHiP, and LDMX in Fig. 8 give experimental colleagues concrete curves. The appendix cross sections are explicit and standard, and the paper is honest: the note added in proof acknowledges the overlapping study by Heeba-Kahlhoefer, and the citation pattern to earlier B-L and freeze-in work is fine.\n\nThe numerical inconsistencies are real and should have been caught before v4. The abstract and the conclusions say the freeze-out lower bound is m_zeta >~ 200 GeV, but the Fig. 4 caption says 100 GeV. Similarly, the conclusions list case (B2) for m_zeta >~ 1.5 TeV while the abstract says 2.5 TeV, and 2.5 TeV is the value used elsewhere as the sequential freeze-in crossover. These are headline numbers, not typos in an unimportant plot, and the 2.5 TeV threshold is presented as one of the paper's new results. A referee should ask for one consistent set.\n\nThe larger worry is kinetic mixing. The paper says in Sec. I that mixings between Z_BL and the SM gauge bosons are ignored for simplicity, with the justification that they are loop suppressed. That argument covers only the radiative contribution. A tree-level operator (epsilon/2) F_mu nu Z_BL^{mu nu} is not forbidden in a general U(1)_{B-L} extension, and it enters direct detection, DM annihilation, and beam-dump production at the same order as the g_BL effects being plotted. So the freeze-in relations and the Fig. 8 reach are conditional on epsilon being negligibly small. This does not kill the paper, and it is a common assumption in the B-L literature, but it should be stated quantitatively in the model section rather than as an aside. If a tree-level epsilon is present at order g_BL, the XENON1T boundary and the experimental sensitivity lines shift and the open window may close or move.\n\nWho is this for: model builders working on U(1) portals and lifetime frontier proposers. I would send it to a serious referee. With the consistency issues fixed and a clearer kinetic mixing discussion, this is an accept; as is, it is a conditional accept.","headline":"Competent, useful B-L DM constraints paper with a defensible central relic-density argument, but it needs a consistency pass on headline bounds and a quantitative statement on kinetic mixing before the freeze-in reach is fully trusted.","tokens_in":19522,"tokens_out":2916,"would_cite":true,"duration_ms":31294,"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":"Thermal $B-L$ dark matter survives only above $m_\\zeta \\gtrsim 200$ GeV and $M_{Z_{BL}} \\gtrsim 10$ GeV.","keywords":["B-L gauge boson","dark matter","freeze-in","freeze-out","direct detection","lifetime frontier","U(1) extension","Dirac fermion dark matter"],"falsifier":"A future direct detection experiment could falsify the thermal window by observing a spin-independent $\\zeta$-nucleon cross section larger than the value implied by $g_{BL}g_\\zeta$ from the allowed parameter curves at a claimed $m_\\zeta$; equivalently, a displaced-vertex search like FASER finding a $Z_{BL}$ with parameters in the region Fig. 8 excludes, or a confirmed thermal relic with $m_\\zeta<200$ GeV in this exact model, would disprove the paper's central bounds.","tokens_in":18395,"feed_emoji":"🔭","tokens_out":13825,"duration_ms":114413,"temperature":0.7,"pith_summary":"The paper asks whether a dark matter particle that carries a new $B-L$ gauge charge can account for the observed relic abundance when the mediating gauge boson is light and weakly coupled. For thermal freeze-out, reproducing the observed relic density forces the dark-matter coupling to $g_\\zeta \\equiv g_{BL}Q \\simeq 0.016\\sqrt{m_\\zeta[{\\rm GeV}]}$, and thermal equilibrium with the Standard Model plasma requires $g_{BL} \\gtrsim 2.7\\times10^{-8}\\sqrt{m_\\zeta[{\\rm GeV}]}$. These conditions, combined with direct-detection, CMB, and cosmic-ray constraints, leave only $m_\\zeta \\gtrsim 200$ GeV and $M_{Z_{BL}} \\gtrsim 10$ GeV. For the weaker-coupling freeze-in regime the same relic condition becomes a set of compact algebraic constraints on $g_\\zeta$ and $g_{BL}$, which translate into definite predictions for lifetime-frontier experiments such as FASER, Belle II, SHiP, and LDMX.","feed_headline":"200 GeV sets the B-L dark matter floor","feed_subtitle":"Thermal freeze-out rules out lighter dark matter; weaker couplings open the freeze-in window.","key_machinery":"The load-bearing object is the $Z_{BL}$ gauge boson with two distinct couplings: $g_{BL}$, common to all Standard Model fermions through their $B-L$ charges, and $g_\\zeta = Qg_{BL}$, the dark fermion's coupling. The analysis runs on the Boltzmann equation for the DM yield, Eq. (8), with thermally averaged annihilation cross sections for $\\zeta\\bar\\zeta \\to f\\bar f$ (s-channel $Z_{BL}$ exchange) and for $\\zeta\\bar\\zeta \\to Z_{BL}Z_{BL}$ (t/u-channel exchange), listed in the appendix. The observed relic density $\\Omega_{DM}h^2 = 0.12$ converts the Boltzmann solution into the compact coupling relations quoted above, while the equilibrium conditions in Eqs. (5) and (7) decide whether the freeze-out or freeze-in regime applies. In the freeze-in case the same machinery is used with a vanishing initial DM abundance, and Eq. (23) implements the sequential freeze-in of $\\zeta$ from an intermediate $Z_{BL}$ population.","core_discovery":"The paper's central claim is that the four-parameter $B-L$ model with masses $m_\\zeta$, $M_{Z_{BL}}$ and couplings $g_{BL}$, $g_\\zeta$ collapses onto narrow, calculable curves when the dark matter is required to match the observed relic abundance. In the thermal freeze-out case the relic condition fixes $g_\\zeta \\simeq 0.016\\sqrt{m_\\zeta[{\\rm GeV}]}$ for $M_{Z_{BL}}^2 \\ll m_\\zeta^2$, and equilibrium with the Standard Model plasma requires $g_{BL} \\gtrsim 2.7\\times10^{-8}\\sqrt{m_\\zeta[{\\rm GeV}]}$. Direct detection then imposes a floor on $M_{Z_{BL}}$, while CMB and AMS-02 constraints rule out $\\zeta$ masses below roughly 200 GeV. In the freeze-in regime the DM starts with zero abundance and is produced either from the thermal plasma or through an intermediate $Z_{BL}$ population, yielding three distinct conditions: case (A) with $g_\\zeta^2 g_{BL}^2 + (0.82/1.2)g_\\zeta^4 \\simeq 8.2\\times10^{-24}$; case (B1) with $g_\\zeta^2 g_{BL}^2 \\simeq 8.2\\times10^{-24}(m_\\zeta/2.5\\,{\\rm TeV})$ for $m_\\zeta \\lesssim 2.5$ TeV through sequential freeze-in; and case (B2) with the same mass-independent product for heavier dark matter.","pith_inferences":["Beyond the paper: if kinetic mixing between $Z_{BL}$ and the SM photon or $Z$ boson is added, the DM annihilation and scattering cross sections acquire mixing-angle-dependent terms, which would reopen parts of the low-mass parameter space the paper closes; measuring the mixing would require a dedicated two-mediator analysis.","Beyond the paper: because the DM charge $Q$ is arbitrary and enters only through $g_\\zeta = Qg_{BL}$, the derived relic-density constraints apply to any $U(1)'$ portal with a vector-like dark fermion, so the quantitative results here are a template for a broader class of models.","Beyond the paper: the sharp transition between direct freeze-in (case B2) and sequential freeze-in (case B1) at $m_\\zeta \\simeq 2.5$ TeV could be probed by searching for a turn in the coupling-mass relation; a future measurement of the DM mass and its annihilation cross section would indicate which production mechanism operated.","Beyond the paper: improved CMB bounds on energy injection from DM annihilation would strengthen the indirect constraints used here and could push the thermal freeze-out mass floor above 200 GeV, narrowing the window further."],"forward_implications":["Thermal $B-L$ dark matter cannot be light: the combined constraints force $m_\\zeta \\gtrsim 200$ GeV and $M_{Z_{BL}} \\gtrsim 10$ GeV, so any low-mass signal in this model must come from the freeze-in regime.","In freeze-in case A the relic condition fixes a definite combination of couplings, $g_\\zeta^2 g_{BL}^2 + (0.82/1.2)g_\\zeta^4 \\simeq 8.2\\times10^{-24}$, making the model predictive rather than merely constrained.","Sequential freeze-in dominates for $m_\\zeta \\lesssim 2.5$ TeV in case B, so the coupling product scales with $m_\\zeta$; this is a distinctive signature that separates the two freeze-in production routes.","Lifetime-frontier experiments (FASER, FASER2, Belle II, SHiP, LDMX) are projected to reach part of the surviving parameter space, so the model's freeze-in window can be experimentally tested rather than remaining purely theoretical.","For $M_{Z_{BL}} \\lesssim 50$ MeV, the direct detection bound is satisfied only for $g_{BL}g_\\zeta \\lesssim 1.5\\times10^{-12}$, which effectively closes the low-mass region to thermal production."],"supporting_citations":[{"why":"Defines the freeze-in production mechanism that the paper uses for all small-coupling cases.","marker":"[20]"},{"why":"Provide the standard asymptotic formula for thermal relic abundance used to convert the Boltzmann solution into Eq. (12).","marker":"[37, 38]"},{"why":"Sets the observed dark matter abundance $\\Omega_{DM}h^2 = 0.12$ that all relic constraints are matched to.","marker":"[39]"},{"why":"Gives the XENON1T upper bound on spin-independent scattering used in Fig. 3 and Fig. 4.","marker":"[40]"},{"why":"Supplies the CMB and AMS-02 indirect detection bounds that exclude the low-mass thermal region.","marker":"[30]"},{"why":"Provides the spin-independent DM-nucleon cross section formula quoted in Eq. (14).","marker":"[36]"},{"why":"Identifies sequential freeze-in, which the paper applies in Case B1.","marker":"[49]"},{"why":"Defines the FASER geometry and sensitivity used for the lifetime-frontier projections in Fig. 8.","marker":"[50]"},{"why":"Earlier freeze-in analysis with a Majorana DM and $g_\\zeta=g_{BL}$; the paper states its results agree with this Dirac generalization.","marker":"[24]"},{"why":"Provides CMB-based lower bounds on dark matter mass and dark-photon portal constraints, used to take $m_\\zeta \\geq 1$ GeV.","marker":"[29]"}],"fun_headline_variants":["B-L DM floor at 200 GeV; freeze-in widens probe","Freeze-out sets B-L DM floor; freeze-in maps lighter zone","B-L DM: 200 GeV floor, freeze-in window","B-L dark matter: freeze-out floor, freeze-in frontier","B-L DM: 200 GeV freeze-out limit, freeze-in scope"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes the new $Z_{BL}$ boson does not mix with the photon or the $Z$ boson of the Standard Model; if that mixing is sizable, the relic-density, direct-detection, and collider constraints all change.","fun_headline_variants_meta":{"raw":{"variants":["B-L DM floor at 200 GeV; freeze-in widens probe","Freeze-out sets B-L DM floor; freeze-in maps lighter zone","B-L DM: 200 GeV floor, freeze-in window","B-L dark matter: freeze-out floor, freeze-in frontier","B-L DM: 200 GeV freeze-out limit, freeze-in scope"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001269,"raw_usage":{"total_tokens":5441,"prompt_tokens":1440,"completion_tokens":4001,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1056,"completion_tokens_details":{"reasoning_tokens":3909}},"tokens_in":1056,"tokens_out":4001,"duration_ms":29700,"temperature":1.0,"reasoning_tokens":3909,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:18:59.287705+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future direct detection experiment could falsify the thermal window by observing a spin-independent $\\zeta$-nucleon cross section larger than the value implied by $g_{BL}g_\\zeta$ from the allowed parameter curves at a claimed $m_\\zeta$; equivalently, a displaced-vertex search like FASER finding a $Z_{BL}$ with parameters in the region Fig. 8 excludes, or a confirmed thermal relic with $m_\\zeta<200$ GeV in this exact model, would disprove the paper's central bounds.","supporting_citations":[{"cited_title":"Cirelli, P","cited_arxiv_id":null,"evidence_quote":"Sets the observed dark matter abundance $\\Omega_{DM}h^2 = 0.12$ that all relic constraints are matched to."},{"cited_title":"Bernal, M","cited_arxiv_id":null,"evidence_quote":"Supplies the CMB and AMS-02 indirect detection bounds that exclude the low-mass thermal region."},{"cited_title":"Biswas and A","cited_arxiv_id":null,"evidence_quote":"Provides the spin-independent DM-nucleon cross section formula quoted in Eq. (14)."},{"cited_title":"Aprile et al","cited_arxiv_id":null,"evidence_quote":"Identifies sequential freeze-in, which the paper applies in Case B1."},{"cited_title":"Agnes et al","cited_arxiv_id":null,"evidence_quote":"Defines the FASER geometry and sensitivity used for the lifetime-frontier projections in Fig. 8."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides CMB-based lower bounds on dark matter mass and dark-photon portal constraints, used to take $m_\\zeta \\geq 1$ GeV."}],"review_version":1}