text:
# Quantized bulk scalar fields in the Randall-Sundrum brane model

## 1 Introduction

The search for a unified field theory led from the very early days after Einstein’s theory of general relativity to the suggestion that there could be extra spatial dimensions beyond the familiar three. The original idea of Kaluza  was to incorporate electromagnetism into a 5-dimensional theory of gravity. In order to explain why the extra spatial dimension has not been observed, Klein  suggested that it was compactified to a circle of extremely small radius, typically taken to be of the order of the Planck length. The generalization to more than one extra spatial dimension  allowed the incorporation of non-abelian gauge fields as part of the higher dimensional metric. These higher dimensional theories are usually referred to by the generic term Kaluza-Klein theories.

In the early 1980’s an immense amount of interest was generated in Kaluza-Klein theories by Witten’s  observation that seven extra spatial dimensions was the minimum required for incorporation of an $`SU(3)\times SU(2)\times U(1)`$ gauge symmetry, and Nahm’s proof  that this was the maximum allowed by supergravity. This focused attention on 11-dimensional theories with various compactifications. (See Ref.  for a review.) Further interest in higher dimensional theories was, and continues to be, generated by superstring theory.

A new angle on Kaluza-Klein theories was pointed out recently . In contrast to the standard belief that extra dimensions must be associated with extremely small length scales, it was noted that the extra dimensions could be much bigger than the Planck scale and still remain as yet undetected. This has the effect of bringing the fundamental Planck scale much closer to the electroweak scale, and thus offers a possible solution to the gauge hierarchy problem. The basic idea behind this  is to have a $`(4+n)`$-dimensional spacetime, with the $`n`$-dimensional space having a volume $`V_n`$. The relationship between the 4-dimensional Planck scale $`M_{Pl}`$ and the fundamental mass scale $`M`$ of the higher dimensional space is $`M_{Pl}V_nM^{n+2}`$. For large volumes $`V_n`$, we can have $`M`$ of the order of the electroweak scale $`(1`$ TeV) and still have $`M_{Pl}10^{19}`$ GeV. Unfortunately this scenario does not provide a completely satisfactory solution to the gauge hierarchy problem since it trades a large ratio between the Planck scale and the electroweak scale for a large ratio between the compactification scale $`V^{1/n}`$ and the electroweak scale.

Randall and Sundrum  suggested a model that overcomes the problem just described. This model is based on a 5-dimensional spacetime with the extra spatial dimension having an orbifold compactification. Two 3-branes with opposite tensions sit at the orbifold fixed points. The line element is

$$ds^2=e^{2kr|\varphi |}\eta _{\mu \nu }dx^\mu dx^\nu r^2d\varphi ^2$$
(1.1)
with $`x^\mu `$ the usual 4-dimensional coordinates, $`|\varphi |\pi `$ with the points $`(x^\mu ,\varphi )`$ and $`(x^\mu ,\varphi )`$ identified. The factor of $`e^{2kr|\varphi |}`$ in the line element is often referred to as the warp factor, since it results in the spacetime not being the direct product of 4-dimensional spacetime and the extra dimension. The 3-branes sit at $`\varphi =0`$ and $`\varphi =\pi `$. $`k`$ is a constant of the order of the Planck scale (the natural scale for the theory), and $`r`$ is an arbitrary constant associated with the size of the extra dimension.

The interesting feature of the Randall-Sundrum model is that it can generate a TeV mass scale from the Planck scale in the higher dimensional theory. A field with a mass $`m_0`$ on the $`\varphi =\pi `$ brane will have a physical mass of $`me^{\pi kr}m_0`$. By taking $`kr12`$, and $`m_010^{19}`$ GeV, we end up with $`m1`$ TeV. In the original version of the Randall-Sundrum model all of the standard model particles are supposed to be confined on the brane, with only gravity in the bulk (5-dimensional) spacetime. Alternatives to confining particles to the brane, and allowing fields to live in the bulk spacetime have been investigated. See for example . It is the possibility of bulk matter fields that we are interested in this paper.

The radius $`r`$ is assumed to be the vacuum expectation value of a scalar field, called the radion. In the original Randall-Sundrum model the value of $`r`$ is undetermined. In order to make the theory physically acceptable a mechanism for determining $`r`$ must be found. In Ref.  a classical method for stabilizing $`r`$ was suggested.

The similar problem of fixing the size of the extra dimensions was found in the older Kaluza-Klein theories that were based on direct products of Minkowski spacetime and a homogeneous space, often chosen to be a sphere. It was necessary to have the extra dimensions of the order of the Planck scale, but there was nothing in the classical theory that determined this. It was recognized by Candelas and Weinberg  that quantum effects from matter fields or gravity could fix the size of the extra dimensions in a natural way. It therefore is sensible to ask if a similar procedure can fix the value of $`r`$ in the Randall-Sundrum model. This problem has received some recent attention  with some differing conclusions. The main purpose of the present paper is to provide details of what we believe to be the correct and complete calculation of the vacuum energy density for bulk scalar fields, and to calculate the induced gravity term in the effective action  since it is known  that this can play an important role in self-consistent solutions and can give different conclusions from merely considering the vacuum energy density. We also clarify the two inequivalent field configurations which can arise as a result of the boundary conditions that can be imposed. Finally, we analyze in detail the requirement that the Randall-Sundrum spacetime be a solution to the quantum corrected Einstein equations. It is shown that this is not just the condition that the Casimir energy density have a minimum.

## 2 Quantum corrections to the classical theory

In order to calculate the quantum corrections to the classical theory we will follow the procedure used in Ref. . This consists of first expanding the higher dimensional field in terms of a complete set of modes appropriate to the extra dimensions. The dependence on the extra dimensions can then be integrated out of the action functional to leave an equivalent 4-dimensional theory, but one with an infinite number of fields. The masses of the fields are usually quantized in some way. This is the procedure used in our computation of the vacuum energy density  for the minimally coupled scalar field.

If we wish to calculate the induced gravity contribution to the effective action we must generalize the metric (1.1) to allow curvature in the 4-dimensional submanifold representing the branes. We will take

$$ds^2=e^{2\sigma (y)}g_{\mu \nu }(x)dx^\mu dx^\nu dy^2$$
(2.1)
where $`y`$ is related to the angle $`\varphi `$ in (1.1) by

$$y=r\varphi $$
(2.2)
and the function $`\sigma (y)`$ is given by

$$\sigma (y)=k|y|.$$
(2.3)
We will allow the submanifold specified by $`y=0`$ to have an arbitrary dimension $`D`$. (The Randall-Sundrum model is then specified by $`D=4`$ and $`g_{\mu \nu }=\eta _{\mu \nu }`$.)

The Einstein-Hilbert gravitational action is

$$S_G=(16\pi G)^1d^Dx𝑑y|\widehat{g}|^{1/2}(\widehat{R}2\mathrm{\Lambda })$$
(2.4)
where we use a caret to denote a $`(D+1)`$-dimensional object. (So $`\widehat{g}_{\widehat{\mu }\widehat{\nu }}`$ is the $`(D+1)`$-dimensional metric found from (2.1), $`\widehat{g}`$ is its determinant, and $`\widehat{R}`$ is the scalar curvature constructed from it.) There are also terms  in the total action functional coming from the two branes at $`y=0,\pi r`$ which represent the hidden and visible branes respectively:

$$S_{brane}=d^Dx𝑑y|\widehat{g}|^{1/2}\left\{V_v\delta (y\pi r)+V_h\delta (y)\right\}.$$
(2.5)
We will return to this part of the action in Sec. 3.

In order that the metric in (2.1) with $`g_{\mu \nu }=\eta _{\mu \nu }`$ satisfy the field equations coming from varying (2.4) and (2.5) with $`\sigma (y)`$ given in (2.3) along with the $`_2`$ identification $`(x^\mu ,y)(x^\mu ,y)`$, it follows that

$`\sigma ^{\prime \prime }`$ $`=`$ $`2k\left[\delta (y)\delta (y\pi r)\right],`$ (2.6)
$`V_h=V_v`$ $`=`$ $`{\displaystyle \frac{(D1)}{4\pi G}}k,`$ (2.7)
$`\mathrm{\Lambda }`$ $`=`$ $`{\displaystyle \frac{1}{2}}D(D1)k^2.`$ (2.8)
These reduce to the results of  when $`D=4`$.

For later use we record the curvature scalar for (2.1):

$`\widehat{R}`$ $`=`$ $`e^{2\sigma }R+D\left[2\sigma ^{\prime \prime }(D+1)(\sigma ^{})^2\right]`$ (2.9)
$`=`$ $`e^{2\sigma }R+4kD\left[\delta (y)\delta (y\pi r)\right]D(D+1)k^2.`$ (2.10)
We will consider a single real scalar field $`\mathrm{\Phi }`$ of mass $`m`$ in the spacetime (2.1):

$$S=\frac{1}{2}d^Dx𝑑y|\widehat{g}|^{1/2}\left\{\widehat{g}^{\widehat{\mu }\widehat{\nu }}_{\widehat{\mu }}\mathrm{\Phi }_{\widehat{\nu }}\mathrm{\Phi }m^2\mathrm{\Phi }^2\xi \widehat{R}\mathrm{\Phi }^2\right\}.$$
(2.11)
We allow a general non-minimal coupling to the curvature with dimensionless coupling constant $`\xi `$. The minimally coupled scalar field has $`\xi =0`$. The conformally coupled scalar field has $`\xi =(D1)/(4D)`$. We will keep $`\xi `$ arbitrary here.

By considering the $`(D+1)`$-dimensional field equation for $`\mathrm{\Phi }`$, following the analysis of Ref. , leads us to define $`f_n(y)`$ as a solution to

$`e^{(D2)\sigma }_y\left[e^{D\sigma }_yf_n\right]+e^{2\sigma }\overline{m}^2f_n`$
$`+4\xi Dke^{2\sigma }[\delta (y)\delta (y\pi r)]f_n=m_n^2f_n,`$ (2.12)
normalized by

$$\underset{\pi r}{\overset{+\pi r}{}}𝑑ye^{(2D)\sigma }f_n(y)f_n^{}(y)=\delta _{nn^{}}.$$
(2.13)
We have defined

$$\overline{m}^2=m^2D(D+1)k^2\xi .$$
(2.14)
We first obtain the boundary conditions satisfied by $`f_n(y)`$ by integrating (2.12) first about $`y=0`$ and then about $`y=\pi r`$. We assume that $`f_n(y)`$ is continuous at these two points. This gives the two results

$`\underset{ϵ0}{lim}\left\{{\displaystyle \frac{f_n(y)}{y}}|_{y=ϵ}{\displaystyle \frac{f_n(y)}{y}}|_{y=ϵ}\right\}`$ $`=`$ $`4\xi kDf_n(0),`$ (2.15)
$`\underset{ϵ0}{lim}\left\{{\displaystyle \frac{f_n(y)}{y}}|_{y=\pi rϵ}{\displaystyle \frac{f_n(y)}{y}}|_{y=\pi rϵ}\right\}`$ $`=`$ $`4\xi kDf_n(\pi r),`$ (2.16)
(In the second result we have used the $`yy`$ identification.)

Next we solve (2.12) for $`y>0(y\pi r)`$ and $`y<0(y\pi r)`$. Call the solutions in these two regions $`f_n^\pm (y)`$. It is easy to show that

$$f_n^\pm (y)=e^{\pm \frac{D}{2}ky}\left\{A_n^\pm J_\nu \left(\frac{m_n}{k}e^{\pm ky}\right)+B_n^\pm Y_\nu \left(\frac{m_n}{k}e^{\pm ky}\right)\right\}$$
(2.17)
where $`A_n^\pm `$ and $`B_n^\pm `$ are constants of integration. The order of the Bessel functions $`\nu `$ is given by

$$\nu =\sqrt{\frac{\overline{m}^2}{k^2}+\frac{D^2}{4}}.$$
(2.18)
At this stage we must choose how $`f_n^+(y)`$ is related to $`f_n^{}(y)`$ by the $`_2`$ identification on $`S^1`$. The situation is like that used some time ago  in considering fields on non-simply connected spacetimes. For a real scalar field there are only two possibilities. We have either

