Lattice QCD Method To Study Proton Radius Puzzle

TL;DR

This paper develops a lattice QCD approach to calculate the proton's charge radius from first principles, aiming to address the proton radius puzzle caused by experimental discrepancies.

Contribution

It formulates a non-perturbative lattice QCD method to compute the proton charge radius directly from QCD, providing a new tool for resolving experimental disagreements.

Findings

Derived a non-perturbative formula for proton charge radius

Established a lattice QCD framework for the proton radius

Provides a first-principles calculation method

Abstract

Recently there has been disagreement between various experiments about the value of the proton radius which is known as the proton radius puzzle. Since the proton is not a point particle the charge radius of the proton depends on the charge distribution (the form factor) of the partons inside the proton. Since this form factor is a non-perturbative quantity in QCD it cannot be calculated by using the perturbative QCD (pQCD) method but it can be calculated by using the lattice QCD method. In this paper we formulate the lattice QCD method to study the charge radius of the proton. We derive the non-perturbative formula of the charge radius of the proton from the first principle in QCD which can be calculated by using the lattice QCD method.

Equations

$$R_P^2=-\frac{6}{G_E(0)}\frac{d{G}_E(Q^2)}{d{Q}^2}{|}_{Q^2=0}$$

$${\mathcal{O}}_P(x)={ϵ}_{kln}{U}_k^T(x)C{\gamma }_5{D}_l(x)U_n(x)$$

$$<0|{\mathcal{O}}_P^{†}(x^{\prime }){\mathcal{O}}_P(0)|0>=\frac{1}{Z[0]}\int [dA][d\bar{U}][dU][d\bar{D}][dD]\times {\mathcal{O}}_P^{†}(x^{\prime }){\mathcal{O}}_P(0)\times \det [\frac{\delta B_f^s}{\delta {\omega }^b}]$$

$$\times \exp [i\int d^{4}x[-\frac{1}{4}{F}_{\sigma \lambda }^s(x)F^{\sigma \delta s}(x)-\frac{1}{2\alpha }[B_f^s(x){]}^2+{\bar{U}}_k(x)[{\delta }^{kn}(i/\partial -m_U)+g{T}_{kn}^{s}A{/}^s(x)]U_n(x)$$

$$+{\bar{D}}_k(x)[{\delta }^{kn}(i/\partial -m_D)+g{T}_{kn}^{s}A{/}^s(x)]D_n(x)]]$$

$$Z[0]=\int [dA][d\bar{U}][dU][d\bar{D}][dD]\times \det [\frac{\delta B_f^s}{\delta {\omega }^b}]\times \exp [i\int d^{4}x[-\frac{1}{4}{F}_{\sigma \lambda }^s(x)F^{\sigma \delta s}(x)-\frac{1}{2\alpha }[B_f^s(x){]}^2$$

$$+{\bar{U}}_k(x)[{\delta }^{kn}(i/\partial -m_U)+g{T}_{kn}^{s}A{/}^s(x)]U_n(x)+{\bar{D}}_k(x)[{\delta }^{kn}(i/\partial -m_D)+g{T}_{kn}^{s}A{/}^s(x)]D_n(x)]]$$

$$F_{\sigma \lambda }^s(x)={\partial }_{\sigma }{A}_{\lambda }^s(x)-{\partial }_{\lambda }{A}_{\sigma }^s(x)+g{f}^{scd}{A}_{\sigma }^c(x)A_{\lambda }^d(x).$$

$${\mathcal{O}}_P(t,\vec{r})=e^{-iHt}{\mathcal{O}}_P(0,\vec{r})e^{iHt}$$

$$\sum _{l^{\prime \prime }}|H_{l^{\prime \prime }}><H_{l^{\prime \prime }}|=1.$$

