Pion and kaon electromagnetic and gravitational form factors

A unified set of predictions for pion and kaon elastic electromagnetic and gravitational form factors is obtained using a symmetry-preserving truncation of each relevant quantum field equation. A key part of the study is a description of salient aspects of the dressed graviton + quark vertices. The calculations reveal that each meson's mass radius is smaller than its charge radius, matching available empirical inferences; and meson core pressures are commensurate with those in neutron stars. The analysis described herein paves the way for a direct calculation of nucleon gravitational form factors.


Introduction
A new era is heralded by construction and planning of high-luminosity, high-energy facilities [1][2][3][4][5][6][7][8][9][10].Following a one-hundred year focus on the proton, such facilities will enable exploration of the structure of many other states produced by the strong interaction Hamiltonian.Principal amongst them are pions (π) and kaons (K), viz.Nature's most fundamental (would-be) Nambu-Goldstone bosons.These states are special because, absent Higgs boson (HB) couplings into quantum chromodynamics (QCD), they are all massless and identical.The Higgs inserts quark current-masses into the QCD Hamiltonian; yet at realistic current-mass values, π-and K mesons remain unnaturally light, with masses that seemingly belie their categorisation as hadrons, e.g., the π mass is similar to that of the µ-lepton, despite pions being strong-interaction bound-states.
Measurements relating to π and K generalised parton distributions (GPDs) are also likely [3-5, 7, 10].Such data would be significant because GPDs provide access to hadron gravitational form factors [36]; so, could open doors to comparisons between the electromagnetic and gravitational structure of π-and K-mesons.This would enable insights to be drawn into the impacts of Nature's two known mass generating mechanisms on the structure of Nambu-Goldstone bosons.It is thus imperative for theory to deliver sound, unifying predictions for π and K electromagnetic and gravitational form factors.
Steps toward a unified set of Poincaré invariant statements about the distributions of charge, mass, and pressure inside pions and kaons are described in Refs.[25,26].Those analyses exploit the GPD overlap representation, with light-front wave functions constrained by π and K valence-quark DFs.Profiting from such studies, Ref. [37] delivered a data driven prediction for the pion mass distribution, showing that the pion's mass radius is ≈ 20% smaller than its charge radius.
Herein, we adopt a less phenomenological approach.Namely, using continuum Schwinger function methods (CSMs) [38,39], we solve the probe+meson scattering problem at leading-order of the truncation scheme introduced in Refs.[40,41].Therewith, we also unify π and K electromagnetic and gravitational distributions with predictions for an array of hadron structural properties that range over systems with up-to three heavy quarks, e.g., Refs.[42][43][44][45][46].
The presentation is arranged as follows.Section 2 introduces an approximation to the probe+meson scattering problem.It is followed, in Sec. 3, by a discussion of dressed probe+quark vertices and their properties.Section 4 details the quark+antiquark scattering kernel which underlies all concrete calculations herein.Numerical results for the dressed probe+quark vertices are described in Sec. 5.This is followed by an explanation of algebraic Ansätze for all elements of the calculation, which are subsequently exploited to provide ultraviolet completions of the results.Predictions for pion and kaon elastic electromagnetic and gravitational form factors are discussed in Sec. 7. The gravitational form factors are used to calculate the pressure and shear force profiles drawn in Sec. 8. Section 9 provides a summary and perspective.