$$f_n(y)=f_n(y),$$
(2.19)
corresponding to an untwisted field, or else

$$f_n(y)=f_n(y),$$
(2.20)
corresponding to a twisted field. The first choice corresponds to the boundary conditions adopted in Ref.  and Ref. . We examine each possibility.

For the untwisted field (2.19) leads to

$$A_n^{}=A_n^+,B_n^{}=B_n^+$$
(2.21)
in which case we can write

$$f_n(y)=e^{\frac{D}{2}k|y|}\left\{A_nJ_\nu \left(\frac{m_n}{k}e^{k|y|}\right)+B_nY_\nu \left(\frac{m_n}{k}e^{k|y|}\right)\right\}$$
(2.22)
for either positive or negative values of $`y`$ (but $`y0,\pi r`$), as given in Ref. . Applying the boundary conditions (2.15) and (2.16) results in the eigenvalues $`m_n`$ determined by solutions to the transcendental equation

$$F_\nu \left(\frac{m_n}{ka}\right)=0,$$
(2.23)
with

$$F_\nu (z)=\frac{2}{\pi }\left[y_\nu (az)j_\nu (z)j_\nu (az)y_\nu (z)\right],$$
(2.24)
where $`j_\nu (z)`$ and $`y_\nu (z)`$ are expressed in terms of Bessel functions by

$`j_\nu (z)`$ $`=`$ $`{\displaystyle \frac{1}{2}}D(14\xi )J_\nu (z)+zJ_\nu ^{}(z),`$ (2.25)
$`y_\nu (z)`$ $`=`$ $`{\displaystyle \frac{1}{2}}D(14\xi )Y_\nu (z)+zY_\nu ^{}(z).`$ (2.26)
We have written

$$a=e^{k\pi r}.$$
(2.27)
For $`\xi =0`$ and $`D=4`$ these expressions reduce to the results of Ref. .

For the twisted field (2.20) leads to

$$A_n^{}=A_n^+,B_n^{}=B_n^+$$
(2.28)
instead of (2.21). This time we have

$$f_n(y)=\frac{y}{|y|}e^{\frac{D}{2}k|y|}\left\{A_nJ_\nu \left(\frac{m_n}{k}e^{k|y|}\right)+B_nY_\nu \left(\frac{m_n}{k}e^{k|y|}\right)\right\}$$
(2.29)
as the solution for $`y>0(y\pi r)`$ and $`y<0(y\pi r)`$. Because $`f_n(y)`$ for the twisted field is an odd function of $`y`$, using the $`yy`$ identification, we must have

$$f_n(0)=0,f_n(\pi r)=0.$$
(2.30)
These last two results are consistent with the general boundary conditions (2.15) and (2.16) since the derivative of an odd function is even, resulting in the left hand sides of the two equations vanishing. The boundary conditions (2.30) lead to a much simpler analysis, with the eigenvalues $`m_n`$ given as solutions to

$$\stackrel{~}{F}_\nu \left(\frac{m_n}{ka}\right)=0$$
(2.31)
with

$$\stackrel{~}{F}_\nu (z)=\frac{2}{\pi }\left[J_\nu (z)Y_\nu (az)J_\nu (az)Y_\nu (z)\right].$$
(2.32)
Thus there can be no conflict with the results of Ref.  that used untwisted boundary conditions. As we discuss in Appendix B, there are some slight differences in the technicalities between the twisted and untwisted scalars.

We now expand the $`(D+1)`$-dimensional scalar field $`\mathrm{\Phi }(x,y)`$ in terms of the $`f_n(y)`$ as

$$\mathrm{\Phi }(x,y)=\underset{n}{}\phi _n(x)f_n(y).$$
(2.33)
This expansion can be substituted into (2.11) and the integration performed over $`y`$. The result is

$$S=\frac{1}{2}\underset{n}{}d^Dx|g|^{1/2}\phi _n(\mathrm{}+m_n^2+\xi R)\phi _n.$$
(2.34)
This shows that we may view the theory as one in $`D`$ dimensions, with an infinite number of fields whose masses are quantized according to (2.23).

At this stage we can work out the one-loop effective action in the usual way: simply compute the contribution from the $`n^{th}`$ mode, and then sum over all modes. Because we are after the vacuum energy, as well as the induced gravity term, it is advantageous to use a heat kernel method, and we follow Ref.  here. (See  for a review.)

The one-loop effective action for the $`n^{th}`$ mode is

$$\mathrm{\Gamma }_n^{(1)}=\frac{i}{2}d^Dx|g|^{1/2}\underset{0}{\overset{\mathrm{}}{}}\frac{ds}{s}e^{ism_n^2}K(s;x,x)$$
(2.35)
where

$`i{\displaystyle \frac{}{s}}K(s;x,x^{})`$ $`=`$ $`(\mathrm{}+\xi R)K(s;x,x^{}),`$ (2.36)
$`K(0;x,x^{})`$ $`=`$ $`\delta (x,x^{}).`$ (2.37)
The result for $`\mathrm{\Gamma }_n^{(1)}`$ is infinite and requires regularization. We choose dimensional regularization here with $`D`$ viewed as $`D+ϵ`$ with $`ϵ0`$ understood. We now make use of the asymptotic expansion (see Ref.  for example)

$$K(s;x,x)i\mathrm{}^ϵ(4\pi is)^{D/2}\underset{k=0}{\overset{\mathrm{}}{}}(is)^ka_k,$$
(2.38)
where

$`a_0`$ $`=`$ $`1,`$ (2.39)
$`a_1`$ $`=`$ $`\left({\displaystyle \frac{1}{6}}\xi \right)R.`$ (2.40)
Higher order terms in (2.38) involve higher powers of the curvature tensor; so the $`a_0`$ and $`a_1`$ terms contain the leading order terms.

The two leading terms in $`\mathrm{\Gamma }_n^{(1)}`$ are found to be

$`\mathrm{\Gamma }_n^{(1)}`$ $``$ $`{\displaystyle \frac{1}{2}}(4\pi )^{D/2}\mathrm{}^ϵ{\displaystyle }d^Dx|g|^{1/2}\{\mathrm{\Gamma }(D/2)(m_n^2)^{D/2}`$ (2.41)
$`+({\displaystyle \frac{1}{6}}\xi )\mathrm{\Gamma }(1D/2)(m_n^2)^{D/21}R\}.`$
If we write

$`\mathrm{\Gamma }^{(1)}`$ $`=`$ $`{\displaystyle \underset{n}{}}\mathrm{\Gamma }_n^{(1)}`$ (2.42)
$`=`$ $`{\displaystyle d^Dx|g|^{1/2}\left\{_\mathrm{\Lambda }^{(1)}+_R^{(1)}+\mathrm{}\right\}}`$
we have

$`_\mathrm{\Lambda }^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}(4\pi )^{D/2}\mathrm{}^ϵ\mathrm{\Gamma }(D/2){\displaystyle \underset{n}{}}(m_n^2)^{D/2},`$ (2.43)
$`_R^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}(4\pi )^{D/2}\mathrm{}^ϵ\left({\displaystyle \frac{1}{6}}\xi \right)R\mathrm{\Gamma }(1D/2){\displaystyle \underset{n}{}}(m_n^2)^{D/21},`$ (2.44)
The term in $`_\mathrm{\Lambda }^{(1)}`$ represents a contribution to the cosmological constant in $`D`$ dimensions induced by quantum corrections to the classical theory. $`_\mathrm{\Lambda }^{(1)}`$ is recognized as the negative of the Casimir, or vacuum energy density, and corresponds to what we calculated earlier  in the case $`\xi =0`$ and $`D=4`$. The term in $`_R^{(1)}`$ represents an induced gravity term which can renormalize the Newtonian gravitational constant found in the classical gravitational action. In the special case $`\xi =1/6`$, which corresponds to the conformal value for the coupling constant when $`D=3`$, we have $`_R^{(1)}=0`$. In this case the lowest order correction to the classical gravitational action would result in an $`R^2`$ theory. For the original 5-dimensional Randall-Sundrum metric, $`D=4`$, and the conformal value is $`\xi =3/16`$, so there is an induced gravity term in the effective action.

## 3 Renormalization

The expressions for $`_\mathrm{\Lambda }^{(1)}`$ and $`_R^{(1)}`$ obtained in (2.43) and (2.44) may be evaluated using the method described in Appendix A or B. With $`\omega (s)`$ defined in (A.1) we have

$`_\mathrm{\Lambda }^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}(4\pi )^{D/2}\mathrm{}^ϵ(ka)^D\omega (D/2),`$ (3.1)
$`_R^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}(4\pi )^{D/2}\mathrm{}^ϵ\left({\displaystyle \frac{1}{6}}\xi \right)R(ka)^{D2}\omega (1D/2).`$ (3.2)
The pole parts of these expressions will determine the form of the required counterterms. If we use $`\mathrm{PP}\{\mathrm{}\}`$ to denote the pole part of any expression, from (A.45) and (A.46) we find

$`\mathrm{PP}\left\{_\mathrm{\Lambda }^{(1)}\right\}`$ $`=`$ $`{\displaystyle \frac{1}{ϵ}}(4\pi )^{D/2}{\displaystyle \frac{k^D}{\mathrm{\Gamma }(D/2)}}\left[a^D+(1)^D\right]d_D,`$ (3.3)
$`\mathrm{PP}\left\{_R^{(1)}\right\}`$ $`=`$ $`{\displaystyle \frac{1}{2ϵ}}(4\pi )^{D/2}\left({\displaystyle \frac{1}{6}}\xi \right)R{\displaystyle \frac{(D2)k^{D2}}{\mathrm{\Gamma }(D/2)}}`$ (3.4)
$`\times \left[a^{D2}+(1)^D\right]d_{D2}.`$
The first few coefficients $`d_k`$ are given in (A.27A.30) for untwisted fields, and in (B.5B.10) for twisted fields. In both cases higher order coefficients may be calculated in a routine manner, but because they become lengthy we will not give the results.

The renormalization can be effected by modifying the procedure used by Ref. , which used the brane tensions. Because we allow curvature on the branes, it is necessary to augment the brane tension terms with curvature or else it will not be possible to remove the pole terms coming from (3.4) when $`D`$ is even. As noted in Ref. , the key feature is the $`a`$-dependence of the results. When $`D`$ is odd, (3.3) and (3.4) show that the pole terms are proportional to $`a^D1`$ and $`a^{D2}1`$ respectively. Such pole terms can be dealt with by a renormalization of the bulk gravitational and cosmological terms. The sign difference between the factors of $`a`$ and the 1 comes about from the integration over the extra $`y`$ coordinate to obtain the effective 4-dimensional theory. When $`D`$ is even ($`D=4`$ for the Randall-Sundrum model) the pole terms are proportional to $`a^D+1`$ and $`a^{D2}+1`$ respectively. It is not possible to use the bulk terms to remove these divergences. The situation is easily understood from the viewpoint of heat-kernel methods and dimensional regularization, and was pointed out a long time ago . In an odd spacetime dimension (as occurs for $`D=4`$) the only possible counterterms come from the boundary of the spacetime. For the Randall-Sundrum model this corresponds to the two branes.

The metric on the visible brane is

$$g_{\mu \nu }^v(x)=\widehat{g}_{\mu \nu }(x,y=\pi r)=a^2g_{\mu \nu }.$$
(3.5)
The metric on the hidden brane is

$$g_{\mu \nu }^h(x)=\widehat{g}_{\mu \nu }(x,y=0)=g_{\mu \nu }.$$
(3.6)
These give rise to

$`R_v`$ $`=`$ $`a^2R,`$ (3.7)
$`R_h`$ $`=`$ $`R`$ (3.8)
as the curvature scalars on the visible and hidden branes We will take

$$S_{brane}=S_{brane}^v+S_{brane}^h$$
(3.9)
where

$`S_{brane}^v`$ $`=`$ $`{\displaystyle d^Dx|g_v|^{1/2}(\kappa _vR_vV_v)},`$ (3.10)
$`S_{brane}^h`$ $`=`$ $`{\displaystyle d^Dx|g_h|^{1/2}(\kappa _hR_hV_h)},`$ (3.11)
as the brane action. $`\kappa _v,\kappa _h,V_v,V_h`$ are all bare quantities. It follows that

