Two-current correlations in the pion in the Nambu and Jona-Lasinio model

We present an analysis of two-current correlations for the pion in the Nambu--Jona-Lasinio model, with Pauli--Villars regularization. We provide explicit expressions in momentum space for two-current correlations corresponding to the zeroth component of the vector Dirac bilinear in the quark vertices, which has been evaluated on the lattice. The numerical results show a remarkable qualitative agreement with recent lattice data. The factorization approximation into one-body currents is discussed.


Introduction
Due to the high partonic densities reached, processes with more than two partons from the two colliding protons participating in the actual scattering -the so called multiple parton interactions-are likely to happen at the LHC [1]. Double parton scattering (DPS), the simplest form of multiple parton interactions, involves two simultaneous hard collisions. It has been indeed observed at the LHC (see, e.g., Ref. [2]). The DPS cross section is expressed in terms of double parton distribution functions (dPDFs) [3,4]. The latter are related to the number density of two partons located at a given transverse separation in coordinate space with given longitudinal momentum fractions. The information encoded in these distributions therefore complements that accessed through electromagnetic probes in terms of generalized parton distributions (GPDs) [5,6], one of the flagships of the future Electron Ion Collider (EIC) [7]. If measured, dPDFs would therefore represent a novel tool to study the threedimensional hadron structure [8][9][10][11][12]. Indeed, they are sensitive to two-parton correlations not accessible via one body distributions, e.g. Parton Distribution Functions (PDFs) and GPDs -see Ref. [1,13] for a recent report. Since they are non perturbative objects, dPDFs have been evaluated within models of the hadron structure -estimates exist at low momentum scales for the proton dPDFs [14][15][16][17][18][19][20]. In order to match theoretical predictions with future experimental analyses, the results of these calculations are then evolved using the proper QCD evolution equations to reach the high momentum scale of the data. Developments on the evolution properties of dPDFs can be found in, e.g., [21].

Definitions and isospin decomposition
Experimental measurements of pion dPDFs appear rather challenging at the moment. Nevertheless, from a complementary viewpoint, lattice evaluations of matrix elements that can be connected to the Mellin moments of double parton distributions are possible [4]. A relevant step towards this goal has been done in Ref. [31], where lattice data of correlation functions of two local currents in a pion, have been presented. In the expression above, the pion charges -k, k = +, −, 0-may be different in the bra and ket state, while the four-momentum P is the same in both, and the two currents are applied at the same time but at different spatial positions, i.e., one has y µ = (0, y). It differs from the dPDFs which are defined from a similar product of two one body operators taken at a distance y µ = y + = 0, y − , y ⊥ . Hence, the currents of interest here are where q and q are u or d quark fields and Γ n are the Dirac bilinears 1, iγ 5 , γ µ , γ µ γ 5 for n = S, P, V µ , A µ , respectively. In the present work we will only consider the time component of the vector current, i.e.
and we will compare our model results with Lattice data in the case of two equal currents of this kind. As shown in Ref. [31], the matrix elements (1) are not all independent due to constraints from isospin symmetry and discrete symmetries. We follow the treatment of Ref. [31] defining the isosinglet current: as well as the isotriplet currents O a n =Qτ a Γ n Q , a = 1, 2, 3 where τ a are the Pauli matrices and Q = (u, d) is the isodoublet of quark fields. As a consequence, the following isospin decomposition is obtained Here the real functions F and G -quantities evaluated in what follows-depend on the pion momentum P and on the current indices n, n [31].
3 Two-current correlations in the NJL model.

