Unpolarized transverse-momentum dependent distribution functions of a quark in a pion with Minkowskian dynamics

The unpolarized twist-2 (leading) and twist-3 (subleading), T-even, transverse-momentum dependent quark distributions in the pion are evaluated for the first time by using the actual solution of a dynamical equation in Minkowski space. The adopted theoretical framework is based on the 4D homogeneous Bethe–Salpeter integral equation with an interaction kernel given by a one-gluon exchange, featuring an extended quark-gluon vertex. The masses of quark and gluon as well as the interaction-vertex scale have been chosen in a range suggested by lattice-QCD calculations, and calibrated to reproduce both pion mass and decay constant. The sum rules to be fulfilled by the transverse-momentum dependent distributions are carefully investigated, particularly the leading-twist one, that has to match the collinear parton distribution function, and hence can be scrutinized in terms of existing data as well as theoretical predictions. Noteworthy, the joint use of the Fock expansion of the pion state facilitates an in-depth analysis of the content of the pion Bethe–Salpeter amplitude, allowing to calculate the gluon contribution to the quark average longitudinal fraction, that results to be ∼6%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 6\%$$\end{document}. The current analysis highlights the role of the gluon exchanges through quantitative analysis of collinear and transverse-momentum distributions, showing, e.g. for both leading and subleading-twists, an early departure from the widely adopted exponential fall-off, for |k⊥|2>m2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\textbf{k}_\perp |^2> m^2$$\end{document}, with the quark mass ∼ΛQCD\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim {\varLambda }_{QCD}$$\end{document}.

It should be recalled that the proton is the elective target of much experimental (see, e.g., Refs.[9][10][11]) and theoretical research (see, e.g., Refs.[48][49][50] and references therein).While the pion, given the experimental challenges its study poses, has surely attracted less efforts in spite of its intriguing double-nature, being both a Goldstone boson (and hence fundamental for investigating the dynamical chiral-symmetry breaking) and a quark-antiquark bound system (i.e. the simplest bound system in QCD).In particular, a first extraction of the pion unpolarized leading-twist TMD from Drell-Yan data can be found in Ref. [51], where the results of the E615 Collaboration [52] has been used, and in Ref. [53], where both the previous data and the E537 Collaboration cross-sections [54] have been included.As to the phenomenological calculations, a broad overview, embracing different approaches, can be gained from Refs.[35,38,40,41,[44][45][46][47][55][56][57][58] (see also Ref. [59] for the generalized TMDs in a spin-0 hadron).
As a conclusion to the above schematic introduction, it has to be emphasized that the vast amount of nowadays theoretical studies on TMDs finds its strong motivation in the very accurate measurements that will come from forthcoming electron-ion colliders, that promise to achieve greatly expected milestones in the experimental investigation of non-perturbative QCD, given the planned high energy and luminosity [60,61].
Our aim is to obtain, for the first time, T-even leading-and subleading-twist unpolarized TMDs (uT-MDs) of the pion, by solving a dynamical equation directly in Minkowski space, namely relying on a genuinely quantum-field theory framework based on the 4D homogeneous BSE [42,43].The homogeneous BSE is an integral equation and therefore suitable for dealing with the fundamentally non-perturbative nature of bound states.One should not get confused by the use of an interaction kernel expressed in a perturbative series, since an integral equation has a peculiar feature of infinitely many times iterating the boson exchanges contained in each term of the kernel, just what one needs for obtaining a pole in the relevant Green's function.In our approach (see Ref. [62] for details and references therein), based on the 4D homogeneous BSE in Minkowski space and the Nakanishi integral representation (NIR) of the BS-amplitude [63,64], the interaction kernel is given by the exchange of a massive vector boson in the Feynman gauge, with three input parameters, inferred from lattice QCD (LQCD) calculations (see, e.g.Refs.[65][66][67]): (i) the constituent-quark and gluon masses, and (ii) a scale parameter featuring the extended quark-gluon vertex.It should be pointed out that the ladder kernel, i.e. the first term in a perturbative series, can be a reliable approximation to evaluate the pion bound state, as suggested by the suppression of the non-planar contributions for N c = 3 within the BS approach in a scalar QCD model [68], and the presence of massive quarks and gluons, featuring the confinement effects in a relatively large system (r ch ∼ 0.66 fm).There is another important consequence stemming from the use of the BS-amplitude.Although in the definition of the q q-pair BS-amplitude there is a simple dependence upon two interacting fermionic fields, one ends up dealing with an infinite content of Fock states (the use of the Fock space allows one to recover a probabilistic language within the BS framework).In particular, by exploiting the Fock expansion of the pion state, one can establish a formal link between the LF-projected BS-amplitude (see, e.g., Refs.[69][70][71]), and the amplitude of the Fock component of the pion state with the lowest number of constituents.Therefore, in our approach, it is natural to call the LF-projected BS-amplitude LF valence wave function (LFWF), to be distinguished from the valence wave function, when a SU (3)-flavor language is adopted.In the latter case, the pion is composed by only two fermionic constituents, suitably dressed.One should keep in mind that within our framework, the pion LFWF contributes only with 70% [62] of the normalization, and consequently a significant role of the higher Fock components has to be highlighted, and possibly analyzed in-depth, as illustrated in what follows.Finally, we would emphasize that the first evaluation of the uTMDs strengthens the reliability of our approach and makes sound the ground for the next step, already in progress, i.e. taking into account the self-energy of the quarks (see Refs. [72] and [73,74]).
Indeed, in spirit, our approach is similar to the one developed in Ref. [45] for evaluating the leading-twist uTMD, where it was also taken into account the selfenergy of the quark propagator (solving the gap equation) and a confining interaction, but in Euclidean space.In this case, one resorts to a suitable method (based on the moments and a parametrization of the Euclidean BS-amplitude) to get the Minkowski-space distribution function.Differently, in our approach the NIR of the BS-amplitude allows one to successfully deal with the analytic structure of the BS-amplitude itself, obtaining an integral equation formally equivalent to the initial 4D homogeneous BSE, but more suitable for the numerical treatment.Many and relevant applications of our approach to the pion, such as the electromagnetic form factor [75], the PDF [76] and the 3D imaging [62], have confirmed its reliability and encouraged to broad the scope of our investigation.It should be pointed out that (it will become clear in what follows) the evaluation of quantities that depend not only upon the longitudinal dof but also the transverse ones leads to sharply increase the sensitivity to the dynamical content of a given phenomenological description of the pion, namely to increase its predictive power.Furthermore, the joint use of the Fock expansion, meaningful in the Minkowski space, allows one to resolve the gluonic content of the pion state.
The paper outline is as follows.In Sect.II, the general formalism and the notations are introduced, highlighting the ingredients of our dynamical approach, namely i) the Bethe-Salpeter amplitude, solution of the 4D homogeneous Bethe-Salpeter equation, and ii) the Nakanishi integral representation of the BSamplitude.In Sect.III, the expressions of leadingand subleading-twist uTMDs are given in terms of the Bethe-Salpeter amplitude of the pion.In Sec.IV and V, the leading and subleading-twist uTMDs are shown and compared with outcomes from other approaches.Finally, in Sect.IV, the conclusions are drawn, and the perspectives of our approach are presented.

