{"id":"d92b44ca-748d-4576-9a3d-b40b882216bf","arxiv_id":"2608.06731","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Solving the in-medium T-matrix in the laboratory frame rather than the center-of-mass frame increases the charm quark drag coefficient by 25-40% at low momenta in a quark-gluon plasma model.","lead":"This paper calculates heavy quark scattering in the quark-gluon plasma directly in the medium's rest frame instead of the usual center-of-mass frame, and finds that the two frames give scattering rates that differ substantially at low momentum. If correct, heavy quark drag coefficients, used to read off plasma transport properties from experiment, could be revised upward by 25-40%.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CM-frame comparison omits the energy mapping: T_CM must be evaluated at the invariant mass sqrt(E^2-P^2), not the lab-frame E, or the 25-40% drag correction may be an artifact of an energy shift.","rationale":"The reader's weakest assumption, the constant width Gamma, affects the absolute scale of A(p) but the relative lab-vs-CM correction is nearly the same for Gamma=100 and 200 MeV in Fig. 2, so it is unlikely to overturn the central comparison. The CM-energy mapping, by contrast, is a binary choice that determines whether the claimed discrepancy is a genuine frame effect or an energy shift. The paper's own text in Section III is silent on this point: it gives the boost for momenta but not the energy argument of Eq. (6). Because Eq. (14) is linear in the squared amplitude and the phase space favors moderate P, a systematic error of about 0.7 GeV in the baseline energy can plausibly produce corrections of the reported size. The P=0 check only validates the code in the rest frame and does not constrain the mapping at P>0. A clarifying calculation with both energy assignments would settle the matter. This is consistent with the reader's conditional verdict; the condition should explicitly include the energy mapping.","tokens_in":9395,"tokens_out":13338,"duration_ms":124786,"concrete_test":"Rerun the CM-frame curves in Figs. 1-2 with T_CM explicitly evaluated at E_CM = sqrt(E^2-P^2) for each lab event, and as a control at E_CM = E. Add the mapping formula to Section III. If the 25-40% drag enhancement is unchanged under the invariant-mass mapping, the concern is resolved; if the ratio moves toward 1 or changes sign, the headline correction is an energy-mapping artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III defines the lab-frame on-shell amplitudes via Eq. (12) and then Lorentz-transforms the lab momenta to obtain q_CM and theta_CM, but it never states the energy argument at which the CM T-matrix of Eq. (6) is evaluated. For T_Lab(E,P,...) and T_CM to describe the same physical event, the CM calculation must use the total energy in the pair rest frame, E_CM = sqrt(E^2-P^2). The text specifies only beta=P/E and gamma for the momentum boost; no equivalent statement is made for the energy. With E fixed at 3.05 GeV in Fig. 1 and P up to 2 GeV, E_CM falls to about 2.3 GeV, a substantial shift. If the CM curves are instead evaluated at the lab-frame E, the baseline entering Eq. (14) describes scattering at too high an invariant mass as P grows, and the reported 25-40% enhancement of the lab-frame drag coefficient would partly or wholly reflect this energy offset rather than the frame dependence of the in-medium T-matrix. This is load-bearing because the central quantitative claim is the comparison against the CM baseline.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a laboratory-frame in-medium two-body T-matrix formalism. Since the medium rest frame picks a preferred direction, and a finite total pair momentum P breaks SO(3) down to SO(2), the authors decompose the Lippmann-Schwinger equation in azimuthal harmonics about the P-axis, obtaining uncoupled equations for the components T_m (Eq. 9). They apply the framework to charm-heavy-light antiquark scattering in the quark-gluon plasma at T=190 and 258 MeV, using a screened Cornell potential and constant quark widths Gamma=100-200 MeV, and compare the resulting amplitudes and drag coefficients with conventional CM-frame calculations. The central numerical claim is that the laboratory-frame drag coefficient exceeds the CM-frame baseline by 25-40% at low momenta (Fig. 2), while the temperature dependence remains similar.","tokens_in":9650,"tokens_out":9983,"duration_ms":84986,"significance":"If the central comparison is correct, this is an important and overdue assessment of a systematic uncertainty in heavy-quark transport: essentially all previous in-medium T-matrix applications were solved in the CM frame, and the frame mismatch has not previously been quantified. The formal decomposition is clean, and the P=0 limit correctly recovers the CM result, providing a nontrivial consistency check. The drag comparison is self-contained rather than a fit: both amplitudes use the same potential, masses, and widths, so the claimed enhancement is an output of the calculation. However, because the energy argument used in the CM baseline is not specified (Section III) and the numerical implementation is not described, the quantitative claim is not yet reproducible.","major_comments":[{"comment":"The CM-frame T-matrix of Eq. (6) has an energy argument E_CM, but the text never states at which energy it is evaluated when comparing to T_Lab(E,P). The momenta are Lorentz-transformed using beta=P/E and gamma, but the total energy in the pair rest frame is a separate quantity; under the relativistic kinematics of Eq. (12) it is E_CM = sqrt(E^2-P^2). With E=3.05 GeV and P up to 2 GeV, E_CM decreases to about 2.3 GeV at the largest P. If the CM amplitude was evaluated at E rather than E_CM, the baseline entering Eq. (14) describes scattering at too high an invariant mass as P grows, and the reported 25-40% enhancement of the drag coefficient could be partly or wholly an energy-shift artifact rather than a frame effect. The authors must state the energy argument, correct the CM calculation if needed, and recompute the drag comparison.","section":"Section III, Eqs. (6) and (12), Fig. 1"},{"comment":"The numerical solution of Eq. (9) is not described. No information is given on the discretization of (q, x''), the momentum cutoff, the number of azimuthal components m retained in Eq. (10), or convergence tests against these parameters. The P=0 check is reassuring but does not establish convergence at P>0, where the m>0 components are essential. Since the central quantitative claim is a 25-40% correction, the authors should provide the numerical scheme and convergence checks.","section":"Section III, Eqs. (9) and (10), Figs. 1 and 2"},{"comment":"The reduction of Eq. (14) to an integral over laboratory-frame variables is described only in words. The text does not specify how the 4-momentum delta functions are used, how the on-shell amplitude is evaluated off the specific configurations shown in Fig. 1, or how the width parameters Gamma_Q and Gamma_q enter the spectral functions used in the drag kernel. Without this mapping the central drag result cannot be reproduced or checked. Please present the reduced phase-space integrals and any interpolations used.","section":"Section IV, Eq. (14), Fig. 2"}],"minor_comments":[{"comment":"The text states p_Q=p_q=1 GeV, but the caption of Fig. 1 quotes E=3.05 GeV; with the stated masses Eq. (12) gives E approximately 2.72 GeV, so one of these numbers should be corrected.","section":"Section III, Fig. 1 caption"},{"comment":"The parameters alpha and sigma of the screened Cornell potential are never given numerically; please list them or state explicitly that they are identical to a specific prior work.","section":"Section III"},{"comment":"The relation V_m=V_{-m} should be spelled out with the phase convention of Eq. (7), since the expansion uses complex exponentials e^{imPhi} while the amplitudes are real.","section":"Section II, Eq. (10)"},{"comment":"There are typos in Section IV ('heayv-light') and Section V ('non-perturabtive'); please proofread the manuscript.","section":"Sections IV and V"},{"comment":"The qualitative statement that the CM approximation is most reliable in the forward region is made on the basis of |T|^2, whereas Fig. 1(b) shows that |M|^2 behaves differently; please clarify which quantity supports each qualitative conclusion.","section":"Section III, Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The energy-mapping issue (major comment 1) is the decisive one. If the CM baseline was evaluated at the lab-frame energy E rather than at the pair rest-frame energy, the 25-40% drag enhancement would not be a genuine frame effect and the paper would need substantial revision. I believe the issue is fixable in a revision, so I do not recommend rejection at this stage. The missing numerical details should be treated as a blocking point for reproducibility, not merely as a stylistic issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nShort version: this is a serious paper with a real technical contribution, but the central numerical comparison has an underspecified energy mapping that could change the headline result. I would send it to a referee, not desk-reject it, but I would ask for a revision that pins down the CM-frame baseline.\n\nWhat's new: the azimuthal decomposition of the in-medium T-matrix about the total-momentum axis is not new machinery—it is standard few-body scattering (refs 57-59)—but applying it to QGP heavy-light scattering and quantifying the lab-vs-CM drag difference is. The derivation is clean, the P=0 limit correctly recovers the CM result, and the framework generalizes to other in-medium problems. The spectral function is admittedly crude (constant width, no self-consistency), but the authors bracket it with Gamma = 100-200 MeV and flag the limitation; that is acceptable for a first pass.\n\nThe soft spot is in Section III. They say they Lorentz-transform lab momenta to get q_CM and theta_CM, but they never state the energy argument used in the CM T-matrix equation (Eq. 6). For the two calculations to describe the same scattering event, T_CM has to be evaluated at E_CM = sqrt(E^2 - P^2), not the lab-frame E. With E=3.05 GeV and P up to 2 GeV, E_CM drops to roughly 2.3 GeV—a big shift. If the CM curves are evaluated at the lab energy, the comparison is partly an energy-shift effect, not a frame effect, and the 25-40% drag correction could be inflated. The stress-test note lands. This is a fixable reporting or calculation issue, but it is load-bearing for the paper's main claim. I would also like numerical details: grid, convergence, truncation of the m-sum. They are absent, which makes the result hard to reproduce as written.\n\nOverall, the thinking is honest and the framework is likely useful. The citation pattern looks appropriate; the self-citation to Ref. [61] is for a standard limit and not a problem. I would not quote the drag number yet. I would send it to a serious referee and ask for the energy mapping to be stated explicitly and the numerics summarized. If the correction survives, this is a meaningful update for heavy-quark diffusion.","headline":"A useful framework for lab-frame in-medium T-matrices, but the headline 25-40% drag correction needs an explicit CM-energy mapping before I would trust the number.","tokens_in":10165,"tokens_out":5217,"would_cite":true,"duration_ms":46138,"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":"Solving the in-medium T-matrix in the medium rest frame raises charm quark drag coefficients by 25-40 percent at low momenta.","keywords":["in-medium T-matrix","laboratory frame","azimuthal decomposition","heavy quark drag","quark-gluon plasma","charm quark","center-of-mass frame","transport coefficients"],"falsifier":"Compute the same charm drag coefficient with a fully self-consistent iterative solution in which the spectral functions are updated together with the T-matrix, and check whether the lab-frame result remains 25-40 percent above the CM-frame result at $p_Q \\le 1$ GeV; alternatively, repeat the on-shell lab-versus-CM amplitude comparison at fixed total momentum with a purely Coulomb potential, since the authors attribute the $\\vec{P}$-dependence to the non-Born part of the kernel.","tokens_in":9162,"feed_emoji":"⚛️","tokens_out":4761,"duration_ms":41663,"temperature":0.7,"pith_summary":"This paper argues that conventional in-medium T-matrix calculations, which solve the scattering equation in the two-particle center-of-mass frame, miss an intrinsic dependence on the total pair momentum that the medium rest frame requires. It develops a practical way to solve the T-matrix directly in the laboratory frame by exploiting the residual azimuthal symmetry around the pair-momentum axis. Applied to charm-light quark scattering in the quark-gluon plasma, this treatment yields drag coefficients that are 25-30 percent larger at T=190 MeV and 30-40 percent larger at T=258 MeV at low momenta, while leaving the temperature dependence essentially unchanged. If correct, existing T-matrix-based extractions of QGP transport properties from heavy-quark data carry a previously unquantified systematic underestimate of tens of percent.","feed_headline":"Lab-frame T-matrix lifts charm drag by up to 40%","feed_subtitle":"Lab-frame scattering adds a 25-40 percent correction to charm drag at low momentum.","key_machinery":"The load-bearing mechanism is an azimuthal decomposition of the laboratory-frame T-matrix about the pair-momentum axis. Because a finite total momentum $\\vec{P}$ breaks the spherical $SO(3)$ symmetry of the scattering state down to $SO(2)$ rotations around $\\vec{P}$, the T-matrix is expanded in azimuthal harmonics $e^{im(\\phi_{\\mathbf{q}'}-\\phi_{\\mathbf{q}})}$. The two-particle propagator is independent of the intermediate azimuthal angle, so the orthogonality of these harmonics decouples the integral equation into independent, uncoupled equations for each component $T_m$. This decomposition makes a five-dimensional, fully $\\vec{P}$-dependent calculation numerically tractable, and it reduces correctly to the ordinary CM-frame T-matrix equation in the $\\vec{P}=0$ limit.","core_discovery":"The central discovery is that the in-medium two-body T-matrix, solved directly in the medium rest frame with full dependence on the total pair momentum $\\vec{P}$ and the complete scattering geometry, differs materially from the conventional CM-frame amplitude, and this difference survives relativistic normalization and thermal phase-space averaging. For charm-light quark scattering, the laboratory-frame drag coefficient exceeds the CM-frame result by 25-40 percent at low heavy-quark momentum, while the temperature dependence of the drag is nearly the same in both frames. The discrepancies trace to the $\\vec{P}$-dependent structure of the two-particle propagator: in the near-forward region the CM approximation stays within about 10 percent of the lab-frame squared amplitude, but in non-forward kinematics the CM-frame result overshoots by roughly 40-150 percent at $P=1$-$1.5$ GeV. After applying the relativistic normalization factor, the forward region also develops large differences of order 100 percent, indicating that the CM-frame shortcut distorts the angular profile of the scattering intensity that feeds transport coefficients.","pith_inferences":["The magnitude of the lab-CM correction is likely model-dependent: the authors trace the effect to the suppression of the $\\vec{P}$-independent Born term relative to the non-perturbative integral term, so a weaker potential or a different screening scale could make the 25-40 percent number larger or smaller.","If the same frame mismatch appears in Sommerfeld-enhanced dark-matter annihilation in the early Universe, where the thermal bath also selects a rest frame, those annihilation rates could carry a similar correction that the CM-frame treatment has not quantified.","A natural testable extension is to compute the lab-frame drag coefficient with momentum-dependent spectral functions obtained from an iterative self-consistent scheme; if the constant-width approximation is the main source of the enhancement, the correction would shrink, whereas if it persists, the pair-momentum dependence itself is robust.","The unphysical displacement of the forward peak in the CM-frame projection suggests that previous T-matrix calculations may have systematically underestimated the contribution of relatively hard, non-forward scatterings to charm quark energy loss, which could matter for azimuthal anisotropy observables."],"forward_implications":["Existing T-matrix-based heavy-quark transport coefficients for the QGP should be revised upward at low momenta by roughly 25-40 percent, with implications for Langevin-type simulations of charm hadron observables.","The same azimuthal-decomposition method can be applied directly to heavy quark-gluon scattering, which is currently missing from the lab-frame calculation and is needed for a complete description of heavy-quark thermalization.","Because the temperature dependence of the drag coefficient is similar in both frames, conclusions about how transport coefficients scale with temperature are preserved, even though the absolute normalization changes.","The pronounced distortion of the forward-scattering profile in the CM-frame projection means that angle-resolved collision kernels entering Boltzmann or Langevin codes should be evaluated with the full pair-momentum dependence rather than a CM-frame amplitude.","The formalism transfers to other many-body systems with a preferred rest frame, such as finite-density nuclear matter, where the laboratory-frame T-matrix likewise dictates the physical observables."],"supporting_citations":[{"why":"Supplies the collision integral definition of the heavy quark drag coefficient used to fold the scattering amplitude into a transport quantity.","marker":"[24]"},{"why":"Provides the constant-width spectral function setup and the CM-frame T-matrix calculation that serves as the baseline for comparison.","marker":"[14]"},{"why":"Establishes the self-consistent T-matrix framework whose constant-width approximation is adopted here, and which is suggested as the future completion of the lab-frame calculation.","marker":"[15]"},{"why":"Supplies the screened Cornell potential whose analytic Fourier transform defines the interaction kernel in the scattering equation.","marker":"[62]"},{"why":"Represents the earlier CM-frame T-matrix application to heavy-quark transport that the lab-frame calculation corrects and extends.","marker":"[18]"},{"why":"Provides the companion T-matrix equation in the $\\vec{P}=0$ limit that the lab-frame framework must recover, serving as a consistency check on the numerical implementation.","marker":"[61]"},{"why":"Documents the strong coupling of partial waves at finite total momentum, motivating the azimuthal decomposition instead of a conventional partial-wave expansion.","marker":"[21]"}],"fun_headline_variants":["Lab-frame T-matrix boosts charm drag by 40%","Scattering in lab frame hikes quark drag up to 40%","CM frame misses up to 40% of charm drag","Quark drag jumps 40% with lab-frame T-matrix","Lab frame corrects charm drag by 40%"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results rest on approximating the quark self-energies by a constant width $\\Gamma=100$-$200$ MeV in the spectral functions rather than solving for them self-consistently, so if the real spectral functions depend strongly on momentum, the reported size of the lab-frame correction could change.","fun_headline_variants_meta":{"raw":{"variants":["Lab-frame T-matrix boosts charm drag by 40%","Scattering in lab frame hikes quark drag up to 40%","CM frame misses up to 40% of charm drag","Quark drag jumps 40% with lab-frame T-matrix","Lab frame corrects charm drag by 40%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000215,"raw_usage":{"total_tokens":1435,"prompt_tokens":956,"completion_tokens":479,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":394}},"tokens_in":572,"tokens_out":479,"duration_ms":4136,"temperature":1.0,"reasoning_tokens":394,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:30:11.226334+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same charm drag coefficient with a fully self-consistent iterative solution in which the spectral functions are updated together with the T-matrix, and check whether the lab-frame result remains 25-40 percent above the CM-frame result at $p_Q \\le 1$ GeV; alternatively, repeat the on-shell lab-versus-CM amplitude comparison at fixed total momentum with a purely Coulomb potential, since the authors attribute the $\\vec{P}$-dependence to the non-Born part of the kernel.","supporting_citations":[{"cited_title":"Svetitsky, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the collision integral definition of the heavy quark drag coefficient used to fold the scattering amplitude into a transport quantity."},{"cited_title":"Riek and R","cited_arxiv_id":null,"evidence_quote":"Provides the constant-width spectral function setup and the CM-frame T-matrix calculation that serves as the baseline for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the self-consistent T-matrix framework whose constant-width approximation is adopted here, and which is suggested as the future completion of the lab-frame calculation."},{"cited_title":"Karsch, M","cited_arxiv_id":null,"evidence_quote":"Supplies the screened Cornell potential whose analytic Fourier transform defines the interaction kernel in the scattering equation."},{"cited_title":"van Hees, M","cited_arxiv_id":null,"evidence_quote":"Represents the earlier CM-frame T-matrix application to heavy-quark transport that the lab-frame calculation corrects and extends."},{"cited_title":"Tiwari and M","cited_arxiv_id":null,"evidence_quote":"Provides the companion T-matrix equation in the $\\vec{P}=0$ limit that the lab-frame framework must recover, serving as a consistency check on the numerical implementation."},{"cited_title":"Schiller, H","cited_arxiv_id":null,"evidence_quote":"Documents the strong coupling of partial waves at finite total momentum, motivating the azimuthal decomposition instead of a conventional partial-wave expansion."}],"review_version":2}