Form factors
Consider a charged pion, π + , built from a u valence quark and a d valence antiquark and suppose that isospin symmetry is exact, so the only difference between these degrees-of-freedom is their electric charge.Then, in rainbow-ladder (RL) truncation [40,41], the five-point Schwinger function that defines the π + (p) → π + (p ′ ) elastic electromagnetic form factor takes the form drawn in the top diagram of Fig. 1.This image translates into L(l, P, Q) = S(l + p)Γ π (l + p/2; p)S(l) where ] is the dressed-quark propagator, calculated in rainbow truncation; Γ π is the RLtruncation pion Bethe-Salpeter amplitude; and Γ γ ν is the RL dressed photon+quark vertex.(For an overview, see, e.g., Refs.[38,39].)Since the gluons that bind the π do not carry electric charge, then, in RL truncation, insofar as Λ γπ ν is concerned, the bottom diagram in Fig. 1 is zero.
Regarding Eq. ( 1), electromagnetic current conservation is guaranteed in RL truncation [47].Thus, Q ν Λ γπ ν (P, Q) ≡ 0; so, and F π (Q 2 ) can be extracted from the projection The expectation value of the energy-momentum tensor in the pion, viz. the π gravitational current, takes the following form: where θ π 2,1 are, respectively, the in-pion mass and pressure distribution form factors, which are also accessible via the pion GPD [36,48].The following relations follow from symmetries: The first identity is a statement of mass normalisation, like F π (Q 2 = 0) = 1 for the electromagnetic form factor; the second is a corollary of EHM, expressed in a soft-pion theorem [49,50]; and the third is a basic consequence of energy-momentum conservation, viz.
In RL truncation, both diagrams in Fig. 1 contribute to the pion gravitational current.The bottom diagram plays a role analogous to Ref. [51, Fig. 3B ′ ], which restores momentum conservation in calculations of pion DFs.Notwithstanding that, it is not necessary herein to develop an explicit expression for this term because it only affects cπ (Q 2 ): Hence, regarding θ 2,1 (Q 2 ), it is sufficient to consider only the upper diagram.This translates into an expression like Eq. ( 1), with the only change being that the photon + quark vertex is replaced by the graviton + quark analogue: Γ γ ν (l + p ′ , l + p) → Γ g µν (l + p ′ , l + p).The gravitational form factors may then be extracted via contractions of the upper/triangle-diagram component of Λ g µν (P, Q) with the following projectors (L µν (P ) = P µ P ν /P 2 ): When considering K form factors, the only difference from the pion is that one must distinguish quark flavours when expressing the currents illustrated by Fig. 1 -see, e.g., Refs.[19,52] for comparative studies of π and K electromagnetic form factors.In the isospin symmetry limit, for any given form factor, F K : F K = F s K +F u K .(For electromagnetic form factors, one must include multiplicative quark electric charge factors.) We work at the hadron scale, whereat quasiparticle degrees of freedom carry all properties of a given hadron.Evolution to higher scales [53][54][55][56], which exposes QCD parton contributions in species-decompositions of hadron structural properties [26,29], is discussed elsewhere [27,57,58].Apart from ensuring the correct anomalous dimensions [42,59], such evolution has no effect on the overall meson form factors, which are our focus.

