Estimation of nucleon D-term in QCD

TL;DR

This paper introduces a method using light-cone sum rules to estimate the nucleon D-term gravitational form factor in QCD, providing preliminary upper bounds and comparing with experimental data and models.

Contribution

It presents a novel approach to estimate the nucleon D-term in QCD using light-cone sum rules at leading order, offering preliminary upper bounds.

Findings

Preliminary upper estimates for the nucleon D-term form factor.

Comparison with experimental data shows consistency within bounds.

Discussion of differences between model predictions and estimated values.

Abstract

Using the light-cone sum rules at leading order, we present an approach to perform the preliminary upper estimation for the nucleon gravitational form factor $D(t)$ ($D$-term contribution). Comparison with the experimental data and with the results of different models is discussed.

Equations

$${\Theta }_q^{\mu \nu }(0)=\frac{i}{2}{\mathcal{R}}_{\tau =2}^{\mu \nu }$$

$$-2i{\Theta }_q^{\mu \nu }(0)$$

$$\displaystyle\widetilde{\cal R}^{\mu\nu}_{\tau=2}=\lim_{y\to 0}\frac{\partial}{\partial y_{\nu}}\,\int\limits_{0}^{1}du\,\frac{\partial}{\partial y_{\mu}}\big{[}\bar{\psi}(0)\,y_{\alpha}\gamma^{\alpha,+}\,[0\,;\,uy]_{A}\,\psi(uy)\big{]}.$$

$${\mathcal{R}}_{\tau =2}^{+-}=$$

$$=$$

$$\displaystyle\langle P^{\prime}|\Theta^{(q)}_{\mu\nu}(0)|P\rangle=\bar{N}(P^{\prime})\Big{[}\mathds{A}(t)\frac{\overline{P}_{\mu}\overline{P}_{\nu}}{m_{N}}+i\mathds{J}(t)\frac{\overline{P}_{\left\{\mu\right.}\sigma_{\left.\nu\right\}\Delta}}{m_{N}}+$$

$$\displaystyle\mathds{D}(t)\frac{\Delta_{\mu}\Delta_{\nu}-g_{\mu\nu}\Delta^{2}}{4m_{N}}+g_{\mu\nu}m_{N}\overline{\mathds{C}}(t)\Big{]}N(P),$$

$$\displaystyle\mathds{J}(t)=\frac{1}{2}\Big{(}\mathds{A}(t)+\mathds{B}(t)\Big{)},\quad a_{\left\{\mu\right.}b_{\left.\nu\right\}}=\frac{1}{2}\big{(}a_{\mu}b_{\nu}+a_{\nu}b_{\mu}\big{)},\quad$$

$$\displaystyle\overline{P}=\frac{1}{2}\big{(}P^{\prime}+P\big{)},\quad\Delta=P^{\prime}-P,\quad\Delta^{2}=-t.$$

$$P=p+\frac{m_N^2}{2}n,\Delta =P\cdot \Delta n+{\Delta }_{\perp }.$$

$$T_{[{\bar{\psi }}_{+}{\psi }_{+}]}^{\mu \nu }(P,\Delta )\equiv ⟨P^{\prime }|{\tilde{\mathcal{R}}}_{\tau =2}^{\mu \nu }|P⟩=\underset{y\to 0}{\lim }\frac{\partial }{\partial y_{\nu }}\underset{0}{\overset{1}{\int }}du\frac{\partial }{\partial w_{\mu }}\int (d^{4}z)e^{-i\Delta \cdot z}⟨0|T\eta (0)$$

$$\displaystyle\times\big{[}\bar{\psi}(w+z)\,w^{\alpha}\gamma^{\alpha,+}\,[w+z\,;\,-w+z]_{A}\,\psi(-w+z)\big{]}|P\rangle$$

$$\displaystyle\eta(0)=\varepsilon^{ijk}\big{[}u^{i}(0)C\gamma_{\alpha}u^{j}(0)\big{]}\gamma_{5}\gamma_{\alpha}d^{k}(0).$$