$$\sum _{\vec{r}}<0|{\mathcal{O}}_P^{†}(t,\vec{r}){\mathcal{O}}_P(0)|0>=\sum _{l^{\prime \prime }}|<H_{l^{\prime \prime }}|{\mathcal{O}}_P(0)|0>{|}^2{e}^{-\int dt{E}_{l^{\prime \prime }}(t)}$$

$$[\sum _{\vec{r}}<0|{\mathcal{O}}_P^{†}(t,\vec{r}){\mathcal{O}}_P(0)|0>{]}_{t\to \infty }=|<P|{\mathcal{O}}_P(0)|0>{|}^2{e}^{-\int dtE(t)}$$

$$E_P=E(t)+E_{BS}(t)$$

$$\frac{d{E}_{BS}(t^{\prime })}{d{t}^{\prime }}=[\frac{{\sum }_{{\vec{r}}^{\prime \prime }}<0|{\mathcal{O}}_P^{†}(t^{\prime \prime },{\vec{r}}^{\prime \prime })[{\sum }_{q,\bar{q},g}\int d^3{r}^{\prime }{\partial }_i{T}^{i0}(t^{\prime },{\vec{r}}^{\prime })]{\mathcal{O}}_P(0)|0>}{{\sum }_{{\vec{r}}^{\prime \prime }}<0|{\mathcal{O}}_P^{†}(t^{\prime \prime },{\vec{r}}^{\prime \prime }){\mathcal{O}}_P(0)|0>}{]}_{t^{\prime \prime }\to \infty }.$$

$$T^{\sigma \lambda }(x)=F^{\sigma \mu s}(x)F_{\mu }^{\lambda s}(x)+\frac{g^{\sigma \lambda }}{4}{F}^{{\sigma }^{\prime }{\mu }^{\prime }s}(x)F_{{\mu }^{\prime }}^{{\lambda }^{\prime }s}(x)++{\bar{U}}_l(x){\gamma }^{\sigma }[{\delta }^{ln}i{\partial }^{\lambda }-ig{T}_{ln}^s{A}^{\lambda s}(x)]U_n(x)$$

$$+{\bar{D}}_l(x){\gamma }^{\sigma }[{\delta }^{ln}i{\partial }^{\lambda }-ig{T}_{ln}^s{A}^{\lambda s}(x)]D_n(x).$$

$$|<P|{\mathcal{O}}_P(0)|0>{|}^2{e}^{-M_{P}t}=[\frac{{\sum }_{{\vec{r}}^{\prime }}<0|{\mathcal{O}}_P^{†}(t^{\prime },{\vec{r}}^{\prime }){\mathcal{O}}_P(0)|0>}{e^{\int d{t}^{\prime }[\frac{{\sum }_{{\vec{r}}^{\prime \prime }}<0|{\mathcal{O}}_P^{†}({\vec{r}}^{\prime \prime },t^{\prime \prime })[{\sum }_{q,\bar{q},g}\int d{t}^{\prime }\int d^3{r}^{\prime }{\partial }_i{T}^{i0}(t^{\prime },{\vec{r}}^{\prime })]{\mathcal{O}}_P(0)|0>}{{\sum }_{{\vec{r}}^{\prime \prime }}<0|{\mathcal{O}}_P^{†}({\vec{r}}^{\prime \prime },t^{\prime \prime }){\mathcal{O}}_P(0)|0>}{]}_{t^{\prime \prime }\to \infty }}}{]}_{t^{\prime }\to \infty }$$

$$\mathcal{M}=\frac{g^{\lambda \delta }}{q^2}[{\bar{u}}_e(k_f){\gamma }_{\lambda }{u}_e(k_i)][ie{\bar{u}}_P(p_f){\Gamma }_{\delta }(p_f,p_i)u_P(p_i)]=\frac{4\pi \alpha }{Q^2}{l}_{\lambda }{J}^{\lambda }$$

