{"id":"f7ea11d1-f2d9-4070-84d3-933b34317528","arxiv_id":"2608.06788","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a magnetized NJL model, the pi0 spectral function develops multiple magnetic-field-induced peaks and the q-qbar scattering phase shift jumps at threshold energies, with anisotropic dependence on finite meson momentum.","lead":"This paper calculates the spectral function and quark-antiquark scattering phase shift of the neutral pion in a background magnetic field, within an effective Nambu-Jona-Lasinio model. The results show multiple magnetic-field-induced peaks and phase-shift jumps, with different behavior for pion momentum along and perpendicular to the field.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Numerical support for the magnetic-field-induced multi-peak structure is not established: the Landau-level truncation and Pauli-Villars implementation used in Eqs. (14)-(16) are unspecified, so the claimed phase-shift jumps and peaks could be regulator artifacts.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing risk: the numerical predictions are presented without convergence or regulator-sensitivity evidence. A good-faith reading shows the analytic setup is coherent, the relation between the phase shift and spectral function in Eqs. (6)-(8) is standard, and the eB=0 limit in the appendix reproduces the expected single-continuum form. The paper also avoids overclaiming experimental observability and explicitly flags that the peak structure relies on numerical calculation. However, the central qualitative novelty, namely that eB generates multiple peaks and abrupt phase-shift jumps, is not derivable from the analytic formulas alone; it depends on numerically resolving threshold singularities that are regularization-sensitive. Because the manuscript omits implementation details (PV parameters, Landau-level cutoff, quadrature rules), a skeptical reader cannot verify whether the divergent threshold behavior in Tables I-III survives PV regularization, nor whether higher Landau levels change the peak count or jump pattern. This is a real concern, but it does not force rejection: it supports the existing CONDITIONAL verdict, since the issue is unverified numerics rather than demonstrated internal inconsistency. No change to the reader's verdict is needed.","tokens_in":23683,"tokens_out":5652,"duration_ms":64249,"concrete_test":"Reproduce the magenta curve in Fig. 2(a2) (T=mu=0, |k_perp|=0.5 GeV, k3=0, eB=20 m_pi^2) with the same G, Lambda, and m0, while varying the two numerical controls: (i) the number of Landau levels N_LL in Eq. (15) from the original value to 2x and 4x, and (ii) the Pauli-Villars subtraction mass and order by +/-20%. Record the number, positions, and heights of peaks between omega=0 and 1.5 GeV, together with the locations and magnitudes of the phase-shift jumps from Eq. (6). If the peak pattern changes by more than one peak, or if any jump moves by more than a few MeV, then the claimed magnetic-field-induced multi-peak and jump structure is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that eB produces multiple Breit-Wigner/non-Breit-Wigner peaks and abrupt phase-shift jumps. The quantitative content of that claim comes from evaluating the polarization function in Eqs. (14)-(16), which contain infinite sums over Landau levels and principal-value integrals whose threshold singularities are divergent in the unregularized 1+1-dimensional kinematics (Tables I-III). The paper states that Pauli-Villars regularization is used, but the numerical implementation is not shown: no PV subtraction mass or order, no cutoff N_LL for the n and l sums, and no convergence test. If the PV regulator broadens or removes the inverse-square-root threshold singularities, the exact zeros of the spectral function and the abrupt jumps of the phase shift at unitary/Landau thresholds could become finite-width artifacts; if N_LL is too small, the higher-omega peak sequence is truncated and the phase-shift pattern changes. The authors themselves flag that peak structure between thresholds 'rel[y] on the numerical calculations' (Sec. III A) and restrict Fig. 2(c2)/3(c2) to omega < 0.6 GeV because of the large number of thresholds, but they do not provide the robustness checks needed to separate physical magnetic-field structure from cutoff/truncation artifacts. The eB=0 appendix is a useful consistency check, but it does not test this structure because there are no divergent thresholds at eB=0.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the q\\bar q scattering phase shift in the π0 channel, Φπ0(ω²,k⊥²,k3²), and the π0 spectral function ρπ0(ω²,k⊥²,k3²) in the two-flavor Nambu-Jona-Lasinio model with a constant external magnetic field and finite meson momentum. The authors derive the RPA polarization function in a Landau-level basis, identify unitary, Landau, and Pauli-blocking thresholds, and present numerical results for the spectral function and phase shift in three regimes: chiral broken phase (T=μ=0), chiral restored phase at finite T, and chiral restored phase at finite μ. The central claims are that the magnetic field generates multiple Breit-Wigner and non-Breit-Wigner peaks, that the phase shift jumps abruptly at threshold points, and that the dependence on k⊥ versus k3 exhibits magnetic-field-induced anisotropy. An appendix contains the eB=0 limit and shows that the familiar NJL polarization function is recovered.","tokens_in":23934,"tokens_out":5264,"duration_ms":59453,"significance":"If the numerical results are robust, this paper would provide a systematic map of π0 spectral and phase-shift structure in a magnetized NJL model, including analytic threshold classifications that could be useful for further Beth-Uhlenbeck-type thermodynamics. The formal part has real strengths: the phase shift and spectral function are consistently defined from the same polarization function, the threshold tables are explicit and match the stated divergence structure, and the eB=0 limit in Appendix A is a meaningful consistency check. The main limitation is that the central quantitative claims—the multiple peaks and the abrupt phase-shift jumps—rest on numerical evaluations whose regularization and truncation details are not documented, so the physical interpretation as magnetic-field-induced structure is not yet fully supported.","major_comments":[{"comment":"The numerical implementation of the polarization function is not specified. The paper states that Pauli-Villars regularization is used, but does not give the PV subtraction masses, the order of subtraction, the cutoff for the Landau-level sums over n and l, or any convergence test. Since the peak structure between thresholds is explicitly said to 'rely on the numerical calculations' (Sec. III A.2, near the end of the T≠0, μ=0 subsection), and since the Landau-level sums and PV regulator directly control the threshold singularities, the multiple-peak and phase-shift-jump claims could be artifacts of the regulator or of a finite truncation. Please provide the implementation details and convergence checks, including variation of the PV parameters and of the Landau-level cutoff.","section":"Section III, Eqs. (14)–(16)"},{"comment":"The analytic statements that ReΠ diverges to +∞ or −∞ and ImΠ diverges to +∞ at the unitary and Landau thresholds appear to be made for the unregularized principal-value integrals in (14)–(15). If Pauli-Villars regularization is applied before solving the gap equation and evaluating the polarization function, the regularized integrands differ from those shown, and the divergence pattern in Tables I–III may be modified. In particular, PV subtraction can render threshold singularities finite, which would soften the exact zeros of the spectral function and the abrupt phase-shift jumps. The paper needs to show the regularized threshold behavior and demonstrate that the jump pattern of Tables I–III survives regularization.","section":"Tables I–III and Eqs. (14)–(16)"},{"comment":"The abstract and summary emphasize that in the large-ω region the phase shift jumps abruptly at the starting and end points of wide spectral peaks. However, for the T=0, μ≠0, k⊥≠0 case the authors restrict the plots to 0<ω<0.6 GeV because of the large number of thresholds (Sec. III A.3). The large-ω jump structure is therefore not actually documented for this central case. Either show the relevant ω region or state explicitly that the large-ω claim is an extrapolation based on the threshold structure; as written, the conclusion exceeds the presented numerical evidence.","section":"Section III A.3 and Fig. 2(c2)/3(c2)"}],"minor_comments":[{"comment":"The relation ρπ0 = 2G sin(2Φπ0)/(1−2G ReΠπ0) gives zero both below and above a bound-state pole, since sin(2π)=0; the delta peak is not represented by this formula. The text should state explicitly that the bound-state contribution is a separate delta-function term arising from the pole of the propagator, and that Eq. (7) describes only the continuum part.","section":"Eq. (7)"},{"comment":"The notation E_{q+k} is used before being defined; please define E_{q+k} = sqrt(2n|Q_f B| + (q3+k3)² + m_q²) immediately after Eq. (14) or before.","section":"Section II, Eq. (14)"},{"comment":"The threshold labels in the figures are written as e.g. 'U^d(0,0)' but the text refers to ω_U^d(n,l). The figures should use the same notation, e.g. ω_U^d(0,0), to avoid ambiguity.","section":"Figures 2 and 3"},{"comment":"There are several typographical errors, including 'coefficient' in the introduction and 'difficult' in Sec. III A.3; a careful proofread is needed.","section":"Throughout"},{"comment":"Reference [48] (Mei et al., Phys. Rev. D 113, 074031 (2026)) appears to cover closely related ground on pion spectral functions in a magnetic field. The introduction should state explicitly what is new in this paper relative to [48], particularly the phase-shift analysis and the finite-momentum dependence.","section":"Introduction and Ref. [48]"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable model study, but its main quantitative claims are not yet supported because the numerical evaluation is underspecified and no convergence tests are shown. The threshold classification is a useful analytic contribution, but the 'caused by external magnetic field' claim would be much stronger if the authors demonstrated robustness to the PV regulator and Landau-level truncation, and if they compared with a moderate eB value rather than only eB=20mπ². The relation to Ref. [48] should also be clarified; at present the novelty is not sharply delineated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe useful thing here is the systematic map of the pi0 spectral function and the q-qbar phase shift as functions of eB, k_perp, and k3 in the two-flavor NJL model. The zero-momentum spectral function was already in the same group's Ref. [48], but the separate k_perp and k3 dependence, the anisotropy, and the phase-shift jumps at unitary and Landau thresholds are new. The analytical framework in Sec. II is coherent: the S-matrix phase shift and the spectral function both follow from the same polarization function, and the eB=0 limit in Appendix A recovers the known NJL polarization function. The threshold classification (unitary, Landau, Pauli-blocking) is careful and consistent with the divergences they derive. That is real work and it holds up.\n\nThe soft spot is the numerical side. The central claims—multiple Breit-Wigner and non-Breit-Wigner peaks and abrupt phase-shift jumps caused by the magnetic field—are quantitative and come from evaluating Eqs. (14)-(16). No code, no data, and no convergence tests are shown. The Pauli-Villars implementation is described only as \"the Pauli-Villars regularization scheme,\" with no subtraction order or cutoff detail beyond Lambda=1.127 GeV. The Landau-level sums are truncated somewhere, but the truncation point is never stated. The paper itself admits in Sec. III A that the peak structure between thresholds relies on the numerical calculations, and it restricts Fig. 2(c2)/3(c2) to omega < 0.6 GeV because there are too many thresholds. That honesty is good, but it does not replace a cutoff-sensitivity check. If the inverse-square-root threshold divergences are regulator-dependent, the exact zeros of the spectral function and the jumps could be artifacts. This is a legitimate concern, not a manufactured one.\n\nThat said, the broad brush is probably right: at eB=0, no such multi-peak structure appears, so the magnetic field is the cause in that sense. The burden is on the authors to show the structure is independent of the PV regulator and the Landau-level truncation. A revision that adds convergence tests and a few sentences about the numerical implementation would remove most of my hesitation.\n\nWho is this for: people working with NJL/RPA models of magnetized QCD matter, and anyone planning Beth-Uhlenbeck thermodynamic calculations in a magnetic background. It deserves a serious referee, and I would send it to one, with the request that the referee push on the numerical details. I would probably cite it if I worked in that niche, but I would not take the specific peak pattern as established physics until the sensitivity tests appear.","headline":"Solid NJL+RPA extension to finite momentum in a magnetic field, with a clean analytical core and honest caveats, but the numerical peak structure needs convergence tests before the field-induced claims can be trusted.","tokens_in":24530,"tokens_out":2622,"would_cite":true,"duration_ms":25848,"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":"An external magnetic field turns the neutral pion's spectral function into a multi-peak structure with sharp phase-shift jumps.","keywords":["NJL model","pi0 spectral function","quark-antiquark scattering phase shift","external magnetic field","Landau levels","unitary and Landau thresholds","Pauli blocking","anisotropy"],"falsifier":"Recompute the polarization function at $eB=20m_\\pi^2$ with the Pauli-Villars cutoff raised by about 20% and with more Landau levels kept: if the number or positions of the peaks, or the jump pattern of the phase shift, changes qualitatively, the multi-peak structure is a regulator artifact rather than a magnetic-field effect. An independent lattice QCD calculation of the $\\pi^0$ spectral function at the same field strength, showing agreement with the predicted threshold zeroes, would confirm the claim.","tokens_in":23446,"feed_emoji":"🧲","tokens_out":12425,"duration_ms":106468,"temperature":0.7,"pith_summary":"The paper computes, in a two-flavor Nambu-Jona-Lasinio (NJL) quark model, how the neutral pion's spectral function and the quark-antiquark scattering phase shift in the pion channel change under an external magnetic field and finite pion momentum. It claims that the magnetic field reorganizes the spectral function into a bound-state delta peak, several resonance peaks, and several additional broad peaks, and that the scattering phase shift jumps abruptly at the energy thresholds where the quark polarization function becomes singular. These jumps and the multi-peak structure are attributed directly to the magnetic field: they are absent in the zero-field limit that the paper presents in its appendix. The results matter because the spectral function and phase shift are the inputs for understanding how magnetized quark matter's thermodynamics and its pion excitations behave, and because the different response to transverse versus longitudinal momentum is a direct signal of the anisotropy the field induces.","feed_headline":"Magnetic field splinters the neutral pion's spectral peaks","feed_subtitle":"In a quark model, a pion's spectrum gains multiple peaks and the phase shift jumps at magnetic thresholds.","key_machinery":"The load-bearing object is the one-loop quark-antiquark polarization function $\\Pi_{\\pi^0}(\\omega^2,\\mathbf{k}_\\perp^2,k_3^2)$ evaluated in the random-phase approximation, with quark propagators built from Landau levels through the Ritus method and with Pauli-Villars regularization. Everything else is read off from its real and imaginary parts: the $S$-matrix element is $S=(1-2G\\Pi(\\omega-i\\epsilon))/(1-2G\\Pi(\\omega+i\\epsilon))=e^{2i\\Phi_{\\pi^0}}$, the phase shift is $\\Phi_{\\pi^0}=\\arctan[2G\\,\\mathrm{Im}\\,\\Pi/(1-2G\\,\\mathrm{Re}\\,\\Pi)]$, and the spectral function is $\\rho_{\\pi^0}=2G\\sin(2\\Phi_{\\pi^0})/(1-2G\\,\\mathrm{Re}\\,\\Pi)$. The singular thresholds of $\\Pi$ do the causal work: unitary thresholds from the non-crossing terms $\\omega=E_q+E_{q+k}$, where $\\mathrm{Re}\\,\\Pi\\to+\\infty$ from below and $\\mathrm{Im}\\,\\Pi\\to+\\infty$ from above; Landau thresholds from crossing terms $\\omega=|E_{q+k}-E_q|$, whose divergence pattern depends on whether $\\mathbf{k}_\\perp$ or $k_3$ is nonzero; and Pauli-blocking thresholds set by the Fermi step when $T=0$ and $\\mu>0$. The magnetic field enters through the Landau-level sums and the Laguerre-polynomial coefficients $B^\\pm_{nlf}(\\mathbf{k}_\\perp^2)$, which is why the results depend on $\\mathbf{k}_\\perp^2$ and $k_3^2$ separately rather than on $\\omega^2+\\mathbf{k}^2$.","core_discovery":"The paper's central claim is that, at field strength $eB=20m_\\pi^2$, the pion spectral function $\\rho_{\\pi^0}(\\omega^2,\\mathbf{k}_\\perp^2,k_3^2)$ is no longer a single delta-plus-continuum curve but a delta peak for the bound pion, several Breit-Wigner-type resonance peaks, and several non-Breit-Wigner peaks (peaks not tied to a pole of the propagator), with positions controlled by unitary thresholds $\\omega_U(n,l)$ and Landau thresholds $\\omega_L(n,l)$; at $T=0$ with finite quark chemical potential, Pauli-blocking thresholds additionally reshape the peak interiors. The companion claim is that the quark-antiquark scattering phase shift $\\Phi_{\\pi^0}$ is discontinuous: it jumps $0$ to $\\pi$ at the bound-state pole, passes through $\\pi/2$ at resonance poles, and jumps abruptly at unitary and Landau thresholds, exactly where the spectral function vanishes. Both features are caused by the external magnetic field, which quantizes quark momenta into Landau levels and breaks transverse-longitudinal symmetry; accordingly, finite transverse momentum $\\mathbf{k}_\\perp$ and finite longitudinal momentum $k_3$ modify the spectral function and phase shift differently, and only $k_3$ acts like a simple Lorentz shift at $T=\\mu=0$.","pith_inferences":["Beyond the paper, the same threshold-and-jump mechanism should appear in other neutral channels of the same model, such as the sigma meson, because it follows from the Landau-level structure of the polarization function rather than from pion-specific dynamics.","Beyond the paper, the predicted anisotropy between $\\mathbf{k}_\\perp$ and $k_3$ could be tested through angular distributions of dileptons or photons from pion decays in heavy-ion collisions with strong magnetic fields, although the paper does not draw that connection.","Beyond the paper, the phase-shift jumps imply that thermodynamic quantities built from the phase shift via a Beth-Uhlenbeck formula acquire nonanalytic dependence on $eB$ at threshold energies, which could affect fluctuation observables near the chiral crossover."],"forward_implications":["In the magnetized NJL model, the $\\pi^0$ spectral function vanishes exactly at every unitary and Landau threshold, so the multiple peaks are confined to intervals whose boundaries are fixed by the magnetic field and the quark mass.","The $q\\bar{q}$ phase shift jumps from $0$ to $\\pi$ at the bound-state pole and from $\\pi$ to a value below or above $\\pi/2$ at unitary thresholds; with transverse momentum, Landau thresholds force a jump to $0$.","Longitudinal pion momentum $k_3$ at zero temperature shifts the spectrum as $\\omega\\to\\sqrt{\\omega^2+k_3^2}$, while transverse momentum $\\mathbf{k}_\\perp$ does not, giving an explicit signal of magnetic-field-induced anisotropy.","At zero temperature and finite quark chemical potential, Pauli-blocking thresholds, rather than unitary thresholds, set where the wide peaks begin and carve dips inside the resonance peaks.","The computed spectral function and phase shift are the ingredients needed to extend Beth-Uhlenbeck-type thermodynamics of the quark-meson system to the magnetized case, which the paper identifies as a next step."],"supporting_citations":[{"why":"Supplies the two-flavor NJL model and the Pauli-Villars regularization scheme used for the quark propagator and polarization function.","marker":"[28]"},{"why":"Establishes the RPA relation between the meson propagator, the quark-antiquark scattering phase shift, and the meson spectral function.","marker":"[33]"},{"why":"Provides the phase-shift and spectral-function formalism in a meson channel at finite temperature and density that this paper extends to magnetic field and finite momentum.","marker":"[34]"},{"why":"Gives the magnetic-field spectral function of pions and the unitary and Landau threshold structure on which the present calculation builds.","marker":"[48]"},{"why":"Provides the Ritus method for charged-fermion propagators in a magnetic field, used to derive the Landau-level polarization function.","marker":"[42]"},{"why":"Supplies the gauge-independent formulation of the Ritus approach used in the quark-propagator expansion.","marker":"[43]"},{"why":"Fixes the model parameters $G$, $\\Lambda$, and $m_0$ by vacuum condensate, pion decay constant, and pion mass, and gives the pseudo-critical temperature and chemical potential.","marker":"[49]"},{"why":"Contributes the unitary and Landau threshold analysis in a hot, chirally imbalanced NJL medium that the paper adapts to the magnetized case.","marker":"[14]"}],"fun_headline_variants":["Pion spectral peaks multiply in magnetic field","Magnetic field cracks pion spectrum into multiple peaks","Pion phase shift snaps at magnetic field thresholds","Magnetic field triggers multi-peak pion spectrum","Pion spectrum and phase shift reshaped by magnetic field"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim assumes that the multiple peaks and phase-shift jumps are real magnetic-field effects rather than artifacts of the ultraviolet cutoff or of truncating the quark energy-level sums, since the paper does not report convergence or cutoff-sensitivity tests.","fun_headline_variants_meta":{"raw":{"variants":["Pion spectral peaks multiply in magnetic field","Magnetic field cracks pion spectrum into multiple peaks","Pion phase shift snaps at magnetic field thresholds","Magnetic field triggers multi-peak pion spectrum","Pion spectrum and phase shift reshaped by magnetic field"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000733,"raw_usage":{"total_tokens":3498,"prompt_tokens":1381,"completion_tokens":2117,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":997,"completion_tokens_details":{"reasoning_tokens":2059}},"tokens_in":997,"tokens_out":2117,"duration_ms":16870,"temperature":1.0,"reasoning_tokens":2059,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:49:19.178289+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the polarization function at $eB=20m_\\pi^2$ with the Pauli-Villars cutoff raised by about 20% and with more Landau levels kept: if the number or positions of the peaks, or the jump pattern of the phase shift, changes qualitatively, the multi-peak structure is a regulator artifact rather than a magnetic-field effect. An independent lattice QCD calculation of the $\\pi^0$ spectral function at the same field strength, showing agreement with the predicted threshold zeroes, would confirm the claim.","supporting_citations":[{"cited_title":"9 polarization function ΠIm(ω2 δ, k2 ⊥,k 2","cited_arxiv_id":null,"evidence_quote":"Establishes the RPA relation between the meson propagator, the quark-antiquark scattering phase shift, and the meson spectral function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the phase-shift and spectral-function formalism in a meson channel at finite temperature and density that this paper extends to magnetic field and finite momentum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the magnetic-field spectral function of pions and the unitary and Landau threshold structure on which the present calculation builds."},{"cited_title":"With increasing ω (0 < ω < 0.14GeV), Φπ0 increases continuously from zero to π at ω = ωd U (0, 0) − 0+, which crosses π/2 at ω = ωBW < ω d U (0, 0)","cited_arxiv_id":null,"evidence_quote":"Provides the Ritus method for charged-fermion propagators in a magnetic field, used to derive the Landau-level polarization function."},{"cited_title":"It should be mentioned that at PB thresholds, no jumps of scattering phase shift happen","cited_arxiv_id":null,"evidence_quote":"Supplies the gauge-independent formulation of the Ritus approach used in the quark-propagator expansion."},{"cited_title":"Usually, we start from a Lagrangian with current quarks and two-body interaction in scalar and pseudo-scalar channels","cited_arxiv_id":null,"evidence_quote":"Fixes the model parameters $G$, $\\Lambda$, and $m_0$ by vacuum condensate, pion decay constant, and pion mass, and gives the pseudo-critical temperature and chemical potential."}],"review_version":1}