$$T_{[{\bar{\psi }}_{+}{\psi }_{+}]}^{\mu \nu (d)}(P,\Delta )=$$

$$\displaystyle\frac{1}{8}\int{\cal D}x_{i}\Big{(}2p_{\mu}\big{[}2x_{1}P+\Delta\big{]}_{\nu}+(\mu\leftrightarrow\nu)\Big{)}T^{(d)\,+}_{em}(x_{i};P,\Delta),$$

$$T_{em}^{(q)+}(x_i;P,\Delta )=$$

$$\displaystyle\Big{\{}m_{N}\,{\cal A}^{+\,(q)}_{em}(x_{i};\Delta^{2},P^{\prime\,2})+\hat{\Delta}_{\perp}\,{\cal B}^{+\,(q)}_{em}(x_{i};\Delta^{2},P^{\prime\,2})\Big{\}}N^{+}(P).$$

$$\displaystyle 2p_{\mu}\big{[}2x_{i}P+\Delta\big{]}_{\nu}+(\mu\leftrightarrow\nu).$$

$$\displaystyle\mathds{F}(t)=2m_{N}\Big{(}m_{N}\overline{\mathds{C}}(t)-\frac{\Delta^{2}_{\perp}}{4m_{N}}\mathds{D}(t)\Big{)}.$$

$$\displaystyle\frac{\lambda_{1}m_{N}}{m^{2}_{N}-P^{\prime\,2}}\Big{(}m^{2}_{N}\mathds{A}(t)+\Big{[}\frac{\Delta^{2}_{\perp}}{4}+P\cdot\Delta\Big{]}\mathds{A}(t)+\frac{\Delta^{2}_{\perp}}{4}\mathds{B}(t)$$

$$\displaystyle+\mathds{F}(t)\Big{)}N^{+}(P)=m_{N}\Big{(}\bar{T}^{+-}_{1\,[\bar{\psi}_{+}\psi_{+}]}+\bar{T}^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}\Big{)}N^{+}(P)$$

$$\displaystyle m_{N}\Big{(}\bar{T}^{+-}_{1\,[\bar{\psi}_{+}\psi_{+}]}+\bar{T}^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}\Big{)}=m_{N}\bar{T}^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}+$$

$$\displaystyle\frac{m_{N}}{\pi}\int_{0}^{s_{0}}\frac{ds}{s-P^{\prime\,2}}\text{Im}\Big{\{}\frac{m^{2}_{N}}{4}\mathcal{A}^{++}(s,t)+\frac{P\cdot\Delta}{4}\mathcal{A}^{+}_{em}(s,t)\Big{\}}.$$

$$\displaystyle\mathds{F}(t)=T^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}-\int\hat{d}\mu(s)\Big{\{}\frac{m^{2}_{N}}{4}\mathcal{A}^{++}-$$

$$\displaystyle\frac{\Delta^{2}_{\perp}}{4}\Big{(}\frac{1}{2}\mathcal{A}^{++}+\mathcal{B}^{++}-\frac{1}{2}\mathcal{A}^{+}_{em}\Big{)}\Big{\}},$$

$$\displaystyle\int\hat{d}\mu(s)\big{\{}\mathcal{F}\big{\}}=\frac{1}{\lambda_{1}\pi}\int_{0}^{s_{0}}ds\,e^{(m_{N}^{2}-s)/M^{2}}\,\text{Im}\big{\{}\mathcal{F}(s,t)\big{\}}.$$

$$\displaystyle\mathds{F}(t)=2\,\bar{T}^{+-}_{4\,[\bar{\psi}_{-}\psi_{-}]}-\int\hat{d}\mu(s)\Big{\{}\frac{m^{2}_{N}}{2}\Big{(}\mathcal{A}^{++}-\mathcal{B}^{++}\Big{)}-$$

$$\displaystyle\frac{\Delta^{2}_{\perp}}{4}\Big{(}2\mathcal{B}^{++}-\mathcal{B}^{+}_{em}\Big{)}\Big{\}}.$$