$$S_{brane}=d^Dx|g|^{1/2}\{(\kappa _h+a^{D2}\kappa _v)R(V_h+a^DV_v)\}.$$
(3.12)
All bare quantities can be expressed in terms of renormalized ones and counterterms in the usual way :

$`\kappa _{v,h}`$ $`=`$ $`\kappa _{v,h}^R+\delta \kappa _{v,h},`$ (3.13)
$`V_{v,h}`$ $`=`$ $`V_{v,h}^R+\delta V_{v,h},`$ (3.14)
This gives a counterterm brane action $`\delta S_{brane}`$ from (3.12) in an obvious way. Because the structure of the pole terms in (3.3) and (3.4) is the same as that of $`\delta S_{brane}`$, there is no problem absorbing them into the counterterms.

There is always the freedom to perform finite renormalizations since the renormalization counterterms should be fixed by imposing some renormalization conditions on the effective action. (See Ref.  for example.) To this end we can write

$`\delta \kappa _{v,h}`$ $`=`$ $`{\displaystyle \frac{\delta \kappa _{v,h}^1}{ϵ}}+\delta \kappa _{v,h}^0,`$ (3.15)
$`\delta V_{v,h}`$ $`=`$ $`{\displaystyle \frac{\delta V_{v,h}^1}{ϵ}}+\delta V_{v,h}^0,`$ (3.16)
where $`\delta \kappa _{v,h}^1,\delta \kappa _{v,h}^0,\delta V_{v,h}^1,\delta V_{v,h}^0`$ are all independent of $`ϵ`$. Noting that the factors of $`a^D`$ in (3.12) are really $`a^{D+ϵ}`$, we find

$`\delta _{brane}`$ $`=`$ $`[{\displaystyle \frac{1}{ϵ}}(\delta \kappa _h^1+a^{D2}\delta \kappa _v^1)`$ (3.17)
$`(\delta \kappa _h^0+a^{D2}\delta \kappa _v^0+a^{D2}\delta \kappa _v^1\mathrm{ln}a)]R`$
$`[{\displaystyle \frac{1}{ϵ}}(\delta V_h^1+a^D\delta V_v^1)+\delta V_h^0`$
$`+a^D\delta V_v^0+a^D\delta V_v^1\mathrm{ln}a]`$
as the Lagrangian density counterterms on the branes. $`D`$ is now fixed to be the physical spacetime dimension (eg 4). In order that the poles in (3.3) and (3.4) be removed we must fix

$`\delta V_v^1`$ $`=`$ $`\left({\displaystyle \frac{k^2}{4\pi }}\right)^{D/2}{\displaystyle \frac{d_D}{\mathrm{\Gamma }(D/2)}},`$ (3.18)
$`\delta V_h^1`$ $`=`$ $`(1)^D\delta V_v^1,`$ (3.19)
$`\delta \kappa _v^1`$ $`=`$ $`{\displaystyle \frac{1}{2}}\left(\xi {\displaystyle \frac{1}{6}}\right){\displaystyle \frac{(D2)k^{D2}}{\mathrm{\Gamma }(D/2)}}d_{D2},`$ (3.20)
$`\delta \kappa _h^1`$ $`=`$ $`(1)^D\delta \kappa _v^1.`$ (3.21)
If we use the results contained in the appendices for either twisted or untwisted fields, (see (A.44,A.45)), it can be seen that $`_\mathrm{\Lambda }^{(1)}`$ involves either terms independent of $`a`$, proportional to $`a^D`$ or $`a^D\mathrm{ln}a`$, apart from the term we have defined as $`A_1`$. Apart from $`A_1`$, the $`a`$-dependence is the same form as the finite counterterms $`\delta V_h^0+a^D\delta V_v^0+a^D\delta V_v^1\mathrm{ln}a`$ appearing in $`\delta _{brane}`$. This means that all terms in $`_\mathrm{\Lambda }^{(1)}`$ apart from $`A_1`$ can be absorbed by finite counterterms. A similar conclusion holds for the curvature dependent part of $``$. The renormalized expression for $`_\mathrm{\Lambda }`$ becomes

$$_\mathrm{\Lambda }=V_h^Ra^DV_v^R+\frac{1}{2}(4\pi )^{D/2}(ka)^DA_1(D/2).$$
(3.22)
The renormalized expression for $`_R`$ becomes

$$_R=\left[\kappa _h^R+a^{D2}\kappa _v^R+\frac{1}{8\pi }\left(\frac{1}{6}\xi \right)(4\pi )^{1D/2}(ka)^{D2}A_1(1D/2)\right]R.$$
(3.23)
The function $`A_1(s)`$ was defined in (A.37) with (A.21) used for $`g_\nu (z)`$ in the untwisted field case, and (B.3) used for twisted fields.

In the next section we look at the possibility of finding self-consistent solutions arising from quantum corrections to the classical theory.

## 4 Self-consistent solutions?

Having obtained the vacuum, or Casimir, energy for the effective $`D`$-dimensional theory we can now see if there are any solutions of the Randall-Sundrum type to the quantum corrected field equations. We will argue that the problem is not quite as simple as calling the Casimir energy a potential and looking for its extremum. Rather, we must require that the Einstein field equations hold, and we will see that this has a different requirement.

If $`\mathrm{\Gamma }`$ denotes the full effective action with quantum effects included we must require

$`{\displaystyle \frac{\delta \mathrm{\Gamma }}{\delta g_{\mu \nu }}}`$ $`=`$ $`0,`$ (4.1)
$`{\displaystyle \frac{\delta \mathrm{\Gamma }}{\delta a}}`$ $`=`$ $`0,`$ (4.2)
with the Randall-Sundrum spacetime (2.1) as a solution with $`g_{\mu \nu }=\eta _{\mu \nu }`$. The classical part of $`\mathrm{\Gamma }`$ comes from (2.4) and (2.5), and to these contributions we must add the one-loop correction $`\mathrm{\Gamma }^{(1)}`$. If we perform the integration over the extra dimension $`y`$ using the metric ansatz (2.1) with $`\sigma `$ given by (2.3) we obtain

$`S_G`$ $`=`$ $`(16\pi G)^1{\displaystyle }dv_x\{{\displaystyle \frac{2}{(D2)k}}(1a^{D2})R+4kD(1a^D)`$ (4.3)
$`2k(D+1)(1a^D){\displaystyle \frac{4\mathrm{\Lambda }}{kD}}(1a^D)\}.`$
$`dv_x=d^Dx|g|^{1/2}`$ is the $`D`$-dimensional invariant volume element. Because we set $`g_{\mu \nu }=\eta _{\mu \nu }`$ after differentiation with respect to $`g_{\mu \nu }`$, the term in the curvature $`R`$ will make no contribution to the field equations. The same will be true for the induced gravity part of the one-loop effective action, as well as for higher order curvature terms that we have not calculated. The induced gravity term is important for relating the overall coefficient of $`R`$ in the effective action to the physical value of the gravitational constant, which is irrelevant to the field equations in the present case. We are therefore justified in only considering the vacuum energy in the present calculation. Accordingly we will take

$$\mathrm{\Gamma }=S_G𝑑v_xF(a)$$
(4.4)
where $`F(a)`$ includes the brane tension terms from (2.5) as well as the Casimir energy.

Variation of $`\mathrm{\Gamma }`$ with respect to $`g_{\mu \nu }`$ followed by taking $`g_{\mu \nu }=\eta _{\mu \nu }`$ results in

$$k(D1)(1a^D)\frac{2\mathrm{\Lambda }}{kD}(1a^D)=8\pi GF(a).$$
(4.5)
Variation of $`\mathrm{\Gamma }`$ with respect to $`a`$ yields

$$k(D1)\frac{2\mathrm{\Lambda }}{kD}=\frac{8\pi G}{D}a^{1D}F^{}(a),$$
(4.6)
after setting $`g_{\mu \nu }=\eta _{\mu \nu }`$. Combining (4.5) and (4.6) gives us a requirement for $`F(a)`$ :

$$0=(1a^D)aF^{}(a)+Da^DF(a).$$
(4.7)
As a check on this last result, if we set

$$F(a)=V_va^D+V_h,$$
(4.8)
corresponding to the brane tension terms in (2.5), we find the condition

$$V_v+V_h=0.$$
(4.9)
This reproduces the result on the balancing of the two brane tensions needed in the Randall-Sundrum model. The remaining requirements (2.7) and (2.8) follow upon using (4.9) in (4.5) and (4.6).

If we now include the vacuum energy by taking

$$F(a)=V_va^D+V_h+f(a),$$
(4.10)
in place of the classical expression (4.8), we must have

$$0=Da^D(V_v+V_h)+(1a^D)af^{}(a)+Da^Df(a)$$
(4.11)
holding. This is clearly not just a simple minimization condition on the potential. From (3.22) we have

$$f(a)=\frac{1}{2}(4\pi )^{D/2}k^Da^DA_1(D/2)$$
(4.12)
with $`A_1`$ defined in (A.37). If we take

$$I(a)=A_1(D/2),$$
(4.13)
then we must have

$$0=\frac{2D(4\pi )^{D/2}}{k^D}(V_v+V_h)DI(a)a(1a^D)I^{}(a).$$
(4.14)
This condition must be fulfilled if the Randall-Sundrum spacetime is to be a self-consistent solution to the effective 4-dimensional Einstein equations with quantum effects included.

In order to satisfy the self-consistency relation (4.14), we can use the freedom to tune the brane tensions. A possible choice is to require that the balancing condition (4.9) between the brane tensions still holds. For this choice it is easy to see that there are no solutions to (4.14). To see this is quite simple. It is easy to show that $`I(a)`$ is negative for $`0<a<1`$, and is monotonically decreasing. The right hand side of (4.14) can never vanish. However, there is no need to require that $`V_v+V_h=0`$ when quantum effects are included<sup>1</sup><sup>1</sup>1We would like to thank J. Garriga for having pointed this out to us., and relaxing this condition allows solutions to (4.14). The interesting thing is that it is possible to find self-consistent solutions for values of $`a`$ which solve the hierarchy problem only at the price of severe fine tuning of the brane tensions<sup>2</sup><sup>2</sup>2This has been noted in  for the case of massless minimally coupled and conformally coupled fields. This is a more general result than  who only were interested in self-consistent solutions which had values for $`a`$ consistent with the solution to the gauge hierarchy problem (ie $`a<<1`$), and who only examined the Casimir, or vacuum, energy.

## 5 Massless, conformally coupled fields

The case of massless, conformally coupled scalars deserves special mention. Taking $`m^2=0`$ and $`\xi =(D1)/(4D)`$ in the order of the Bessel functions determined in (2.18) gives $`\nu =1/2`$. Because Bessel functions of half-integer order can be given in terms of elementary functions, the analysis becomes much simpler than for general $`\nu `$. If we examine the twisted fields first, and use

$`J_{1/2}(z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{2}{\pi z}}}\mathrm{sin}z,`$ (5.1)
$`Y_{1/2}(z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{2}{\pi z}}}\mathrm{cos}z,`$ (5.2)
it is easy to see that $`\stackrel{~}{F}_{1/2}(z)`$ defined in (2.32) is given by

$$\stackrel{~}{F}_{1/2}(z)=\frac{4\mathrm{sin}(1a)z}{\pi ^2\sqrt{a}z}.$$
(5.3)
We therefore have the mass eigenvalues given by

$$m_n=\frac{\pi ka}{(1a)}n,\mathrm{for}n=1,2,3,\mathrm{}.$$
(5.4)
We can evaluate the sum over modes in $`_\mathrm{\Lambda }^{(1)}`$ directly using the Riemann $`\zeta `$-function, without using the contour integral method of Appendix B. It is easy to see that