Symmetry-preserving vertices
A key factor in calculations of probe+meson form factors is the dressed probe+quark vertex.Following thirty years of study [47,[60][61][62], the photon+quark vertex is well understood.It satisfies a Ward-Green-Takahashi (WGT) identity (l Thus, four terms in Γ γ ν are fixed by the dressed-quark propagator [63][64][65], that associated with σ µν l (′) +µ l +ν being zero.This is called the Ball-Chiu (BC) part of the vertex.Eight independent nonzero terms remain.They can be obtained by solving the inhomogeneous vectorchannel Bethe-Salpeter equation using the same RL formulation as employed for S, Γ π in Eq. (1).
Combining this information, one can decompose the vertex as follows: Here, [47, Eq. ( 26)]: where is a matrix-valued regular function (free of kinematic singularities for real arguments) and the second term in Eq. ( 7) exhibits, inter alia, timelike-Q 2 poles, one at the mass-squared of each neutral vector meson supported by the interaction.The lightest such state is the ρ 0 -meson.Equation ( 7) is neither an approximation nor simplification.Instead, it merely introduces a physical understanding into the calculation of Γ γ ν .The new element herein is the dressed graviton + quark vertex, Γ g µν , which has not hitherto been studied using a realistic gluon + quark interaction.It is first worth recalling that Γ g µν satisfies its own WGT identity, viz. a straightforward generalisation of that for scalar fields in Ref. [66]: Next, recall that Γ g µν is an isoscalar (I = 0) J P C = 2 ++ vertex.Hence, by analogy with Γ γ ν , Γ g µν must exhibit timelike-Q 2 poles (or resonance-like structures) at the mass-squared positions of each I = 0 meson for which a 2 ++ tensor coupling can be constructed.Naturally, this includes the tower of isoscalar tensor mesons, but that does not exhaust the range of possibilities.Poles associated with isoscalar-scalar mesons also emerge -see, e.g., Refs.[67][68][69].The lowest-mass singularity in Γ g µν (k, Q) occurs in this channel.It is associated with the tensor structure Q 2 T µν (Q) [67] and therefore -see Eq. ( 3) -contributes significantly to θ π 1 (Q 2 ) without much affecting θ π 2 (Q 2 ).The trajectory of I = 0 axialvector mesons can also contribute to Γ g µν (k, Q).In this case, the available tensor structures are more complex; e.g., the following form is admissible: 1,2 are associated with {µν}-symmetric tensors in the graviton + meson vertex and this term is antisymmetric, then axialvector mesons are not relevant in the present analysis.
Capitalising on these observations, the following minimal Ansatz is evidently admissible as a resolution of Eq. ( 9): where Γ 2 αβ (k; Q) is the tensor + quark vertex generated by the inhomogeneity ) so as to ensure the absence of kinematic singularities.Notwithstanding this, dynamical singularities do appear in Γ 2 µν ; namely, one at the pole position of each I = 0 tensor meson.Given the structure of P θ2 µν in Eq. (5a), one should expect the lowest-mass singularity in this vertex to affect θ π 2 .The next pole lies much deeper in the timelike region, so must play a much lesser role.Since Γ 2 µν possesses eight independent Dirac matrix valued tensor structures, then this Ansatz involves fourteen such nonzero terms.
Of course, Eq. ( 10) is not complete.Considering the available four-vectors and tensor structures, and recognising that fermion on-shell conditions cannot be used, Γ g µν may contain a large number of additional independent terms.Nevertheless, in general, one may write where represents all possible transverse structures not already included.Plainly, like Γ 2 0µν , Γ gT µν does not contribute to resolving Eq. ( 9).It may be determined by solving the appropriate Bethe-Salpeter equation.In doing so, one sees the emergence of isoscalar scalar mesons in the graviton+quark vertex.
Again, Eq. ( 12) is not an approximation: it just introduces a physical understanding into the calculation of Γ g µν .

Dressed vertices
Having specified the Bethe-Salpeter kernel, every element in Eq. ( 1) and its analogue for the graviton+meson current, Eq. ( 3), can be computed.For instance, the photon+quark vertex is the solution of an inhomogeneous Bethe-Salpeter equation: where Λ dl represents a symmetry-preserving regularisation of the four-dimensional integral, with Λ being the regularisation scale; S is the propagator of the interacting quark, obtained from the analogous rainbow gap equation; and Z 2 (ζ, Λ) is the quark wave function renormalisation constant, determined as part of solving that gap equation.
The solution for Γ γ ν has eleven independent Dirac matrix valued terms and can readily be projected into the form of Eq. (7).Good numerical methods for solving sets of coupled gap and Bethe-Salpeter equations are described, e.g., in Refs.[75,76].
The remaining element in Eq. ( 3) is the bound-state amplitude of the meson being probed, which can be obtained from an analogous homogeneous Bethe-Salpeter equation [75].
Consider now the graviton+quark vertex, the RL result for which may be obtained by solving where Z 0 m is the chiral-limit mass renormalisation constant and the WGT identity, Eq. ( 9), ensures that it is Z 2 which appears here.
As already noted, in general, Γ g µν possesses a large number of independent terms.On the other hand, we have also emphasised that, physically, one may expect just a few contributions to be important, viz.those parts which saturate the WGT identity, and pieces associated with an f 2 tensor meson pole and an analogous scalar meson resonance.We therefore continue with a truncated form of Γ g µν .Specifically, the vertex given by Eq. ( 12), with Γ 2 µν obtained by solving the tensor analogue of Eq. ( 16) defined via the Eq. ( 11) inhomogeneity, and where Γ I (k; Q), which has four independent Dirac matrix valued structures, D j=1,4 Here, Equation ( 19) is an inhomogeneous Bethe-Salpeter equation.Its solution exhibits a pole at the mass of each scalar meson generated by the interaction in Eq. (13).In practice and analogous to the tensor vertex, the lightest scalar dominates in the calculation of spacelike gravitational form factors because the first excitation lies roughly 1 GeV higher in mass [73], i.e., well into the timelike region.
Solving Eq. ( 17) and the tensor analogue of Eq. ( 16), then the lowest mass isoscalar-scalar and -tensor mesons are expressed in Γ g µν as pole terms of the form: where Γ S,T are their associated canonically normalised bound-state amplitudes [75,78].The calculated masses and residue coefficients are (in GeV): Here we have listed results for the probe+light-quark and probe+s-quark vertices.The latter are required in the calculation of kaon form factors.The photon+quark vertices, obtained from Eq. ( 16), display analogous features with, e.g., pole contributions from neutral vector mesons at timelike momenta [60][61][62].
All elements required for calculation of π and K electromagnetic and gravitational form factors are now available: dressed quark propagators (2 one-variable scalar functions for each quark), meson Bethe-Salpeter amplitudes (4 two-variable functions for each meson), and probe+quark vertices (2 one-variable and 8 two-variable functions for the electromagnetic interaction of each quark and 2 one-variable and 12 two-variable functions for the graviton+quark interaction of each quark).The scalar functions are stored as arrays of numbers, wherewith the integrand in Eq. ( 1) -and its analogues for other systems and probes -can be formed and the associated integral evaluated using standard quadrature and interpolation schemes.

Algebraic Ansätze
We choose to complement our fully numerical work with calculations that employ minimal effective algebraic inputs for each of the functions.Such representations of quark propagators and meson Bethe-Salpeter amplitudes were used in Refs.[25,26] to deliver results for π and K GPDs and, therefrom, electromagnetic and gravitational form factors for these mesons.Therein, Our variant of the model defining parameters is listed in Table 1: the dressed-quark masses are somewhat larger than those in Refs.[26, Table I] and Λ K is smaller.The other values are unchanged.(n P =π,K are computed normalisation constants, which ensure F P (0) = 1 = θ P 2 (0).)Using Eq. (23b), then Eq. ( 8) simplifies to Γ BC ν = γ ν ; hence, we use Γ γ ν = γ ν as the algebraic representation.
Somewhat more care must be taken with the graviton + quark vertex owing to the importance of scalarand tensor-meson poles.Considering Eqs. ( 10) -( 12), ( 18), (21), we are led to the following Ansätze where, on Extension to Q 2 ≳ 5 GeV 2 is discussed in Sec. 8. Fixing the masses in Eq. ( 26) to be those computed in RL truncation, listed in Eq. ( 22), then each graviton+quark vertex Ansatz has two parameters: κ q , r Sq .
Regarding κ q , we desire to implement a mass scale, κ q m Tq , in the form factor that is associated with the tensor channel.Having already used the ground state mass for the leading pole term, we chose κ u = κ s =: κ = the ratio of masses of the first radial excitation and ground state in the light-quark channel.Using the listings in Ref. [74,RPP], one finds a value of 1.13 for this ratio; so, we proceed with κ = 1.13, a fixed value, viz. it is not subsequently varied.Using this approach, one finds that the residues of the tensor poles in P T u,s are (in GeV): f Tu = 0.046, f Ts = 0.059.Given the simplicity of the Ansätze, these values are a fair match with the results in Eq. (22).
Again, given the simplicity of the Ansätze, the f Su,s values compare tolerably with the RL results in Eq. (22).
A common way to estimate uncertainties in RL predictions is to reevaluate all results with ±5% variations of ω in Eq. (13).Given the correlation between ω and derived mass scales, we translate this approach into an estimation of uncertainties via simultaneous ±5% variations of each mass scale in the algebraic Ansätze, viz.M q , m Sq , m Tq .The resulting change in a given value is the uncertainty listed in each instance below.

Electromagnetic and gravitational form factors
Using the RL results, described in Secs.4, 5, for all elements in the photon+meson current, Eq. ( 1), and its analogue for the graviton+meson current, one can directly proceed to deliver a unified set of predictions for π and K electromagnetic and gravitational form factors.The first studies of F π,K (Q 2 ) were completed in Ref. [52].Using brute-force numerical techniques, as therein, then in the calculation of each form factor one encounters moving singularities in the complex-l 2 domain sampled by the bound-state equations [75] such that there is a maximum value of Q 2 beyond which the evaluation of integrals like that in Eqs. ( 1) is no longer possible with conventional algorithms.
More advanced methods have been developed [19,42,59], exploiting the perturbation theory integral representation (PTIR) [80].However, constructing accurate PTIRs is time consuming.This is especially true in our case because one would need to build PTIRs for each quark propagator, Bethe-Salpeter amplitude, photon+quark and graviton+quark vertex considered herein, i.e., roughly 200 scalar functions.We therefore continue with a straightforward RL approach, computing all form factors on the directly accessible domain and then using the algebraic Ansätze, detailed in Sec.6, to assist in defining their ultraviolet completions.
As evident in Fig. 2A, using simple algorithms, the π RL calculation fails on Q 2 ≳ 2 GeV 2 .Thereupon, the algebraic Ansätze, Sec. 6, become valuable.On Q 2 ≲ 2 GeV 2 , they deliver results in agreement with the RL truncation.This justifies their use in developing ultraviolet completions of the RL predictions.
The results in Ref. [26] were obtained using algebraic Ansätze for pion and kaon generalised parton distributions (GPDs), constrained entirely by hadron-scale π and K valence-parton DFs.Owing to the D-term ambiguity [82] encountered, e.g., in connecting GPDs with gravitational form factors, reliable results for θ P 1 (0) could not be obtained in Ref. [26].Herein, there is no D-term ambiguity: as already seen for the pion, the RL truncation supplies a definite prediction for θ P 1 (t).The CSM predictions for kaon elastic electromagnetic and gravitational form factors are drawn in Figs.2B,  3.These figures show that, for the K, the RL calculation fails on Q 2 ≳ 1.7 GeV 2 .The algebraic Ansätze, Sec. 6, are valuable on this domain: since they deliver results in agreement with RL truncation on Q 2 ≲ 1.7 GeV 2 , they can be used to assist in defining ultraviolet completions of the RL predictions.
Fig. 3 Kaon gravitational form factors. Panel A. Mass, Legend.CSM predictions -solid curves; and results from algebraic Ansätze -dashed curves.Kaon total -purple; s quark in K -red; and u quark in Kblue.Pion comparison curves -green.
Figure 2B shows that the u quark in K elastic electric form factor is almost indistinguishable from the u quark in π form factor, F π .This weak environmental sensitivity is typical of RL results for light pseudoscalar meson elastic electromagnetic form factors [52].On the other hand, as seen elsewhere [19], there is a marked difference between the s in K and u in K electromagnetic form factors: , respectively.A similar distinction is expressed in the ratio of s in K and u in K valence parton DFs [24,84,85].
These comparisons confirm the character of the results for θ Ku,s 2 (0) in Fig. 3A; namely, the accentuated magnitude of flavour splitting is an artefact of RL truncation.The issue can be traced to the fact that an efficacious RL truncation must express all effects of EHM in the form of G (s) -see Eq. ( 13), achieving a description of hadron phenomena by overmagnifying the gaugesector interaction strength on s ≲ Λ 2 I [71].Moving beyond RL truncation, EHM is also manifested in corrections to the gluon+quark vertex [87], thereby enabling an explanation of observables with a far less infrared enhanced form of G ; hence, more realistic expressions of EHM+HB interference.
Further regarding the gravitational form factors, the ordering of kaon radii is the same as that found for the pion: r θ1 K > r F K > r θ2 K -see Table 2.Moreover, in each case, the net kaon radius is smaller than the kindred pion radius.In fact, averaging over all form factors: r K /r π = 0.85 (6).Here, matching the electromagnetic form factor pattern, the θ 1,2 K u quark in K radii are very like the analogous u in π radii: r θ1 Ku ≈ 1.09r θ1 π ; r θ2 π ≈ 1.03r θ2 Ku .

Pressure profiles
Working with the Poincaré-invariant elastic form factors described above, one may compute Breit frame charge, mass, pressure, and shear force distributions via appropriate three dimensional Fourier transforms [79].
Other transforms are possible [88]; but since the input function in any case is always the same, then whatever type of transform is chosen, it is just a mathematical operation on the same input object.Thus, judiciously interpreted, all outputs are qualitatively equivalent.Calculation of a Fourier transform requires knowledge of the subject form factor on Q 2 ∈ [0, ∞).However, our direct CSM predictions do not extend beyond Q 2 = 2 GeV 2 .The algebraic Ansätze discussed in Sec.6 can be used to assist in developing their ultraviolet completions.Those results are readily obtained on Q 2 ∈ [0, 5] GeV 2 ; but in order to arrive at realistic pictures, they must be augmented by inclusion of additional ln Q 2 suppression at ultraviolet momenta.This is necessary to express the fact that, in four spacetime dimensions, a quantum field theory treatment of hadron form factors introduces such scaling violations.Neglect-  4].
The required completions can be achieved by fitting the form factor results obtained using the algebraic Ansätze with the following function: where F 0 := F (y = 0).The coefficients are listed in Table 3.
Charge and mass distribution Breit frame density profiles, calculated from F = F P , θ P 2 , respectively, are fairly straightforward.In all cases, r 2 ρ F (r) vanishes at r = 0; and with increasing r, grows to a maximum value (at r ∼ 0.2 − 0.4), whose location is closer to r = 0 and magnitude is greater for form factors with smaller radii; then approaches zero from above at a rate that increases as the form factor radius becomes smaller.Naturally, consistent with the results in Table 2, the density profile of a meson's heavier valence quark is more compact than that of the lighter valence quark.
Pseudoscalar meson pressure and shear force distributions may be defined as follows [79]: where isospin symmetry is assumed, q = u, s, ∆ = √ Q 2 , E(∆) 2 = m 2 P + ∆ 2 /4 and j 2 (z) is a spherical Bessel function, For bound states, ∞ 0 dr r 2 p q P (r) ≡ 0. As revealed by Fig. 4, these densities possess more interesting structure than the charge and mass profiles.Semiquantitatively, our ultraviolet completed RL predictions yield the same pictures obtained using the GPD models detailed in Refs.[25,26].
Consider first the pressure profiles, drawn in Fig. 4A.Compared with the GPD model results [25,26], these profiles are more compact because the RL predictions for θ As noted elsewhere [25,26,89], profiles like those in Fig. 4A can be drawn for neutron stars.They indicate r ≃ 0 pressures therein of roughly 0.1 GeV/fm [90].Plainly, the form factors computed herein yield core pressures in Nambu-Goldstone bosons that are similar in magnitude to that of a neutron star.
The K profile is more compact, so the associated peak core pressure is higher than in the π: the ratio is ≈ 1.5.Further, the s-quark contributes more of the total K pressure than its partner u-quark and its peak/trough intensities are greater.
The shear pressures, drawn in Fig. 4B, are an indicator of the strength of deformation forces within the meson.Evidently, these forces are maximal in the neighbourhood upon which the pressure changes sign.Within these neighbourhoods, the forces driving the quark and antiquark apart are overwhelmed by attractive confinement pressure.It is interesting to compare the total π and K shear forces: ( The π result is greater owing to the long tail of s π (r): on r ∈ [0, 0.9], the integrated strengths are equal.The ordering in Eq. ( 29) is opposite to that found using the GPD models in Refs.[25,26].This difference owes largely to (a) the failure of those models to express the suppression of θ K 1 (0) -see Table 2; and (b) a low value for r θ1 Ku .

Summary and perspective
Using continuum Schwinger function methods, we delivered predictions for π and K elastic electromagnetic, F π,K , and gravitational, θ π,K 1,2 , form factors, thereby unifying them with numerous additional properties of these Nambu-Goldstone bosons and other hadrons.The calculations were completed using a symmetry-preserving formulation of each quantum field equation relevant to the calculation of the bound-state wave functions and interaction currents.In particular, we described and calculated the structure of dressed graviton + quark vertices [Secs.3,5].The analysis highlights that, just as the low momentum transfer behaviour of hadron electromagnetic form factors is sensitive to the properties of the lightest neutral vector meson in a given quark+antiquark channel, so their mass and pressure form factors are sensitive, respectively, to the analogous tensor and scalar mesons.
For both π and K, the pressure distribution, θ 1 , is softer than F , which is than the mass distribution, θ 2 ; and each K form factor is harder than its kindred π form factor [Sec. 7].Examining form factor flavour separations, the sizes of the effects are typical of those found when considering interference between Nature's two known sources of mass, i.e., emergent hadron mass and Higgs-boson couplings into QCD.
Regarding density and pressure profiles -Sec.8, defined as form factor Fourier transforms, each follows the expected pattern, e.g., K profiles are more compact than those of the π, and the K core pressure is higher.Our direct calculations, which do not suffer from the D-term ambiguity encountered in GPD-based analyses, confirm that the pressure at the core of a Nambu-Goldstone boson is of the same magnitude as that in the heart of a neutron star.Considering the pressure profiles, one might define a confinement radius as the location at which these functions cross zero.We find this radius to be ≈ 0.39 fm for the pion.The kaon value is roughly 33% smaller, with the s in K radius being ≈ 15% smaller than that of the partner u quark.Contrary to earlier phenomenological results, the total shear pressure in the π exceeds that in the K.
The analysis described herein should enable a direct calculation of nucleon gravitational form factors using either a quark+diquark picture [91] or a three-body approach [44,92].

Fig. 1
Fig. 1 Probe+pion interaction in RL truncation.Both panels.Solid lines -dressed quarks, S; orange shaded circlespion Bethe-Salpeter amplitude, Γ π ; and blue shaded circles -dressed probe+quark vertex, Γ P .Lower panel.Shaded C region -gluon-binding contribution to the probe+pion interaction, which vanishes for F π but cancels a Q-longitudinal contribution from the top diagram when computing the gravi-ton+pion coupling, thereby ensuring cπ ≡ 0 in Eq. (3).Analogous diagrams define the kaon form factors.

Fig. 4
Fig. 4 Pressure (Panel A) and shear (Panel B ) distributions obtained using Eqs.(28) as explained in Sec. 8. Legend.πsolid green curves; K -purple curves; s in K -dashed red curves; and u in K -dashed blue curves.

Pq 1
fall faster with increasing ∆.The qualitative features of the profiles in Fig.4Asuggest a physical interpretation.Namely, the meson pressures are positive and large on the neighbourhood r ≃ 0, whereupon the meson's dressed-valence constituents are pushing away from each other.With increasing separation, the pressure switches sign, indicating a transition to the domain wherewithin confinement forces exert their influence on the pair.The zeros occur at the following locations (in fm): r π c = 0.39(1), r K c = 0.26(1), r Ku c

Table 3
Inserted into Eq.(27),these parameters define realistic ultraviolet completions of the π, K electromagnetic and gravitational form factors.