$$\displaystyle R(t)=\frac{4}{\Delta^{2}_{\perp}}\,\Big{\{}2\bar{T}^{+-}_{4\,[\bar{\psi}_{-}\psi_{-}]}-\bar{T}^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}\Big{\}}$$

$$\displaystyle R(t)=\int\hat{d}\mu(s)\Big{\{}\frac{\mathcal{A}^{++}}{2}-\mathcal{B}^{++}\Big{\}}\Big{(}1+\frac{2m^{2}_{N}}{\Delta^{2}_{\perp}}\Big{)}-\int\hat{d}\mu(s)\Big{\{}\frac{\mathcal{A}^{+}_{em}}{2}-\mathcal{B}^{+}_{em}\Big{\}}.$$

$$2\underset{t\to \infty }{\lim }{\bar{T}}_{4[{\bar{\psi }}_{-}{\psi }_{-}]}^{+-}\approx \underset{t\to \infty }{\lim }{\bar{T}}_{1[{\bar{\psi }}_{-}{\psi }_{-}]}^{+-}.$$

$$\displaystyle\bar{T}^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}\stackrel{{\scriptstyle t\to\infty}}{{\approx}}\int\hat{d}\mu(s)\Big{\{}\frac{1}{2}\mathcal{A}^{+}_{em}-\frac{1}{2}\mathcal{A}^{++}\Big{(}1-\frac{2m^{2}_{N}}{t}\Big{)}\Big{\}}.$$

$$\displaystyle\mathds{D}(t)\approx-\int\hat{d}\mu(s)\Big{\{}\frac{1}{2}\mathcal{B}^{++}+\frac{m^{2}_{N}}{t}\mathcal{A}^{++}\Big{\}}.$$

$$\text{mathds}{D}^{\text{fit}}(t)=\frac{D}{(1+et{)}^n},$$

$$D=-2.7,e=0.96{\text{GeV}}^{-2},n=3.43$$

$$p_0=-\frac{1}{24{\pi }^2{m}_N}{\int }_0^{\infty }dtt\sqrt{t}\text{mathds}{D}^{\text{fit}}(t)=0.79\text{GeV}{\text{fm}}^{-3}$$

Full text

Gravitational form factors, QCD, Electromagnetic form factors, Light-cone sum rules, $D$-term

††preprint: Contribution to the Special Issue “QCD and Hadron Structure” devoted to the memory of Prof. G.V.  Efimov

I Introduction

It is well-known that the hadron matrix element of energy-momentum tensor (EMT) can provide information on fundamental characteristics of particles such as mass and spin Polyakov:2002yz ; Polyakov:2018zvc ; Polyakov:2018exb ; Belitsky:2005qn ; Diehl:2003ny . One of the EMT form factors has been related to the $D$-term Kivel:2000fg which can be called as the last unknown fundamental hadron characteristic determining the spatial deformations as well as defining the mechanical properties of hadrons Polyakov:2002yz ; Polyakov:2018zvc ; Polyakov:2018exb . The analogy of $D$-term with the vacuum cosmological constant has been studed in Teryaev:2013qba . The $D$-term has also been calculated in the framework of the dispersion relation approach where the good agreement with the chiral quark-soliton model has been achieved Pasquini:2014vua . Recently, the results of Polyakov:2018zvc have been extended to the different frames where the nucleon has the non-vanishing average momentum Lorce:2018egm .

Since the energy-momentum tensor is the interplay between the gravitation as a external field and the matter fields (in full analogy with the gauge fields which interact with particles through the corresponding electromagnetic current), we adopt the light-cone sum rules (LCSRs) developed for the different nucleon electromagnetic form factors in Anikin:2013aka ; Anikin:2015ita ; Anikin:2016teg .