$`_\mathrm{\Lambda }^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}(4\pi )^{D/2}\left({\displaystyle \frac{\pi ka}{1a}}\right)^D\mathrm{}^ϵ\mathrm{\Gamma }(D/2)\zeta (D)`$ (5.6)
$`=`$ $`{\displaystyle \frac{1}{2\sqrt{\pi }}}(4\pi )^{D/2}\left({\displaystyle \frac{ka}{1a}}\right)^D\mathrm{\Gamma }\left({\displaystyle \frac{D+1}{2}}\right)\zeta (D+1).`$
In the last line we have let $`ϵ0`$, since the result is analytic at $`ϵ=0`$, and used some of the relations of the $`\mathrm{\Gamma }`$\- and $`\zeta `$-functions to simplify the result. Although no infinte renormalizations are required, it is still necessary to do finite ones.

As a check on this result, we can now apply the method of Appendix B. We need to know $`A_1(D/2)`$ which is given by

$`A_1(D/2)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1+D/2)}}{\displaystyle \underset{0}{\overset{\mathrm{}}{}}}𝑑zz^D{\displaystyle \frac{d}{dz}}\mathrm{ln}\stackrel{~}{g}_{1/2}(z)`$ (5.8)
$`=`$ $`{\displaystyle \frac{2}{\mathrm{\Gamma }(D/2)}}{\displaystyle \underset{0}{\overset{\mathrm{}}{}}}𝑑zz^{D1}\mathrm{ln}\left[{\displaystyle \frac{1e^{2(1a)z}}{1e^{2z}}}\right].`$
(In the second line we have integrated by parts and made use of the Bessel functions (5.1,5.2) for $`\nu =1/2`$. ) The integral is easily evaluated in terms of the Riemann $`\zeta `$-function with the result

$$A_1(D/2)=\frac{1}{\sqrt{\pi }}\mathrm{\Gamma }\left(\frac{D+1}{2}\right)\zeta (D+1)\left[(1a)^D1\right].$$
(5.9)
Note that when expanded in powers of $`a`$, $`A_1`$ begins with order $`a`$ as found in the case of general $`\nu `$. The result in (5.9) gives rise to the same result as found by a direct summation of the mass eigenvalues, after imposing the same renormalization conditions, and provides a useful check on the method.

The untwisted conformally coupled scalar field case was also dealt with in Ref.  who exploited the conformal invariance to relate the problem to one in flat spacetime with boundaries present. It is worth emphasizing that this requires the theory to have no conformal anomaly, or else extra terms will arise when transforming back to the original spacetime. The fact that no such terms do arise is guaranteed by our method which makes no use of conformal transformations. In effect we have shown that there is no conformal anomaly for conformally coupled scalars in the Randall-Sundrum spacetime in any dimension. Vanishing of the conformal anomaly is guaranteed by the vanishing of $`\stackrel{~}{d}_k`$ when $`\nu =1/2`$. (This is easily checked to be true for the coefficients listed in (B.5B.10).)

We now look at the case of untwisted fields. The functions $`j_\nu `$ and $`y_\nu `$ defined in (2.25,2.26) are easily seen to be

$`j_{1/2}(z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{2z}{\pi }}}\mathrm{cos}z,`$ (5.10)
$`y_{1/2}(z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{2z}{\pi }}}\mathrm{sin}z,`$ (5.11)
We then find

$$F_{1/2}(z)=\frac{4\sqrt{a}}{\pi ^2}z\mathrm{sin}(1a)z,$$
(5.12)
where $`F_{1/2}(z)`$ was defined in (2.24). The mass eigenvalues are given by the same expression (5.4) as we had in the twisted field case. A direct summation over the mass eigenvalues again leads to (5.6).

We can also check the result by using the contour integral method of Appendix A. The integral $`A_1(D/2)`$ becomes

$$A_1(D/2)=\frac{1}{\sqrt{\pi }}\mathrm{\Gamma }\left(\frac{D+1}{2}\right)\zeta (D+1)\left[(1a)^D+12^D\right].$$
(5.13)
The only difference with the twisted field result is the constant term. The $`(1a)^D`$ dependence is the same. Note that as a consequence, when $`a`$ is small, $`A_1`$ is approximately a constant, rather than proportional to $`a`$. This would seem to contradict the general results established in Appendix A. The root cause of this lies in the behaviour of $`F_{1/2}(z)`$ as $`z0`$. Clearly $`F_{1/2}(z)z^2`$ as $`z0`$, which does not satisfy the assumption we made in Appendix A that $`F_{1/2}(0)`$ is finite and non-zero. In this particular case we must remove the zero of $`F_{1/2}(z)`$ at $`z=0`$ and take instead $`z^2F_{1/2}(z)`$ which does satisfy the conditions of being finite and non-zero as $`z0`$ allowing the deformation of contour used in Appendix A. If we compare $`z^2F_{1/2}(z)`$ with $`\stackrel{~}{F}_{1/2}(z)`$ it is obvious that the same result will now be obtained for $`A_1(D/2)`$ for the untwisted field as we found for the twisted field, and which is in agreement with the direct summation over mass modes. In any case, the difference between the twisted and untwisted fields for $`A_1`$ results in a term proportional to $`a^D`$ in the effective action which can be dealt with by renormalization.

For either twisted or untwisted field configurations it is easy to show that the only solutions to (4.14) is $`a=0`$ for $`V_v+V_h=0`$. When this last condition is relaxed, it is possible to find self-consistent solutions, but if we demand the hierarchy problem to be solved a severe fine tuning of the brane tensions is required.

## 6 Discussion and conclusions

In the present paper we have examined the one-loop radiative corrections to the effective action arising from the quantization of a massless scalar field non-minimally coupled to the curvature of the background spacetime proposed by Randall and Sundrum . Although it is not possible to obtain an exact result for general, curved versions of the Randall-Sundrum spacetime, it is possible to obtain an expansion in powers of the curvature. The vacuum (Casimir) energy is the first term in this expansion. The vacuum energy has received attention  due to its possible role in stabilising the radius of the extra dimension in a self-consistent way, as an alternative to the mechanism proposed in . Because of the different conclusions and inconsistent results obtained in previous analyses, we have presented a careful and complete calculation here. Differences and comparisons of our results are made with Refs.  at appropriate stages.

One important clarification concerns the boundary conditions adopted for the scalar field. We emphasized that there were two inequivalent field configurations (conventionally referred to as twisted and twisted ). Either can be used, and in general the results will be different.

Another clarification concerns the renormalisation of the theory. We presented a full calculation for both field configurations for the generally coupled massive scalar field with all relevant counterterms indicated. To do this it is necessary to perform a full analytic continuation of relevant integrals as outlined in Appendices A and B. As we emphasized, the same method can be used for both twisted and untwisted fields, but the precise details of the analytic continuation are slightly different due to the different asymptotic behaviour of the Bessel functions between the two field configurations. Ref.  did not perform the detailed analytic continuation, and really only analyzed the untwisted field configuration. Nevertheless, as we showed above, the essential conclusions of their paper concerning the renormalisation are correct in either case. One reason for presenting the detailed calculation of the counterterms is that the pole term in dimensional regularization is related to the conformal anomaly for the scalar field.

As noted in the older Kaluza-Klein models based on spaces with homogeneous extra dimensions , a full analysis of self-consistency can only be achieved if the quantum correction to the gravitational part of the effective action is known. We presented the necessary analysis for this above. In order to renormalize the theory when the branes were curved, it was necessary to add on additional contributions to the brane tensions involving the curvature. This suggest that the brane sector of the model has to be viewed as a series, with the brane tensions as the first term in the expansion.

Another clarification concerns the conditions for self-consistency of the Randall-Sundrum spacetime. With flat branes we showed that the induced gravity contribution was not important, in contrast to the situation in . Additionally, it was argued that the requirement of self-consistency is not obtained simply by looking at the minimum of the Casimir energy. By looking at the quantum corrected Einstein equations for the effective 4-dimensional (or more generally $`D`$-dimensional) theory we obtained a consistency condition which the radius of the extra dimension had to satisfy. We found no solutions to this if we demand $`V_v+V_h=0`$, meaning that it is not possible to stabilize the radius by quantum effects coming from bulk scalar fields in this case. This conclusion can be altered by relaxing the balancing condition $`V_v+V_h=0`$, even if we demand that the gauge hierarchy problem also be solved, but only at the price of fine tuning the brane tensions.

We also examined massless conformally coupled scalars fields as a limiting case of our general results. This provides a useful check on the results, since it is possible to perform the calculations in a different, and much simpler, way. It also provides a comparison with Ref.  where a conformal transformation was used to relate the calculation to one in a flat spacetime with boundaries. A curious feature of the massless conformally coupled case is that the imposition of twisted or untwisted boundary conditions makes no difference to the final result. As we mentioned above, our calculation provides a proof that there is no conformal anomaly for the Randall-Sundrum spacetime in any number of dimensions.

An interesting point concerns the relation between viewing the quantized scalar field as an effective 4-dimensional theory, and a 5-dimensional theory. From the 5-dimensional viewpoint, self-consistency must be satisfied from the 5-dimensional Einstein field equations. To do this, it is necessary to evaluate not just the total energy, but rather the energy density. (More specifically we must know the vacuum expectation value of the stress-energy tensor.) Due to the boundary conditions and the nature of the spacetime, the vacuum energy density for the 5-dimensional scalar field will not be expected to be constant in general. It would be interesting to study this problem, although the presence of Bessel functions in the mode eigenfunctions would make this rather difficult.

Final comments concern the whole question of using quantum effects to obtain a self-consistent result for the radius of the extra dimensions. The conclusions of Ref.  were rather negative concerning the possibility of simultaneously solving the gauge hierarchy problem and self-consistency. It is still interesting to ask if it is possible to obtain self-consistent solutions at all, regardless of the size of the extra dimensions. All that we currently know about this is that it is possible to do so with bulk scalar fields, but, fine tuning of the brane tensions is required when we demand that the hierarchy problem is simultaneously solved.

Whether or not a theory with more realistic matter content can do so is an open question. Finally, even if it is possible to find a self-consistent solution there still remains the problem of determining if it is stable.

## Acknowledgements

A. Flachi is grateful to the University of Newcastle upon Tyne for the award of a Ridley Studentship. We are grateful to J. Garriga for comments concerning Ref. , and to I. G. Moss for discussions on the conformally coupled case.

## Appendix A Evaluation of the sums for untwisted fields

In this appendix we outline the evaluation of the sum over modes needed to calculate the vacuum energy and induced gravity terms. Begin by defining

$$\omega (s)=\mathrm{\Gamma }(s)\underset{n=1}{\overset{\mathrm{}}{}}x_n^{2s},$$
(A.1)
where $`x_n`$ is the $`n^{\mathrm{th}}`$ positive solution to $`F_\nu (x_n)=0`$ with $`F_\nu (z)`$ defined in terms of Bessel functions in (2.24). For large $`n`$ we know $`x_nn`$ so that (A.1) is defined initially for $`\mathrm{}(s)>1/2`$ where the sum converges. We want to analytically continue $`\omega (s)`$ to $`s=D/2`$ and $`1D/2`$. The basic method for doing this is to convert the sum into a contour integral. A simple application of the residue theorem shows that

$$\omega (s)=\frac{\mathrm{\Gamma }(s)}{2\pi i}_𝒞𝑑zz^{2s}\frac{d}{dz}\mathrm{ln}F_\nu (z),$$
(A.2)
where $`𝒞`$ is any contour which encloses the positive zeros of $`F_\nu (z)`$ defined in (2.24). A convenient choice is shown in Fig. 1.

This procedure is basically the same as that used in Refs.  for studying the Casimir effect in a spherical shell. Ref.  uses an integration by parts to remove the derivative in (A.2), and it is this basic starting point that Refs.  follow. It is more convenient to stick with (A.2) for our purposes, although either method could be used. A difference between our approach and that of Ref.  is that we do not need to use the uniform asymptotic expansion of the Bessel functions, making our calculation slightly simpler.

