{'fraction_non_alphanumeric': 0.04235826686785945, 'fraction_numerical': 0.02137673349141338, 'mean_word_length': 4.3170964660936, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Tomography of light mesons in light-cone quark model

We investigate the tomographical structure of pion and kaon in light cone quark model (LCQM). In particular, we study the parton distribution amplitude (PDA) of pion and kaon. We obtain the parton distribution function (PDF) and the generalized parton distributions (GPDs) of the pion and kaon. The valence quark PDA and PDF of pion, after QCD evolution, are found to be consistent with the data from the E791 and the E615 experiments at Fermilab, respectively. Further, we investigate the transverse momentum distributions (TMDs) of pion and kaon. We also discuss the unpolarized TMD evolution for pion and kaon in this model.


INTRODUCTION
The nonperturbative structure of hadron is well described by the distribution of partons inside the hadron in both position and momentum space. The distribution amplitudes (DAs) are among the most basic quantities which not only provide important information on bound states in QCD but also play an essential role in describing the various hard exclusive processes [1,2] of QCD via the factorization theorem [3] analogous to parton distributions in inclusive processes. DAs are the longitudinal projection of the hadronic wave functions obtained by integrating out the transverse momenta of the constituents of the hadron [4,5]. The lowest moments of the hadronic DAs for a quark and an antiquark inside a meson also give us the knowledge of decay constants and transition form factors [6][7][8]. The parton distribution functions (PDFs) [9,10], which are accessible in hard inclusive processes such as deep inelastic scattering (DIS) or Drell-Yan processes, encode the distribution of longitudinal momentum and polarization carried by the constituents. The generalized parton distributions (GPDs) [11][12][13][14] reveal the parton distribution in the direction transverse to the hadron motion providing the spatial distribution. Unlike the PDFs which are function of longitudinal momentum fraction carried by the active parton x only, GPDs being function of x, the longitudinal momentum transferred ζ, and the total momentum transferred from initial state to final state of the hadron t, provide us the knowledge of 3-dimensional (3D) spatial structure of hadron. One can obtain the form factors, charge distributions, PDFs etc. from GPDs under certain conditions [15][16][17]. The momentum tomography of hadrons is described by the transverse momentum-dependent parton distributions (TMDs) [18][19][20]. The TMDs are function of longitudinal momentum fraction x and transverse momentum possessed by the parton (k ⊥ ).
From the theoretical point of view, the pion DAs have been studied using Chernyak-Zhitnitsky (CZ) [74] and 75] approaches. Further, the pion and the kaon DAs have been theoretically calculated using the Dyson-Schwinger equations [76,77] as well as from Poincaré-covariant Bethe-Salpeter wave functions [78]. The pion PDFs have been studied in Nambu-Jona-Lasinio (NJL) model [79]. The valence-quark pion and kaon PDFs incorporated with gluon contributions have been studied in Ref. [80]. The pion and the kaon DAs and PDFs have been discussed in light-front constituent quark model by using the symmetric quark-bound state vertex function [81] and chiral constituent quark model [82]. The pion PDF has also been the subject of detailed analyses in the phenomenological models in Refs. [83,84], also including anti-de Sitter (AdS)/QCD models [85][86][87][88] and the chiral quark model [89]. The pion PDFs have also been studied within lattice QCD [90][91][92][93][94][95]. Additionally, the first global fit analysis of PDFs in the pion has been reported in Ref. [96]. Although the PDFs are expected to be universal, tension exists regarding the behavior of the pion valence PDF. The large-x behavior of the pion valence PDF is expected to fall off linearly or slightly faster from the analyses of the Drell-Yan data [31,33]. This is supported by the constituent quark models [83,84], the Nambu-Jona-Lasinio (NJL) model [84], and duality arguments [97]. However, this observation disagrees with perturbative QCD, where the behavior of the same function has been predicted to be (1 − x) 2 [98][99][100][101], a behavior further supported by the Bethe-Salpeter equation (BSE) approach [102,103]. Meanwhile, the reanalysis of the data for the Drell-Yan process [35] including the next-to-leading logarithmic threshold resummation effects shows a considerably softer valence PDF at high x when compared to the NLO analysis [31,33]. The pion GPDs have been attempted using covariant and light-front constituent quark models [104], two chiral quark models-NJL model and spectral quark model [89]. Recently, the GPDs of pion for zero skewedness [85] and for non-zero skewedness [105] have been studied with AdS/QCD approach. Furthermore, the studies on pion and kaon GPDs are explained in Poincaré covariant Bethe-Salpeter constituent quark model by considering the support parameter of pseudoscalar mesons [106].
TMDs contain important information on the 3D internal structure of hadrons, especially the spin-orbit correlations of quarks within them [107]. The pion TMDs have been studied by considering the Drell-Yan process with pion beams [108]. The TMD fragmentation functions of elementary particles in pion and kaon from the NJL-jet model have been explained in Ref. [109] using Monte-Carlo approach. The pion TMDs are also evaluated in NJL model with Pauli-Villars regularization [110]. The transverse structure of the pion in momentum space inspired by light-front holography has been reported in Refs. [111,112].
To understand the relativistic effects of the motion of quarks and gluons in the hadrons, light-cone formalism is used which is a convenient frame to study the applications to the exclusive processes. The Wigner-rotation is taken into account, when one transforms a composite system from one reference frame to another. The advantage of using the light-front dynamics is that the Wigner rotation related to the spin states is unity in different frames under the Lorentz transformation. The light-cone quark model (LCQM) finds application in QCD low-scale regime. The pion has chiral symmetry constraints, particularly the explicit chiral and spontaneous symmetry breaking, leading to the pion structure being the simplest valence-quark substructure to study. The LCQM is successful in explaining the electromagnetic form factors of pion and kaon. The results of electromagnetic form factors (EFFs) have already been compared with the experimental data available at low-energy scale in Ref. [113,114], and is found to be consistent with the results evaluated in this model. The decay constants and charge radii for both particles have also been predicted in the LCQM. In light of the progress, it therefore becomes necessary to enhance this model to study several distributions of partons in the mesons.
In the present work, we have implemented the Melosh-Wigner transformation to derive the light-cone spinflavor wave functions of pion and kaon. We have investigated the DAs of quark in pion and kaon at the model initial energy scale µ 0 . The QCD evolution is applied on the low-energy scale model calculated DAs and compared with asymptotic result. The moments at different evolution scales corresponding to pion and kaon are also compared with different available theoretical predictions. Further, we evaluate the PDFs for both pseudoscalar mesons in this model. To compare the model results with the available experimental data of PDFs, the nextto-next-leading order (NNLO) DGLAP evolution is put into consideration by choosing the appropriate scales. We have studied the GPDs and the TMDs of pion and kaon from the overlap of light-cone wave functions (LCWFs). As the pseudoscalar mesons in their ground states possess a simple valence quark substructure, one is left with only one GPD and one TMD. In addition to this, the scale evolution of unpolarized pion and kaon TMDs is also discussed.
The manuscript is organized as follows. In Section II, we discuss the details of light-cone quark model. The wave functions and distribution amplitudes for pion and kaon are detailed in Section III. Section IV includes the parton distribution functions for both pion and kaon mesons and its DGLAP evolution to the higher energy scales. In Section V, a detailed description and graphical interpretations of the GPD are given for u-quark in pion and kaon. In Section VI, we present the unpolarized quark TMD of both pseudoscalar mesons. In this section, the details of TMD scale evolution is also presented. Finally, the results are concluded in Section VII.