$$J^{\lambda }={\bar{u}}_P(p_f)[{\gamma }^{\lambda }{F}_1(Q^2)+\frac{{\sigma }^{\lambda \delta }{q}_{\delta }}{2M_P}{F}_2(Q^2)]u_P(p_i)=<P|\frac{2}{3}\bar{U}{\gamma }^{\lambda }U-\frac{1}{3}\bar{D}{\gamma }^{\lambda }D|P>.$$

$$G_E(Q^2)=F_1(Q^2)-\frac{Q^2}{4M_P^2}{F}_2(Q^2).$$

$$\sum _{\vec{r}}{e}^{\vec{p}\cdot \vec{r}}<0|{\mathcal{O}}_P^{†}(t,\vec{r}){\mathcal{O}}_P(0)|0>=\sum _{l^{\prime \prime }}|<H_{l^{\prime \prime }}(\vec{p})|{\mathcal{O}}_P(0)|0>{|}^2{e}^{-\int dt{E}_{l^{\prime \prime }}(\vec{p},t)}$$

$$[\sum _{\vec{r}}{e}^{\vec{p}\cdot \vec{r}}<0|{\mathcal{O}}_P^{†}(t,\vec{r}){\mathcal{O}}_P(0)|0>{]}_{t\to \infty }=|<P(\vec{p})|{\mathcal{O}}_P(0)|0>{|}^2{e}^{-\int dtE(\vec{p},t)}$$

$$j_{\lambda }^q(x)=\frac{2}{3}\bar{U}(x){\gamma }_{\lambda }U(x)-\frac{1}{3}\bar{D}(x){\gamma }_{\lambda }D(x)=\sum _f{e}_f{\bar{q}}_f(x){\gamma }_{\lambda }q(x)$$

$$<0|{\mathcal{O}}_P^{†}(x^{\prime \prime })j_0^q(x^{\prime }){\mathcal{O}}_P(0)|0>=\frac{1}{Z[0]}\int [dA][d\bar{U}][dU][d\bar{D}][dD]\times {\mathcal{O}}_P^{†}(x^{\prime \prime })j_0^q(x^{\prime }){\mathcal{O}}_P(0)\times \det [\frac{\delta B_f^s}{\delta {\omega }^b}]$$

$$\times \exp [i\int d^{4}x[-\frac{1}{4}{F}_{\sigma \lambda }^s(x)F^{\sigma \delta s}(x)-\frac{1}{2\alpha }[B_f^s(x){]}^2+{\bar{U}}_k(x)[{\delta }^{kn}(i/\partial -m_U)+g{T}_{kn}^{s}A{/}^s(x)]U_n(x)$$

$$+{\bar{D}}_k(x)[{\delta }^{kn}(i/\partial -m_D)+g{T}_{kn}^{s}A{/}^s(x)]D_n(x)]].$$

$$\sum _{{\vec{r}}^{\prime \prime },{\vec{r}}^{\prime }}{e}^{i\vec{q}\cdot ({\vec{r}}^{\prime \prime }-{\vec{r}}^{\prime })}<0|{\mathcal{O}}_P(t^{\prime \prime },{\vec{r}}^{\prime \prime })j_0^q(t^{\prime },{\vec{r}}^{\prime }){\mathcal{O}}_P(0)|0>=\sum _{n^{\prime \prime },n^{\prime }}<0|{\mathcal{O}}_P|H_{n^{\prime \prime }}(\vec{p})>$$

$$<H_{n^{\prime \prime }}(\vec{p})|j_0^q|H_{n^{\prime }}><H_{n^{\prime }}|{\mathcal{O}}_P|0>e^{-[\int d{t}^{\prime \prime }{E}_{n^{\prime \prime }}(\vec{p},t^{\prime \prime })-\int d{t}^{\prime }{E}_{n^{\prime \prime }}(\vec{p},t^{\prime })]}{e}^{-\int d{t}^{\prime }{E}_{n^{\prime }}(t^{\prime })}$$

