text:
# I Introducción

## I Introducción

El método canónico de cuantización permite establecer una correspondencia entre las cantidades clásicas observables y los operadores cuánticos. El hecho de que al emplearlo para cuantizar campos aparezcan energías de vacío divergentes es una consecuencia de la ambigüedad en la elección del orden usado para escribir los términos del hamiltoniano, que son productos de operadores que no conmutan entre sí. Los hamiltonianos escritos con diferentes órdenes difieren en términos proporcionales al operador identidad . Restar estos términos para pasar de un ordenamiento a otro es equivalente a trasladar el nivel cero de energía. Ya que solamente las diferencias de energías entre estados son observables, la elección del cero de energía no afecta la dinámica de los campos. Así, en la teoría de campos libres, la ambigüedad se resuelve con la prescripción de orden normal que conduce a un valor medio de vacío nulo para el operador hamiltoniano.

En 1948, H. B. G. Casimir mostró que la presencia de dos placas conductoras neutras conduce a un valor medio de vacío no nulo para el operador hamiltoniano de un campo electromagnético.

Por extensión, se conoce como efecto Casimir a la variación que sufren los valores medios de vacío de observables de un campo al imponérsele condiciones de contorno externas o al interactuar con un campo externo generado por determinada configuración de fuentes. La diferencia entre las energías de vacío de estas dos configuraciones se define como Energía de Casimir. Así, se tiene

$$E_{Cas}=E_0E_0^l.$$
(1)
donde $`E_0^l`$ es la energía de vacío del campo libre y $`E_0`$ es la energía de vacío del campo distorsionada por condiciones de contorno o la interacción con un campo externo.

Los métodos de evaluación de energías de vacío pueden reunirse en dos categorías : el método de suma directa de modos, aplicable toda vez que el espectro clásico sea explícitamente determinable, y el método local, que requiere el conocimiento de la función de Green del problema. En este trabajo, usaremos el método de suma de modos. Según éste, la energía de vacío se obtiene de sumar los autovalores de energía de los infinitos modos del problema particular. Para el caso de un campo de Dirac como el que trataremos, la energía de vacío resulta

$$E_0=\frac{1}{2}\underset{k}{}g_k|E_k|,$$
(2)
donde $`E_k`$ son los autovalores del espectro clásico de energías, correspondientes a autofunciones de cuadrado integrable, y $`g_k`$ la degeneración de cada uno de ellos.

La suma en (2) resulta, en general, divergente. A fin de dar una interpretación a las divergencias y poder calcular (1) debe recurrirse a alguna técnica de regularización. En este trabajo, usaremos el método de regularización de la función zeta .

La función zeta de un operador $`\widehat{A}`$ es definida, cuando son conocidos sus autovalores $`\lambda _n`$, como

$$\zeta _A=\underset{n}{}\lambda _n^s.$$
(3)
Esta suma resulta convergente sólo para $`Re(s)`$ mayor que cierta cantidad, que depende de la dimensión del espacio y del orden del operador . Sin embargo, la función puede extenderse analíticamente a todo el plano complejo, como una función meromorfa con polos en un conjunto discreto de valores de $`s`$.

En nuestro caso, definiremos la suma (2) como

$$E_0=\left(\frac{\mu ^{s+1}}{2}\underset{k}{}g_k|E_k|^s\right)_{s=1},$$
(4)
donde $`\mu `$ es un parámetro de escala con unidades de energía.

Es decir: con el objeto de evaluar las energías de vacío que aparecen en la definición (1) de la energía de Casimir, realizaremos la extensión analítica de la función zeta correspondiente al hamiltoniano clásico, evaluándola posteriormente en el valor $`s=1`$ de su extensión analítica. Al hacerlo, pueden ocurrir dos cosas: que el procedimiento arroje un resultado finito o bien que la extensión presente un polo en el valor de $`s`$ de interés. Cuando esto último ocurre, debe realizarse una renormalización de los parámetros en el hamiltoniano clásico, a fin de dar una interpretación física al resultado. En lo que sigue, veremos ejemplos de ambas situaciones.

El contenido del presente trabajo es el siguiente:

En la sección II, determinamos el operador hamiltoniano clásico para fermiones en un campo magnético uniforme en la dirección del eje $`z`$, de acuerdo al gauge elegido para el potencial electromagnético y la representación elegida de la matrices de Dirac.

En la sección III, obtenemos los autovalores del espectro clásico de energía, en caso en el que los fermiones están confinados en un plano perpendicular al campo magnético, hallando las autofunciones de cuadrado integrable del hamiltoniano clásico con impulso $`p_z=0`$.

En la sección IV, obtenemos las energías de vacío y de Casimir del problema en $`2+1`$ dimensiones, mediante la técnica de regularización de la función zeta.

En la sección V, determinamos el espectro clásico de energía cuando todo el espacio tridimensional es accesible para los fermiones. Lo hacemos aplicando una transformación de Lorentz sobre las autofunciones obtenidas en la sección III.

En la sección VI, obtenemos las energías de vacío y de Casimir del problema en $`3+1`$ dimensiones, por medio de una regularización zeta y una renormalización de la carga eléctrica.

En la sección VII, consideraremos la variación en la energía de vacío por efecto conjunto del campo magnético y condiciones de contorno en la dirección del campo. Aquí nuevamente relizamos una regularización zeta y una renormalización de la carga eléctrica.

## II Generalidades

Estudiaremos el problema de un campo de Dirac en interacción con un campo electromagnético, caracterizado por su potencial $`A_\mu (x)`$. La ecuación de Dirac en el espacio de Minkowski de $`3+1`$ dimensiones y en unidades naturales, $`\mathrm{}=c=1`$, se escribe entonces

$$(i\gamma ^\mu _\mu e\gamma ^\mu A_\mu m)\mathrm{\Psi }=0.$$
(5)
Trataremos el caso de un campo magnético uniforme $`\stackrel{}{B}`$ en la dirección del eje $`z`$. Sabemos que la elección de un potencial no es única y que el operador hamiltoniano depende de esta elección. Eligiendo el potencial electromagnético $`A_\mu =(0,\frac{B}{2}y,\frac{B}{2}x,0)`$ los operadores hamiltoniano y componente $`z`$ del momento angular total conmutan y, por lo tanto, tendremos autofunciones simultáneas de ambos operadores.

Utilizaremos la siguiente representación de las matrices de Dirac.

$$\gamma ^0=\left(\begin{array}{cc}\sigma _3& 0\\ 0& \sigma _3\end{array}\right)\gamma ^1=\left(\begin{array}{cc}i\sigma _2& 0\\ 0& i\sigma _2\end{array}\right)\gamma ^2=\left(\begin{array}{cc}i\sigma _1& 0\\ 0& i\sigma _1\end{array}\right)\gamma ^3=\left(\begin{array}{cc}0& i\sigma _1\\ i\sigma _1& 0\end{array}\right),$$
(6)
donde $`\sigma _i`$ son las matrices de Pauli.

A fin de determinar el espectro clásico de energías, factorizamos la dependencia temporal del espinor de Dirac y lo escribimos en términos de espinores de dos componentes

$$\mathrm{\Psi }(x)=\psi _E(\stackrel{}{x})e^{iEt}=\left(\begin{array}{c}\varphi (\stackrel{}{x})\\ \chi (\stackrel{}{x})\end{array}\right)e^{iEt}.$$
(7)
Obtenemos entonces

$$\widehat{H}\left(\begin{array}{c}\varphi (\stackrel{}{x})\\ \chi (\stackrel{}{x})\end{array}\right)=E\left(\begin{array}{c}\varphi (\stackrel{}{x})\\ \chi (\stackrel{}{x})\end{array}\right),$$
(8)
donde el operador hamiltoniano está dado, en el caso de campo eléctrico nulo, por

$$H=\gamma ^0\stackrel{}{\gamma }(i\stackrel{}{}e\stackrel{}{A})+m\gamma ^0=\left(\begin{array}{cc}H_+& i\sigma _2_z\\ i\sigma _2_z& H_{}\end{array}\right),$$
(9)
con

$$H_\pm =\left(\begin{array}{cc}m& i_x\pm _y\frac{eB}{2}(y\pm ix)\\ i_x_y\frac{eB}{2}(yix)& m\end{array}\right)$$
(10)
o, en coordenadas cilíndricas

$$H_\pm =\left(\begin{array}{cc}m& ie^{i\theta }\left(_r\frac{i}{r}_\theta \frac{eB}{2}r\right)\\ ie^{\pm i\theta }\left(_r\pm \frac{i}{r}_\theta \pm \frac{eB}{2}r\right)& m\end{array}\right).$$
(11)
En lo siguiente estudiaremos, en primer lugar, el problema de un campo de Dirac confinado en un plano perpendicular al campo magnético. A continuación analizaremos el caso en el que todo el espacio es accesible a los fermiones y, finalmente, estudiaremos el efecto de condiciones de contorno en la dirección del campo.

## III Espectro en $`2+1`$ dimensiones

Si los electrones están confinados en un plano perpendicular al campo magnético, entonces las autofunciones resultan independientes de $`z`$

$$\mathrm{\Psi }(x)=\left(\begin{array}{c}\varphi (r,\theta )\\ \chi (r,\theta )\end{array}\right)e^{iEt},$$
(12)
y el operador hamiltoniano resulta diagonal en bloques

$$H=\left(\begin{array}{cc}H_+& 0\\ 0& H_{}\end{array}\right),$$
(13)
con $`H_\pm `$ igual que en la ecuación (11).

Cada uno de estos bloques corresponde a una de las representaciones no equivalentes de las matrices de Dirac en $`2+1`$ dimensiones .

Estudiaremos, en primer lugar, el problema

$$H_+\varphi (r,\theta )=\left(\begin{array}{cc}m& ie^{i\theta }\left(_r\frac{i}{r}_\theta \frac{eB}{2}r\right)\\ ie^{i\theta }\left(_r+\frac{i}{r}_\theta +\frac{eB}{2}r\right)& m\end{array}\right)\varphi (r,\theta )=E\varphi (r,\theta ).$$
(14)
Dado que vamos a obtener una base común de autofunciones de los operadores hamiltoniano y componente $`z`$ del momento angular total, proponemos a $`\varphi `$ autofunción de este último

$$\varphi _k=\left(\begin{array}{c}f_k(r)e^{ik\theta }\\ g_k(r)e^{i(k+1)\theta }\end{array}\right),\text{con}k𝐙.$$
(15)
Entonces, se satisface el siguiente sistema de ecuaciones de primer orden acoplado

$$i\left(\frac{d}{dr}+\frac{k+1}{r}\frac{eB}{2}r\right)g_k(r)=(Em)f_k(r)$$
(16)
$$i\left(\frac{d}{dr}\frac{k}{r}+\frac{eB}{2}r\right)f_k(r)=(E+m)g_k(r),$$
(17)
o bien, definiendo la variable adimensional $`\rho =\left(\frac{eB}{2}\right)^{1/2}r`$

$$i\left(\frac{eB}{2}\right)^{1/2}\left(\frac{d}{d\rho }+\frac{k+1}{\rho }\rho \right)g_k(\rho )=(Em)f_k(\rho )$$
(18)
$$i\left(\frac{eB}{2}\right)^{1/2}\left(\frac{d}{d\rho }\frac{k}{\rho }+\rho \right)f_k(\rho )=(E+m)g_k(\rho ).$$
(19)
Analizaremos, en primer lugar, el problema $`Em`$. En este caso deberemos resolver

