{"id":"0fa99760-65bc-4cc6-850d-ea5165afe80d","arxiv_id":"2412.02027","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"An anomaly-free two-component secluded WIMP model with a dark photon, dark Higgs, and scotogenic Dirac neutrino masses can account for the observed relic density while evading cosmological bounds.","lead":"This paper constructs a complete dark matter model containing two dark matter particles plus a dark photon and dark Higgs, and shows a region of its parameters matches the observed dark matter abundance while also giving neutrinos tiny masses. It matters because the model is testable by next-generation direct detection experiments such as DARWIN, and does not require kinetic mixing between the dark and visible sectors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Early-universe thermal equilibrium between the dark sector and SM is asserted in Sec. IV, not derived; Sec. V cites decoupling conditions but never shows the same tiny portals (epsilon <= 1e-6, theta < 1e-3, m_Xi > 1 TeV) both equilibrate at high T and decouple before BBN.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing concern: the early-universe thermal equilibrium between the dark sector and the Standard Model is assumed rather than derived. My stress-test agrees. The central claim that the model explains the relic abundance even without kinetic mixing rests on this assumption. If the tiny portals never thermalized the dark sector with the SM at high temperature, then the freeze-out equations (18) are not the correct description and the reported abundance and BBN compatibility are not established. This is not a demonstrated internal inconsistency; it is a missing calculation, so the appropriate verdict remains conditional rather than accepted or rejected. The proposed concrete test, a numerical integration of the coupled temperature and abundance evolution, would settle whether the assumption holds. The paper has some independent support: the one-component limit in Fig. 5 is checked against micrOMEGAs, and Appendix C provides an analytic annihilation cross-section, but neither addresses the thermalization history. Thus the reader's conditional verdict is unchanged.","tokens_in":21786,"tokens_out":4691,"duration_ms":52037,"concrete_test":"Numerically integrate the coupled Boltzmann equations for the dark-sector temperature and number densities (DM species, Z', h2, and right-handed neutrinos) using the SARAH/SPheno-generated couplings, starting from a high reheating temperature, for representative benchmark points with epsilon = 0, theta in [1e-6, 1e-3], and m_Xi > 1 TeV. Verify whether the dark and visible sectors reach kinetic and chemical equilibrium before DM freeze-out, and whether Z' and h2 have decayed before BBN. If no equilibrium is reached for allowed benchmarks, Section IV's 'Since the dark sector is in thermal equilibrium' fails and the relic-density calculation needs to be redone with the actual initial abundance.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim, including the 'even without kinetic mixing' compatibility statement, depends on the assumption that the dark sector and the right-handed neutrinos were in thermal equilibrium with the Standard Model in the early universe. Section IV opens with 'Since the dark sector is in thermal equilibrium,' and Section V asserts that the visible and dark sectors are in equilibrium unless the inter-sector interaction rate is much smaller than Hubble, with the small-portal regime (epsilon ~ 1e-6, theta < 1e-3, m_Xi > 1 TeV) cited as satisfying both equilibrium and decoupling, with reference to Ref. [24] but without showing a thermal-history calculation. For epsilon = 0, the only communications with the SM are the tiny scalar mixing theta and the scotogenic Yukawa sector in Eq. (6), so it is not self-evident that the dark sector ever reached the same temperature as the SM plasma, nor that the right-handed neutrinos equilibrated before decoupling. If the sectors never equilibrated, the freeze-out equations (18) use the wrong initial condition, the relic-density curves in Figs. 6 and 7 lose their basis, and the BBN and Delta Neff constraints in Appendix F, which assume a thermalized right-handed neutrino population, are not the relevant ones. The paper does not derive or numerically integrate the coupled dark-visible sector temperature evolution.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a UV-complete, anomaly-free extension of the Standard Model with a secluded U(1)_D gauge symmetry. The new fields (Table I) comprise a dark Higgs S that breaks U(1)_D, a dark photon Z', a dark Higgs h2, nine chiral fermions including three massless right-handed neutrinos charged under U(1)_D, and a decoupled inert scalar sector (eta, Phi). Two Dirac fermions, Psi1 and the lighter of the two mixed states Psi_2^1, are stable dark matter candidates; the right-handed neutrinos plus the inert scalars generate Dirac neutrino masses at one loop (Fig. 2), with the couplings y_nR solved from the measured neutrino masses and PMNS matrix (App. A). The relic density is computed in a one-component WIMP limit (Eqs. (24)-(25), App. C) and with micrOMEGAs 6.0.3 for the full two-component system (Eqs. (18), Figs. 6-7) over the scan ranges of Table III, selecting points with Omega h^2 equal to about 0.12. The paper shows direct-detection prospects (Figs. 9-10) and claims compatibility with invisible-decay, BBN, and Delta N_eff constraints, including for zero kinetic mixing (epsilon = 0).","tokens_in":22271,"tokens_out":29719,"duration_ms":262785,"significance":"The paper's main asset is completeness: it is the first phenomenological exploration of the model class catalogued in Ref. [23], with an explicit anomaly-free charge assignment, one-loop Dirac neutrino masses, and relic densities computed with standard tools. The analytic one-component calculation is cross-checked against micrOMEGAs (Fig. 5), and the annihilation cross-sections in App. C reproduce known vector and axial-vector limits. The fits are transparently labelled: y_nR is fixed 'by construction' (App. A) to reproduce neutrino data, and the scan is filtered on Omega h^2 = 0.12, so the results establish existence of viable parameter regions and falsifiable direct-detection targets (Fig. 9) rather than sharp predictions. The significance is tempered by the two gaps identified below: the assumed early-universe thermal history of the dark sector, and the epsilon = 0 Delta N_eff claim; both are load-bearing for the headline 'even without kinetic mixing' statement.","major_comments":[{"comment":"The stress-test concern about the thermal history is warranted. Section IV opens with 'Since the dark sector is in thermal equilibrium,' and the same section asserts that the visible and dark sectors decouple when epsilon ~ 10^-6, theta < 10^-3, and m_Xi > 1 TeV, citing Ref. [24]; neither statement is derived for the model's own couplings. This assumption is load-bearing because the Boltzmann equations (18), the analytic relic density (24), and the Delta N_eff treatment of Appendix F all presuppose a single thermal bath at temperature T with standard equilibrium initial abundances. At epsilon = 0 the only channels connecting the dark sector to the SM are the scalar mixing theta and the y_nL Yukawa couplings to inert scalars with m_Xi > 1 TeV, and it is a parametric question (depending on the Table III ranges) whether these equilibrate the dark sector at T larger than about m_Psi/25 and later decouple it before BBN. If the dark sector never equilibrated, the correct computation would be freeze-in or a two-temperature Boltzmann system, and Figs. 5-7 would not determine the relic abundance. I request a rate-versus-Hubble comparison for the dominant DS-SM equilibration processes for the benchmark of Fig. 5 and for representative epsilon = 0 points, or a coupled T_DS/T_SM evolution.","section":"Sec. IV; App. F"},{"comment":"The epsilon = 0 Delta N_eff claim is not supported by the presented calculation. Appendix F computes Gamma_nuR only for the kinetic-mixing-mediated coupling (Eqs. (F2)-(F3)), and the text concludes that the blue epsilon = 0 region in Fig. 9 has 'zero contribution to the relativistic degrees of freedom at tree level.' But at epsilon = 0 the dark sector is still populated (the freeze-out mechanism requires it), and after freeze-out the DM annihilation energy is transferred to Z' and h2, which decay predominantly into the massless right-handed neutrinos; that nu_R population is a decoupled dark-radiation component whose contribution to Delta N_eff is not computed anywhere in the paper. The magenta-region constraint in Fig. 9 therefore reflects only the kinetic-mixing channel, and the abstract and Sec. VI claims of compatibility with Delta N_eff 'even without kinetic mixing' need either a computation of this channel or a quantitative demonstration that it is negligible (e.g., early decoupling and entropy dilution).","section":"App. F; Fig. 9"},{"comment":"The paper's headline statements are existence claims, and the reader should be able to judge how representative the viable regions are. The scan filters on Omega h^2 = 0.12 (Eq. (21)) and fixes y_nR from neutrino data (Eqs. (A13)-(A14)); the manuscript should report the fraction of randomly scanned points that pass Eq. (21) and all other constraints, and how the viable regions shift under the Table III bounds. Note also a quantitative inconsistency: Eq. (21) imposes Omega h^2 = 0.1200 +- 0.0012, while Fig. 6 and its text use a 10% band; the actual filter used to produce Figs. 6-10 should be stated once, precisely.","section":"Eq. (21); Table III; Figs. 6-10"}],"minor_comments":[{"comment":"Equation (F1) is dimensionally inconsistent as printed: Gamma_nuR has dimension of energy while 3H n_nuR has dimension of energy to the fourth power; the intended condition is presumably Gamma_nuR > 3H, so the 'n_nuR' should be removed.","section":"Eq. (F1)"},{"comment":"The notation for the second dark matter fermion is hard to parse: the mass eigenstates are written as m_Psi1_2 and m_Psi2_2 in the Fig. 5 caption and as Psi_2^1, Psi_2^2 elsewhere; a consistent notation such as m_{Psi_2^{(1)}} and m_{Psi_2^{(2)}} would remove the ambiguity with the first DM candidate Psi1.","section":"Figs. 5-7"},{"comment":"The text says the relic abundance is 'almost always dominated by annihilation processes directly to the SM particles,' but in Table II the classes 1100 and 2200 include final states with Z' and h2, and the secluded regime is dominated by annihilation into dark mediators; the phrase should be reconciled with the classification.","section":"Sec. V (Fig. 7 discussion)"},{"comment":"Please state whether the SARAH/SPheno/micrOMEGAs model files and the scan-filtering code will be made available; the reproducibility of Figs. 6-10 depends on these.","section":"Sec. V"},{"comment":"Minor typos and infelicities: 'Bolztmann equations' (Sec. V), 'the branching ratio of Standard model particles into invisible' (abstract), and 'that is below the current limit of XENONnT for the two DM components' (Sec. V) should be reworded.","section":"Abstract; Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of a phenomenological journal, and its differentiation from the overlapping Ref. [66] (chiral dark matter with scotogenic Majorana neutrino masses) is clearly stated and genuine: here the right-handed neutrinos are massless and the scotogenic masses are Dirac. My main request is a thermal-history section; without it, the two headline claims (relic abundance and Delta N_eff compatibility at epsilon = 0) rest on an unverified initial condition. I would also ask the editor to encourage the authors to include scan statistics (number of accepted points and passing fractions) and to report the epsilon = 0 benchmark in detail. The dimensional typo in Eq. (F1) suggests that the Delta N_eff numerics in Fig. 9 deserve a careful recheck."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First numerical phenomenology for a specific anomaly-free chiral U(1)D secluded WIMP model with two-component DM and scotogenic Dirac neutrino masses. It does what it sets out to do: the internal calculations are standard, the one-component limit matches micrOMEGAs, and the 2211 conversion process is highlighted for the first time in this benchmark. The p-wave annihilation cross-section in Eq. (C4) is new as far as I know.\n\nCredit where due: the model is anomaly-free, the mass matrices and cross sections are explicit, and the Fig. 5 benchmark agrees with the numerical code. The direct-detection and Delta Neff plots give a useful map of the viable parameter space. The paper is also honest about the fitting nature: couplings ynR are solved from neutrino data in Appendix A, and the relic density is imposed by scanning. That is not hidden.\n\nSoft spots, in order. No public code or benchmark files, so the scan cannot be independently reproduced. The abstract's \"explains the relic abundance\" overstates things; the model fits Omega h^2 by construction. That is normal for the field but it is not a prediction. The thermal-equilibrium assumption is the one I'd push on. Section IV opens with \"Since the dark sector is in thermal equilibrium\" and Section V asserts decoupling for tiny epsilon and theta without showing the rate calculation or thermal history. For epsilon = 0, the only links to the SM are the scalar mixing and the scotogenic Yukawa couplings, and it is not self-evident that those keep the dark sector at SM temperature down to freeze-out while still passing Delta Neff and BBN. The citation to Ref. [24] may cover it, but the paper should reproduce the key estimate or phrase the assumption more carefully. That is a legitimate referee question, not a fatal flaw: the portals are there, so a successful thermal region plausibly exists.\n\nThe paper is for model-builders working on secluded DM, Dirac neutrino masses, or Abelian dark extensions. It deserves a serious referee; the major comments should be about thermalization, benchmark parameters, and releasing the code.","headline":"A solid first phenomenology scan of a specific chiral secluded WIMP benchmark; the thermal-equilibrium assumption is the main thing to push on.","tokens_in":22780,"tokens_out":5215,"would_cite":true,"duration_ms":53456,"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":"The paper builds a complete dark-sector model in which two WIMP dark matter candidates, a dark photon, and a dark Higgs arise from one extra Abelian gauge symmetry, and claims this sector alone can produce the observed relic abundance…","keywords":["secluded WIMP dark matter","dark photon","dark Higgs","Dirac neutrino mass","scotogenic mechanism","two-component dark matter","U(1)D gauge symmetry","anomaly-free chiral fermions"],"falsifier":"Measure the dark photon's invisible decay to neutrinos at a beam-dump or fixed-target experiment over the 1-100 GeV mass range; if a $Z'$ with the required coupling is absent, or its lifetime exceeds one second so it survives to Big Bang nucleosynthesis, the model's central compatibility claim is ruled out.","tokens_in":21582,"feed_emoji":"🌌","tokens_out":11024,"duration_ms":106838,"temperature":0.7,"pith_summary":"The paper constructs a complete, renormalizable extension of the Standard Model in which dark matter is not one particle but two weakly interacting massive fermions that annihilate inside an isolated 'secluded' dark sector. The central claim is that one extra Abelian gauge symmetry, broken by a dark Higgs, leaves a dark photon and a dark Higgs as the only portals, and that this sector can reproduce the observed relic abundance $\\Omega h^{2}\\approx 0.12$ even when the kinetic mixing between the dark photon and the Standard Model is set to zero. The same sector, through two inert scalars that never acquire vacuum values, generates small Dirac masses for the three Standard Model neutrinos at one loop, matching current oscillation data by construction. A sympathetic reader would care because it shows a UV-complete example of a class of models—chiral, anomaly-free, multi-component dark sectors with radiative neutrino masses—that survives cosmological and laboratory constraints and leaves some parameter space within reach of next-generation direct-detection experiments.","feed_headline":"Two dark particles explain relic dark matter and neutrino masses","feed_subtitle":"A secluded WIMP sector with a dark photon and dark Higgs reproduces the measured relic density, with zero kinetic mixing still allowed.","key_machinery":"The central object is the anomaly-free charge vector $D=[(1,-10),(1,-10),(-4,-5),9,9,9]$ under $U(1)_D$, with dark Higgs charge $S=9$. This vector solves the anomaly equations, gives three massive Dirac pairs after symmetry breaking (one of which is a two-flavor system whose lightest mass eigenstate is the first dark matter candidate, and a one-generation Dirac fermion is the second), and leaves three massless charge-9 right-handed neutrinos. A dark photon $Z'$ and a dark Higgs $h_2$ mediate annihilation; two inert scalars, a doublet $\\eta$ and a singlet $\\Phi$, close the one-loop 'scotogenic' Dirac neutrino mass diagram, meaning neutrinos remain massless at tree level. The load-bearing equation is the coupled Boltzmann system, whose conversion term $\\langle\\sigma_{2211}v\\rangle$ couples the two number densities and is what makes the two-component nature change the relic density.","core_discovery":"The central claim is that a single self-consistent charge assignment $D=[(1,-10),(1,-10),(-4,-5),9,9,9]$ under $U(1)_D$, with a dark Higgs of charge $9$, yields all the ingredients at once: two stable Dirac fermions (the lightest state of a two-flavor pair and a single one-generation Dirac fermion), a dark photon that decays directly into three right-handed neutrinos, and a radiative Dirac neutrino mass from a decoupled inert scalar sector. Solving the two coupled Boltzmann equations with the conversion process $\\Psi_2\\Psi_2\\to\\Psi_1\\Psi_1$, the paper finds parameter regions where the two components together saturate $\\Omega h^{2}=0.12$, with dark matter masses from about 1 GeV to a few TeV. The neutrino decay mode of the dark photon opens a channel that keeps the model consistent with Big Bang nucleosynthesis and $\\Delta N_{\\text{eff}}$ even at zero kinetic mixing, and a scan under all constraints leaves viable points above the projected sensitivity of future liquid-xenon direct-detection experiments.","pith_inferences":["An extension the authors do not pursue: because the dark photon decays predominantly to neutrinos, its most direct experimental signature may be missing-energy events in beam-dump and fixed-target searches rather than dilepton resonances.","The charge assignment studied here is one benchmark from a family of about a thousand anomaly-free solutions, and the paper's own comparison suggests the same two-mediator, two-component phenomenology can be realized by other charge vectors; whether the relic-density behavior is identical for the rest of the family is a testable extension.","The near-decoupled dark Higgs, with mixing angle below $10^{-3}$ and mass near 550 GeV, could still be probed by precision measurements of the 125 GeV Higgs couplings at future colliders, giving a complementary handle on the viable points that direct detection cannot reach.","Quantifying the cascade annihilation $\\Psi\\Psi\\to Z'Z'\\to 4\\nu$ at neutrino telescopes could make the model's otherwise dark indirect signal visible; the paper notes the neutrino-flux enhancement but leaves the calculation for future work."],"forward_implications":["With zero kinetic mixing, the relic abundance can still reach $\\Omega h^{2}=0.12$, because the dark photon's decay into right-handed neutrinos provides the connection back to the Standard Model without needing $\\epsilon$.","Whenever the two dark matter candidates are close in mass and the dark Yukawa coupling $y_c$ is above about 0.1, the conversion process $\\Psi_2\\Psi_2\\to\\Psi_1\\Psi_1$ affects the heavier component's abundance, so isolated single-component freeze-out calculations are not reliable in this model.","Spin-independent direct detection is dominated by the vector portal and scales as $\\epsilon^{2}$, so the $\\epsilon=0$ region is currently invisible; nonetheless a slice of the viable parameter space lies above the projected sensitivity of the next generation of xenon-based detectors.","Annihilation into the Standard Model is dominated by neutrino final states, making the model effectively dark for gamma-ray telescopes; the strongest experimental probes are therefore direct detection, collider searches, and possibly neutrino telescopes.","The neutrino spectrum is Dirac and radiatively generated, so the model predicts no neutrinoless double beta decay, in contrast with Majorana seesaw constructions."],"supporting_citations":[{"why":"Defines the secluded-WIMP dark matter paradigm, where the relic abundance is set by annihilation into lighter dark mediators that later decay to the Standard Model; the whole dark-sector freeze-out analysis rests on this mechanism.","marker":"[14]"},{"why":"Establishes the linear and cubic anomaly cancellation conditions that fix which chiral fermion charge sets are allowed under the extra Abelian symmetry.","marker":"[22]"},{"why":"Supplies the effective Dirac neutrino mass operator and the enumeration of roughly one thousand anomaly-free charge solutions from which the benchmark charge set is taken.","marker":"[23]"},{"why":"Provides the inert scalar sector that realizes one-loop scotogenic Dirac neutrino masses and the decoupling condition (heavy inert scalars) used in the scan.","marker":"[24]"},{"why":"Gives the Dirac scotogenic construction with anomaly-free Abelian symmetries that links the right-handed neutrinos, the two-flavor dark fermions, and the inert scalars in the neutrino-mass loop.","marker":"[25]"},{"why":"Numerically solves the coupled two-component Boltzmann equations including the dark matter conversion term, which is how the relic density values are actually computed.","marker":"[40]"},{"why":"Provides the observed relic abundance $\\Omega h^{2}\\approx 0.12$ that the parameter scan is required to reproduce.","marker":"[41]"},{"why":"Sets the current spin-independent direct-detection bound used to filter the scanned models.","marker":"[71]"}],"fun_headline_variants":["Secluded WIMP duo sets relic density and neutrino mass","Two-component dark matter: relic density and Dirac neutrinos","Dark photon and Higgs: relic density and neutrino mass","Zero kinetic mixing with two WIMP dark matter components","Secluded WIMPs: two components fill relic density and neutrino mass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The dark sector and the Standard Model plasma came into thermal equilibrium in the early Universe, and the same tiny couplings later decoupled it before Big Bang nucleosynthesis; the paper takes this equilibrium as a starting point rather than deriving it from a full thermal-history calculation.","fun_headline_variants_meta":{"raw":{"variants":["Secluded WIMP duo sets relic density and neutrino mass","Two-component dark matter: relic density and Dirac neutrinos","Dark photon and Higgs: relic density and neutrino mass","Zero kinetic mixing with two WIMP dark matter components","Secluded WIMPs: two components fill relic density and neutrino mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00052,"raw_usage":{"total_tokens":2495,"prompt_tokens":902,"completion_tokens":1593,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":518,"completion_tokens_details":{"reasoning_tokens":1511}},"tokens_in":518,"tokens_out":1593,"duration_ms":10694,"temperature":1.0,"reasoning_tokens":1511,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:55:39.439281+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the dark photon's invisible decay to neutrinos at a beam-dump or fixed-target experiment over the 1-100 GeV mass range; if a $Z'$ with the required coupling is absent, or its lifetime exceeds one second so it survives to Big Bang nucleosynthesis, the model's central compatibility claim is ruled out.","supporting_citations":[{"cited_title":"Linkage of Dirac Neutrinos to Dark U(1) Gauge Symmetry","cited_arxiv_id":"2101.12138","evidence_quote":"Numerically solves the coupled two-component Boltzmann equations including the dark matter conversion term, which is how the relic density values are actually computed."}],"review_version":1}