$$[\sum _{{\vec{r}}^{\prime \prime },{\vec{r}}^{\prime }}{e}^{i\vec{q}\cdot ({\vec{r}}^{\prime \prime }-{\vec{r}}^{\prime })}<0|{\mathcal{O}}_P(t^{\prime \prime },{\vec{r}}^{\prime \prime })j_0^q(t^{\prime },{\vec{r}}^{\prime }){\mathcal{O}}_P(0)|0>{]}_{t^{\prime \prime }>>>t^{\prime },t^{\prime }\to \infty }=<0|{\mathcal{O}}_P|P(\vec{p})>$$

$$<P(\vec{p})|j_0^q|P><P|{\mathcal{O}}_P|0>e^{-[\int d{t}^{\prime \prime }E(\vec{p},t^{\prime \prime })-\int d{t}^{\prime }E(\vec{p},t^{\prime })]}{e}^{-\int d{t}^{\prime }E(t^{\prime })}$$

$$[\sum _{{\vec{r}}^{\prime }}<0|{\mathcal{O}}_P(t^{\prime },{\vec{r}}^{\prime }){\mathcal{O}}_P(0)|0>{]}_{t^{\prime }\to \infty }=|<P|{\mathcal{O}}_P|0>{|}^2{e}^{-\int d{t}^{\prime }E(t^{\prime })}.$$

$$[\sum _{{\vec{r}}^{\prime \prime }}{e}^{i\vec{p}\cdot {\vec{r}}^{\prime \prime }}<0|{\mathcal{O}}_P({\vec{r}}^{\prime \prime },t^{\prime \prime }-t^{\prime }){\mathcal{O}}_P(0)|0>{]}_{t^{\prime \prime }>>>t^{\prime },t^{\prime }\to \infty }$$

$$=|<P(\vec{p})|{\mathcal{O}}_P|0>{|}^2{e}^{-[\int d{t}^{\prime \prime }E(\vec{p},t^{\prime \prime })-\int d{t}^{\prime }E(\vec{p},t^{\prime })]}.$$

$$<P(\vec{p})|j_0^q|P>=\sqrt{|<P(\vec{p})|{\mathcal{O}}_P|0>{|}^2}\times \sqrt{|<P|{\mathcal{O}}_P|0>{|}^2}$$

$$\times [\frac{{\sum }_{{\vec{r}}^{\prime \prime \prime },{\vec{r}}^{\prime }}{e}^{i\vec{q}\cdot ({\vec{r}}^{\prime }-{\vec{r}}^{\prime \prime \prime })}<0|{\mathcal{O}}_P(t^{\prime },{\vec{r}}^{\prime })j_0^q(t^{\prime \prime \prime },{\vec{r}}^{\prime \prime \prime }){\mathcal{O}}_P(0)|0>}{[{\sum }_{{\vec{r}}^{\prime \prime \prime }}<0|{\mathcal{O}}_P(t^{\prime \prime \prime },{\vec{r}}^{\prime \prime \prime }){\mathcal{O}}_P(0)|0>][{\sum }_{{\vec{r}}^{\prime }}{e}^{i\vec{p}\cdot {\vec{r}}^{\prime }}<0|{\mathcal{O}}_P(t^{\prime }-t^{\prime \prime \prime },{\vec{r}}^{\prime }){\mathcal{O}}_P(0)|0>]}{]}_{t^{\prime }>>>t^{\prime \prime \prime },t^{\prime \prime \prime }\to \infty }.$$

Full text

I Introduction

By using the spectroscopic method involving the electronic hydrogen (the ordinary hydrogen atom consisting of proton and electron) the charge radius of the proton is measured to be 0.8768 $\times 10^{-15}$ meter. Similarly by using the electron-proton scattering method the charge radius of the proton is measured to be 0.8775 $\times 10^{-15}$ meter which is consistent with the spectroscopic method. The CODATA-2014 world average value of the charge radius of the proton by using the electrons, i. e., by using the above two methods, is 0.8751 $\times 10^{-15}$ meter pr .