$$\left[\frac{1}{\rho }\frac{d}{d\rho }\left(\rho \frac{d}{d\rho }\right)+\frac{k^2}{\rho ^2}+\rho ^2\right]f_k(\rho )=\left(\frac{2}{eB}\right)\left[E^2m^2+eB(k+1)\right]f_k(\rho ),$$
(20)
$$g_k(\rho )=\frac{i}{E+m}\left(\frac{eB}{2}\right)^{1/2}\left(\frac{d}{d\rho }\frac{k}{\rho }+\rho \right)f_k(\rho ).$$
(21)
A fin de resolver la ecuación de segundo orden, hacemos el siguiente cambio de variables

$$\xi =\rho ^2d\xi =2\rho d\rho \frac{1}{\rho }\frac{d}{d\rho }=2\frac{d}{d\xi }\text{y}\rho \frac{d}{d\rho }=2\xi \frac{d}{d\xi },$$
(22)
y proponemos soluciones de la forma

$$f_k\left(\rho (\xi )\right)=e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}F(\xi ).$$
(23)
Al reemplazarlas en la ecuación diferencial, obtenemos

$`\left({\displaystyle \frac{2}{eB}}\right)\left[E^2m^2+eB(k+1)\right]e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}F(\xi )`$ $`=`$ $`\left[4{\displaystyle \frac{d}{d\xi }}\left(\xi {\displaystyle \frac{d}{d\xi }}\right)+{\displaystyle \frac{k^2}{\xi }}+\xi \right]e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}F(\xi )`$ (24)
$`=`$ $`\left(4\xi {\displaystyle \frac{d^2}{d\xi ^2}}4{\displaystyle \frac{d}{d\xi }}+{\displaystyle \frac{k^2}{\xi }}+\xi \right)e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}F(\xi ).`$ (26)
Dado que

$$\frac{d}{d\xi }\left(e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}F(\xi )\right)=e^{\frac{\xi }{2}}\left(\frac{\xi ^{\frac{k}{2}}F(\xi )}{2}+\frac{k\xi ^{\frac{k}{2}1}F(\xi )}{2}\xi ^{\frac{k}{2}}\frac{dF}{d\xi }(\xi )\right),$$
(27)
$$\frac{d^2}{d\xi ^2}\left(e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}F(\xi )\right)=e^{\frac{\xi }{2}}\xi ^{\frac{k}{2}}\left[\frac{d^2F}{d\xi ^2}+\left(\frac{k}{\xi }1\right)\frac{dF}{d\xi }+\left(\frac{1}{4}\frac{k}{2\xi }+\frac{\frac{k}{2}\left(\frac{k}{2}1\right)}{\xi ^2}\right)\right].$$
(28)
$`F(\xi )`$ satisface la ecuación diferencial

$$\xi \frac{d^2F}{d\xi ^2}+(k+1\xi )\frac{dF}{d\xi }+\frac{E^2m^2}{2eB}F(\xi )=0,$$
(29)
que puede identificarse con la ecuación de Kummer .

$$z\frac{d^2w}{dz^2}+(bz)\frac{dw}{dz}aw=0,$$
(30)
donde $`b=k+1`$ y $`a=\frac{m^2E^2}{2eB}`$.

Las soluciones de esta ecuación, estudiadas detalladamente en el Apéndice A, conducen a autofunciones de cuadrado integrable sólo en los casos en que $`a`$ es un entero negativo o cero, y $`ba+1`$. Estas soluciones son polinomios generalizados de Laguerre. Cuando $`a=n`$ los autovalores permitidos para la energía deben satisfacer

$$E^2=m^2+2eBn\text{con}n=0,1,2,\mathrm{}$$
(31)
mientras que las soluciones $`F(\xi )`$ resultan