The normalization of f
where F q π (t) is the quark contribution to the electromagnetic (em) form factor of the pion.The latter results to be equal to F π (t) = e q F q π (t) + e q F q π (t), with t = (P − P ) 2 , and is related to the matrix element of the four-current by N c P | ψq (0)γ µ ψ q (0)|P = 2P µ F π (t = 0).Finally, it should be pointed that inserting a complete basis in Eq. (1) and exploiting the good and bad components of the fermionic field one can easily demonstrate that f q 1 (γ, ξ) ≥ 0 (see Ref. [77]).
In order to describe the pion by taking into account at some extent the QCD dynamics in the nonperturbative regime, it is useful to resort to the Mandelstam framework [78], where the interacting quarkpion vertex is expressed in terms of the (reduced) BSamplitude, i.e. the solution of the 4D homogeneous BSE, and defined by where the fermionc field fulfills the Poincaré translation ψ(x) = e i P •x ψ(0)e −i P •x (recall that only the component P − is interacting in the LF dynamics, see, e.g., Ref. [79]).Thus, by using the Feynman-like diagrammatic picture inherent to the Mandelstam framework (see, e.g., Ref. [75] for the application to the em form factor), one can write the following expression for f q 1 (γ, ξ) where For the sake of completeness, let us write the BSE in ladder approximation, i.e. the one we are adopting for the numerical calculations, viz.
where quark and antiquark momenta are off-shell, i.e. p 2 q(q) = (±k + P 2 ) 2 = m 2 , and q = k − k is the gluon four-momentum.In Eq. ( 6), the fermion propagator, the gluon propagator in the Feynman gauge and the quark-gluon vertex, dressed through a simple form factor, are where g is the coupling constant, µ the mass of the exchanged vector-boson and Λ is a scale parameter, featuring the extension of the color distribution in the interaction vertex of the dressed constituents.Moreover, in Eq. ( 6), one has Γ ν (q) = C Γ T ν (q) C −1 , where C = iγ 2 γ 0 is the charge-conjugation operator.The normalization of the BS-amplitude reads (cf.Refs.[80] and [62] for details) The antiquark uTMD is given by where the minus sign results from the property of the normal-ordered em current to be odd under the action of the charge conjugation operator.It is noteworthy that in Appendix C 1, it is proven the identity of the normalization condition, Eq. ( 8), and the half sum of Eqs.
(1) and (9).Within a SU (3)-flavor symmetry framework, one describes a pion as a bound system of a massive q q pair.This leads to introduce the so-called valencequark PDF in the pion, that is charge symmetric (once the isospin breaking is disregarded [81]) as well as fulfills the charge conjugation.From those properties one deduces that the SU (3)-valence PDFs in the charged pions must verify: In our BS framework, in addition to the fermionic dof (still massive) one introduces also gluonic dof, by adding an explicit dynamical description of the binding.This amounts to the ladder exchange of infinite number of massive gluons.Therefore, at the initial scale, the quark and anti-quark longitudinal-momentum fraction distributions are not expected to be symmetric with respect to ξ = 1/2 (as it follows from the charge symmetry), given the gluon-momentum flow in the composite pion (see Sect.IV).The symmetric combination of quark and anti-quark contribution allows one to fulfill the charge symmetry, and hence it is relevant in the comparison with experimental data (see Ref. [76]).In what follows, in addition to the quark distributions, symmetric and anti-symmetric combinations are introduced for all the uTMDs we are going to analyze.
The half sum (difference) of the quark and antiquark contributions, Eqs.(1) and ( 9), yields the following charge-symmetric (anti-symmetric) expression for the leading-twist uTMD inside a π + meson f Analogously to Eq. (1), one can define the T-even subleading quark uTMDs, starting from the decomposition of the pion correlator [6,82].To be specific, one has two twist-3 uTMDs (see, e.g., Ref. [35] for the pion case) In analogy to Eq. ( 2), one gets for the twist-3 e q (ξ) (see Refs. [77,82] and Ref. [35] for the pion in phe- where the matrix element P | ψq (0) 1 ψ q (0)|P has to be proportional to the pion sigma term, once a QCD framework is adopted.As a matter of fact, one gets where m cur is the quark current mass and σ π is the pion sigma term, that becomes σ π = M/2, in the leading order of the chiral expansion, i.e. the Gell-Mann-Oakes-Renner relation [83].It should be pointed that recent LQCD calculations [84] confirm, with high accuracy, the Gell-Mann-Oakes-Renner relation in the range of the explored pion masses.Indeed, the QCD equations of motion gives a decomposition of the collinear PDF e(ξ) = dγ e(γ, ξ) in three terms.Among them, there is a singular term proportional to the pion sigma term, that reads (see, e.g., Ref. [85]) while the other two terms, one is due to quarkantiquark-gluon correlations and the other is proportional to the quark mass, do not contribute to Eq. ( 14) (see Ref. [85], where the issue is analyzed, taking the nucleon as actual case).In our phenomenological model the strength is distributed over the whole range of ξ (as in Ref. [35]), without the singularity at ξ = 0, as it will be shown in Sect.V.Moreover, one has for the first moment [85] 1 where the singular term and the gluonic contribution vanish, and only the term proportional to the quark mass contributes.
From the equations of motion of a free-quark model, one deduces the following relations between the above uTMDs (see,e.g., Ref [35,85,86]) where the uTMDs with a tilde are the gluonic contributions.The relevant point is the dependence of all the subleading-twist uTMDs from only the leading one, modulo the gluonic terms.In our fully interacting framework, one can anticipate that the relations are not recovered, and rather heavily broken.For a derivation of the first line of Eq. ( 17), fully consistent with QCD, one could apply the formalism presented in Ref. [85].Following Eq. ( 10), one readily writes down chargesymmetric and the anti-symmetric combinations for the subleading TMDs.One has to take care how the scalar and vector operators behave under the charge conjugation that impose a different combination of signs (cf.below Eq. ( 9)).Namely, one gets A. The BS-amplitude and its Nakanishi integral representation It is useful to briefly recall some features of our approach for obtaining the actual solution of the ladder BSE given in Eq. (6).The basic ingredient is the NIR of the BS-amplitude (see Ref. [64] for the general introduction, and Refs.[62,76,[87][88][89] for the application to a two-fermion case), but let us first introduce the general decomposition of the BS-amplitude, Φ(k; P ), for a 0 − bound state, viz.[87,90] Φ(k; P ) = S 1 (k; P )φ 1 (k; P ) + S 2 (k; P )φ 2 (k; P ) +S 3 (k; P )φ 3 (k; P ) + S 4 (k; P )φ 4 (k; P ) , where φ i 's are unknown scalar functions, that depend upon the kinematical scalars at disposal (k 2 , k • P and P 2 ), and S i 's are suitable Dirac structures, given by The functions φ i must be even for i = 1, 2, 4 and odd for i = 3, under the change k → −k, as dictated by the anti-commutation rules of the fermionic fields, and they can be written in terms of the NIR as follows where The real functions g i (γ , z ; κ 2 ), the unknowns of the problem under scrutiny, are the Nakanishi weight functions (NWFs), and assumed to be unique, following the uniqueness theorem from Ref. [64].The properties of the scalar functions φ i under the exchange k → −k translate to properties of the NWFs, but under the exchange z → −z .
Finally, it should be mentioned that NWFs are determined by solving a system of integral equation, so that one is able to non-perturbatively embed dynamical information that characterize the BS interaction kernel.The system of integral equations is formally deduced from the initial BSE, by exploiting the analytic structure of the scalar functions φ i , made explicit by means of the NIR.In fact, after inserting Eqs.(20) and (22) in the BSE, Eq. ( 6), and performing both the Dirac traces and a LF projection, i.e. the integration over the k − = k 0 − k 3 component of the relative momentum, one gets a coupled system of integral equations for the NWFs (see details in Ref. [89]).Once the NWFs are known, the BS-amplitude can be fully reconstructed through an inverse path, i.e.Eqs.(22) and (20).

III. THE UNPOLARIZED TMDS AND THE PION BS-AMPLITUDE
The evaluation of the leading-and subleading-twist uTMDs, given in Eqs. ( 10), ( 18) and ( 19), can be performed by inserting the decomposition of the BSamplitude in Eq. ( 20), obtaining where A new variable z is defined as z = −2k + /P + and the three quantities: i) the functions T i (γ, ξ), ii) the operators O i and iii) the phase η