However, the muonic hydrogen (hydrogen atom consisting of proton and muon) experiment in the year 2010 found that the charge radius of the proton is 0.84087 $\times 10^{-15}$ meter mh ; mh1 . The CODATA-2018 value of the charge radius of the proton is 0.8414 $\times 10^{-15}$ meter cod . This disagreement between various experiments about the value of the charge radius of the proton is known as the proton radius puzzle which remains an unsolved problem in science.

It is well known that the proton is not a point particle but it is a composite particle consisting of quarks and gluons which are the fundamental particles of the nature. The up quark has the fractional electric charge $\frac{2e}{3}$ and the down quark has fractional electric charge $-\frac{e}{3}$ where $e$ is the magnitude of the charge of the electron. Hence the charge radius of the proton depends on the charge distribution (the form factor) of the partons inside the proton.

The electric charge radius $R_{P}$ of the proton $P$ is given by rp

[TABLE]

where $G_{E}(Q^{2})$ is the electric form factor and $Q^{2}$ is the momentum transfer square of the virtual photon in the lepton-proton scattering.

The interaction between the quarks and gluons inside the proton is described by the quantum chromodynamics (QCD) ymj which is a fundamental theory of the nature. The short distance partonic cross section can be calculated by using the perturbative QCD (pQCD) due to the asymptotic freedom in QCD gwj . Using the factorization theorem in QCD fcj ; fcj1 ; fcj2 the hadronic cross section can be calculated from the partonic cross section at the high energy colliders by using the experimentally extracted parton distribution function (PDF) and fragmentation function (FF).

The formation of the proton from the quarks and gluons is a long distance phenomenon in QCD. Due to the asymptotic freedom the QCD coupling becomes large at the large distance where the pQCD is not applicable. Hence the formation of the proton from the quarks and gluons cannot be studied by using the pQCD. The non-perturbative QCD is necessary to study the formation of the proton from the quarks and gluons.

However, the analytical solution of the non-perturbative QCD is not known yet because of the presence of the cubic and quartic power of the gluon fields in the QCD lagrangian inside the path integration in the generating functional in QCD [see section II for details]. The path integration in QCD can be performed numerically in the Euclidean time by using the lattice QCD method. Hence the lattice QCD provides the first principle method to study the formation of the proton from the quarks and gluons.

Since the electric form factor $G_{E}(Q^{2})$ of the partons inside the proton in eq. (1) is a non-perturbative quantity in QCD it cannot be calculated by using the perturbative QCD (pQCD) method but it can be calculated by using the lattice QCD method.

Recently we have formulated the lattice QCD method to study the proton formation from the quarks and gluons pj and to study the proton spin crisis psj and to study the proton decay pd by implementing the non-zero boundary surface term in QCD which arises due to the confinement of quarks and gluons inside the finite size proton nkbsj . In this paper we extend this to formulate the lattice QCD method to study the proton radius puzzle. We derive the non-perturbative formula of the charge radius of the proton from the first principle in QCD which can be calculated by using the lattice QCD method.

The paper is organized as follows. In section II we describe the lattice QCD method to study the proton formation from quarks and gluons by implementing the non-zero boundary surface term in QCD due to confinement. In section III we formulate the lattice QCD method to study the proton radius puzzle and derive the non-perturbative formula of the charge radius of the proton from the first principle in QCD which can be calculated by using the lattice QCD method. Section IV contains conclusions.

II Formation of proton from quarks and gluons Using lattice QCD Method

The partonic operator for the proton ($P$) formation is given by

[TABLE]

where $U_{k}(x)$ is the up quark field and $D_{k}(x)$ is the down quark field, $C$ is the charge conjugation operator and $k,l,n=1,2,3$ are the color indices.

In the path integral formulation of QCD the vacuum expectation value of the non-perturbative partonic correlation function of the type $<0|{\cal O}^{\dagger}_{P}(x^{\prime}){\cal O}_{P}(x^{\prime\prime})|0>$ is given by

