{"id":"746c9d93-ab62-4a97-b495-c7a81bc60c80","arxiv_id":"2608.12079","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A Sudakov resummation of dark-sector radiation converts the divergent M_X^2 distribution in invisible dark photon production into an integrable line shape and shifts recast BaBar limits by up to about 14%.","lead":"Dark photons produced at electron-positron colliders can radiate extra invisible dark-sector particles before decaying, and this paper calculates how that radiation reshapes the missing-mass signal used in dark photon searches. The new resummed prediction is integrable where the fixed-order one diverges, and it changes the recast BaBar limit on the kinetic mixing parameter by up to about 14%.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The resummed line shape and the claimed ~14% limit shift depend on setting Q_hard^2 ≈ s; under the BaBar LowM cut the invisible branch is kinematically bounded by q_cut^2 ≈ 40 GeV^2, so the hard-scale choice needs a band.","rationale":"I read the paper as a standard Sudakov treatment of a collinear endpoint, and I did not find an algebraic inconsistency in the primary exponentiation: Eq. (79) is the derivative of the no-emission survival factor and integrates to a finite result. The secondary Sudakov factor and the matching/beta_peak construction are internally consistent, and the authors flag the sum-veto simplification in a footnote. The one condition that is both load-bearing and under-specified is the hard scale Q_hard^2. Because the BaBar LowM cuts restrict q to q_cut^2 ≈ 40 GeV^2, while the paper motivates Q_hard^2 ~ s, the recast numbers in Table IV and Fig. 5 have an unquantified scale sensitivity. This is exactly the kind of condition a conditional acceptance should require, so the reader's CONDITIONAL verdict stands. I set agreement_with_reader to partial because I focus on the second of the reader's two flagged premises rather than the mass hierarchy of Eq. (71), which is satisfied for the massless benchmarks used in the recast.","tokens_in":29906,"tokens_out":18199,"duration_ms":183104,"concrete_test":"Recompute Eq. (117) and Table IV with Q_hard^2 set to q_cut^2 = s - 2 E_gamma sqrt(s) ≈ 40 GeV^2 instead of s, plus a scale band Q_hard^2 in {s/4, s, 4s}, for the three charge benchmarks and alpha' = 0.5. If any epsilon'/epsilon_0 entry moves by more than about 0.02, the 14% claim needs an explicit scale-dependence statement; if all entries move by less, the hard-scale choice is not load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the value of Q_hard^2 in the Sudakov exponent. Equation (71) states Q_hard^2 ~ s, and Eq. (78) defines the primary factor Δ_pri = exp[-a_pri ln(Q_hard^2/M_X^2)]. The resummed distribution Eq. (79) is a_pri/M_X^2 (M_X^2/Q_hard^2)^{a_pri}, and the matched r_resum(q) of Eq. (116) inherits this scale. But in the BaBar LowM selection used in Sec. VI, E_gamma > 3 GeV and sqrt(s) = 10.58 GeV bound the invisible branch invariant mass by q <= q_cut^2 = s - 2 E_gamma sqrt(s) ≈ 40 GeV^2, not by s ≈ 112 GeV^2. Since the matched distribution is normalized through beta_peak in Eq. (122), using s instead of q_cut^2 changes the low-q weight and the detector-level shape. For alpha' = 0.5, a_pri ≈ 0.066, so the effect is moderate but not negligible: the integrated resummed weight below q_2 scales like (q_2/Q_hard^2)^{a_pri}, which differs by about 6% between Q_hard^2 = 112 and 40 GeV^2. Table IV quotes epsilon'/epsilon_0 up to 1.137 with no Q_hard uncertainty band, so the claim that dark-sector radiation weakens BaBar limits by up to ~14% is not yet robust.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies e+e- -> gamma A' production in a dark Abelian Higgs model with a light dark fermion chi, computing NLO virtual and real corrections from the dark Higgs and chi sectors. It shows that the inclusive NLO cross section is IR-safe after resonance subtraction, identifies a quasi-collinear 1/M_X^2 enhancement in the differential missing-mass-squared distribution, and resums it with a primary plus secondary Sudakov treatment to obtain an integrable, normalized line shape. The matched particle-level distribution is then convolved with a Crystal-Ball detector response and used in a simplified recast of the BaBar mono-photon search, with the result that dark-sector radiation can modify the inferred epsilon limit by up to about 14% for alpha'=0.5.","tokens_in":30281,"tokens_out":16562,"duration_ms":177572,"significance":"If the Sudakov treatment is accepted, the paper provides the first analytic resummed M_X^2 prediction for invisible dark-photon production in this model, together with an explicit inclusive IR-safety check and a transparent matching/recast pipeline. The derivations in Eqs. (44)-(53), the collinear coefficient in Eq. (67), and the normalized endpoint distribution in Eq. (79) are internally consistent, and the paper does not fit parameters to data. The main phenomenological claim is, however, sensitive to the choice of hard scale and to the treatment of intermediate dark-photon masses, so the numerical conclusion is not yet fully robust.","major_comments":[{"comment":"The hard scale Q_hard^2 is identified with s (Eq. (71)) without a specified value or uncertainty band, but the kinematic endpoint of the BaBar LowM selection is qmax = s - 2 E_cut^gamma sqrt(s) ~ 48 GeV^2 at sqrt(s)=10.58 GeV, not s ~ 112 GeV^2. The primary-resummed distribution r_resum(q) = C_sec (a_pri/q)(q/Q_hard^2)^{a_pri} and the normalization beta_peak in Eq. (122) depend on Q_hard^2. For alpha'=0.5, a_pri ~ 0.066, and changing Q_hard^2 from 112 to 48 GeV^2 changes the integrated Sudakov weight below q2 by roughly (48/112)^{a_pri} ~ 0.945, a 5.5% effect of the same order as the claimed limit shifts in Table IV. The manuscript should either adopt qmax as the hard scale or provide a scale-variation band for epsilon'/epsilon_0 before the ~14% statement can be considered robust.","section":"Section V.A and Section VI.D (Eqs. (78), (116), (122))"},{"comment":"The recast uses the massless matched spectrum only for mA' <= 0.1 GeV and the finite-mass NLO distribution for mA' >= 0.5 GeV. For mA' >= 0.5 GeV, m_A'^2 >= 0.25 GeV^2, so the hierarchy m_A'^2 << M_X^2 in Eq. (71) is not satisfied in the low-M_X tail, and the fixed-order quasi-collinear 1/M_X^2 divergence is left un-resummed in exactly the mass region displayed in Fig. 5. Since the low-M_X tail feeds the binned likelihood through Eq. (128), the plotted limit ratios for mA' ≳ 0.5 GeV rest on an un-resummed distribution. Please either extend the resummation with a mass-dependent low-scale cutoff or restrict the resummed limit claim to the massless regime and clearly label the high-mass extension as an un-resummed estimate.","section":"Section VI.D (Figs. 4-5, Eq. (121))"}],"minor_comments":[{"comment":"For the BaBar LowM selection, qmax = s - 2 E_cut^gamma sqrt(s) is approximately 48.5 GeV^2 at sqrt(s)=10.58 GeV, not 40 GeV^2; the numerical values used in the text and Table IV should be made consistent.","section":"Section VI.A"},{"comment":"With sqrt(s)=10.58 GeV, q1=0.01 s = 1.12 GeV^2 and q2=0.10 s = 11.2 GeV^2; if the rounded values 1 and 10 GeV^2 are intended, this should be stated explicitly.","section":"Eq. (119)"},{"comment":"The statement that the final conclusion is insensitive to q_det is not demonstrated; a short scan over q_det or an argument based on the detector resolution would support this claim.","section":"Section VI.C"},{"comment":"The recast uses digitized BaBar background expectations from a fitted curve without propagating the digitization uncertainty; a validation against the official bin-by-bin workspace or an estimate of the digitization error would improve the robustness of S_fit and epsilon'/epsilon_0.","section":"Section VI.D"},{"comment":"The leading-log treatment of the secondary Sudakov factor replaces the sum constraint on t_sec,i/zeta_i by individual constraints (Eq. (90)) and relies on angular ordering; the numerical precision quoted for C_sec in Table I should be accompanied by an estimate of the uncertainty from these approximations.","section":"Section V.C"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. The paper gives the first analytic Sudakov-resummed M_X^2 distribution for invisible dark photon production in a dark Abelian Higgs model, including a light dark fermion channel. That part is a genuine extension of the authors' earlier fixed-order and shower work. The headline numerical claim—dark-sector radiation can weaken BaBar limits by up to about 14%—is not yet robust, because it rests on identifying the hard scale Q_hard^2 with s, while the BaBar LowM cut actually bounds the invisible branch at q_cut^2 ≈ 48 GeV^2.\n\nThe core physics argument is sound. The fixed-order differential distribution has a quasi-collinear 1/M_X^2 divergence when dark-sector masses are small; the primary Sudakov factor converts that into an integrable, normalized endpoint distribution. The IR-safety analysis in Sec. III is consistent, and the collinear coefficient in Eq. (67) is right. I also find the secondary-radiation analysis convincing: the dipole cancellation means secondary emissions contribute at order a few percent at most, so the primary Sudakov factor is a good approximation.\n\nThe soft spot is the hard scale. The resummed weight below q2 scales like (q2/Q_hard^2)^(a_pri); for a_pri ≈ 0.066 that differs by about 6% between Q_hard^2 = 112 and 48 GeV^2. Since β_peak is fixed by K_NLO minus the integral of the matched distribution, this shifts the detector-level template and the ε'/ε0 ratios in Table IV by a few percent. No uncertainty band is given. The matching scales q1, q2 are set to 0.01s and 0.10s without a sensitivity scan, and the BaBar background is digitized and treated as fixed. These are fixable, but they should be addressed.\n\nI don't see a load-bearing flaw. The exponentiation is standard, the self-citations to [70,83] are normal, and the central claim that the divergence is Sudakov-suppressed survives. The paper would benefit from a scale-variation estimate; with that, it would be a solid contribution.\n\nThis is for dark-photon and dark-sector phenomenologists, particularly people working on mono-photon recasts. It deserves a serious referee. I'd send it to review and ask for a hard-scale variation and a scan over the matching scales.","headline":"Sudakov resummation of the dark photon missing-mass distribution is a real advance, but the claimed 14% BaBar shift needs a hard-scale uncertainty band.","tokens_in":30836,"tokens_out":7311,"would_cite":true,"duration_ms":69100,"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 establishes that the quasi-collinear $1/M_X^2$ singularity in the invisible dark photon missing-mass distribution is resummed by a primary Sudakov factor into an integrable, normalized jet-mass line shape, and that including…","keywords":["dark photon","dark Abelian Higgs model","invisible dark photon production","missing mass distribution","jet mass","Sudakov resummation","collinear divergence","mono-photon search"],"falsifier":"Measure the mono-photon missing-mass spectrum in a high-statistics $e^+e^-$ run at $\\sqrt{s}=10$ GeV with $m_{A'}\\ll\\sqrt{s}$ and known $\\alpha'$: the resummed prediction makes the cumulative distribution scale as $(M_X^2/Q_{\\rm hard}^2)^{a_{\\rm pri}}$, with $a_{\\rm pri}=(\\alpha'/\\pi)(Q_S^2/12+Q_\\chi^2/3)$, whereas the fixed-order prediction would show the unintegrated $a_{\\rm pri}/M_X^2$ power law down to the detector resolution.","tokens_in":29655,"feed_emoji":"📡","tokens_out":13908,"duration_ms":121761,"temperature":0.7,"pith_summary":"This paper asks what happens to the $M_X^2$ distribution of mono-photon events when the dark photon lives in a dark Abelian Higgs model with a light dark Higgs and a light dark fermion. It shows that the total NLO cross section is infrared-safe after resonance subtraction, but the differential distribution develops a quasi-collinear $1/M_X^2$ divergence in the region where all dark-sector masses are small. It then proves that a primary Sudakov factor, the no-radiation probability of the final-state dark photon, converts this divergence into an integrable, normalized line shape that is the jet mass of the dark photon branch. The result is a resummed prediction that, after matching to fixed order and detector response, changes the inferred kinetic-mixing limits by up to about 14%.","feed_headline":"Resumming dark radiation shifts dark photon limits by up to 14%","feed_subtitle":"The missing-mass singularity becomes an integrable jet-mass line shape, changing recast searches at e+e- colliders.","key_machinery":"The central object is the squared missing mass $M_X^2$, interpreted as the jet mass of the invisible dark-photon branch. The load-bearing mechanism is the primary Sudakov factor $\\Delta_{\\rm pri}$, the no-radiation probability between the hard scale $Q_{\\rm hard}^2\\sim s$ and the measured scale $M_X^2$, which multiplies the fixed-order singular distribution and converts $(M_X^2)^{-1}$ into an integrable $(M_X^2)^{-1+a_{\\rm pri}}$ endpoint enhancement. The secondary Sudakov factor treats soft radiation from the daughters of the primary splitting, with an angular-ordering constraint that enforces the dipole cancellation removing the would-be double-log contribution. An additive matching formula, fixed order plus resummed minus the singular overlap, joins the low-$M_X^2$ resummed region to the high-$M_X^2$ hard region.","core_discovery":"The central discovery is that the fixed-order quasi-collinear NLO distribution, $$ \\frac{1}{\\sigma_0}\\frac{d\\sigma_{\\rm NLO}}{$dM_X^{2}$} \\simeq \\frac{a_{\\rm pri}}{$M_X^{2}$}, \\qquad a_{\\rm pri} = \\frac{\\$\\alpha$'}{\\pi}\\left(\\frac{$Q_S^{2}$}{12}+\\frac{Q_\\$chi^{2}$}{3}\\right), $$ is converted by the primary Sudakov factor $\\Delta_{\\rm pri}(Q_{\\rm hard}^2,M_X^2)=\\exp[-a_{\\rm pri}\\ln(Q_{\\rm hard}^2/M_X^2)]$ into the integrable, normalized line shape $$ \\frac{1}{\\sigma_0}\\frac{d\\sigma_{\\rm pri}^{\\rm resum}}{$dM_X^{2}$} = \\frac{a_{\\rm pri}}{$M_X^{2}$}\\,\\exp\\!\\left[-a_{\\rm pri}\\ln\\frac{Q_{\\rm hard}^2}{$M_X^{2}$}\\right]. $$ The same coefficient $a_{\\rm pri}$ controls the singularity and the exponent because the primary splitting kernels $P_{T\\to Ls}(z)=z(1-z)$ and $P_{T\\to\\chi\\bar\\chi}(z)=z^2+(1-z)^2$ have no soft endpoint pole. Secondary radiation from the daughter system is suppressed by a dipole angular-ordering cancellation, so the secondary factor $C_{\\rm sec}$ is within about 2% of unity even at $\\alpha'=0.5$. After additive matching to the fixed-order hard region and convolution with detector response, the predicted spectrum shifts recast kinetic-mixing limits by up to about 14%.","pith_inferences":["A high-statistics measurement of the mono-photon $M_X^2$ spectrum at a future $e^+e^-$ collider could extract the combination $Q_S^2/12+Q_\\chi^2/3$ from the shape alone, separating the dark-sector coupling structure from the overall kinetic-mixing scale.","The factorization $d\\sigma_{\\gamma+X}\\simeq d\\sigma_{\\gamma+A'^*}\\times dP_{A'\\to X}$ should apply to other invisible dark-sector final states with a collinear tail, such as invisible dark Higgs or $Z'$ production, so the machinery developed here is a template for resumming those line shapes.","The size of the recast shift depends on identifying $Q_{\\rm hard}^2$ with $s$; varying $Q_{\\rm hard}^2$ over a plausible range would bracket an additional systematic uncertainty on the 14% estimate that the paper does not quantify."],"forward_implications":["The $M_X^2$ distribution is integrable and normalized after the primary Sudakov resummation, so every finite bin of the missing-mass spectrum is a well-defined prediction even though the fixed-order distribution diverges as $1/M_X^2$.","Secondary soft radiation modifies the line shape by less than about 2% even at $\\alpha'=0.5$, so the leading-log primary-resummed result is already a good approximation for the endpoint region.","Dark-sector radiation broadens the reconstructed missing-mass signal, which generally weakens the mono-photon constraint on the kinetic-mixing parameter: up to about 10% for equal dark charges and up to about 14% for charge assignments $(Q_S,Q_\\chi)=(1/4,1)$ at $\\alpha'=0.5$.","The matched particle-level distribution, which combines the resummed low-$q$ region with the fixed-order high-$q$ region and subtracts the singular overlap, connects smoothly to the hard region and can be convolved with detector response.","The same primary radiation coefficient $a_{\\rm pri}$ controls both the fixed-order singularity and the resummation exponent, so the shape and the rate correction are tied to the same combination of dark-sector charges."],"supporting_citations":[{"why":"Supplies the dark final-state radiation splitting kernels and the mono-photon recast framework used for the detector-level analysis.","marker":"[70]"},{"why":"Establishes the infrared-safe inclusive NLO cross section in the same model and provides the dark-Higgs virtual self-energy expression that the resummation builds on.","marker":"[83]"},{"why":"Provides the experimental mono-photon search dataset, background expectations, and detector-response parameters used in the recast.","marker":"[44]"},{"why":"Introduces the Sudakov method used to resum the collinear logarithms into an integrable endpoint distribution.","marker":"[86]"},{"why":"States the mass-singularity cancellation theorem used to show that the inclusive NLO rate is finite after real and virtual contributions are combined.","marker":"[84]"},{"why":"Supplies the companion degeneracy theorem that justifies the same inclusive cancellation.","marker":"[85]"},{"why":"Provides the general final-state resummation and matching principles used by the additive matching formula that joins resummed and fixed-order spectra.","marker":"[93]"}],"fun_headline_variants":["Sudakov resummation tames dark photon missing-mass singularity","Dark sector radiation tamed: limits shift by up to 14%","Jet-mass line shape from dark photon production","Resummed dark photon spectra alter recast limits by 14%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the strong mass hierarchy $m_{A'}^2, m_s^2, m_\\chi^2 \\ll M_X^2 \\ll Q_{\\rm hard}^2 \\sim s$, which lets every dark-sector particle be treated as massless in the collinear region; if the dark Higgs or dark photon is not light compared with the measured missing mass, the $1/M_X^2$ enhancement stops controlling the distribution and the resummed line shape no longer applies.","fun_headline_variants_meta":{"raw":{"variants":["Sudakov resummation tames dark photon missing-mass singularity","Dark sector radiation tamed: limits shift by up to 14%","Jet-mass line shape from dark photon production","Resummed dark photon spectra alter recast limits by 14%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1552,"prompt_tokens":1053,"completion_tokens":499,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":669,"completion_tokens_details":{"reasoning_tokens":428}},"tokens_in":669,"tokens_out":499,"duration_ms":5368,"temperature":1.0,"reasoning_tokens":428,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:18:12.995790+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the mono-photon missing-mass spectrum in a high-statistics $e^+e^-$ run at $\\sqrt{s}=10$ GeV with $m_{A'}\\ll\\sqrt{s}$ and known $\\alpha'$: the resummed prediction makes the cumulative distribution scale as $(M_X^2/Q_{\\rm hard}^2)^{a_{\\rm pri}}$, with $a_{\\rm pri}=(\\alpha'/\\pi)(Q_S^2/12+Q_\\chi^2/3)$, whereas the fixed-order prediction would show the unintegrated $a_{\\rm pri}/M_X^2$ power law down to the detector resolution.","supporting_citations":[{"cited_title":"Zheng, Y","cited_arxiv_id":null,"evidence_quote":"Establishes the infrared-safe inclusive NLO cross section in the same model and provides the dark-Higgs virtual self-energy expression that the resummation builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the companion degeneracy theorem that justifies the same inclusive cancellation."}],"review_version":1}