S(AS) i
are given by Finally, the integrand F i j in Eq. ( 23) reads where a i j (k − , γ, ; S(AS)) are polynomial in k − (up to the cubic power) and can be found in Appendix A for each uTMDs, we are considering.
By exploiting the NIR, Eq. ( 22), one can perform the integration on k − .This integration amounts to restrict the LF-time to x + = 0, and it is also known as LF-projection (see, e.g., Refs.[70,71,91]).After carrying out the k − -integration, the expression of each T
The symmetric and anti-symmetric combinations of the T-even leading-twist uTMD, f S(AS) 1 (γ, ξ), allow us to address the evaluation of both quark and anti-quark contributions, f q(q) 1 (γ, ξ), that in the BS framework plus the Fock expansion of the pion state have interesting features, distinct from the ones of f S(AS) 1 (γ, ξ).After integrating the leading-twist f q(q) 1 (γ, ξ) on γ, one gets the quark PDF u q (ξ), while the symmetric combination provides the charge-symmetric PDF u S (ξ), i.e. the one is expected to have relevance at the valence scale (see, e.g., Ref. [81]).Indeed, in the Mandelstam approach the quark and antiquark PDFs do not have in general a symmetry with respect to ξ = 1/2, since each receives contributions from states containing an infinite number of gluons, as a consequence of the ladder-interaction kernel.But if we restrict to the contribution from the first Fock component in the expansion of the pion state, one gets the LF-valence u LF val (ξ), that is given by the BSamplitude projected onto the null plane [79] and is fully compliant with the charge symmetry (see below the discussion on the differences among u q (ξ), u S (ξ) and u LF val (ξ)).To illustrate general features and relations, in this Section we give some details, referring to Appendix C for a more complete discussion.
The symmetric and anti-symmetric leading-twist uTMDs, can be decomposed as follow where the non-vanishing symmetric contributions are given by Eqs.(C2), (C3), (C4) and I 3d (γ, ξ; S) = 0, respectively.The anti-symmetric quantities are shown in Eqs.(C5), (C6), (C7) and (C8), respectively.Two comments are in order.The symmetry properties of the above quantities with respect to the transformation z → −z are demonstrated in Appendix C, and can be translated into the symmetry with respect to ξ → 1 − ξ (that implements the charge-symmetry).A relevant feature is given by the presence in the expressions of I d,2d,3d of the partial derivatives ∂ n /∂z n , that should be considered dual of the n-th moment in k − of the relevant functions, generated by the formal step of the LF-projection (cf Eq. ( 26)).This is not a surprise since the variable z is proportional to k + .
A first consistency check of our formalism has been carried out in Appendix C 1, where it is shown that, within the Mandelstam approach, f S 1 (γ, ξ) and in turn f q 1 (γ, ξ) are normalized to 1, as naturally follows from the canonical BS-amplitude normalization [80,92], performed according to Eq. ( 8) (see also Ref. [62]).In particular, the integral on γ and ξ of I N (γ, ξ; S) saturates the normalization, while the other two terms provide vanishing contributions.Hence, one gets It should be recalled that all the calculated uTMDs vanish outside the interval 0 ≤ ξ ≤ 1, as dictated by the conservation of the plus components of the fourmomenta of both pion and constituents (cf.Eq. ( 10)).
It is understood that the integral of f AS 1 (γ, ξ) is vanishing, given the antisymmetry with respect to ξ → 1 − ξ.
with the normalization that follows from Eq. ( 29) and the vanishing result of the double integration of f AS 1 (γ, ξ).Finally, the quark and anti-quark PDFs are evaluated through with the normalization still given by Eq. ( 29).Within the SU(3)-flavor symmetry, one has to implement the charge symmetry (see, e.g.Ref. [81]) at the initial scale, and therefore u S (ξ) is the PDF to be compared, after the proper evolution, with the experimental data, as it has been shown in Ref. [76].
In the left panel of Fig. 1, u S (ξ) and its three contributions (see Eqs. (D4), (D5) and (D6)) are shown.The calculation has been carried out by adopting the BS-amplitude obtained by using the solution of the BSE as described in Ref. [62], using the following values of the three input parameters: m = 255 MeV, µ = 637.5 MeV and Λ = 306 MeV, able to reproduce the pion decay constant f P DG π = 130.50(1)(3)(13) MeV [93] (recall that the pion charge radius results to be r ch = 0.663 fm [75], in excellent agreement with r P DG ch = 0.659 ± 0.004 fm [94]).A remarkable cancellation among the contributions takes place, and this represents a common feature for all the integrated quantities generated by the uTMDs we are considering.In the right panel, one can see the comparison between the quark PDF, u S(AS) (ξ) and the LF-valence PDF, resulting from the one-to-one relation between the LF-projected BS amplitude and the valence amplitude of the Fock expansion of the pion state.In particular, the LF-valence PDF (see Refs. [62,76]), is given by where ξ = (1 − z)/2, ψ ↑↓ (γ, z) is the anti-aligned component of the LF-valence amplitude and ψ ↑↑ (γ, z) the aligned one (of purely relativistic nature having an orbital angular momentum equal to 1).These amplitudes are suitable combinations of the LF-projected scalar functions φ i (k; P ), Eq. ( 22).The integral on ξ of LF-valence PDF gives the probability of the valence state in the Fock expansion and amounts to The striking feature shown in the left panel is the shift toward low ξ of the quark PDF, so that for this quantity the symmetry ξ → 1 − ξ is slightly violated.