II LIGHT-CONE QUARK MODEL (LCQM)
The hadron eigenstate |M (P + , P ⊥ , S z ) in connection with multi-particle Fock eigenstates |n is defined as [115] |M (P + , P ⊥ , S z ) = n,λi where P = (P + , P − , P ⊥ ) is considered as the total momentum of meson and S z is the longitudinal spin projection. The momenta of meson having mass M and its constituents having masses m 1 and m 2 in light-cone frame are defined as The multi-particle Fock states containing n constituents where the ith constituent holding the light-cone longitudinal momentum fraction x i = k + i P + , the transverse momentum k ⊥i and helicity λ i are normalized as The light-cone wave function in LCQM is written as where ϕ and χ correspond to the momentum space and spin wave functions respectively and superscript F stands for the front form. The LCWF of pion (or kaon) can be obtained through the transformation of the instant-form SU(6) wave functions using Melosh-Wigner rotation. The spin wave function of the pseudoscalar meson in the instant form (T ) can be written as [114,116] where χ ↑,↓ i is the two-component Pauli spinor. One can relate the light-cone spin states |J, λ F and the ordinary instant-form spin states |J, s T as follows where U J is the Melosh-Wigner rotation operator. The spin space wave function of pseudoscalar meson can obtained in the infinite momentum frame by implementing the transformation Eq. (8) in Eq. (7). The Melosh-Wigner transformation is used to connect the instant form spin states and light-front form spin states as Here we take the instant-form four-momenta for 2 quarks as: (9) and (10), A meson is a bound state of a quark (Q) and an antiquark (Q) viz. QQ, where the masses of two partons are denoted as m 1 and m 2 . For pion, we consider m 1 = m 2 = m whereas for kaon, we take m 1 = m 2 because of its composition.
The light-cone spin wave function of pseudoscalar 'P' meson, which can be pion or kaon depending on their composition, has the form with S z and λ being the spin projection of pion (or kaon) and quark helicity, respectively. Since for pion (kaon) having masses m (m 1 and m 2 ), the z-component of spin is zero (S z = 0), therefore, the component coefficients κ F Sz=0 (x, k ⊥ , λ 1 , λ 2 ) in spin wave function are indicated as Here x i (i = 1, 2) is the light-cone quark momentum fraction in light-front dynamics with the constraint This leads to x 1 + x 2 = 1 for the case of meson. If we assume the momentum fraction of one parton x 1 = x, then for the other parton it becomes The component coefficients κ F Sz=0 (x, k ⊥ , λ 1 , λ 2 ) in spin wave function given in Eq. (12) must satisfy the following normalization conditions for pion (or kaon) The momentum space wave functions ϕ π(K) (x, k ⊥ ) in Eq. (6) are adopted using Brodsky-Huang-Lepage (BHL) prescription. For pion we have, and for kaon we have where A π and A K are the normalization constants for pion and kaon, respectively.
The two-particle Fock state expansion can be described in terms of LCWFs ψ Sz (x, k ⊥ , λ 1 , λ 2 ) defined in Eq. (6) and we have