By using the behaviour of the Bessel functions for small $`z`$ it can be shown that the contribution from the semicircle of radius $`r`$ behaves like $`r^{22s}`$ for small $`r`$. This means that if we restrict ourselves to the strip $`\frac{1}{2}<\mathrm{}(s)<1`$ the contribution to the contour coming from this portion of the contour will vanish as $`r0`$. (The fact that $`F_\nu (0)`$ is finite and non-zero is crucial to this result.) For the semicircle of radius $`R`$, the integrand behaves like $`R^{2s}`$ for large $`R`$, which follows from using the large argument expansion for the Bessel functions in $`F_\nu (z)`$. The contribution from this portion of the contour also vanishes if we let $`R\mathrm{}`$ and restrict ourselves to the strip $`\frac{1}{2}<\mathrm{}(s)<1`$. The only contributions to the integral for $`\frac{1}{2}<\mathrm{}(s)<1`$ come from the portions of the contour along the positive and negative imaginary axis. After a short calculation it can be shown that

$$\omega (s)=\frac{1}{\mathrm{\Gamma }(1s)}\underset{0}{\overset{\mathrm{}}{}}𝑑zz^{2s}\frac{d}{dz}\mathrm{ln}P_\nu (z),$$
(A.3)
with

$$P_\nu (z)=\frac{2}{\pi }\left[i_\nu (az)k_\nu (z)i_\nu (z)k_\nu (az)\right].$$
(A.4)
and

$`i_\nu (z)`$ $`=`$ $`{\displaystyle \frac{D}{2}}(14\xi )I_\nu (z)+zI_\nu ^{}(z),`$ (A.5)
$`k_\nu (z)`$ $`=`$ $`{\displaystyle \frac{D}{2}}(14\xi )K_\nu (z)+zK_\nu ^{}(z).`$ (A.6)
We now wish to analytically continue the expression (A.3) for $`\omega (s)`$ out of the strip $`\frac{1}{2}<\mathrm{}(s)<1`$ to $`\mathrm{}(s)<\frac{1}{2}`$. The simplest way to do this is to note that the impediment to convergence of (A.3) for $`\mathrm{}(s)<\frac{1}{2}`$ comes from the behaviour of $`P_\nu (z)`$ at large $`z`$. We will define

$`I_\nu (z)`$ $`=`$ $`{\displaystyle \frac{e^z}{\sqrt{2\pi z}}}\mathrm{\Sigma }^{(I)}(z),`$ (A.7)
$`K_\nu (z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{\pi }{2z}}}e^z\mathrm{\Sigma }^{(K)}(z),`$ (A.8)
which gives

$`i_\nu (z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{z}{2\pi }}}e^z\mathrm{\Sigma }_\nu ^{(I)}(z),`$ (A.9)
$`k_\nu (z)`$ $`=`$ $`\sqrt{{\displaystyle \frac{\pi z}{2}}}e^z\mathrm{\Sigma }_\nu ^{(K)}(z),`$ (A.10)
from (A.5) and (A.6) where

$`\mathrm{\Sigma }_\nu ^{(I)}(z)`$ $`=`$ $`\mathrm{\Sigma }^{(I)}(z)+{\displaystyle \frac{d}{dz}}\mathrm{\Sigma }^{(I)}(z)+{\displaystyle \frac{1}{2z}}[D(14\xi )1]\mathrm{\Sigma }^{(I)}(z),`$ (A.11)
$`\mathrm{\Sigma }_\nu ^{(K)}(z)`$ $`=`$ $`\mathrm{\Sigma }^{(K)}(z){\displaystyle \frac{d}{dz}}\mathrm{\Sigma }^{(K)}(z){\displaystyle \frac{1}{2z}}[D(14\xi )1]\mathrm{\Sigma }^{(K)}(z).`$ (A.12)
For large $`z`$ the asymptotic expansions for the Bessel functions show that

$$\mathrm{\Sigma }^{(I)}(z)\underset{k=0}{\overset{\mathrm{}}{}}\alpha _kz^k,$$
(A.13)
where

$$\alpha _k=\frac{(1)^k\mathrm{\Gamma }(\nu +k+\frac{1}{2})}{2^kk!\mathrm{\Gamma }(\nu k+\frac{1}{2})},$$
(A.14)
and

$$\mathrm{\Sigma }^{(K)}(z)\mathrm{\Sigma }^{(I)}(z).$$
(A.15)
Because the divergent behaviour of $`\omega (s)`$ comes from the large $`z`$ behaviour of the integrand, we will split up

$$\omega (s)=\omega _1(s)+\omega _2(s),$$
(A.16)
where

$`\omega _1(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}P_\nu (z),`$ (A.17)
$`\omega _2(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}P_\nu (z).`$ (A.18)
(The choice of 1 as the upper limit in (A.17) is not important; any value greater than 0 will be equally good.) Analysis of the terms in $`P_\nu (z)`$ responsible for divergence when $`\mathrm{}(s)<\frac{1}{2}`$ leads us to define

$$P_\nu (z)=\frac{1}{\pi }a^{1/2}ze^{(1a)z}\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)g_\nu (z),$$
(A.19)
where

$`g_\nu (z)`$ $`=`$ $`1e^{2(1a)z}{\displaystyle \frac{\mathrm{\Sigma }_\nu ^{(K)}(z)\mathrm{\Sigma }_\nu ^{(I)}(az)}{\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)}}`$ (A.20)
$`=`$ $`1{\displaystyle \frac{k_\nu (z)i_\nu (az)}{i_\nu (z)k_\nu (az)}}.`$ (A.21)
Substitution of (A.19) into (A.17) results in

$`\omega _1(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}\{{\displaystyle \frac{1}{2s}}+{\displaystyle \frac{1a}{12s}}\}+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}dzz^{2s}\{{\displaystyle \frac{d}{dz}}\mathrm{ln}g_\nu (z)`$ (A.22)
$`+{\displaystyle \frac{d}{dz}}\mathrm{ln}\left[\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)\right]\}`$
where the result has been analytically continued to $`\mathrm{}(s)<0`$. Substitution of (A.19) into (A.18), assuming initially $`\frac{1}{2}<\mathrm{}(s)<1`$, and then analytically continuing to $`\mathrm{}(s)<0`$, results in

$`\omega (s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}g_\nu (z)`$ (A.23)
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\left[\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)\right]`$
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\left[\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)\right].`$
The first two terms are analytic for $`\mathrm{}(s)<0`$; however, the last term has poles at various values of $`s`$ which we need to determine if we are to perform the required analytic continuation.

We have constructed $`\mathrm{\Sigma }_\nu ^{(I)}(z)`$ so that $`\mathrm{\Sigma }_\nu ^{(I)}(z)1+𝒪(z^1)`$ as $`z\mathrm{}`$. This means that we can write

$$\mathrm{ln}\mathrm{\Sigma }_\nu ^{(I)}(z)\underset{k=1}{\overset{\mathrm{}}{}}d_kz^k$$
(A.24)
for large $`z`$, with some coefficients $`d_k`$. The coefficients $`d_k`$ will determine the poles of $`\omega (s)`$. Because of (A.15) we have

$$\mathrm{ln}\left[\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)\right]\underset{k=1}{\overset{\mathrm{}}{}}\sigma _kz^k$$
(A.25)
with

$$\sigma _k=\left[1+\frac{(1)^k}{a^k}\right]d_k.$$
(A.26)
Using (A.13) and (A.14) in (A.24) it is easy to calculate the coefficients $`d_k`$. The first four coefficients are found to be

$`d_1`$ $`=`$ $`\alpha {\displaystyle \frac{3}{8}}{\displaystyle \frac{1}{2}}\nu ^2,`$ (A.27)
$`d_2`$ $`=`$ $`{\displaystyle \frac{1}{4}}\nu ^2{\displaystyle \frac{3}{16}}+{\displaystyle \frac{1}{2}}\alpha {\displaystyle \frac{1}{2}}\alpha ^2,`$ (A.28)
$`d_3`$ $`=`$ $`{\displaystyle \frac{1}{2}}\alpha \nu ^2{\displaystyle \frac{1}{2}}\alpha ^2+{\displaystyle \frac{3}{8}}\alpha +{\displaystyle \frac{23}{48}}\nu ^2+{\displaystyle \frac{1}{24}}\nu ^4{\displaystyle \frac{21}{128}}+{\displaystyle \frac{1}{3}}\alpha ^3,`$ (A.29)
$`d_4`$ $`=`$ $`\alpha \nu ^2{\displaystyle \frac{1}{2}}\alpha ^2+{\displaystyle \frac{3}{8}}\alpha +{\displaystyle \frac{13}{16}}\nu ^2{\displaystyle \frac{1}{8}}\nu ^4`$ (A.30)
$`+{\displaystyle \frac{1}{2}}\alpha ^3+{\displaystyle \frac{1}{2}}\alpha ^2\nu ^2{\displaystyle \frac{1}{4}}\alpha ^4{\displaystyle \frac{27}{128}},`$
where

$$\alpha =\frac{D}{2}(14\xi ).$$
(A.31)
(For $`\alpha =2`$, which is the case for $`D=4`$ and $`\xi =0`$, the results reduce to those given in .) By adding and subtracting terms to the last integrand of (A.23) we may analytically continue the result to $`\mathrm{}(s)<0`$. The more negative $`\mathrm{}(s)`$ becomes, the more terms we need to subtract. In this way we obtain