B. Analysing the shift and the gluon content
The PDF calculations based on the BS-amplitude are able to capture an explicit gluonic effect, to be taken distinct from the one responsible for the effective mass of the constituents.In particular, the difference between the two symmetric PDFs, i.e. u S (ξ) and u LF val (ξ) (recall that has P val = 0.7), can be traced back to the non negligible probability of the higher Fock states (HFS), where a q q pair interacts by exchanging any number of gluons.Interestingly, the difference can be effectively described only by a factor, since it turns out that u LF val (ξ)/P val largely overlaps u S (ξ).Finally, also the small, but relevant, shift of the quark PDF with respect to u S (ξ) has to be ascribed to the presence of HFS, as discussed in what follows.
To get a qualitative view, we remind that the pion state can be, in principle, decomposed in Fockcomponents, which are schematically written in ladder approximation as Due to the charge symmetry, each Fock-component is invariant by q ↔ q, and hence the valence state |q q provides a symmetric contribution to u q (ξ), identified with u LF val (ξ).The following terms contain gluons up to infinity.In our model, the gluon has an effective mass about twice the quark mass, so that the HFS cumulative effect results in a small shift of the u q (ξ) peak at ξ < 1/2, as shown in the right panel of Fig. 1.Actually, a similar effect, related to the increasing mass of the remnant, can be also recognized in the nucleon, where one has a valence parton distribution with a peak around 1/3 due to the presence of the other two constituent quarks.In the case of the pion, the effect is small since the valence component |q q has 70% of probability (as generated by our dynamical calculation), and hence is largely dominant.
To become more quantitative and illustrate this effect, we schematically write the quark PDF by using the Fock expansion of the pion state, Eq. (34), and inserting LF variables [79], one has where ξ 1(2) is the longitudinal-momentum fraction of the quark (antiquark) in each Fock state, composed by a q q pair and n − 2 gluons, generated by the iteration of the one-gluon exchange.Moreover, is the probability amplitude of the corresponding Fock component and fulfills a normalization condition that follows from the one of the pion state.In the n-th state one has Since ξ i > 0 for massive particles, the average value of ξ 1 starts to decreases while the number of gluons increases, as quantitatively shown in what follows.
Looking at the right panel of Fig. 1, one can realize that while the valence term, with probability P val = 0.7, has a peak at the HFS shift the peak to ξ 1 < 1/2, and decrease the tail, due to the constraint of the overall normalization.This is reflected in the evaluation of the first moment (recall ξ q ≡ ξ 1 ) where P n is the probability of the n-th Fock state beyond the valence one.The first term in Eq. ( 37) is equal to 0.35, since 1/2 is weighted by P val , and the rest is weighted by 0.3.Notice that for each HFS, normalized to 1, one has where the gluon bosonic nature leads to the factor n − 2.
The actual value of the first moment of u q (ξ) is that amounts to an average of ξ q HF S equal to 0.40.We can further analyse ξ q HF S , aiming at extracting a quantitative estimate of the exchanged-gluon contribution, ξ g .From the momentum sum rule Eq. (37), and recalling Eq. ( 36), we get where Dashed line: DSE calculation from Fig. 5 of Ref. [95].Dotted line: BLFQ result at 4.0 GeV [96].Shaded area: LQCD calculation extracted via Mellin moments from Ref. [97].Full squares: reanalyzed data by using the ratio between the fit 3 of Ref. [98], evolved to 5.2 GeV, and the experimental data [52], at each data point (see Ref. [76] for details).
Moreover, since each Fock component fulfills the charge symmetry, i.e. q ↔ q, the corresponding quark and antiquark momentum densities are equal and hence for the Mellin moments one has ξ k q HF S = ξ k q HF S (this property does not imply the charge symmetry of the total density, given the presence of the gluon contribution, cf.Eq. ( 37)).From Eq. ( 40), it follows that the gluon contribution reads Then, in our model one has ξ g = 0.2.We should note that i) ξ q > ξ q HF S , as it should be, and ii) the massive gluons carry 20% of the HFS momentum fraction and contribute to the total longitudinal fraction by 6% (recalling that P HF S = 0.3).This result indicates that the exchanged gluons in the pion are not soft (differently from the ones considered in Ref. [99] where the subtraction of the effect due to soft gluons is advocated for getting a symmetric PDF from the LF projected BS amplitude).It has to be emphasized that the above analysis, made transparent by the adopted LF variables, is valid in any gauge (both covariant gauges or the light-cone one), and the only difference is the amount of the shift one gets.The possibility to regain the full gauge-invariance by taking into account the additional gluon exchanges that could affect the interaction between the knocked-out quark and the spectator one (see, e.g., the analysis of the gauge-invariance and the hand-bag contribution in Ref. [14,100]) will be explored elsewhere.
The real test of the longitudinal dof is obviously the comparison between the PDF and the experimental data [52].As it is shown in Ref. [76], after evolving u S (ξ) from an assigned initial scale of 360 MeV (suggested by the inflection point of the effective running charge α s (Q 2 )) to the scale of 5.2 GeV, as given by the reanalysis in Ref. [101], the result compares very satisfactorily with the experimental data extracted by taking into account logarithmic resummation effects in the hard part of the Drell-Yan cross-section, as performed in Ref. [98].Moreover, we have achieved a nice agreement with other dynamical calculations, such as the Dyson-Schwinger result of Ref. [95], the basis light-front quantization calculation of Refs.[102,103], and also the recent LQCD outcomes of Ref. [97].In particular, both the overall shape and, importantly, the tail for ξ → 1, gives great confidence in our formalism, and encourages the further steps we have undertaken in this work.Normalized pion transverse distribution function, Eq. ( 43), vs γ/m 2 .The normalization is given by D ⊥ (0) = 22.945 GeV −2 .Thick solid line: full calculation.Dashed line: the same as the full line, but times (γ/m 2 ) 4 .Dash-dotted line: the same as the full line, but times (γ/m 2 ) 2 .Dash-double-dotted line: exponential form e −γ/(m 0.42) 2 , with the parameter from Table 1 of Ref. [35], corresponding to a Gaussian Ansatz for f1(γ, ξ) (see text).Right Panel.Pion unpolarized transverse-momentum distribution f S 1 (γ, ξ), Eq. ( 10), for ξ = 0.5.Solid line: full calculation as in Fig. 4. Dashed line: LF constituent quark model [35,56].Dash-dotted line: LF wave function from DSE calculations [45].Dash-Double-dotted line: NJL model [38].The adopted quark mass m = 255 MeV.
In Fig. 2, one can observe a further comparison, involving the product ξ u(ξ), that sheds more light on the link between the shift of the peak and the gluon dynamics taken explicitly into account in the ladder kernel of the BSE.In particular, we get a scale of 389 MeV for u q (ξ), the solid line in the left panel of Fig. 1, by backward-evolving its first moment, ξ q = 0.471 (cf.Eq. ( 39)), to 0.5, the first moment of u S (ξ), that has an assigned hadronic scale of 360 MeV, as above mentioned and thoroughly discussed in Ref. [76].In the right panel, the comparison at 360 MeV between ξ u S (ξ) and the backward-evolved ξ u q (ξ) shows that the effect of the interaction taken into account in the ladder BSE is reproduced at large extent by applying a leading-order DGLAP evolution with an effective running charge as suggested in Ref. [104] and already applied to our PDF in Ref. [76].This is not surprising once we remind that the dressing of the quarkgluon vertex, as expressed by the effective charge, is governed by the same interaction kernel present in the BSE (i.e. the q q amputated T-matrix).The left panel shows the comparison at 5.2 GeV between the evolved u S (ξ), starting from the scale of 360 MeV, and the evolved u q (ξ), starting from the scale of 389 MeV.Nicely, the difference is even smaller.

