{"id":"b00e6aef-a55e-4c85-8f1c-fc0470b4c02d","arxiv_id":"2507.06025","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a nonlocal quark model, the pion gravitational form factors are calculated, with bubble diagrams contributing about 8% of A(0) and 69% of D(0).","lead":"The paper computes the gravitational form factors of the pion in a nonlocal quark model, finding the D-term close to the predicted -1 and showing that bubble diagrams contribute strongly. It matters because it estimates how quark and gluon-like contributions to the pion's mass and pressure distribution enter in this model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central 69%-bubble claim is not gauge invariant and is even flagged in the paper as not well-defined; absent a scheme-independent definition, it is a representation artifact.","rationale":"The paper's most salient new observation is the 8%/69% decomposition, but that decomposition is not a gauge-invariant or representation-independent quantity. The paper states this limitation itself in Section IV: the separate bubble contribution is not well-defined because only the sum of all diagrams gives a conserving energy-momentum tensor. This is not a disagreement with an external convention; it is an internal property of the calculation. The reader's weakest_assumption focused on the survival of the nonlocal vertex shape after removal of the background field B, which is a legitimate model-input concern about the total amplitude. However, even if that assumption is granted, the headline percentage is still ambiguous because the split between vertex-expansion and propagator-expansion diagrams depends on how the equivalence principle is implemented in the nonlocal vertex. The reader's rationale also mentions that the diagram decomposition is not gauge invariant, so there is partial agreement, but the reader's formal weakest_assumption is not the same as this representation-dependence. A concrete re-computation under an alternative covariantization would settle whether the 69% number is stable; if it is not, the claim should be withdrawn or explicitly reframed as scheme-dependent. Because the original reader already assigned CONDITIONAL and the total form-factor calculation remains a valid conditional model result, the verdict does not need to change.","tokens_in":9357,"tokens_out":4630,"duration_ms":56059,"concrete_test":"Recompute the isolated (b)+(c) contribution to D(0) using an alternative but equivalent covariantization of the nonlocal vertex: rewrite the exponent entering Eq. (9) with a different ordering of the covariant derivatives (e.g., symmetrize before exponentiating, or integrate one derivative by parts before varying the tetrad), while keeping the total action and the total D(0) unchanged. If the fraction of D(0) coming from diagrams (b)+(c) shifts by more than about 10 percentage points from 69%, the claimed bubble contribution is representation-dependent and should not be reported as a physical result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's abstract highlights that bubble diagrams contribute approximately 8% of A(0) and 69% of D(0). Yet this decomposition into bubble diagrams (Fig. 2b,c) and impulse diagrams (Fig. 2d,e) is not a gauge-invariant or representation-independent quantity. The paper itself concedes in Section IV that their separate contribution is not well-defined, because only the sum of all diagrams yields a conserving energy-momentum tensor. The split is fixed by how the minimal substitution in Eq. (7) is distributed among the two derivative operators inside the nonlocal vertex in Eq. (9). Different but physically equivalent implementations—for example, reordering the exponentiated covariant derivatives in Eq. (5), moving spin-connection terms by integration by parts, or choosing a different symmetric ordering before covariantizing—will move contributions between the vertex-expansion diagrams (b,c) and the propagator-expansion diagrams (d,e), leaving the total ⟨Tμν⟩ unchanged. Since only the total is fixed by the Ward identity, the 8%/69% figures are properties of a particular diagrammatic bookkeeping rather than robust model predictions. Even granting the model entirely, the central numerical claim has no objective meaning unless a scheme-independent definition of 'bubble contribution' is supplied.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the gravitational form factors A(t) and D(t) of the pion in a nonlocal quark model. Gravity is introduced by promoting ordinary derivatives to covariant derivatives in the nonlocal quark-meson vertex via the equivalence principle, and the energy-momentum tensor is defined as the vielbein variation of the effective action. The calculation is performed in a one-loop approximation and organized into the diagrams of Fig. 2: a volume-element diagram, two diagrams from the nonlocal vertices, and two impulse-approximation diagrams from the quark propagator. The author reports A(0) close to 1, D(0) close to -1 consistent with soft-pion theorems, a pion mass radius of 1.12 fm, and a substantial contribution from the bubble diagrams (approximately 8% of A(0) and 69% of D(0)).","tokens_in":9626,"tokens_out":4137,"duration_ms":47662,"significance":"If the result is robust, the paper introduces a workable method for coupling nonlocal quark models to gravity and provides an explicit model estimate of gluonic-type contributions to the pion energy-momentum tensor that are absent in the usual impulse approximation. The derivation is explicit, the momentum integrals are finite because of the vertex falloff, and the parameters are fixed from the pion mass and the pi->gamma gamma width rather than from the form factors themselves, which is a genuine strength. The prediction D(0) approximately -1 is also a nontrivial check of the model. However, the headline quantitative claim about the bubble diagram contribution is not gauge invariant and is in fact self-described in the paper as not a well-defined quantity, so the central interpretive claim needs either a scheme-independent definition or substantial reframing.","major_comments":[{"comment":"The headline claim that bubble diagrams contribute approximately 8% of A(0) and 69% of D(0) is not gauge invariant. The split into diagrams (b),(c) and (d),(e) depends on how the minimal substitution in Eq. (7) is distributed among the derivative operators inside the nonlocal vertex of Eq. (9), and in particular on the ordering of the exponentiated covariant derivatives before expansion. A physically equivalent reordering or integration by parts moves contributions between the vertex-expansion and propagator-expansion diagrams without changing the total energy-momentum tensor. The paper itself concedes in the same paragraph that the separate contribution 'is not a well-defined quantity.' Since these percentages are the main advertised result, the author must either supply a scheme-independent definition of the bubble contribution or remove/reframe this claim.","section":"Section IV, paragraph beginning 'An important observation'; Eqs. (7), (9), (15)"},{"comment":"The Ward-identity consistency check A(0)=1 is not explicitly quantified. The figure suggests A(0) is close to 1, but no numerical value, numerical tolerance, or error estimate is given for A(0), D(0), or the derived radius sqrt(<r^2>)=1.12 fm. Because A(0)-1 is a direct measure of the numerical and truncation accuracy of the calculation, please state the computed A(0) and D(0) and the precision of the Gaussian/proper-time integrals used to obtain them.","section":"Section IV and Fig. 3; Eq. (6)"},{"comment":"The relative weight of the bubble and impulse contributions is controlled by the assumption that the nonlocal vertex shape and the parameter Lambda, originally derived from the gluon propagator in a homogeneous self-dual background field, remain unchanged when the background is removed and gravity is coupled through minimal substitution. This is the central model assumption and it is not tested. A concrete test would be to compute the pion electromagnetic form factor in the same framework and compare with experiment or lattice data, or to vary the ordering/covariantization prescription and show that the bubble fraction changes; without such a test, the 69% figure should be presented as scheme-dependent.","section":"Section II and III, Eqs. (4), (5), (9)"}],"minor_comments":[{"comment":"The sentence 'only the sum of all diagrams in Fig. 3 yields conserving energy-momentum tensor' should refer to Fig. 2, because Fig. 3 shows the form factors rather than the diagrams.","section":"Section IV, paragraph on the separate contribution"},{"comment":"In the integrands for diagrams (b), (c), (d), and (e), the second integration variable is written as d4x in both factors; it should be d4x' to match the momentum labels exp(-ipx) and exp(ip'x').","section":"Appendix B, equations for diagrams (b)-(e)"},{"comment":"The notation 'JaJ' in the definition of the current J_Q is unexplained and appears to involve a typographical duplication; please clarify the index structure.","section":"Eq. (3)"},{"comment":"The mean square radius formula is presented without derivation; adding a brief explanation of the terms, especially the role of the 1/(4M_pi^2) correction, would improve readability.","section":"Section IV, radius formula"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for hep-ph and represents a serious calculation, not a superficial model exercise. The main problem is interpretive: the total form factors are well-defined and the D(0) result is interesting, but the bubble-diagram decomposition that the abstract emphasizes is not a gauge-invariant observable and the author's own caveat in Section IV already concedes this. This is fixable by reframing the paper's claims or by defining an invariant measure of the bubble contribution, so I do not recommend rejection. I also note that the paper would be strengthened by reporting numerical values with uncertainties for A(0), D(0), and the radius."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a competent model calculation that does something new – it applies the gauging of nonlocal Lagrangians to gravity and computes the pion's A(t) and D(t). The total form factors look plausible, and D(0) ≈ −1 is a nice check of the soft pion theorem. But the paper's lead claim, that bubble diagrams contribute 8% of A(0) and 69% of D(0), is not a gauge-invariant statement. The paper itself admits in Section IV that the split is not well-defined because only the sum of diagrams yields a conserved energy-momentum tensor. The stress-test note is right: you can move contributions between bubble and impulse diagrams by reordering the covariantized derivatives or choosing a different symmetric ordering. So those percentages are bookkeeping artifacts, not model predictions.\n\nWhat's genuinely useful is the total A(t), D(t), and the mechanical radius sqrt(<r²>) = 1.12 fm. The derivation is careful and the use of the equivalence principle is standard. The nonlocal vertex shape is imported from the background-field model, which is an assumption, but a reasonable one for this class of models. The numerical evaluation is not fully reproducible – no error bars, no code, and A(0)=1 consistency is not shown explicitly – but that's typical for this kind of paper.\n\nThe main problem is framing. The abstract sells the bubble contribution as the finding, but the paper's own caveat undercuts it. A referee should ask the author to either give a scheme-independent definition of the bubble part (unlikely) or to reframe the paper around the total form factors and the model's prediction for the mass radius, presenting the bubble decomposition as an illustrative diagnostic. The D(0) result and the full t-dependence are worth having on the record.\n\nBottom line: it deserves a serious referee, but the central numerical claim as stated should not be taken at face value. If you cite it, cite the total form factors, not the 69%.","headline":"Competent calculation of pion gravitational form factors in a nonlocal model, but the headline bubble-diagram percentages are not gauge-invariant and should not be taken as physical.","tokens_in":10121,"tokens_out":2426,"would_cite":true,"duration_ms":25431,"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":"In a nonlocal quark model, coupling gravity through the equivalence principle generates bubble diagrams from the nonlocal pion vertex that supply about 8% of $A(0)$ and 69% of $D(0)$, with $D(0)$ ending up near $-1$.","keywords":["gravitational form factors","pion","nonlocal quark model","energy-momentum tensor","D-term","bubble diagrams","equivalence principle","soft pion theorems"],"falsifier":"Recompute diagrams (b) and (c) with a different nonlocal vertex profile, for example replacing the exponential $1/p^2$ regulator by a Gaussian regulator, re-fitting $m$ and $\\Lambda$ to the pion mass and $\\pi\\to\\gamma\\gamma$ width; if the bubble share of $D(0)$ moves away from 69% by more than a few percent, the result is tied to the specific vertex shape rather than to the equivalence-principle coupling.","tokens_in":9162,"feed_emoji":"⚛️","tokens_out":7706,"duration_ms":80553,"temperature":0.7,"pith_summary":"The paper aims to show that a consistently coupled gravitational field in a nonlocal quark model generates a specific set of gravitational form factors for the pion, and that the nonlocal vertex—the piece representing gluon effects—contributes much more than the usual quark impulse diagrams for the D-term. In particular, bubble diagrams are found to give about 8% of $A(0)$ and 69% of $D(0)$, while the total D-term remains near $-1$ as required by soft-pion theorems. If this is right, any model estimate of the pion's gravitational form factors that stops at the impulse approximation misses the dominant mechanism behind the D-term. The paper also extracts a gravitational radius of 1.12 fm, larger than the pion charge radius, suggesting the mass distribution is more extended than the charge distribution.","feed_headline":"Bubble diagrams carry 69% of the pion's gravitational D-term","feed_subtitle":"Coupling gravity to nonlocal quark vertices adds a large gluonic contribution; D(0) stays near -1.","key_machinery":"The load-bearing object is the nonlocal quark-pion vertex $$V_a = i\\gamma_5 M^a \\$int_0^{1}$ dt \\exp\\!\\left(\\frac{\\partial_\\$leftrightarrow^{2}$}{\\$Lambda^{2}$}t\\right), \\quad \\partial_\\leftrightarrow = \\xi \\partial_\\leftarrow - \\xi' \\partial_\\rightarrow ,$$ whose shape and scale come from the gluon propagator in a background self-dual field. Gravity is inserted by minimal substitution into this vertex and the quark kinetic term, using tetrads and spin connection; varying the effective action with respect to the tetrad gives the Hilbert energy-momentum tensor. Expansion of the nonlocal vertex in the metric fluctuation produces diagrams (b) and (c) — the \"bubble\" terms — alongside the usual impulse diagrams from the quark propagator. The vertex's $1/p^2$ ultraviolet behavior keeps all momentum integrals finite, and Schwinger proper-time representations turn them into Gaussian integrals that are evaluated analytically before numerical proper-time integration.","core_discovery":"Starting from the generating functional of a nonlocal quark model for the pion, the paper couples an external gravitational field through the equivalence principle, then defines the flat-space energy-momentum tensor by varying the effective action with respect to the tetrad. In the one-loop approximation the matrix element of $T_{\\mu\\nu}$ contains five diagram classes; besides the standard impulse diagrams, the expansion of the nonlocal vertex yields bubble diagrams. The paper's central numerical finding is that these bubble diagrams contribute roughly 8% of $A(0)$ and about 69% of $D(0)$, while the full $D(0)$ comes out near $-1$, consistent with soft-pion theorems. The corresponding gravitational mean-square radius of the pion is about 1.12 fm, larger than the charge radius of 0.659 fm, indicating that mass is distributed less compactly than charge in this model.","pith_inferences":["One could test the robustness of the 69% figure by repeating the calculation with a different nonlocal vertex profile, such as a Gaussian regulator fitted to the same pion observables; a large change would show the result is sensitive to the specific vertex shape rather than the equivalence-principle mechanism itself.","The same minimal-substitution construction should apply to other hadrons, where an analogous bubble contribution would modify the gluonic part of the energy-momentum tensor; lattice results for nucleon gravitational form factors would be a natural external check.","If the bubble contribution dominates the D-term at low $t$, then extractions of the pion D-term from generalized parton distributions will need to include gluonic and nonlocal-vertex contributions even at small momentum transfer.","The near equality of $D(0)$ to $-1$ at the physical pion mass suggests the model's chiral-limit corrections are small; computing $D(0)$ as a function of $M_\\pi$ would give a concrete estimate of those corrections."],"forward_implications":["The pion's D-term in this model is controlled by bubble contributions from the nonlocal vertex, so impulse-approximation-only calculations in separable or nonlocal quark models systematically under-predict the gluonic part of $T_{\\mu\\nu}$.","The model reproduces $D(0)\\approx -1$ without tuning, which supports the soft-pion theorem and suggests the bubble terms are the finite-mass corrections that bring the D-term from the chiral-limit value to the physical one.","The gravitational radius of the pion (1.12 fm) exceeds its charge radius (0.659 fm) in this model, so the spatial distribution of energy is predicted to be broader than the distribution of charge.","Because the Hilbert definition yields a symmetric, conserved energy-momentum tensor, the Ward identity $A(0)=1$ is satisfied by the full set of diagrams, which would not happen if only the quark-current part were kept.","Gluonic (bubble) contributions to gravitational form factors are much larger than analogous nonlocal-vertex corrections to electroweak form factors, since gravity couples directly to the gluonic content encoded in the vertex."],"supporting_citations":[{"why":"Establishes the gravitational Ward identity that fixes $A(0)=1$ and forces $\\bar{c}(t)=0$, the constraint the model must satisfy.","marker":"[7]"},{"why":"Soft-pion theorems predicting $D(0)=-1$, the benchmark the numerical D-term is checked against.","marker":"[8–10]"},{"why":"Lattice QCD computation of pion gravitational form factors that provides the existing nonperturbative comparison.","marker":"[13]"},{"why":"Rainbow-ladder Bethe–Salpeter calculation of pion electromagnetic and gravitational form factors, the impulse-approximation result to which bubble diagrams are added.","marker":"[14]"},{"why":"Show that consistent gauging of nonlocal or bound-state Lagrangians generates interaction and bubble currents, the basis for diagrams (b) and (c).","marker":"[15–18]"},{"why":"Derives the nonlocal quark current and vertex shape from hadronization in a background self-dual gluon field, the origin of the model.","marker":"[19]"},{"why":"Determines the scale $\\Lambda$ and Regge-spectrum phenomenology that fix the model parameters used here.","marker":"[22]"},{"why":"Relates the nonlocal currents to Bethe–Salpeter amplitudes, supplying the operator basis for the vertex.","marker":"[24]"},{"why":"Particle Data Group values for the pion mass and $\\pi\\to\\gamma\\gamma$ width used to fix the quark mass and coupling $h$.","marker":"[29]"}],"fun_headline_variants":["Bubble diagrams dominate pion's gravitational D-term at 69%","Nonlocal quark model: bubbles drive pion's gravitational form factor","Pion's gravitational radius exceeds charge radius in nonlocal model","Gravity meets quarks: bubble loops shape pion's mass distribution","Equivalence principle in quark model yields pion's gravitational shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The nonlocal vertex shape taken from the gluon propagator in a background self-dual field is assumed to stay exactly the same once that background is removed and gravity is coupled by minimal substitution, and that fixed shape and the fixed scale $\\Lambda$ determine the relative weight of the bubble and impulse diagrams.","fun_headline_variants_meta":{"raw":{"variants":["Bubble diagrams dominate pion's gravitational D-term at 69%","Nonlocal quark model: bubbles drive pion's gravitational form factor","Pion's gravitational radius exceeds charge radius in nonlocal model","Gravity meets quarks: bubble loops shape pion's mass distribution","Equivalence principle in quark model yields pion's gravitational shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000195,"raw_usage":{"total_tokens":1262,"prompt_tokens":758,"completion_tokens":504,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":374,"completion_tokens_details":{"reasoning_tokens":417}},"tokens_in":374,"tokens_out":504,"duration_ms":5126,"temperature":1.0,"reasoning_tokens":417,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:11:53.095089+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute diagrams (b) and (c) with a different nonlocal vertex profile, for example replacing the exponential $1/p^2$ regulator by a Gaussian regulator, re-fitting $m$ and $\\Lambda$ to the pion mass and $\\pi\\to\\gamma\\gamma$ width; if the bubble share of $D(0)$ moves away from 69% by more than a few percent, the result is tied to the specific vertex shape rather than to the equivalence-principle coupling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the gravitational Ward identity that fixes $A(0)=1$ and forces $\\bar{c}(t)=0$, the constraint the model must satisfy."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the nonlocal quark current and vertex shape from hadronization in a background self-dual gluon field, the origin of the model."},{"cited_title":"Meson masses within the model of induced nonlocal quark currents","cited_arxiv_id":"hep-ph/9601344","evidence_quote":"Determines the scale $\\Lambda$ and Regge-spectrum phenomenology that fix the model parameters used here."},{"cited_title":"Regge spectra of excited mesons, harmonic confinement and QCD vacuum structure","cited_arxiv_id":"1603.01447","evidence_quote":"Relates the nonlocal currents to Bethe–Salpeter amplitudes, supplying the operator basis for the vertex."},{"cited_title":"Dipole polarizabilities of light pseudoscalar mesons within the Domain Model of QCD vacuum","cited_arxiv_id":"2208.00253","evidence_quote":"Particle Data Group values for the pion mass and $\\pi\\to\\gamma\\gamma$ width used to fix the quark mass and coupling $h$."}],"review_version":1}