$`\omega (s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}g_\nu (z)`$ (A.32)
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\left[\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)\right]`$
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\left\{\mathrm{ln}\left[\mathrm{\Sigma }_\nu ^{(I)}(z)\mathrm{\Sigma }_\nu ^{(K)}(az)\right]{\displaystyle \underset{k=1}{\overset{N}{}}}\sigma _kz^k\right\}`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{k=1}{\overset{N}{}}}{\displaystyle \frac{k\sigma _k}{(k+2s)}}.`$
Provided that $`N`$ is chosen so that $`\mathrm{}(N+1+2s)>0`$ the analytic continuation to any value of $`s`$ with $`\mathrm{}(s)<0`$ may be performed. All of the possible poles can be found in the last term of (A.32). Some procedure like the one we have described is essential to provide the proper analytic continuation of $`\omega (s)`$ to $`\mathrm{}(s)<0`$.

We are interested in the value of $`\omega (s)`$ at $`s=pD/2`$ for $`p=0,1`$. There is a simple pole at this value of $`s`$ coming from the $`k=D2p`$ term in the sum in the last term of (A.32). We can write

$$\omega (s)=\omega _{\mathrm{pole}}(s)+\omega _{\mathrm{reg}}(s)$$
(A.33)
where

$$\omega _{\mathrm{pole}}(s)=\frac{1}{\mathrm{\Gamma }(1s)}\frac{(D2p)\sigma _{D2p}}{(D2p+2s)}$$
(A.34)
contains the pole term (as well as a finite term found by expanding about the pole), and $`\omega _{\mathrm{reg}}(s)`$ is regular (ie analytic) at $`s=pD/2`$ and is given by an expression like (A.32) but with the term $`k=D2p`$ omitted from the sum in the last term. We are free to simply set $`s=pD/2`$ in $`\omega _{\mathrm{reg}}(s)`$, where $`D`$ is the physical dimension (ie 4 for the Randall-Sundrum model).

The most important thing about the regular part $`\omega _{\mathrm{reg}}(s)`$ is the dependence on $`a`$. By splitting up the logarithms in (A.32) it can be seen that

$$\omega _{\mathrm{reg}}(s)=C(s)+A_1(s)+A_2(s)$$
(A.35)
where

$`C(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\mathrm{\Sigma }_\nu ^{(I)}(z)`$ (A.36)
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\left\{\mathrm{ln}\mathrm{\Sigma }_\nu ^{(I)}(z){\displaystyle \underset{k=1}{\overset{N}{}}}d_kz^k\right\}`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{k=1}{\overset{N}{}}^{}}{\displaystyle \frac{kd_k}{(k+2s)}}`$
is independent of $`a`$, and

$`A_1(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}g_\nu (z)`$ (A.37)
$`A_2(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\mathrm{\Sigma }_\nu ^{(K)}(az)`$ (A.38)
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\left\{\mathrm{ln}\mathrm{\Sigma }_\nu ^{(K)}(az){\displaystyle \underset{k=1}{\overset{N}{}}}{\displaystyle \frac{d_k(1)^k}{(az)^k}}\right\}`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{k=1}{\overset{N}{}}^{}}{\displaystyle \frac{k(1)^kd_k}{(k+2s)a^k}}`$
contain the $`a`$-dependence. The prime on the sums in (A.36) and (A.38) denote that the term $`k=D2p`$, which gives rise to the pole, is omitted.

The $`a`$-dependence of (A.38) can be made more explicit by a change of variable $`za^1z`$ in both integrals :

$`A_2(s)`$ $`=`$ $`{\displaystyle \frac{a^{2s}}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{a}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\mathrm{\Sigma }_\nu ^{(K)}(z)`$ (A.39)
$`+{\displaystyle \frac{a^{2s}}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{a}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\left\{\mathrm{ln}\mathrm{\Sigma }_\nu ^{(K)}(z){\displaystyle \underset{k=1}{\overset{N}{}}}{\displaystyle \frac{d_k(1)^k}{z^k}}\right\}`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{k=1}{\overset{N}{}}^{}}{\displaystyle \frac{k(1)^kd_k}{(k+2s)a^k}}.`$
Now split up the integration range of the second integral as $`{\displaystyle \underset{a}{\overset{\mathrm{}}{}}}={\displaystyle \underset{a}{\overset{1}{}}}+{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}`$. After a short calculation we are left with

$`A_2(s)`$ $`=`$ $`a^{2s}A(s){\displaystyle \frac{a^{2s}}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{a}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}{\displaystyle \underset{k=1}{\overset{N}{}}}{\displaystyle \frac{(1)^kd_k}{z^k}}`$ (A.40)
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{k=1}{\overset{N}{}}^{}}{\displaystyle \frac{k(1)^kd_k}{(k+2s)a^k}}`$
where

$`A(s)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{0}{\overset{1}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\mathrm{ln}\mathrm{\Sigma }_\nu ^{(K)}(z)`$ (A.41)
$`+{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{1}{\overset{\mathrm{}}{}}}𝑑zz^{2s}{\displaystyle \frac{d}{dz}}\left\{\mathrm{ln}\mathrm{\Sigma }_\nu ^{(K)}(z){\displaystyle \underset{k=1}{\overset{N}{}}}{\displaystyle \frac{d_k(1)^k}{z^k}}\right\}`$
is independent of $`a`$ (and is analytic for $`\mathrm{}(s)<0`$). Finally we may perform the integration in (A.40) to find that

$$A_2(s)=a^{2s}\overline{A}(s)\frac{(D2p)(1)^D}{\mathrm{\Gamma }(1+D/2p)}d_{D2p}a^{2pD}\mathrm{ln}a$$
(A.42)
with

$$\overline{A}(s)=A(s)\frac{1}{\mathrm{\Gamma }(1s)}\underset{k=1}{\overset{N}{}}^{}\frac{k(1)^kd_k}{(k+2s)}.$$
(A.43)
We can conclude that for $`s=pD/2`$

$`\omega _{\mathrm{reg}}(pD/2)`$ $`=`$ $`C(pD/2)+A_1(pD/2)+a^{2pD}\overline{A}(pD/2)`$ (A.44)
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(1s)}}{\displaystyle \underset{k=1}{\overset{N}{}}^{}}{\displaystyle \frac{k(1)^kd_k}{(k+2s)}}`$
where $`C(s),\overline{A}(s)`$ are independent of $`a`$.

Now turn to a similar analysis for the simpler pole part of $`\omega (s)`$ given in (A.34). Because there is a pole, we will use dimensional regularisation and set $`DD+ϵ`$, with $`s=p(D+ϵ)/2`$. Note that it is only in the $`s`$-dependent terms that we replace $`DD+ϵ`$, since the factor of $`(D2p)\sigma _{D2p}`$ came from a term in the asymptotic expansion which involved only integral powers. Because $`\sigma _kd_k`$, and $`d_k`$ depends on the order of the Bessel functions $`\nu `$, there can be a $`D`$-dependence hidden inside $`\sigma _{D2p}`$. However it is easy to see that this dependence can only affect terms which are constant, or proportional to $`a^{2pD}`$. Such terms can be absorbed into redefined functions $`C(pD/2)`$ and $`\overline{A}(pD/2)`$ respectively, which we show in Sec. 3 can be removed by finite renormalizations. There is no point in giving the explicit form of such terms, and we can safely ignore the $`ϵ`$-dependence of $`\sigma _{D2p}`$. It is possible to write

$$\omega _{\mathrm{pole}}(s)=\frac{\alpha _1}{ϵ}+\alpha _0+\alpha _1a^{2pD}$$
(A.45)
for

$$\alpha _1=\frac{(D2p)}{\mathrm{\Gamma }(1+D/2p)}\left[1+\frac{(1)^D}{a^{D2p}}\right]d_{D2p}$$
(A.46)
with $`\alpha _0`$ and $`\alpha _1`$ independent of $`a`$. (As explained above, it is a straightforward matter to obtain explicit results for $`\alpha _0`$ and $`\alpha _1`$ at $`s=pD/2`$ by expanding (A.34) about the pole.)

## Appendix B Evaluation of the sums for twisted fields

In this appendix we look at the calculation of $`\omega (s)`$ using the solution to (2.31). The details are very similar to those of the preceding appendix, so we will be brief. $`\omega (s)`$ is expressed as a contour integral as in (A.2) where $`\stackrel{~}{F}_\nu (z)`$ appears in place of $`F_\nu (z)`$. (A.3) still holds, but with $`P_\nu (z)`$ now given by

$$P_\nu (z)=\frac{2}{\pi }\left[I_\nu (z)K_\nu (az)I_\nu (az)K_\nu (z)\right].$$
(B.1)
Using (A.7) and (A.8) results in

$$P_\nu (z)=\frac{1}{\pi }a^{1/2}z^1e^{(1a)z}\mathrm{\Sigma }^{(I)}(z)\mathrm{\Sigma }^{(K)}(az)g_\nu (z)$$
(B.2)
where

$$g_\nu (z)=1\frac{K_\nu (z)I_\nu (az)}{K_\nu (az)I_\nu (z)}.$$
(B.3)
The difference between the twisted and untwisted cases lies in the asymptotic behaviour of the integrand of the contour integral for large $`|z|`$. (Compare (A.19) with (B.2).) In the untwisted case the leading term in the asymptotic expansion comes from the derivative of the Bessel functions in (A.5,A.6). These terms are absent from the twisted field expressions. Because of this different asymptotic behaviour, which lead to divergent integrals dealt with by adding and subtracting terms as discussed in the previous appendix, it is difficult to deal with both twisted and untwisted fields in one single calculation.

The steps in the analytic continuation of $`\omega (s)`$ follow much as before with slightly altered expressions. In place of (A.23) we have the simpler expansion

$$\mathrm{ln}\mathrm{\Sigma }^{(I)}(z)\underset{k=1}{\overset{\mathrm{}}{}}\stackrel{~}{d}_kz^k$$
(B.4)
where the first few coefficients are given by

$`\stackrel{~}{d}_1`$ $`=`$ $`{\displaystyle \frac{1}{8}}(14\nu ^2),`$ (B.5)
$`\stackrel{~}{d}_2`$ $`=`$ $`{\displaystyle \frac{1}{16}}(14\nu ^2),`$ (B.6)
$`\stackrel{~}{d}_3`$ $`=`$ $`{\displaystyle \frac{1}{384}}(25104\nu ^2+16\nu ^4),`$ (B.7)
$`\stackrel{~}{d}_4`$ $`=`$ $`{\displaystyle \frac{1}{128}}(1356\nu ^2+16\nu ^4),`$ (B.8)
$`\stackrel{~}{d}_5`$ $`=`$ $`{\displaystyle \frac{1}{5120}}(10734748\nu ^2+1840\nu ^464\nu ^6),`$ (B.9)
$`\stackrel{~}{d}_6`$ $`=`$ $`{\displaystyle \frac{1}{192}}(103465\nu ^2+216\nu ^416\nu ^6).`$ (B.10)
The remaining analysis of Appendix A follows as before with these new coefficients in place of the $`d_k`$, and (B.3) used for $`g_\nu (z)`$. (This changes the result we called $`A_1(s)`$ in (A.37).)


format: text

text:
# Photon Structure Functions from Quenched Lattice QCD

## Abstract

We calculate the first moment of the photon structure function, $`x^\gamma =_0^1𝑑xF_2^\gamma (x,Q^2)`$, on the quenched lattices with $`\beta =6.0`$ using the formalism developed by the authors recently. In this exploratory study, we take into account only the connected contractions. The result is compared with the experimental data as well as model predictions.

preprint: UMD PP#01-145preprint: DOE/ER/40762-229

The parton (quark and gluon) distributions in hadrons, which can be measured from lepton-hadron deep inelastic scattering and other hard processes, have played a crucial role in understanding high-energy scattering and the hadron structure. Recently it has become possible to calculate moments of these distributions from first principles using lattice QCD .

There are many other physical observables which are yet to be calculated from lattice QCD. A special class of of these involve either virtual or real photons. Until recently, the conditions under which these quantities can be calculated using lattice QCD remained unclear. The reason is that the photon is not an eigenstate of QCD. Rather, a “photon” state in nature is a superposition of the U(1) gauge boson and quark-gluon configurations which are suppressed by the electromagnetic coupling. The matrix elements of QCD operators in the photon states are time-dependent correlations which are defined in the Minkowski space, and the standard method used for calculating hadron matrix elements on the lattice is not directly applicable .

In our earlier paper , we showed that the matrix element of a quark-gluon operator between photon states can be evaluated using lattice QCD. The relevant expression is

$$\gamma (p\lambda ^{})|𝒪(0)|\gamma (p\lambda )=e^2d^4xd^4ye^{\omega (x_4y_4)}e^{i\stackrel{}{p}(\stackrel{}{x}\stackrel{}{y})}0|T_Eϵ^{}(\lambda ^{})J(x)𝒪(0)ϵ(\lambda )J(y)|0$$
(1)
where every quantity has been expressed in the Euclidean space: $`T_E`$ is the Euclidean time-ordering; the electromagnetic current $`J_\mu (\mu =1,2,3,4)`$ is defined as $`_fe_f\overline{\psi }_f\gamma _\mu \psi _f`$ with Euclidean $`\gamma `$ matrices $`\gamma _i=i\gamma _M^i(i=1,2,3)`$ and $`\gamma _4=\gamma ^0`$; the Euclidean photon polarization vector is defined as $`ϵ_i=ϵ_M^i`$ and $`ϵ_4=iϵ^0(ϵ^4=iϵ^0)`$; $`i\omega `$ is the Euclidean-space photon energy. \[In Ref. , we did not continue the electromagnetic current and the polarization vector to Euclidean space and superficially the result there differs by an overall sign.\] The above expression is the same as the naive analytic continuation of the matrix element from the Minkowski space. However, Eq. (1) is only valid when there exists an energy gap between the photon energy and the lowest hadronic state of the same quantum number. Otherwise, the integral becomes divergent.

One of the quantities we can study using the above formalism is the photon structure functions, which can be measured from collisions between real or virtual photons with highly virtual ones, achievable in $`e^+e^{}e^+e^{}+`$hadrons . While a photon is normally considered a structureless particle, it can fluctuate into a charged fermion-antifermion pair or more complicated hadronic states which can be revealed through interactions with a highly virtual photon. In fact, the photon structure functions can be defined from its hadronic tensor and the moments of the structure functions can be obtained from operator product expansion, much in the same way as for the hadron structure functions. There exists a substantial amount of experimental data for unpolarized structure functions already ; and future experiments from HERA and $`e^+e^{}`$ are expected to measure the polarized structure functions as well . (For a recent review on theoretical and experimental progress on this subject, refer to Ref. .)

Here we report on the first lattice study of the structure function $`F_2^\gamma `$. More specifically, we measure the first moment of $`F_2^\gamma `$ which corresponds to the average fractional momentum carried by the quarks (weighed with the fractional charge squared for each flavor) in a real photon. The measurement is done on $`\beta =6.0`$ lattices, corresponding to $`Q^2a^2(2.4\text{GeV})^2`$, in the quenched approximation for three degenerate flavors of quarks. Also, we choose to ignore the contributions from disconnected diagrams in which the quark lines do not connect any of the external vertices.

We begin by briefly reviewing the definition of the photon structure functions and their relation to the cross sections. We closely follow  for the definition of $`F_2^\gamma (x,Q^2)`$. For the deep inelastic $`e+\gamma `$ process in a frame where the incoming electron and the photon are collinear, we write

$`e^\pm (k)+\gamma (p)e^\pm (k^{})+X(p_X),`$ (2)
$`k^\mu =(E,0,0,E),`$ (3)
$`k_{}^{}{}_{}{}^{\mu }=(E^{},0,E^{}\mathrm{sin}\theta ,E^{}\mathrm{cos}\theta ),`$ (4)
$`p^\mu =(E_\gamma ,0,0,E_\gamma ).`$ (5)
The relevant kinetic variables are $`Q^2=(kk^{})^2`$, $`\nu =pq`$, $`q=kk^{}`$, $`x=Q^2/2(pq)`$, $`y=qp/kp`$. The photon structural information is contained in the following tensor,

$`W^{\mu \nu }`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle \frac{1}{2}}{\displaystyle \underset{X,\lambda }{}}\gamma (p\lambda )|J^\mu (0)|XX|J^\nu (0)|\gamma (p\lambda )(2\pi )^4\delta ^4(p_Xpq)`$ (6)
$`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle \frac{1}{2}}{\displaystyle \underset{\lambda }{}}{\displaystyle \gamma (p\lambda )|[J^\mu (x),J^\nu (0)]|\gamma (p\lambda )e^{iqx}d^4x}`$ (7)
$`=`$ $`W_1^\gamma (\nu ,Q^2)\left[g^{\mu \nu }+{\displaystyle \frac{q^\mu q^\nu }{q^2}}\right]+W_2^\gamma (\nu ,Q^2)\left[p^\mu {\displaystyle \frac{pq}{q^2}}q^\mu \right]\left[p^\nu {\displaystyle \frac{pq}{q^2}}q^\nu \right].`$ (8)
Here $`J_\mu (x)=_fe_f\overline{\psi }_f(x)\gamma _\mu \psi _f(x)`$ is the electromagnetic current and $`e_f`$ is the fractional charge of the quark of flavor $`f`$. The covariant normalization has been used for the photon state, $`\gamma (p\lambda )|\gamma (p^{}\lambda ^{})=(2\pi )^32\omega _\stackrel{}{p}\delta (\stackrel{}{p}\stackrel{}{p^{}})\delta _{\lambda ,\lambda ^{}}`$.

Since we are interested in the unpolarized photon properties here, we have averaged over the photon polarization $`\lambda `$ in the above equations. \[The polarized photon structure functions are interesting by themselves and can be studied in a similar approach .\]

Using these variables, the differential cross section is given by

$$\frac{d\sigma }{dE^{}d\mathrm{cos}\theta }=\frac{8\pi \alpha _{\mathrm{em}}^2E_{}^{}{}_{}{}^{2}}{Q^4E_\gamma }\left[2E_\gamma ^2W_2^\gamma (\nu ,Q^2)\mathrm{cos}^2\left(\frac{\theta }{2}\right)+W_1^\gamma (\nu ,Q^2)\mathrm{sin}^2\left(\frac{\theta }{2}\right)\right],$$
(9)
or, in terms of $`x`$ and $`y`$ and scaling functions,

$$\frac{d\sigma }{dxdy}=\frac{16\pi \alpha _{\mathrm{em}}^2EE_\gamma }{Q^4}\left[(1y)F_2^\gamma (x,Q^2)+xy^2F_1^\gamma (x,Q^2)\right],$$
(10)
where the dimensionless scaling functions are defined as $`F_2^\gamma (x,Q^2)=\nu W_2^\gamma (\nu ,Q^2)`$ and $`F_1^\gamma (x,Q^2)=W_1^\gamma (\nu ,Q^2)`$.

In the parton model, the photon scaling function $`F_2^\gamma (x)`$ is related to the quark distributions $`q_f^\gamma (x)`$ in the photon—the probability of finding a quark with momentum fraction $`x`$,

$$F_2^\gamma (x)=x\underset{f}{}e_f^2(q_f^\gamma (x)+\overline{q}_f^\gamma (x)).$$

Taking into account QCD radiative corrections, the moments of the photon structure functions can be written in the following factorized form,

$$_0^1x^{n2}F_2^\gamma (x,Q^2)dx=\underset{f}{}C_{2n,f}(\frac{Q^2}{\mu ^2},g(\mu ))x^{n1}_f^\gamma ,(n=2,4,\mathrm{},)$$

where $`C_{2n,f}`$ are the perturbative coefficient functions and $`x^{n1}_f^\gamma `$ are the nonperturbative moments of the parton (quark and gluon) distributions ($`f`$ here refers to gluons as well). For quarks, the moments can be calculated as the matrix elements of local operators,

$$\frac{1}{2}\underset{\lambda }{}\gamma (p\lambda )|𝒪_{\{\mu _1\mathrm{}\mu _n\},f}|\gamma (p\lambda )=2x^{n1}_f^\gamma [p_{\mu _1}\mathrm{}p_{\mu _n}\mathrm{traces}],$$
(11)
$$𝒪_{\mu _1\mathrm{}\mu _n,f}=\left(\frac{i}{2}\right)^{n1}\overline{\psi }_f\gamma _{\mu _1}\underset{\mu _2}{\overset{}{D}}\mathrm{}\underset{\mu _n}{\overset{}{D}}\psi _f,$$
(12)
where $`\underset{\mu }{\overset{}{D}}=\underset{\mu }{\overset{}{D}}\underset{\mu }{\overset{}{D}}`$ is the covariant derivative and $`\{\mathrm{}\}`$ denotes the symmetrization of indices.

In the following discussion, we ignore the radiative corrections to the coefficient functions and take $`C_{2n,f}=e_f^2`$ for quarks and 0 for gluons. We then define $`x^\gamma =_0^1F_2^\gamma (x,Q^2)𝑑x`$, which can be calculated as the matrix elements of the following operator;

$$𝒪_{\mu \nu }^{(\mathrm{cont})}(x)=\underset{f}{}e_f^2\frac{i}{2}\overline{\psi }_f(x)\gamma _{\{\mu }\underset{\nu \}}{\overset{}{D}}\psi _f(x).$$
(13)
In lattice QCD, we need to define the corresponding lattice operator in Euclidean space:

$$𝒪_{\mu \nu }^{(\mathrm{latt})}(x)=\underset{f,y}{}e_f^2\left[\overline{\psi }_f(x)\mathrm{\Gamma }_{\mu \nu }^𝒪(x,y)\psi _f(y)+\overline{\psi }_f(y)\mathrm{\Gamma }_{\mu \nu }^𝒪(y,x)\psi _f(x)\right],$$
(14)
$$\mathrm{\Gamma }_{\mu \nu }^𝒪(x,y)=\frac{i}{8}[\gamma _\mu (U_\nu (x)\delta _{y,x+\widehat{\nu }}U_\nu ^{}(y)\delta _{y,x\widehat{\nu }})+(\mu \nu )],$$
(15)
where the lattice spacing($`a`$) is set to 1 and the lattice covariant derivatives are defined as

$$\underset{\mu }{\overset{}{D}}\psi (x)=\frac{1}{2}(U_\mu (x)\psi (x+\widehat{\mu })U_\mu ^{}(x\widehat{\mu })\psi (x\widehat{\mu })),$$

$$\overline{\psi }(x)\underset{\mu }{\overset{}{D}}=\frac{1}{2}(\overline{\psi }(x+\widehat{\mu })U_\mu ^{}(x)\overline{\psi }(x\widehat{\mu })U_\mu (x\widehat{\mu })).$$

The corresponding Euclidean matrix element to Eq. (12) is

$$\frac{1}{2}\underset{\lambda }{}\gamma (p\lambda )|𝒪_{\mu \nu }^{(\mathrm{latt})}|\gamma (p\lambda )=2ix^\gamma p_{\mu E}p_{\nu E},$$
(16)
where the photon Euclidean four-momentum is defined as $`p_{\mu E}=(\stackrel{}{p},i\omega )`$.

The photon matrix element in lattice QCD is

$$\begin{array}{c}\gamma (p\lambda ^{})|𝒪_{\mu \nu }(0)|\gamma (p\lambda )^{(\mathrm{latt})}=e^2ϵ_\alpha ^{}(\lambda ^{})ϵ_\beta (\lambda )\hfill \\ \hfill \times \underset{x,y}{}e^{\omega (x_4y_4)}e^{i\stackrel{}{k}(\stackrel{}{x}\stackrel{}{y})}0|T_EJ_\alpha (x)𝒪_{\mu \nu }^{(\mathrm{latt})}(0)J_\beta (y)|0.\end{array}$$
(17)
The Euclidean Green’s function can be calculated as follows,

$`0|T_EJ_\alpha (x)𝒪_{\mu \nu }^{(\mathrm{latt})}(0)J_\beta (y)|0={\displaystyle \underset{z,f,f^{},f^{\prime \prime }}{}}`$ $`e_fe_f^{}e_{f^{\prime \prime }}^2`$ (18)
$`\times 0|T_E\overline{\psi }_f(x)\gamma _\alpha \psi _f(x)[\overline{\psi }_{f^{\prime \prime }}(0)\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)\psi _{f^{\prime \prime }}(z)`$ $`+\overline{\psi }_{f^{\prime \prime }}(z)\mathrm{\Gamma }_{\mu \nu }^𝒪(z,0)\psi _{f^{\prime \prime }}(0)]\overline{\psi }_f^{}(y)\gamma _\beta \psi _f^{}(y)|0`$ (19)
$`={\displaystyle \underset{z,f}{}}e_f^4(2\kappa )^3[\mathrm{Tr}[\gamma _\alpha S(x,0)\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)S(z`$ $`,y)\gamma _\beta S(y,x)]_g`$ (20)
$`+\mathrm{Tr}[\gamma _\beta S(y,0)\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)S(z,x)\gamma _\alpha S(x,y)]_g`$ $`+\mathrm{Tr}[\gamma _\alpha S(x,z)\mathrm{\Gamma }_{\mu \nu }^𝒪(z,0)S(0,y)\gamma _\beta S(y,x)]_g`$ (21)
$`+\mathrm{Tr}`$ $`[\gamma _\beta S(y,z)\mathrm{\Gamma }_{\mu \nu }^𝒪(z,0)S(0,x)\gamma _\alpha S(x,y)]_g]+\mathrm{},`$ (22)
where we have shown explicitly the connected contraction with one trace and the ellipsis represents disconnected diagrams with more than one trace, and $`\mathrm{}_g`$ denotes averaging over gauge configurations. The propagator $`S(x,y)`$ is solved from $`_xD(y,x)S(x,z)=\delta _{y,z}`$ where

$$D(x,y)=\delta _{x,y}\kappa \underset{\mu }{}[(1\gamma ^\mu )U_\mu (x)\delta _{y,x+\widehat{\mu }}+(1+\gamma ^\mu )U_\mu ^{}(y)\delta _{x,y+\widehat{\mu }}],$$
(23)
is the scaled Wilson-Dirac operator. The Dirac and color indices are suppressed for readability and the traces are over those indices.

In this exploratory study, we will ignore the disconnected contribution. For the connected contraction, we can show, using the charge conjugation property of the Wilson-Dirac operator, that the average traces for the left and the right derivatives are the same and, hence, we only need to evaluate one of them. The Wilson-Dirac operator $`D(x,y)`$ is invariant under the charge conjugation:

$$\psi (x)C\overline{\psi }(x)^T,\overline{\psi }(x)\psi (x)^TC^1,U_\mu (x)(U_\mu (x))^{},$$
(24)
where the matrix $`C`$ satisfies $`C^1\gamma _\mu C=\gamma _\mu ^T`$ and $`T`$ denotes transpose operation. Using Eq. (24), we can show

$$S(x,y|U)^T=C^1S(y,x|U^{})C,$$
(25)
where we explicitly denote the dependence of the propagator on the gauge field.

Now we can relate the traces which arise from the left derivative of $`\mathrm{\Gamma }_{\mu \nu }^𝒪`$ with those from the right derivative. For example,

$`\mathrm{Tr}\left[\gamma _\alpha S(x,0|U)\gamma _\nu U_\mu (0)S(\widehat{\mu },y|U)\gamma _\beta S(y,x|U)\right]`$ (26)
$`=\mathrm{Tr}\left[S(y,x|U)^T\gamma _\beta ^TS(\widehat{\mu },y|U)^TU_\mu ^T(0)\gamma _\nu ^TS(x,0|U)^T\gamma _\alpha ^T\right]`$ (27)
$`=\mathrm{Tr}\left[S(x,y|U^{})\gamma _\beta S(y,\widehat{\mu }|U^{})(U^{})_\mu ^{}(0)\gamma _\nu S(0,x|U^{})\gamma _\alpha \right].`$ (28)
Using the previous equation and averaging over gauge configurations, we get

$$\mathrm{Tr}[\gamma _\alpha S(x,0)\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)S(z,y)\gamma _\beta S(y,x)]_g=\mathrm{Tr}[\gamma _\beta S(y,z)\mathrm{\Gamma }_{\mu \nu }^𝒪(z,0)S(0,x)\gamma _\alpha S(x,y)]_g,$$
(29)
and similarly for $`xy`$ and $`\alpha \beta `$. Finally, we arrive at the following expression used for actual calculation:

$$\begin{array}{c}\gamma (p\lambda ^{})|𝒪_{\mu \nu }(0)|\gamma (p\lambda )_{\mathrm{conn}.}^{(\mathrm{latt})}=e^2ϵ_\alpha ^{}(\lambda ^{})ϵ_\beta (\lambda )\underset{x,y}{}e^{\omega (x_4y_4)}e^{i\stackrel{}{p}(\stackrel{}{x}\stackrel{}{y})}\underset{z,f}{}e_f^4(2\kappa )^3\times \hfill \\ \hfill 2\mathrm{Tr}[\gamma _\alpha S(x,0)\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)S(z,y)\gamma _\beta S(y,x)]+\mathrm{Tr}[\gamma _\beta S(y,0)\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)S(z,x)\gamma _\alpha S(x,y)]_g.\end{array}$$
(30)
The above matrix element is evaluated by summing over the propagators starting from point sources placed at the position of $`𝒪`$ instead of one of the electromagnetic vertices. The reason for this strategy, different from hadronic structure functions , is that, as evident from Eq. (20), the electromagnetic currents at both source and sink should be summed over all the time slices with proper momentum dependence, in contrast to the hadron structure functions where the time coordinates of the hadron operators are fixed. \[Because of the translational invariance in the time direction, one can still fix the time coordinate of one of the electromagnetic vertices instead of the operator observable. The answer shall not be much different if the time direction is sufficiently large.\]

In practical simulations, two sequential propagators are generated from each of the $`12=(4\times 3)`$ point sources

$`M_\alpha (\stackrel{}{p},\omega ,y)={\displaystyle \underset{x}{}}S(y,x)\gamma _\alpha S(x,0)e^{i\stackrel{}{p}\stackrel{}{x}+\omega x_4},`$ (31)
$`M_\beta ^{}(\stackrel{}{p},\omega ,x)={\displaystyle \underset{y}{}}S(x,y)\gamma _\beta S(y,0)e^{+i\stackrel{}{p}\stackrel{}{y}\omega y_4}.`$ (32)
The complementary propagator generated from the operator $`𝒪`$ is

$`N_{\mu \nu }^𝒪(x)={\displaystyle \underset{z}{}}\mathrm{\Gamma }_{\mu \nu }^𝒪(0,z)S(z,x),`$which can be calculated from the point sources as $`N_{\mu \nu }^𝒪(x)=_z\gamma _5S(x,z)\gamma _5\mathrm{\Gamma }_{\mu \nu }^𝒪(z,0)`$. In terms of the above propagators,

$$\begin{array}{c}x_{\mathrm{conn}.}^\gamma /\alpha _{\mathrm{em}}=\frac{4\pi }{2p_\alpha p_\beta }\underset{f}{}e_f^4ϵ_\mu ^{}ϵ_\nu (2\kappa )^3\hfill \\ \hfill \times 2\left[\underset{y}{}\mathrm{Tr}[N_{\mu \nu }^𝒪(y)\gamma _\beta M_\alpha (y)]e^{i\stackrel{}{p}\stackrel{}{y}\omega y_4}+\underset{x}{}\mathrm{Tr}[N_{\mu \nu }^𝒪(x)\gamma _\alpha M_\beta ^{}(x)]e^{i\stackrel{}{p}\stackrel{}{x}+\omega x_4}\right],\end{array}$$
(33)
where $`\alpha `$ and $`\beta `$ are fixed indices (without summation) and the polarization vector can be chosen as a linear one in the $`x`$ direction, $`ϵ_\mu =(1,0,0,0)`$. For 3 degenerate quarks, $`_fe_f^4`$ = $`(1/3)^4+(2/3)^4+(1/3)^4`$ = 2/9. \[Alternatively, one can generate doubly sequential propagators:

$`M_{(2)\alpha \beta }(z)={\displaystyle \underset{x,y}{}}S(z,y)\gamma _\beta S(y,x)\gamma _\alpha S(x,0)e^{i\stackrel{}{p}\stackrel{}{y}\omega y_4}e^{i\stackrel{}{p}\stackrel{}{x}+\omega x_4},`$ (34)
$`M_{(2)\alpha \beta }^{}(z)={\displaystyle \underset{x,y}{}}S(z,x)\gamma _\alpha S(x,y)\gamma _\beta S(y,0)e^{i\stackrel{}{p}\stackrel{}{y}\omega y_4}e^{i\stackrel{}{p}\stackrel{}{x}+\omega x_4},`$ (35)
which would eliminate the need for separate inversions for different operators. However, the successive inversions are then not desirable numerically when high precision is needed.\]

We use $`\beta =6.0,16^3\times 32`$ quenched configurations generated by Kilcup et al. , available from NERSC lattice archive. The lattice spacing is determined by taking the $`\rho `$ meson mass at the chiral limit: $`a^1`$ 2.4 GeV . To decide on the boundary condition for the Wilson fermion in time direction, we made a test run with $`\beta =5.7,8^2\times 16\times 32`$ lattices. For the antiperiodic boundary condition, the finite size effect gives a noisy $`x^\gamma `$ at a $`m_q`$ 200 MeV $`(\kappa =0.16)`$. On the other hand, the effect is under much better control for a $`m_q100`$ MeV $`(\kappa =0.165)`$ for the Dirichlet boundary condition. For heavier quark masses, both boundary conditions give consistent results. Therefore, we use a Dirichlet boundary condition to obtain the main result of the paper.

The photon momentum fraction $`x^\gamma `$ is measured by evaluating $`\gamma (p\lambda )|𝒪_{34}(0)|\gamma (p\lambda )`$ for the momentum $`\stackrel{}{k}=(0,0,2\pi /16),\omega =2\pi /16`$. We use hopping parameters $`\kappa =`$ 153, 0.154, 0.155, which corresponds to pion masses $`m_\pi `$ 0.424, 0.364, 0.301 in lattice units, respectively . We also evaluate $`x^\gamma `$ for $`\stackrel{}{k}=(0,0,\pi /16),\omega =\pi /16`$ ($`\kappa =`$ 154, 0.155), on the same Monte Carlo lattices, which is possible by employing the antiperiodic boundary condition for the pseudofermion fields in the $`z`$ direction.

To be close to the real world, the moment is calculated for three flavors of quarks with the degenerate masses. In the quenched approximation, the “$`\rho `$ meson” is the lowest hadronic state with photon quantum number. With these parameters, $`m_\rho `$ = 0.426 (in lattice units) for the lightest quark mass ($`\kappa =0.155`$) and $`\sqrt{m_\rho ^2+\stackrel{}{k}^2}`$ = 0.577. This gives $`L_t\times \left(\sqrt{m_\rho ^2+\stackrel{}{k}^2}\omega \right)5`$.

To test if this is large enough, we measured $`x^\gamma `$ with 3 different time positions of the operator $`𝒪`$. Figure 1 shows the value of the first moment of $`F_2^\gamma `$ for $`\kappa =0.153,0.154`$ with $`T_𝒪`$ = 16,17,18 (time coordinate runs from 0 to 31). The numerical values of $`x^\gamma `$ for different $`T_𝒪`$’s are consistent with each other (Fig. 1), which shows the temporal size of the box is large enough for the mass range. Note that the lattice result has been converted to a continuum renormalization scheme ($`\overline{MS}`$ in this case) by

$$𝒪_{\mu \nu }^{\overline{\mathrm{MS}}}(Q^2)=Z_𝒪𝒪_{\mu \nu }^{\mathrm{LATT}}(a^2)$$
(36)
$$Z_𝒪=1+\frac{g_0^2}{16\pi ^2}\frac{N_c^21}{2N_c}\left(\gamma ^{\overline{\mathrm{MS}}}\mathrm{log}(Q^2a^2)(B^{\mathrm{LATT}}B^{\overline{\mathrm{MS}}})\right)$$
(37)
Here we use the renormalization constant calculated perturbatively for $`Q=a^1`$ , $`Z_𝒪=0.9892`$. (The mixing between the operator $`𝒪_{\mu \nu }`$ (Eq. (13)) and gluon operator $`\mathrm{Tr}[F_{\mu \rho }F_{\rho \nu }]`$ is absent in the quenched approximation.)

Figure 2 shows the first moment of $`F_2^\gamma `$ for 3 degenerate quark flavors. The values of $`x^\gamma /\alpha _{\mathrm{em}}`$ for both momenta ($`\omega =2\pi /16,\pi /16`$) are well within the error bar. A linear fit to the chiral limit ($`m_\pi =0`$) gives $``$ 0.72(8), which is significantly larger than existing experimental results . Also, existing theoretical models such as the GRV model  and the Quark Parton Model (QPM)  predicts $`x^\gamma /\alpha _{\mathrm{em}}`$ $``$ (0.3$``$0.4) for the value of $`Q^2`$ investigated here.

We notice that the previous lattice studies of hadron structure functions  showed a consistent overestimation of the first moment of the structure functions ($`x`$). It has been suspected that the absence of sea quarks enhances the valence quark contribution, resulting in larger values of $`x`$. However, recent results from unquenched structure function studies  show no significant deviation from quenched results, which suggests that much smaller quark masses and/or a larger lattice size, quenched or unquenched, may be necessary to approach the continuum value. So some amount of overestimation for $`x^\gamma `$ is not unexpected. Still, the difference is more significant than those for hadron structure functions. This may suggest that $`F_2^\gamma `$ is actually larger at the extremes of $`x`$, where the experimental data is not available.

Another source of discrepancy is that the contribution from disconnected diagrams may be large. We are currently investigating this. \[It should be noted that for 3 degenerate (up, down and strange) quarks, the trace of electromagnetic operator, $`\mathrm{Tr}[J_\beta (x)]`$, is zero since $`_fe_f=0`$; we would only need to evaluate the diagrams similar to disconnected insertion (DI) diagrams studied for quantities, such as the strangeness magnetic moment of the nucleon .\]

In summary, the formalism necessary to compute the moments of the photon structure function is presented. We then compute the first moment of the unpolarized photon structure function $`F_2^\gamma `$ on quenched $`\beta =6.0`$ ($`Qa^1`$ 2.4 GeV) lattice configurations. The result is somewhat larger than the existing theoretical models and experimental results. The discrepancy is likely due to the quenched approximation and the disconnected contribution. If the result persists, it may indicate $`F_2^\gamma (x)`$ is significantly higher at near $`x=0`$ or $`x=1`$ than the theoretical models show. However, the systematic errors present in the measurement need to be understood to make more meaningful comparison.

This work is supported by funds provided by the U.S. Department of Energy (DOE) under grant No. DOE-FG02-94ER-40762. The numerical calculation reported here was performed on the Calico Alpha Linux Cluster at the Jefferson Laboratory, Virginia.