{"id":"a62fe7e1-e3f9-4e7d-9b4a-2238ee296da8","arxiv_id":"2608.04598","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"First-order phase transition dynamics in a chiral fluid model produce separable early pole-mass and late low-mass dilepton enhancements, strongest at low beam energies.","lead":"Simulations of heavy-ion collisions predict that a first-order QCD phase transition would leave two distinct signatures in low-mass dilepton spectra, an early bump near the vector-meson peaks and a late bump at lower masses, which survive after averaging over the collision. Why read it? The signatures give HADES, CBM, and the RHIC Beam Energy Scan a specific observable to look for when hunting for the QCD phase transition.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The COV baseline changes m_pi to 300 MeV while m_pi enters the self-energy in Eqs. (8), (13), (15), (16); the manuscript never specifies that the rate calculation is exempted, so the claimed isolation of EoS dynamics is not established.","rationale":"The paper's contribution is a controlled comparison intended to separate EoS dynamics from hadronic rescattering. The control is implemented by changing m_pi while claiming the microscopic self-energies remain identical. Equations 8, 13, 15, and 16 show that m_pi enters the rate calculation directly, so the control is not clean unless a specific undocumented choice was made. Using m_pi = 300 MeV in the rates gives the COV baseline a narrower vacuum rho width and different pion kinematics; using m_pi = 138 MeV for the rates while the EoS is generated with 300 MeV introduces a hybrid that is internally inconsistent. In either reading, the Section 7 statement that all FOPT–COV differences can be attributed directly to the EoS is not supported by the described calculation. The time-resolved decomposition itself—early pole-mass emission from the supercooled stage and late low-mass emission from reheating—remains a plausible and interesting mechanism, and the paper deserves credit for constructing the same self-energy functional for both scenarios. The concern is therefore about the comparison baseline, not about the physical mechanism. A single recomputation or an explicit statement of which pion mass enters the self-energy would resolve it, so a conditional verdict is appropriate rather than an outright rejection.","tokens_in":14492,"tokens_out":14597,"duration_ms":210687,"concrete_test":"Recompute the COV (m_pi = 300 MeV) trajectory with the self-energies in Eqs. (8), (13), (15), (16) evaluated at m_pi = 138 MeV while keeping the same T(tau), mu_B(tau) history, and compare the integrated yields in the [276,700] and [700,850] MeV windows with the version that uses m_pi = 300 MeV. If the two COV spectra differ by an amount comparable to the FOPT–COV separation in Figs. 6–8, the isolation claim fails; if they are essentially identical, the authors should state the exemption explicitly and only the proxy concern remains.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The isolation claim that FOPT–COV differences are purely EoS effects requires the two scenarios to use the same vector-meson self-energy. The COV baseline is produced by changing the pion mass (Sec. 2.1), and m_pi appears in the rate calculation: q_cm(s) in Eq. (8), threshold behavior of resonance widths in Eq. (13), the thermal pion distribution through omega_k = sqrt(k^2 + m_a^2) in Eq. (15), and the vacuum rho self-energy in Eq. (16). If m_pi = 300 MeV is used in those expressions, the COV scenario has a different vacuum width and different in-medium kinematics, so part of the FOPT–COV separation is a direct m_pi effect rather than EoS dynamics. If instead the self-energy is kept at m_pi = 138 MeV, the thermodynamic evolution and the emission rates use inconsistent pion masses. The paper does not state which choice was made, so Section 7's attribution of all differences to the EoS is unsupported. This is the reader's proxy concern, but stronger: the control variable itself enters the supposedly controlled rates.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies dilepton emission from heavy-ion collisions at sqrt(s_NN) = 2.20–6.20 GeV using non-equilibrium chiral fluid dynamics (NχFD). The authors construct a first-order phase transition scenario with the physical pion mass and a crossover scenario by raising m_pi to 300 MeV, then compute in-medium rho and omega spectral functions from resonance-driven forward scattering amplitudes with pions and nucleons. They report two temporally separated FOPT signatures: an early enhancement in the 700–850 MeV pole-mass window and a late enhancement in the 276–700 MeV low-mass continuum, both of which survive integration over the fireball history and produce a FOPT–COV separation in the excitation function, strongest at 2.20 GeV. The abstract and Section 7 conclude that these differences can be attributed directly to equation-of-state dynamics because the microscopic vector-meson self-energies are kept identical between the two scenarios.","tokens_in":14885,"tokens_out":5151,"duration_ms":59483,"significance":"The study is a useful proof-of-principle forward simulation that makes falsifiable predictions for HADES, CBM, and the RHIC Beam Energy Scan. Its strengths are that no data are fitted, the vector-meson self-energies use a fixed PDG/UrQMD resonance table without tuning to the FOPT/COV split, and the time-resolved versus integrated analysis clearly exposes the proposed mechanism. The main result—an early pole-mass enhancement versus a late low-mass enhancement—is internally consistent and nontrivial. However, the claimed attribution of all FOPT–COV differences to equation-of-state dynamics is undercut by the control-variable issue described below; the quantitative significance of the predictions depends on resolving that issue.","major_comments":[{"comment":"The central isolation claim is not established because the control variable used to create the COV scenario also enters the supposedly controlled rate calculation. The COV scenario is generated by changing m_pi from 138 to 300 MeV in Section 2.1, but m_pi appears in q_cm(s) in Eq. (8), in the resonance-width threshold behavior in Eq. (13), in the thermal pion distribution through omega_k = sqrt(k^2 + m_a^2) in Eq. (15), and in the vacuum rho self-energy in Eq. (16). The manuscript nowhere states whether the COV rate calculation uses m_pi = 300 MeV or keeps m_pi = 138 MeV. If the former, part of the FOPT–COV separation is a direct pion-mass effect on vacuum and in-medium self-energies; if the latter, the thermodynamic evolution and the emission rates use inconsistent pion masses. Either way, the repeated statements in Sections 3 and 7 that all differences are \"directly attributed\" to EoS dynamics are unsupported as written. The authors should either specify and justify the choice, or rerun the COV scenario with the rate calculation explicitly decoupled from the pion-mass variation.","section":"Sec. 2.1 and Sec. 3"},{"comment":"The phrase \"identical microscopic vector-meson self-energies\" is potentially misleading. The self-energies in Eq. (15) depend on T and mu_B through the distribution functions n_a(omega_k; T, mu_B), so the numerical self-energy values differ between FOPT and COV trajectories even if the functional form and parameters are the same. What is identical is the model, not the self-energy at a given time. The attribution of all differences to EoS dynamics should be phrased as \"identical functional form for the self-energy as a function of T and mu_B,\" and the paper should acknowledge explicitly that the different trajectories therefore produce different in-medium spectral functions. This clarification matters for the interpretation of Figures 5–8, where part of the observed separation is the expected consequence of evaluating the same self-energy formula at different thermodynamic points.","section":"Sec. 3 and Sec. 4"},{"comment":"The claim that the two signatures \"remain identifiable after integration over the full evolution\" is made on the basis of a one-dimensional Bjorken expansion with a fixed transverse radius R_T = 6.5 fm and a freeze-out criterion d^2 sigma / d tau^2 = 0. The late-time low-mass enhancement in Figure 5 is integrated over precisely the period where the 1D expansion is least physically reliable, and the freeze-out condition is itself EoS-dependent. I therefore ask the authors to provide a quantitative check that the integrated FOPT–COV separation is not controlled by the chosen freeze-out definition or by the fixed R_T. Since the paper already describes itself as a proof of principle, this check would substantially strengthen the claim that the separation is robust.","section":"Sec. 7"}],"minor_comments":[{"comment":"The caption states \"The squared speed of sound c_s^2(tau) (top) and proper time evolution of the order parameter sigma(tau) (bottom),\" but the text in Section 2.2 and the panel layout indicate that sigma(tau) is the top panel and c_s^2(tau) the bottom panel. Please reverse the order in the caption.","section":"Fig. 2 caption"},{"comment":"References 4 and 5 appear to be the same citation (Borsanyi et al., JHEP 09, 073 (2010)); please check and remove the duplicate.","section":"References"},{"comment":"In the low-mass window the FOPT–COV separation is visible only at 2.20 GeV after integration, whereas the pole-mass window shows separation over 2.20–3.50 GeV. The discussion in Section 6 should state this asymmetry more prominently, since it qualifies the general statement that both signatures survive integration.","section":"Sec. 6 and Fig. 8"},{"comment":"The choice of the lower edge of the low-mass window at 276 MeV is not motivated. Please state whether this value is related to the rho-to-dilepton threshold, an experimental acceptance cut, or another criterion, since the time-profile and excitation-function results depend on the window boundaries.","section":"Sec. 4"}],"recommendation":"major_revision","confidential_remarks":"This is a competent model study from an established group, and the internal logic of the simulation is consistent. My recommendation of major revision is driven by one crucial ambiguity: the parameter used to generate the COV scenario is the same pion mass that enters the vector-meson self-energies, so the claimed isolation of EoS dynamics is not yet demonstrated. If the authors can show that the rate calculation is insensitive to this choice, or explicitly decouple the two uses of m_pi, the paper would be suitable for publication. I do not see a deeper circularity: the temporal ordering of the enhancements is a genuine forward-simulation output. The paper fits the journal scope and would be of interest to the heavy-ion dilepton community, but the central attribution claim must be fixed first."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the main claim is new and the simulation logic is clean. The paper shows, within its model, that a first-order transition gives an early pole-mass enhancement and a late low-mass enhancement, and that both survive integration. The full UrQMD/PDG resonance set for the self-energies is a real extension over Eletsky et al., and the excitation function benchmark is a useful experimental target. The temporal decomposition is the genuinely new part.\n\nThe soft spot that matters most: Section 7 attributes all FOPT–COV differences to EoS dynamics, but the COV scenario is generated by changing m_pi to 300 MeV, and m_pi appears in the self-energy: q_cm in Eq. (8), resonance widths in Eq. (13), thermal pion distribution in Eq. (15), and the vacuum rho width in Eq. (16). The paper never states whether m_pi=300 MeV was used in those expressions. If it was, part of the separation is a direct pion-mass effect on the rates; if it wasn't, the thermodynamic evolution and rate calculations use inconsistent pion masses. Either way, the claim that the microscopic rates are identical across scenarios is not established on the text as written. This is a load-bearing ambiguity, not a nitpick. The authors should either use the physical m_pi in all rate formulas and say so, or perform a check with m_pi held fixed.\n\nOther soft spots are minor. The 1D Bjorken geometry with fixed R_T is a known simplification and the authors call it a proof of principle. The artificial crossover is a modeling choice; it is more defensible if the rate kinematics are held fixed. No code or data are shipped, so the quantitative predictions are not independently reproducible.\n\nWhat is solid: the internal logic is consistent, the time profiles support the early/late interpretation, and the excitation function trend is a clear, falsifiable benchmark. The paper is worth a serious referee, but the m_pi ambiguity should be fixed before publication. If you handle it, ask the authors to state explicitly which pion mass enters Eqs. (8)–(16) in each scenario, and ideally rerun with m_pi fixed.","headline":"Solid phenomenological study with a clean temporal separation argument, but the isolation claim is undermined by an unstated choice of pion mass in the self-energy formulas.","tokens_in":15278,"tokens_out":2172,"would_cite":false,"duration_ms":24170,"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":"A first-order QCD phase transition leaves two temporally separate dilepton signatures, an early pole-mass enhancement and a late low-mass one, that survive integration over the full fireball evolution.","keywords":["heavy-ion collisions","QCD phase transition","dileptons","equation of state","vector meson spectral functions","non-equilibrium chiral fluid dynamics","first-order phase transition","low beam energies"],"falsifier":"Measure the dilepton yield in the pole-mass window $M\\in[700,850]$ MeV as a function of beam energy at $\\sqrt{s_{NN}}\\approx2.2$–$3.5$ GeV: if the measured excitation function tracks the crossover baseline rather than showing the predicted growing excess toward low energies, the paper's central separation claim is falsified.","tokens_in":14259,"feed_emoji":"⚛️","tokens_out":7635,"duration_ms":83418,"temperature":0.7,"pith_summary":"The paper aims to show that lepton pairs emitted from heavy-ion collisions can distinguish a first-order QCD phase transition from a smooth crossover, at low beam energies, even after all time information is integrated away. The authors build both scenarios inside the same non-equilibrium chiral fluid dynamics model, with identical microscopic in-medium vector-meson self-energies, so any spectral difference is traced to the equation of state. Their central result is that a first-order transition imprints two temporally separated signatures: an early enhancement near the $\\rho$ and $\\omega$ pole masses, signaling the onset of equation-of-state softening, and a later enhancement in the low-mass continuum, produced by reheating and a longer-lived fireball. Both pieces survive the full evolution integral, and the pole-mass window is the cleanest discriminator, with the largest separation at $\\sqrt{s_{NN}}=2.20~\\mathrm{GeV}$. A sympathetic reader would care because this gives a concrete, testable benchmark for upcoming low-energy dilepton measurements.","feed_headline":"Two dilepton signatures can tell a phase transition from a crossover","feed_subtitle":"Early pole-mass and late low-mass dilepton enhancements survive the full fireball, separating first-order from crossover scenarios.","key_machinery":"The load-bearing object is the non-equilibrium chiral fluid dynamics evolution, in which the chiral condensate $\\sigma$ obeys a Langevin equation coupled back into the fluid via energy-momentum exchange, producing supercooling, reheating, and delayed freeze-out when the trajectory crosses a first-order boundary. The second ingredient is the in-medium vector-meson propagator: the retarded self-energy $\\Sigma_V$ is built from the forward scattering amplitude $f^{Va}_{\\rm c.m.}$ summed over a broad resonance set, and the thermal dilepton rate is obtained from the imaginary part of the electromagnetic current correlator $\\mathrm{Im}\\,\\Pi^R_{\\rm em}$ evaluated along the non-equilibrium trajectory. What this machinery does is attach a physically identical hadronic in-medium mechanism to two different phase structures, so that any difference between the first-order and crossover scenarios in the spectra is attributable to the equation of state alone.","core_discovery":"Within a quark-meson effective theory with non-equilibrium dynamics, the chiral order parameter is propagated by a Langevin equation while exchanging energy and momentum with an expanding quark fluid. Shifting the explicit symmetry-breaking term through the pion mass produces either a first-order transition (physical $m_\\pi=138$ MeV) or a crossover ($m_\\pi=300$ MeV) while leaving the hadronic rescattering content unchanged. The in-medium $\\rho$ and $\\omega$ spectral functions are computed from resonance-driven forward scattering amplitudes with pions and nucleons, so the two scenarios differ only through the bulk evolution. In the first-order case the trajectory overshoots into a mechanically unstable (spinodal) region, supercools, reheats, and freezes out later; the paper shows that this produces an early emission excess in the pole-mass window $M\\in[700,850]$ MeV and a late emission excess in the low-mass window $M\\in[276,700]$ MeV. After integrating the emission rate over the whole fireball history, both excesses remain visible in the invariant-mass spectra and in the beam-energy excitation function, establishing that equation-of-state dynamics can be separated from hadronic in-medium broadening in dilepton observables.","pith_inferences":["A ratio of low-mass to pole-mass yields as a function of beam energy would likely amplify the first-order-versus-crossover separation beyond the two individual windows, since the two signatures have opposite time order; the paper does not compute this ratio, but the time-resolved results imply it.","The artificial crossover via $m_\\pi=300$ MeV is the main interpretive risk; if the real QCD crossover in this baryon-density range is steeper or located elsewhere, the quantitative size and beam-energy range of the separation could shift, even though the qualitative two-stage mechanism should persist.","The same two-window, two-time scheme could be applied to other first-order-transition candidates where a softening equation of state and a reheating stage are expected, turning the time-resolved signatures into a generic diagnostic rather than a QCD-specific prediction.","Because the low-mass enhancement is a chronometer of reheating, combining dilepton spectra with photon or interferometry measurements at the same beam energies could measure the duration of phase conversion directly."],"forward_implications":["In the $M\\in[700,850]$ MeV window, the first-order-transition enhancement survives integration over the full evolution for all beam energies in $\\sqrt{s_{NN}}=2.20$–$3.50$ GeV, making the pole-mass excitation function the strongest discriminator between a first-order transition and a crossover.","In the low-mass window $M\\in[276,700]$ MeV, the first-order-transition enhancement after integration appears mainly at the lowest beam energy, because the late reheating emission accumulates over a long period.","If confirmed, the predicted low-energy enhancement would indicate that the collision trajectories cross or approach a first-order phase boundary, constraining the QCD phase diagram at high baryon density.","Combining the dilepton enhancement with two-particle interferometry radii, which also encode the space-time evolution of the fireball, could provide a cross-check of the delayed freeze-out predicted by the first-order scenario."],"supporting_citations":[{"why":"supplies the forward-scattering formalism from which the in-medium vector-meson self-energies are constructed.","marker":"[38]"},{"why":"provides the non-equilibrium chiral fluid dynamics model coupling the sigma field to the expanding fluid.","marker":"[35]"},{"why":"defines the shock-adiabat initial conditions that set each beam-energy trajectory.","marker":"[43]"},{"why":"gives the thermal dilepton emission rate through the imaginary part of the electromagnetic current correlator.","marker":"[44]"},{"why":"provides the comprehensive hadronic resonance table that fixes the scattering amplitudes and branching ratios.","marker":"[39, 40]"}],"fun_headline_variants":["Dilepton timing separates phase transition from crossover","Early and late dilepton excesses flag first-order transition","Two dilepton windows reveal QCD phase transition signature","Low-mass dileptons distinguish EoS softening from rescattering","Phase transition leaves twin dilepton fingerprints in spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The crossover is simulated by raising the pion mass to 300 MeV, and the paper's attribution of all first-order-versus-crossover spectral differences to equation-of-state dynamics assumes this artificial change is a faithful stand-in for the real QCD crossover and does not secretly alter the self-energies or the meaning of the thermodynamic trajectory.","fun_headline_variants_meta":{"raw":{"variants":["Dilepton timing separates phase transition from crossover","Early and late dilepton excesses flag first-order transition","Two dilepton windows reveal QCD phase transition signature","Low-mass dileptons distinguish EoS softening from rescattering","Phase transition leaves twin dilepton fingerprints in spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00074,"raw_usage":{"total_tokens":3304,"prompt_tokens":945,"completion_tokens":2359,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":561,"completion_tokens_details":{"reasoning_tokens":2282}},"tokens_in":561,"tokens_out":2359,"duration_ms":15916,"temperature":1.0,"reasoning_tokens":2282,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T21:01:44.281003+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the dilepton yield in the pole-mass window $M\\in[700,850]$ MeV as a function of beam energy at $\\sqrt{s_{NN}}\\approx2.2$–$3.5$ GeV: if the measured excitation function tracks the crossover baseline rather than showing the predicted growing excess toward low energies, the paper's central separation claim is falsified.","supporting_citations":[{"cited_title":"Eletsky, M","cited_arxiv_id":null,"evidence_quote":"supplies the forward-scattering formalism from which the in-medium vector-meson self-energies are constructed."},{"cited_title":"Nahrgang, S","cited_arxiv_id":null,"evidence_quote":"provides the non-equilibrium chiral fluid dynamics model coupling the sigma field to the expanding fluid."},{"cited_title":"Bumnedpan, J","cited_arxiv_id":null,"evidence_quote":"defines the shock-adiabat initial conditions that set each beam-energy trajectory."},{"cited_title":"Rapp, Adv","cited_arxiv_id":null,"evidence_quote":"gives the thermal dilepton emission rate through the imaginary part of the electromagnetic current correlator."}],"review_version":2}