In the present paper, we develop the light-cone sum rules for the preliminary upper estimation $M$ of the nucleon gravitational form factor $\mathds{D}(t)$, $|\mathds{D}(t)|\leq M$. Unfortunately, the full information on the gravitation form factor $\mathds{D}(t)$ ($D$-term form factor) cannot be obtained within the LCSRs due to the lack of the sum rules that use the minus and perpendicular light-cone projections of electromagnetic current. To estimate $\mathds{D}(t)$ we study the valence quark contributions in nucleon that originate from the leading collinear twist-$2$ combination of the electromagnetic current. In order to approach the region of small $t$ , where the LCSRs approach is really useless, we first approximate the form factor derived for the region of large $t$ by the appropriate multipole function and, then, we do the extrapolation of fitted form factor to the region of small $t$. Our final preliminary prediction demonstrates reasonably good agreement with the first experimental data Hagler:2007xi ; Burkert:2018bqq together with the chiral quark-soliton and Skyrme model results for $D$-term contributions Polyakov:2018zvc .

II Energy-momentum tensor

For the quark contribution, the Belinfanté improved energy-momentum tensor is nothing but the local geometrical twist-$2$ operator which reads

[TABLE]

and it can be expressed through the non-local operator as (cf. Balitsky:1987bk )

[TABLE]

where $y=(y^{+},y^{-},\vec{\bf y}_{\perp})$. From now on, we do not show the trace subtractions. In the paper, we are restricted by the bilinear quark combinations with the spin projection $s_{a}=+1$. Hence, we deal with the collinear twist-$2$ quark combination $[\bar{\psi}_{+}\,\psi_{+}]$ of ${\cal R}^{\mu\nu}_{\tau=2}$. In this connection, we also introduce the local geometrical twist-2 operator which reads

[TABLE]

Notice that this operator actually differs from the operator ${\cal R}^{\mu\nu}_{\tau=2}$ but $\widetilde{\cal R}^{++}_{\tau=2}={\cal R}^{++}_{\tau=2}.$

Moreover, the operator ${\cal R}^{+-}_{\tau=2}$ can be written as

[TABLE]

where the first term with the quark combination involving $[\bar{\psi}_{+}\,\psi_{+}]$ defines $\widetilde{\cal R}^{+-}_{\tau=2}$ and the second term with $[\bar{\psi}_{-}\,\psi_{-}]$ is kept intact because it is beyond the direct computations within our approach. Despite, the $[\bar{\psi}_{-}\,\psi_{-}]$ contribution to the projection $\Theta^{+-}$, which we eventually need, remains unavailable for the explicit LCSRs calculations, we however can still make some reasonable estimation for this kind of contribution.

III Gravitational form factors

We now remind the standard definitions for the relevant gravitational form factors. The quark contribution to the hadron matrix element of the Belinfanté improved energy-momentum tensor operator reads (see, for instance, Polyakov:2018exb )

[TABLE]

where

[TABLE]

Notice that the Belinfanté improved energy-momentum or angular momentum tensor includes the contribution from the spin momentum tensor. As usual, the hadron momenta can be expanded as Anikin:2013aka

[TABLE]

To calculate the form factor $\mathds{D}(t)$ we have to concentrate on a projection ${\cal R}^{+-}_{\tau=2}$. Since our exact computations are restricted by the consideration of $[\bar{\psi}_{+}\,\psi_{+}]$ combination, we are dealing with the corresponding projections of $\widetilde{\cal R}^{\mu\nu}_{\tau=2}$ rather than the projections of ${\cal R}^{\mu\nu}_{\tau=2}$. With respect to the $[\bar{\psi}_{-}\,\psi_{-}]$ combination, we propose a reliable recipe for estimations of this contribution.

In analogy with the nucleon electromagnetic form factors Anikin:2013aka ; Anikin:2015ita ; Anikin:2016teg , we define the amplitude which corresponds to the hadron matrix element of the energy-momentum tensor as

