Nucleon’s energy-momentum tensor form factors in light-cone QCD

We use the energy-momentum tensor (EMT) current to compute the EMT form factors of the nucleon in the framework of the light cone QCD sum rule formalism. In the calculations, we employ the most general form of the nucleon’s interpolating field and use the distribution amplitudes (DAs) of the nucleon with two sets of the numerical values of the main input parameters entering the expressions of the DAs. The directly obtained results from the sum rules for the form factors are reliable at Q ≥ 1 GeV : To extrapolate the results to include the zero momentum transfer squared with the aim of estimation of the related static physical quantities, we use some fit functions for the form factors. The numerical computations show that the energy-momentum tensor form factors of the nucleon can be well fitted to the multipole fit form. We compare the results obtained for the form factors at Q = 0 with the existing theoretical predictions as well as experimental data on the gravitational form factor d 1 (0). For the form factors M 2 (0) and J(0) a consistency among the theoretical predictions is seen within the errors: Our results are nicely consistent with the Lattice QCD and chiral perturbation theory predictions. However, there are large discrepancies among the theoretical predictions on d 1 (0). Nevertheless, our prediction is in accord with the JLab data as well as with the results of the Lattice QCD, chiral perturbation theory and KM15-fit. Our fit functions well define most of the JLab data in the interval Q ∈ [0, 0.4] GeV , while the Lattice results suffer from large uncertainties in this region. As a by-product, some mechanical properties of the nucleon like the pressure and energy density at the center of nucleon as well as its mechanical radius are also calculated and their results are compared with other existing theoretical predictions.