Leading order
We describe now the calculation of the two-current correlations, and in particular of the functions F i , with i = 0, 1, 2, and G l , with l = 1, 2, in the NJL framework. This framework is the most realistic model for the pseudoscalar mesons based on a local quantum field theory built only with quarks [32]. It respects the realization of chiral symmetry and gives a good description of low energy properties. Mesons are described as bound states, in a fully covariant fashion, using the Bethe-Salpeter amplitude, in a field theoretical framework. In this way, Lorentz covariance is preserved. The NJL model is a non-renormalizable field theory and therefore a regularization procedure has to be implemented. We have performed our calculations in the Pauli-Villars regularization scheme -a well established method. The NJL model, together with its regularization procedure, can be regarded as an effective theory of QCD. Some basic features of the NJL model and details on the regularization scheme are reported in Appendix A.1. Model calculations of meson partonic structure within this approach have a long story of successful predictions [33][34][35][36][37][38][39][40][41][42][43]. Collinear parton distributions evaluated within a low energy model can be associated to a RGE scale denoted by Q 0 . The determination of that hadronic scale reflects the degrees of freedom of the model -here the NJL model contains only valence quarks [44,45]. The identification of such a low lying scale enables predictions from the model to measured quantities, through perturbative QCD evolution. Let us describe the main steps of our calculation. We define, following Eq. (6), with α, β = s, 1, 2, 3. In order to proceed to the calculation in the NJL model, we perform the Fourier transform to momentum space We use in this expression the following definition of the state for the pion i with φ π i the quark-pion vertex function for a π i . In the NJL model the amplitude φ π i (k, P ) is independent on the relative and total quark-antiquark momenta, k and P, respectively, and we have φ π i (k, P ) = ig πqq iγ 5 τ i , where g πqq is the quark-pion coupling constant and τ i is the isospin matrix associated to the corresponding pion π i . Inserting Eq. (9) in the definition (8), we obtain contributions associated to different Feynman diagrams. The fact that G ij αβ ( y) is defined at y 0 = 0 implies that where T is the Feynman amplitude corresponding to the diagram. The two contributions to be considered at leading order are shown in Fig. 1. In the following, we will call them rhombus and trapezoid diagrams, respectively. Therefore one has where G (r)ij αβ ( q) and G (t)ij αβ ( q) stand for contributions from the rhombus and trapezoid diagrams, respectively, which read Here τ α = (τ s = 1, τ ) , with τ the standard Pauli matrices. According to their isospin structure, which is determined by φ i and τ α terms, the two expressions above yield four types of contributions once the isospin traces are evaluated -a straightforward lengthy procedure. For the rhombus contribution, they read with In the equations above, m is the quark mass and we have assumed that Γ n = Γ n = Γ. For the trapezoid contribution, they read and By comparing Eq. (15) and (17) with Eq. (6) we obtain the relations, with In Ref. [31] it is shown that iF 3 ( q) is different from zero only if Γ n = Γ n , and charge conjugation constraints imply that G ( q) = 0 for the cases Γ n = Γ n . We observe that these results are automatically reproduced if I r Γ (m, q) and I t i Γ (m, q) are functions of | q| only. We will see later on that this fact actually occurs in our calculation.
The integrals over q 0 and k 0 in the expressions I r,t i Γ ( q) can be easily evaluated using the poles of the propagators. For instance, when Γ = V 0 = γ 0 we obtain where Now, the angular integration can be performed, arriving to where k = | k| and q = | q| . In a similar way we obtain After regularization, the same integrals becomē One can see that these integrals do not depend on the direction of q, as anticipated, and therefore Eqs. (21) become We conclude this section observing that the above equations define the leading order contributions to the functions (20) in the Γ = γ 0 = V 0 case:

Dressing the vertex and other contributions
In addition to the contribution we have just discussed, one has to consider that, in the present NJL scheme, the local vertex can be dressed by the chiral interaction through the diagrams depicted in Fig. 2. This dressing, which could produce corrections to the results obtained in the previous subsection, implies to change the bare vertex τ α Γ as follows, for the isoscalar and isovector vertices respectively where q µ = p µ 1 − p µ 2 = (0, q) . In the vector current case, only the isoscalar current can be dressed in principle. However, it is also easily proved that such a contribution vanishes, as  Figure 3: Other type of contribution originated by the dressing of the currents.
Therefore, in the particular case of the vector current under investigation here, there is no dressing of the vertex. Finally, we consider the diagram depicted in Fig. 3. This contribution can be written in the form -(p) stands for pole contribution, with and The direct calculation gives, for M ij , with I 2 (q 2 ) and I 3 (q 2 ) defined in App. A.1. The two quark intermediate state is described by the quantity, and for N αβ we have with Putting all these results together, we obtain By comparing to the structures for the isospin amplitudes defined in Eqs. (6), it is observed that the pole contributions add to the bare leading-order results in Eqs. (20) in the following manner with In the case Γ = γ 0 , shown in the following as an example, I p Γ (m, q) is a function of q = | q|,