C. Transverse degree of freedom
In the left panel of Fig. 3, it is shown the transverse distribution defined by It has to be pointed out that the integration on ξ eliminates the anti-symmetric term f AS 1 (γ, ξ), and therefore one gets the same transverse distribution also by using f q 1 (γ, ξ).In order to emphasize the analysis of the general pattern, we have presented D ⊥ (γ)/D ⊥ (0), so that the widely adopted exponential or power-like fall-off can be readily compared to our result.
In addition, in the left panel of Fig. 3 one can find: i) an exponential form D ⊥ (γ)/D ⊥ (0) = e −γ/(m 0.42) 2 , with the parameter given in Table 1 of Ref. [35], corresponding to the so-called Gaussian Ansätz (recall γ = |k ⊥ | 2 ), amounting to a factorized form for f S 1 (γ, ξ) ∼ u S (ξ)e −γ/(0.42 2 ) very often adopted in phenomenological studies; ii) our full results multiplied by (γ/m 2 ) 2 and iii) our full results multiplied by (γ/m 2 ) 4 .Not surprisingly, a power-like shape provides a better approximation to the dynamically calculated D ⊥ (γ), but the proper power is different from the one expected by the action of only a one-gluon exchange, that should govern the ultraviolet behavior and lead to a (γ/m 2 ) 2 (as suggested by a generalized counting rule in Ref. [105]).Indeed, the adopted form-factor featuring the extension of the quark-gluon interaction vertex (cf.Eq. ( 7)) generates a different power-like fall-off, namely (γ/m 2 ) 4 , as already pointed out in Refs.[88,89].Finally, it is worth noticing that, unlike the PDF, the two terms in f S 1 (γ, ξ) containing derivatives of the delta-function do not contribute, as it is discussed at the end of Appendix C 1.
In the right panel of Fig. 3, it is presented the quantitative comparison between f S 1 (γ, ξ) at ξ = 0.5 and some phenomenological outcomes from i) the approach based on the LF wave function obtained by using the DSE calculation in Ref. [45]; ii) the LF constituent quark-model of Ref. [35,56]; iii) the NJL model with Pauli-Villars regulator as given in Ref. [38].For γ/m 2 → 0, there are remarkable differences that, indeed, are present also on the tails.This last feature impacts the value of γ/m 2 1 2 , as shown in Table I, where, for the sake of completeness, the value of u S (ξ = 0.5) and the pion charge radius are also presented.As can be expected, the larger the average transverse moment, the smaller the radius of charge.The current model has the smaller γ/m 2 1 2 (of the order of the infrared scale Λ QCD , effectively incorporated in the QCD-inspired choice of our parameters) which in turn leads to a larger charge radius, in agreement with the experimental value.TABLE I.The average value γ/m 2 (with m = 0.255 MeV), u S (ξ = 0.5) and the pion charge radius are presented for: i) f S 1 (γ, ξ = 0.5) from the present approach (NIR+BSE); ii) the outcome from the LF wave function obtained by using DSE calculation [45] (LFDSE); iii) the LF constituent quark-model of Ref. [35,56] (LFCQM) and iv) the NJL with Pauli-Villars regulator [38].(Recall that the most recent PDG value of the charge radius is r P DG ch = 0.659 ± 0.004 fm [94]).γ/m 2 1 2 u S (ξ = 0.In Fig. 4, the uTMD f S 1 (γ, ξ) is shown in full, in order to appreciate the main features, i.e. i) the peak at ξ = 0.5 for running γ/m 2 , ii) the vanishing val- (Color online) Pion unpolarized transversemomentum distribution f S 1 (γ, ξ), Eq. ( 10), at the initial scale.The normalization is ues at the end-points and iii) the order of magnitude fall-off already for γ/m 2 > 2. Comparing to other approaches, one can notice the sharp difference with the results from the LF constituent model in Ref. [56] and the LF holographic framework, like in Refs.[40,57] where a double-humped structure is found due to the ξ-dependence in the holographic wave functions.Also the value at ξ = 0.5 and small γ/m 2 is substantially lower than ours (almost an order of magnitude less).Differently, the shape of our f S 1 (γ, ξ) is more similar, i.e. without any double-humped structure, to the one obtained in Ref. [45], where the pion LF-wave function is determined from a beyond rainbow-ladder Dyson-Schwinger equations (DSE) in Euclidean space, by exploiting the γ-dependent moments in ξ and a suitable parametrization of the BS-amplitude.