$$F(\xi )=\{\begin{array}{cc}C_{k,n}.L_n^k(\xi )\hfill & \text{si }k0\hfill \\ C_{k,n}.\xi ^kL_{n+k}^k(\xi )\hfill & \text{si }nk<0\hfill \end{array},$$
(32)
donde $`C_{k,n}`$ es una constante compleja que se ajusta en cada caso para que los estados resulten de norma 1. Entonces, las soluciones de la ecuación (20) son

$$f_k(\rho )=\{\begin{array}{cc}C_{k,n}.e^{\frac{\rho ^2}{2}}\rho ^kL_n^k(\rho ^2)\hfill & \text{si }k0\hfill \\ C_{k,n}.e^{\frac{\rho ^2}{2}}\rho ^kL_{n+k}^k(\rho ^2)\hfill & \text{si }nk<0\hfill \end{array}.$$
(33)
Usando relaciones de recurrencia entre los polinomios generalizados de Laguerre , es fácil obtener las soluciones de la ecuación (21), para $`n0`$

$$g_k(\rho )=\{\begin{array}{cc}iC_{k,n}\frac{(2eB)^{\frac{1}{2}}}{E+m}e^{\frac{\rho ^2}{2}}\rho ^{k+1}L_{n1}^{k+1}(\rho ^2)\hfill & \text{si }k0\hfill \\ iC_{k,n}\frac{(2eB)^{\frac{1}{2}}}{E+m}n.e^{\frac{\rho ^2}{2}}\rho ^{k1}L_{n+k}^{k1}(\rho ^2)\hfill & \text{si }nk<0\hfill \end{array},$$
(34)
y $`g_k(\rho )=0`$ para $`n=0`$.

Ahora analizaremos el caso $`E=m`$. En este caso, la ecuación (19) queda desacoplada

$$\left(\frac{d}{d\rho }\frac{k}{\rho }+\rho \right)f_k(\rho )=0,$$
(35)
y su solución es simplemente

$$f_k(\rho )=C_k.e^{\frac{\rho ^2}{2}}\rho ^k,$$
(36)
donde $`k`$ solamente puede tomar valores enteros positivos para que la componente superior de $`\varphi `$ resulte de cuadrado integrable, y $`C_k`$ es una constante de normalización.

Por su parte, la ecuación (18) es una ecuación de primer orden inhomogénea

$$\left(\frac{d}{d\rho }+\frac{k+1}{\rho }\rho \right)g_k(\rho )=\frac{2miC_k}{\left(\frac{eB}{2}\right)^{1/2}}e^{\frac{\rho ^2}{2}}\rho ^k.$$
(37)
Proponemos soluciones de la forma

$$g_k(\rho )=G(\rho )e^{\frac{\rho ^2}{2}}\rho ^{(k+1)},$$
(38)
por lo que $`G(\rho )`$ debe satisfacer

$$\frac{dG}{d\rho }=\frac{2miC_k}{\left(\frac{eB}{2}\right)^{1/2}}e^{\rho ^2}\rho ^{2k+1},$$
(39)
$$G(\rho )=\frac{2miC_k}{\left(\frac{eB}{2}\right)^{1/2}}𝑑\rho e^{\rho ^2}\rho ^{2k+1},$$
(40)
y haciendo el cambio de variables $`\rho ^2=x`$

$`G(\rho (x))`$ $`=`$ $`{\displaystyle \frac{miC_k}{\left(\frac{eB}{2}\right)^{1/2}}}{\displaystyle 𝑑xe^xx^k}`$ (41)
$`=`$ $`{\displaystyle \frac{miC_k}{\left(\frac{eB}{2}\right)^{1/2}}}e^x{\displaystyle \underset{n=0}{\overset{k}{}}}{\displaystyle \frac{k!}{(kn)!}}x^{kn}+Cte.`$ (42)
Por lo tanto, la solución general de (37) es la solución general de su homogénea correspondiente más una solución particular de la inhomogénea

$$g_k(\rho )=\frac{miC_k}{\left(\frac{eB}{2}\right)^{1/2}}e^{\frac{\rho ^2}{2}}\underset{n=0}{\overset{k}{}}\frac{k!}{(kn)!}\rho ^{k2n1}+De^{\frac{\rho ^2}{2}}\rho ^{(k+1)},$$
(43)
la cual conduce, para todo valor de de la constante $`D`$, a autofunciones de la componente inferior de $`\varphi `$ que no son de cuadrado integrable.

Por lo tanto, los autovalores de $`H_+`$ que corresponden a autofunciones de cuadrado integrable son los siguientes:

$$E_+=\{\begin{array}{c}m\hfill \\ \pm \sqrt{m^2+2eBn}\text{con}n=1,2,3,\mathrm{}\hfill \end{array}$$
(44)
El problema

$$H_{}\chi (r,\theta )=\left(\begin{array}{cc}m& ie^{i\theta }\left(_r+\frac{i}{r}_\theta +\frac{eB}{2}r\right)\\ ie^{i\theta }\left(_r\frac{i}{r}_\theta \frac{eB}{2}r\right)& m\end{array}\right)\chi (r,\theta )=E\chi (r,\theta ),$$
(45)
se resuelve en forma enteramente análoga.

En efecto, factorizando la dependencia angular, obtenemos un sistema de ecuaciones de primer orden acopladas similar al de (18) y (19), que puede desacoplarse en una ecuación de segundo orden para el caso $`Em`$. Nuevamente, los autovalores de la energía satisfacen

$$E^2=m^2+2eBn\text{con}n=0,1,2,\mathrm{}$$
(46)
y las autofunciones pueden expresarse en términos de polinomios generalizados de Laguerre.

Por último, al estudiar las soluciones del caso $`E=m`$, se encuentra que este autovalor no tiene asociadas autofunciones de cuadrado integrable.

Finalmente, nos quedan los siguientes autovalores de $`H_{}`$ correspondientes a autofunciones de cuadrado integrable:

$$E_{}=\{\begin{array}{c}m\hfill \\ \pm \sqrt{m^2+2eBn}\text{con}n=1,2,3,\mathrm{}\hfill \end{array}$$
(47)
## IV Energía de vacío en $`2+1`$ dimensiones

Como se explicó en la sección I, la energía de vacío está dada por la ecuación (2)

$$E_0=\frac{1}{2}\underset{n}{}g_n.|E_n|,$$
(48)
donde los $`E_n`$ son los autovalores dados por las ecuaciones (44) y (47). El factor de degeneración $`g_n`$ proviene de los distintos valores que puede tomar $`k`$ (cualquier entero mayor o igual que n). Éste puede calcularse contando el número de estados que se encuentran en el problema de la partícula libre con módulo cuadrado del impulso lineal entre $`2eBn`$ y $`2eB(n+1)`$, si pensamos que al encender el campo todos estos niveles colapsan en uno. Integrando la densidad de estados de la partícula libre, resulta

$$g_n=\frac{L_xL_y}{(2\pi )^2}𝑑p_x𝑑p_y=\frac{A}{(2\pi )^2}2\pi _{\sqrt{2eBn}}^{\sqrt{2eB(n+1)}}p𝑑p=\frac{A}{2\pi }eB,$$
(49)
donde $`A`$ es el área perpendicular al campo magnético.

Conviene, entonces, considerar la energía de vacío por unidad de área

$$\epsilon _0=\frac{E_0}{A}=\frac{eB}{4\pi }\underset{n}{}|E_n|=\frac{eB}{4\pi }2\left(m+2\underset{n=1}{\overset{\mathrm{}}{}}\sqrt{m^2+2eBn}\right).$$
(50)
Como es usual, este resultado es sólo formal y para darle sentido debemos regularizarlo. Hacemos esto con el método de regularización de la función zeta, descrito en la sección I. Debemos multiplicar a la suma por un parámetro para que conserve unidades de energía por unidad de área. Usamos $`\mu `$ con unidades de masa

$`\epsilon _0`$ $`=`$ $`{\displaystyle \frac{2eB}{4\pi }}\left\{\mu ^{s+1}\left[m^s+2{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left(m^2+2eBn\right)^{s/2}\right]\right\}_{s=1}`$ (51)
$`=`$ $`{\displaystyle \frac{2eB}{2\pi }}\left\{\mu ^{s+1}\left[{\displaystyle \frac{m^s}{2}}(2eB)^{s/2}{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}\left(n+{\displaystyle \frac{m^2}{2eB}}\right)^{s/2}\right]\right\}_{s=1}.`$ (53)
Es conveniente introducir los parámetros adimensionales

$$\{\begin{array}{c}b=\frac{2eB}{\mu ^2}\hfill \\ c=\frac{m^2}{\mu ^2}\hfill \end{array},$$
(54)
de modo que

$$\epsilon _0=\frac{b}{2\pi }\mu ^3\left[\frac{c^{s/2}}{2}b^{s/2}\underset{n=0}{\overset{\mathrm{}}{}}\left(n+\frac{c}{b}\right)^{s/2}\right]_{s=1},$$
(55)
Para $`Re(s)>2`$, la serie en la ecuación (55) converge y coincide con la función zeta Hurwitz.

$$\zeta _H(s,a)=\underset{n=0}{\overset{\mathrm{}}{}}(n+a)^s,$$
(56)
que es convergente para $`Re(s)>1`$, pero se extiende analíticamente a una función meromorfa con un polo simple en $`s=1`$. Por lo tanto, la densidad de energía de vacío puede definirse, por extensión analítica, como

$$\epsilon _0=\frac{b}{2\pi }\mu ^3\left[\frac{c^{s/2}}{2}(b)^{s/2}\zeta _H(s/2,\frac{c}{b})\right]_{s=1},$$
(57)
Dado que $`\zeta _H(s/2,\frac{c}{b})`$ tiene un único polo simple en $`s=2`$, la densidad de energía de vacío por unidad de área resulta finita e igual a

$$\epsilon _0=\frac{2eB}{2\pi }\left[\frac{m}{2}(2eB)^{1/2}\zeta _H(1/2,\frac{m^2}{2eB})\right].$$
(58)
Vale la pena señalar que, en el caso de un campo puramente magnético, la energía de vacío debe coincidir con la densidad lagrangiana efectiva al orden 1 loop del sistema. En efecto, el resultado de la ecuación (58) coincide con la densidad lagrangiana efectiva obtenida en la referencia , a partir del cálculo con regularización zeta del determinante del operador de Dirac.

Para calcular la densidad de energía de Casimir, debemos conocer la densidad de energía de vacío del problema sin campo. A fin de evaluar (58) para $`B=0`$, usaremos un desarrollo asintótico de la función $`\zeta _H(s,a)`$, válido cuando $`|a|\mathrm{}`$ .

$`\zeta _H(s,a)`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Gamma }(s)}}[a^{(s1)}\mathrm{\Gamma }(s1)+{\displaystyle \frac{1}{2}}a^s\mathrm{\Gamma }(s)+`$ (61)
$`+{\displaystyle \underset{n=1}{\overset{m}{}}}{\displaystyle \frac{B_{2n}}{(2n)!}}a^{(2n+s1)}\mathrm{\Gamma }(2n+s1)+O(a^{(2m+s+1)})],`$
donde $`B_{2n}`$ son los coeficientes de Bernoulli.

Al escribir la densidad de energía de vacío en términos de este desarrollo asintótico resulta, para $`eBm^2`$

$$\epsilon _0=\frac{m^3}{3\pi }\frac{m^3}{2\pi }\frac{1}{\mathrm{\Gamma }(1/2)}\underset{n=1}{\overset{m}{}}\frac{B_{2n}}{(2n)!}\left(\frac{2eB}{m^2}\right)^{2n}\mathrm{\Gamma }(2n3/2)+O\left[\left(\frac{2eB}{m^2}\right)^{2(m+1)}\right].$$
(62)
Vemos que el primer término, independiente del campo magnético, es $`\frac{m^3}{3\pi }`$, y por lo tanto corresponde a la densidad de energía de vacío del problema sin campo. En el Apéndice B obtenemos el mismo resultado sumando directamente los autovalores del espectro clásico de energías de la partícula libre.

Resumiendo, la densidad de energía de Casimir resulta, para valores arbitrarios del campo

$`\epsilon _{Cas}`$ $`=`$ $`\epsilon _0\epsilon _0^l`$ (63)
$`=`$ $`{\displaystyle \frac{2eB}{2\pi }}\left[{\displaystyle \frac{m}{2}}(2eB)^{1/2}\zeta _H(1/2,{\displaystyle \frac{m^2}{2eB}})\right]{\displaystyle \frac{m^3}{3\pi }}.`$ (65)
Nótese que, para campos pequeños, el segundo término del desarrollo (61) cancela, en la ecuación (58), el término proporcional al campo. De este modo, resulta para la energía de Casimir un desarrollo en potencias pares del campo magnético.

## V Espectro en $`3+1`$ dimensiones

En lugar de resolver explícitamente (8) podemos aplicar una transformación de Lorentz sobre las soluciones independientes de $`z`$. En particular, considaremos un boost en la dirección del eje $`z`$, desde el sistema de referencia en el que $`p_z=0`$ hasta el sistema en el cual $`p_z^{}0`$

$$\mathrm{\Lambda }^\mu \nu =\left(\begin{array}{cccc}\mathrm{cosh}\omega & 0& 0& \mathrm{sinh}\omega \\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ \mathrm{sinh}\omega & 0& 0& \mathrm{cosh}\omega \end{array}\right),\text{donde}\mathrm{cosh}\omega =\frac{1}{\sqrt{1v^2}}.$$
(66)
Entonces

$$\{\begin{array}{cc}E^{}& =\mathrm{cosh}\omega E\\ p_z^{}& =\mathrm{sinh}\omega E\end{array},$$
(67)
por lo tanto

$$E_{}^{}{}_{}{}^{2}=\mathrm{cosh}^2\omega E^2=(1+\mathrm{sinh}^2\omega )E^2=E^2+p_{z}^{}{}_{}{}^{2}.$$
(68)
Por su parte, los espinores se transforman según

$$\mathrm{\Psi }^{}=S\mathrm{\Psi },$$
(69)
donde

$$S=\mathrm{exp}\left(i\frac{\omega }{2}\frac{i}{2}[\gamma ^0,\gamma ^3]\right)=\mathrm{cosh}\left(\frac{\omega }{2}\right)\mathrm{𝟏}_4+\mathrm{sinh}\left(\frac{\omega }{2}\right)\gamma ^0\gamma ^3.$$
(70)
En la representación elegida para las matrices de Dirac resulta

$$S=\left(\begin{array}{cc}\mathrm{cosh}\left(\frac{\omega }{2}\right)\mathrm{𝟏}_2& \mathrm{sinh}\left(\frac{\omega }{2}\right)\sigma _2\\ \mathrm{sinh}\left(\frac{\omega }{2}\right)\sigma _2& \mathrm{cosh}\left(\frac{\omega }{2}\right)\mathrm{𝟏}_2\end{array}\right).$$
(71)
A partir de las soluciones

$$\mathrm{\Psi }=\left(\begin{array}{c}\varphi (r,\theta )\\ 0\end{array}\right)e^{iE_+t},$$
(72)
donde $`\varphi `$ es una de las autofunciones del bloque $`H_+`$, se obtiene

$$\mathrm{\Psi }^{}=\left(\begin{array}{c}\mathrm{cosh}\left(\frac{\omega }{2}\right)\varphi (r,\theta )\\ \mathrm{sinh}\left(\frac{\omega }{2}\right)\sigma _2\varphi (r,\theta )\end{array}\right)e^{iE_+t},$$
(73)
o bien, en términos de las coordenados del sistema primado,

$$\mathrm{\Psi }^{}=\left(\begin{array}{c}\mathrm{cosh}\left(\frac{\omega }{2}\right)\varphi (r^{},\theta ^{})\\ \mathrm{sinh}\left(\frac{\omega }{2}\right)\sigma _2\varphi (r^{},\theta ^{})\end{array}\right)e^{ip_z^{}z^{}}e^{iE_+^{}t^{}}.$$
(74)
Estas autofunciones corresponden a los autovalores de la energía

$$E_+^{}=\{\begin{array}{c}\sqrt{m^2+p_{z}^{}{}_{}{}^{2}}\hfill \\ \pm \sqrt{m^2+p_{z}^{}{}_{}{}^{2}+2eBn}\hfill \end{array}\text{con}p_z^{}𝐑\text{y}n=1,2,3,\mathrm{}$$
(75)
Similarmente, a partir de las autofunciones $`\chi `$ del bloque $`H_{}`$, se obtienen los siguientes espinores en el sistema de referencia primado

$$\mathrm{\Psi }^{}=\left(\begin{array}{c}\mathrm{sinh}\left(\frac{\omega }{2}\right)\sigma _2\chi (r^{},\theta ^{})\\ \mathrm{cosh}\left(\frac{\omega }{2}\right)\chi (r^{},\theta ^{})\end{array}\right)e^{ip_z^{}z^{}}e^{iE_{}^{}t^{}},$$
(76)
que corresponden a autovalores

$$E_{}^{}=\{\begin{array}{c}\sqrt{m^2+p_{z}^{}{}_{}{}^{2}}\hfill \\ \pm \sqrt{m^2+p_{z}^{}{}_{}{}^{2}+2eBn}\hfill \end{array}\text{con}p_z^{}𝐑\text{y}n=1,2,3,\mathrm{}$$
(77)
## VI Energía de vacío en $`3+1`$ dimensiones

En este caso, la energía de vacío se obtiene a partir de sumar los autovalores dados en la ecuaciones (75) y (77)

$$E_0=g\frac{L_z}{2\pi }_{\mathrm{}}^{\mathrm{}}𝑑p_z\left(\sqrt{m^2+p_z^2}+2\underset{n=1}{\overset{\mathrm{}}{}}\sqrt{m^2+p_z^2+2eBn}\right),$$
(78)
donde el factor de degeneración $`g`$ es nuevamente $`\frac{A}{2\pi }eB`$.

Redefinimos el índice de la suma para que ésta comience desde 0 y definimos la densidad de energía de vacío por unidad de volumen

$$\epsilon _0=\frac{E_0}{AL_z}=\frac{eB}{(2\pi )^2}_{\mathrm{}}^{\mathrm{}}𝑑p_z\left(\sqrt{m^2+p_z^2}2\underset{n=0}{\overset{\mathrm{}}{}}\sqrt{m^2+p_z^2+2eBn}\right).$$
(79)
Nuevamente, la expresión que encontramos para la densidad de energía de vacío es formalmente divergente. Repitiendo el procedimiento antes utilizado, la extendemos analíticamente con la técnica de regularización zeta. Como antes debemos agregar algún parámetro para que la densidad de energía de vacío mantenga las unidades apropiadas. Usando $`\mu `$ con unidades de masa, resulta

$$\epsilon _0=\frac{eB}{(2\pi )^2}\left\{\mu ^{s+1}_{\mathrm{}}^{\mathrm{}}𝑑p_z\left[(m^2+p_z^2)^{s/2}2\underset{n=0}{\overset{\mathrm{}}{}}\left(m^2+p_z^2+2eBn\right)^{s/2}\right]\right\}_{s=1}.$$
(80)
Si usamos la transformada de Mellin

$$Z^s=\frac{1}{\mathrm{\Gamma }(s)}_0^{\mathrm{}}𝑑tt^{s1}e^{Zt},$$
(81)
la ecuación (80) se escribe

$`\epsilon _0`$ $`=`$ $`{\displaystyle \frac{eB}{(2\pi )^2}}\{\mu ^{s+1}{\displaystyle _{\mathrm{}}^{\mathrm{}}}dp_z[{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+p_z^2)t]`$ (84)
$`2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+p_z^2+2eBn)t]]\}_{s=1}.`$
El orden de integración puede intercambiarse para $`Re(s)`$ suficientemente grande, de modo que

$`\epsilon _0`$ $`=`$ $`{\displaystyle \frac{eB}{(2\pi )^2}}\{\mu ^{s+1}[{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}(m^2t){\displaystyle _{\mathrm{}}^{\mathrm{}}}dp_z\mathrm{exp}(p_z^2t)`$ (87)
$`2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+2eBn)t]{\displaystyle _{\mathrm{}}^{\mathrm{}}}dp_z\mathrm{exp}(p_z^2t)]\}_{s=1}.`$
Usando el resultado de la siguiente integral definida

$$_{\mathrm{}}^{\mathrm{}}𝑑xe^{ax^2}=\sqrt{\frac{\pi }{a}},$$
(88)
resolvemos las integrales a lo largo de $`p_z`$. Así obtenemos

$`\epsilon _0`$ $`=`$ $`{\displaystyle \frac{eB}{(2\pi )^2}}\{\mu ^{s+1}[{\displaystyle \frac{\sqrt{\pi }}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}(m^2t)`$ (91)
$`2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\sqrt{\pi }}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}[(m^2+2eBn)t]]\}_{s=1}.`$
Luego de realizar las integrales en $`t`$, resulta

$$\epsilon _0=\frac{eB}{(2\pi )^2}\left\{\mu ^{s+1}\left[\sqrt{\pi }\frac{\mathrm{\Gamma }(\frac{s1}{2})}{\mathrm{\Gamma }(s/2)}(m^2)^{\frac{s1}{2}}2\sqrt{\pi }\frac{\mathrm{\Gamma }(\frac{s1}{2})}{\mathrm{\Gamma }(s/2)}\underset{n=0}{\overset{\mathrm{}}{}}(m^2+2eBn)^{\frac{s1}{2}}\right]\right\}_{s=1},$$
(92)
Nuevamente, la serie en el segundo término define, para $`Re(s)>3`$, una zeta de Hurwitz. Por extensión analítica, obtenemos

$$\epsilon _0=\frac{2eB}{(2\pi )^2}\sqrt{\pi }\left\{\mu ^{s+1}\frac{\mathrm{\Gamma }(\frac{s1}{2})}{\mathrm{\Gamma }(s/2)}\left[\frac{(m^2)^{\frac{s1}{2}}}{2}(2eB)^{\frac{s1}{2}}\zeta _H(\frac{s1}{2},\frac{m^2}{2eB})\right]\right\}_{s=1}.$$
(93)
A diferencia de lo que ocurría en $`2+1`$ dimensiones, la extensión de la función zeta del hamiltoniano clásico presenta, en $`3+1`$, un polo simple en el valor $`s=1`$, debido al comportamiento de la función $`\mathrm{\Gamma }(s)`$ en los enteros negativos.

Usando nuevamente los parámetros adimensionales definidos en (54), (93) puede escribirse

$`\epsilon _0`$ $`=`$ $`{\displaystyle \frac{\sqrt{\pi }}{(2\pi )^2}}b\mu ^4\left\{{\displaystyle \frac{\mathrm{\Gamma }(\frac{s1}{2})}{\mathrm{\Gamma }(s/2)}}\left[{\displaystyle \frac{c^{\frac{s1}{2}}}{2}}b^{\frac{s1}{2}}\zeta _H({\displaystyle \frac{s1}{2}},{\displaystyle \frac{c}{b}})\right]\right\}_{s=1}`$ (94)
$`=`$ $`{\displaystyle \frac{\sqrt{\pi }}{\pi ^2}}b\mu ^4\left\{{\displaystyle \frac{\mathrm{\Gamma }(\frac{s+3}{2})}{\mathrm{\Gamma }(s/2)(s+1)(s1)}}\left[{\displaystyle \frac{c^{\frac{s1}{2}}}{2}}b^{\frac{s1}{2}}\zeta _H({\displaystyle \frac{s1}{2}},{\displaystyle \frac{c}{b}})\right]\right\}_{s=1},`$ (96)
o bien

$$\epsilon _0=\frac{\sqrt{\pi }}{\pi ^2}b\mu ^4\frac{Z(s)}{s+1}|_{s=1},$$
(97)
donde

$$Z(s)=\frac{\mathrm{\Gamma }(\frac{s+3}{2})}{\mathrm{\Gamma }(s/2)(s1)}\left[\frac{c^{\frac{s1}{2}}}{2}b^{\frac{s1}{2}}\zeta _H(\frac{s1}{2},\frac{c}{b})\right].$$
(98)
A fin de identificar las partes singular y finita de la densidad de energía de vacío, hacemos un desarrollo alrededor de $`s=1`$.

$$\epsilon _0=\frac{\sqrt{\pi }}{\pi ^2}b\mu ^4\left[\frac{Z(1)}{s+1}+Z^{}(1)\right]_{s=1}.$$
(99)
Para la parte singular debemos calcular

$$Z(1)=\frac{\mathrm{\Gamma }(1)}{(2)\mathrm{\Gamma }(1/2)}\left[\frac{c}{2}b\zeta _H(1,\frac{c}{b})\right],$$
(100)
y usando la extensión analítica de la función Zeta de Hurwitz

$$\zeta _H(1,a)=\frac{1}{2}\left(a^2a+\frac{1}{6}\right),$$
(101)
$$Z(1)=\frac{b}{8\sqrt{\pi }}\left(\frac{c^2}{b^2}+\frac{1}{6}\right).$$
(102)
La parte finita está dada por

$`Z^{}(1)`$ $`=`$ $`{\displaystyle \frac{b}{16\sqrt{\pi }}}\left(\psi (1)\psi (1/2)+1\right)\left[{\displaystyle \frac{c}{b}}2\zeta _H(1,{\displaystyle \frac{c}{b}})\right]`$ (105)
$`{\displaystyle \frac{b}{16\sqrt{\pi }}}\left({\displaystyle \frac{c}{b}}\mathrm{ln}c2\mathrm{ln}b\zeta _H(1,{\displaystyle \frac{c}{b}})+2\zeta _H^{}(1,{\displaystyle \frac{c}{b}})\right)`$
$`=`$ $`{\displaystyle \frac{b}{16\sqrt{\pi }}}\left(\psi (1)\psi (1/2)+1\right)\left({\displaystyle \frac{c^2}{b^2}}+{\displaystyle \frac{1}{6}}\right)`$ (109)
$`{\displaystyle \frac{b}{16\sqrt{\pi }}}\left[{\displaystyle \frac{c}{b}}\mathrm{ln}c+\mathrm{ln}b\left({\displaystyle \frac{c^2}{b^2}}{\displaystyle \frac{c}{b}}+{\displaystyle \frac{1}{6}}\right)+2\zeta _H^{}(1,{\displaystyle \frac{c}{b}})\right].`$
Reemplazando (102) y (105) en (99), la densidad de energía de vacío resulta

$`\epsilon _0`$ $`=`$ $`\mu ^4\left[{\displaystyle \frac{c^2}{8\pi ^2(s+1)}}+{\displaystyle \frac{\frac{b^2}{6}}{8\pi ^2(s+1)}}\right]_{s=1}+`$ (114)
$`+{\displaystyle \frac{\mu ^4}{16\pi ^2}}[(\psi (1)\psi (1/2)+1\mathrm{ln}c)(c^2+{\displaystyle \frac{b^2}{6}})+`$
$`+(c^2cb+{\displaystyle \frac{b^2}{6}})\mathrm{ln}\left({\displaystyle \frac{c}{b}}\right)2b^2\zeta _H^{}(1,{\displaystyle \frac{c}{b}})].`$
Vale la pena señalar, en este punto, que la parte finita de la expresión anterior coincide, a menos de términos proporcionales a $`b^2`$ (y, como veremos, renormalizables) con la densidad lagrangiana efectiva calculada en  a partir del determinante del operador de Dirac. La diferencia entre ambos resultados está de acuerdo con la discusión general presentada en .

De la expresión (114), vemos que el residuo en el polo presenta un término independiente del campo magnético, y otro proporcional a su cuadrado.

A fin de dar una interpretación física a este resultado verificaremos, en primer lugar, que el término independiente de $`b`$ se cancela con uno idéntico que aparece en la densidad de energía de vacío a $`B=0`$.

Para hallar la densidad de energía de vacío sin campo $`\epsilon _0^l`$, podemos considerar un desarrollo de asintótico de la derivada de la función zeta en el límite de los campos pequeños, y quedarnos con los términos del resultado final independientes del campo (en el Apéndice C la reobtendremos directamente a partir de la suma de los autovalores del espectro clásico de la partícula libre).

El siguiente desarrollo para la derivada de la función zeta de Hurwitz es válido en el límite $`a\mathrm{}`$ .

$$\zeta _H^{}(1,a)=\frac{1}{2}\left(a^2a+\frac{1}{6}\right)\mathrm{ln}a\frac{1}{4}a^2+\frac{1}{12}+\underset{n=1}{\overset{m1}{}}\frac{B_{2n+2}a^{2n}}{(2n)(2n+1)(2n+2)}+O(a^{2m}).$$
(115)
Haciendo uso del mismo resulta, para $`bc`$

$`Z^{}(1)`$ $`=`$ $`{\displaystyle \frac{b}{16\sqrt{\pi }}}\{(\psi (1)\psi (1/2)+1\mathrm{ln}c)({\displaystyle \frac{c^2}{b^2}}+{\displaystyle \frac{1}{6}})+{\displaystyle \frac{1}{2}}{\displaystyle \frac{c^2}{b^2}}{\displaystyle \frac{1}{6}}`$ (118)
$`2{\displaystyle \underset{n=1}{\overset{m1}{}}}{\displaystyle \frac{B_{2n+2}\left(\frac{c}{b}\right)^{2n}}{(2n)(2n+1)(2n+2)}}+O\left[\left({\displaystyle \frac{b}{c}}\right)^{2m}\right]\}`$
y, por lo tanto,

$$\epsilon _0^l=\mu ^4\left[\frac{c^2}{8\pi ^2(s+1)}\right]_{s=1}+\frac{\mu ^4c^2}{16\pi ^2}\left[\psi (1)\psi (1/2)+3/2\mathrm{ln}c\right].$$
(119)
De este modo, la densidad de energía de Casimir resulta

$`\epsilon _{Cas}`$ $`=`$ $`\epsilon _0\epsilon _0^l`$ (120)
$`=`$ $`\mu ^4\left[{\displaystyle \frac{\frac{b^2}{6}}{8\pi ^2(s+1)}}\right]_{s=1}+{\displaystyle \frac{\mu ^4}{16\pi ^2}}[(\psi (1)\psi (1/2)+1\mathrm{ln}c){\displaystyle \frac{b^2}{6}}{\displaystyle \frac{c^2}{2}}+`$ (124)
$`+(c^2cb+{\displaystyle \frac{b^2}{6}})\mathrm{ln}\left({\displaystyle \frac{c}{b}}\right)2b^2\zeta _H^{}(1,{\displaystyle \frac{c}{b}})].`$
Esta expresión presenta aún una divergencia, cuya dependencia funcional con el campo es igual a la que aparece en la densidad clásica de energía

$$\epsilon _{cl\stackrel{´}{a}}=\frac{1}{2}B^2=\frac{\mu ^4b^2}{8e^2}$$
(125)
A fin de eliminar la singularidad que persiste en la densidad de energía de Casimir, renormalizamos la carga $`e`$. Para ello, consideramos a los términos proporcionales al cuadrado del campo como la corrección cuántica de la densidad de energía electromagnética clásica.

De esta manera, definimos la carga renornalizada $`e_R`$ a través de

$$\frac{(\epsilon _{cl\stackrel{´}{a}}+\epsilon _0)}{b^2}|_{b^2=0}=\frac{\mu ^4}{8e_R^2}$$
(126)
Evaluamos esta derivada a partir de la ecuación (94) y obtenemos

$`{\displaystyle \frac{(\epsilon _{cl\stackrel{´}{a}}+\epsilon _0)}{b^2}}`$ $`=`$ $`{\displaystyle \frac{\mu ^4}{8e^2}}+{\displaystyle \frac{\sqrt{\pi }}{(2\pi )^2}}\mu ^4\{{\displaystyle \frac{\mathrm{\Gamma }(\frac{s1}{2})}{\mathrm{\Gamma }(s/2)}}[\left({\displaystyle \frac{s3}{4}}\right)b^{\frac{s+1}{2}}\zeta _H({\displaystyle \frac{s1}{2}},{\displaystyle \frac{c}{b}})`$ (129)
$`\left({\displaystyle \frac{s1}{4}}\right)cb^{\frac{s+3}{2}}\zeta _H({\displaystyle \frac{s+1}{2}},{\displaystyle \frac{c}{b}})+{\displaystyle \frac{c^{\frac{s1}{2}}b^1}{4}}]\}_{s=1}.`$
donde hemos usado la siguiente propiedad de la función zeta de Hurwitz

$$\frac{d\zeta _H(s,a)}{da}=s\zeta _H(s+1,a)$$
(130)
Podemos evaluar las zetas de Hurwitz en el límite $`b^20`$ con el desarrollo (61). Finalmente, resulta

$$\frac{(\epsilon _{cl\stackrel{´}{a}}+\epsilon _0)}{b^2}|_{b^2=0}=\frac{\mu ^4}{8e^2}\frac{\sqrt{\pi }}{8\pi ^2}\mu ^4B_2c^{\frac{s+1}{2}}\frac{\mathrm{\Gamma }(\frac{s+1}{2})}{\mathrm{\Gamma }(s/2)}|_{s=1}.$$
(131)
Dado que esta expresión es singular en el límite $`s1`$, igual que como hicimos antes, separamos las partes finita y singular mediante un desarrollo alrededor de $`s=1`$, con lo cual obtenemos

$`{\displaystyle \frac{(\epsilon _{cl\stackrel{´}{a}}+\epsilon _0)}{b^2}}|_{b^2=0}`$ $`=`$ $`{\displaystyle \frac{\mu ^4}{8e_R^2}}`$ (132)
$`=`$ $`{\displaystyle \frac{\mu ^4}{8e^2}}+{\displaystyle \frac{\mu ^4}{96\pi ^2}}\left[{\displaystyle \frac{2}{s+1}}+\psi (1)\psi (1/2)\mathrm{ln}c\right]_{s=1},`$ (134)
y de aquí se tiene

$`e_R^2`$ $`=`$ $`{\displaystyle \frac{e^2}{1+\frac{e^2}{12\pi ^2}\left[\frac{2}{s+1}+\psi (1)\psi (1/2)\mathrm{ln}c\right]}}`$ (135)
$``$ $`e^2\left\{1{\displaystyle \frac{e^2}{12\pi ^2}}\left[{\displaystyle \frac{2}{s+1}}+\psi (1)\psi (1/2)\mathrm{ln}c\right]\right\}.`$ (137)
Entonces, el parámetro que aparece en el hamiltoniano debe ser renormalizado según:

$$e^2=e_R^2+\frac{e_R^4}{12\pi ^2}\left[\frac{2}{s+1}+\psi (1)\psi (1/2)\mathrm{ln}c\right]+O(e_R^5).$$
(138)
Finalmente, para la densidad de energía al orden 1 loop:

$`\epsilon _{Cas}`$ $`=`$ $`\epsilon _{cl\stackrel{´}{a}}+\epsilon _0\epsilon _0^l`$ (139)
$`=`$ $`{\displaystyle \frac{\mu ^4b^2}{8e_R^2}}+{\displaystyle \frac{\mu ^4}{16\pi ^2}}\left[{\displaystyle \frac{b^2}{6}}{\displaystyle \frac{c^2}{2}}+\left(c^2cb+{\displaystyle \frac{b^2}{6}}\right)\mathrm{ln}\left({\displaystyle \frac{c}{b}}\right)2b^2\zeta _H^{}(1,{\displaystyle \frac{c}{b}})\right],`$ (141)
y en términos de los parámetros originales

$`\epsilon _{Cas}`$ $`=`$ $`{\displaystyle \frac{B^2}{2}}+{\displaystyle \frac{1}{16\pi ^2}}\{{\displaystyle \frac{(2e_RB)^2}{6}}[\mathrm{ln}\left({\displaystyle \frac{m^2}{2e_RB}}\right)+1]+m^4[\mathrm{ln}\left({\displaystyle \frac{m^2}{2e_RB}}\right){\displaystyle \frac{1}{2}}]`$ (144)
$`m^22e_RB\mathrm{ln}\left({\displaystyle \frac{m^2}{2e_RB}}\right)2(2e_RB)^2\zeta _H^{}(1,{\displaystyle \frac{m^2}{(2e_RB)}})\}.`$
Nótese que, una vez renormalizada la carga, la densidad de energía de Casimir es independiente del parámetro de escala $`\mu `$.

## VII Energía de vacío en $`3+1`$ dimensiones con condiciones de contorno

En esta sección, analizaremos el efecto de tamaño finito en la dirección del campo sobre la energía de vacío. Impondremos sobre los campos de Dirac condiciones de contorno no locales, del tipo

$$\psi (r,\theta ,L)=e^{i\alpha }\psi (r,\theta ,0).$$
(145)
Vemos que esta condición se reduce a la condición de periodicidad para $`\alpha =0`$, mientras que de $`\alpha =\pi `$ se obtiene el problema con condiciones de contorno antiperiódicas.

Nuevamente tenemos autovectores simultáneos de la energía y del impulso en la dirección del campo, sólo que ahora el espectro de autovalores $`p_z`$ pasa a ser discreto.

$$\begin{array}{ccccc}\text{Si}& e^{ip_zL}=e^{i\alpha }& & p_z=\frac{2k\pi +\alpha }{L}& \text{donde }k\text{ es un entero}\end{array}.$$
(146)
Entonces, en la expresión para la energía de vacío, debemos ahora sumar sobre el espectro discreto de autovalores $`p_z`$

$`E_0^L`$ $`=`$ $`{\displaystyle \frac{eBA}{2\pi }}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\left\{\left[m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2\right]^{1/2}+2{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left[m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2+2eBn\right]^{1/2}\right\}`$ (147)
$`=`$ $`{\displaystyle \frac{eBA}{2\pi }}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\left\{2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}\left[m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2+2eBn\right]^{1/2}\left[m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2\right]^{1/2}\right\}.`$ (149)
Una vez regularizada mediante la técnica de la función zeta, (147) se escribe

$`E_0^L`$ $`=`$ $`{\displaystyle \frac{eBA}{2\pi }}(\mu ^{s+1}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\{2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}[m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2+2eBn]^{s/2}`$ (152)
$`[m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2]^{s/2}\})_{s=1}.`$
Usamos la transformada de Mellin (81) para luego poder hacer la suma en $`n`$

$`E_0^L`$ $`=`$ $`{\displaystyle \frac{eBA}{2\pi }}(\mu ^{s+1}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\{{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{2}{\mathrm{\Gamma }(s/2)}}\times `$ (157)
$`\times {\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2+2eBn)t]`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+\left({\displaystyle \frac{2\pi k+\alpha }{L}}\right)^2)t]\})_{s=1}`$
$`=`$ $`{\displaystyle \frac{eBA}{2\pi }}(\mu ^{s+1}\{{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{2}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+\left({\displaystyle \frac{\alpha }{L}}\right)^2+2eBn)t]\times `$ (165)
$`\times {\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\mathrm{exp}({\displaystyle \frac{4\pi \alpha }{L^2}}tk{\displaystyle \frac{4\pi ^2}{L^2}}tk^2)`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}\mathrm{exp}[(m^2+\left({\displaystyle \frac{\alpha }{L}}\right)^2)t]\times `$
$`\times {\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\mathrm{exp}({\displaystyle \frac{4\pi \alpha }{L^2}}tk{\displaystyle \frac{4\pi ^2}{L^2}}tk^2)\})_{s=1}.`$
Podemos expresar este resultado en términos de la función theta de Jacobi, que queda definida por

$$\theta (x,y)=\underset{k=\mathrm{}}{\overset{\mathrm{}}{}}e^{2\pi kx}e^{\pi k^2y}.$$
(166)
Por lo tanto

$$\underset{k=\mathrm{}}{\overset{\mathrm{}}{}}\mathrm{exp}\left(\frac{4\pi \alpha }{L^2}tk\frac{4\pi ^2}{L^2}tk^2\right)=\theta (\frac{2\alpha t}{L^2},\frac{4\pi t}{L^2}).$$
(167)
Considerando la siguiente propiedad de inversión de la función theta de Jacobi

$$\theta (x,y)=\frac{1}{\sqrt{y}}\mathrm{exp}\left(\frac{\pi x^2}{y}\right)\theta (\frac{x}{iy},\frac{1}{y}),$$
(168)
resulta

$`{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\mathrm{exp}\left({\displaystyle \frac{4\pi \alpha }{L^2}}tk{\displaystyle \frac{4\pi ^2}{L^2}}tk^2\right)`$ $`=`$ $`\sqrt{{\displaystyle \frac{\pi }{t}}}\left({\displaystyle \frac{L}{2\pi }}\right)\mathrm{exp}\left({\displaystyle \frac{\alpha ^2t}{L^2}}\right)\theta ({\displaystyle \frac{i\alpha }{2\pi }},{\displaystyle \frac{L^2}{4\pi t}})`$ (169)
$`=`$ $`\sqrt{{\displaystyle \frac{\pi }{t}}}\left({\displaystyle \frac{L}{2\pi }}\right)\mathrm{exp}\left({\displaystyle \frac{\alpha ^2t}{L^2}}\right){\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}\mathrm{exp}\left(i\alpha k{\displaystyle \frac{L^2}{4t}}k^2\right).`$ (171)
Reemplazando en la expresión para la energía de vacío

$`E_0^L`$ $`=`$ $`{\displaystyle \frac{eBA}{2\pi }}\sqrt{\pi }{\displaystyle \frac{L}{2\pi }}(\mu ^{s+1}\{2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}[(m^2+2eBn)t]\times `$ (176)
$`\times {\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}e^{ik\alpha }\mathrm{exp}({\displaystyle \frac{L^2k^2}{4t}})`$
$`{\displaystyle \frac{1}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}(m^2t){\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}e^{ik\alpha }\mathrm{exp}({\displaystyle \frac{L^2k^2}{4t}})\})_{s=1}.`$
Nuevamente, calculamos la densidad de energía de vacío por unidad de volumen. Escribiendo por separado el término correspondiente a $`k=0`$ podemos identificar a la densidad de energía de vacío sin condiciones de contorno.

$`\epsilon _0^L`$ $`=`$ $`{\displaystyle \frac{eB}{(2\pi )^2}}(\mu ^{s+1}{\displaystyle \frac{\sqrt{\pi }}{\mathrm{\Gamma }(s/2)}}\{{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}(m^2t)`$ (183)
$`2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}[(m^2+2eBn)t]\})_{s=1}+`$
$`+{\displaystyle \frac{eB}{(2\pi )^2}}(\mu ^{s+1}{\displaystyle \frac{\sqrt{\pi }}{\mathrm{\Gamma }(s/2)}}\{{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}(m^2t)2{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }\mathrm{exp}({\displaystyle \frac{L^2k^2}{4t}})`$
$`2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}[(m^2+2eBn)t]2{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }\mathrm{exp}({\displaystyle \frac{L^2k^2}{4t}})\})_{s=1}.`$
Efectivamente, los dos primeros términos corresponden a la densidad de energía de vacío $`\epsilon _0`$ del problema en todo el espacio estudiado en la sección previa (ver ecuación (91))

$`\epsilon _0^L`$ $`=`$ $`\epsilon _0{\displaystyle \frac{2eB}{(2\pi )^2}}\{{\displaystyle \frac{\sqrt{\pi }\mu ^{s+1}}{\mathrm{\Gamma }(s/2)}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}dtt^{\frac{s1}{2}1}\mathrm{exp}(m^2t)\times `$ (186)
$`\times [2{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}(e^{2eBt})^n1]\mathrm{exp}({\displaystyle \frac{L^2k^2}{4t}})\}_{s=1}.`$
A fin de estudiar la variación de energía debida al tamaño finito, realizamos la suma geométrica, cuyo resultado es

$$\underset{n=0}{\overset{\mathrm{}}{}}x^n=\frac{1}{1x}\text{para}x<1,$$
(187)
y evaluamos para $`s=1`$

$$\epsilon _0^L=\epsilon _0+\frac{eB}{(2\pi )^2}\underset{k=1}{\overset{\mathrm{}}{}}e^{ik\alpha }_0^{\mathrm{}}𝑑tt^2\mathrm{exp}(m^2t)\mathrm{coth}(eBt)\mathrm{exp}\left(\frac{L^2k^2}{4t}\right).$$
(188)
Cabe señalar aquí que, para $`\alpha =0,\pi `$, este resultado coincide con el encontrado, mediante el esquema de regularización del tiempo propio en .

A fin de obtener una expresión más explícita para la densidad de energía de vacío, usamos el siguiente desarrollo para la cotangente hiperbólica

$$x\mathrm{coth}(x)=1+\frac{x^2}{3}2x^4\underset{q=1}{\overset{\mathrm{}}{}}\frac{1}{\pi ^2q^2(x^2+\pi ^2q^2)}.$$
(189)
La contribución a (188) del primer término del desarrollo es

$`=`$ $`{\displaystyle \frac{1}{(2\pi )^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑tt^3\mathrm{exp}\left(m^2t{\displaystyle \frac{L^2k^2}{4t}}\right)`$ (190)
$`=`$ $`{\displaystyle \frac{1}{(2\pi )^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑xx\mathrm{exp}\left({\displaystyle \frac{m^2}{x}}{\displaystyle \frac{L^2k^2x}{4}}\right).`$ (192)
La integral puede resolverse con el siguiente resultado , el cual es válido para $`Re(\beta )>0`$ y $`Re(\gamma )>0`$

$$_0^{\mathrm{}}𝑑xx^{\nu 1}\mathrm{exp}\left(\frac{\beta }{x}\gamma x\right)=2\left(\frac{\beta }{\gamma }\right)^{\nu /2}K_\nu \left(2\sqrt{\beta \gamma }\right),$$
(193)
donde $`K_\nu `$ son las funciones de Bessel modificadas. La contribución de este término, independiente del campo, es

$$\epsilon _0^c=\frac{2m^2}{\pi ^2L^2}\underset{k=1}{\overset{\mathrm{}}{}}\frac{e^{ik\alpha }}{k^2}K_2(mLk).$$
(194)
Este término se reduce, para $`\alpha =0`$, al resultado bien conocido para la energía de Casimir de un campo de Dirac en una caja con condiciones periódicas .

El segundo término del desarrollo (189) dará una contribución cuadrática en los campos

$`=`$ $`{\displaystyle \frac{1}{3}}{\displaystyle \frac{(eB)^2}{(2\pi )^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑tt^1\mathrm{exp}\left(m^2t{\displaystyle \frac{L^2k^2}{4t}}\right)`$ (195)
$`=`$ $`{\displaystyle \frac{2}{3}}\left({\displaystyle \frac{eB}{2\pi }}\right)^2{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }K_0(mLk).`$ (197)
donde nuevamente hemos usado (193) para resolver la integral.

Podemos analizar la convergencia de estos dos primeros términos a partir del comportamiento asintótico de las funciones de Bessel modificadas para argumentos grandes . Para $`z\mathrm{}`$

$$K_\nu (z)=\sqrt{\frac{\pi }{2z}}e^z\left[1+O\left(\frac{1}{z}\right)\right],$$
(198)
independientemente de $`\nu `$. Vemos que el factor $`e^z`$ garantiza la convergencia de las sumas (194) y (195)

Para evaluar el último término en (189) sumamos y restamos el término que correspondería a $`k=0`$ para poder extender la suma desde $`\mathrm{}`$ hasta $`+\mathrm{}`$ y poder reconstruir la función theta de Euler.

$`=`$ $`2{\displaystyle \frac{(eB)^4}{(2\pi )^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑tt{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\mathrm{exp}\left(m^2t\frac{L^2k^2}{4t}\right)}{\pi ^2q^2[(eBt)^2+\pi ^2q^2]}}`$ (199)
$`=`$ $`{\displaystyle \frac{(eB)^4}{(2\pi )^2}}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑tt{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\mathrm{exp}\left(m^2t\frac{L^2k^2}{4t}\right)}{\pi ^2q^2[(eBt)^2+\pi ^2q^2]}}+`$ (203)
$`+{\displaystyle \frac{(eB)^4}{(2\pi )^2}}{\displaystyle _0^{\mathrm{}}}𝑑tt{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\mathrm{exp}(m^2t)}{\pi ^2q^2[(eBt)^2+\pi ^2q^2]}}.`$
En el primer término podemos formar el siguiente factor y reescribirlo usando la definición de la función $`\theta (x,y)`$ y su propiedad de inversión (168)

$$\underset{k=\mathrm{}}{\overset{\mathrm{}}{}}e^{ik\alpha }e^{\frac{L^2k^2}{4t}}=\sqrt{\frac{t}{\pi }}\left(\frac{2\pi }{L}\right)\underset{k=\mathrm{}}{\overset{\mathrm{}}{}}\mathrm{exp}\left[\left(\frac{2k\pi +\alpha }{L}\right)^2t\right].$$
(204)
La contribución a la densidad de energía de vacío es

$`=`$ $`{\displaystyle \frac{(eB)^2}{2\pi ^{3/2}L}}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{\pi ^2q^2}}{\displaystyle _0^{\mathrm{}}}𝑑t.t^{3/2}{\displaystyle \frac{\mathrm{exp}\left[m^2t\left(\frac{2k\pi +\alpha }{L}\right)^2t\right]}{\left[t^2+\left(\frac{q\pi }{eB}\right)^2\right]}}+`$ (207)
$`+\left({\displaystyle \frac{eB}{2\pi }}\right)^2{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{\pi ^2q^2}}{\displaystyle _0^{\mathrm{}}}𝑑t.t{\displaystyle \frac{\mathrm{exp}(m^2t)}{\left[t^2+\left(\frac{q\pi }{eB}\right)^2\right]}}.`$
Las integrales las resolvemos usando la siguiente solución que es válida para $`Re(\beta )>0`$, $`Re(\mu )>0`$ y $`Re(\nu )>1`$ .

$`{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{x^\nu \mathrm{exp}(\mu x)}{x^2+\beta ^2}}`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Gamma }(\nu )}{2}}\beta ^{\nu 1}[\mathrm{exp}(i\mu \beta +i{\displaystyle \frac{(\nu 1)\pi }{2}})\mathrm{\Gamma }(1\nu ,i\mu \beta )+`$ (210)
$`+\mathrm{exp}(i\mu \beta i{\displaystyle \frac{(\nu 1)\pi }{2}})\mathrm{\Gamma }(1\nu ,i\mu \beta )],`$
donde $`\mathrm{\Gamma }(a,z)`$ es la función gamma incompleta .

La contribución a la densidad de energía es

$``$ $`{\displaystyle \frac{(eB)^{3/2}}{8\pi ^{5/2}L}}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^{3/2}}}\{\mathrm{exp}({\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]+{\displaystyle \frac{i\pi }{4}})\times `$ (213)
$`\times \mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{eB}}\left[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2\right])+`$
$`+`$ $`\mathrm{exp}({\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]{\displaystyle \frac{i\pi }{4}})\mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2])\}+`$ (215)
$`+`$ $`{\displaystyle \frac{(eB)^2}{8\pi ^4}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^2}}\left[\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{eB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{eB}})+\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{eB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{eB}})\right].`$ (217)
Para verificar la convergencia de las series, consideramos el comportamiento asintótico de las funciones gamma incompletas . Para $`z\mathrm{}`$ en $`|\mathrm{arg}(z)|<\frac{3\pi }{2}`$.

$$\mathrm{\Gamma }(a,z)z^{a1}e^z\left[1+\frac{a1}{z}+\frac{(a1)(a2)}{z^2}+\mathrm{}\right].$$
(218)
Para $`k`$ y/o $`q`$ suficientemente grandes, los sumandos se comportan de la siguiente manera: en ambos casos los primeros términos del comportamiento asintótico se cancelan entre los dos términos de cada sumando. Para el sumando de la primera suma (doble en este caso)

$$\frac{3}{q^4}.\frac{e^{i\pi /4}(eB)^{5/2}}{\left\{i\pi \left[m^2+\left(\frac{2k\pi +\alpha }{L}\right)^2\right]\right\}^{5/2}}.$$
(219)
Para el segundo sumando, que sólo se suma en $`q`$

$$\frac{(eB)^2}{\pi ^2m^2q^4}.$$
(220)
Concluimos que las sumas de estos dos últimos términos de la densidad de energía de vacío son convergentes. La expresión que finalmente damos para la densidad de energía de vacío del campo de fermiones en interacción con un campo magnético uniforme y con una condición de contorno en la dirección del campo es la siguiente

$`\epsilon _0^L`$ $`=`$ $`\epsilon _0+\epsilon _0^c+{\displaystyle \frac{2}{3}}\left({\displaystyle \frac{eB}{2\pi }}\right)^2{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }K_0(mLk)`$ (229)
$`{\displaystyle \frac{(eB)^{3/2}}{8\pi ^{5/2}L}}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^{3/2}}}\{\mathrm{exp}({\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]+{\displaystyle \frac{i\pi }{4}})\times `$
$`\times \mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{eB}}\left[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2\right])+`$
$`+\mathrm{exp}({\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]{\displaystyle \frac{i\pi }{4}})\mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2])\}+`$
$`+{\displaystyle \frac{(eB)^2}{8\pi ^4}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^2}}\left[\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{eB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{eB}})+\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{eB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{eB}})\right].`$
En esta expresión de la densidad de energía de vacío hemos aislado una contribución, $`\epsilon _0`$, independiente de las condiciones de contorno y otra, $`\epsilon _0^c`$, debida exclusivamente al tamaño finito e independiente del campo. Si definimos la densidad de energía de Casimir como la diferencia entre $`\epsilon _0^L`$ en (229) y su valor para $`B=0`$ y $`L\mathrm{}`$,obtenemos (usando $`\epsilon _0`$ de (114))

$`\epsilon _{Cas}`$ $`=`$ $`\epsilon _0^L\epsilon _0^l`$ (230)
$`=`$ $`\left[{\displaystyle \frac{\frac{(2eB)^2}{6}}{8\pi ^2(s+1)}}\right]_{s=1}+{\displaystyle \frac{1}{16\pi ^2}}\{[\psi (1)\psi (1/2)+1+\mathrm{ln}\left({\displaystyle \frac{m^2}{2eB}}\right)\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)]{\displaystyle \frac{(2eB)^2}{6}}+`$ (244)
$`+m^4[\mathrm{ln}\left({\displaystyle \frac{m^2}{2eB}}\right){\displaystyle \frac{1}{2}}]m^22eB\mathrm{ln}\left({\displaystyle \frac{m^2}{2eB}}\right)2(2eB)^2\zeta _H^{}(1,{\displaystyle \frac{m^2}{2eB}})\}+`$
$`+\epsilon _0^c+{\displaystyle \frac{2}{3}}\left({\displaystyle \frac{eB}{2\pi }}\right)^2{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }K_0(mLk)`$
$`{\displaystyle \frac{(eB)^{3/2}}{8\pi ^{5/2}L}}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^{3/2}}}\{\mathrm{exp}({\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]+{\displaystyle \frac{i\pi }{4}})\times `$
$`\times \mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{eB}}\left[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2\right])+`$
$`+\mathrm{exp}({\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]{\displaystyle \frac{i\pi }{4}})\mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{eB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2])\}+`$
$`+{\displaystyle \frac{(eB)^2}{8\pi ^4}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^2}}\left[\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{eB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{eB}})+\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{eB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{eB}})\right].`$
Nuevamente, los términos proporcionales a $`B^2`$ pueden eliminarse por renormalización de la carga eléctrica. En este caso, debido al tamaño finito, aparecerán nuevas contribuciones a la carga renormalizada (por ejemplo, el tercer término en (229)). Para ver que, en efecto, éste es el único término que agrega una nueva contribución a la carga renormalizada, la calculamos explícitamente a partir de su definición (126) y de (188). Para ello, debemos mostrar que la serie derivada en (188) es uniformemente convergente. El segundo término en (188) se escribe, al desarrollar la cotangente hiperbólica,

$$\frac{1}{(2\pi )^2}\underset{k=1}{\overset{\mathrm{}}{}}e^{ik\alpha }_0^{\mathrm{}}𝑑tt^3\left(1+2\underset{q=1}{\overset{\mathrm{}}{}}\frac{(eBt)^2}{(eBt)^2+\pi ^2q^2}\right)\mathrm{exp}\left(m^2t\frac{L^2k^2}{4t}\right).$$
(245)
La serie derivada resulta

$$\frac{2e^2}{(2\pi )^2}\underset{k=1}{\overset{\mathrm{}}{}}e^{ik\alpha }_0^{\mathrm{}}𝑑tt^1\underset{q=1}{\overset{\mathrm{}}{}}\frac{\pi ^2q^2}{[(eBt)^2+\pi ^2q^2]^2}\mathrm{exp}\left(m^2t\frac{L^2k^2}{4t}\right),$$
(246)
y el módulo de sus términos puede acotarse con los términos independientes de $`B`$ de una serie convergente

$`\left|{\displaystyle \frac{2e^2}{(2\pi )^2}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑tt^1{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\pi ^2q^2}{[(eBt)^2+\pi ^2q^2]^2}}\mathrm{exp}\left(m^2t{\displaystyle \frac{L^2k^2}{4t}}\right)\right|`$ (247)
(248)
$`{\displaystyle \frac{2e^2}{(2\pi )^2}}{\displaystyle \frac{\zeta _R(2)}{\pi ^2}}{\displaystyle _0^{\mathrm{}}}𝑑tt^1\mathrm{exp}\left(m^2t{\displaystyle \frac{L^2k^2}{4t}}\right)={\displaystyle \frac{e^2}{\pi ^4}}\zeta _R(2)K_0(mLk).`$ (249)
Concluimos que la serie derivada en (188) es uniformemente convergente y, por lo tanto, podemos evaluar su derivada término a término. Así, obtenemos

$`{\displaystyle \frac{}{B^2}}\left[{\displaystyle \frac{1}{(2\pi )^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }{\displaystyle _0^{\mathrm{}}}𝑑tt^3\left(1+2{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(eBt)^2}{(eBt)^2+\pi ^2q^2}}\right)\mathrm{exp}\left(m^2t{\displaystyle \frac{L^2k^2}{4t}}\right)\right]_{B^2=0}=`$ (250)
(251)
$`={\displaystyle \frac{e^2}{6\pi ^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }K_0(mLk).`$ (252)
De esta manera, se agrega un término en la carga renormalizada que depende del tamaño finito en la dirección del campo

$`{\displaystyle \frac{\mu ^4}{8e_R^2}}`$ $`=`$ $`{\displaystyle \frac{\mu ^4}{8e^2}}+{\displaystyle \frac{\mu ^4}{96\pi ^2}}\left[{\displaystyle \frac{2}{s+1}}+\psi (1)\psi (1/2)\mathrm{ln}c\right]_{s=1}+`$ (255)
$`+{\displaystyle \frac{\mu ^4}{24\pi ^2}}{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}e^{ik\alpha }K_0(mLk).`$
La expresión que finalmente nos queda para la densidad de energía de Casimir

$`\epsilon _{Cas}`$ $`=`$ $`{\displaystyle \frac{B^2}{2}}+{\displaystyle \frac{1}{16\pi ^2}}\{{\displaystyle \frac{(2e_RB)^2}{6}}[\mathrm{ln}\left({\displaystyle \frac{m^2}{2e_RB}}\right)+1]+m^4[\mathrm{ln}\left({\displaystyle \frac{m^2}{2e_RB}}\right){\displaystyle \frac{1}{2}}]`$ (266)
$`m^22e_RB\mathrm{ln}\left({\displaystyle \frac{m^2}{2e_RB}}\right)2(2e_RB)^2\zeta _H^{}(1,{\displaystyle \frac{m^2}{(2e_RB)}})\}+\epsilon _0^c`$
$`{\displaystyle \frac{(e_RB)^{3/2}}{8\pi ^{5/2}L}}{\displaystyle \underset{k=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^{3/2}}}\{\mathrm{exp}({\displaystyle \frac{iq\pi }{e_RB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]+{\displaystyle \frac{i\pi }{4}})\times `$
$`\times \mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{e_RB}}\left[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2\right])+`$
$`+\mathrm{exp}({\displaystyle \frac{iq\pi }{e_RB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2]{\displaystyle \frac{i\pi }{4}})\mathrm{\Gamma }(1/2,{\displaystyle \frac{iq\pi }{e_RB}}[m^2+\left({\displaystyle \frac{2k\pi +\alpha }{L}}\right)^2])\}+`$
$`+{\displaystyle \frac{(e_RB)^2}{8\pi ^4}}{\displaystyle \underset{q=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{q^2}}\left[\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{e_RB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{e_RB}})+\mathrm{exp}\left({\displaystyle \frac{im^2q\pi }{e_RB}}\right)\mathrm{\Gamma }(0,{\displaystyle \frac{im^2q\pi }{e_RB}})\right].`$
Otra vez debe destacarse que, una vez renormalizada la carga, la densidad de energía de Casimir no depende del parámetro de escala $`\mu `$.

## VIII Resumen de los resultados obtenidos

En este trabajo, hemos calculado las correcciones cuánticas a la energía de vacío (energías de Casimir) para campos de Dirac en presencia de un campo magnético uniforme de background $`\stackrel{}{B}`$. Hemos analizado las correcciones debidas exclusivamente a la presencia del campo de background y al efecto combinado de dicha presencia y la imposición de condiciones de contorno. En todos los casos hemos utilizado el método de suma de modos, combinado con la regularización zeta .

En primer lugar, hemos evaluado la energía de Casimir para campos de Dirac confinados en un plano perpendicular a la dirección de $`\stackrel{}{B}`$ (problema en $`2+1`$ dimensiones). En este caso, la energía de Casimir, definida por extensión analítica, resulta finita y no requiere, por lo tanto, una renormalización de los parámetros de la teoría.

A continuación, hemos tratado el problema en $`3+1`$ dimensiones, cuando todo el espacio tridimensional es accesible a los fermiones. En este caso, la corrección cuántica a la energía de vacío resulta divergente, con residuos en los polos de dos tipos: independientes de $`B`$ y proporcionales a $`B^2`$. La primera de estas divergencias desaparece en el cálculo de la energía de Casimir, dado que está presente también en ausencia del campo magnético. Los términos proporcionales a $`B^2`$ pueden eliminarse por renormalización de la carga eléctrica.

Finalmente, hemos tratado el efecto combinado de campo de background y condiciones de contorno no locales (”twisted”), cuando el espacio accesible a los fermiones tiene tamaño finito en la dirección de $`\stackrel{}{B}`$. En esta situación, hemos podido separar las contribuciones a la energía de vacío debidas a uno y otro factor externo exclusivamente de aquella debida a la acción combinada de ambos. Aquí, nuevamente, debimos renormalizar la carga eléctrica a fin de obtener un resultado finito para la energía de Casimir

En los casos sin borde, hemos verificado la consistencia con resultados previos de otros autores  para la acción efectiva (ya que para campo eléctrico nulo debe coincidir con la energía de vacío). En el caso con condiciones ”twisted” hemos mostrado que, en los límites particulares de condiciones periódicas y antiperiódicas, nuestra energía de vacío, expresada en forma integral, coincide con resultados anteriores . En este trabajo, hemos encontrado, para un parámetro de ”periodicidad” arbitrario, una expresión más explícita de la energía de Casimir, como serie de funciones especiales.

Agradecimientos: Quisiera expresar mi agradecimiento a Horacio Falomir, a Gabriela Beneventano, a Pablo Pisani y a Karin Révora por sus valiosos comentarios y sugerencias. También quisiera agradecer a la Fundación Antorchas por el apoyo recibido, fundamental para la realización de este trabajo. Por último, agradezco a Mariel Santangelo la confianza depositada en mí, su inestimable dedicación y su guía durante estos últimos dos años.

## Apéndice A: Soluciones de la ecuación de Kummer

Las soluciones independientes de la ecuación de Kummer

$$z\frac{d^2w}{dz^2}+(bz)\frac{dw}{dz}aw=0$$
(267)
se pueden expresar en términos de

$$M(a,b,z)=1+\frac{az}{b}+\frac{a(a+1)z^2}{b(b+1)2!}+\mathrm{}+\frac{a(a+1)\mathrm{}(a+n1)z^n}{b(b+1)\mathrm{}(b+n1)n!}+\mathrm{},$$
(268)
definida sólo para $`b`$ distinto de un entero negativo o cero y

$$U(a,b,z)=\frac{\pi }{\mathrm{sin}\pi b}\left[\frac{M(a,b,z)}{\mathrm{\Gamma }(1+ab)\mathrm{\Gamma }(b)}z^{1b}\frac{M(1+ab,2b,z)}{\mathrm{\Gamma }(a)\mathrm{\Gamma }(2b)}\right],$$
(269)
si bien estas dos funciones son soluciones de la ecuación de Kummer, no siempre resultan linealmente independientes.

Para estudiar los distintos casos consideraremos las cuatro siguientes soluciones

$$\begin{array}{c}y_1(a,b,z)=M(a,b,z)\hfill \\ y_2(a,b,z)=z^{1b}M(1+ab,2b,z)\hfill \\ y_3(a,b,z)=U(a,b,z)\hfill \\ y_4(a,b,z)=e^zU(ba,b,z)\hfill \end{array}$$
(270)
Estudiaremos los comportamientos asintóticos de las soluciones en los distintos casos, para $`z\mathrm{}`$ y para $`z0`$, y veremos para que valores de $`a`$ y $`b`$ conducen a autofunciones de cuadrado integrable, recordando que en nuestro problema $`b`$ sólo toma valores enteros. En particular consideraremos el comportamiento de $`M(a,b,z)`$ para $`z`$ grande, cuando $`a`$ es distinto de un entero negativo,

$$M(a,b,z)=\frac{\mathrm{\Gamma }(b)}{\mathrm{\Gamma }(a)}e^zz^{ab}\left[1+O(|z|^1)\right]\text{para}z\mathrm{},Re(z)>0$$
(271)
y el comportamiento de $`U(a,b,z)`$ en un entorno del origen

$$U(a,b,z)=\{\begin{array}{cc}\frac{\mathrm{\Gamma }(b1)}{\mathrm{\Gamma }(a)}z^{1b}+O(|z|^{b2})& \text{si }b>2\\ & \\ \frac{1}{\mathrm{\Gamma }(a)}z^1+O(|\mathrm{ln}z|)& \text{si }b=2\\ & \\ \frac{1}{\mathrm{\Gamma }(a)}[\mathrm{ln}z+\psi (a)]+O(|z\mathrm{ln}z|)& \text{si }b=1\\ & \\ \frac{1}{\mathrm{\Gamma }(1+a)}+O(|z\mathrm{ln}z|)& \text{si }b=0\\ & \\ \frac{\mathrm{\Gamma }(1b)}{\mathrm{\Gamma }(1+ab)}+O(|z|)& \text{si }b<0\end{array}$$
(272)
1. *$`an`$ (negativo o cero) y $`b`$ entero $`b1`$*

La solución general es

$$w(z)=Ay_1(a,b,z)+By_3(a,b,z),$$
(273)
ya que su wronskiano es no nulo

$$W(y_1,y_3)=\frac{\mathrm{\Gamma }(b)}{\mathrm{\Gamma }(a)}z^be^z$$
(274)
Si $`z\mathrm{}`$, entonces $`y_1e^zz^{ab}`$ y, por lo tanto, las autofunciones que se obtienen a partir de esta solución divergen exponencialmente.

Si $`z0,y_3z^{1b}`$ para $`b2`$ y entonces las autofunciones, que divergen como $`\rho ^k`$ para $`k1`$ en el origen, no son de cuadrado integrable. Mientras que para $`b=1`$ la divergencia es logarítmica y en consecuencia esta componente del espinor sí es de cuadrado integrable. Sin embargo, al derivar la componente restante de (19) la divergencia en el origen va como $`k^1`$.
2. *$`an`$ (negativo o cero) y $`b`$ entero $`b0`$*

En este caso, la solución general es

$$w(z)=Ay_2(a,b,z)+By_3(a,b,z),$$
(275)
dado que son soluciones linealmente independientes

$$W(y_2,y_3)=\frac{\mathrm{\Gamma }(2b)}{\mathrm{\Gamma }(1+ab)}z^be^z$$
(276)
Si $`z\mathrm{}`$, entonces $`y_2e^zz^{ab}`$ y ,por lo tanto, las autofunciones que se obtienen a partir de esta solución divergen exponencialmente.

Si $`z0,y_21`$ y entonces las autofunciones divergen como $`\rho ^k`$ con $`k1`$ en el origen y no son de cuadrado integrable.
3. *$`a=n`$ (negativo o cero) y $`b`$ entero $`b1`$*

Escribimos la solución general como

$$w(z)=Ay_1(a,b,z)+By_4(a,b,z),$$
(277)
para las cuales

$$W(y_1,y_4)=\frac{\mathrm{\Gamma }(b)}{\mathrm{\Gamma }(ba)}e^{iϵ\pi b}z^be^z$$
(278)
Entonces $`y_1=M(n,b,z)=\frac{n!}{(b+1)(b+2)\mathrm{}(b1+n)}L_n^{b1}`$, polinomio de Laguerre generalizado que conduce a autofunciones de cuadrado integrable.

Mientras que $`y_4=e^zU(b+n,b,z)e^z(z)^{bn}`$ diverge exponencialmente, para $`z\mathrm{}`$.
4. *$`a=n`$ (negativo o cero) y $`b`$ entero $`n+1b0`$*

La solución general es ahora

$$w(z)=Ay_2(a,b,z)+By_4(a,b,z),$$
(279)
y el wronskiano de estas soluciones es

$$W(y_2,y_4)=\frac{\mathrm{\Gamma }(2b)}{\mathrm{\Gamma }(1a)}z^be^z$$
(280)
Tenemos $`y_2=z^{1b}M(n+1b,2b,z)=z^{1b}\frac{(n+b1)!}{(2b)(2b+1)\mathrm{}(2b+n1)}L_{n+b1}^{b+1}`$, con las que resultan autofunciones de cuadrado integrable.

Por el contrario, las soluciones $`y_4=e^zU(b+n,b,z)e^z(z)^{bn}`$ divergen exponencialmente,para $`z\mathrm{}`$.
5. *$`a=n`$ (negativo o cero) y $`b`$ entero $`bn`$*

Para la solución general nuevamente proponemos la combinación

$$w(z)=Ay_2(a,b,z)+By_3(a,b,z),$$
(281)
ya que wronskiano de estas soluciones es

$$W(y_2,y_3)=\frac{\mathrm{\Gamma }(2b)}{\mathrm{\Gamma }(1+ab)}z^be^z$$
(282)
Si $`z\mathrm{}`$, entonces $`y_2e^zz^{ab}`$ y, por lo tanto, las autofunciones que se obtienen a partir de esta solución divergen exponencialmente.

Si $`z0,y_21`$ y entonces las autofunciones divergen como $`\rho ^k`$ con $`k1`$ en el origen y no son de cuadrado integrable.

## Apéndice B: Campo libre en $`2+1`$ dimensiones

La energía de vacío del campo libre se obtiene de integrar las energías de la partícula libre en el espacio de impulsos, multiplicadas por la densidad de estados de la partícula libre. Esto es

$$\epsilon _0^l=\frac{E_0}{A}=22\frac{1}{2}\frac{1}{(2\pi )^2}d^2p\left(m^2+p^2\right)^{\frac{1}{2}}.$$
(283)
Emplearemos el método de la función zeta para regularizar esta integral, el mismo que utilizamos en la sección IV. El parámetro $`\mu `$ con unidades de masa asegura que la expresión tenga unidades de densidad de energía por unidad de área.

$$\epsilon _0^l=\frac{2}{(2\pi )^2}\left[\mu ^{s+1}_{\mathrm{}}^{\mathrm{}}𝑑p_x_{\mathrm{}}^{\mathrm{}}𝑑p_y(m^2+p_x^2+p_y^2)^{s/2}\right]_{s=1}.$$
(284)
Usamos la transformada de Mellin (81)

$$\epsilon _0^l=\frac{2}{(2\pi )^2}\left[\frac{\mu ^{s+1}}{\mathrm{\Gamma }(s/2)}_0^{\mathrm{}}𝑑tt^{s/21}e^{m^2t}_{\mathrm{}}^{\mathrm{}}𝑑p_xe^{p_x^2t}_{\mathrm{}}^{\mathrm{}}𝑑p_ye^{p_y^2t}\right]_{s=1},$$
(285)
y evaluando las integrales de los impulsos con (88)

$`\epsilon _0^l`$ $`=`$ $`{\displaystyle \frac{2}{(2\pi )^2}}\left[{\displaystyle \frac{\mu ^{s+1}}{\mathrm{\Gamma }(s/2)}}\pi {\displaystyle _0^{\mathrm{}}}𝑑tt^{s/211}e^{m^2t}\right]_{s=1}`$ (286)
$`=`$ $`{\displaystyle \frac{1}{(2\pi )}}\left[\mu ^{s+1}{\displaystyle \frac{\mathrm{\Gamma }(s/21)}{\mathrm{\Gamma }(s/2)}}(m^2)^{(s/21)}\right]_{s=1}.`$ (288)
Esta expresión para la densidad de energía de vacío está regularizada ya que la función gamma sólo presenta singularidades en los enteros negativos y cero.

$$\epsilon _0^l=\frac{1}{(2\pi )}\frac{\mathrm{\Gamma }\left(3/2\right)}{\mathrm{\Gamma }\left(1/2\right)}m^3=\frac{m^3}{3\pi },$$
(289)
y el resultado coincide con el término independiente del campo del desarrollo (62).

## Apéndice C: Campo libre en $`3+1`$ dimensiones

De igual modo, expresamos la energía de vacío del problema de campo libre en $`3+1`$ dimensiones,

$$\epsilon _0^l=\frac{E_0}{L_xL_yL_z}=22\frac{1}{2}\frac{1}{(2\pi )^3}d^3p\left(m^2+p^2\right)^{\frac{1}{2}}.$$
(290)
Esta expresión es divergente. La regularizamos nuevamente por medio del método de la función zeta, e introducimos un parámetro para que conserve las unidades apropiadas

$$\epsilon _0^l=\frac{2}{(2\pi )^3}\left[\mu ^{s+1}_{\mathrm{}}^{\mathrm{}}𝑑p_x_{\mathrm{}}^{\mathrm{}}𝑑p_y_{\mathrm{}}^{\mathrm{}}𝑑p_z(m^2+p_x^2+p_y^2+p_z^2)^{s/2}\right]_{s=1}.$$
(291)
Hacemos uso de la transformada de Mellin (81) y obtenemos

$`\epsilon _0^l`$ $`=`$ $`{\displaystyle \frac{2}{(2\pi )^3}}[{\displaystyle \frac{\mu ^{s+1}}{\mathrm{\Gamma }(s/2)}}{\displaystyle _0^{\mathrm{}}}dtt^{s/21}e^{m^2t}\times `$ (294)
$`\times {\displaystyle _{\mathrm{}}^{\mathrm{}}}dp_xe^{p_x^2t}{\displaystyle _{\mathrm{}}^{\mathrm{}}}dp_ye^{p_y^2t}{\displaystyle _{\mathrm{}}^{\mathrm{}}}dp_ze^{p_z^2t}]_{s=1}.`$
Resolvemos las integrales con (88) y resulta

$`\epsilon _0^l`$ $`=`$ $`{\displaystyle \frac{2}{(2\pi )^3}}\left[{\displaystyle \frac{\mu ^{s+1}}{\mathrm{\Gamma }(s/2)}}\pi ^{3/2}{\displaystyle _0^{\mathrm{}}}𝑑tt^{\frac{s3}{2}1}e^{m^2t}\right]_{s=1}`$ (295)
$`=`$ $`{\displaystyle \frac{1}{4\pi ^{3/2}}}\left[\mu ^{s+1}{\displaystyle \frac{\mathrm{\Gamma }(\frac{s3}{2})}{\mathrm{\Gamma }(s/2)}}(m^2)^{\frac{s3}{2}}\right]_{s=1}.`$ (297)
Nuevamente usamos el parámetro adimensional $`c`$ definido por (54), y factorizamos la singularidad que la función $`\mathrm{\Gamma }(\frac{s3}{2})`$ tiene en $`s=1`$,

$$\epsilon _0^l=\frac{2\mu ^4}{\pi ^{3/2}}\left[\frac{\mathrm{\Gamma }(\frac{s+3}{2})}{\mathrm{\Gamma }(s/2)(s3)(s1)(s+1)}c^{\frac{s3}{2}}\right]_{s=1},$$
(298)
y separamos en partes singular y finita haciendo un desarrollo alrededor de la singularidad

$$\epsilon _0^l=\frac{2\mu ^4}{\pi ^{3/2}}\left[\frac{Z(1)}{s+1}+Z^{}(1)\right]_{s=1},$$
(299)
donde la función $`Z(s)`$ queda definida por

$$Z(s)=\frac{\mathrm{\Gamma }(\frac{s+3}{2})c^{\frac{s3}{2}}}{\mathrm{\Gamma }(s/2)(s3)(s1)}.$$
(300)
Para la parte singular

$$Z(1)=\frac{\mathrm{\Gamma }(1)c^2}{\mathrm{\Gamma }(1/2)(4)(2)}=\frac{c^2}{16\sqrt{\pi }}$$
(301)
Para el término finito

$$Z^{}(1)=\frac{c^2}{32\sqrt{\pi }}\left[\psi (1)\psi (1/2)+\frac{3}{2}\mathrm{ln}c\right].$$
(302)
Finalmente, para la densidad de energía de vacío del campo libre en $`3+1`$ dimensiones obtenemos

$$\epsilon _0^l=\frac{m^4}{8\pi ^2(s+1)}+\frac{m^4}{16\pi ^2}\left[\psi (1)\psi (1/2)+\frac{3}{2}\mathrm{ln}\frac{m^2}{\mu ^2}\right],$$
(303)
que coincide con el resultado (119) obtenido en la sección V.