Gluon distributions and mass decompositions of the pion and kaon

We present the gluon distribution functions of pion and kaon in small-$x$ and large-$x$ regions, and compare them with the results obtained from lattice QCD and continuum Schwinger function methods. Whether in the small-$x$ region or the large-$x$ region, our gluon distribution of pion is consistent with the results of lattice QCD and continuum Schwinger function method. In addition, the first four moments of gluon distributions of pion and kaon at different $Q^{2}$ scales are calculated. Furthermore, we present the mass decompositions of pion and kaon with the dynamical parton distribution functions calculated by the DGLAP equation with parton-parton recombination corrections. The mass structures of pion and kaon are totally different from that of the proton.


I. INTRODUCTION
The pion, the lightest bound state in quantum chromodynamics (QCD), plays an important role in nuclear physics as it is a Nambu-Goldstone boson [1,2] with dynamical chiral symmetry breaking (DCSB).For the kaon meson, although it was discovered in the 1950s [3], little is known about its structure.The studies of the structures of pion and kaon can reflect the physics of DCSB, which is helpful to reveal the relative influence of DCSB on the chiral symmetry breaking of the quark masses, and it is of great significance to understand the nonperturbative QCD.Studying the pion and kaon parton distribution functions (PDFs) is not only for characterizing their structures but also for further understanding of DCSB and non-perturbative QCD.At present, people know less about the pion and kaon PDFs than the nucleon PDFs, for the experimental data are scarce, especially for the gluon distributions.
Most of the global analyses of pion PDFs mostly rely on Drell-Yan data.In our previous global QCD analysis [4], studies of pion and kaon PDFs were based mostly on pion-induced Drell-Yan data [5][6][7][8], and the leadingneutron deep inelastic scattering (DIS) data of e − p collision at HERA [9,10].There are some recent studies, including the work of Bourrely and Soffer [11] that extract the pion parton distributions from the Drell-Yan π+W data with quantum statistical approach.JAM Collaboration [12,13] performed the first Monte Carlo global QCD analysis of the Drell-Yan πA data and the leading-neutron DIS data from HERA to reach the low-x region.This Monte Carlo global QCD analysis of pion PDFs by the JAM Collaboration reveals that the gluons carry significantly higher momentum fraction than that only inferred from Drell-Yan data.The xFitter developers' team [14] used their open source QCD fitting framework to promote a PDF extraction from the Drell-Yan πA and photon production data.It is found that these data can well constrain the valence distribution, but they are not sensitive enough to the accurate determination of sea and gluon distributions.Chang et al. [15] show that the pion-induced J/ψ production data can impose useful additional constraints on pion PDFs, especially for the gluon distribution of the pion in the large-x region.
As we know, the studies of pion and kaon gluon distributions from a theoretical perspective can provide useful information for experiments.Most theoretical model calculations only predict the valence quark distribution of pion [16][17][18][19][20][21][22][23][24][25][26][27][28], while gluon and sea PDFs are predicted by the Dyson-Schwinger equation (DSE) continuum method with QCD evolution [18].The description of the pion gluon PDF that based on the Rainbow-Ladder truncation of DSE is consistent with the pion gluon PDF results of JAM Collaboration [12,13] within two sigma.Recently, Fan and Lin present the x-dependent gluon distribution of pion from lattice-regularised QCD (LQCD) with the pseudo-PDF approach [29] for the first time.
In this work, we present our x-dependent pion gluon distribution from a global QCD analysis of the available experimental data [4] and compare it with the results of LQCD calculations, continuum Schwinger function methods (CSMs), and other theoretical calculations.The xdependent gluon distribution of kaon is also presented under a certain scale Q 2 .In addition, the first four moments of gluon distributions of pion and kaon at different Q 2 scales are calculated and compared with some model results of DSE [30], JAM [12,13], xFitter [14] and LQCD [31] respectively.Furthermore, according to the hadron mass decomposition which is based on the structure of the QCD energy-momentum tensor and the pion and kaon PDFs from our global QCD analysis, we study the mass structures of the pion and kaon.
dicted the pion's gluon distribution with CSMs [30].In this work, we present pion and kaon PDFs with uncertainties from a global QCD analysis of the experimental data within the framework of dynamical parton model [4].Figure 1 displays the normalized pion gluon distribution xg π (x, Q 2 )/ ⟨xg π ⟩ at Q 2 = 4 GeV 2 from us, compared with the LQCD [29] and DSE [30,32] calculations.The top pad in Fig. 1 shows the gluon distribution in the small-x region, and the bottom pad in Fig. 1 shows the gluon distribution in the large-x region.There shows agreements among the three results of gluon distribution on the entire depicted domain.Judged by the agreement between our work and other theoretical predictions (LQCD and CSMs) for the pion gluon distribution in Fig. 1, our pion gluon distribution obtained from global QCD analysis is also reliable.[30,32] and LQCD [29] calculation results for the pion gluon distribution in the small-x region.(Bottom) Comparisons of our result to the CSMs [30,32] and LQCD [30,32] for the pion gluon distribution in the large-x region.
The pion gluon distribution xg π (x, Q 2 ) we calculated at Q 2 = 0.61 GeV 2 is shown in Fig. 2, compared with the gluon distribution from the dressed quark [33] in the pion.By comparison, one sees that our result based on the dynamical parton model is consistent with Freese's result.The dressing-gluon is derived from the constituent quarks, thus when x<0.1, the distribution of dressinggluon decreases.For our result, most of the gluons generated under the dynamical parton model are from parton splitting process, thus when x<0.1, our dynamical parton model still has the gluons which are not related to the dressed quarks.Figure 3 shows our pion gluon distribution xg π (x, Q 2 ) at Q 2 = 10 GeV 2 , compared to the pion gluon distribution from the framework of a statistical model fitted to the pion-induced Drell-Yan data and J/ψ production data of the pion [34].
Figure 4 shows our normalized gluon distribution xg K (x, Q 2 )/ xg K of the kaon in the large-x region at Q 2 = 4 GeV 2 .By comparing the kaon gluon distri-butions, one sees that our result is consistent with the DSE'20 [30] in the M S scheme and the lattice QCD calculation [31].Our kaon gluon distribution xg K (x, Q 2 ) at Q 2 = 10 GeV 2 from a global QCD analysis is shown in Fig. 5, compared with Bourrely's result from an analysis of meson-induced Drell-Yan and quarkonium production data within a statistical model [35].
In order to compare with more other theoretical calculations, we calculate the first four moments of the pion and kaon gluon distributions at different Q 2 .The moments of the gluon distribution is defined as follows, where g(x, Q 2 ) is the gluon distribution.
Figure 6 shows the first four moments of pion and kaon gluon distributions as a function of Q 2 scale.The blue solid and red dashed curves show the evolutions of the first four moments of pion and kaon gluon distributions, respectively.The (a), (b), (c), and (d) pads present the first-order moment, second-order moment, third-order moment, and fourth-order, respectively.The  g / ⟨x⟩ g , x 3 g / ⟨x⟩ g for pion and kaon, which are in good agreements with the corresponding results from DSE [30], JAM [12,13], xFitter [14] and LQCD [31] approaches.

