{"id":"228af7f2-5da5-48b0-85ff-c2fdd2f88093","arxiv_id":"2411.15108","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"A feebly interacting dark matter particle can become a pseudo-FIMP by staying in equilibrium through conversion with a strongly self-interacting SIMP partner; the paper demonstrates this with coupled Boltzmann equations and a two-scalar example.","lead":"Dark matter could be two particles at once: one with very weak ordinary interactions and one with strong self-interactions, and the pair can stay in thermal balance through a conversion process. This paper studies that two-particle setup, called pseudo-FIMP plus SIMP, and builds the simplest two-scalar model that produces it.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Modified-equilibrium formulas (3.3)-(3.4) are asserted without derivation and are dimensionally inconsistent, so the pFIMP freeze-out behavior at large conversion is not established.","rationale":"The reader's weakest assumption concerns kinetic equilibrium at the scan value lambda_chiH ~ 1e-3, which is a legitimate concern. However, the more load-bearing issue is the modified equilibrium formulas (3.3)-(3.4): they are presented as the basis for interpreting the numerical solutions as pFIMP freeze-out, but they are (a) not derived, and (b) dimensionally inconsistent as written. Even if the units are a typographical artifact, the lack of a derivation means the central mechanism---that the feeble component reaches and tracks an equilibrium curve before freeze-out---is asserted rather than demonstrated in the model-independent Sec. 3. This is distinct from and more central than the kinetic-equilibrium check, because it directly underwrites the pFIMP identification and the modified freeze-out behavior that produces the headline mass reach. I therefore agree partially with the reader: the kinetic-equilibrium mismatch is real, but the equilibrium-yield derivation deserves priority. The paper otherwise has independent support in the form of explicit cross-section expressions and numerical comparisons, and the existence of some pFIMP-like behavior for intermediate lambda_chiPH is plausible; I would not move the verdict to REJECT because the issue can likely be repaired with a derivation or corrected formulas. The verdict should remain CONDITIONAL, with the requested revision now prioritized on eqs. (3.3)-(3.4) and the numerical code/scan reproducibility.","tokens_in":21990,"tokens_out":1756,"duration_ms":15101,"concrete_test":"Re-derive the modified equilibrium conditions by setting dY_w/dx = 0 and dY_s/dx = 0 in eqs. (3.1)-(3.2) simultaneously, and check whether the resulting algebraic relations reduce to eqs. (3.3) and (3.4) under the hierarchy gamma_sw >> gamma_ww, with all terms written with consistent units (e.g., multiplying eq. (3.4) by the appropriate power of s). If the algebra yields different expressions, recompute the relic-density points in fig. 4 using the corrected equilibrium yields; if the pFIMP region III and the headline mass reach (~50 MeV) shift materially, the central claim needs revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that a feeble DM component w reaches and follows thermal equilibrium via SS->ww conversion, with \"modified equilibrium\" yields given by eqs. (3.3) and (3.4). Yet these formulas are stated without derivation in Sec. 3, and they do not pass dimensional inspection. In eq. (3.3), n_eq'_w = n_eq_w (n_s / n_eq_s) has units of (number density)^2, not number density unless the ratio is interpreted as dimensionless, which the notation does not support. In eq. (3.4), the factor (n_w / n_eq_w)^2 is dimensionless, but inside the bracket it multiplies <sigma v>_{ss->ww} (units cm^3 s^-1) while the other two terms have units cm^3 s^-1; after the square root the bracket has units (cm^3 s^-1)^{1/2}, which cannot be added to <sigma v>_{ss->SM SM} (cm^3 s^-1) inside the numerator. The cBEQ source terms in eqs. (3.1)-(3.2) are rate equations that do not by themselves imply the algebraic modified-equilibrium conditions. Because the pFIMP identification (region III) and the statement that the FIMP \"follows equilibrium before freeze-out\" rest on these formulas, the quantitative freeze-out dynamics in the pFIMP regime is not demonstrated. The numerical solutions in figs. 2-3 may be correct, but the analytic interpretation that the particle follows the modified equilibrium curve is unsupported without a derivation that clarifies the units and approximations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a two-component dark matter setup in which one component is a SIMP and the other has only feeble couplings to the visible sector but a sizeable coupling to the SIMP, making the latter a pseudo-FIMP (pFIMP). After reviewing the single-component SIMP, the authors present a model-independent analysis based on solving the coupled Boltzmann equations (3.1)-(3.2), with the DM-DM conversion rate varied from negligible to large. They identify four regimes, including a pFIMP regime in which the feeble component reaches thermal equilibrium through conversion and freezes out. They then construct a concrete two-scalar model with Z2 x Z3 symmetry, solve the coupled equations, scan the parameter space under relic density, unitarity, and self-interaction constraints, and discuss detection prospects through a vector-like lepton extension. The paper claims that the SIMP mass range is extended up to about 50 MeV when the pFIMP is heavier.","tokens_in":22422,"tokens_out":9053,"duration_ms":91517,"significance":"If correct, the paper extends the pFIMP mechanism from a WIMP partner to a SIMP partner and provides the simplest scalar realization. The numerical solution of coupled Boltzmann equations is a credible method, and the paper includes useful appendices: a semi-analytic SIMP solution compared with numerical results (Appendix A), cross-section formulas (Appendix B), self-interaction expressions (Appendix C), and a kinetic-equilibration estimate (Appendix D). The concrete model and the parameter-space scan, with self-interaction and unitarity constraints, give falsifiable predictions for SIMP masses and couplings. However, the central pFIMP identification rests on the modified-equilibrium formulas (3.3)-(3.4), which are presented without derivation, and the kinetic-equilibrium assumption is verified only for parameters different from those used in the main scan. These issues are load-bearing and require a major revision.","major_comments":[{"comment":"The pFIMP regime is characterized by the statement that the feeble component 'follows equilibrium before freeze out,' but the modified equilibrium number densities in Eqs. (3.3) and (3.4) are asserted without derivation. Please derive these expressions from the coupled Boltzmann equations in the limit of large conversion rate, stating all approximations (e.g., neglect of SM production/destruction terms, steady-state condition dY/dx ~ 0), and validate them by comparing with the numerical solutions in Fig. 2. Note that, contrary to a purely dimensional objection, the formulas are dimensionally consistent if n_s and n_w denote number densities: in Eq. (3.4) the numerator second term (n_s^2/n_eq_s)<sigma v2> has units cm^3 s^-1, matching the other terms. However, the asymmetric appearance of n_s^2/n_eq_s in the numerator and n_s in the denominator needs clarification, as does the meaning of a 'modified equilibrium' for a species that is not itself in chemical equilibrium. Without this derivation, the quantitative identification of region III and the claim that the pFIMP tracks the modified equilibrium are not established.","section":"Sec. 3, Eqs. (3.3)-(3.4)"},{"comment":"The coupled Boltzmann equations assume that the dark and visible sectors share a single temperature. Appendix D demonstrates that the SIMP kinetic-equilibration condition Gamma_{chi f -> chi f} > H is satisfied for lambda_chiH ~ 0.1, whereas the relic-density scan in Sec. 4.2 fixes lambda_chiH ~ 1e-3. Since the elastic scattering rate scales approximately as lambda_chiH^2, a reduction by two orders of magnitude may invalidate the kinetic-equilibrium assumption at the benchmark points used in the scan. Please compute Gamma_{chi f -> chi f}/H at the freeze-out temperature for representative scan points with lambda_chiH = 1e-3, or impose the kinetic-equilibrium condition in the scan. This is load-bearing because a dark-sector temperature different from the SM bath temperature changes the form of the Boltzmann equations and the freeze-out conclusions.","section":"Sec. 4.2 and Appendix D"},{"comment":"The model-independent analysis relies on hand-picked numerical values for <sigma v>_{ss->SM SM}, <sigma v>_{ww->SM SM}, <sigma v2>_{3s->2s}, and the conversion cross-section, with no exploration of how the four-region classification and the pFIMP threshold depend on these inputs. Since the paper claims a model-independent conclusion, please show the robustness of the pFIMP regime under order-of-magnitude variations of these cross-sections, or clearly state that the conclusions are illustrative. At minimum, specify how the conversion cross-section maps to the rate ratios gamma_sw/gamma_3s->2s used to define regions III and IV.","section":"Sec. 3, Fig. 2"},{"comment":"In the concrete model, the conversion term has a factor 1/4 in the pFIMP equation (4.2) and a factor 1/2 in the SIMP equation (4.3). This is consistent with the definition Y_s = 2Y_chi and the process chi chi* -> phi phi, but the reasoning is not stated. Please spell out the connection between Y_s, Y_chi, and the symmetry factors so that the density-balance between the two equations is transparent to the reader.","section":"Sec. 4.1, Eqs. (4.2)-(4.3)"}],"minor_comments":[{"comment":"The symbol mu_sw in the definitions of H(x), s, and Y_eq is not defined. Since the horizontal axis in Fig. 2 is labeled mu_sw/T, please define mu_sw explicitly and clarify the convention for x in a two-mass system.","section":"Sec. 3, after Eq. (3.2)"},{"comment":"In several places the notation 'sigma^T_{ss->ww}' or 'sigmaT' appears without definition; use a consistent notation for thermally averaged cross-sections.","section":"Sec. 3"},{"comment":"The text 'pFIMP-SMIP model' in the conclusions is a typo; it should read 'pFIMP-SIMP model.'","section":"Sec. 5"},{"comment":"The sentence 'allows lambda_chi ~ x 10^-2' is incomplete; it should read 'lambda_chi ~ 10^-2' or give the explicit numerical value.","section":"Sec. 4.2"},{"comment":"The phrase 'The ncecessary condition' is a typo; it should be 'The necessary condition.'","section":"Sec. 2"},{"comment":"The text says 'SIMP mass is allowed up to ~ 50 MeV when m_phi > m_chi' in one place and 'Delta m <= 100 MeV' in the conclusions; please ensure the mass-separation statements are consistent and clearly defined.","section":"Sec. 4.2 and Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a development of the authors' earlier pFIMP-WIMP work; the SIMP partner is a new application. The dimensional consistency objection raised in one reading does not, on inspection, invalidate Eqs. (3.3)-(3.4) if number densities are meant, but the missing derivation and the kinetic-equilibrium parameter mismatch are genuine issues that should be fixed before publication. The paper fits the journal's scope and, after the requested revision, could be suitable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper does what it says—demonstrates that a feebly coupled DM component can be brought into equilibrium by a SIMP partner, both in a model-independent setup and in a concrete two-scalar Z2 x Z3 model. The existence proof holds up. The main weaknesses are a missing derivation of the modified-equilibrium formulas (3.3)-(3.4) and a kinetic-equilibrium assumption checked at a different coupling than the scan uses. Both are fixable.\n\nThe novelty is real. pFIMP was previously studied with a WIMP partner; this is the first application to SIMP, and the coupled Boltzmann treatment plus the Z2 x Z3 benchmark give concrete results, including the 50 MeV SIMP mass reach when the pFIMP is heavier. The paper is honest about using fixed numbers in the model-independent section, and the self-interaction constraints are properly applied.\n\nI disagree with the stress-test note I was handed: it claims eqs. (3.3) and (3.4) are dimensionally inconsistent. That is not right. In (3.3), n_s/neq_s is dimensionless, so neq'_w has number-density units. In (3.4), the terms inside the numerator and denominator each have cm^3/s units—the (n_w/neq_w)^2 factor is dimensionless, and n_s <sigma v^2> has cm^3/s because <sigma v^2> carries cm^6/s for a 3->2 process. So the ratio is dimensionless and the square root is well-defined. The formulas may be questionable, but on dimensional grounds they are fine.\n\nThe real problem is that (3.3)-(3.4) are asserted without derivation. They are the basis for identifying the pFIMP regime: the claim that the FIMP follows a modified equilibrium curve before freeze-out rests on these equations. A referee should ask the authors to derive them from the cBEQs in the limit where conversion dominates, or at least to state the approximations clearly. Without that, the analytic interpretation is under-supported, though the numerical curves in figs. 2-3 stand on their own.\n\nThe other substantive gap is kinetic equilibrium. Appendix D shows a dark-SM common temperature holds for lambda_chiH ~ 0.1, but the relic scan fixes lambda_chiH ~ 1e-3. The paper acknowledges this in Sec. 5, but does not verify it at the scan value. That matters because the freeze-out conclusions assume T_dark = T_SM, and the fix should be easy.\n\nWho is this for? Phenomenologists working on multi-component DM, especially SIMP models. It is a useful existence proof and a clean benchmark. It deserves a serious referee; it is not ready for acceptance as-is, but the fixes are straightforward. I would send it to review with a request for the derivation and the kinetic-equilibrium check.","headline":"A legitimate, mostly sound extension of pFIMP to SIMP partners: the existence proof holds, but the modified-equilibrium formulas need derivation and the kinetic-equilibrium check is done at a different coupling than the scan uses.","tokens_in":22988,"tokens_out":4317,"would_cite":true,"duration_ms":36655,"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 shows that a feebly interacting dark-matter component can become a pseudo-FIMP and freeze out like a thermal relic when its dark-sector partner is a strongly interacting SIMP.","keywords":["pseudo-FIMP","SIMP dark matter","two-component dark matter","coupled Boltzmann equations","relic density","dark-matter self-interaction","Z2 x Z3 scalar dark matter","freeze-out"],"falsifier":"Compute $\\sum_f \\Gamma_{\\chi f\\to\\chi f}/H(T)$ at $T=m_\\chi/25$ for a relic-allowed benchmark with $\\lambda_{\\chi H}=10^{-3}$ and $m_\\chi$ in the $10$\\textendash$50$ MeV range; if this ratio drops below about one, dark and visible temperatures decouple and the Boltzmann solutions that assume $T_{\\rm dark}=T_{\\rm SM}$ miscompute freeze-out.","tokens_in":21793,"feed_emoji":"🌌","tokens_out":8644,"duration_ms":78299,"temperature":0.7,"pith_summary":"The paper argues that a dark-matter component with only feeble couplings to ordinary matter can nevertheless be a thermal relic if it lives alongside a strongly interacting massive particle (SIMP): once the conversion process that turns two SIMPs into two of the feeble particles is fast enough, the feeble component reaches equilibrium and freezes out like a WIMP, becoming what the paper calls a pseudo-FIMP. The claim is established first model-independently, by solving the coupled Boltzmann equations (3.1)–(3.2) and scanning the conversion rate, and then in a concrete model with a real scalar pFIMP and a complex scalar SIMP stabilised by a $\\mathbb{Z}_2\\otimes\\mathbb{Z}_3$ symmetry. In that model the relic-density-allowed region includes SIMP masses up to about $50$ MeV when the pFIMP is the heavier component, and the dominant constraint is dark-matter self-interaction. If correct, this widens the class of viable dark-matter production mechanisms beyond single-component SIMP or WIMP–FIMP setups, and it ties the detectability of the feeble component to the dark-sector partner's access to the visible sector.","feed_headline":"SIMP partner can make a feeble dark-matter component a thermal relic","feed_subtitle":"Coupled Boltzmann equations show the feeble component thermalizes through dark-sector conversion, expanding the viable relic parameter…","key_machinery":"The engine of the argument is the pair of coupled Boltzmann equations (3.1)–(3.2) tracking the yields of the weak component $w$ and the SIMP $s$, with the cross-component conversion term $\\langle\\sigma v\\rangle_{ss\\to ww}$ mediating energy and number exchange. A SIMP is defined by a number-changing $3\\to2$ self-annihilation in the dark sector, so the SIMP equation also carries the $\\langle\\sigma v^2\\rangle_{3s\\to2s}$ term; when the conversion rate is large enough, the feeble component follows a modified equilibrium distribution before freeze-out, with the form depending on whether the SIMP or the pFIMP is heavier. The concrete model supplies the same physics through a real scalar $\\phi$ (pFIMP) and a complex scalar $\\chi$ (SIMP) interacting through the portal coupling $\\lambda_{\\chi\\phi}$; varying this single coupling moves the system from pure FIMP behaviour through pFIMP freeze-out to the regime where the heavier component depletes into the lighter one.","core_discovery":"The central discovery is that SIMP dark matter can host a pseudo-FIMP. In the model-independent treatment, the yield $Y_w$ of the weakly coupled component and $Y_s$ of the SIMP obey coupled Boltzmann equations whose conversion term $\\langle\\sigma v\\rangle_{ss\\to ww}$ is the knob. For negligible conversion the weak component is an ordinary freeze-in FIMP; once the conversion rate $\\gamma_{sw}$ becomes comparable to the SIMP's self-annihilation rate $\\gamma_{3s\\to2s}$ it tracks equilibrium and freezes out with a density locked to the SIMP's, with the modified equilibrium yields given by eqs. (3.3) and (3.4) for the two mass hierarchies. The concrete $\\mathbb{Z}_2\\otimes\\mathbb{Z}_3$ scalar model realizes this with $\\phi$ as the pFIMP and $\\chi$ as the SIMP; after imposing relic density, unitarity, perturbativity, vacuum stability and self-interaction bounds, the surviving parameter space allows SIMP masses up to about $50$ MeV when the pFIMP is heavier, and it is the self-interaction bounds, not relic density alone, that most tightly fix the portal couplings $\\lambda_\\phi$ and $\\lambda_{\\chi H}$.","pith_inferences":["If the shared-temperature assumption fails at the scan value $\\lambda_{\\chi H}\\sim10^{-3}$, the freeze-out temperatures and relic abundances computed in Sec. 4.2 would need to be redone with two independent dark-sector and visible-sector temperatures; a dedicated kinetic-equilibration scan over the whole relic-allowed region would settle this.","The same coupled-equation structure should apply to pFIMP partners other than scalars, such as fermionic SIMPs or dark vector mesons, whenever a $3\\to2$ process sets the bath density.","In the large-conversion regime IV the 'SIMP' stops being defined by its own $3\\to2$ freeze-out, so one should expect its self-interaction phenomenology to be diluted; comparing halo-shape predictions between regions III and IV could serve as a model-independent test.","A measurement of the dark-matter momentum distribution, or of dark radiation, at MeV scales could distinguish a pFIMP from a freeze-in FIMP even when the total relic density is fixed, because the two production histories give different phase-space and temperature evolutions."],"forward_implications":["A feebly coupled dark-matter candidate does not have to be produced by freeze-in; in multicomponent models it can freeze out after equilibrating through dark-sector conversion with a SIMP partner.","In the two-component model the relic-allowed SIMP mass range reaches about $50$ MeV when the pFIMP is heavier, and the self-interaction bound, rather than relic density alone, sets the strongest limits on the parameter space.","Self-interaction constraints from Bullet and Abell clusters become the deciding phenomenological test, constraining couplings like $\\lambda_\\phi$ and $\\lambda_{\\chi H}$ that barely affect the relic abundance.","Direct and indirect detection of the pFIMP is hard unless the SIMP communicates with the visible sector through a light mediator; the vector-like-lepton extension discussed in the paper opens electron-scattering and annihilation channels.","The four dynamical regions identified by the conversion-rate ratio give a classification scheme: pure SIMP plus FIMP, converted FIMP, pFIMP freeze-out, and conversion-dominated depletion."],"supporting_citations":[{"why":"Defines the pseudo-FIMP concept for a thermal partner; this work transfers that construction to a SIMP partner.","marker":"[36]"},{"why":"Supplies the minimal complex-scalar $\\mathbb{Z}_3$ SIMP model, its Lagrangian, and the semi-analytic freeze-out solution that the two-component model builds on.","marker":"[19]"},{"why":"Introduces the SIMP mechanism in which $3\\to2$ dark-sector annihilation sets the relic density.","marker":"[14]"},{"why":"Provides a concrete SIMP model and the freeze-out dynamics that the paper adopts for $\\chi$.","marker":"[15]"},{"why":"Gives a multi-component SIMP model with $\\mathbb{Z}_2\\times\\mathbb{Z}_3$ structure and the effective self-interaction cross-section formula used for the two-component case.","marker":"[33]"},{"why":"Source of the dark-matter self-interaction framework and the Bullet and Abell cluster constraints that dominate the parameter-space scan.","marker":"[17]"},{"why":"Supplies the observed relic density $\\Omega_{\\rm DM}h^2=0.1200\\pm0.0012$ that the scan targets.","marker":"[6]"},{"why":"Underpins the kinetic-equilibration and self-interaction treatment for SIMP dark matter, including the DM-fermion elastic scattering rate.","marker":"[41]"},{"why":"Establishes that pFIMP detection proceeds through a thermal-DM loop, the strategy the paper extends to the SIMP case with a light mediator.","marker":"[37]"}],"fun_headline_variants":["SIMP can turn a feeble dark matter component into a thermal relic","pFIMP thermalizes via SIMP conversion in two-component dark matter","SIMP conversion yields a thermal pFIMP","SIMP hosts a pseudo-FIMP unlocking new relic density paths","Coupled SIMP-pFIMP dynamics expand viable dark matter parameter space"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the dark sector and the Standard Model bath share a single temperature throughout freeze-out; the paper's kinetic-equilibration check uses a Higgs-portal coupling of order $0.1$, while the relic scan fixes $\\lambda_{\\chi H}\\sim10^{-3}$, so the shared-temperature condition is not demonstrated at the value used in the scan.","fun_headline_variants_meta":{"raw":{"variants":["SIMP can turn a feeble dark matter component into a thermal relic","pFIMP thermalizes via SIMP conversion in two-component dark matter","SIMP conversion yields a thermal pFIMP","SIMP hosts a pseudo-FIMP unlocking new relic density paths","Coupled SIMP-pFIMP dynamics expand viable dark matter parameter space"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001137,"raw_usage":{"total_tokens":4750,"prompt_tokens":1001,"completion_tokens":3749,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":617,"completion_tokens_details":{"reasoning_tokens":3661}},"tokens_in":617,"tokens_out":3749,"duration_ms":26916,"temperature":1.0,"reasoning_tokens":3661,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:29:18.214384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $\\sum_f \\Gamma_{\\chi f\\to\\chi f}/H(T)$ at $T=m_\\chi/25$ for a relic-allowed benchmark with $\\lambda_{\\chi H}=10^{-3}$ and $m_\\chi$ in the $10$\\textendash$50$ MeV range; if this ratio drops below about one, dark and visible temperatures decouple and the Boltzmann solutions that assume $T_{\\rm dark}=T_{\\rm SM}$ miscompute freeze-out.","supporting_citations":[{"cited_title":"SIMPler realisation of Scalar Dark Matter","cited_arxiv_id":"1904.07562","evidence_quote":"Supplies the minimal complex-scalar $\\mathbb{Z}_3$ SIMP model, its Lagrangian, and the semi-analytic freeze-out solution that the two-component model builds on."},{"cited_title":"A multi-component SIMP model with $U(1)_X \\rightarrow Z_2 \\times Z_3$","cited_arxiv_id":"2103.05956","evidence_quote":"Gives a multi-component SIMP model with $\\mathbb{Z}_2\\times\\mathbb{Z}_3$ structure and the effective self-interaction cross-section formula used for the two-component case."}],"review_version":1}