{"id":"46925a0a-eb8c-4881-aac3-e09240f68e29","arxiv_id":"2506.11185","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new calculation of the Landau-Pomeranchuk-Migdal effect for scalar dark matter production shows that it can add up to roughly 27 percent to the predicted relic abundance.","lead":"This paper calculates how dark matter is produced in the early universe when many soft collisions in the hot plasma boost the production rate, an effect previously studied only for fermionic dark matter. It derives this effect for a scalar dark matter candidate and estimates that it can change the predicted dark matter abundance by up to about 27 percent.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1–27% LPM contribution is conditional on f_ζ: the paper's own switch-off comparison shows a 30% spread at δ=0.1, comparable to the effect it claims, and no scheme is valid at T∼mF0.","rationale":"The reader's weakest-assumption analysis and my stress-test converge on the same point: Eq. (3.25) defines a phenomenological switch-off that is unproven in the T∼mF0 regime, and the paper explicitly says so. This matters because the freeze-in relic density is an integral over temperature, and the switch-off is active where the integral can receive significant contributions, especially for larger mass splittings where the decay peak sits at z∼O(1). The paper's own comparison of three switch-off prescriptions gives a spread comparable to the headline LPM effect, so the headline '1–27%' should not be read as a controlled prediction of the scalar LPM contribution alone. I do not see a stronger internal inconsistency: the redefinition leading to Eq. (3.11) is algebraically consistent, the Born-limit subtraction is a sensible way to avoid double counting, and the comparison with HTL and semi-classical approaches is a useful addition. The transparency about the limitation is a genuine credit, and the proposal of a follow-up calculation with vacuum masses is appropriate. For that reason I would keep the reader's CONDITIONAL verdict rather than move to ACCEPT or REJECT: the paper delivers a first scalar LPM equation and a clear framework, but the central quantitative range is conditional on an interpolation whose systematic error is not yet controlled.","tokens_in":25846,"tokens_out":11591,"duration_ms":136077,"concrete_test":"Compute the z-integral in Eq. (4.3) for the d_R realization at (G=1.2, δ=0.1) and (G=1.2, δ=10), splitting the LPM correction term (γ_LPM − γ_LPM_Born) f_ζ(mF0/T) into a z<1 part, where f_ζ≈1 and the collinear LPM treatment is valid, and a z>1 part, where the switch-off is active and unvalidated. If the z>1 part contributes more than ~10% of the total LPM correction at either benchmark, the headline 1–27% range is not robust against the switch-off prescription; if the z>1 part is negligible, the switch-off is not load-bearing for those benchmarks.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim (Sec. 4.2, Fig. 10) is computed with Eq. (3.22): γ = γ_1PI + (γ_LPM − γ_LPM_Born) f_ζ(mF0/T), where f_ζ is defined in Eq. (3.25) as the normalized derivative of the fermionic thermal-loop function dJ_F/dz². That object controls the real part of the scalar self-energy in the static limit; it is used as a proxy for the suppression of the imaginary part (the LPM rate) without a derivation. The LPM power counting (3.2) assumes p∥∼T, p⊥∼gT, and masses ∼gT, so the collinear equation (3.11) is not valid for mF0∼T. The freeze-in integrand is non-negligible in exactly this regime, and the paper concedes (Sec. 3.4, Sec. 5) that none of the available switch-off schemes is strictly valid there. The quantitative impact is not small: switching between f_ζ, f_κ and the integrand-level interpolation changes the final relic density by up to ~30% at δ=0.1 and ~9% at δ=10 (Sec. 3.4), while the claimed LPM effect in the same regions is 27% and 8%. The headline range is therefore not a measurement of the LPM effect but a statement conditioned on one conservative interpolation; the uncertainty associated with the switch-off is of the same order as the signal.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the 1PI-resummed freeze-in calculation for scalar dark matter of Ref. [1] by including the Landau-Pomeranchuk-Migdal (LPM) effect. It derives a new integral equation for the LPM rate of a scalar particle, Eq. (3.11), computes its Born limit, and combines the LPM rate with the previous 1PI rate using a phenomenological switch-off function f_zeta defined in Eq. (3.25). On this basis the authors report that the LPM contribution changes the scalar DM relic density by 1--27%, with larger effects for larger gauge couplings and smaller mass splittings, and they compare their results with HTL-based and semi-classical Boltzmann approaches.","tokens_in":26154,"tokens_out":3275,"duration_ms":35146,"significance":"If the central numerical claim is accepted, the paper would be a useful step beyond the 1PI-resummed treatment of Ref. [1], providing the first explicit LPM resummation equation for a scalar FIMP and a detailed comparison of common approximations. The paper is transparent about the ad hoc nature of the switch-off and quantifies its uncertainty by comparing three schemes. It also provides fit functions and a scan over model realizations, which are valuable for phenomenological use. However, the central quantitative result, the 1--27% range in Sec. 4.2, is conditioned on a switch-off function whose theoretical status is not derived, and the paper's own comparison shows that the associated uncertainty is comparable to the claimed effect. The derivation of the LPM equation is sketched rather than fully self-contained, and the regime T ~ m_F0 where freeze-in is most sensitive is precisely the regime where the collinear LPM approximation is not valid.","major_comments":[{"comment":"The switch-off function f_zeta is defined as the normalized derivative of the fermionic thermal loop function, which controls the real part of the scalar self-energy in the static limit, yet it is used to suppress the imaginary part that constitutes the LPM rate. No derivation or justification is given for why this real-part suppression is an adequate proxy for the LPM imaginary-part suppression. The paper's own comparison (Sec. 3.4, Fig. 5) shows that switching between f_zeta, f_kappa and the integrand-level interpolation changes the final relic density by about 30% at delta=0.1 and about 9% at delta=10, while the claimed LPM effect in the same regions is 27% and 8%. Therefore the headline range in Sec. 4.2 is not a robust measurement of the LPM contribution; it is a statement conditioned on one specific interpolation, and the central quantitative claim should either be accompanied by an explicit systematic uncertainty band or be reformulated as a range that accounts for the switch-off ambiguity.","section":"Sec. 3.4, Eq. (3.25); Sec. 4.2, Fig. 10"},{"comment":"The LPM integral equation is derived under the power-counting assumption p_parallel ~ T, p_perp ~ gT, and masses ~ gT, as stated in Eq. (3.2). This assumption is not satisfied when m_F0 ~ T, and the freeze-in integrand is non-negligible in exactly this intermediate regime, as the paper itself emphasizes at the beginning of Sec. 1 and in Sec. 3.4. The paper concedes that none of the three switch-off schemes is strictly valid in the T ~ m_F0 region most relevant for freeze-in. Without a controlled treatment of that regime, the quantitative impact on the relic density, including the claimed 8--27% effects, is not established to the precision that the abstract and conclusion suggest. The authors should clearly restrict the domain of validity of the numerical results or supply an additional systematic check, such as a comparison with a calculation that retains vacuum masses in the LPM kernel, even in a simplified model.","section":"Sec. 3.1, Eq. (3.11); Sec. 3.4"},{"comment":"The Born-limit subtraction in Eq. (3.22) is an essential step: it removes the collinear decay contribution from the LPM rate to avoid double counting with the 1PI rate. The paper states in Sec. 3.3 that the LPM rate and its Born limit coincide in the non-relativistic regime, and treats this as a consistency check. However, both quantities are computed from the same collinear approximation, so this agreement does not independently validate the subtraction in the transition region. The authors should clarify what independent check, if any, exists for the subtraction at intermediate temperatures, and how sensitive the final relic density is to the exact form of the subtraction as opposed to the switch-off function.","section":"Sec. 3.2, Eq. (3.18); Sec. 3.3, Eq. (3.22)"}],"minor_comments":[{"comment":"In the sentence 'such as the LMP effect' the acronym should be 'LPM', not 'LMP'.","section":"Sec. 2"},{"comment":"In the text following Eq. (3.22), the switch-off function is written as f(m_F,0F/T); this appears to be a typo for f(m_F,0/T).","section":"Sec. 3.3"},{"comment":"The caption states that the plots consider a mediator of type d_R with G=1.2 and mass splitting delta=10, but the left and right panels correspond to different mass splittings (delta=0.1 and delta=10, respectively, as discussed in the text). The caption should be clarified.","section":"Fig. 5"},{"comment":"The smooth interpolation method is attributed to Ref. [34] in Sec. 3.4 and to Ref. [56] in Sec. 5; please verify that these citations refer to the correct and distinct works, and make the attribution consistent.","section":"Sec. 3.4, Sec. 5"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this paper if you care about scalar FIMP freeze-in rates: it gives the first scalar LPM integral equation and combines it with the 1PI-resummed framework from Ref. [1]. The derivation is transparent: collinear power counting, recursion relation adapted from the fermionic case, Born-limit consistency check. The comparison of three switch-off schemes is genuinely useful, and the paper is admirably explicit that none is valid at T~M.\n\nThe central quantitative claim, a 1-27% LPM contribution to the relic density, is conditional on the chosen switch-off f_zeta. The stress-test note is on target: Eq. (3.25) uses the derivative of the fermionic thermal loop function, which controls the real part of the scalar self-energy in the static limit, as a proxy for the suppression of the imaginary part; there is no derivation for that step. The power counting (3.2) breaks down for mF0~T, exactly where freeze-in is most relevant. And the paper's own numbers show a ~21% spread between f_zeta and f_kappa on the LPM-only rate, and ~30% difference between f_zeta and smooth interpolation for the complete rate at delta=0.1. That is the same size as the headline 27% effect.\n\nWhat saves this from being a fatal problem is honesty. Sec. 3.4 states in plain terms that these are interpolation schemes, not a valid calculation in the relativistic regime; Sec. 5 says a full mass-dependent LPM treatment is follow-up work. The paper does not oversell. The 1-27% range should be read as \"given our conservative f_zeta choice,\" and a careful abstract would say so more clearly.\n\nOther soft spots are minor: the gauge dependence of the 1PI rate is discussed but only estimated at the percent level; the fit function for the LPM rate is limited to the grid scanned; and the comparison with Boltzmann approaches is useful but inherits all the switch-off ambiguity.\n\nVerdict: solid, worthwhile, should go to peer review. A referee should ask for the 27% to be reported together with the switch-off spread, and for a clear statement that the T~M region is interpolated, not computed. But the new scalar LPM equation is a real contribution, and the paper's uncertainty accounting is better than most.","headline":"First scalar LPM equation combined with 1PI-resummed freeze-in rate; the 1-27% effect is real but conditional on a switch-off whose uncertainty at small splitting is comparable to the effect—and the authors say so.","tokens_in":26719,"tokens_out":2472,"would_cite":true,"duration_ms":26217,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Including multiple soft plasma scatterings changes scalar dark-matter freeze-in relic densities by up to 27 percent.","keywords":["scalar dark matter","freeze-in","Landau-Pomeranchuk-Migdal effect","1PI resummation","thermal field theory","relic density","multiple soft scattering","FIMP"],"falsifier":"Compute the LPM production rate with the physical vacuum masses of the mediator and dark matter kept in the $T\\sim M$ regime, without invoking a switch-off function; if the resulting relic density differs from the thermal-function result by more than the 30 percent spread among switch-off schemes, the paper's quoted 1 to 27 percent LPM contribution would not hold in that regime. This is precisely the calculation the paper identifies as the missing step.","tokens_in":25616,"feed_emoji":"🌌","tokens_out":9027,"duration_ms":93031,"temperature":0.7,"pith_summary":"This paper improves the calculation of scalar dark matter production by freeze-in, the process in which dark matter accumulates from rare decays and scatterings without ever reaching thermal equilibrium. It claims that a previously missing piece of the production rate, the Landau-Pomeranchuk-Migdal (LPM) effect of coherent multiple soft scatterings in the plasma, is a leading-order contribution of size $g^2 T$ and must be added to the earlier calculation based on one-particle-irreducible (1PI) resummed propagators. The paper derives for the first time an integral equation for the LPM rate of a scalar particle, Eq. (3.11), and combines it with the 1PI result. The resulting relic density is 1 to 27 percent larger than the 1PI-only result, with the largest enhancement for larger gauge couplings and smaller mass splittings. The result matters because most freeze-in predictions so far omitted this leading-order effect, and the paper calibrates simpler approximations against the improved rate.","feed_headline":"Soft-scattering effect lifts scalar dark matter yield by up to 27%","feed_subtitle":"New freeze-in rate includes coherent multiple plasma scatterings; effect peaks for large gauge couplings and small mass splittings.","key_machinery":"The load-bearing object is the scalar LPM integral equation, Eq. (3.11), for the function $\\chi(\\vec{k}_\\perp)$, together with the spectral self-energy Eq. (3.10) that it feeds. The equation balances the pole-distance term $\\epsilon(\\vec{k}_\\perp)$, which contains the vacuum and thermal masses of the mediator and the Standard-Model fermion, against a collision integral with kernel $K_i(\\vec{q}_\\perp)=1/\\vec{q}_\\perp^2 - 1/(\\vec{q}_\\perp^2 + m_{D i}^2)$, so that arbitrarily many soft scatterings are summed at the same parametric order $g^2 T$ as the one-loop rate. The companion switch-off function $f_\\zeta(m_{F,0}/T)$, defined from the derivative of a fermionic thermal loop function, suppresses the LPM contribution as the temperature approaches the vacuum mass scale, where the collinear scaling $p_\\perp\\sim gT$ used in the derivation fails.","core_discovery":"The central claim is that the LPM effect contributes at leading order to scalar dark matter production and therefore changes the predicted relic density. The paper's new object is Eq. (3.11), the scalar analogue of the LPM resummation: an integral equation for $\\chi(\\vec{k}_\\perp)$ in which the free-propagation energy difference $\\epsilon(\\vec{k}_\\perp)$ is balanced against a collision term integrating over soft transverse gauge-boson momenta with a Debye-screened kernel. The spectral self-energy built from $\\chi$ in Eq. (3.10) resums an infinite ladder of soft scatterings between the two fermion lines of the production diagram. Because this LPM derivation assumes collinear ultra-relativistic kinematics, the paper attaches a phenomenological switch-off function $f_\\zeta(m_{F,0}/T)$ and subtracts the LPM Born limit to avoid double counting the decay. For five charge assignments, over $G\\in[0.4,1.6]$ and $\\delta\\in[0.1,10]$, the LPM contribution increases the relic density by 1 to 27 percent, and the paper reports that the choice of switch-off scheme introduces up to 30 percent uncertainty for small mass splittings.","pith_inferences":["If the switch-off ambiguity is the dominant residual uncertainty, then freeze-in relic densities in the small-splitting, large-coupling region should be quoted as a band rather than a single number.","A full LPM calculation that keeps the physical mediator and dark matter masses in the $T\\sim M$ regime, which the paper leaves to future work, could move the central relic density outside the quoted 1 to 27 percent range.","The thermal switch-off function, tied to the scalar thermal mass, might also apply to other scalar self-energy processes such as axion-like particle production from a thermal plasma.","The model dependence through gauge charges suggests that the same resummation could be used to rank mediator representations by production efficiency before scanning couplings."],"forward_implications":["Scalar singlet dark matter produced from gauge-charged fermions should have its freeze-in relic density revised upward by 1 to 27 percent, with the correction largest for strongly coupled, nearly degenerate dark sectors.","The scalar LPM integral equation provides a reusable tool for any scalar production process mediated by fermions in a plasma, not only dark matter freeze-in.","Common semi-classical treatments are now benchmarked against the improved rate: vacuum-mass decays stay within about 10 percent for large mass splittings, while thermal-mass decays plus scatterings can deviate by roughly -30 to +25 percent.","The HTL-resummed 1PI rate without LPM happens to match the full result to about 10 percent across the parameter scan, because its overestimate is offset by the missing LPM enhancement.","The LPM rate depends on the specific gauge charges of the mediator, unlike the 1PI one-loop rate, so relic-density predictions must be made separately for each realization such as $q_L$, $e_R$, or $d_R$."],"supporting_citations":[{"why":"Supplies the 1PI-resummed production rate to which the LPM contribution is added.","marker":"[1]"},{"why":"Provides the fermionic DM freeze-in framework with LPM and the susceptibility switch-off used as a comparison benchmark.","marker":"[12]"},{"why":"Supplies the recursion relation for ladder diagrams from which the scalar LPM integral equation is derived.","marker":"[36]"},{"why":"Establishes the power counting showing that multiple soft scatterings contribute at leading order.","marker":"[28]"},{"why":"Provides the susceptibility-based switch-off prescription and the neutrino LPM context.","marker":"[35]"},{"why":"Motivates the thermal loop function whose derivative defines the new switch-off function.","marker":"[40]"},{"why":"Supplies the smooth interpolation scheme for dilepton rates used as one of the switch-off alternatives.","marker":"[34]"},{"why":"Provides the phase-space and scattering expressions used in the HTL-based comparison rate.","marker":"[33]"}],"fun_headline_variants":["LPM effect boosts scalar DM freeze-in by up to 27%","Multiple soft scatterings raise scalar DM relic density","Scalar DM freeze-in gets up to 27% boost from LPM","LPM correction boosts scalar DM production by 27%","Soft scattering LPM effect adds up to 27% to scalar DM yield"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the phenomenological switch-off function used to turn off the LPM contribution at low temperature reliably interpolates between the ultra-relativistic regime, where the LPM calculation is valid, and the non-relativistic regime, where the 1PI rate should take over, even though the paper states that no available switch-off method is strictly valid in the $T\\sim M$ region most relevant for freeze-in and that the choice among schemes changes the final relic density by up to 30 percent for small mass splittings.","fun_headline_variants_meta":{"raw":{"variants":["LPM effect boosts scalar DM freeze-in by up to 27%","Multiple soft scatterings raise scalar DM relic density","Scalar DM freeze-in gets up to 27% boost from LPM","LPM correction boosts scalar DM production by 27%","Soft scattering LPM effect adds up to 27% to scalar DM yield"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00099,"raw_usage":{"total_tokens":4272,"prompt_tokens":1099,"completion_tokens":3173,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":715,"completion_tokens_details":{"reasoning_tokens":3083}},"tokens_in":715,"tokens_out":3173,"duration_ms":22775,"temperature":1.0,"reasoning_tokens":3083,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:12:01.649880+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the LPM production rate with the physical vacuum masses of the mediator and dark matter kept in the $T\\sim M$ regime, without invoking a switch-off function; if the resulting relic density differs from the thermal-function result by more than the 30 percent spread among switch-off schemes, the paper's quoted 1 to 27 percent LPM contribution would not hold in that regime. This is precisely the calculation the paper identifies as the missing step.","supporting_citations":[{"cited_title":"Biondini and J","cited_arxiv_id":null,"evidence_quote":"Provides the fermionic DM freeze-in framework with LPM and the susceptibility switch-off used as a comparison benchmark."},{"cited_title":"Besak and D","cited_arxiv_id":null,"evidence_quote":"Supplies the recursion relation for ladder diagrams from which the scalar LPM integral equation is derived."},{"cited_title":"Arnold, G.D","cited_arxiv_id":null,"evidence_quote":"Establishes the power counting showing that multiple soft scatterings contribute at leading order."},{"cited_title":"Ghiglieri and M","cited_arxiv_id":null,"evidence_quote":"Provides the susceptibility-based switch-off prescription and the neutrino LPM context."},{"cited_title":"Ghisoiu and M","cited_arxiv_id":null,"evidence_quote":"Supplies the smooth interpolation scheme for dilepton rates used as one of the switch-off alternatives."},{"cited_title":"Besak and D","cited_arxiv_id":null,"evidence_quote":"Provides the phase-space and scattering expressions used in the HTL-based comparison rate."}],"review_version":1}