V. THE SUBLEADING-TWIST UTMDS
In this Section we present the numerical results for (T-even) uTMDs beyond the leading-twist.The detailed expressions can be found in the Appendix E, but it is useful to recall that the decomposition in symmetric and antisymmetric combinations adopted for f 1 (γ, ξ) remains still valid, as well as the relations  46), ii) e S (ξ) (dashed line) and e AS (ξ) (dotted line), Eqs.(45).It is also shown e q EoM (ξ) (dash-dotted line), Eq. ( 47).Right panel.Quark unpolarized collinear PDFs: ξ e q (ξ).Solid line: full calculation as in the left panel.Dashed line: m/M u q (ξ), with u q (ξ) shown in the right panel of Fig. 1.Double-dot-dashed line: ξ e q EoM (ξ), Eq. ( 47).
with the quark and anti-quark contributions.
As introduction to the outcomes of our dynamical approach, it is worth anticipating that the comparison between full calculations and naive estimates one can infer from Eq. ( 17) by using a valence approximation of the leading-twist f 1 (γ, ξ), highlights the inspiring statement one can read in Ref. [77]: the higher-twist distributions are naturally related to multiparton distributions.The role of the exchanged gluons becomes definitely clear through a remarkable shift of the peak in all the sub-leading uTMD we have analyzed, as already discussed in the previous Section, as well as through the sharp difference with the naive estimates, which exclude the effect of the one-gluon exchange.

Longitudinal degree of freedom
In the left panel of Fig. 5, the following collinear PDFs are shown and e q (ξ) = e S (ξ) + e AS (ξ) .
Moreover, in the spirit of Ref. [35], we also present the collinear PDF, e q EoM (ξ), obtained by integrating the first line in Eq. ( 17), but disregarding the gluon contribution, viz where u LF val (ξ)/P val , normalized to 1 (cf.Eq. ( 33)), approximates the integral of f q 1;EoM (γ, ξ).The large difference between our e q (ξ) and (m/M ξ)u LF val (ξ)/P val indicates the sizable role of the gluon contribution from the HFS generated by our dynamical model.In addition, one should point out that the strength of e q (ξ) is spread out on the whole range of ξ, and not concentrated at the end-point ξ = 0 as QCD investigations indicate.The latter feature leads to the singular  44), for ξ = 0.5.Solid line: full calculation.Dashed line: LF constituent quark model of Ref. [35,56], but multiplied by m/(M 0.5) (cf.Eq. ( 47)).Dash-dotted line: the same as the dashed line but with the LF wave function from DSE calculations [45].Dash-Double-dotted line: the same as the dashed line but with the NJL model [38].The adopted quark mass m = 255 MeV.
contribution given in Eq. ( 15) (see, e.g., Ref. [85], for a detailed discussion, but notice that the focus is on the nucleon).
In the right panel of Fig. 5, the comparison between ξ e(ξ) and the other two approximations: i) (m/M )f q 1 (ξ) and ii) (m/M )u LF val (ξ) (cf.Eq. ( 47)) is carried out.The relevance of such a comparison is given by the possibility of more directly assessing the gluon role, since the factor ξ eliminates the singular term present in the QCD analysis of e(ξ), and one remains with the mass contribution (m/M )f q 1 (ξ) and the term from the quark-gluon-antiquark correlator.
Still within the QCD framework (see, e.g., Ref. [85]), the moments ξ n e q , for n ≤ 2, read as follow and, for n > 2, they receive contributions not only from the (n − 1)-th moment of f q 1 (γ, ξ) , but also from the n-th moment of the twist-3 contribution pertaining to the quark-gluon-antiquark correlator.Given TABLE II.The moments ξ n e q of the quark twist-3 e q (γ, ξ) for n < 4, and the ratio R(n, e q , f q 1 ) = ξ n e q / ξ n−1 , and the values of ξ 1 f q 1 = 0.471 and ξ 2 f q 1 = 0.266 have been numerically evaluated).
2.190 0.814 0.445 0.292 R(n, e q , f q 1 ) 2.190 0.814 0.943 1.10 the highly non trivial dynamical content of the e q (ξ) moments, it is interesting to show the results obtained with our dynamical model.In Table II, both the moments up to n = 3 and the ratio R(n, e q , f q 1 ) = ξ n e q / ξ n−1 f q 1 are presented.In particular, as to the first two moments, to get rid of the dependence upon m cur it is helpful to compare the result obtained by multiplying the zero-th and the first moment, (cf.Eq. ( 48)), with final outcome σ π /M .The estimate of σ π at the leading order of the chiral expansion leads to σ π /M = 1/2, as satisfactorily confirmed by the lattice calculations in Ref. [84], where σ lat π = 78.2±4.2MeV, for M = 149.5±1.3MeV and m cur ∼ 4.9 MeV.Eliminating the current quark mass, that is outside our approach, through the above product, we get σ π /M = 1.78, instead of ∼ 0.5.Such a conspicuous difference is surely influenced by the different distribution of the e q (ξ) strength, as already mentioned, and points to a needed enrichment of the gluon dynamics in our approach.However, it is worth mentioning that for a simple non relativistic constituent quark model one has σ N R π = 2m, so that σ N R π /M = 3.64 (with m = 255 MeV the constituent mass), almost twice the result obtained in the BS framework.
In QCD, the ratios R(1, e q , f q 1 ) and R(2, e q , f q 1 ) are equal and amount to m cur /M (see Eq. ( 48)), while R(3, e q , f q 1 ) = m cur /M + ∆ 3 g , where ∆ 3 g contains the contribution from the twist-3 gluonic contribution.In our calculation, the ratios for n = 1, 2 are almost equal, but different from the naive expectation m/M = 1.82 with the adopted m = 255 MeV.The difference with the third ratio indicates the onset of the contribution from the twist-3 gluonic term.