III. MASS STRUCTURES OF THE PION AND KAON
The mass structure of pion is very interesting.According to Goldstone theorem, the pion is essentially different from other ordinary hadrons, which is a collective mode in QCD vacuum.We will see that the mass structure does reflect that.The QCD structure of hadron mass can be decomposed into four different parts: in which M q is the quark energy contribution, M g is the gluon energy contribution, M m is the quark mass contribution and M a is the trace anomaly contribution.
These four mass terms mainly depend on the momentum fraction carried by all the quarks from our analysis [4], QCD trace anomaly parameter b calculated via chiral perturbation theory [38] and the quark mass anomalous dimension γ m .The quark mass anomalous dimension γ m [39] is predicted by perturbative QCD as follows: where a s = α s /π = g 2 /(4π 2 ), and g denotes the renormalized strong coupling constant.It is worth noticing that the anomalous dimension γ m defined by Ji [38] is -2 times [40] that defined in Ref. [39].For the running coupling constant α s , we adopted Cui's latest parametrization of the effective charge from lattice QCD [41].The adopted running coupling constant α s including the gluon mass effect of dynamical chiral symmetry breaking is written as: In general, the matrix element a(Q 2 ) in Eq. ( 3) can be deduced from the quark distributions that are measured in π-N Drell-Yan process.The matrix element a(Q 2 ) is the momentum faction carried by the quarks in the pion, which is given by, Here q f (x, Q 2 ) and qf (x, Q 2 ) are the quark and antiquark distributions with flavor f in the pion.In this work, we use the pion PDFs from our recent global QCD analysis of the experimental data to evaluate the matrix element a(Q 2 ).
For the trace anomaly parameter b in Eq. ( 3), it can be derived from the pion mass in the chiral perturbation theory.The first-order chiral perturbation theory gives a clean prediction of the trace anomaly parameter to be b π = 0.5 [38].Here we provide a brief explanation of the criteria for estimating the trace anomaly parameter b.First, according to the original derivation of mass decomposition and the requirement of Lorentz covariance [38], b is viewed to be a scale-invariant quantity which has been discussed in many previous studies [38,40,[42][43][44].Two estimates of trace anomaly parameter b of proton in the limits of chiral SU(3) and heavy strange quark are b p = 0.17 and b p = 0.11 provided by Ji [38].Wang et al. firstly present an extraction of the QCD trace anomaly parameter to be b p = 0.07 ± 0.17 [40] from the experimental data of nearthreshold J/ψ photoproduction cross section from GlueX at JLab.Recently, Kou et al. analyze the updated differential cross-section data from Hall C and GlueX Collaborations [43,45], and obtain a new value of trace anomaly parameter b p = 0.52 ± 0.09 [44], which is close to the trace anomaly parameter for the pion estimated in the chiral perturbation theory [46].But this does not necessarily prove that b is independent of the hadrons.Second, the mass decomposition of different hadrons should be different.Therefore, the trace anomaly parameter b could be different for different hadrons.Third, we assume that in the chiral limit, the first-order chiral perturbation theory under SU(3) gives the same mass formula for the pion and the kaon.This is due to the fact that the vacuum condensations of quarks with three flavors are treated equally [46].Therefore the parameter b of the kaon is also taken to be b K = 0.5 in this work.
Using the trace anomaly parameter b explained above, we calculated the mass decompositions of pion and kaon at different Q 2 scales.Figure 7 shows the pion mass decomposition as a function of Q 2 , where the different markers indicate the percentages of different contributions to the pion mass.By comparing to the mass decomposition of the proton [40], it is found that the proportion of the quark mass contribution in the pion is smaller than that in the proton, while the proportion of the trace anomaly contribution in the pion is larger than that in the proton.The reason for the large difference between the mass decompositions of the pion and the proton is that the trace anomaly parameter b p ≈ 0.1 of the proton [38,40] is much smaller than that of the pion.Quark mass plays a more important role in the pion mass comparing to the proton, as the chiral-limit mass of the pion is zero while the chiral-limit mass of the proton is huge.Moreover, the proportion of quark energy M q in the pion is smaller than that in the proton, especially at high Q 2 scale.The quark energy M q is the sum of kinetic energy and potential energy.The reason for the smaller proportion of M q in the pion is that the quark potential energy in the pion is lower compared to that in the proton.In addition to the small proportion of M a to the pion mass, the small proportion of M q may be another reason why the pion mass is so light.Figure 8 shows the kaon mass decompositions at different Q 2 scales, where the different markers indicate the percentages of different contributions to the kaon mass.The Q 2 -dependences of different mass origins in the pion or kaon mass decomposition are discussed.Firstly, we emphasize that the trace anomaly is scale-invariant, which is governed by the symmetry of QCD energymomentum tensor.In other words, the trace anomaly parameter b is scale-invariant and it is associated with the hadron's scalar charge [38].Secondly, the quark mass term M m depends on the anomalous dimension γ m and this gives the Q 2 -dependence of M m .Thirdly, the Q 2dependence of the quark energy term M q depends on the momentum fraction carried by all quarks a(Q 2 ), and the traditional DGLAP evolution tells us that a(Q 2 ) decreases as Q 2 increases.This is the main reason why the quark energy term M q is negatively correlated with the Q 2 scale.Fourthly, the gluon energy term M g increases with increasing Q 2 , and it is also governed by the DGLAP evolution equation.With the increasing Q 2 , the parton splitting processes radiate more and more gluons, resulting in a decrease of quark kinetic energy, and an increase of gluon kinetic energy, which is in Fig. 6.In whole description of the hadron mass decomposition, the anomalous dimension γ m only affects the terms related to quarks, as the gluon is massless.The gluon contributions to the hadron mass only depends on the trace anomaly and the momentum fraction carried by the gluons.

IV. SUMMARY
The pion and kaon gluon distributions from our previous global QCD analysis applying DGLAP equations with parton recombination corrections are basically consistent with the DSE and LQCD calculations in both small-x and large-x regions.The future experimental measurements will provide excellent chances to verify these predictions.Additionally, the first four moments of the gluon distributions of pion and kaon at different Q 2 scales are predicted, which agree well with some theoretical models.
We also give the mass structures of pion and kaon at different Q 2 scales, according to the mass decomposition of QCD energy-momentum tensor.The mass of the pion or the kaon can be decomposed into four parts in this work.For the pion and the kaon, the quark mass contribution M m is the largest one contributed to the meson mass at high Q 2 .Compared to the proton, the trace anomaly contribution of pion or kaon is small.These four origins of the hadron mass depend on the momentum fraction a by the quarks, the trace anomaly parameter b, and the quark mass anomalous dimension γ m .In this work, we find some apparent differences between the mass decompositions of the pion and the proton [40].Besides the trace anomaly contribution M a , we notice that the quark energy contribution M q , the gluon energy contribution M g , and the quark mass contribution M m all have obvious dependence on Q 2 .
The US Electron Ion Collider (EIC) [47][48][49] under the construction and the proposed Electron-Ion Collider in China (EicC) [50][51][52] provide great opportunities to deepen our understanding on the meson structures through the measurements of Sullivan process [48,53,54].These high-energy projects will surely increase our knowledge about the pion and kaon gluon distributions and the emergence of the meson mass.

FIG. 1 .
FIG. 1. (Color online) (Up) Comparisons of our result to theCSMs[30,32] and LQCD[29] calculation results for the pion gluon distribution in the small-x region.(Bottom) Comparisons of our result to the CSMs[30,32] and LQCD[30,32] for the pion gluon distribution in the large-x region.

FIG. 3 .
FIG. 3. (Color online) Our obtained pion gluon distribution from a global QCD analysis compared with the result from a statistical model constrained by pion-induced Drell-Yan data and J/ψ production data at Q 2 = 10 GeV 2 [34].

2 Q
FIG. 4. (Color online) Comparisons of our kaon gluon distribution with the DSE[30] and LQCD[31] results in the large-x region.

FIG. 6 .
FIG. 6. (Color online)The first four moments of the pion and kaon gluon distribution functions as a function of the square of four momentum transfer Q 2 .The (a), (b), (c), and (d) panels show the first-order moment, second-order moment, third-order moment, and fourth-order moment, respectively.The blue solid and red dashed curves show the first four moments of pion and kaon gluon distributions, respectively.

FIG. 7 .
FIG. 7. (Color online) The QCD decompositions of the pion mass under different Q 2 scales, where the different markers indicate the percentages of different contributions to the mass as a function of Q 2 .

FIG. 8 .
FIG. 8. (Color online) The decompositions of the kaon mass under different Q 2 scales, where the different markers indicate the percentages of different contributions to the kaon mass.

TABLE I .
Comparisons of our result from a global QCD analysis in the dynamical parton model with some other model calculations for ⟨x⟩ π g and ⟨x⟩ K g at Q 2 = 4, 27 GeV 2 .This work 0.38 4 first moments of the gluon distributions of pion and kaon are listed in TableI.Our ⟨x⟩ predictions agree well with the other model calculations.It is worth noting that LQCD and DSE are two main nonperturbative approaches derived directly from the QCD theory.TableII π g and ⟨x⟩ K g

TABLE II .
Comparisons of our result from a global QCD analysis in the dynamical parton model with some other model calculations for x 2 g / ⟨x⟩ g , x 3 g / ⟨x⟩ g of pion and kaon at Q 2 = 4 GeV 2 .