[TABLE]

where $w=uy$. $\eta(0)$ stands for the Ioffe interpolation current defined as

[TABLE]

Below, the colour indices $i,j,k,..$ are omitted.

By the direct calculation of (III), using the distribution amplitude parametrization (see for example Anikin:2013aka ), one concludes that in order to extend the preceding calculations performed for electromagnetic form factors to the case of gravitational form factors we have to weight the electromagnetic form factors with the certain tensor structure. So, the exact expression for the amplitude for $d$-quark contribution is given by

[TABLE]

where

[TABLE]

and ${\cal A}^{+\,(q)}_{em}$, ${\cal B}^{+\,(q)}_{em}$ have been taken from Anikin:2013aka . In a similar manner we can write down the $u$-quark contribution.

From eqn. (10), we conclude that the tensor structure of the amplitude (III) is entirely determined by the tensor

[TABLE]

Since the Ioffe current $\eta(0)$ has been used as the interpolation current, we focus on the valence quark contributions to nucleon. Notice that all amplitudes and form factors should be understood as the objects where the summation over $u$- and $d$-flavours have been implemented: $T_{\mu\nu}=T^{(u)}_{\mu\nu}+T^{(d)}_{\mu\nu}$ etc.

IV Estimation of the form factor $\mathds{D}(t)$

We are now going over to the plus-minus light-cone projection of the amplitude which is $T^{+-}(P,\Delta)$. For the sake of shortness, we introduce the combination given by

[TABLE]

For our goal, we first focus on the $T_{1}$-type (or $\mathcal{A}$-type) of amplitudes, we have the following relations (for the LCSR notations, see Anikin:2013aka )

[TABLE]

where $\mathds{A}(t)$ and $\mathds{B}(t)$ can be expressed through the plus-plus projection of the amplitudes, $\mathcal{A}^{++}$ and $\mathcal{B}^{++}$ respectively, and

[TABLE]

Notice that unlike the plus-plus light-cone projection, the eqn. (IV) do not form the closed system of equations. Hence, there is no way to compute the form factors $\mathds{D}(t)$ directly. Instead, we shall try to estimate the form factor $\mathds{D}(t)$.

The next step is to perform the Borel transforms of eqn. (IV). Afterwards, we derive the first representation for $\mathds{F}(t)$ which reads

[TABLE]

where

[TABLE]

In the same manner, working with the $T_{4}$-type (or $\mathcal{B}$-type) of amplitudes, we derive the the second equivalent representation for $\mathds{F}(t)$ which is

[TABLE]

Combining (IV) and (IV), we obtain the relation as

[TABLE]

where, by introduction,

[TABLE]

The behavior of $R$ as functions of $t$ has been presented in Fig. 1. The numerical analysis, based on our work Anikin:2013aka , shows that $R\approx 0$, therefore we are able to estimate the unknown LCSRs contributions as

[TABLE]

As one can see in Fig. 1, this estimation works with a good precision for $t\geq 4\,\text{GeV}^{2}$.

Further, since the amplitudes of $\mathcal{A}$- and $\mathcal{B}$-types are nothing but the invariant amplitudes of $T_{1}$- and $T_{4}$-types, having focused on $\bar{T}^{+-}_{1\,[\bar{\psi}_{-}\psi_{-}]}$ only, we impose that

[TABLE]

We are now going over to the analysis of the $D$-term. Using the QCD equations of motion, we can express $\overline{\mathds{C}}(t)$ through the nucleon matrix elements of the quark-gluon operator (for the quark contribution) or the gluon-gluon operator (for the gluon contribution) Tanaka:2018wea . Therefore, the calculation of $\overline{\mathds{C}}(t)$ is not available at the leading order at all. Moreover, this form factor meets the condition $\sum_{a}\overline{\mathds{C}^{a}}(t)=0$ which shows that the total (quark and gluon) energy-momentum tensor is conserved. So, at the moment, the consideration of $\overline{\mathds{C}}(t)$ is not presented in the present paper and is postponed to the forthcoming work. Neglecting $\overline{\mathds{C}}$-term from $\mathds{F}$-combination at the large $t$ and using the approximation (22), we estimate the form factor $\mathds{D}(t)$ by (see, the Eqn. (IV))