and, after regularization, the last integral becomes
Putting all these results together we have finally and the final result for the functions under investigation, including all the contributions evaluated so far, can be summarized as follows

Comparison with lattice data
The comparison of our results with the lattice evaluation of Ref. [31] in the vector-vector case is presented in Figs [31], respectively. A qualitative agreement is clearly found. A similar conclusion holds for F 2 , given in Eqs. (43) and shown in Fig. 5.
Besides, the results of the evaluations of the rhombus and trapezoid diagrams in our approach are related, in the calculation scheme of Ref. [31], to lattice contractions called C 1 and C 2 . The correspondence formally reads     Fig. 6, but in momentum space. The latter quantities can also be directly compared to the findings of Ref. [31], which will come of useful for the interpretation of contributions from diagrams. Our results for these functions are shown in Figs. 6 and 7, in coordinate and momentum space, respectively. In particular the former results show an encouraging agreement with the lattice data reported, in coordinate space, in Figs. 11 (a) of Ref. [31]. Overall, the NJL scheme seems to tentatively reproduce the trend of the available lattice results.
Following an argument discussed in Ref. [29], it is observed that, in the absence of correlations, two-current distributions could be computed from the single-current ones. This factorization hypothesis is comparable to an argument that has often been used more recently for dPDFs -see, e.g., Ref. [4]. This hypothesis is examined in Ref. [31]. In particular, for the isospin combination π + (P )| O uu n (y) O dd n (0) |π + (P ) , inserting a complete set of intermediate states between the two currents and keeping, among them, only the pion ground state, it can be shown that the two-current observable can be expressed in terms of single-current Form Factors.
In the most recent evaluation, the the full lattice result for the C 1 contraction -as it dominates the relevant isospin combination-is compared to single-current expressions using known results for the charges -i.e. Form Factors at zero momentum transfer-as well as a common power law fall-off for the momentum transfer dependence of the Form Factor. While the results at zero momentum transfer are acceptable, the attempt is found to generally fail to a large extent.
Based on our previous experience with dPDFs [24], we have evaluated diagrams of the type shown in Fig. 8, where the intermediate two quark state can have the quantum numbers of a σ particle or a pion, as allowed in the present calculation scheme. After a tedious evaluation detailed in the Appendix A.3 for Γ = γ 0 and assuming dominance of the pion pole, in agreement with the procedure followed in Ref. [31], it is found that the only contribution arising from these diagrams can be written where F em is the pion electromagnetic form factor and E π = m 2 π + q 2 . This expression -after using the relevant isospin combination from Eq. (66)-agrees with that adopted in Ref. [31] to approximate the full result given by the function C 1 . In the NJL scheme used here the contributions to C 1 arise from the rhombus diagram only; such diagram cannot accommodate for any contribution of the type of Fig. 8. Therefore no comparison between the result of Eq. (45) with the contributions coming from C 1 -as suggested in Ref. [31]-can be made according to the spirit of the present paper. Accordingly we do not discuss relations based on that factorization hypothesis any further.

Conclusions
Two-partons correlations represent an elusive non perturbative information, theoretically encoded in two-current correlations as well as double parton distribution functions. The latter quantities are under theoretical investigation for the proton. There are perspectives to observe them during the LHC operation in the next years. While a possible access to the same observable for the pion appears challenging, lattice results have already been obtained for two-current correlations in the pion, a quantity that could be related to double parton distribution functions [31]. A consistent field-theoretical approach, based on the Nambu-Jona-Lasinio model with Pauli-Villars regularization, has been used in this paper for a systematic analysis of two-current correlations in the pion. We have given explicit expressions for the time-time component of the vector-vector two-current correlations in momentum space. We have compared our results for two-current correlations with lattice data, obtaining an overall good qualitative agreement. As happens in model calculation of the related quantities double parton distributions, it is found that, at the low scale relative to the model description, these functions encode novel non-perturbative information, not present in one-body quantities. This conclusion is highlighted in the consideration of a possible factorization of the two-current correlation into one-body Form Factors, which is found to be an unconnected contribution in our low-energy model. It will be interesting to compare our results for pion double distribution functions with forth-coming lattice data, whose evaluation is presently in progress [46].
Effectively, in the large N c approximation, the quark condensate is defined by Pion and sigma masses are defined by the relations with the scalar polarization and the pseudoscalar polarization The pion-quark and sigma-quark coupling constant are respectively defined by The pion decay constant is The NJL model is a non-renormalizable field theory and a regularization procedure has to be defined for the calculation of I 1 (m) and I 2 (m, q 2 ) . We will introduce now the Pauli-Villars regularization method for the NJL model.

