The transverse polarization of Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document} hyperons in e+e-→Λ↑hX\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow \Lambda ^\uparrow h X$$\end{document} processes within TMD factorization

We investigate the transverse polarization of the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document} hyperon in the processes e+e-→Λ↑π±X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow \Lambda ^\uparrow \pi ^\pm X$$\end{document} and e+e-→Λ↑K±X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow \Lambda ^\uparrow K^\pm X$$\end{document} within the framework of the transverse momentum dependent (TMD) factorization. The transverse polarization is contributed by the convolution of the transversely polarizing fragmentation function (PFF) D1T⊥\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{1T}^\perp $$\end{document} of the lambda hyperon and the unpolarized fragmentation function D1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_1$$\end{document} of pion/kaon. We adopt the spectator diquark model result for D1T⊥\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{1T}^{\perp }$$\end{document} to numerically estimate the transverse polarization in e+e-→Λ↑hX\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow \Lambda ^\uparrow h X$$\end{document} process at the kinematical region of Belle Collaboration. To implement the TMD evolution formalism of the fragmentation functions, we apply two different parametrizations on the nonperturbative Sudakov form factors associated with the fragmentation functions of the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}, pion and kaon. It is found that our prediction on the polarization in the Λπ+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda \pi ^+$$\end{document} production and Λ¯π-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{\Lambda }} \pi ^-$$\end{document} is consistent with the recent Belle measurement in size and sign, while the model predictions on the polarizations in Λπ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda \pi ^-$$\end{document} and ΛK±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda K^\pm $$\end{document} productions show strong disagreement with the Belle data. The reason for the discrepancies is discussed and possible approaches to improve the calculation in the future are also discussed.