Transverse degree of freedom
The transverse dof can be analyzed globally by introducing the following transverse distribution function, as already accomplished with the leading-twist uTMD, viz.
In the left panel of Fig. 6, it is presented our calculation and the ratio D ⊥ (γ)/D ⊥ (0) to show the similar fall-off, as generated from gluon dynamics and the form-factor featuring the quark-gluon vertex.
A more close view of the decreasing as a function of γ/m 2 is provided by the right panel of Fig. 6, where it is shown the comparison between our calculation of e(γ, ξ = 0.5) and the outcomes obtained by using Eq. ( 47) with i) the LF wave function from the constituent quark model of Ref. [35,56], ii) the LF wave function from DSE calculations [45] and iii) the PDF from the NJL model [38].The differences again point to the role of the interaction in the various approaches, and highlight the relevance of an experimental investigation of the transverse dof.
In Fig. 7, the full dependence of e S (γ, ξ) is presented, displaying a double-hump shape that for larger γ/m 2 becomes smoother and smoother.
In an analogous way, for i = 2 and using Table VIII, one gets the twist-3 f ⊥S(AS) (γ, ξ), with the following decomposition where the above functions are given in Appendix E. Notice that in this case P 2d (γ, ξ; S(AS)) = 0.

Longitudinal degree of freedom
In the left panel of Fig. 8, the following collinear PDFs are shown and the corresponding quark combination.As a reference, it is also presented f ⊥q EoM (ξ), obtained from the The same as in Fig. 5, but for f ⊥q (ξ), f ⊥S (ξ) and f ⊥AS (ξ), Eq. ( 51), and f ⊥q EoM (ξ) as given in Eq. ( 52).Right panel.Quark unpolarized collinear PDFs ξ f q⊥ (ξ).Solid line: full calculation as in left panel.Dashed line: ξ f q⊥ (ξ) obtained by using the second line in Eq. ( 17) and our f q 1 (ξ).Double-dot-dashed line: the same as the dashed line but using the valence approximation of the PDF, u LF val (ξ), with norm equal to 1.  50), for ξ = 0.5.Solid line: full calculation.Dashed line: by using f1(γ, ξ = 0.5) in Fig. 3 from the LF constituent quark model of Ref. [35,56] (cf. the second line in Eq. ( 17), without the gluonic term).Dash-dotted line: the LF wave function from DSE calculations [45].Dash-Double-dotted line: the NJL model [38].The adopted quark mass m = 255 MeV.
second line of Eq. ( 17), without the gluon term, as follows For the sake of completeness, in the right panel of Fig. 8, the product ξ f q⊥ (ξ) is compared to f q 1 (ξ) and u LF val (ξ) that represents the approximation to f ⊥q EoM (ξ) as given in Eq. (52).Also for f q⊥ (ξ), the full calcu-lation substantially differs from approximated evaluations, prompting further investigation of the gluon contributions.

Transverse degree of freedom
Also for f ⊥ (γ, ξ), we introduce the transverse distribution function, viz.
In the left panel of Fig. 9, a comparison, built with the same spirit as in the left panel of Fig. 6, is shown for the ratio P ⊥ (γ)/P ⊥ (0).
A more detailed view of the fall-off can be gained from the right panel of Fig. 9, where f S⊥ (γ, ξ = 0.5) is compared with the results obtained by using i) the LF constituent quark model of Ref. [35,56] (cf. the second line in Eq. ( 17), without the gluonic term).ii) the LF wave function from DSE calculations [45] and iii) the NJL model [38].Finally in Fig. 10, the full dependence of f S⊥ (γ, ξ) is shown.Also in this uTMD, the double-hump shape decreases when γ/m 2 increases.
To summarize, a coherent view of the tail in γ is plainly provided by Figs. 3, 6 and 9. Namely, the interaction taken into account in the ladder kernel together with the extended structure of the quark-gluon vertex governs the fall-off of both the leading and subleading-twist uTMDS.Therefore, the quantitative estimates obtained through our dynamical model, in Minkowski space, is shown to be in a favorable position to provide insights into the interplay between transverse dof and the role of gluons.