Understanding of the inner structures of the nucleons based on the quarks and gluons degrees of freedom is one of the most prominent research subjects of the hadron physics. A powerful instrument to probe the nucleon's structure is to investigate the various form factors (FFs) of the nucleon as they carry direct information on the nucleon' substructure and geometric shape. Indeed, the electromagnetic FFs of the nucleon unveil the way of distributions of the charge and magnetization of quarks inside the nucleon. The scalar and axial-vector FFs encompasses information on specific viewpoints of the nucleon's substructure nature such as chiral and flavor symmetries and their breakdown. For these reasons, extensive research has been carried out on different form factors of the nucleon for decades. However, the gravitational or energy-momentum tensor form factors (EMTFFs) of the nucleon have recently been received considerable attention both in theory and experiment, despite they were recommended by Pagels in 1966 [1]. These form factors cannot be extracted directly from the experiment: They can be obtained from hard exclusive reactions using the Mellin moments of certain generalized parton distributions (GPDs). These form factors give us a tool for systematic studies of the properties of the nucleon and calculate different related observables like energy, angular momentum and pressure distributions inside the nucleon, etc.
The matrix element of the EMT current between the nucleon states is characterized by four form factors and parametrized as follows [2,3] whereP = (p ′ + p)/2, ∆ = p ′ − p, σ µν = i 2 [γ µ , γ ν ], Q 2 = −∆ 2 and u(p, s) is the spinor of the nucleon with mass m N . Here, M 2 (Q 2 ), J(Q 2 ), d 1 (Q 2 ) andc(Q 2 ) are the EMTFFs of the nucleon. The M 2 (Q 2 ) form factor gives knowledge on the fractions of the momenta carried by the quark and gluon constituents of the nucleon. It is also related to the energy density distribution inside the nucleon. The J(Q 2 ) form factor gives instruction about how the total angular momenta of quarks and gluons form the nucleon's spin. The third form factor, d 1 (Q 2 ) (called the D-term), provides information on the distribution and stabilization of strong force in the nucleon. It can be obtained by the beam charge asymmetry in deeply virtual Compton scattering (DVCS). The negative sign of this form factor at zeromomentum transfer obtained from various theoretical studies is thought to be in connection with the spontaneous chiral symmetry breaking [4][5][6][7]. The quark and gluon parts of the EMT current are not conserved separately, but their sum is conserved. The form factorc(Q 2 ) characterizes the order of the non-conservation of the quark part of EMT current. This form factor is substantial to specify the distributions of the pressure forces inside the nucleon separately for quarks and gluons. This is also used to study the forces among quarks and gluons inside the nucleon. Hence, the EMTFFs provide new perspectives to the internal structure of the nucleon. More details can be found in a recent paper [8].
The EMTFFs of the nucleon have been investigated in the framework of Lattice QCD [9][10][11][12][13][14][15][16], chiral perturbation theory (χPT) [17][18][19][20][21][22], instant and front form (IFF) [23], Skyrme model [24,25], chiral quark soliton model (χQSM) [7,[26][27][28][29][30][31][32][33][34][35], light-cone QCD sum rules at leading order(LCSR-LO) [36], dispersion relation (DR) [37] and instanton picture (IP) [38]. In Ref. [9], Hagler et al. calculated the quark part of the EMTFFs of the nucleon by means of the Lattice QCD. They obtained J q ∼ 0.34 ± 0.04 and M q 2 ∼ 0.68 ± 0.07 at re-normalization scale of µ 2 = 4 GeV 2 . In Ref. [10], Mathur et al. calculated the quark total angular momentum of the nucleon from the quark EMTFFs on the Lattice QCD and they found J q = 0.30 ± 0.07 at re-normalization scale of µ = 1.74 GeV. In Ref. [11], Gockeler et al. performed a quenched Lattice computation of the first moment of twist-two GPDs of the proton, and assessed the total quark contribution to the spin of the proton. They obtained J q = 0.33 ± 0.07, M q 2 = 0.55 ± 0.11 and d q 1 = −1.0 ± 0.05 at re-normalization scale of µ = 2 GeV. In Ref. [12], Bratt et al. presented their predictions for the substructure of the nucleon from a mixed-action computation using 2+1 flavors of asqtad sea and domain wall valence fermions. They carried out extrapolations of their data based on various chiral effective field theory pattern at re-normalization scale of µ 2 = 4 GeV 2 . In Ref. [13] Hagler et al. presented a exhaustive study of the lowest moments of nucleon GPDs in 2 + 1 Lattice QCD by the help of domain wall valence quarks and refined staggered sea quarks without including the disconnected diagrams. They performed extrapolations of their results based on different chiral effective field theory schemes at re-normalization scale of µ 2 = 4 GeV 2 . In Ref. [14], Brommel et al. reported on a calculation of nucleon's GPDs based on simulations with two dynamical non-perturbatively improved Wilson quarks with pion masses down to 350 MeV. They found J q = 0.226 ± 0.013 and M q 2 = 0.572 ± 0.012 at re-normalization scale of µ 2 = 4 GeV 2 . In Ref. [16], Deka et al. reported a comprehensive computation of the quark and gluon momenta in the nucleon. The computations encompass the contributions of quarks coming from both the connected and disconnected inclusions at re-normalization scale of µ = 2 GeV. In Ref. [22] Dorati et al. evaluated the basic properties related to the structures of baryons at low energies, by the help of the method of the covariant chiral perturbation theory in the baryon sector at leading-one-loop order. They investigated the quark-mass dependence of the isoscalar moments in the forward limit and estimated the contributions of quarks to the total spin of the nucleon at re-normalization scale of µ = 2 GeV. In Ref. [23] Lorce et al. evaluated, in details, the distributions of energy, radial pressure and tangential pressure inside the nucleon in the framework of both the instant form and the front form of dynamics at re-normalization scale of µ = 2 GeV. In Ref. [36] the author developed a method based on the light-cone sum rules at the leading order of α s to compute the gravitational form factors for the valence quark combinations in a nucleon at re-normalization scale of µ 2 = 2 GeV 2 . In Refs. [24,25], the EMTFFs of the nucleon are studied via the Skyrme and in-medium modified Skyrme models and they discuss how medium effects act on the form factors. In Refs. [7,[26][27][28][29][30][31][32][33][34][35], the EMTFFs of the nucleon are studied by means of the chiral quark-soliton model. It should be noted here that the Skyrme and χQSM models show the total form factors which are re-normalization scale independent. In Ref. [37], Pasquini et al. presented a depiction of the D-term form factor for hard exclusive reactions, making use of unsubtracted t-channel dispersion relations at re-normalization scale of µ 2 = 4 GeV 2 . In Ref. [39], Burkert et al. presented an analysis of JLab data where an experimental information on the quark contribution to the D-term was obtained at re-normalization scale of µ 2 = 1.5 GeV 2 . In Ref. [40] the distributions of pressure and shear forces inside the proton are discussed via Lattice QCD computations of the EMTFFs at re-normalization scale of µ = 2 GeV.
In the present study, we compute the quark parts of the EMTFFs of the nucleon by the help of the light-cone QCD sum rules (LCSR) method, as one of the powerful and successful nonperturbative methods in hadron physics [41][42][43]. The LCSR method is based on the operator product expansion (OPE) near the light-cone and expansion is carried out over the twists of the operators and the features of the hadrons under study are stated with respect to the features of the vacuum and the light-cone distribution amplitudes of the hadrons. Since the form factors are quantities with respect to the features of the vacuum and distribution amplitudes of the hadrons, any uncertainties in these parameters are reflected to the uncertainties of the estimations of the form factors. This method is quite accomplished in determining the baryonic form factors at high Q 2 (see e.g. [44][45][46][47][48][49][50][51][52][53][54]).
This manuscript is organized as follows. In section II, we formulate and derive the light-cone QCD sum rules for the nucleon EMTFFs. In section III, we present our numerical results for the nucleon EMTFFs. In section IV, we discuss the mechanical structure of the nucleon using the EMTFFs. Section V is reserved for the conclusions on the obtained results. A remark on different definitions and notations for the EMTFFs is given in Appendix A. The explicit expressions of the EMTFFs are moved to the Appendix B.

II. FORMALISM
In order to calculate the EMTFFs of the nucleon within LCSR, we begin our calculations with the subsequent correlation function: where q = p ′ − p and T µν is the energy-momentum tensor current. The quark and gluon parts of the EMT current are defined as The second part of Eq. (3) can be rewritten as [8] where γ m is the anomalous dimension of the mass operator. It should be noted here that we work in the chiral limit (m u = m d = 0) and ignore the gluon fields contributions. Hence, in the present study, we will deal only with the the first part of the quark part of the EMT current in Eq. (3). The covariant derivative ← → D µ is defined as and A µ is the gluon field. In the correlation function above, J N (0) is the nucleon's interpolating current. In this study, we decide on the most general form of the interpolating current for nucleon, which is written as where a, b, c are the color indices, t is an arbitrary mixing parameter, and C is the charge conjugation operator. Choosing t = −1 reduces the above current to the famous Ioffe current. In order to calculate sum rules for EMTFFs, we need to evaluate the correlator in two different languages. First, it is computed with respect to the QCD degrees of freedom: In terms of the parameters of the quarks and gluons and their non-perturbative interactions with the QCD vacuum.This representation is called the QCD or theoretical representation of the correlation function and it is obtained by the help of OPE in deep Euclidean space. In the second representation, the correlation function is calculated in terms of the hadronic parameters like the mass, residue, form factors and other hadronic degrees of freedom. This representation of the correlation function is called the physical or hadronic representation. Equating the coefficients of various Lorentz structures from two different representations of the same correlation function and carrying out a Borel transformation with the aim of eliminating the contributions of the continuum and higher states, we obtain sum rules for the EMTFFs of the nucleon. To further suppress the unwanted contributions and enhance the ground state contribution we apply the continuum subtraction procedure with accompany of the quark hadron duality assumption.
First we focus on the calculation of the hadronic side of the correlation function. To this end, we saturate the correlation function with a complete set of the nucleon state, the integration over four-x leads to where dots represent the unwanted contributions coming from the continuum and higher states. The above relation is further simplified by introducing the following definition: where λ N is the nucleon overlap amplitude or its residue. Inserting Eqs. (1) and (8) into Eq. (7), and performing summation over the spins of the Dirac spinors, we obtain the hadronic side of the correlation function in terms of the hadronic properties as well as different Lorentz structures as On QCD side, we insert the explicit forms of the interpolating current for the nucleon and the EMT current into the correlation function and perform the required contractions using the Wick theorem. The resultant expression is in terms of the light quark propagator as well as the matrix elements of the quark fields sundwiched between the vacuum and nucleon states. The latter will be defined in terms of the nucleons DAs later. As a result we get where I is the unit matrix and S(x) represents the up/down quark propagator which is given, in the limit m q = 0, as where, qq is the quark condensate and m 2 0 is specified with respect to the mixed quark-gluon condensate as m 2 0 ≡ qg s G µν σ µν q / qq . Since the expressions proportional to the gluon field strength tensor (G µν ) are related to the four and five-particle distribution amplitudes, the contributions of these terms are expected to be small [55] and, therefore, these contributions will be neglected in our calculations. Furthermore, the terms proportional to qq are killed and they do not contribute after applying the Borel transformations. Hence, only the first term of the propagator survives in the calculations.
As it is clear from Eq. (10), to proceed in the calculations, we need to know the matrix elements of the quark operators sandwiched between the vacuum and nucleon states, i. e.
where a 1 , a 2 and a 3 are some real numbers. These matrix elements are parameterized in terms of the nucleon's distributions amplitudes of different twists in the basis of the QCD conformal partial wave expansion approach [47]: where N φ is the spinor of the nucleon. The "calligraphic" functions, leaving aside the terms proportional to x 2 which contain V M 1 , A M 1 and T M 1 , can be denoted in terms of the functions of the specific twist as: where A i , P i , S i , T i and V i are axialvector, pesudoscalar, scalar, tensor, and vector distribution amplitudes, respectively. The expansion of the matrix element is essentially an expansion in increasing twists of the distribution amplitudes. The distribution amplitudes V 1 , A 1 and T 1 have twist three, are functions of a i px, can be described as where x i with i = 1, 2, 3 are equivalent to the longitudinal momentum fractions carried by the quarks inside the nucleon. The explicit representations of the nucleon distribution amplitudes (A i , P i , S i , T i , V i ) can be found in Ref. [47]. The desired LCSRs for the EMTFFs are acquired by matching the coefficients of various structures from both the hadronic and QCD sides of the correlation function in the momentum space. We use the structures p ′ µ q ν , p ′ µ p ′ ν q /, q µ q ν and g µν to find the sum rules for the form factors M q 2 (Q 2 ), J q (Q 2 ), d q 1 (Q 2 ) andc q (Q 2 ), respectively. For the EMTFFs of the nucleon we obtain: andc The explicit forms of the various F functions that come into view in Eqs. (14) to (17) are presented in Appendix B with respect to the distribution amplitudes of the nucleon. For the sake of brevity, in Appendix B, only the results for the M q 2 (Q 2 ) form factor are presented, explicitly. The last step is to apply the Borel transformation with respect to the variable p ′2 as well as the continuum subtraction with the aim of suppression of the contributions of the higher states and continuum. These steps are performed by the help of the subsequent replacement rules (see e.g. [47]): where, M 2 is the Borel mass squared parameter and x 0 is the solution of the quadratic equation for s = s 0 : with s 0 being the continuum threshold.
One of the main input parameters in the expressions of the sum rules for EMTFFs is the nucleon's residue, λ N . We use the expression of this parameter, in terms of the hadronic and QCD degrees of freedom as well as the auxiliary parameters entering the calculations, calculated via mass two-point sum rules [50]. It is given as where x = s 0 /M 2 , and

III. NUMERICAL RESULTS
The present section encompasses the numerical analyses of nucleon EMTFFs. In order to obtain the numerical results of the form factors, expressions of the distribution amplitudes for nucleon are needed. We borrow them from Ref. [47]. These distribution amplitudes include eight nonperturbative hadronic parameters, which are obtained in the framework of different models. In further numerical computations we take into account two different sets of these parameters: 1) QCD sum rules (QCDSR) based distribution amplitudes, where corrections to the distribution amplitudes are considered and the parameters in distribution amplitudes are obtained from QCDSR (Set-I), 2) The condition that the next to leading conformal spin contributions vanish, fixes five of the eight parameters, and remaining three parameters (f N , λ 1 , λ 2 ) are borrowed from QCDSR. This set is called asymptotic set (Set-II). The values of these parameters in two different sets are given in Table I. In addition we use: m u = m d = 0, m N = 0.94 GeV, qq = (−0.24 ± 0.01) 3 GeV 3 and m 2 0 = 0.8 ± 0.1 GeV 2 [56]. There are three auxiliary parameters of the QCDSR to be fixed: the continuum threshold s 0 , the Borel mass parameter M 2 and the mixing parameter t. The continuum threshold starts from the point, where the excited states and continuum contribute to the correlation function. We use the continuum threshold in the range s 0 ≃ (2.25 − 2.40) GeV 2 , which is pretty much fixed in the literature from the nucleon spectrum analyses. The working region of M 2 is decided to be in the interval 1.0 GeV 2 ≤ M 2 ≤ 2.0 GeV 2 . As it can be seen from Fig. 1 (as an example), the results of the form factors are roughly independent of the Borel parameter in the interval 1.0 GeV 2 ≤ M 2 ≤ 2.0 GeV 2 . We include into the final results the errors coming from the variations of the physical observables with respect to the auxiliary parameters, which remain below the limits accepted by the sum rules computations. The next step is to specify the optimal mixing parameter t. Our numerical calculations indicate that the form factors are not sensitive   to cosθ (with t = tanθ) when it varies in the region -0.2 ≤ cosθ ≤ -0.4. We see that the famous Ioffe current for the nucleon, which corresponds to the choice cosθ ≃ −0.71 remains out of the reliable the working region. Figure 2 shows the changes of the energy-momentum tensor form factors of the nucleon with respect to the momentum transfer squared Q 2 . As is also clear from this figure, the sum rules for the EMTFFs give reliable results for Q 2 ≥ 1. To extrapolate the results to the smaller points as well as Q 2 = 0, which enables us to compute the static properties of the nucleon, we use some fit functions such that the results of fit functions coincide with the LCSRs predictions at the region Q 2 ≥ 1. Our numerical computations show that, the EMTFFs of the nucleon can be well described by the multipole fit functions defined as [36] The values of the fit parameters, i. e. the form factors at Q 2 = 0, m p and p for different EMTFFs obtained from the sum rules analyses are shown in Table II. The errors in the presented results are due to the variations in the computations of the working regions of M 2 , s 0 and mixing parameter t as well as the uncertainties in the values of the input parameters and the nucleon distribution amplitudes. Although the central values of the form factors at Q 2 = 0 obtained via two sets of distribution amplitudes differ slightly from each other but they are consistent within the presented errors.
The individual quark and gluon EMTFFs are not re-normalization scale independent. The numerical values of distribution amplitudes are used at the scale µ 2 = 1 GeV 2 in Ref. [57]. In the present study our estimations correspond to µ 2 = 1 GeV 2 , as well. Different sources use different scales to calculate the EMTFFs. In order to compare the results, we should bring them in the same re-normalization scale. For this purpose, we use where n f is the number of flavors, µ i is the initial re-normalization scale and Approaches M q 2 (0) J q (0) d q 1 (0)c q (Q 2 ) Lattice QCD [9] 0.89 ± 0.08 0.45 ± 0.04 --Lattice QCD [10] 0.76 ± 0.14 0.38 ± 0.07 --Lattice QCD [11] 0.73 ± 0.07 0.43 ± 0.07 -1.33 ± 0.06 -Lattice QCD [12] 0.73 ± 0.06 0.32 ± 0.03 -2.15 ± 0.10 -Lattice QCD [13] 0.64 ± 0.03 0.34 ± 0.02 -2.78 ± 0.36 -Lattice QCD [14] 0.76 ± 0.01 0.30 ± 0.01 --Lattice QCD [15] 0.80 ± 0.03 0.40 ± 0.02 --Lattice QCD [16] 0.79 ± 0.06 0.43 ± 0.05 --χPT [22] 0.70 ± 0.02 0.32 ± 0.06 -2.38 ± 0.08 -IFF [23] 0.72 0.32 -2.11 -0.14 Skyrme [24] 1 0.   with Λ = 0.2 GeV. Table III displays a comparison of our results at Q 2 = 0 with those of the various theoretical models, Lattice QCD and existing experimental data for d q 1 (0) at re-normalization scale µ 2 = 1 GeV 2 . For the M 2 (0) form factor, almost all approaches give, more or less, similar predictions. For the J q (0) form factor, our estimations are in good agreements, within the errors, with the predictions of Refs. [10-16, 22, 23]. Note that in the χQSM and Skyrme models 2 J q (0) = M q 2 = 1, since there are only quarks and antiquarks to carry the nucleon's angular momentum and they must carry 100% of it. Despite all of the results for d q 1 (0) from different sources have the same sign, there are large discrepancies among the results. Our predictions for both of the DAs are in accord with the JLab data. Our estimations are also consistent, within the errors, with the predictions of Refs. [12,13,22,23,32,36,37,58], but they differ from the other predictions presented in the table. The negative sign of the d q 1 (0) form factor indicates a profound connection with the spontaneous chiral symmetry breaking in QCD (see also [4][5][6]) as well as an attractive relation with the criterion of stability of the nucleon [32]. The values obtained for thec q (0) form factor using both of the DAs in the present study as well as the prediction of IP [38] are quite small, telling us that the quark and gluon subsystems inside the nucleon interact weakly. Our predictions forc q (0), however, differ with the predictions of IFF [23], substantially.
At the end of this section, we would like to compare the behaviors of the M q 2 (Q 2 ), J q (Q 2 ) and d q 1 (Q 2 ) form factors at small values of (Q 2 ) with the Lattice predictions as well as JLab data for the d q 1 (Q 2 ) gravitational form factor. To this end we plot Fig. 3. The Lattice data points in this figure are taken from the LHPC Collaboration [13] and the JLab data form Ref. [39]. For the form factor M q 2 (Q 2 ), our predictions are consistent with those of the Lattice QCD. In the case of J q (Q 2 ) and d q 1 (Q 2 ) form factors the Lattice results suffer from large uncertainties at small values of Q 2 . Our predictions, especially obtained via set-I DAs reproduces most of the JLab data at small values of Q 2 . This can be considered as a good assurance for the behaviors of the EMTFFs with respect to Q 2 at all regions and in particular at Q 2 = 0.  , J q (Q 2 ) and d q 1 (Q 2 ) as functions of Q 2 at lower values of Q 2 and at the scale µ 2 = 1 GeV 2 . The Lattice data points are taken from the LHPC Collaboration [13] and JLab data from [39].