Introduction
Understanding the internal parton structure of hadrons and the fragmentation mechanism of the parton into hadrons is one of the main goals in QCD and high energy physics. Particularly, the production of a polarized hyperon from unpolarized pp collisions has been observed [1,2] and it formed a long-standing challenge [3,4] in hadron physics and spin physics. It is suggested [5] that a polarizing fragmentation a e-mail: xiaoyuwang@zzu.edu.cn b e-mail: zhunlu@seu.edu.cn (corresponding author) function (PFF) [6], denoted by D ⊥ 1T , can account for the polarization of the production. As a time-reversal-odd and transverse momentum dependent (TMD) fragmentation function, D ⊥ 1T describes the fragmentation of an unpolarized quark to a transversely polarized hadron, and it is usually viewed as the analog of the Sivers function [7,8] which gives the azimuthal asymmetry in the distribution of unpolarized quarks inside a transversely polarized nucleon. Furthermore, D ⊥ 1T may play an important role in the spontaneous polarization, such as: q → ↑ X [9]. Thus, the study of the polarized production could also provide the information on the spin structure of the hyperon. This is intriguing since the hyperon can not serve as a target in high energy scattering processes.
As D ⊥ 1T is a chiral-even function, in principle it may be accessed directly without any unknown, chiral-odd, counterparts. However, significant signal on the transverse polarization of the hyperon has not been observed in the single inclusive e + e − annihilation (SIA) experiment performed by OPAL at LEP [10]. Alternative to SIA, e + e − → ↑ + h + X [11][12][13] and semi-inclusive deep inelastic scattering (SIDIS) p → + ↑ + X were also suggested [11] to study the polarization. Those measurements could provide further understanding of the origin of the sizable transverse polarisation of hyperons observed in different processes [2,5,6,[14][15][16][17]. Recently, nonzero transverse polarization of production in SIA and semi-inclusive e + e − → (¯ ) + K ± (π ± ) + X process was measured by the Belle Collaboration [13], making the extraction [18,19] of the polarized fragmentation function of the possible. It is worth pointing out that, since there are no hadrons in the initial state, electronpositron annihilation is recognized as a rather clean process to access TMD fragmentation functions. Motivated by the recent Belle data, the authors in Ref. [20] applied an approach that assumes a kinematical configuration between the PFF and the unpolarized fragmentation function, which allows the comparison between numerical calculations and experimental data in the SIA process. The comparison showed the estimate results are in agreement with the data. On the other hand, model calculations may also provide an approach to acquire knowledge of the transverse polarization of the hyperon. A calculation of D ⊥ 1T for light flavors based on a spectator-diquark model has been performed in Ref. [21] and was used to make predictions on physical observables.
The main purpose of this work is to apply the TMD factorization (P hT /z h Q) [22][23][24][25][26] to estimate the transverse polarization in the transversely polarized production e + e − → (¯ ) + π ± (K ± ) + X process. In the last decades, TMD factorization has been applied in various high energy processes [26][27][28][29][30][31][32][33][34][35][36]. In the TMD formalism, the differential cross section in the region P hT /z h Q can be expressed as the convolution of the hard scattering factors and the well-defined TMD distributions or fragmentation functions. The TMD formalism embeds evolution information of those functions, of which the energy evolution (or the scale dependence) are governed by the so-called Collins-Soper equation [22,23,26,37]. The solution of the evolution equation indicates that the changes of TMDs from a initial scale to another scale may be determined by an exponential form of the Sudakov-like form factor [23,26,29,38], which can be separated to the perturbative part and nonperturbative part. The former one is perturbatively calculable, while the later one can not be calculated directly and is usually obtained by phenomenological extraction from experimental data. In Refs. [29,33,35,39] the authors extracted the nonperturbative Sudakov form factor corresponding to the unpolarized fragmentation function. In the e + e − → ↑ h ± X process, the convolution of the PFF for the hyperon and the unpolarized fragmentation function for pion/kaon can give rise to the transverse polarization of hyperon. Particularly, we will take into account the TMD evolution for both the PFF D ⊥ 1T and the unpolarized fragmentation function D 1 . As a comparison, we adopt two different parameterizations on the nonperturbative part for the fragmentation functions [33,35].
The remaining content of the paper is organized as follows. In Sect. 2, we present the formalism of the transverse polarization contributed by the convolution of D ⊥ 1T and the unpolarized fragmentation function D 1 in the e + e − → ↑ +h+X process within TMD factorization. In Sect. 3, we investigate the evolution effect for the unpolarized and transversely polarized fragmentation functions at leading order and present our choice on the nonperturbative Sudakov form factors associated with the fragmentation functions in details. In Sect. 4, we numerically estimate the transverse polarization at the energy √ s = 10.58 GeV which is accessible at Belle. We also compare the results calculated from different choices of the nonperturbative ingredients associated with the TMD evolution. Finally, We summarize the paper in Sect. 5.
2 Transverse polarization in e + e − → ↑ + h + X process In this section, we will present the detailed framework of transverse polarization in e + e − annihilation process with a transversely polarized hyperon and a light meson h (pion or kaon) produced in the final state In the e + e − annihilation process, the electron and the positron annihilate into a virtual photon, which then decays into a quark-antiquark pair, the latter one fragments into final-state transversely polarized hyperon and the light meson h (pion or kaon in this work), respectively. and stand for the four-momenta of the initial-state positron and electron, q is the four-momentum of the virtual photon with q = + , which is time-like (q 2 = Q 2 > 0). k and p are the four-momenta of the quark and antiquark, respectively, while K and P are the four-momenta of the final-state and pion/kaon meson. The center of mass energy for the process can be written as s = ( + ) 2 = Q 2 .
In the ideal case, the transversely polarized and the pion/kaon meson should be produced completely back-toback. However, the radiation of the gluon in the fragmentation process and the transverse momentum dependence make the hadrons deviate from the ideal back-to-back state. The TMD factorization can be used to describe the imbalance from the back-to-back state as well as to calculate the differential cross section. There are two experimental methods to define the reference frame in the e + e − annihilation process in literature [27,[40][41][42]. In this work, we adopt the secondhadron momentum frame, which means that the momentum direction of the second hadron-pion/kaon meson is defined as the z axis. Thexz plane is determined by the lepton and the pion/kaon meson momentum directions, while the hadron plane is determined by z axis and the momentum direction of the hadron. Hence φ is the azimuthal angle of the hadron plane relative to the lepton plane, while φ S represents the azimuthal angle of the hyperon polarization vector S ⊥ in the lepton frame. k T and p T are the transverse momenta of the two fragmenting quarks, which are related to the transverse momenta of the final hadrons through K ⊥ = −z 1 k T and P ⊥ = −z 2 p T . Finally, P ⊥ = −z 1 q T is the transverse momentum of hadron with respect to the transverse momentum of the virtual photon q T , which is also recognised as the transverse momentum of the process. The invariants z 1 and z 2 are defined as z 1 = 2K ·q Q 2 and z 2 = 2P·q Q 2 , which can be identified as the light-cone momentum fraction (neglecting term of order (q 2 T /(Q 2 ))) of the and the pion/kaon meson in the fragmenting quark and antiquark, respectively.
The transverse polarization P T in e + e − → ↑ h X can be defined as [19,20] where d σ = 1 2 (σ (S ⊥ ) − σ (−S ⊥ )), and dσ in the denominator is the unpolarized differential cross section. The transverse-spin dependent differential cross section can be expressed as [19] dσ (S ⊥ ) in which the azimuthal angle of the lepton plane φ l has been integrated out. Here, M is the mass of the , and the unit vectorĥ is defined asĥ = P ⊥ | P ⊥ | = q T q T [28,43]. The notation F represents the convolution of the corresponding fragmentation functions in the transverse momentum space with ω( p T , k T ) an arbitrary function of p T and k T . Since it is convenient to deal with the TMD evolution effect in the b space, which is conjugated to the k T space, we perform the Fourier transformation for the delta function to obtain the denominator in Eq. (2) as where J 0 is the Bessel function of the zeroth order. Thus, the unpolarized fragmentation function in b space (hereafter the tilde terms represent the ones in b space, with b=|b ⊥ |) is the Fourier transformation of the fragmentation function in with P ⊥ the transverse momentum of the hadron with respect to the fragmenting quark, and P ⊥ = |P ⊥ |. Similarly, the numerator can be expressed as: where the PFF of the hyperon in b space is defined as The energy dependence of the fragmentation functions will be discussed in details in the following section.

The TMD evolution of fragmentation functions
In this section, we will set up the formalism of the TMD evolution for both the PFF and the unpolarized fragmentation function.
The key purpose of the TMD evolution is to deal with the energy dependence of the TMD fragmentation functions in the b space. The advantage of the b space is that the differential cross section can be written as a simple product instead of complicated convolution of the fragmentation functions in the transverse momentum space. There are two energy dependencies of the TMD fragmentation functions D(z, b; μ, ζ D ) in b space, one is the renormalization scale μ related to the collinear fragmentation functions, and the other one ζ D is the scale related to the cutoff in the operator definition of the TMD function to regularize the singularity [22,23,26,29,31,39]. Hereafter, we set μ = √ ζ D = Q for simplicity, and the TMD fragmentation functions can be written asD(z, b; Q) (Here, D is a shorthand for any fragmentation function, such as the PFF D ⊥ 1T and the unpolarized fragmentation function D 1 ). Therefore, the main focus turns to solve the energy evolution equations ofD(z, b; Q) [26,44] defined in Eqs. (7) and (9).

Solving the TMD evolution equation for the fragmentation function in b-space
The energy dependence for the ζ D is encoded in the Collins- while the μ dependence is derived from the renormalization group equation as withK the evolution kernel, and γ K , γ D the anomalous dimensions. The solutions of these evolution equations were studied in details in Refs. [26,37,45]. Here, we will only discuss the final result. After solving the above evolution equations, the overall structure of the solutions are identical to each other in all TMD factorization schemes, and the evolution effects are incorporated into the exponential form factors [22][23][24][25][26]45] as Here, D(Q) is the hard scattering factor which can be calculated through perturbative QCD, and S(Q, b) is the Sudakovlike form factor. Eq. (13) shows that the TMD fragmentation functions D at an arbitrary scale Q can be evolved from an initial scale μ i through the evolution encoded by the exponential form exp(−S (Q, b)).
Studying the behavior of the TMD fragmentation function in the b space is quite important since it can determine the transverse momentum dependence of the physical observables through the inverse Fourier transformation. To do this, one need the information of S(Q, b) in the entire b region. In the small b region 1/Q b 1/ , the b-dependence is perturbative, while it turns to be non-perturbative in the large b region. Thus, to combine the information between the two regions, a matching procedure should be adopted. In the original CSS approach [23,[46][47][48], a parameter b max is introduced as the boundary between the two regions, which allows a smooth transition from perturbative region to nonperturbative region as well as to avoid hitting on the Landau is always in the perturbative region. In the original CSS approach it has the following form [23,34] Apart from the above form, there are also several different choices on the form of b * (b) in literature [35,49].

Sudakov form factor
After introducing the b max and the b * prescription, the Sudakov form factor can be separated into two parts: the perturbative part S P and the non-perturbative part S NP . It is important to keep in mind thatD (z, b; Q) follow exactly the same QCD evolution effect in the perturbative region, which means that S P (Q, b) is universal and is the same for different kinds of fragmentation functions (namely, S P is spin-independent). The details of S P have been studied in literature, it is usually written in the following general form [31,33,39,50,51]: where the coefficients A and B in Eq. (15) can be expanded as the series of α s /π : In this work, we adopt A (n) up to A (2) and B (n) up to B (1) in the accuracy of next-to-leading-logarithmic (NLL) order [23,29,31,46,48,50]: with C F = 4/3, C A = 3 and T R = 1/2. On the other hand, the non-perturbative part of the Sudakov form factor S NP can not be calculated perturbatively, it is usually parameterized and extracted from experimental data. There are several different approaches to parameterize S NP , we will discuss two of them in details.
One of the widely used non-perturbative Sudakov form factor S NP associated with an unpolarized TMDFFs has the following form (Echevarria-Idilbi-Kang-Vitev (EIKV) parametrization) [48,[52][53][54]: Since the information of the nonperturbative Sudakov form factor associated with the polarized fragmentation function of the hyperon still remains unknown, we assume it to be the same as the one for the unpolarized fragmentation NP . In Eq. (21), g i (b) are functions of the impact parameter b. Particularly, g 2 (b) contains the information on the large b behavior of the evolution kernelK , while g ff 1 contains the information about the intrinsic nonperturbative transverse motion of bound partons, i.e., it depends on the type of the hadron and quark flavor. It might also depend on the momentum fraction of the hadrons in the fragmenting quark z [55]. One should note that g 2 (b) is universal among different types of TMDs and does not depend on the particular process, which is one of the important predictions of TMD factorization [26, 29,33,34]. On the other hand, the parameter g ff 1 was assumed to be related to the intrinsic transverse momentum squared p 2 T Q 0 for TMD fragmentation functions at the initial energy scale Q 0 as In Ref. [33], the authors show that the Sudakov factor with the following parameters (and Q 0 = √ 2.4 GeV) leads to a reasonably good description of all experimental data on SIDIS, DY lepton pair and W/Z boson production, Besides the EIKV parametrization, several other forms for S NP have been also proposed [26, 29,33,39,50,51] recently. Particularly, we will discuss the evolution formalism from Ref. [35], which is called as Bacchetta-Delcarro-Pisano-Radici-Signori (BDPRS) parametrization. In this approach, the fragmentation functions evolving from the initial energy scale to the final energy scale also has the exponential form as where g K = −g 2 b 2 /2, following the choice in Refs. [48,52,56].D a→h  1NP (z, b 2 ) is the intrinsic nonperturbative part of the fragmentation function parameterized as with β, γ , δ and N 3,4 ≡ g 3,4 (ẑ) withẑ = 0.5 are free parameters fitted to the available data from SIDIS, Drell-Yan, and Z boson production processes yielding β = 1.65, γ = 2.28, δ = 0.14, λ F = 5.50 GeV −2 , g 2 = 0.13 GeV 2 , N 3 = 0.21 GeV 2 , N 4 = 0.03 GeV 2 . Furthermore, in Ref. [35], the new b * prescription different from Eq. (14) was also proposed as Again, b max is the boundary of the nonperturbative and perturbative b-space region fixed by b max = 2e −γ E GeV −1 ≈ 1.123 GeV −1 , with γ E ≈ 0.577 the Euler's constant [22]. Besides, the authors in Ref. [35] also chose to saturate b * at the minimum value b min ∝ 2e −γ E /Q. In this work, we will adopt the both the EIKV evolution formalism and the BDPRS evolution formalism to calculate the transverse polarization of to investigate the impact of the different evolution formalisms on the polarization.

The evolved fragmentation functions
In the perturbative region 1/Q b 1/ , the TMD fragmentation functions can be expressed as the convolution of the perturbatively calculable coefficients and the corresponding collinear counterparts of the TMD fragmentation functions at the fixed energy scale μ b , which is a dynamic scale q←i (α s /π ) n is the perturbatively calculable coefficient function with i summing over the quark and antiquark flavors. Here, we will adopt the leading order (LO) result, i.e. C (0) q←i = δ iq δ (1−z). In other words, Here, D 1T [57] as: where D

⊥(1)
1T is the first transverse moment of D ⊥ /q 1T . It is straightforward to rewrite the scale-dependent TMD fragmentation functions of the pion/kaon meson and the hyperon in b spacẽ The factor of 1 2 in front of S P comes from the fact that S P of quarks and antiquarks satisfies the relation [58] With all the above ingredients, we can explicitly write out the evolved TMD fragmentation functions as Thus, the fragmentation function in the transverse momentum space can be obtained by performing the inverse Fourier transformation 1T .
with J 1 being the Bessel function of the first order.

Numerical calculation
In this section, we numerically estimate the transverse polarization in the process e + e − → (¯ ) + K ± (π ± ) + X using the framework set up above and compare the estimation with recent experimental data measured by the Belle Collaboration [13]. In order to obtain the numerical results, one needs to utilize the corresponding collinear parts of the TMD fragmentation functions as the inputs of the TMD evolution effects. For the unpolarized collinear fragmentation function D 1 (z) for pion and kaon, we adopt the leading order set of DSS parametrization [59].
The other important input is the twist-3 collinear correlation function D 1T , for which we apply the spectator diquark model results [21], where the contributions from both the scalar diquark and the axial-vector diquark spectators are included. Assuming the SU(6) spin-flavor symmetry, the fragmentation functions of the hyperon for light flavors satisfy the relations between different quark flavors and diquark types where u, d and s denote the up, down and strange quarks, respectively. One should notice that in this model only the valence quarks contribute to the fragmentation function, while the sea quark contribution is zero. D and D ⊥(s) 1T represent the contribution from the axial-vector diquark and scalar diquark, and have the form 1 where k 2 can be written as z , with m q , m D , M the masses of the parent quark, the spectator diquark and fragmenting hyperon, respectively. The 2 has the general form 2 = λ 2 z α (1 − z) β . At one loop level, there are four diagrams that can generate imaginary phases contributing to the nonzero polarized fragmentation function. For the scalar case, the four nonzero contributions are while for the axial-vector case, the four contributions are 2 .
Here, the expressions for the functions of A, B, C P − , D P − , I i have been given in Ref. [21].
where the values of the last four parameters are fixed, and the coupling constant α s is chosen as 0.817 at the model scale. Figure 1  1T (z, z, μ b ) at the model scale can be obtained through Eq. (32). For consistency, we apply the unpolarized fragmentation function of the hyperon D /q 1 (z) calculated from the same model. Furthermore, to solve the energy dependence of the collinear counterparts for the fragmentation functions, we apply the QCDNUM evolution package [61] to evolve the unpolarized fragmentation function D 1T (z, z, μ b ) from Ref. [62], which has the same form as that for the unpolarized fragmentation function. When applying the model results forD ⊥(3) 1T and D 1 as the nonperturbative input to calculating the polarization, some attention should be paid, since the results are obtained as a low scale μ 2 0 = 0.23 GeV 2 , which is beyond the limit of applicability of the perturbative QCD evolution equations. This is the common feature of many model calculations for TMD distributions and fragmentation functions because these models are usually valid at a hadronic scale. As the evolutions forD and D 1 in the nonperturbative region are unknown, we assume they follow the same evolutions as those in the perturbative region. This is of course an approximation to the real situation. However, we expect the assumption will not qualitative change our result since P T is the ratio between the polarized cross section and the unpolarized cross section, so that there is some cancelation of the effects between the numerator and the denominator.
Note that the z-variable defined in the Belle experiments is the energy fraction of the hadron with respect to the fragmenting quark, denoted by z h = 2E h / √ s (h = , π or K ), while our analysis in the previous sections and the input collinear functionD 1T (z, z, μ b ) apply the light-cone momentum fraction of the hadron for z. The two kinds of variables are connected by To achieve a better comparison between our calculation and the Belle data, we apply the transformation in Eq. (45) to change the light-cone momentum fraction z to the energy fraction z h and calculate the polarization P T . The Belle Collaboration measured the data in z bins with boundaries at z = 0.2, 0.3, 0.4, 0.5, 0.9 [13]. In Fig. 2, we plot the numerical results of the transverse polarization in the process e + e − → ↑ π + X at Q = 10.58 GeV as functions of z π for different z bins. Since the positivity bound of D ⊥ 1T in the model calculation was violated at large z region (z > 0.75) [21], the z bins are shown in Fig. 2 up to 0.5, i.e. the analysis is performed in z 1 bins with boundaries at z = 0.2, 0.3, 0.4, 0.5. To make the TMD factorization valid in the kinematic region P h⊥ /z h Q, the integration over the transverse momentum P h⊥ is performed in the region of 0 < P h⊥ < 0.7 GeV. The solid lines show the transverse polarization calculated from the BDPRS evolution formalism [35] and the b * prescription in Eq. (27), while the dashed lines show the transverse polarization using the EIKV parametrization [33] for the nonperturbative Sudakov form factor for comparison, with the b * prescription in Eq. (14). The solid squares show the experimental data measured by the Belle Collaboration [13], with the error bars including both the systematic error and the statistical error.
As shown in Fig. 2, the transverse polarization of in the process e + e − → ↑ π + X is negative, in agreement with most of the Belle data. In this case, the polarization is dominated by the term D , as in our model there is no sea contribution to fragmentation functions and D is negative in our model lead to the prediction of negative polarization in π + production . Our estimates also show that the size of the polarization increases with increasing z π , and the results from the EIKV parameterization are close to those from the BDPRS parametrization. We can conclude that different methods dealing with the non-perturbative evolution lead to the similar results for the polarization P T . We also find that using different definition for z (momentum fraction or energy fraction) will slightly reduce the size of the numerical result. In our calculation, using the energy fraction z h can describe the Belle data better than using the momentum fraction z.
In Fig. 3, we plot the the transverse polarization of vs z π in different z bins, but in the process e + e − → ↑ π − X . In this case the polarization from our model calculation is firmly negative in each z bin. The reason is that in ↑ π − production our model indicates the term D This result is similar to our model prediction for the polarization in π + production, since in our model D ⊥u 1T = D ⊥d 1T and D ⊥s 1T is much smaller than D ⊥u 1T , as shown in Fig. 1. However, the polarization from the Belle data for ↑ π − production is positive in the 0.2 < z < 0.3 bin, and turn to be negative in the 0.4 < z < 0.5 region with relatively smaller size. In the intermediate bin 0.3 < z < 0.4, the polarization is consistent with zero. This shows that our calculation for ↑ π − production strongly disagrees with the Belle data. In addition, the charge-conjugation symmetry implies P T (¯ h − ) = P T ( h + ). In these respects, the calculated results for P T (¯ π − ) are consistent with experimental data of Belle. Fig. 2 The transverse polarization as function of z π calculated in the process e + e − → ↑ + π + + X . The red solid lines correspond to the results from the BDPRS parametrization [35] [Eqs. (24) and (25)] on the nonperturbative form factor while the dashed lines correspond to the results calculated from the EIKV parametrization [33] [Eq. (21)] on the nonperturbative Sudakov form factor. The solid squares represent the Belle data for comparison Finally, we also calculate the transverse polarization of in the processes e + e − → ↑ K + X and e + e − → ↑ K − X . The results are shown in Figs. 4 and 5, respectively. We conclude that the overall signs of the polarization in both the ↑ K + and ↑ K − productions are negative, consistently with the sign of the Belle data for the ↑ K + production and the ↑ K − production in the lower z K region. However, our predicted polarization in the ↑ K + production is rather small, which is significantly different from the Belle data; while our prediction in the ↑ K − production is several times larger than the Belle data. The reason for a poor model result in ↑ K + production is because the D ⊥ /s 1T ⊗ D K + /s 1 term provides the main contributions, while in our model D ⊥ /s 1T is much smaller than that of the u and d quarks. In the case of ↑ K − production, the dominant contribution is the D term, leading to a rather large polarization. In addition, the charge-conjugation symmetry implies that P T (¯ K − ) = P ( K + ) and P T (¯ K + ) = P T ( K − ).
Some comments are in order. Although our prediction for the polarization in ↑ π + production agrees with the Belle data, there are large discrepancies between the Belle data and our predictions for ↑ π − production as well as ↑ K ± production. From the available parameterizations on D ⊥ 1T from Refs. [18] and [19] which can satisfactorily describe the polarizations in ↑ π ± and ↑ K ± productions, one can find that a positive D and substantial sea contribution are crucial for describing the Belle data in all cases. For example, positive polarizations in ↑ π − production requires a positive D ⊥ /u 1T ; a negative sea PFF D ⊥ /ū 1T will lead to a negative nonzero polarization in ↑ K + [18,19]. These features are different from our model input, in which D ⊥ /u 1T is negative and is equal to D comes from the assumption of SU (6) symmetry in the model calculation of the fragmentation function. Besides, the spectator model does not incorporate sea contents of the fragmentation functions. For a further study, a model that can allocate substantial SU(6) breaking effect may lead to very different result for the PFFs of the u and d quarks, particularly in sign. This will improve the model description of the polarization in ↑ π − production at Belle. Incorporating the sea quarks in the model could also provide a better description on the ↑ K ± production. Finally, in our calculation we have not considered the Y -term in the cross section. The inclusion of the Y -term could change the size of the Polarization in the large momentum region.

Conclusion
The e + e − annihilation with a transversely polarized hyperon produced in the final state is a useful tool to study the hyperon spin structure and the non-perturbative fragmenta-tion mechanism. In this work, we have studied the transverse polarization of in the processes e + e − → (¯ ) + π ± + X and e + e − → (¯ )+ K ± + X at Belle by applying the TMD factorization formalism. The polarization was contributed by the convolution of the T-odd PFF D ⊥ 1T , which describes the fragmentation of an unpolarized quark into a transversely , and the unpolarized fragmentation function of pion or kaon. To evolve the TMD fragmentations of the pion, kaon and from the initial energy to the experimental energy, we have taken into account the TMD evolution effects, in which the Sudakov form factor plays an important role. The Sudakov form factor is separated into perturbative part and non-perturbative part, for the former one, we adopted the results from the perturbative QCD at NLL accuracy, while for the latter one, we considered two different non-perturbative TMD evolution formalism for comparison. As the nonperturbative Sudakov form factor associated with the PFF is still unknown, we assume that it has the same form as that of the unpolarized fragmentation function. The hard coefficients associated with the corresponding collinear functions in the TMD evolution formalism are kept at the leading-order accuracy. For the fragmentation functions of , we have chosen the results from the diquark spectator model. It has been shown that different choices of nonperturbative Sudakov form factors in the TMD evolution formalism lead to similar results for transverse polarization of in process e + e − annihilation. Within the framework of TMD evolution, our prediction for the polarization in ↑ π + production agrees with the Belle data. However, there are large discrepancies between the Belle data and our predictions for ↑ π − production as well as ↑ K ± production. A comparison with our model input for the PFF with the available parameterizations Refs. [18] and [19] indicates that a positive D ⊥ /u 1T and sizable sea PFF are essential to describe the many aspects of the Belle data. Possible improvements of the model calculation include allocating substantial SU(6) breaking effect to lead to a opposite sign of the u quark PFF compared with the d quark PFF, incorporating the sea quarks in the model, and taking into the Y -term in the large momentum region. Further studies are needed to understand the transverse polarization of the Lambda production in e + e − annihilation from the model aspect.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .