{"id":"6dcb5526-794a-4c64-9cb1-a90507d5395a","arxiv_id":"2501.03893","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Regeneration keeps K*/K falling linearly with time, giving hadronic-stage lifetimes 2 to 4 times longer than the exponential-decay estimate used by STAR.","lead":"This paper uses a set of coupled rate equations to show that the K*/K ratio in heavy-ion collisions falls off linearly, not exponentially, because resonance decays are counterbalanced by regeneration. This changes the inferred duration of the hadronic stage, making it two to four times longer than earlier STAR estimates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section II's rate equations omit the advertised collisional loss L_coll; adding K* dissociation could steepen the slope in Eq. (12) and shrink the claimed 2-4x lifetimes.","rationale":"The reader's weakest assumption concerns the universality of the slope across centralities and energies. That is a real and important gap, but I see a more fundamental problem: the slope is computed from a network whose only K*-changing process is h* ↔ h + X (Eq. 4), even though the introduction explicitly advertises 'loss rates' L_coll in Eq. (2). This is an internal inconsistency in the manuscript. At BES energies, where baryon densities are high, K* collisional dissociation is a standard inelastic channel in hadronic transport and could contribute comparably to the vacuum decay width. Adding it would steepen the slope and lower the equilibrium K* abundance, acting on both the numerator and denominator of Eq. (12). The benchmark against STAR data (Fig. 1) is not a substitute for including these channels, because the kinetic freeze-out time is inferred from the truncated evolution itself; the agreement only demonstrates that the model can find a time at which the data is reproduced. The proposed test—adding Γ_coll and re-extracting the slope—would settle whether the central factor-2-4 claim survives. Because the omission is concrete and testable, the paper should not be accepted as a final quantitative result without either implementing the missing loss terms or showing that they are numerically negligible. This supports the reader's CONDITIONAL verdict rather than changing it.","tokens_in":8632,"tokens_out":12795,"duration_ms":136833,"concrete_test":"Extend the Section II network by adding an explicit collisional loss term -Γ_coll N_K* to Eq. (4), with Γ_coll estimated from K*+N and K*+π inelastic cross sections taken from transport models or measured dissociation data. Recompute d(K*/K)/dt and Δt_hadronic from Eq. (12) for central Au+Au at √sNN = 7.7 and 39 GeV. If the slope steepens by ≳30% (e.g., from -0.02 fm^-1 to -0.03 fm^-1 or beyond), the inferred lifetimes shrink proportionally and the factor-2-4 improvement over STAR must be revised; also rerun with K*+baryon channels included to bracket the effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central relation Eq. (12) rests entirely on the numerically extracted slope d(K*/K)/dt ≈ -0.02 fm^-1. But the equations actually solved in Section II, Eq. (4) and the 23-species network, contain only the decay/regeneration channel h* ⇌ X + h. The 'loss rate due to collisions' L_coll announced in Eq. (2) is never implemented. In the dense hadronic medium at RHIC-BES energies, K* can also be dissociated by inelastic collisions such as K* + N → K + N or K* + π → K_1 → ...; these channels add extra loss and lower the equilibrium K* abundance, making the net slope more negative than -0.02 fm^-1. Since Eq. (12) divides the measured suppression by this slope, omission of collisional losses makes the inferred Δt_hadronic too large. The quantitative claim that lifetimes are 2-4 times the STAR estimates therefore depends on an unverified completeness assumption. The good reproduction of STAR K*/K in Fig. 1 does not resolve this, because the model's kinetic freeze-out time is itself read off from the same incomplete evolution; matching the data at some time only shows the truncated network can reach the measured ratio, not that the missing channels are negligible.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the time evolution of the K*/K ratio during the hadronic phase in Au+Au collisions at RHIC-BES energies using a network of 23 coupled rate equations with partial chemical equilibrium and detailed balance. The authors report that the K*/K ratio decreases approximately linearly rather than exponentially, because regeneration counteracts decay, and they extract an essentially energy-independent slope d(K*/K)/dt ≈ -0.02 fm^-1 (Section IV, Fig. 3). They propose an improved relation, Eq. (12), between the K*/K suppression and the hadronic-phase lifetime, and apply it to STAR data to obtain lifetimes that are 2-4 times larger than the STAR estimates based on the exponential formula in Eq. (1).","tokens_in":8948,"tokens_out":4179,"duration_ms":43126,"significance":"If the central claims hold, the paper offers a more realistic method for extracting the hadronic-stage duration from resonance suppression and identifies a systematic bias in previous exponential-decay estimates. The explicit rate-equation framework with detailed balance, the reproduction of centrality and energy dependence of the STAR K*/K data in Fig. 1, and the clear linear behavior shown in Fig. 3 are genuine strengths. The proposed relation Eq. (12) is simple and directly usable by STAR, NA61, and future FAIR experiments. However, the quantitative validity of the result depends on the completeness of the loss terms and on the universality of the slope, both of which need additional support before the factor-2-4 statement can be considered robust.","major_comments":[{"comment":"The loss rate due to collisions, L_coll, announced in Eq. (2) is never implemented in the equations actually solved. The evolution equation Eq. (4) contains only the decay/regeneration channel h* ⇌ X + h, and no collisional dissociation term is defined or discussed in Section II. In the dense hadronic medium at RHIC-BES energies, channels such as K* + N → K + N or K* + π → K_1 → ... would add an extra loss mechanism that makes d(K*/K)/dt more negative than -0.02 fm^-1. Since Eq. (12) divides the measured suppression by this slope, omitting collisional losses could make the inferred Δt_hadronic systematically too large, and the claimed factor of 2-4 relative to STAR would rest on an unverified completeness assumption. The authors should either implement L_coll in the network or provide a quantitative estimate of its magnitude and demonstrate that it is negligible for the slope.","section":"Section II, Eqs. (2) and (4)"},{"comment":"The slope s = d(K*/K)/dt ≈ -0.02 fm^-1 is extracted from central Au+Au collisions in Fig. 3, but Eq. (12) is then applied to all centralities and energies in Fig. 4 without demonstrating that the slope is centrality independent. The slope also depends on the volume parametrization Vch in Eq. (9), the chemical freeze-out time tch taken from UrQMD, and the transverse expansion time t_perp in Eq. (11); none of these dependencies is quantified. The paper provides no uncertainty on s and no sensitivity analysis, although any change in s linearly rescales all extracted lifetimes. The statement in Section II that t_perp has only minor influence is not supported by a calculation. Please show the slope as a function of centrality, provide an uncertainty band for s, and propagate that band into Fig. 4.","section":"Section IV, Eq. (12), Figs. 3 and 4"},{"comment":"The benchmarking against STAR data in Fig. 1 does not independently validate the slope used in Eq. (12). The kinetic freeze-out point in Fig. 3 is read off when the model's K*/K matches the same STAR data that enter the numerator of Eq. (12), and the slope is taken from the same model run. This consistency is necessary but not sufficient to certify the slope as a universal constant. An independent cross-check, such as comparing the model's K*/K evolution against a transport-code calculation with controlled resonance processes, would substantially strengthen the claim that the slope is a reliable input for lifetime extraction.","section":"Section III, Fig. 1"}],"minor_comments":[{"comment":"The caption states that filled symbols with bands denote the rate-equation calculations, but the meaning of the band is not defined in the text; please specify whether it reflects numerical uncertainty, centrality-bin width, or a model variation.","section":"Fig. 1 caption"},{"comment":"The vertical axis extends to negative values and some STAR points appear at negative hadronic lifetimes, but the text does not explain how a negative duration is to be interpreted; please address this for clarity.","section":"Fig. 4"},{"comment":"The paper refers to a network of 23 rate equations and to 'the specific rate equations described in [24]', but the included species and the relevant decay/regeneration channels are not itemized; a table listing the 23 species and their principal channels would improve reproducibility.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The comparison to 'other estimates' in Section IV cites refs. [10,42,43], several of which share authors with the present paper; this is not disqualifying, but the editor may wish to ensure that the agreement is not presented as an independent cross-check. The missing L_coll term in Eq. (4) is the main technical issue and should be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper makes a real correction to the standard resonance-suppression estimate. STAR's exponential-decay formula misses regeneration; with the rate-equation network the K*/K ratio falls approximately linearly and the extracted hadronic stage is 2-4 times longer, bringing it into line with HBT/transport estimates. The qualitative point survives.\n\nThe new result is not the rate-equation framework itself — that is from the authors' 2022 paper — but the application to RHIC-BES energies, the linear-decrease finding, and the improved lifetime relation Eq. (12). Fig. 1's reproduction of the STAR centrality and energy dependence is a genuine benchmark, and the authors are clear that regeneration counteracts decay.\n\nSoft spots. First, the stress-test note is correct: Eq. (2) advertises a collision loss term L_coll, but the equations actually solved, Eq. (4) plus the 23-species network, contain only the decay/regeneration channel h* ⇌ X + h. K* + N or K* + π inelastic dissociation is not in the network. Missing loss would make the slope more negative than -0.02 fm^-1, so the claimed factor 2-4 is an upper bound, not a robust number. This is a real omission.\n\nSecond, the universal slope is asserted with no uncertainty and no centrality check. It is read off central collisions in Fig. 3 and applied to all centralities via Eq. (12), even though the volume expansion, chemical freeze-out time, and V_ch parametrization all enter. The circularity concern also has bite: the model is benchmarked against STAR K*/K and then used to extract the lifetimes from the same data. Not fatal, but it means the 2-4x number is not an independent measurement.\n\nThird, the 'now in line with other estimates' comparison cites several papers with overlapping authors; not disqualifying, but the agreement is less independent than it looks. And no code or data are released, so sensitivity analysis is impossible.\n\nWho is this for? Practitioners analyzing K*/K at RHIC-BES, NA61, or FAIR. They should read it, and they should not take 0.02 fm^-1 as a universal constant until the collisional loss and centrality dependence are quantified. I would send it to a serious referee: the linear-regeneration point is important enough to put on record, and the missing L_coll needs to be either implemented or explicitly argued away. I would not desk-reject it.","headline":"Regeneration makes K*/K decay roughly linear and implies a 2-4x longer hadronic stage, but Eq. (12)'s universal slope is asserted, not demonstrated, and the advertised collision-loss term is absent from the solved network.","tokens_in":9469,"tokens_out":4190,"would_cite":true,"duration_ms":35180,"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 hadronic stage in heavy-ion collisions lasts 2-4 times longer than STAR's estimate because resonance regeneration makes the K*/K ratio fall linearly rather than exponentially.","keywords":["hadronic stage duration","resonance suppression","K*/K ratio","rate equations","partial chemical equilibrium","heavy-ion collisions","RHIC beam energy scan","regeneration"],"falsifier":"Measure the $K^*/K$ ratio in a centrality or system where the hadron phase duration is independently known (e.g., from two-particle HBT radii) and check whether Eq. (12) with the universal slope reproduces that duration; a systematic mismatch beyond the stated uncertainties would show that the slope is not universal.","tokens_in":8449,"feed_emoji":"⏱️","tokens_out":6508,"duration_ms":53246,"temperature":0.7,"pith_summary":"This paper argues that the duration of the hadronic stage in central Au+Au collisions at RHIC-BES energies has been substantially underestimated. By running a full set of coupled rate equations that include both resonance decay and regeneration, the authors find that the $K^*/K$ ratio falls essentially linearly in time, not exponentially as previously assumed. The slope is nearly energy-independent, $d(K^*/K)/dt \\approx -0.02~\\mathrm{fm}^{-1}$, giving a direct formula for the hadron phase lifetime from the ratio difference between chemical and kinetic freeze-out. The resulting lifetimes are 2-4 times larger than STAR's estimates and fall in line with independent methods.","feed_headline":"Hadron phase lives 2-4 times longer than STAR's estimate","feed_subtitle":"Regeneration makes K*/K fall linearly, not exponentially, giving a universal clock for the hadronic stage.","key_machinery":"The machinery is a network of 23 coupled rate equations for hadron species in partial chemical equilibrium, where each resonance's multiplicity evolves under detailed balance between its decay $h^* \\to h + X$ and its regeneration $h + X \\to h^*$. Equilibrium multiplicities are fixed at chemical freeze-out using a standard parametrization of $T_{\\rm ch}$ and $\\mu_{B,\\rm ch}$, and the fireball expansion enters through a parametrized volume $V(t)$ with conserved entropy, baryon number, and strangeness. The load-bearing output is the nearly energy-independent slope of the $K^*/K$ ratio versus time, which converts a measured ratio difference into a hadron phase lifetime through Eq. (12).","core_discovery":"The central finding is that the time evolution of the $K^*/K$ ratio after chemical freeze-out is controlled by a counterbalance between resonance decay and regeneration, and this balance makes the ratio decrease approximately linearly at a slope that is essentially independent of collision energy. Because regeneration compensates for decays, the ratio does not follow the simple exponential decay law used in earlier estimates. The difference between the $K^*/K$ ratio at chemical and kinetic freeze-out, divided by the universal slope of about $0.02~\\mathrm{fm}^{-1}$, yields the lifetime of the hadronic stage, which the authors find to be 2-4 times larger than the STAR collaboration's previous extraction.","pith_inferences":["The same detailed-balance mechanism should apply to other resonance ratios such as $\\rho/\\pi$ or $\\Delta/N$; if those ratios also evolve linearly, the universal-slope method could be extended to multi-resonance consistency checks.","One could test the universality claim directly by comparing central and peripheral collisions: if the slope changes with system size or centrality, Eq. (12) would need a centrality-dependent slope.","Because the slope is anchored to the assumed volume expansion law and the chemical freeze-out time from UrQMD, a direct measurement of the time dependence of $K^*/K$ (for instance through femtoscopy) could validate the linear evolution independently."],"forward_implications":["The hadronic phase at RHIC-BES energies is roughly 2-4 times longer than the STAR collaboration's estimate, resolving the tension between resonance suppression and other lifetime probes.","Analyses that treat resonance suppression as a pure exponential decay will systematically underestimate hadron phase durations, so the linear formula provides a corrected tool for future measurements.","The energy independence of the slope means a single constant can be used across the beam-energy scan, simplifying extraction of lifetimes at NA61, STAR, and FAIR.","The agreement with independent 4-8 fm/c estimates suggests that the rate-equation network captures the dominant physics of the late hadronic stage, at least for the $K^*$ channel."],"supporting_citations":[{"why":"Supplies the partial-chemical-equilibrium rate equation framework used for the resonance time evolution.","marker":"[24]"},{"why":"Provides the measured $K^*/K$ ratios and the STAR estimates that the new lifetimes are compared against.","marker":"[26]"},{"why":"Gives the parametrization of $T_{\\rm ch}$ and $\\mu_{B,\\rm ch}$ used to set initial conditions.","marker":"[38]"},{"why":"Supplies the time-dependent volume parametrization used in the rate equations.","marker":"[40]"},{"why":"Fixes the chemical freeze-out time $t_{\\rm ch}$ from UrQMD simulations.","marker":"[41]"},{"why":"Provides independent previous estimates of the hadron phase lifetime (4-8 fm/c) with which the new results agree.","marker":"[10]"},{"why":"Converts HBT radii to an emission duration, used as an independent lifetime estimate.","marker":"[44]"}],"fun_headline_variants":["Hadronic phase lasts 2-4 times longer than thought","K* regeneration revises hadron lifetime upward","Hadron phase clock: K*/K slope gives new lifetimes","Resonance regeneration extends hadronic stage","STAR's hadron lifetime underestimated by factor 2-4"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The slope of about $-0.02~\\mathrm{fm}^{-1}$ is extracted for central Au+Au collisions and then used as a universal constant for all centralities and energies, with no demonstration that it is centrality-independent and no estimate of its uncertainty.","fun_headline_variants_meta":{"raw":{"variants":["Hadronic phase lasts 2-4 times longer than thought","K* regeneration revises hadron lifetime upward","Hadron phase clock: K*/K slope gives new lifetimes","Resonance regeneration extends hadronic stage","STAR's hadron lifetime underestimated by factor 2-4"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000552,"raw_usage":{"total_tokens":2575,"prompt_tokens":831,"completion_tokens":1744,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":447,"completion_tokens_details":{"reasoning_tokens":1664}},"tokens_in":447,"tokens_out":1744,"duration_ms":11754,"temperature":1.0,"reasoning_tokens":1664,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:44:36.687943+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $K^*/K$ ratio in a centrality or system where the hadron phase duration is independently known (e.g., from two-particle HBT radii) and check whether Eq. (12) with the universal slope reproduces that duration; a systematic mismatch beyond the stated uncertainties would show that the slope is not universal.","supporting_citations":[{"cited_title":"Status of Chemical Freeze-Out","cited_arxiv_id":"hep-ph/0607164","evidence_quote":"Gives the parametrization of $T_{\\rm ch}$ and $\\mu_{B,\\rm ch}$ used to set initial conditions."},{"cited_title":"Pan and S","cited_arxiv_id":null,"evidence_quote":"Supplies the time-dependent volume parametrization used in the rate equations."},{"cited_title":"Probing chemical freeze-out criteria in relativistic nuclear collisions with coarse grained transport simulations","cited_arxiv_id":"2007.06440","evidence_quote":"Fixes the chemical freeze-out time $t_{\\rm ch}$ from UrQMD simulations."},{"cited_title":"Timescales in heavy ion collisions","cited_arxiv_id":"1607.06188","evidence_quote":"Converts HBT radii to an emission duration, used as an independent lifetime estimate."}],"review_version":1}