{"id":"bccce650-0a9c-4fa4-a4c4-c47fcb1ea8ba","arxiv_id":"2412.02536","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Using PHSD, the authors compute dark photon contributions to dilepton spectra and obtain epsilon^2(M_U) limits that scale with an assumed 0.3% to 10% surplus over the Standard Model yield.","lead":"This paper adds dark photon decay channels to a heavy-ion collision model and derives upper limits on the dark photon mixing parameter from the condition that their dilepton yield stays below an assumed extra fraction of known Standard Model decays. The result is a set of sensitivity curves showing what precision future SIS, HADES and RHIC experiments would need to compete with existing dark photon searches.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed upper limit on epsilon^2 is set by the hand-chosen surplus C_U=0.3%, tuned to match BaBar14, not derived from the heavy-ion dilepton data or its uncertainties.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: C_U is an input chosen to match BaBar14 rather than derived from the heavy-ion data, and since Eq. (15) is linear in C_U, any error or arbitrariness in C_U propagates directly into the claimed bound. My stress-test adds the concrete quantitative point that the chosen 0.3% level is far below the actual precision of the HADES/STAR dilepton spectra, so those data cannot exclude a surplus of 0.3% and therefore cannot support the quoted limit. The paper is internally consistent and transparent about its preliminary nature, and the PHSD shape of the bound is a valid sensitivity prediction. However, the central claim of an independent upper limit near the world constraint is overstated; the work should be reframed as a precision-requirements forecast. This is exactly the reader's CONDITIONAL verdict, so no change to the verdict is needed. The concern is load-bearing because if C_U is corrected to a data-derived value, the headline constraint shifts by orders of magnitude, changing the scientific message from 'heavy-ion data constrain dark photons at current world limits' to 'heavy-ion data would need 0.3% precision to reach those limits.'","tokens_in":13085,"tokens_out":5483,"duration_ms":56282,"concrete_test":"Use the published HADES and STAR dilepton spectra with their full systematic and statistical uncertainties to compute, for each mass bin, a 90% CL upper bound on (N_meas - N_SM)/N_SM; set C_U to this data-derived bound in Eq. (15) and recompute epsilon^2(M_U). If the resulting curve lies above ~10^-5 across 0.2 e+e- decay. For each invariant-mass bin, Eq. (15) converts an assumed fractional surplus C_U of the U-boson yield over the summed Standard Model dilepton yield into an upper limit on the kinetic mixing parameter epsilon^2(M_U). Using HADES data at SIS18 energies and STAR data at RHIC energies, the authors report that with C_U = 10% the extracted epsilon^2(M_U) is compatible with BaBar09 limits, and that with C_U = 0.3% the result matches the BaBar14 exclusion and the global compilation value epsilon^2 ~ 3 x 10^-7 for 0.2 < M_U < 1.2 GeV. The paper frames this as a procedure for constraining dark photon kinetic mixing with heavy-ion dilepton data.","tokens_in":13389,"tokens_out":3235,"duration_ms":32612,"significance":"If the derived limit were genuinely anchored to the precision of the heavy-ion dilepton measurements, the paper would offer an independent, energy-frontier probe of dark photon parameter space complementary to fixed-target and collider searches. The study's strengths are its use of the PHSD model, which reproduces the measured Standard Model dilepton spectra across several collision systems, the inclusion of multiple new U-boson production channels beyond earlier work, and the transparent scaling relation Eq. (15) that makes the sensitivity to C_U explicit. The paper is also honest about the missing QGP and charm sources at RHIC. However, the central quantitative claim---that heavy-ion data can exclude epsilon^2 down to ~3 x 10^-7---depends entirely on the ad hoc choice C_U = 0.3%, which is selected to match BaBar14 rather than determined from the experimental uncertainties of the HADES or STAR data. As a sensitivity forecast for future measurements, the study is useful; as a derived upper limit, the present form is not yet supported.","major_comments":[{"comment":"The central extraction of epsilon^2(M_U) is linearly proportional to the arbitrary surplus parameter C_U: Eq. (15) is epsilon^2 = C_U * (dN_sumSM/dM)/(dN_sumU@eps=1/dM). The value C_U = 0.3% is introduced in Sec. 4 (right panel of Fig. 3) specifically to make the resulting curve coincide with the BaBar14 bound, and a 10% surplus is used in the left panel. Since the paper does not determine C_U from the statistical and systematic uncertainties of the measured dilepton spectra, the claimed agreement with the world limit epsilon^2 ~ 3 x 10^-7 is imposed by input assumption rather than inferred from the heavy-ion data. The authors should either quote limits for a range of C_U values tied to the known experimental uncertainties, or perform a proper statistical comparison to the HADES/STAR data points (e.g., a chi-square or CLs procedure) with the uncertainties propagated onto epsilon^2.","section":"Sec. 4, Eq. (15)"},{"comment":"The statement that a dark photon contribution 'must remain within the uncertainties of the experimental data' is not operationalized. The manuscript does not provide the statistical or systematic uncertainties of the HADES and STAR dilepton spectra used in the comparison, nor does it assign an uncertainty band to the extracted epsilon^2(M_U) curves. Published HADES and STAR dilepton measurements have point-to-point uncertainties of several percent or more, far larger than the C_U = 0.3% used to claim sensitivity to epsilon^2 ~ 3 x 10^-7. Without an explicit mapping from experimental precision to C_U, the quoted limit is a statement of model sensitivity, not a constraint derived from the data.","section":"Sec. 4, paragraph after Eq. (15)"},{"comment":"The paper acknowledges that for Au+Au collisions at 19 and 200 GeV the mass region above 1 GeV is dominated by QGP and charm contributions that are not included in the model, and that the extracted epsilon^2 deviates in this region. Since the paper extends the claimed constraints up to M_U ~ 2 GeV using exactly these RHIC systems, the high-mass portion of the derived limit is not a reliable exclusion. The claim of constraints in the 1-2 GeV range should either be restricted to the SIS18 systems where these contributions are negligible, or the missing QGP and charm dilepton sources should be incorporated before any limit is quoted for the RHIC data.","section":"Sec. 4, Fig. 3 (left panel), discussion of RHIC region"}],"minor_comments":[{"comment":"The axis label for the Ar+KCl panel reads '3.5 AGeV', while the text and caption state 1.76 A GeV. This inconsistency should be corrected.","section":"Fig. 1 and text, Sec. 4"},{"comment":"The branching ratio expression contains a term (1 + 2 m_mu^2/M_U) that is dimensionally inconsistent; it should be (1 + 2 m_mu^2/M_U^2), as in the standard two-body decay width.","section":"Eq. (13)"},{"comment":"In the left panel of Fig. 3, both p+p at 3.5 AGeV and p+Nb at 3.5 AGeV are labeled with blue lines, which makes the legend difficult to read. Different colors or line styles would improve clarity.","section":"Figs. 1-3"},{"comment":"The text refers to 'dark matter (DM) sources' when the intended meaning is dark photon / U-boson contributions. This terminology is confusing because dark photons are mediators, not dark matter particles, and should be changed throughout.","section":"Sec. 3 and Sec. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its limitations and the PHSD model comparisons with dilepton data are valuable. However, the central quantitative claim is currently a tautology: the C_U = 0.3% value is adjusted to match BaBar14, so the resulting epsilon^2 ~ 3 x 10^-7 is not an independent constraint. Whether this is rescuable within the manuscript's scope depends on whether the authors can anchor C_U to the actual experimental uncertainties or reformulate the result as a sensitivity study for future higher-precision dilepton measurements. If such a reformulation is made, the paper would be a useful contribution to the heavy-ion dark-photon search program."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe honest one-line: this is a solid sensitivity study, not a measurement. The paper extends the group's earlier PHSD dark-photon work with additional production channels (omega Dalitz, direct rho/omega/phi decays, K+), pushes the mass range to 2 GeV, and adds RHIC energies. That part is real, new work. The PHSD baseline is checked against HADES and STAR dilepton data and looks good. The mass-dependent shape of the U-boson yield is a genuine calculation.\n\nThe problem is the headline “upper limit.” Eq. (15) defines epsilon^2 = C_U * (SM yield)/(U yield at epsilon=1). C_U is the allowed dark-photon surplus over the SM. The paper first uses C_U=10%, gets limits tens of times weaker than BaBar, then “adjusts” C_U to 0.3% so that the curve lines up with BaBar14. That makes the agreement with BaBar14 a tautology. The heavy-ion data themselves have several-percent systematic uncertainties; C_U=0.3% is not derived from them. With C_U=10%, the same PHSD calculation gives epsilon^2 ~ 1e-5 to 1e-6, which is not competitive. So the quoted 3e-7 is an assumption about required precision, not a constraint from the data.\n\nAlso, the high-mass part (M>1 GeV) at RHIC is weakened by the admitted omission of QGP and charm sources for U production. For SIS energies the shape is more reliable, but the normalization still depends on the arbitrary C_U.\n\nWhat the paper does well: the PHSD framework is appropriate, the acceptance treatment is sensible, and the work is transparent about its preliminary nature. If reframed as a forecast of the precision needed to probe epsilon^2 ~ 3e-7 with future dilepton measurements—with C_U derived from realistic systematic budgets and error bars propagated onto epsilon^2—it would be a useful paper. The channel list is a genuine step beyond the earlier papers.\n\nWho should read it: heavy-ion dilepton experimentalists planning upgrades, and phenomenologists working on dark photon searches in nuclear collisions. I would send it to review, but with the clear expectation that the authors reframe the central claim and anchor C_U to data uncertainties rather than to BaBar14.\n\nCheers,\n[Your name]","headline":"A useful sensitivity forecast for dark photons in heavy-ion dileptons, but the quoted epsilon^2 limits are set by a hand-chosen surplus parameter tuned to BaBar14, not by the data.","tokens_in":14021,"tokens_out":2359,"would_cite":true,"duration_ms":24031,"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 claims that the measured dilepton spectra in heavy-ion collisions, combined with the PHSD transport model, can constrain the dark photon kinetic mixing parameter $\\epsilon^2$ down to about $3\\times10^{-7}$ in the mass range…","keywords":["dark photons","kinetic mixing parameter","U-boson","heavy-ion collisions","PHSD transport model","dilepton spectra","HADES","STAR"],"falsifier":"A high-statistics dilepton measurement in the invariant-mass window $0.2$–$1.2$ GeV with point-to-point systematic uncertainties below $0.3\\%$ would either expose a narrow dark-photon peak at a specific $M_U$ or show that any dark-photon contribution is smaller than $0.3\\%$ of the Standard Model yield, contradicting the paper's assumed $C_U$.","tokens_in":12839,"feed_emoji":"🔭","tokens_out":14226,"duration_ms":113873,"temperature":0.7,"pith_summary":"This paper claims that dilepton spectra measured in heavy-ion collisions, from SIS to RHIC energies, can be used to place upper limits on the dark photon kinetic mixing parameter $\\epsilon^2(M_U)$. Using the PHSD transport model, which reproduces the Standard Model dilepton background, the authors extend the model to include dark photon production via Dalitz decays of pions, etas, omegas, $\\Delta$ resonances, and direct decays of vector mesons and kaons. By requiring that the dark photon contribution not exceed a specified fraction $C_U$ of the Standard Model yield in each mass bin, they extract $\\epsilon^2(M_U)$ from the existing HADES and STAR data; with $C_U = 0.3\\%$ the result matches the global upper limit near $\\epsilon^2 \\sim 3 \\times 10^{-7}$ for masses between 0.2 and 1.2 GeV. The authors note that for the highest RHIC energies the mass region above about 1 GeV is complicated by unmodeled quark-gluon plasma and charm contributions.","feed_headline":"Heavy-ion dilepton data bound dark photons to 3e-7","feed_subtitle":"PHSD model turns HADES and STAR spectra into kinetic-mixing limits near 3e-7 for masses 0.2–1.2 GeV.","key_machinery":"The central object is the dark ($U$) boson, a hypothetical massive vector particle that kinetically mixes with the Standard Model photon through the Lagrangian term $\\mathcal{L} \\sim \\frac{\\epsilon}{2} F_{\\mu\\nu} F'^{\\mu\\nu}$. The argument is carried by Eq. (15) above, which converts the assumed surplus fraction $C_U$ into a mass-dependent bound on $\\epsilon^2$ by taking the ratio of the Standard Model dilepton rate to the $U$-boson rate at $\\epsilon = 1$. The $U$-boson production rates are assembled from partial-width formulas for the Dalitz decays $\\pi^0,\\eta,\\omega \\to \\gamma U$ and $\\Delta \\to N U$, the decays $V \\to U$ for $V = \\rho,\\phi,\\omega$, the Dalitz decay $\\omega \\to \\pi^0 U$, and the kaon decay $K^+ \\to \\pi^+ U$. The mechanism that makes the procedure work is that the $U$-boson contribution scales linearly with $\\epsilon^2$ while the Standard Model background is fixed by the PHSD description, so the ratio in Eq. (15) turns a measured or assumed excess fraction into an $\\epsilon^2$ bound.","core_discovery":"The central claim is that a microscopic transport calculation, calibrated to reproduce Standard Model dilepton production in $p+p$, $p+A$, and $A+A$ collisions, can convert the measured absence of an unexplained dilepton excess into bounds on the kinetic mixing of a hypothetical $U$-boson with the photon. The bound is obtained from Eq. (15), $$\\$epsilon^{2}$(M_U) = C_U \\cdot \\frac{dN_{\\rm sum\\;SM}/dM}{dN_{\\rm sum\\;U,\\epsilon=1}/dM},$$ which equates the kinetic mixing parameter to the allowed surplus fraction $C_U$ times the ratio of the summed Standard Model dilepton yield to the summed $U$-boson dilepton yield computed at $\\epsilon = 1$. The paper identifies which production channels dominate at different masses: pion Dalitz decay below $m_{\\pi^0}$, and vector meson decays ($\\rho,\\omega,\\phi$) and $\\Delta$ resonances above about 1.2 GeV. The authors show that with $C_U = 10\\%$ the extracted curve matches the BaBar09 exclusion, and with $C_U = 0.3\\%$ it matches BaBar14 and the global compilation at $\\epsilon^2 \\sim 3\\times10^{-7}$ in the mass window $0.2$–$1.2$ GeV.","pith_inferences":["If the $C_U = 0.3\\%$ tuning reflects the true systematic precision of the data, a null search in these systems would push the allowed kinetic mixing below $3\\times10^{-7}$ in the $0.2$–$1.2$ GeV window, a region where dedicated dark photon searches currently have comparable sensitivity; heavy-ion experiments could thereby serve as an independent confirmation route.","The same ratio method could be applied to $\\mu^+\\mu^-$ pairs at higher invariant masses, or to future high-luminosity heavy-ion runs with better understood backgrounds, to extend the reach below $3\\times10^{-7}$.","Because the paper omits Drell-Yan and Bremsstrahlung production of $U$-bosons, the constraints above 1 GeV are likely conservative; including those channels could strengthen the high-mass bound."],"forward_implications":["Existing HADES and STAR dilepton measurements can already constrain dark photon kinetic mixing to the level of the current global limit near $\\epsilon^2 \\sim 3\\times10^{-7}$ for masses $0.2$–$1.2$ GeV.","The inclusion of vector meson and kaon decay channels extends the reach of heavy-ion dilepton data to dark photon masses up to about 2 GeV.","The extracted bound is proportional to the allowed surplus $C_U$, so any improvement in the precision of dilepton spectrum measurements translates directly into proportionally tighter dark photon limits.","The consistency with BaBar09 and BaBar14 suggests that heavy-ion experiments can serve as a complementary cross-check of dark photon searches at beam-dump and collider facilities."],"supporting_citations":[{"why":"Supplies the surplus-ratio method for extracting $\\epsilon^2$ from dilepton spectra that the present analysis extends to new channels and energies.","marker":"Schmidt et al., 2021"},{"why":"Provides the PHSD calculation of the Standard Model dilepton sources used as the background.","marker":"E. L. Bratkovskaya, Cassing, and Mosel (2001)"},{"why":"Contributes the HADES dilepton data and the relation between U-boson yield and virtual photon coupling.","marker":"Agakishiev et al., 2014"},{"why":"Provides the partial-width formulas for $\\pi^0$ and $\\eta$ Dalitz decays into a U-boson.","marker":"Batell, Pospelov, and Ritz (2009a)"},{"why":"Provides the $U \\to e^+e^-$ branching ratio entering the yield formula.","marker":"Batell, Pospelov, and Ritz (2009b)"},{"why":"Provides the $K^+ \\to \\pi^+ U$ partial width used at RHIC energies.","marker":"Pospelov (2009)"},{"why":"Extends the $U \\to e^+e^-$ branching ratio to masses up to 2 GeV.","marker":"Liu, Weiner, & Xue (2015)"},{"why":"Supplies the STAR Au+Au dilepton data at 200 GeV used for the high-energy comparison.","marker":"Adamczyk et al., 2015"},{"why":"Supplies the STAR Au+Au dilepton data at 19.6 GeV used for the intermediate-energy comparison.","marker":"Han (2024)"},{"why":"Provides the BaBar14 result used to calibrate $C_U = 0.3\\%$.","marker":"Lees et al., 2014"}],"fun_headline_variants":["Heavy-ion dileptons tighten dark photon mixing to 3e-7","Dark photon bounds from heavy-ion data: 3e-7 limit","PHSD model uses heavy-ion spectra to cap dark photon mixing","Dark photon kinetic mixing constrained by heavy-ion collisions","Heavy-ion data yield dark photon mixing bound near 3e-7"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole bound rests on choosing the surplus $C_U$ — the share of the dilepton yield a dark photon may add — and setting it to $0.3\\%$ because that value makes the result match BaBar14, rather than deriving it from the statistical and systematic uncertainties of the HADES and STAR data.","fun_headline_variants_meta":{"raw":{"variants":["Heavy-ion dileptons tighten dark photon mixing to 3e-7","Dark photon bounds from heavy-ion data: 3e-7 limit","PHSD model uses heavy-ion spectra to cap dark photon mixing","Dark photon kinetic mixing constrained by heavy-ion collisions","Heavy-ion data yield dark photon mixing bound near 3e-7"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000221,"raw_usage":{"total_tokens":1515,"prompt_tokens":1075,"completion_tokens":440,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":691,"completion_tokens_details":{"reasoning_tokens":361}},"tokens_in":691,"tokens_out":440,"duration_ms":4460,"temperature":1.0,"reasoning_tokens":361,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:21:57.618081+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-statistics dilepton measurement in the invariant-mass window $0.2$–$1.2$ GeV with point-to-point systematic uncertainties below $0.3\\%$ would either expose a narrow dark-photon peak at a specific $M_U$ or show that any dark-photon contribution is smaller than $0.3\\%$ of the Standard Model yield, contradicting the paper's assumed $C_U$.","supporting_citations":[{"cited_title":", Bratkovskaya, E","cited_arxiv_id":null,"evidence_quote":"Supplies the surplus-ratio method for extracting $\\epsilon^2$ from dilepton spectra that the present analysis extends to new channels and energies."},{"cited_title":", Cassing, W","cited_arxiv_id":null,"evidence_quote":"Provides the PHSD calculation of the Standard Model dilepton sources used as the background."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contributes the HADES dilepton data and the relation between U-boson yield and virtual photon coupling."},{"cited_title":"APACrefauthors \\ 2009 , Phys","cited_arxiv_id":null,"evidence_quote":"Provides the $K^+ \\to \\pi^+ U$ partial width used at RHIC energies."},{"cited_title":", Weiner, N","cited_arxiv_id":null,"evidence_quote":"Extends the $U \\to e^+e^-$ branching ratio to masses up to 2 GeV."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the STAR Au+Au dilepton data at 200 GeV used for the high-energy comparison."},{"cited_title":"APACrefauthors \\ 2024 , EPJ Web Conf","cited_arxiv_id":null,"evidence_quote":"Supplies the STAR Au+Au dilepton data at 19.6 GeV used for the intermediate-energy comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the BaBar14 result used to calibrate $C_U = 0.3\\%$."}],"review_version":1}