[TABLE]

where $B_{f}^{s}(x)$ is the gauge fixing term with color index $s=1,...,8$, the $A_{\sigma}^{s}(x)$ is the gluon field with Lorentz index $\sigma=0,1,2,3$, the $\alpha$ is the gauge fixing parameter, $m_{U}$ is the mass of the up quark, $m_{D}$ is the mass of the down quark and $Z[0]$ is the generating functional in QCD given by

[TABLE]

with

[TABLE]

In eq. (3) we do not have ghost fields as we directly work with the ghost determinant ${\rm det}[\frac{\delta B_{f}^{s}}{\delta\omega^{b}}]$.

The time evolution of the partonic operator in the Heisenberg representation is given by

[TABLE]

where the QCD Hamiltonian of the partons is given by $H$.

The complete set of proton energy-momentum eigenstates is given by

[TABLE]

Using eqs. (6) and (7) in (3) we find in the Euclidean time

[TABLE]

where $\int dt$ is an indefinite integration. In the large time limit we neglect the higher energy level contribution to find

[TABLE]

where $E(t)$ is the energy of all the partons inside the proton in its ground state and $|P>$ is the energy-momentum eigenstate of the proton in its ground state.

Due to non-vanishing boundary surface term $E_{BS}(t)$ in QCD due to the confinement of partons inside the finite size proton we find nkbsj

[TABLE]

where $E_{P}$ is the energy of the proton, $E(t)$ is energy of all the partons inside the proton and $E_{BS}(t)$ is the non-vanishing boundary surface term in QCD given by

[TABLE]

In eq. (11) the energy-momentum tensor density $T^{\sigma\lambda}(x)$ of the partons in QCD is given by

[TABLE]

Using eqs. (10) and (11) in (9) we find

[TABLE]

which can be calculated by using the lattice QCD method where $M_{P}$ is the mass of the proton and $\int dt^{\prime}$ is indefinite integration.

Eq. (13) is the non-perturbative formula to study the proton formation from quarks and gluons using the lattice QCD method by implementing non-zero boundary surface term in QCD due to confinement.

III Lattice QCD Method to study charge radius of the proton

In the previous section we have formulated the lattice QCD method to study the formation of the proton from the quarks and gluons by implementing the non-zero boundary surface term in QCD due to the confinement. In this section we will extend this to formulate the lattice QCD method to study the proton radius puzzle. We will derive the non-perturbative formula of the charge radius of the proton from the first principle in QCD which can be calculated by using the lattice QCD method.

This method is also used to study various quantities in QCD in vacuum psj ; pd ; avj and in QCD in medium amj to study the quark-gluon plasma at RHIC and LHC qk ; qk1 ; qk2 .

In the single (virtual) photon exchange the amplitude ${\cal M}$ for the electron-proton ($eP$) elastic scattering process is given by

[TABLE]

where $k_{i}(k_{f})$, $p_{i}(p_{f})$ are the initial (final) momentum of the electron, proton, $q^{\mu}=p^{\mu}_{f}-p^{\mu}_{i}$ is the momentum of the virtual photon, $u_{e},u_{P}$ are the electron, proton spinors respectively, $l^{\mu},J^{\mu}$ are the leptonic, hadronic (electromagnetic) currents respectively and $q^{2}=-Q^{2}$. The general form of the vertex function $\Gamma^{\lambda}(p_{f},p_{i})$ satisfying the relativistic invariance and current conservation is given by

[TABLE]

The electric form factor $G_{E}(Q^{2})$ is given by rp

[TABLE]

Using eqs. (6) and (7) in (3) we find in the Euclidean time

[TABLE]

where ${\vec{p}}$ is the momentum of the proton and $\int dt$ is an indefinite integration. In the large time limit we neglect the higher energy level contribution in eq. (17) to find

[TABLE]