A.2 Pauli-Villars regularization scheme
In Section 2, we have used the Pauli-Villars regularization in order to render the occurring integrals finite. The way to proceed in this method is: (1) remove from the numerator all the powers of the integrated momentum, which will be replaced by external momenta, and the mass of the constituent quark, m; (2) for each resulting integral, which is of the form make the substitution with M 2 j = m 2 + j Λ 2 , c 0 = c 2 = 1 and c 1 = −2. Following this procedure, we obtain for the momentum integral of one propagator and for the one of two propagators With the conventional values ūu = −(0.250 GeV) 3 , f π = 0.0924 GeV and m π = 0.140 GeV, we get m = 0.238 GeV, Λ=0.860 GeV and m 0 = 5.4 MeV. For the pion-quark coupling constant we get g 2 πqq = 6.279. We can obtain the chiral limit taking m 0 = 0, without changing Λ and m. In that case ūu and f π do not change but one has m π = 0 and g 2 πqq = 6.625. For a proper comparison with the lattice data of Ref. [31], we have applied the same model with Λ = 1.022 GeV, m = 0.242 GeV and m 0 = 20 MeV, which lead to a massive pion with m π = 0.300 GeV, f π = 0.1 GeV and ūu = −(0.285 GeV) 3 , as used in Ref. [31]. In this case we have g 2 πqq = 4.527.

A.3 Intermediate pion states
Let us consider now the evaluation of the calculation of the process depicted in Fig. 8. There are eight different diagrams of this type. This contribution can be written in the following form The first (second) line of this equation describes an intermediate two quark state with the quantum numbers of a σ(π) particle. In the case of the intermediate pion the additional isospin index, c, corresponds to its isospin. The quantities introduced in Eq. (60) are Performing the isospin traces, we obtain the following structures (we disregard the G (F )ij as and G (F )ij sa cases): Using Eqs. (51), (52), (53) and (54), the intermediate two quark states appearing in Eq. (66) can be approximated by the quantities, Let us consider, as always in this paper, the Γ = γ 0 case. Performing the traces present in Eq. (67) we have so that X F S V 0 ( q) = X F + P S V 0 ( q) = Z F S Γ ( q) = 0. After a tedious but straightforward calculation we obtain, where F em (q 2 ) and F em (q 2 ) are F em q 2 = 2N c g 2 πqq P 2 P 2 − (P · P ) 2 P 2 P · P − P 2 I 2 P 2 + P 2 P · P − P 2 I 2 P 2 + (P · P ) 2 P · P − P 2 − P 2 I 2 (P − P ) 2 − P 2 P 2 2 P · P − P 2 − P 2 I 3 (P, P ) , F em q 2 = 2N c g 2 πqq P 2 P 2 − (P · P ) 2 P 2 P · P − P 2 I 2 P 2 − P 2 P · P − P 2 I 2 P 2 + (P · P ) P 2 − P 2 I 2 (P − P ) 2 − P 2 P 2 P 2 − P 2 I 3 (P, P ) , with P µ = P µ + q µ . These two form factors correspond to the standard electromagnetic form factors, π (P )| J µ em |π (P ) = (P + P ) µ F em q 2 + (P − P ) µ F em q 2 .
When both pions are on-shell, F em (q 2 ) vanishes and F em (q 2 ) reads F em q 2 = 4N c g 2 πqq m 2 π + (P · P ) −m 2 π I 2 m 2 π − (P · P ) I 2 (P − P ) 2 + m 4 π I 3 (P, P ) Putting together all these ingredients we have The q 0 integral will be dominated by the pole contribution and we obtain the final result which is the expression reported at the end of Section 4.