[TABLE]

We emphasize that the form factor $\mathds{D}(t)$ is independent one. Eqn. (23) shows that the estimated $\mathds{D}(t)$ as a function of $t$ behaves like a combination of the form factors $\mathcal{A}^{++}$ and $\mathcal{B}^{++}$.

According to the strategy of Anikin:2013aka , the gravitational form factor $\mathds{D}(t)$ is not reachable for a region of small $t$ within LCSRs approach. However, we can do approximate fit for this form factor in the region of large $t$ and, then, perform an extrapolation of the obtained fitting functions to the region of small $t$. Finally, the preliminary upper estimation of $\mathds{D}(t)$ can be described by (cf. Tanaka:2018wea )

[TABLE]

In the region of small $t$, for all sets of LCSR-parametrizations, this approximating function matches, see Fig. 2, with the first experimental data Burkert:2018bqq ; Hagler:2007xi together with the results of the chiral quark-soliton and Skyrme models presented in Polyakov:2018zvc ; Goeke:2007fp ; Goeke:2007fq ; Cebulla:2007ei . Our function $\mathds{D}^{\text{fit}}(t)$ lies below than the quark-soliton prediction and above than the result of Skyrme model.

V Mechanical properties of nucleon

Having calculated the form factor $\mathds{D}(t)$, one computes the pressure in the center of nucleon and the hadron mechanical radius Polyakov:2018zvc . We have

[TABLE]

for the pressure $[\text{cf.}\,\,p_{0}^{\text{\cite[cite]{\@@bibref{Authors Phrase1YearPhrase2}{Polyakov:2018zvc}{\@@citephrase{(}}{\@@citephrase{)}}}}}=0.23\,\text{GeV}\text{fm}^{-3}]$, and

[TABLE]

for the mechanical radius.

VI Conclusions

The main challenge in the QCD description of any form factors is the calculation of soft contributions which correspond to the Feynman mechanism to transfer the rather large momentum. As demonstrated in many papers (see, for example, Balitsky:1986st ; Balitsky:1989ry ; Braun:2006hz ; Braun:2001tj ) the LCSRs approach is attractive due to the fact that the soft contributions are calculated in terms of the same DAs that the pQCD calculations include. Thus, the LCSRs can be positioned as one of the most direct relations of the different kind of form factors and hadron DAs that is available at present.

In the paper, we have presented the first preliminary estimation for the valence quark contributions to the gravitational form factor $\mathds{D}(t)$ (the so-called $D$-term) implemented due to the LCSRs techniques at the leading order. To the region of large $t$ where the LCSRs approach can be used, we have observed that the estimated form factor $\mathds{D}$ as a function of $t$ has the similar behaviour as the suitable combination of the well-known gravitational form factors $\mathds{A}(t)$ and $\mathds{B}(t)$. Regarding the small $t$ region, where the first experimental data are available Hagler:2007xi ; Burkert:2018bqq , the estimated $\mathds{D}(t)$ can be approximated by the fitting multipole function $\mathds{D}^{\text{fit}}(t)$. With the upper-estimated $\mathds{D}^{\text{fit}}(t)$, a few quantities characterizing the so-called mechanical properties of the nucleon have been calculated. We emphasize that our result is consistent with the experimental data as well as with the results of the chiral quark-soliton and Skyrme models presented in Polyakov:2018zvc ; Goeke:2007fp ; Goeke:2007fq ; Cebulla:2007ei .

We are grateful to V. Braun, M. Deka, A. Manashov, M. Polyakov and P. Schweitzer for useful discussions, comments and remarks.