where $E({\vec{p}},t)$ is the energy of all the partons inside the proton of momentum ${\vec{p}}$ in its ground state and $|P({\vec{p}})>$ is the energy-momentum eigenstate of the proton $P$ of momentum ${\vec{p}}$ in its ground state.

From eq. (15) the electromagnetic current operator $j^{\lambda}_{q}(x)$ of the quarks inside the proton is given by

[TABLE]

where $q_{f}(x)$ and $e_{f}$ are the quark field and the fractional electrical charge of the quark of the flavor $f$ respectively.

In the path integral formulation of QCD the vacuum expectation value of the non-perturbative 3-point partonic correlation function of the type $<0|{\cal O}_{P}(x^{\prime\prime})j^{q}_{0}(x^{\prime}){\cal O}_{P}(0)|0>$ is given by

[TABLE]

In this paper we assume that the initial proton is at rest, i. e., ${\vec{p}}_{i}=0$. This means ${\vec{p}}_{f}={\vec{q}}={\vec{p}}$. Using eqs. (6) and (7) in (20) we find in the Euclidean time

[TABLE]

where ${\vec{p}}$ is the momentum of the final state proton and ${\vec{q}}={\vec{p}}$.

In the limit $t^{\prime\prime}>>>t^{\prime},~{}~{}t^{\prime}\rightarrow\infty$ we find by neglecting the higher energy level contributions

[TABLE]

where $E({\vec{p}},t)$ is the energy of all the partons inside the proton of momentum ${\vec{p}}$, the $E(t)$ is the energy of all the partons inside the proton at rest, $|P({\vec{p}})>$ is the energy-momentum eigenstate of the proton of momentum ${\vec{p}}$ and $|P>$ is the energy-momentum eigenstate of the proton at rest.

From eq. (9) we find for the proton at rest

[TABLE]

Similarly from eq. (18) we find for the proton with momentum ${\vec{p}}$

[TABLE]

From eqs. (22), (23) and (24) we find

[TABLE]

For the initial proton at rest with energy-momentum eigenstate $|P>$ and the final proton of momentum ${\vec{p}}$ with energy-momentum eigenstate $|P({\vec{p}})>$ we find that the electric form factor $G_{E}(Q^{2})$ of the partons inside the proton is related to the proton matrix element $<P({\vec{p}})|j^{q}_{0}|P>$ of the quark electromagnetic current $j^{q}_{0}$ via the relation eff

[TABLE]

where $E_{P}({\vec{p}})$ is the energy of the proton of momentum ${\vec{p}}={\vec{q}}$ and $Q^{2}=-q^{2}$.

From eq. (13) we find for the proton at rest

[TABLE]

where $M_{P}$ is the mass of the proton. Similarly for the proton with momentum ${\vec{p}}$ we find

[TABLE]

where $E_{P}({\vec{p}})$ is the energy of the proton with momentum ${\vec{p}}$.

Using eqs. (27), (28) and (25) in (26) we find

[TABLE]

where the initial proton is at rest and the final proton momentum is ${\vec{p}}={\vec{q}}$ with $Q^{2}=-q^{2}$.

The eq. (29) is the non-perturbative formula of the electric form factor $G_{E}(Q^{2})$ of the proton derived from the first principle in QCD which can be calculated by using the lattice QCD method. The charge radius $R_{P}$ of the proton is obtained from this electric form factor $G_{E}(Q^{2})$ in eq. (29) by using eq. (1).

IV Conclusions

Recently there has been disagreement between various experiments about the value of the proton radius which is known as the proton radius puzzle. Since the proton is not a point particle the charge radius of the proton depends on the charge distribution (the form factor) of the partons inside the proton. Since this form factor is a non-perturbative quantity in QCD it cannot be calculated by using the perturbative QCD (pQCD) method but it can be calculated by using the lattice QCD method. In this paper we have formulated the lattice QCD method to study the charge radius of the proton. We have derived the non-perturbative formula of the charge radius of the proton from the first principle in QCD which can be calculated by using the lattice QCD method.