III WAVE FUNCTIONS AND DISTRIBUTION AMPLITUDES (DAS)
The constituent quark masses and the harmonic scale β are only two input parameters required to compute the pion and the kaon distribution funtions. The parameters used in the present work are listed in Table I. Comparison between the momentum space wavefuntion of the pion ϕ π (x, k ⊥ ) and the kaon ϕ K (x, k ⊥ ) is shown in Fig.  1. The pion wave function shows symmetry over x = 0.5 whereas, due to dissimilar quark masses, kaon wave function appears to be asymmetrical along x. LFWFs give unique access to light cone distributions by integrating out the transverse momentum [4]. Among those, the DAs control the exclusive processes at large momentum transfer. In light-front formalism, the leading-twist DAs for pseudoscalar mesons are defined through the correlation [117][118][119][120]  where P represents pseudoscalar meson and the decay constant is denoted with f P . Substituting the pseudoscalar meson states and the quark field operators lead to the definition of DAs in terms of LCWFs, one has [117,121] with the normalization condition at any scale: Using the LFWFs given in Eq. (6), we compute the DAs of the pion and the kaon at the model scale. Next, the LO QCD evolution of the DAs is carried using the Efremov-Radyushkin-Brodsky-Lepage (ERBL) equations [4,5]. In a Gegenbauer basis, one has [122] φ with a n (µ) = 2(2n + 3) 3(n + 1)(n + 2) where C 3 2 n (2x − 1) is a Gegenbauer polynomial. The strong coupling constant α(µ) is given by The factor γ (0) n 2β0 defines the anomalous dimensions and where c A = 3 and n F correspond to the number of active flavors. The color factor c F = 4 3 and Λ QCD = 0.226 GeV [86].
In Fig. 2(a) and fig. 2(b), we show the evolution of the pion and the kaon DAs, respectively from the initial scales µ 2 0 to µ 2 = 10 GeV 2 which is the scale relevant to the E791 data [123]. As can be seen, the pion DA is close to the asymptotic DA already at µ 2 = 1 GeV 2 and the DA approaches towards the asymptotic DA with increasing evolution scale, however, the effect is small. The evolved pion DA in this model shows a good agreement with the E791 data. Unlike the pion, the evolved kaon DA is still distinct from the asymptotic DA at µ 2 = 1 GeV 2 or even at µ 2 = 10 GeV 2 . However, the trend of evolution is same as the pion DA.
It is useful to compute the moments in order to quantitatively compare with other theoretical predictions. The n-th moment is defined as, where z can be ξ = (2x − 1) or x −1 . Our predictions for the moments of the pion DA are compared to other theoretical approaches in Table II. From the table, it is important to note that the LCQM DAs have their moments larger than their asymptotic values in case of pion, and the value goes on decreasing when evolve towards the latter from below. LCQM moments show a similar behavior as the moments obtained in other theoretical approaches except LF holographic model, and our moments for the pion DA are closer to the asymptotic results compared to other approaches. The moments of the kaon DA are compared to other theoretical predictions in Table III. In this case, the even moments in LCQM model are lower than their asymptotic values, while the predicted odd moments are greater than zero, which is their asymptotic value. Further, the inverse moment shows the higher value compared to the asymptotic one. We notice a similar trend as observed in LF holograpic model [87] where the exception lies in inverse moment. It can be noticed that our moments for the pion and kaon DAs are close to the asymptotic results when the DAs are evolved to a large scale ∼ 10 4 GeV.

