{"id":"63f347e6-0629-4442-9cf2-bac159d473b3","arxiv_id":"2507.07523","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"A light-front quark model with Gaussian wavefunctions, calibrated on decay constants, reproduces pion and kaon form factors and predicts systematic quark-mass-asymmetry trends for heavier meson radii.","lead":"This paper calculates the electromagnetic form factors and charge radii of pions, kaons, and heavy mesons using a light-front quark model. The results match measured pion and kaon data at low momentum transfer and predict how quark mass differences shape the charge distributions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed proportionality between peak Q2F and Δm, and the claimed monotonic decrease of charge radii with meson mass, are contradicted by Table III: the pion has a nonzero peak at Δm=0, and B+ has a larger radius than D+.","rationale":"The Reader's formal weakest_assumption concerned the one-body current/zero-mode approximation and the Gaussian wavefunction. While that is a legitimate model-dependence caveat, the more decisive issue is internal consistency: the paper's headline trend claims are contradicted by its own Table III. The proportionality of peak Q2F to Δm cannot accommodate the pion, and the claimed monotonic decrease of charge radii with meson mass fails for B+ versus D+. The Reader did note the loosely quantified proportionality in the rationale, but did not make it the primary concern. The underlying numerical computation is a standard LFQM calculation, and the π/K charge radii agree with experiment, which is real evidence in favor of the model. The concern does not invalidate the computation; it requires the authors to soften or re-derive the trend claims in the abstract and conclusions. Therefore the verdict should remain CONDITIONAL, i.e., UNCHANGED relative to the Reader's decision.","tokens_in":8607,"tokens_out":12843,"duration_ms":147051,"concrete_test":"Take Table III's columns 2 and 3 for the six charged mesons (π+, K+, D+, D_s, B+, B_c) and perform a linear regression of peak Q2F on Δm, both with a free intercept and with the intercept fixed to zero. Report the best-fit intercept, its uncertainty, and the chi-squared per degree of freedom; if the free intercept is not consistent with zero (or the zero-intercept fit has chi2/dof >> 1), the 'approximately proportional' claim is refuted. Additionally, check the monotonicity claim by comparing the B+ and D+ entries: B+ has a larger meson mass but a larger charge radius, directly contradicting 'decrease with increasing meson mass' unless the Δm dependence is explicitly separated in a multivariate fit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and conclusions state two central trends: (i) for charged mesons, the peak values of Q2F_P(Q2) are approximately proportional to the constituent-quark mass difference Δm; (ii) charge radii decrease with increasing meson mass and decreasing Δm. Table III does not support either statement as written. The pion has Δm=0 but a peak Q2F of 0.361, so a proportional relation cannot pass through the origin; among the other charged mesons, the ratio peak/Δm varies from about 1.3 (B_c) to about 2.0 (D_s, D+), roughly a 50% spread. The radius trend also fails: B+ (m≈5.28 GeV, sqrt<r2>=0.564 fm) has a larger charge radius than D+ (m≈1.87 GeV, sqrt<r2>=0.411 fm), despite the larger meson mass. In Section III.B the authors attempt a heuristic formula, but it is not derived, not fitted, and not used to qualify the abstract's 'approximately proportional' claim. These trend claims are load-bearing because they are the paper's main phenomenological message beyond the π/K comparison, and they are stated as definitive findings without a quantitative criterion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript computes the electromagnetic form factors and charge radii of pseudoscalar mesons (π, K, D, D_s, B, B_s, B_c) in the light-front quark model with Gaussian wavefunctions, Eq. (6). The model parameters (constituent quark masses and Gaussian widths β) are fixed by mesonic decay constants rather than by form-factor data. The authors report good agreement with low-Q² pion and kaon data, compare their radii with lattice QCD and other models, and propose two global trends: the peak of Q²F_P(Q²) is approximately proportional to the constituent mass difference Δm, and charged-meson charge radii decrease with increasing meson mass. They also claim that neutral-meson radii are governed by the heavy-quark charge and call for future experimental tests.","tokens_in":8867,"tokens_out":6440,"duration_ms":70522,"significance":"If the central results hold, the paper provides a unified set of predictions for heavy-meson form factors and charge radii from a simple, transparent model whose parameters are not fitted to the form-factor data themselves. The honest comparison with NA7, JLab, and recent lattice data for π and K, and the explicit agreement of the pion charge radius with the PDG value, are genuine strengths. The paper is also careful to list parameter uncertainties and to compare with several other approaches. The significance is diminished, however, by the overstatement of the global trend claims, which are not supported by the table that supposedly demonstrates them; these claims are the paper's main phenomenological message beyond the π/K comparison.","major_comments":[{"comment":"The statement that the peak values of Q²F_P(Q²) are approximately proportional to Δm is contradicted by Table III as written. The pion has Δm=0 but a peak Q²F of 0.361, so a proportionality through the origin cannot hold. Among the other charged mesons, the ratio peak/Δm varies from about 1.32 (B_c) to 2.01 (D_s), a spread of roughly 50%. The same issue affects the neutral-meson statement in the Conclusions. The authors should either fit a quantitative relation (e.g., with an intercept or a power law) and quote its quality, or soften the claim to a qualitative ordering.","section":"Abstract and Section III.B, Table III"},{"comment":"The claim that charged-meson charge radii decrease with increasing meson mass is not supported by Table III: B⁺ (m≈5.28 GeV) has √<r²>=0.564 fm, while D⁺ (m≈1.87 GeV) has 0.411 fm. The proposed phenomenological relation √<r²> ∝ e_h Δm_l / β^{3/2} m_s^h in Section III.B is not derived, not fitted, and not used quantitatively, and its symbols are undefined. The authors should replace the global monotonic statement with the actual dependencies on β, quark charges, and mass asymmetry, or present a quantitative fit that is actually compared with the table.","section":"Section III.B and Table III"},{"comment":"The claim of 'obvious asymptotic behavior earlier than LQCD' is misleading. With the Gaussian wavefunction Eq. (6), the overlap integral in Eq. (4) acquires an exp(-x̄² Q²/(4β²)) suppression at large q⊥, which does not reproduce the pQCD 1/Q² tail that the authors themselves cite from Ref. [67]. Footnote 1 concedes a possible higher-twist explanation 'in preparation'. Until that analysis is provided, the high-Q² behavior should be described as model-dependent rather than asymptotic.","section":"Section III.A, Fig. 1 and footnote 1"},{"comment":"The calculation keeps only the one-body valence diagram in the q⁺=0 frame and neglects zero-mode and pair contributions. This approximation is known to be delicate for neutral mesons and at endpoint regions x→0,1. Since the charged/neutral endpoint differences are a central result of the paper, the authors should either provide a consistency check (e.g., comparison with the covariant light-front approach or with an alternative current component) or explicitly state this limitation in the discussion of the endpoint behavior.","section":"Section II, Eqs. (1)–(4)"}],"minor_comments":[{"comment":"The phrase 'parameters derived from the confinement of mesonic decay constants' should be 'parameters fitted to mesonic decay constants' or 'constrained by mesonic decay constants'.","section":"Abstract"},{"comment":"The sentence 'It is straightly to extract FP(Q²) from Eq.(1)' should read 'It is straightforward to extract FP(Q²) from Eq. (1)'.","section":"Section II"},{"comment":"The uncertainties in the Gaussian parameters β are quoted in Table II, but the text does not explain how these uncertainties are propagated into the charge radii in Table III; a brief statement would improve reproducibility.","section":"Table II and Table III"},{"comment":"The symbol '/s32' appearing in the figure legend and caption appears to be a rendering artifact; the meson labels should be cleanly typeset.","section":"Fig. 2"},{"comment":"The sentence 'the one ⟨r²⟩_{K0} = −0.090 fm²' should be 'our ⟨r²⟩_{K0} = −0.090 fm²'.","section":"Section III.A"},{"comment":"The empirical relation √<r²>_P ∝ e_h Δm_l / β^{3/2} m_s^h uses undefined symbols and appears dimensionally inconsistent as written; the authors should define all symbols and state the intended dimensionality, or remove the formula.","section":"Section III.B"},{"comment":"Reference [59] is a self-citation to a preprint by the same authors; its status (published or in preparation) should be clarified in the bibliography.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The numerical implementation appears sound, and the π/K comparisons with data and lattice QCD are the paper's strongest part. The main reason for major revision is not the calculation itself but the mismatch between the abstract/conclusions and Table III: the 'proportionality' and 'decrease with meson mass' statements are load-bearing and are not supported by the presented numbers. I see no grounds for rejection, because these issues can be fixed by reframing the trends quantitatively and by softening the language. The self-citation [59] is not problematic in my view."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The honest take: this is a clean, standard light-front quark model calculation of all pseudoscalar meson electromagnetic form factors using a single parameter set fitted to decay constants rather than to the form-factor data. The pion and kaon results agree well with experiment in the low-Q² region, and the charge radii are in line with PDG and other model predictions. The new content is the breadth—a common treatment of π, K, D, B, Bs, Bc—and the systematic comparison with other approaches. As a reference calculation, it is useful.\n\nThe soft spots are real but localized. The abstract and conclusions make two trend claims: (i) peak Q²F is approximately proportional to Δm, and (ii) charge radii decrease with meson mass. Table III contradicts both as stated. The pion has Δm = 0 but a peak of 0.361, so the proportionality cannot pass through the origin. And B+ has a larger radius than D+ (0.564 fm vs 0.411 fm) despite being much heavier. The heuristic formula in Sec. III.B is not derived, not fitted, and not used to qualify anything; it reads as a rough observation, not a result. These claims should be softened, quantified with a fit, or removed. This is a load-bearing issue for the paper's phenomenological message, but it does not invalidate the underlying calculation.\n\nTwo other points. First, the high-Q² behavior: the Gaussian tail matches neither pQCD nor the lattice trend, and the authors concede in a footnote that the asymptotic onset is earlier than LQCD. The abstract's phrase \"similar asymptotic behavior\" is a stretch. Second, the one-body current approximation is standard, but for neutral mesons the zero-mode/pair contributions are known to be delicate, and the imaginary charge radii are a model-dependent way to report negative mean-square radii; the text does not discuss this limitation.\n\nWho is this for? People working on LFQM and hadron structure who want a consistent set of pseudoscalar charge radii and form factor shapes across the whole family. It deserves serious peer review because the calculation is sound and the comparisons are useful, even though the interpretive trends need revision. A good referee should require either a quantitative criterion for the proportionality or removal of the trend statements from the abstract.\n\nSend it to review. It is not a desk reject.","headline":"A solid, standard LFQM calculation of pseudoscalar EMFFs whose two headline trend claims are overstated in the abstract and conclusions.","tokens_in":9438,"tokens_out":1787,"would_cite":true,"duration_ms":22721,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.40.Gp","12.39.Ki"],"model":"deepseek-v4-flash","headline":"The paper argues that a single light-front quark model with Gaussian wavefunctions and decay-constant-fitted parameters reproduces the measured pion and kaon form factors and predicts how quark mass asymmetry controls the charge radii and…","keywords":["light-front quark model","pseudoscalar mesons","electromagnetic form factors","charge radii","quark mass asymmetry","pion form factor","kaon form factor","heavy mesons"],"falsifier":"Measure the D+ or B+ electromagnetic form factor at JLab or the EIC up to $Q^{2}$ of roughly 10 $GeV^{2}$ and compare the peak of $Q^{2}$ F_P($Q^{2}$) and its position with the model's predictions of 3.06 and 9.12, respectively; also measure the K0 charge radius to check the predicted negative mean square radius near -0.090 $fm^{2}$.","tokens_in":8341,"feed_emoji":"⚛️","tokens_out":4299,"duration_ms":46950,"temperature":0.7,"pith_summary":"The paper tries to establish that one relatively simple light-front quark model, with Gaussian wavefunctions and parameters fixed by mesonic decay constants, can describe the electromagnetic size and charge shape of every pseudoscalar meson. It reproduces the measured pion and kaon form factors at low momentum transfer and their charge radii, and it predicts systematic trends: heavier quark content shrinks the charge radius, while the constituent quark mass difference sets the height of the form-factor peak. If the picture is right, the internal charge structure of heavy D and B mesons follows a clean, mass-asymmetry-dominated pattern that upcoming experiments can test.","feed_headline":"One quark model reproduces pion and kaon form factors","feed_subtitle":"Gaussian wavefunctions fitted to decay constants predict how quark mass gaps shrink charge radii across heavy mesons.","key_machinery":"The central object is the light-front quark model expression for the electromagnetic form factor, F_P($Q^{2}$) = e_q I(m_q, m_bar{q}, $Q^{2}$) + e_bar{q} I(m_bar{q}, m_q, $Q^{2}$), where I is an overlap integral over the longitudinal momentum fraction x and transverse momentum k_perp of a Gaussian radial wavefunction. The quark mass asymmetry enters through the invariant mass M_0 and the Jacobian of the variable transformation, and the Gaussian parameter $\\beta$ controls the bound-state size. This machinery converts quark masses, charges, and one width parameter into the full $Q^{2}$ dependence of every pseudoscalar meson's form factor and charge radius.","core_discovery":"The authors claim that the light-front quark model yields electromagnetic form factors whose charged- and neutral-meson endpoint behaviors differ significantly, because the slope at $Q^{2}$ = 0 encodes the spatial separation of quark and antiquark charges; for neutral mesons the outer, lighter quark's charge dominates, so the mean square charge radius becomes negative. They further claim that for charged mesons the peak value of $Q^{2}$ F_P($Q^{2}$) is approximately proportional to the mass difference between constituent quarks, that the charge radius decreases as the meson mass increases, and that for neutral mesons the heavy quark's electric charge primarily determines the radius. These trends are presented as predictions for heavy-flavor mesons, with the pion and kaon results serving as validation.","pith_inferences":["Inference: If the mass-asymmetry pattern holds, the charge radius of any unmeasured pseudoscalar meson could be estimated from its constituent quark masses and charges alone, without a full dynamical calculation.","Inference: A precise measurement of the K0 charge radius would be a sharp test of the charge-separation mechanism, since the model ties its negative sign to the unequal spatial extent of the d and s quark charge distributions.","Inference: The Gaussian wavefunction ansatz may be too rigid to describe the very high-Q^2 tail; the early onset of asymptotic behavior could be an artifact of that ansatz rather than a physical prediction, so lattice results at higher Q^2 would clarify this.","Inference: The same formalism could be extended to other current components or to include zero-mode contributions, which would quantify how much of the predicted endpoint behavior depends on the one-body valence approximation."],"forward_implications":["The pion and kaon form factors in the low-Q^2 region are reproduced, so the model can serve as a benchmark for extracting charge radii from future low-Q^2 experiments.","The predicted form factor shapes for D and B mesons give concrete targets for JLab and EIC measurements at moderate and high momentum transfer.","The neutral kaon's negative mean square charge radius is reproduced, supporting a picture where the lighter, oppositely charged quark sits farther from the center than the heavy quark.","Heavy-meson charge radii decrease with meson mass, with B_c predicted to have the smallest radius among the states considered, a trend that can be checked once heavy-meson radii are measured.","The model reaches the perturbative 1/Q^2 asymptotic behavior earlier than current lattice results, a difference that future high-Q^2 data can discriminate."],"supporting_citations":[{"why":"Supplies the low-Q^2 NA7 pion form factor data against which the model's endpoint behavior is validated.","marker":"[5]"},{"why":"Supplies lattice QCD results for pion and kaon form factors used as intermediate- and high-Q^2 comparisons.","marker":"[22]"},{"why":"Supplies Bethe-Salpeter charge radius predictions that the model's heavy-meson radii are compared with.","marker":"[31]"},{"why":"Provides the light-front quark model formalism from which the form factor overlap integral is taken.","marker":"[60, 61]"},{"why":"Supplies the PDG mesonic decay constants used to fit the quark masses and Gaussian width parameters, and the experimental charge radii for pions and kaons.","marker":"[66]"},{"why":"Supplies the perturbative QCD prediction that F_P(Q^2) falls as 1/Q^2 at high Q^2, used to judge the model's asymptotic behavior.","marker":"[67]"}],"fun_headline_variants":["Meson form factors follow quark mass gaps","Light-front model predicts heavy meson charge radii","Quark mass asymmetry shapes meson charge radii","Pion, kaon form factors fit light-front quark model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the plus-component electromagnetic current in the q^+ = 0 frame is saturated by the one-body valence diagram, so any omitted zero-mode or pair-creation contribution would change the predicted form factors and charge radii.","fun_headline_variants_meta":{"raw":{"variants":["Meson form factors follow quark mass gaps","Light-front model predicts heavy meson charge radii","Quark mass asymmetry shapes meson charge radii","Pion, kaon form factors fit light-front quark model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000241,"raw_usage":{"total_tokens":1523,"prompt_tokens":951,"completion_tokens":572,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":511}},"tokens_in":567,"tokens_out":572,"duration_ms":6658,"temperature":1.0,"reasoning_tokens":511,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:38:36.712240+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the D+ or B+ electromagnetic form factor at JLab or the EIC up to $Q^{2}$ of roughly 10 $GeV^{2}$ and compare the peak of $Q^{2}$ F_P($Q^{2}$) and its position with the model's predictions of 3.06 and 9.12, respectively; also measure the K0 charge radius to check the predicted negative mean square radius near -0.090 $fm^{2}$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the low-Q^2 NA7 pion form factor data against which the model's endpoint behavior is validated."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies lattice QCD results for pion and kaon form factors used as intermediate- and high-Q^2 comparisons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Bethe-Salpeter charge radius predictions that the model's heavy-meson radii are compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the PDG mesonic decay constants used to fit the quark masses and Gaussian width parameters, and the experimental charge radii for pions and kaons."}],"review_version":1}