IV. MECHANICAL PROPERTIES OF NUCLEON
Having calculated the form factors M 2 (Q 2 ) and d 1 (Q 2 ), it is straightforward to calculate the pressure p 0 and the energy density E at the center of nucleon as well as estimate the hadron mechanical radius. The related formulas are given as [8]: where z = Q 2 . Our results on the mechanical quantities of the nucleon compared to other existing theoretical predictions are shown in Table IV. One can see from Table IV that our predictions on p 0 , within the errors, are very close to that of [7], however, they differ with other predictions, considerably. Our results on E using both sets of DAs are close to the predictions of [25] and [32], but demonstrate considerable deviations from other presented predictions. Our predictions on r 2 mech using both sets are in good consistencies with the only existing prediction [36] within the presented uncertainties. The presented results and their comparison with probable future experimental data can be very useful in understanding the structure of the nucleon.

V. SUMMARY AND CONCLUDING REMARKS
The energy-momentum tensor or gravitational form factors of nucleon are basic quantities that carry valuable information on different aspects of the nucleon's structure. These are used to calculate the pressure and energy distributions inside the nucleon as well as quantities related to its geometric shape. The EMTFFs are also sources of information on the fractions of the momenta carried by the quarks and gluons as ingredients of the nucleon. They help us know how the total angular momenta of quarks and gluons form the nucleon's spin. They also provide knowledge on the distribution and stabilization of the strong force inside the nucleon. We extracted the EMTFFs of the nucleon by applying the light-cone QCD sum rule formalism and using two different sets of the parameters inside the nucleon's distribution amplitudes. In the calculations, we used the most general interpolating current of the nucleon in terms of its constituent quark fields. We observed that the results don't depend on the choice of the DAs, considerably and the two sets give close results to each other. We found that the EMTFFs of nucleon are best described by a multipole fit function, helped us to extrapolate the results to the regions that the LCSRs results are not reliable and applicable.
We extracted the numerical values of the EMTFFs at Q 2 = 0 and compared the results with the existing theoretical predictions as well as the results of the Lattice QCD and JLab data. We observed a consistency among the theoretical predictions, within the uncertainties, for the values of the form factors M q 2 (0) and J q (0). Our results on these form factors are nicely consistent with the Lattice QCD and chiral perturbation theory predictions. However, there are large discrepancies among the theoretical predictions on the gravitational form factor d q 1 (0). Nevertheless, our prediction is in accord with the JLab data as well as with the predictions of the Lattice QCD, chiral perturbation theory and KM15-fit. We obtained a very small value for thec q (0) form factor referring to a good conservation of the quark part of the EMT current.
We discussed the behavior of the EMTFFs with respect to Q 2 and observed that all form factors approach to zero at large values of Q 2 . Making use of the multipole fit function, we glanced the behavior of the FFs at small values of Q 2 , where we have some experimental data on d q 1 (Q 2 ) provided by JLab. We saw that the fit function considered in the present study well defines most of the JLab data in the interval Q 2 ∈ [0, 0.4] GeV 2 . The Lattice QCD results for d q 1 (Q 2 ) and J q (Q 2 ) suffer from large uncertainties in this region. The behaviors of the form factor M q 2 (Q 2 ) obtained using two sets of DAs in the present study are well consistent with the Lattice QCD predictions that contain small uncertainties at small values of Q 2 .
Making use of the fit functions of the form factors, we calculated the pressure and energy density distributions at the center of nucleon as well as the mechanical radius of the nucleon and compared with the existing theoretical predictions. Our predictions on p 0 , using two sets of DAs, are very close to that of χQSM [7] within the errors, however, they differ with other predictions presented in table IV, considerably. Our predictions on E at the center of the nucleon are close to the predictions of Skyrme model [25] and χQSM [32], but demonstrate considerable deviations from other presented predictions. The predictions of the present study on r 2 mech are in good consistencies with the only existing prediction provided by LCSR-LO approach [36] within the presented uncertainties.
The presented results in this study together with the predictions of Lattice QCD and other theoretical predictions on the nucleon's EMTFFs may help experimental groups to measure the values of these form factors at a wide range of Q 2 . The good consistency between our predictions and the existing JLab data on d q 1 (Q 2 ) in the interval Q 2 ∈ [0, 0.4] GeV 2 , strengthens this hope. Any experimental data on the energy momentum tensor as well as the electromagnetic, axial and other form factors of the nucleon and their comparison with the theoretical predictions can help us gain valuable knowledge on the internal structures and geometric shapes of the nucleons as building blocks of the visible matter. Such investigations may also help us answer many fundamental questions by means of the quark-gluon structures of the nucleons. By exploring the Gordon equality 2M Nū ′ γ α u =ū ′ (iσ αβ ∆ β + 2P α )u an alternative decomposition of Eq. (1) is obtained: where A(Q 2 ) = M 2 (Q 2 ), A(Q 2 ) + B(Q 2 )) = 2 J(Q 2 ), C(Q 2 ) = 1 5 d 1 (Q 2 ).
Appendix B: Explicit forms of the F functions for the M q 2 (Q 2 ) form factor