{"id":"1ede17c6-da8d-42cf-88b7-1eb8c1ee6c1f","arxiv_id":"2509.09897","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A Fourier-mode method solves the AMY gluon-splitting equation for anisotropic collision kernels, showing the angular-averaged kernel reproduces the anisotropic rate while the common T* isotropic form deviates by over 10%.","lead":"A new computer method calculates how often gluons split inside the hot, squished matter created in heavy-ion collisions. It finds that a simple angular-average of the plasma's interactions works well, but the standard shortcut used in simulations is noticeably off for slower particles.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sub-0.5% angular-averaging result is obtained from an O(xi)-truncated kernel at xi=1.35, where the exact model is strongly anisotropic; the agreement may be an artifact of comparing a model to its own truncation.","rationale":"The reader flagged the same O(xi) truncation as weakest assumption; this stress-test agrees and sharpens it: at xi=1.35 the expansion parameter is not small, and the specific comparison made is between two quantities derived from the same truncated kernel. Because the central claim is explicitly about how well angular averaging approximates the anisotropic rate, a computation using the full model is required before accepting the conclusion. The numerical method itself appears internally consistent: convergence in nfourier is checked, isotropic limits are recovered, and code availability is claimed. The Eq. (23) typo and lack of comparison to Ref. [45] are secondary; the truncation issue is the load-bearing one. Since the concern is checkable and the authors may already have the exact kernel in their solver, the right action is to require this validation rather than reject the work.","tokens_in":15144,"tokens_out":5185,"duration_ms":45053,"concrete_test":"Recompute the xi=1.35, z=0.1 panel without the O(xi) expansion: use m_D^2(phi) from Eq. (7) directly in Eq. (5), Fourier transform numerically to C(b,phi_b), and solve Eq. (32) with nfourier=7; do the same for the exact angular average <C(b)>_phi. If the exact anisotropic rate differs from the exact angular-averaged rate by more than 1%, or if either differs from the O(xi) result by more than a few percent, the reported <0.5% agreement is an artifact of truncating both sides at the same order and the central conclusion requires qualification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's main physical conclusion—the angular-averaged kernel reproduces the anisotropic rate to <0.5%, while the T* form deviates >10%—is based on rates computed from Eq. (11), the first-order expansion of the collision kernel in xi. The displayed results use xi=0.5, 1.0, 1.35, where O(xi^2) terms are not small. In the exact model (7), m_D^2(phi)=\\bar m_D^2(1-2xi/3+xi cos^2 phi); at xi=1.35 and cos^2 phi=0 this is 0.1 \\bar m_D^2, so the perturbation parameter xi(cos^2 phi-2/3) reaches -0.9. The denominator 1/(q^2(q^2+m_D^2)) is then substantially distorted by the truncation. More importantly, both compared quantities are evaluated with the same truncated kernel: the anisotropic rate and the angular-averaged rate are the O(xi) part plus its own phi-average. Their near-equality could simply mean that the first-order angular modulation has a small integrated effect on the b-space rate, not that the exact anisotropic kernel is well approximated by its angular average. The exact angular-averaged kernel contains O(xi^2) pieces (from averaging the full ratio m_D^4(phi)/(q^2(q^2+m_D^2(phi)))), which are omitted; with these included, the difference from the exact anisotropic rate could exceed 0.5%, and the T* comparison could shift as well. Eq. (23) has an apparent typo, but the decisive issue is the evaluation at xi values beyond the validity of Eq. (11).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a numerical method for computing gluon splitting rates in the AMY formalism when the collision kernel C(q_perp) (equivalently, the dipole cross section C(b)) is anisotropic. The method Fourier-decomposes the angular dependence of the impact-parameter-space function g and the dipole cross section, leading to a coupled system of radial ODEs with boundary conditions fixed by the small-b singular behavior. The method is validated against an isotropic thermal kernel and against variations of b_min, b_max, and the number of Fourier modes. It is then applied to a simple anisotropic model in which the Debye mass depends on angle as m_D^2(phi) = (1 - 2 xi/3 + xi cos^2 phi) bar m_D^2, with the collision kernel expanded to first order in xi. The reported numerical findings are: (i) the anisotropic rate agrees with the rate obtained from the angular-averaged kernel to better than 0.5% for the cases shown; (ii) the commonly used isotropic T_* form deviates from the nonequilibrium rate by more than 10% at small energies; and (iii) seven Fourier modes suffice. The code is publicly available.","tokens_in":15558,"tokens_out":4889,"duration_ms":42391,"significance":"If the main quantitative claim survives scrutiny, the result is practically useful: it suggests that, for this class of observables, an angular-averaged collision kernel may be a much better approximation than the standard T_*-based isotropic kernel for far-from-equilibrium plasmas. The paper has notable strengths: the numerical method is general and clearly described, the code is released, isotropic limits and Fourier-mode convergence are checked, and no parameter is fitted to the target rate. The central physical conclusion, however, rests on a first-order-in-xi expansion that is evaluated at xi values where the expansion parameter is not small, so the reported sub-0.5% agreement and the >10% T_* deviation need additional support before the conclusions can be regarded as robust for the exact model.","major_comments":[{"comment":"The rates shown in Figs. 2 and 3 are computed from the collision kernel expanded to first order in xi, Eqs. (11)-(12), but the results are presented at xi = 0.5, 1.0, and 1.35. At xi = 1.35 the expansion parameter xi(cos^2 phi - 2/3) ranges from about -0.9 to +0.45, so omitted O(xi^2) terms are not numerically small. Since both the 'anisotropic' rate and the 'angular-averaged' rate are evaluated with the same O(xi) truncated kernel, their agreement below 0.5% may only demonstrate that the first-order angular modulation has a small integrated effect, rather than that the exact anisotropic kernel is well approximated by its exact angular average. The angular average of the full kernel contains O(xi^2) contributions that are absent from Eq. (12), and those could change both the comparison to the anisotropic rate and the comparison to the T_* form. I ask the author to add results at small xi (e.g., 0.1 and 0.2), to estimate the O(xi^2) sensitivity at xi = 1.35, or to solve the same model with the full undivided angular kernel (or a resummed/Pade-improved version). The conclusions in Secs. V B and VI should be restated only after such a check.","section":"Sec. III, Eq. (11), used in Sec. V B and Figs. 2-3"},{"comment":"Equation (23) defines D(z,b) = -1/2 [C(b) + C((1-z)b) + C((1-z)b)], with the same argument (1-z)b repeated and the expected C(z b) term missing. The symmetric form required by Eq. (2) should be C(b) + C((1-z)b) + C(z b). If the code follows the printed formula, then the z = 0.1 results in Figs. 2 and 3 are solutions to a different equation. Please correct the typo and confirm, preferably in the text, that the implementation uses the corrected expression.","section":"Eq. (23)"}],"minor_comments":[{"comment":"The sentence 'while still deviating significantly from equilibrium' is slightly ambiguous; it should specify that it is the rate from the angular-averaged kernel that still deviates significantly from the thermal rate.","section":"Abstract and Sec. VI"},{"comment":"The denominators cos phi_b and sin phi_b in Eq. (28) may look singular when phi_b = pi/2 or 0; a parenthetical remark explaining that the expression is evaluated as a limit after combining the two terms would help the reader.","section":"Eq. (28)"},{"comment":"Reference [46] gives no arXiv identifier or persistent DOI; since the code is central to the reproducibility of the numerical claims, a stable identifier or repository link with version information should be provided.","section":"Ref. [46]"},{"comment":"The sentence 'Every system is initialized with exactly one nonzero coefficient' is slightly imprecise because, as stated in the preceding sentence, the coefficients with |n| >= 2 are set to zero for every system; rewording would improve clarity.","section":"Sec. IV B, text before Eq. (41)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the Fourier-mode method is a real advance, and the isotropic numerical crosschecks are credible. The headline claim -- that angular averaging the kernel reproduces the anisotropic rate to better than 0.5% -- is less robust than advertised, because both the anisotropic rate and its angular average are computed from the same O(xi)-truncated kernel at xi = 0.5, 1.0, 1.35. At xi = 1.35 the expansion parameter xi(cos^2 phi - 2/3) ranges from -0.9 to +0.45, so O(xi^2) terms are not small. The stress-test note has this right: the near-agreement may simply mean the first-order angular modulation integrates to a small effect on the rate, not that the exact anisotropic kernel is well represented by its angular average. To make the claim stick, the author should either restrict the numerical test to xi where the O(xi) expansion is controlled, or keep the exact angular dependence in the kernel (even numerically) and compare rates from the full model. The T* comparison could shift as well, since it also uses the truncated kernel as the nonequilibrium reference.\n\nWhat is actually new and good: the decomposition of the angular dependence into Fourier modes in impact-parameter space is a clean, workable extension of the Aurenche-Gelis-Moore-Zaraket and Moore-Schlichting-Schlusser-Soudi methods. The bmax/bmin sensitivity tests in Fig. 1 are the right sanity checks, and the n_fourier = 7 vs 11 agreement gives reasonable confidence in the numerics. The author also correctly flags that the angular-averaging result may be observable-dependent. This is honest and shows the work is not overreaching on its own terms.\n\nSofter issues, in rough order: (1) Eq. (23) really does look like a typo -- D(z,b) should contain C(b) + C(zb) + C((1-z)b), not a repeated C((1-z)b). Even if the code is correct, the displayed equation needs fixing. (2) The paper never compares against the small-xi perturbative result of Hauksson and Gale [45], which would be a natural and cheap validation of the method. (3) Reference [46] for the code is just a title and year, with no URL or arXiv ID; for a numerical methods paper, that undercuts the public-availability claim.\n\nWho it is for: people doing AMY-type rates in QCD kinetic theory or jet quenching in anisotropic media. The method itself is worth having even if the model-test conclusions are scaled back. It deserves a serious referee; my recommendation is to send it out, with the main request being that the author either redoes the analysis at safe xi values or handles the exact kernel, and also fixes the typo and the code-reference issue.","headline":"A genuinely useful numerical method for anisotropic AMY rates, wrapped around a headline physics claim that outruns the O(xi) approximation it is computed from.","tokens_in":16023,"tokens_out":3500,"would_cite":true,"duration_ms":30894,"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":"This paper presents a numerical method for computing gluon splitting rates in an anisotropic plasma, and shows that for its model the angular-averaged collision kernel reproduces the anisotropic rate to better than 0.5%, while the…","keywords":["gluon splitting rate","AMY formalism","anisotropic QCD plasma","collision kernel","dipole cross section","Fourier-mode decomposition","QCD kinetic theory","jet energy loss"],"falsifier":"Recompute the rate with the collision kernel kept to second order in $\\xi$ (or with the full angular-dependent screening mass without expansion) at $\\xi=1.35$ and $p\\sim T$; if the anisotropic rate then differs from the angular-averaged rate by more than 0.5 percent, or if the $T_*$-based rate comes within 10 percent of the nonequilibrium rate, the paper's central numerical conclusions fail.","tokens_in":14941,"feed_emoji":"🧵","tokens_out":6942,"duration_ms":58211,"temperature":0.7,"pith_summary":"This paper presents a numerical method for computing the gluon splitting rate $g\\to gg$ in an anisotropic quark-gluon plasma within the AMY effective-kinetic-theory formalism. Earlier rate calculations assumed an isotropic collision kernel; here the angular dependence of the dipole cross section is expanded in Fourier modes, turning the rate problem into a large coupled system of ordinary differential equations. The method is tested on a simple anisotropic model obtained by giving the Debye mass an angular dependence and expanding to first order in the anisotropy parameter $\\xi$. The central numerical finding is that for this model the rate from the angular-averaged kernel agrees with the full anisotropic rate to better than 0.5%, while the commonly used isotropic kernel based on $T_*$ deviates by more than 10% at small parton energies. This matters because QCD kinetic-theory simulations of early heavy-ion collisions use the isotropic $T_*$ form for inelastic processes.","feed_headline":"Anisotropic gluon splitting rates match angular-averaged kernel","feed_subtitle":"The standard isotropic T* approximation misses the nonequilibrium rate by over 10 percent.","key_machinery":"The central object is the dipole cross section $C(b)$, the Fourier transform of the collision kernel $C(q_\\perp)$, which is the sole medium input in the AMY rate equation. The argument is carried by a decomposition of the angular dependence of both the solution $g(b,\\phi_b)$ and the potential $D(b,\\phi_b)$ into Fourier modes $e^{in\\phi_b}$. This reduces the impact-parameter integral equation to a coupled system of second-order ordinary differential equations for the modes $g_n(b)$, with boundary conditions fixed by the singular small-$b$ behavior and regularity at infinity. An anisotropic model kernel is constructed by giving the Debye mass the angular dependence $m_D^2(\\phi)=(1-2\\xi/3+\\xi\\cos^2\\phi)\\bar m_D^2$ and expanding the collision kernel to first order in $\\xi$, so the comparison between anisotropic, angular-averaged, and $T_*$-based rates tests the method and the physics simultaneously.","core_discovery":"On the paper's own terms, the discovery is that the gluon splitting rate in a plasma with angular-dependent screening is nearly insensitive to the angular structure of the collision kernel once that structure is averaged over directions. For anisotropy parameters $\\xi=0.5, 1.0, 1.35$ and splitting fractions $z=0.1, 0.5$, the anisotropic rate computed with 7 or 11 Fourier modes differs from the angular-averaged-kernel rate by less than 0.5 percent across parton energies from about $0.1\\,T$ to $10^4\\,T$. In contrast, the isotropic approximation used in QCD kinetic-theory simulations, which replaces the kernel by a thermal form with effective infrared temperature $T_*$ and anisotropic Debye mass $m_D$, deviates from the nonequilibrium rate by more than 10 percent, most strongly near the temperature scale. The paper also establishes that the Fourier-mode method itself is general: it applies to arbitrary anisotropic collision kernels, including those obtained from numerical kinetic-theory simulations.","pith_inferences":["If the 0.5-percent insensitivity to angular structure persists for realistic collision kernels, the full angular resolution of the splitting kernel could be replaced by an angular average in integrated-rate calculations, freeing computational resources for other aspects of kinetic-theory simulations.","The near-perfect agreement between anisotropic and angular-averaged rates probably reflects that the rate integrates over all transverse momentum and all emission angles; more differential observables, such as the transverse-momentum spectrum of emitted gluons, are likely to show stronger sensitivity to anisotropy than the integrated rate.","Because the numerical results are obtained at $\\xi$ up to 1.35 from a kernel truncated at $O(\\xi)$, a natural test is to extend the model to second order; if $O(\\xi^2)$ terms shift the rate by more than the claimed 0.5 percent, the angular-averaging conclusion would need qualification at strong anisotropy.","The $T_*$ isotropic form fails mainly in the small-energy region where the formation time is long, which is exactly the regime most relevant for energy loss of moderate-energy jets; this may change the expected shape of jet suppression in the early plasma."],"forward_implications":["The paper argues that QCD kinetic-theory simulations using the isotropic $T_*$ form may misestimate collinear gluon splitting rates in the anisotropic early stages of heavy-ion collisions by more than 10 percent near the temperature scale.","For the class of kernels tested, replacing an anisotropic kernel by its angular average changes the integrated splitting rate by less than 0.5 percent, suggesting that angular averaging is a safe simplification for this particular observable.","The same Fourier-mode method applies to $q\\to qg$ and $g\\to q\\bar q$ splitting rates, which satisfy the same integral equation with different color factors.","The method accepts any anisotropic collision kernel, including realistic kernels generated by kinetic-theory simulations, not just the analytic first-order-in-$\\xi$ model used here.","Convergence of the Fourier series at $n_{\\text{fourier}}=7$ shows that only a small number of angular modes is needed for energy-integrated rates in this parameter range."],"supporting_citations":[{"why":"Supplies the AMY rate formula and integral equation for $g\\to gg$ that the whole calculation solves.","marker":"[15]"},{"why":"Provides the simple anisotropic collision-kernel model with angular-dependent Debye mass that the paper uses as test case.","marker":"[45]"},{"why":"Supplies the impact-parameter-space reformulation and singular boundary-condition technique that the Fourier-mode method generalizes.","marker":"[53]"},{"why":"Describes a recent non-perturbative treatment of the same rate equation for isotropic kernels, providing context and comparison for the numerical approach.","marker":"[49]"},{"why":"Defines the QCD effective kinetic theory and the inelastic collision term in which the isotropic $T_*$ approximation is commonly used.","marker":"[34]"},{"why":"Gives the thermal small-$q_\\perp$ collision-kernel expression that underlies both the anisotropic model and the isotropic comparison form.","marker":"[54]"}],"fun_headline_variants":["Gluon splitting rates match angular-averaged kernel","Anisotropic plasma: angular average suffices for gluon splits","Standard T* fails for nonequilibrium gluon rates","Fourier method shows gluon splittings ignore kernel anisotropy","Angular average predicts anisotropic gluon rates well"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The anisotropic collision kernel is used only to first order in the anisotropy parameter $\\xi$, yet the results are evaluated at $\\xi=0.5$, $1.0$, and $1.35$, where neglected $O(\\xi^2)$ terms could be large; if they are, the reported rates and the 0.5% and 10% comparisons would change.","fun_headline_variants_meta":{"raw":{"variants":["Gluon splitting rates match angular-averaged kernel","Anisotropic plasma: angular average suffices for gluon splits","Standard T* fails for nonequilibrium gluon rates","Fourier method shows gluon splittings ignore kernel anisotropy","Angular average predicts anisotropic gluon rates well"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000687,"raw_usage":{"total_tokens":3097,"prompt_tokens":909,"completion_tokens":2188,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":2106}},"tokens_in":525,"tokens_out":2188,"duration_ms":14213,"temperature":1.0,"reasoning_tokens":2106,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:58:10.925479+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the rate with the collision kernel kept to second order in $\\xi$ (or with the full angular-dependent screening mass without expansion) at $\\xi=1.35$ and $p\\sim T$; if the anisotropic rate then differs from the angular-averaged rate by more than 0.5 percent, or if the $T_*$-based rate comes within 10 percent of the nonequilibrium rate, the paper's central numerical conclusions fail.","supporting_citations":[{"cited_title":"Effective kinetic theory for high temper- ature gauge theories,","cited_arxiv_id":null,"evidence_quote":"Defines the QCD effective kinetic theory and the inelastic collision term in which the isotropic $T_*$ approximation is commonly used."}],"review_version":2}