{"id":"8269f81b-deaf-44d3-97b3-ac1ffc005fdd","arxiv_id":"2501.13243","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The first calculation of the valence quark distribution function of the radially excited pion predicts a three-peaked profile with zeroes at x≈0.2 and x≈0.8.","lead":"Physicists used a symmetry-preserving approximation to quantum chromodynamics to calculate parton distribution functions of the pion and its first radial excitation. They predict the excited pion's valence quark distribution has a three-peaked shape with two internal zeroes, a structure never computed before.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The three-peak π1 valence DF is driven by Schlessinger-point-extrapolated Mellin moments m≥6; the paper does not demonstrate that this extrapolation is stable enough to support the endpoint peaks.","rationale":"I read the paper in good faith: it uses a symmetry-preserving CSM framework, reproduces the ground-state pion DF consistently with prior work, and the π1 spectrum and decay constant match empirical expectations, which gives real support to the calculation. The central novel claim, however, is the three-peak hadron-scale π1 valence DF, and the reader correctly identifies the SPM extrapolation of the Mellin moments as the load-bearing assumption. My stress test sharpens this: not only are the high moments extrapolated rather than computed, but the reconstruction parameters, especially a6, change sign and magnitude between the two SPM variants, so the endpoint peaks are evidently controlled by the extrapolated tail. The paper also contains a small internal tension: it says reliable direct π1 moments reach m≤6, but the Table 2 columns are based on SPM from m≤5, which is exactly the order at which the extrapolation begins. A direct computation of m=6 and a refit with or without the SPM-extended moments would settle whether the three peaks survive. This is a robustness concern, not an internal inconsistency of the framework, so conditional acceptance remains the right verdict; I would not reject or fully accept on this basis.","tokens_in":17253,"tokens_out":7194,"duration_ms":69982,"concrete_test":"Using the same quark propagator and Bethe-Salpeter amplitudes, compute the π1 Mellin moments m=6 and m=7 directly from Eq. (14b) with increased numerical precision, bypassing the SPM, and compare with Table 2; then refit Eq. (19) to the moment set m=0,...,5 with the same non-negativity constraint. If the directly computed m=6 differs from the SPM value by more than the quoted 0.1%, or if the refit without SPM-extended moments no longer exhibits zeros near x=0.2,0.8 with secondary peaks near x=0.1,0.9, then the three-peak valence DF is an artifact of the extrapolation rather than a robust prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4 states that standard algorithms provide reliable π1 moments for m≤6, yet the two π1 columns of Table 2 are described as SPM extrapolations based on moments 0≤m≤5. Hence the m=6,...,10 entries, and especially the m≥8 entries that lie above the scale-free moments, come from the Schlessinger point method, not from direct computation. Section 5 reconstructs qπ1 via Eq. (19), whose only shape parameters beyond the baseline are a2,a4,a6; the zeros near x≈0.2,0.8 and the secondary peaks near x≈0.1,0.9 are the mechanism by which the low moments stay below and the high moments above the scale-free DF. Table 3 shows that the two SPM variants give a6=-0.003 and a6=+0.0684, a large relative difference, while the paper claims the resulting DFs are qualitatively identical; this indicates the endpoint structure is controlled by the extrapolated tail. No uncertainty from the SPM is propagated into the reconstructed DF, and the sentence 'in nontrivial test cases, M(z) returns moments out to m=10 whose relative error is <0.1% in magnitude' is not a stability test for this particular five-point input set. If the SPM overestimates the m≥6 moments, the three-peak structure would weaken or disappear, so the paper's headline claim rests on the least-secure step of the analysis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript computes the hadron-scale valence quark distribution functions of the pion ground state (π0) and its first radial excitation (π1) using a symmetry-preserving continuum Schwinger function framework. The valence Mellin moments are obtained from dressed quark propagators and Bethe-Salpeter amplitudes, the pointwise DFs are reconstructed using a parametrized form, and all DFs (valence, glue, sea) are evolved to 3.2 GeV with an all-orders evolution scheme. The central prediction is that the π1 valence DF has a three-peak structure at the hadron scale, with zeros near x≈0.2 and x≈0.8 and secondary peaks near x≈0.1 and x≈0.9, in contrast to the single-peak ground-state pion DF.","tokens_in":17652,"tokens_out":5924,"duration_ms":57132,"significance":"If the three-peak prediction is robust, it is a novel and testable statement about the structure of excited pseudoscalar mesons, and the paper provides useful benchmarks for other nonperturbative frameworks and for future lattice QCD studies. The analysis has several credible ingredients: an EHM-improved Bethe-Salpeter kernel that reproduces the meson spectrum, a reconstruction of the ground-state pion DF consistent with existing results, and an all-orders evolution scheme that is applied to all parton species. The main weakness is that the headline feature rests on Schlessinger-point extrapolation of Mellin moments, and the paper does not demonstrate that this extrapolation is stable for the specific five-point input set used for the π1.","major_comments":[{"comment":"The three-peak structure of the π1 valence DF in Fig. 3 is driven by the m≥6 (especially m≥8) π1 moments being larger than the scale-free moments. These high moments are obtained by SPM extrapolation from the m≤5 computed moments. The text states that standard algorithms provide reliable access to all m≤6 moments for the π1, yet neither the π(a)1 nor the π(b)1 column in Table 2 uses the m=6 moment as an input; the m=6 entry is itself an extrapolated value. The validation sentence in Section 4 ('in nontrivial test cases, M(z) returns moments out to m=10 whose relative error is <0.1% in magnitude') is a test of the method on other functions, not on this particular five-point input set. Since the endpoint peaks in the reconstructed DF are the central claim, the authors should provide a direct stability test: include the actually computed m=6 moment as a check, vary the input set used for the SPM, compare the SPM predictions with the computed moments, and propagate the resulting uncertainty into the reconstructed DF.","section":"Section 4, Table 2"},{"comment":"The two SPM variants (columns π(a)1 and π(b)1) yield reconstruction parameters a6 = -0.003 and a6 = +0.0684, a substantial relative difference in the Gegenbauer coefficient that controls the oscillatory structure in Eq. (19). The paper states that the resulting DFs are 'qualitatively identical,' but this is not quantified, and the physical claim is precisely about the qualitative features: the zeros near x≈0.2, 0.8 and the secondary peaks near x≈0.1, 0.9. Without an uncertainty envelope for the reconstructed DF, or a demonstration that all plausible reconstructions preserve the three-peak structure, the robustness of the headline prediction is not established. Please provide the range of reconstructed DFs obtained from the two moment sets and state explicitly whether the three peaks and the zero locations are stable under this variation.","section":"Section 5, Table 3"},{"comment":"The discussion that dismisses the lattice QCD results of Ref. [81] relies on the ordering of the hadron-scale moments in Table 2: the paper argues that the lattice ordering ⟨x^m⟩_{qπ1} > ⟨x^m⟩_{qπ0} at ζ3 is incompatible with the moments in Table 2 and with DGLAP evolution. This conclusion is contingent on the hadron-scale moments themselves, which for the π1 depend on the SPM extrapolation. If the SPM overestimates the higher π1 moments, the ordering of π0 vs π1 moments at ζ3 could be reversed. The paper should either treat the lattice comparison as an indication rather than a definitive incompatibility, or show explicitly that the SPM uncertainty is too small to affect the ordering conclusion.","section":"Section 6, paragraph after Fig. 4"}],"minor_comments":[{"comment":"Please clarify why the reliable m=6 moment for the π1 is not used as an input to the SPM, given the statement that standard algorithms provide access to all m≤6 moments; this would strengthen the credibility of the extrapolation.","section":"Section 4"},{"comment":"The summary says 'we computed eleven Mellin moments of the π0,1 valence DFs.' For the π1, moments m=6,...,10 are extrapolated rather than directly computed; 'obtained' or 'determined' would be more precise.","section":"Section 7 and abstract"},{"comment":"The caption contains a placeholder '[?]' for the G-parity reference; please supply the full citation.","section":"Table 4 caption"},{"comment":"The two π1 reconstructions are very close; consider adding a magnified inset near the endpoint peaks (x≈0.1 and x≈0.9) so that the secondary peaks are clearly distinguishable.","section":"Figure 3"},{"comment":"The relation qπ0 ∝ |φπ0|^2 is cited but not derived; for the excited state the qualitative use is reasonable, but the sentence locating the π1 DF zeros 'in the vicinity of' the DA zeros should be supported by a quantitative comparison of the zero positions.","section":"Section 5, Eq. (20)"},{"comment":"The parameters ω, D, η, and the current masses are model inputs; the paper would benefit from a statement on the sensitivity of the DF predictions to their variation, or a reference to the earlier study where these values were determined.","section":"Section 3, Eqs. (11)-(13)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is part of an established program on pion structure using continuum Schwinger function methods. The central claim—the three-peak π1 valence DF—is novel and would be of interest to the hadron structure community. However, the claim is directly contingent on the Schlessinger-point extrapolation of the high Mellin moments, and the paper does not currently provide a robust validation of that step for the specific input set. This is fixable within the manuscript's scope by adding stability tests and uncertainty propagation, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the first calculation of the hadron-scale valence quark DF for the first radial excitation of the pion, π1. The prediction of a three-peak structure, with zeros near x=0.2,0.8 and secondary peaks near 0.1,0.9, is striking and is the kind of concrete, falsifiable output that should push lattice and other frameworks. The paper also does several things well: the symmetry-preserving kernel with a dressed-quark anomalous chromomagnetic moment fixes a known failure of rainbow-ladder for excited states, the ground-state pion DF comes out consistent with earlier work, and the authors are appropriately cautious about using RL truncation for π1. The calculation of Mellin moments directly from the Bethe-Salpeter amplitude, with recursion relations checked, is solid.\n\nThe soft spots are real, though. First, the three-peak structure depends on the high Mellin moments, m≥6, and those come from a Schlessinger point method extrapolation based only on m≤5 moments. The paper's validation of SPM is limited to test cases, not a stability study for this specific five-point input. The two SPM variants give a6 = -0.003 vs +0.0684, a large relative difference, and the claim that the resulting DFs are qualitatively identical does not establish that the endpoint peaks survive. No uncertainty from the SPM is propagated into the reconstructed DF. Second, the reconstruction ansatz in Eq. (19) builds in the three-peak shape through the Gegenbauer terms; the zeros are effectively put in by hand, and the connection to the zeros of the distribution amplitude is suggestive but not a derivation. Third, the paper dismisses the exploratory lattice study that finds the opposite ordering of moments, saying only that more lattice work is needed. That is not fatal, but it is a loose end.\n\nNone of this makes the central claim circular—the DF is not fit to external DF data, and the framework has a good track record—but it does mean the headline prediction is only as solid as the least-secure step, and that step is not adequately stress-tested.\n\nWho is this for? Anyone working on hadron structure, especially pion excitations, lattice QCD, or phenomenological parton distributions. It deserves a serious referee: the novelty is clear, the calculation is a proper first, and the prediction is concrete. I would send it to review with a request for a sensitivity analysis of the SPM (vary the input moments, try alternative interpolation schemes, show whether the endpoint peaks survive), plus a more careful discussion of the lattice discrepancy and a propagated uncertainty band on the three-peak structure.","headline":"First hadron-scale valence DF for the pion's first radial excitation, with a novel three-peak structure that is intriguing but rests on SPM-extrapolated high moments whose stability is not demonstrated.","tokens_in":18139,"tokens_out":2378,"would_cite":true,"duration_ms":26317,"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":"This paper predicts that the pion's first radial excitation has a three-peak valence quark distribution at the hadron scale, with zeros near x≈0.2 and x≈0.8.","keywords":["pion distribution functions","radial excitation","Mellin moments","parton distributions","Bethe-Salpeter equation","dynamical chiral symmetry breaking","emergent hadron mass","three-peak structure"],"falsifier":"An independent determination of the hadron-scale $\\pi_1$ valence Mellin moments—for instance from a lattice or light-front calculation—that finds the $m\\ge8$ moments below the scale-free values, or a pointwise reconstruction without zeros near $x\\approx0.2$ and $x\\approx0.8$, would falsify the three-peak claim. In the paper's Table 2, the reconstruction needs $\\langle x^8\\rangle$ at least near the scale-free value $0.0350$.","tokens_in":17077,"feed_emoji":"⚛️","tokens_out":11681,"duration_ms":111075,"temperature":0.7,"pith_summary":"Working from a symmetry-preserving approximation to QCD's bound-state equations, this paper predicts the parton distribution functions of the ground-state pion and its first radial excitation. Its central result is that the excited pion's valence quark distribution at the hadron scale is not a single bell-shaped curve: it has three peaks, with a central maximum near $x=1/2$, zeros near $x\\approx0.2$ and $x\\approx0.8$, and secondary peaks near $x\\approx0.1$ and $x\\approx0.9$. This shape follows from the pattern of chiral symmetry breaking, and it is the kind of structural fact that distinguishes a radially excited bound state from its ground state. If the prediction holds, it gives other nonperturbative approaches a sharp target: the same moments, the same zeros, the same peak locations.","feed_headline":"First excited pion shows three-peak valence quark structure","feed_subtitle":"Unlike the ground state, the radial excitation's valence distribution has three peaks and two zeros.","key_machinery":"The carrier of the argument is the Mellin-moment representation of the hadron-scale valence quark distribution, computed from the light-front projection of the Bethe-Salpeter amplitude and the dressed quark propagator. Symmetry fixes the $m=0,1$ moments; the paper computes moments up to $m=5$ (or $6$) directly, extends the sequence to $m=10$ using the Schlessinger point method, and reconstructs the pointwise distribution from an ansatz that respects the QCD endpoint behaviour $(1-x)^2$. For the excited state, the load-bearing physical ingredient is the Bethe-Salpeter kernel built from a process-independent effective charge and a dressed-quark anomalous chromomagnetic moment; this kernel is what fixes the failure of rainbow-ladder truncation, which violates the Cauchy-Schwarz stability condition for a non-negative distribution. Evolution to $\\zeta=3.2\\,$GeV is performed with an all-orders scheme, which also produces the glue and sea distributions.","core_discovery":"On the paper's own terms, the discovery is that the hadron-scale valence distribution function of the pion's first radial excitation, $q_{\\pi_1}(x;\\zeta_H)$, has three peaks. The $m=0$ and $m=1$ Mellin moments of $\\pi_0$ and $\\pi_1$ are forced equal by baryon-number and momentum conservation, while every moment with $m\\ge2$ is larger for the ground state. The computed low-order moments are followed by a Schlessinger-point extension to $m\\le10$ and a pointwise reconstruction, which produces a central peak at $x=1/2$ that is taller and narrower than the scale-free profile $30x^2(1-x)^2$, zeros at $x\\approx0.2$ and $x\\approx0.8$, and secondary peaks near $x\\approx0.1$ and $x\\approx0.9$. The paper interprets this as the momentum-space analogue of a first radial excitation: central momentum is favoured, the extreme all-or-nothing fractions are suppressed, and the support that would sit near $x\\approx0.2,0.8$ is pushed outward to the endpoint shoulders.","pith_inferences":["If confirmed, the three-peak structure would give searches a cleaner discriminator between ground-state and excited pion structure than mass spectra alone, since two Hamiltonians with similar spectra can have very different momentum-space wave functions.","The endpoint peaks are controlled by the extrapolated moments $m\\ge8$; a future calculation that resolves these moments directly would settle whether the peaks are real or an artifact of the analytic continuation.","A similar three-peak pattern could be expected for other radially excited pseudoscalar mesons, because the mechanism is the node structure of the excited-state wave function combined with the symmetries of the pion-like system.","The paper's argument implies that the $\\pi_1$ valence distribution at the hadron scale is effectively the modulus-squared of its light-front amplitude; if one measures or computes the amplitude's zero positions accurately, the DF zeros should coincide with them."],"forward_implications":["A radially excited pseudoscalar meson can have a valence distribution that is not bell-shaped: the $\\pi_1$ DF possesses three peaks, with the outer peaks separated from the central maximum by zeros near $x\\approx0.2$ and $x\\approx0.8$.","The moment ordering is robust under evolution: $m=0,1$ moments remain equal and every $m\\ge2$ moment of $\\pi_0$ stays larger than the $\\pi_1$ partner at $\\zeta=3.2\\,$GeV, so the valence difference survives to measurable scales.","Glue and sea distributions of $\\pi_0$ and $\\pi_1$ at $3.2\\,$GeV are nearly identical, so the radial excitation's fingerprint is carried almost entirely by its valence distribution.","The relation $q(x)\\propto|\\varphi(x)|^2$, already known for the ground-state pion, holds at least qualitatively for the radial excitation: the zeros of the $\\pi_1$ distribution amplitude line up with the zeros of its valence DF.","Because the reconstruction respects non-negativity and the $(1-x)^2$ endpoint behaviour, the three-peak shape, if correct, is not an artefact of the polynomial expansion but a consequence of the high-order moments lying above the scale-free values."],"supporting_citations":[{"why":"Gives the light-front formula (Eq. 7) used to compute the valence quark DF moments from the Bethe-Salpeter amplitude and dressed quark propagator.","marker":"[32]"},{"why":"Supplies the Mellin-moment methodology, the recursion relations, and the earlier pion DF analysis whose moment extension procedure is followed here.","marker":"[34]"},{"why":"Defines the EHM-improved Bethe-Salpeter kernel with a dressed-quark anomalous chromomagnetic moment, the ingredient that makes excited-state calculations reliable; rainbow-ladder truncation fails for the $\\pi_1$.","marker":"[23]"},{"why":"Provides the distribution amplitudes of $\\pi_0$ and $\\pi_1$, including the two zeroes of the $\\pi_1$ amplitude that the reconstructed DF's zeros are compared with.","marker":"[24]"},{"why":"Establishes the all-orders evolution scheme for DFs and its application to pion parton distributions, used here to evolve from $\\zeta_H$ to $3.2$ GeV.","marker":"[33]"},{"why":"Gives the process-independent effective charge whose screening mass fixes the hadron scale $\\zeta_H = 0.331(2)$ GeV used in the evolution.","marker":"[48]"},{"why":"Introduces the logarithmic reconstruction ansatz for the ground-state pion DF that is extended with a Gegenbauer expansion for the excited state.","marker":"[29, 79]"},{"why":"Provides the Schlessinger point method used to extend the moment sequence from $m\\le5$ to $m\\le10$.","marker":"[74-77]"},{"why":"The lattice exploratory study of $\\pi_0$ versus $\\pi_1$ valence structure whose moment ordering the present calculation argues is incompatible with its results.","marker":"[81]"}],"fun_headline_variants":["Excited pion's valence quarks show three peaks","First radial pion excitation has triple-humped quark distribution","Pion's first excited state reveals three-peak valence structure","Three-peak valence distribution marks pion's first radial excitation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the analytic extension from the five directly computed Mellin moments to the tenth is reliable; the outer peaks only appear if the extrapolated high-order moments stay at or above the scale-free baseline.","fun_headline_variants_meta":{"raw":{"variants":["Excited pion's valence quarks show three peaks","First radial pion excitation has triple-humped quark distribution","Pion's first excited state reveals three-peak valence structure","Three-peak valence distribution marks pion's first radial excitation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000258,"raw_usage":{"total_tokens":1666,"prompt_tokens":1114,"completion_tokens":552,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":730,"completion_tokens_details":{"reasoning_tokens":486}},"tokens_in":730,"tokens_out":552,"duration_ms":6098,"temperature":1.0,"reasoning_tokens":486,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:19:50.077528+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent determination of the hadron-scale $\\pi_1$ valence Mellin moments—for instance from a lattice or light-front calculation—that finds the $m\\ge8$ moments below the scale-free values, or a pointwise reconstruction without zeros near $x\\approx0.2$ and $x\\approx0.8$, would falsify the three-peak claim. In the paper's Table 2, the reconstruction needs $\\langle x^8\\rangle$ at least near the scale-free value $0.0350$.","supporting_citations":[{"cited_title":"Chang, C","cited_arxiv_id":null,"evidence_quote":"Gives the light-front formula (Eq. 7) used to compute the valence quark DF moments from the Bethe-Salpeter amplitude and dressed quark propagator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Mellin-moment methodology, the recursion relations, and the earlier pion DF analysis whose moment extension procedure is followed here."},{"cited_title":"Xu, Z.-Q","cited_arxiv_id":null,"evidence_quote":"Defines the EHM-improved Bethe-Salpeter kernel with a dressed-quark anomalous chromomagnetic moment, the ingredient that makes excited-state calculations reliable; rainbow-ladder truncation fails for the $\\pi_1$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the distribution amplitudes of $\\pi_0$ and $\\pi_1$, including the two zeroes of the $\\pi_1$ amplitude that the reconstructed DF's zeros are compared with."},{"cited_title":"Cui, J.-L","cited_arxiv_id":null,"evidence_quote":"Gives the process-independent effective charge whose screening mass fixes the hadron scale $\\zeta_H = 0.331(2)$ GeV used in the evolution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The lattice exploratory study of $\\pi_0$ versus $\\pi_1$ valence structure whose moment ordering the present calculation argues is incompatible with its results."}],"review_version":1}