VI. CONCLUSIONS
The twist-2 (leading) and twist-3 (subleading) unpolarized (T-even) transverse-momentum dependent parton distribution functions have been calculated for the pion within an approach based on the solution of the Bethe-Salpeter equation in Minkowski space, namely, within a genuinely relativistic quantum-field theory framework.We achieved this goal by exploiting the Nakanishi integral representation of the BSamplitude in order to get actual solution of the homogeneous BSE, in ladder approximation, through a system of integral equations that determine the Nakanishi weight functions relevant for the problem under scrutiny [62,88,89].After obtaining the pion electromagnetic form factor [75], and the pion PDF [76], we extended the yield of our approach by exploring the dependence of the parton distributions upon the transverse momentum.This additional step has itsown importance in view of the planned experimental efforts to achieve a fully three-dimensional investigation of hadrons (mainly of the nucleon and, more challenging, the pion).
The relevant message one gets from our calculations is given by the essential role of the gluon exchange, that cannot be captured by purely phenomenological model.The joint use of the Fock expansion of the pion state, allows us to shed light on the gluonic content of the quark PDF obtained through the BS amplitude, even determining a quantitative estimate, ∼ 6% of the average longitudinal momentum fraction ξ q .Moreover, the latter analysis explains also the source of the small, but theoretically relevant, shift between the u q (ξ) and the PDF that fulfills the charge symmetry (an issue already investigated within the Dyson-Schwinger approach, e.g., in Ref. [99], where a different interpretation was proposed).As to the transverse degree of freedom, a power-like fall-off of the transverse distributions, obtained by integrating on ξ the uTMDs, is supported by the one-gluon exchange interaction that governs the ultraviolet region, according to our calculations.This outcome could suggest to reconsider an exponential or Gaussian Ansatzes for describing the high-momentum content (γ >> m 2 ) of the uTMDS.Summarizing, our approach can be placed among those in which the dynamics can be studied in Minkowski space and in some detail.Moreover, the additive construction of the interaction kernel allows one to address step-by-step recognized effects, achieving an implementation of the dynamics in a controlled way.In this Appendix the traces in Eqs.(10), (18) and (19), are explicitly evaluated, presenting the expressions of the functions a i ,j (k − , γ, z; S(AS)) and b i n; ,j (γ, z; S(AS)).For the sake of convenience, let us rewrite the generic trace entering Eq. ( 23) with O i and η

S(AS) i
given by To proceed one has to insert the expression of the BS-amplitude, Eq. ( 20), and the definitions, Eqs. ( 7) and ( 21), in Eq. (A1).Then one gets the results shown in Tables III, IV and V, for a i j (k − , γ, z; S(AS)).It is also useful to organize the functions a i j in powers of k − for preparing the integration on such a variable (cfr Appendix B), i.e.
In Tables VI, VII and VIII, one can find the expressions for b i n; j (γ, z; S(AS)) where the quantities a i j (k − , γ, z; S(AS)) and b i n; j (γ, z; S(AS)) are given in Appendix A. The first step (see also Refs.[62,75,76]) is to introduce the NIR of φ (k, P ), Eq. ( 22), and then apply the Feynman parametrization as follows φ (k, P )φ j (k, P ) = 30 where Hence, the following general expression, one can straightforwardly deduce from well-known result in Ref. [106] (corresponding to the case m = 0) is useful for performing the relevant integrals where δ (m) (α) = ∂ m δ(α)/∂α m .Finally, combining the results in Eqs.(B2) and (B4), one gets where and the derivatives of the delta function is with respect to α. Recalling that ∂ α/∂z = −1, one can also write Finally, by using where the coefficient c mk can be obtained by repeatedly applying the Leibniz rule for the product of functions and f (m−k) indicates the (m − k)-th derivative (with f (0) (z) ≡ f (z)), one recasts Eq. (B7) in a form more suitable for the further elaboration.In practice, one trades derivatives on the delta functions with derivatives on the functions b i n;ij (γ, z; S(AS)), getting Hence, by taking into account the expressions of b i 2; j (γ, z; S(AS)) and b i 3; j (γ, z; S(AS)), given in Tables VI, VII and VIII, one can drop some derivatives.In particular the second derivative of b i 2; j (γ, z; S(AS)) and all the derivatives of b i 3; j (γ, z; S(AS)), obtaining F i j (γ, z; S(AS)) = −24i F i 0 j (γ, z; S(AS)) + F i 1; j (γ, z; S(AS)) + F i 2; j (γ, z; S(AS)) where ∂ ∂z b i 2; j (γ, z; S(AS)) and Collecting the above results, Eq. ( 23) becomes F i 0; j (γ, z; S(AS)) + F i 1; j (γ, z; S(AS)) + F i 2; j (γ, z; S(AS)) + F i 3; j (γ, z; S(AS)) .(B15) It is also useful for getting more explicit expressions to perform the integral on v in the Eqs.(B11), (B12) and (B13).This can be accomplished by using the following result The combination of the theta-functions implements the constraint 0 ≤ v 0 ≤ 1.Moreover, notice that i) simultaneously changing the signs of z, z and z the function ∆(z, z , z ) does not change sign, this reflects the symmetry with respect ξ = 0.5, as implemented through the charge symmetry in Eq. ( 23); ii) ∆(z, z , z ) is even under the exchange z → z and in the limit z − z = → 0 one has so that the singularity can be addressed without particular problems.By taking into account Eq. (B16), and the symmetries with respect to the transformation z → z and γ → γ , one gets ∂ ∂z b i 2; j (γ, z; S(AS)) and where the functions b i n; j (γ, z; S(AS)) are given in the Tables of Appendix A and Also notice that for a bound state one has and therefore no poles are associated to such a quantity.It should be pointed out that the presence of the theta-functions, that ensure 0 ≤ v 0 ≤ 1, prevents singular behaviors, shrinking the area of integration in the space z ⊗ z , when z → z .Interestingly, in the Appendix C 1, it is shown that only F i0 j (γ, z), i.e. without derivative of the delta-function, contributes to the norm of the twist-2 uTMD f 1 (γ, ξ).
Appendix C: The leading-twist uTMD f In this Appendix, the symmetric and anti-symmetric combinations of the quark and antiquark leading-twist uTMDs are explicitly given and their relevant features discussed.
By specializing the expressions in Eq. ( 27), one can write where the four contributions are obtained from Eqs. (B19), (B20), (B21) and (B22), respectively.Inserting the functions b 0 n; j (γ, z; S) listed in in the first three columns of Table VI in Appendix A, one gets the following non vanishing symmetric contributions, viz. and with β 0 (z 2 ), Ḡ j and v 0 given in Eq. (B23).The symmetry property under the transformation z → −z can be easily demonstrated, recalling also that under the exchange z → −z and γ → γ the functions Ḡ j (γ , z ; γ , z ) do not change, since the NWFs g i (γ, z; κ 2 ) are even for i = 1, 2, 4 and odd for i = 3.Moreover, under z → −z and z → −z one also has v 0 → (1 − v 0 ), so that β 0 (z 2 ) remains unchanged, as well as Θ(z − z) Θ(z − z )/(z − z ).The anti-symmetric combinations are The anti-symmetry with respect to the transformation z → −z can be easily shown by using the properties listed below Eq. (C4).