IV PARTON DISTRIBUTION FUNCTIONS (PDFS)
The pion (kaon) PDF gives the probablity of finding the quark in the pion (kaon) where the quark carries a longitudinal momentum fraction x = k + /P + . At fixed light-front time, the PDF can be expressed as [125] f P (x) = 1 2 Since the spin is zero in both the cases S = 0, we deal with unpolarized parton distribution function which comes from the above relation by substituting Γ = γ + . The overlap form of PDF is defined by putting the pion (kaon) states from Eq. (18). We have Using the LCWFs given in Eq. (6), we evaluate the quark distribution functions, f π(K) (x) at the initial scale and plot them as a function of x in Fig. 3. Due to equal constituent quark (antiquark) mass, the distribution in the pion appears to be symmetric over x = 0.5, while in kaon, the light quark distribution is maximum at a slightly lower value of quark momentum fraction in the longitudinal direction. The peak is broader in the case of the pion as compared to the kaon. The distribution peak has higher amplitude in case of the kaon as compared to  We now have our PDFs for the light mesons at scales relevant to constituent quark masses which are several hundred MeV. At the model scales, both PDFs for the valence quark (antiquark) is normalized to 1: Meanwhile, the following momentum sum rule: is a consequence of Eq. (30). This states that the valence quark and antiquark together carry the entire light front momentum of the meson, which is appropriate to a lowresolution model.  [73] and blue colored points represents the data for modified FNAL-E615 data [35]. The red solid curve denotes the valence quark PDF obtained from BLFQ [145].

(A) QCD Evolution for Pion PDF
The valence quark distributions at high µ 2 scale can be determined with the initial input by performing the QCD evolution. We adopt the NNLO DGLAP equations [138][139][140] of QCD, to evolve our PDFs from our model scales to higher scales µ 2 needed for the comparison with experiment. The scale evolution allows quarks to emit and absorb gluons while the emitted gluons allowed to generate quark-antiquark pairs as well as additional gluons. In this picture, the higher scale reveals the gluon and sea quark components of the constituent quarks through QCD.
We explicitly evolve our initial PDFs from the LCQM model for the pion to the relevant experimental scales µ 2 = 16 GeV 2 using the Higher Order Perturbative Par-ton Evolution toolkit to numerically solve the NNLO DGLAP equations [141]. While applying the DGLAP equations numerically, we impose the condition that the running coupling α s (µ 2 ) saturates in the infrared at a cutoff value of max α s = 1 [142][143][144][145].
In Fig. 4, we compare our result for the valence quark PDF of the pion with the FNAL-E-0615 [73] and FNAL-E-0615 modified data [35]. Also, the theoretical result of the valence quark PDF evaluated by using the basis light front quantization (BLFQ) [145] is taken for comparison with our result. We notice that our result does not fit to FNAL-E-0615 experiment. Meanwhile, our result fits to the modified E615 data after being reanalyzed to take into account soft gluon resummation [35]. We are able to fit the reanalyzed data using an initial scale µ 2 0 = 0.246 GeV 2 (similar to the initial scale used in Refs [86,87,122]). At the initial scale, this model consists only valence quarks, no sea quark or gluon contributes. The behavior of the pion PDF at large x is still an unresolved issue. However, our observation at large x agrees with perturbative QCD, where the behavior of the PDF has been predicted to be (1 − x) 2 [98][99][100][101], a behavior further supported by the Bethe-Salpeter equation (BSE) approach [102,103].

V GENERALIZED PARTON DISTRIBUTIONS (GPDS)
We calculate the GPDs of quark for pion and kaon in LCQM. GPDs have support region x [−1, 1] [11,146]. However, for present calculations we restrict ourself to only DGLAP region i.e. ζ < x < 1. One can define the correlation to evaluate unpolarized GPD H(x, ζ = 0, t) as [11] H P (x, 0, t) = 1 2 By inserting the initial and the final states of pion (π + (P ) and π + (P )) and kaon (K + (P ) and K + (P )) from Eq. (18) in above equation, we obtain the unpolarized quark GPD H(x, 0, t) in overlap form of LCWFs as where the final state struck quark momentum is written as By substituting the respective wave functions for pion and kaon from Eqs. (6), (12), (16) and (17), we get the explicit expressions of GPDs. For the case of pion we have with in initial and final states, respectively.
For the case of kaon we have with in the initial and the final states, respectively. Here, t = −q 2 ⊥ is denoted as the total momentum transferred to the meson. The detailed discussion on graphical representation of GPD in case of kaon has been already explained in Ref. [151].
We use the parameters mentioned in Table I to calculate the GPDs H(x, 0, t) of u-quark in pion (π + meson).
To understand the dependence of the u-quark GPD on x and −t, in Figs. 5(a) and 5(b), we give the 3-D graphical representation of H in case of pion and kaon respectively. The quark distribution in pion with respect to the longitudinal momentum fraction x is maximum at the central value (x = 0.5) when the momentum transferred to the pion is zero. As the value of momentum transferred −t increases, the peak shifts towards higher values of x and the magnitude of distribution becomes lower. Unlike pion GPD, the kaon GPD has the maximum at lower x (< 0.5) when t = 0. This is because of the presence of strange quark having larger mass. While with increasing −t, the peaks along x get shifted to larger values of x same as pion GPD. This is a model indepent behavior of GPDs which has been observed in other phenomenological models for pion [105] as well as for nucleon [147][148][149][150].

VI TRANSVERSE MOMENTUM-DEPENDENT PARTON DISTRIBUTIONS (TMDS)
TMDs provide the distribution of partons in momentum space and are functions of longitudinal momentum fraction x = k + /P + and transverse momentum k ⊥ carried by the struck quark. To evaluate the pion and the kaon TMDs, the unintegrated quark-quark correlator can be defined as [108,152] Since the spin is zero for the pion and the kaon (S = 0), we are left with only the TMD which comes from the combination where meson is taken to be unpolarized. After putting the states of respective mesons and Γ = γ + , we get the explicit expressions of unpolarized pion TMD f π 1 (x, k 2 ⊥ ) and unpolarized kaon TMD f K 1 (x, k 2 ⊥ ). The overlap form of unpolarized TMD f 1 (x, k 2 ⊥ ) is given as For pion, the TMD evaluated is given as and for kaon, we have In Figs. 6(a) and 6(b), to get the combined information, we show the 3-D picture of distribution xf 1 of a unpolarized quark in the unpolarized pion and kaon with respect to the longitudinal momentum fraction and squared of the quark transverse momentum respectively. The probability of finding the quark in pion is more as compared to kaon, if the momentum fraction carried by that quark is higher in the longitudinal direction. As we increase the k 2 ⊥ , the momentum distribution starts lowering down in both cases. The distribution decreases when the transverse momentum carried by quark increases. The probability to find the quark in pion and kaon starts decreasing and eventually becomes zero by the increase in the quark transverse momentum.
In Figs. 7 and 8, we show the unpolarized pion and kaon TMD f 1 (x, k ⊥ ) in momentum plane (k x , k y ) by choosing different values of x. We observe that the symmetric distribution peak appear to be narrow in transverse momentum plane, while having the lower quark momentum fraction in case of pion and kaon.

(A) TMD evolution
The unpolarized TMD evolution is factorized in the framework of Collins-Soper-Sterman (CSS) formalism [153]. It includes the perturbative affects in the larger energies and momentum transfer regimes, in addition to the non-perturbative affects and the regime of low transverse momentum. The unpolarized TMD scale evolution is executed in b ⊥ space which can be done by taking the Fourier transformation of f 1 (x, k 2 ⊥ ) [111,[154][155][156][157]. We havef The TMD evolution off 1 (x, b 2 ⊥ ) is given as where the µ 2 -evolution operation is directed by the nonperturbative Sudakov factor g k (b ⊥ ) and TMD evolution factor R(µ, µ 0 , b ⊥ ). We have where g 2 is a free parameter and can be extracted from the experimental data. For the present work, the parameter g 2 has been taken from Ref. [158] and is given as g 2 = 0.13 GeV 2 . Further, we have with γ K and γ F being the anomalous dimensions, which are given as [153] and The parameter µ b used in Eq. (46) can be expressed in terms of parton impact parameter b ⊥ as Here C 1 is a constant and we choose its value as C 1 = 2e −γ E [153] with γ E = 0.577 being the Euler constant. In Fig. 9(a), we provide the graphical representation of unpolarized pion and kaon TMD evolution k ⊥ f 1 w.r.t x by choosing the different scale values µ 2 = 1 and 10 GeV 2 . It can be clearly seen that when the TMD evolution is implemented, the distribution peaks become broader and the magnitude goes on decreasing with increase in µ 2 . It would be important to mention here that for the case of pion and kaon, the TMD is evolved from the model scale µ 2 0 = 0.246 GeV 2 . There is negligible effect of evolution at µ = µ 0 , it simply provide the TMD function multiplied by k ⊥ . With growing µ 2 , the width of TMD peaks increases and its values experiences the rapid decrease in magnitude. In Fig. 9 (b), we observe the broad peak in case of pion as compared to kaon. With the evolution of f 1 , the magnitude of distribution decreases. Further, we see the clear asymmetry in case of kaon because of heavy spectator antiquark mass.

VII CONCLUSION
We present various quark distributions in the pion and the kaon in light-cone quark model for the valence quarks suitable for low-resolution properties. The lightcone wave functions in this model have been obtained by transforming the instant-form wave functions through the Melosh-Wigner rotation. We have obtained reasonable agreement with the experimental data for the pion DA which is also very close to asymptotic DA after LO QCD evolution following ERBL equation. Due to unequal quark masses, we observed distinctly different behavior of kaon DA from the asymptotic DA. The initial scale PDFs have been evaluated using the overlaps of the LCWFs. We then applied QCD evolution to our initial pion PDF in order to incorporate degrees of freedom relevant to higher-resolution probes which allows us to compare our QCD-evolved PDF with experimental data. The pion PDF at higher scale relevant to E-615 experiment has been computed based on the NNLO DGLAP equations. The initial low-resolution scale is the only adjustable parameter involved in QCD scale evolution. A good agreement with the reanalysis E-615 data has been observed when we evolved the pion PDF from the initial scale µ 2 0 = 0.246 GeV 2 . Further, we have evaluated GPDs in DGLAP region for zero skewedness i.e. 0 < x < 1, which provide us three-dimensional structure of hadron. For both the pseudoscalar mesons, depending upon the total momentum transferred to the composite system, we observe the change in distribution with respect to active quark longitudinal momentum fraction. The transverse structure of pion and kaon has also been examined. To observe the combined effect, we have presented the 3D picture of TMD with respect to x and k 2 ⊥ . We observed that as the active quark carries larger longitudinal momentum, the broadening in the transverse momentum plane also increases. Furthermore, we have presented the effect of µ 2 dependence on unpolarized pion and kaon TMDs. We have observed that the magnitudes of the distributions decrease and became wider as µ 2 increases. the black and red curves represent the unpolarized pion TMD k ⊥ f π 1 and the unpolarized kaon TMD k ⊥ f K 1 respectively, evolved from initial scale µ 2 0 = 0.246 GeV 2 to different scales µ 2 : µ 2 = 1 GeV 2 (solid curves) and µ 2 = 10 GeV 2 (dotdashed curves) at fixed x = 0.5 w.r.t. k ⊥ . (b) Right panel : the black and red curves represent the unpolarized pion TMD f π 1 and the unpolarized kaon TMD f K 1 evolved at different scales µ 2 : µ 2 = 1 GeV 2 (solid curves) and µ 2 = 10 GeV 2 (dotdashed curves) at fixed k ⊥ = 0.5 GeV w.r.t. x.