\documentclass[a4paper,12pt]{article}
\usepackage{a4wide}

%%%%%%%%%%%%%%%%%%%%%          MACROS FILE      %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%    GREEK ALPHABET  

\def\al{\alpha}
\def\be{\beta}
\def\ga{\gamma}
%\def\de{\delta}
\def\ep{\epsilon}
\def\ze{\zeta}
\def\io{\iota}
\def\ka{\kappa}
\def\la{\lambda}
\def\rh{\rho}
\def\si{\sigma}
\def\om{\omega}
\def\ph{\phi}
\def\th{\theta}
\def\up{\upsilon}
\def\Ga{\Gamma}
\def\De{\Delta}
\def\La{\Lambda}
\def\Si{\Sigma}
\def\Om{\Omega}
\def\Th{\Theta}
\def\Up{\Upsilon}

%    LETTERE E SIMBOLI FISICI


\newcommand{\omn}{\omega_{n}}
\newcommand{\omnd}{\omega_{n}^{2}}
\newcommand{\Omn}{\Omega_{n}}
\newcommand{\Omnd}{\Omega_{n}^{2}}
\newcommand{\mud}{\mu^{2}}
\newcommand{\nablad}{\nabla^{2}}
\newcommand{\piop}{\hat{\pi}}
\newcommand{\phiop}{\hat{\phi}}
\newcommand{\phid}{\phi^{2}}
\newcommand{\phiq}{\phi^{4}}
\newcommand{\phiz}{\phi_{0}}
%\newcommand{\phibar}{\bar{\phi}}
\newcommand{\phibar}{\stackrel{\,\_}{\phi}}
\newcommand{\bra}{\left\langle}
\newcommand{\ket}{\right\rangle} 
\newcommand{\stato}[1]{|\phi_{#1}\!\!>}

%    INTEGRALI E SOMME

\newcommand{\intdtaudx}{\int_{0}^{\beta}d\tau \int d^{3}x }
\newcommand{\intk}{\int\frac{d^{3}k}{(2\pi)^{3}} }
\newcommand{\intnk}{\sum_{n}\int\frac{d^{3}k}{(2\pi)^{3}} }
\newcommand{\sumn}{\sum_{n}}
\newcommand{\mellin}[1]{\frac{1}{\Gamma (s)}\int_{0}^{\infty} t^{#1}}
\newcommand{\dfx}{ {(2\pi)^{-4}} \int_{-\infty}^{\infty} {d^4x}\;}

%    DERIVATE

\newcommand{\de}[1]{{{\partial \over \partial #1}}}
\newcommand{\dde}[2]{{{\partial #1 \over \partial #2}}}
\newcommand{\detau}{\partial_{\tau}}
\newcommand{\detaud}{\partial_{\tau}^{2}}
\newcommand{\deindetau}[1]{\frac{\partial #1}{\partial\tau} }
\newcommand{\deindetaud}[1]{\frac{\partial^{2} #1}{\partial\tau^{2}} }

%    MATRICI
\newcommand{\matrice}[4]{\left(\begin{array}{cc} {#1} & {#2} \\ {#3} & {#4} \end{array} \right) }

%          MATH MODE

\newcommand{\half}{\frac{1}{2}}
\def\beq{\begin{eqnarray}}    %%%  begequation/eqnarray
\def\enq{\end{eqnarray}}      %%%  endequation/eqnarray
%\newcommand{\bea}{\begin{eqnarray}}
%\newcommand{\eea}{\end{eqnarray}}
%\newcommand{\beq}{\begin{equation}}
%\newcommand{\eeq}{\end{equation}}
%\newcommand{\beq}{\begin{displaymath}}
%\newcommand{\enq}{\end{displaymath}}
%\newcommand{\beqn}{\[}
%\newcommand{\eeqn}{\]}
\newcommand{\nnel}{\nonumber \\ {}}
 
%%%  calligraphic

\newcommand{\ca}[1]{{\cal #1}}        
\newcommand{\calB}{{\cal B}}
\newcommand{\calL}{{\cal L}}
\newcommand{\calV}{{\cal V}}
\newcommand{\calW}{{\cal W}}
\newcommand{\calZ}{{\cal Z}}
\newcommand{\Zmath}{{\cal Z}}

\newcommand{\jvec}{\vec{j}}
\newcommand{\phivec}{\vec{\phi}}
\renewcommand{\vec}[1]{{\bf #1}}       %%%  vectors in bold
\def\hs{\qquad}               %%%  horizontal space
\def\nn{\nonumber}            %%%  no number for eqnarray

%    brackets

\def\ap{\left.}               %%%  open bracket
\def\at{\left(}               %%%  open (
\def\aq{\left[}               %%%  open [
\def\ag{\left\{}              %%%  open {
\def\cp{\right.}              %%%  close bracket
\def\ct{\right)}              %%%  close )
\def\cq{\right]}              %%%  close ]
\def\cg{\right\}}             %%%  close }

%%%%%    SPECIAL SYMBOLS
\def\R{{\hbox{{\rm I}\kern-.2em\hbox{\rm R}}}}   %%% real numbers
\def\H{{\hbox{{\rm I}\kern-.2em\hbox{\rm H}}}}   %%% Hilbert space
\def\N{{\hbox{{\rm I}\kern-.2em\hbox{\rm N}}}}   %%% natural numbers
\def\C{{\ \hbox{{\rm I}\kern-.6em\hbox{\bf C}}}} %%% complex numbers
\def\Z{{\hbox{{\rm Z}\kern-.4em\hbox{\rm Z}}}}   %%% integers numbers
\def\ii{\infty}                                  %%% infinit
\def\X{\times\,}                                 %%% times
\newcommand{\fr}[2]{\mbox{$\frac{#1}{#2}$}}      %%% small fraction
\def\Det{\mathop{\rm Det}\nolimits}              %%% Determinant
\def\tr{\mathop{\rm tr}\nolimits}                %%% trace
\def\Tr{\mbox{Tr}\ }                             %%% Trace
%\newcommand{\Tr}{{\rm Tr}}

\def\PP{\mathop{\rm PP}\nolimits}                %%% Finite part
\def\Res{\mathop{\rm Res}\nolimits}              %%% Residue
\renewcommand{\Re}{\mathop{\rm Re}\nolimits}     %%% Real 
\renewcommand{\Im}{\mathop{\rm Im}\nolimits}     %%% Imaginary
\def\lap{\Delta}                                 %%% Laplacian
\begin{document}
\begin{center}

{\LARGE \bf New physics in the} \\ 
{\LARGE \bf charged relativistic Bose gas} \\
{\LARGE \bf using zeta-function regularization?\footnote{Talk given at the ``SEWM2000 - Strong and
Electro-Weak Matter 2000'' Conference, June 2000,
Marseilles, France. Proceedings to be published by World Scientific.}}
\\[4mm]
{\large Antonio Filippi}
\\[4mm]
{\it Theoretical Physics Group, Imperial College, \\ Prince Consort
Road, London SW7 2BZ, United Kingdom\\E-mail: a.filippi@ic.ac.uk}
\\[4mm]
\end{center}

\abstract{The multiplicative anomaly, recently introduced in QFT,
plays a fundamental role in solving some mathematical inconsistencies of the widely
used zeta-function regularization method. Its physical relevance is
still an open question and is here analyzed in the light of a
non-perturbative method. Even in this approach the ``different
physics'' seems to hold and not to be easily removable by renormalization.}

\section{Introduction}

The evaluation of functional determinants of pseudo-differential
operators is often a fundamental point in quantum field theory
calculations.
As these operators are unbounded the corresponding determinant is
undefined, unless a rigorous regularization procedure is applied.
One of the most successful is the zeta-function regularization method \cite{dowcri}-\cite{bcvz}.
This uses the so called ``generalized Riemann zeta function''  $\zeta(s|A)=\mbox{Tr}\, A^{-s}$,
(A a pseudo-differential operator) which is well defined for a sufficiently large real part of $s$ and can be analytically
continued to a function meromorphic in all the plane and analytic at
$s=0$. As such its derivative with respect to $s$ at zero is well defined and
the logarithm of the zeta-function regularized functional determinant
will then be defined by
\begin{eqnarray}
\ln\det\frac{A}{\sigma^2} =-\zeta'(0|A)-\zeta(0|A)\ln \sigma^2 ,
\end{eqnarray}
where $\sigma^2$ is a renormalization scale mass.

The recent introduction \cite{evz}-\cite{regensburg} of the
long-overlooked 
multiplicative anomaly \cite{wod}-\cite{konvis}
in this framework has provoked a lively debate
\cite{byt}-\cite{mcktom2}.
Its mathematical rigor and importance for
quantum field theory has been widely
proved by this and other authors \cite{mcktom} but the question of its physical
relevance is still open.

The multiplicative anomaly can appear when we evaluate a quantity of the
form $\ln \det (A B)$, with $A$ and $B$ two commuting pseudo-differential operators.
The fact is that it is not always true that the equality $\ln \det
(AB)=\ln \det(A)+\ln \det (B)$ holds, as normally assumed in
quantum field theory until recently.
On the contrary, using zeta-function regularization, an additional term $a(A,B)$,
called the `multiplicative anomaly,
may be present on the right hand-side. The anomaly can rarely be computed
directly as the difference of the two sides. Fortunately Wodzicki.
produced a very useful indirect formula for its evaluation (see\cite{wod}).

The fact that the anomaly can add an
``anomalous'' term to the ``standard'' result, the controversial physical relevance
of which is still under scrutiny, will be the subject of this work.

For some of the very simple models studied up to now
the multiplicative anomaly is not always present,
while in others its ``anomalous term'' can be easily renormalized away
at one-loop
by a finite shift in the counterterms \cite{evz}.
Renormalizability is therefore a fundamental issue to analyze to
understand the physical content of the multiplicative anomaly.

To this end I will need to go beyond the one-loop approximation
used until now.
Since we finally want to analyze the symmetry
breaking transition, the
appropriate candidate seems a self-consistent approach similar to the
large $N$ expansion \cite{raybook}.

The relativistic charged Bose gas at finite temperature has long been of
interest on its own\cite{kap1}-\cite{chemical}. The large N limit
of the $0(N)$ interacting field can be found in Haber and
Weldon. None of these remarkable papers adopts a fully regularized
approach as there are
formal manipulations of infinite sums involved.
We showed how the anomaly is crucial for
the consistency of the zeta-function method and how it creates an
unexpected ``anomalous term''. I refer the reader to the relevant
papers \cite{efvz}-\cite{regensburg} for details and results.

\section{The non-perturbative approach}

The inclusion of the chemical potential $\mu$ in a Hamiltonian representation
of the grand-canonical partition function is usually (see \cite{kap1,habwel3,bbd,chemical})
realized by defining an effective Hamiltonian $H=H_o + \mu Q$ where $Q$ is the
charge density operator. Through integration on the momenta the
functional integral can be cast in Lagrangian form.
Although I am
only interested in $N=2$ I will leave the $N$ in the interacting term
for clarity, its sum over
repeated indices $a=1,2$ understood. The starting action is therefore
\begin{equation}
S[\phi]\!\!=\!\!\int_0^\beta\!\!\! d\tau\int\!\!\! d^3 x
\!\left(\partial_{\nu}\phi_{a}\partial^{\nu}\phi_{a}\!\!+\!m^2_{o}\phi^2\!\!+\!\frac{\lambda_{o}}{4N}\phi^4\!\!-\!\mu^2\phi^2\!\!+\!2i\mu(\phi_{2}\partial_{\tau}\phi_{1}\!\!-\!\phi_{1}\partial_{\tau}\phi_{2})\right)
\end{equation}
I can then insert in the functional integral a Gaussian representation
of unity, through an auxiliary field $B(x)$ \cite{raybook}.
Considering constant sources $J_a$ for the $\phi_{a}$ fields,
 the partition function then  becomes
\begin{equation}
Z[J]=\int{\cal D}B{\cal D}\phi_{1}{\cal D}\phi_{2}e^{-\frac{1}{2} \int dx [\phi_{a}(x)A_{ab}(x)\phi_{b}(x)-\frac{N}{\lambda_{o}}B^2(x)]}
\end{equation}
where I defined the matrix valued differential operator
\begin{equation}
A(x)=\matrice{-\partial^{2}_{\tau}-\nabla^2+m^2_{o}+B(x)-\mu^{2}}{-2i\mu\partial_{\tau}}{2i\mu\partial_{\tau}}{-\partial^{2}_{\tau}-\nabla^2+m^2_{o}+B(x)-\mu^{2}} \label{aop}
\end{equation}
 The functional integral in $\phi_a$ can be recast in exponential form to
contribute in an effective action for the $B(x)$ field.
It is also clear that this effective action is of order $N$.
We can therefore apply a saddle point approximation in the field $B(x)$,
such that

\begin{equation}
Z[J]\simeq\int{\cal D}\phi_{1}{\cal D}\phi_{2}e^{-\frac{1}{2} S[\phi,\bar{B}]+ \frac{1}{\beta}J_{a} \int dx \phi_{a}(x)}
\end{equation}
where constant $\bar{B}$ is the large-N saddle-point.

At this point it proves useful to turn to the momentum
representation,
\begin{equation}
Z[J]=e^{\frac{\beta V N}{2\lambda_o}\bar{B}^2+\frac{V}{2\beta}J_{a}{\cal
A}_{0ab}^{-1}J_{b}}e^{-\frac{1}{2}\log\det(A\sigma^2)}
\end{equation}
where I denote the eigenvalues matrix  as ${\cal A}(n,k)$ and  ${\cal
A}(0,0)$ as ${\cal A}_0$. Then
\begin{equation}
W=-\frac{1}{\beta V}\log Z = -\frac{N}{2\lambda_o}\bar{B}^2-\frac{1}{2\beta^2}J_{a}{\cal
A}_{0ab}^{-1}J_{b}+\frac{1}{2\beta V}\log\det(A\sigma^2)
\end{equation}
where ${\bar B}$ is
\begin{equation}
\frac{ N}{\lambda_o}\bar{B}\!+\!\frac{1}{\beta V}\frac{\partial}{\partial
\bar{B}}\left(\!-\!\frac{1}{2}\log\det(A\sigma^2)\right)\!-\!\frac{1}{2\beta^2
}\left(J_{1}({\cal A}_{011}^{-1})^2 J_{1}\!+\!J_{2}({\cal A}_{022}^{-1})^2 J_{2}\right) \!\!=\!\! 0
\end{equation}

Going back to (\ref{aop}) we can see that $\bar{B}$ gives a
contribution to the mass of the field $\phi$, $ M^2=m_o^2 + \bar{B}
\label{effmass} $.
With this substitution the operator $A$ has
the same form as the one for the non-interacting field.
Calculations are then highly simplified and we can resort to previous
results in \cite{efvz}-\cite{regensburg}, for which zeta-function 
regularization was used and in which the multiplicative
anomaly contribution is properly accounted for:
\begin{equation}
-\!\frac{1}{2}\log\det(A\sigma^2)\!\!=\!\!-\!\frac{\beta
V}{32\pi^2}M^4\left(\log(M^2\sigma^2)\!-\!\frac{3}{2}\right)\!+\!S
\underbrace{\!-\!\frac{\beta
V}{16\pi^2}\mu^2\left(M^2\!-\!\frac{\mu^2}{3}\right)} \label{vecchio}
\end{equation}
where $S$ is the standard particles-antiparticles thermal term,
and the underlined part is the one resulting from the multiplicative anomaly.

The constant semi-classical fields are 
\begin{equation}
\bar{\phi}_1=\beta\frac{\partial W}{\partial
J_1}=-\frac{1}{\beta}{\cal A}_{011}^{-1} J_{1}\ \ \ \ \ \ \ \ \ \bar{\phi}_2=\beta\frac{\partial W}{\partial J_2}=-\frac{1}{\beta}{\cal A}_{022}^{-1} J_{2}
\label{semiclfield}
\end{equation} 
so that, using (\ref{vecchio}), we can express $M^2$ in terms of $\bar{\phi^2}$ as
\begin{equation}
M^2 \!\!\!=\!\!\! m_o^2\!+\!\frac{1}{16\pi^2}\frac{\lambda_o}{N}M^2\log(M^2\sigma^2)-\!\frac{\lambda_o}{N\beta
V}\frac{\partial S}{\partial M^2}\!+\!\frac{\lambda_o}{2N}\bar{\phi}_a\bar{\phi}_a
\!+\underbrace{\!\frac{1}{16\pi^2}\frac{\lambda_o}{N}\mu^2} \label{mquadrin}
\end{equation}
and
\beq  
\Gamma[\bar{\phi}]&=&M^2\left[\frac{N}{\lambda_o}m_o^2+\frac{1}{2}\bar{\phi}^2\right]+M^4\left[\frac{1}{32\pi^2}\log(M^2\sigma^2)-\frac{1}{64\pi^2}-\frac{N}{2\lambda_o}\right] \\ \nonumber 
&-&\frac{1}{\beta V}S -\frac{1}{2}\mu^2\bar{\phi}^2-\frac{N}{2\lambda_o}m^4 \underbrace{+\frac{1}{16\pi^2}\mu^2\left(M^2-\frac{\mu^2}{3}\right)}
\enq
where we have redefined the
mass scale $\sigma^2\rightarrow\sigma^2 e$
so as to bring our notation into correspondence with that of Haber and Weldon.
It is now time to put $N=2$ for good.
Performing standard renormalization techniques it seems
that the anomaly cannot be renormalised away. This aspect needs to be
further analyzed, and details will be given elsewhere.


Explicit calculation shows that  the unbroken phase charge density is
\begin{equation}
\rho=-\underbrace{\frac{\mu
M^2}{8\pi^2}+\frac{\mu^3}{12\pi^2}}+\frac{1}{\beta V}\frac{\partial
S}{\partial\mu}
\end{equation}

For the broken phase where, after the
phase transition, $\mu=M$,  we can find the expression for
$\bar{\phi}^2$ and the charge density takes the form
\beq
\rho&\!=\!&\!\mu\frac{4}{\lambda_o}\left[\mu^2\!-\!
m_o^2\!-\!\frac{\lambda_o}{32\pi^2}\mu^2\log(\mu^2\sigma^2)\!+\!\frac{\lambda_o}{2\beta
V}\left.\frac{\partial S}{\partial M^2}\right|_{M^2=\mu^2}\!\!\!
\!-\!\underbrace{\!\frac{\lambda_o}{32\pi^2}\mu^2}\right] \\ \nonumber
&+&\frac{1}{\beta V} \left.\frac{\partial S}{\partial \mu}\right|_{M^2=\mu^2}-\underbrace{\frac{1}{24\pi^2}\mu^3}
\enq
It seems like the anomaly could change the ``standard'' transition
temperature and, for fixed charge $\rho$, we can see a 
difference in physical observables (e.g. pressure) in the two phases.

\section{Comments}

The relevance of zeta-function regularization in QFT can not be easily
dismissed. It is as reliable regularization technique as others.
The multiplicative anomaly is
indisputably necessary for mathematical consistency. These results
seem to show that the ``anomalous term'' it creates in certain
conditions is not trivially removable and could play a role in the
physical measurables of the system. 
To my knowledge this is the first extension of the zeta-function
regularization method in a non-perturbative calculation, within a
neat procedure that could be useful on other models as well.
It is clearly a non-trivial problem, and also
very interesting as it goes to the roots of one of the most used
regularization methods and of the functional integral approach itself.

\section*{Acknowledgments}
This work has been developed in collaboration with R. Rivers, Imperial
College. My thanks also go to E. Elizalde, L. Vanzo, S. Zerbini and T. Evans
for stimulating discussions.
The author wishes to acknowledge financial support from the European Commission
under TMR contract N. ERBFMBICT972020 and
the Foundation Blanceflor Boncompagni-Ludovisi, n\'ee Bildt.

\begin{thebibliography}{99}

\bibitem{dowcri} 
J.S.~Dowker and R.~Critchley,  
Phys.~Rev. {\bf D 13}, 3224 (1976).

\bibitem{haw} 
S.W.~Hawking,  
Commun.~Math.~Phys. {\bf 55}, 133 (1977).

\bibitem{book} 
E.~Elizalde, S.~D.~Odintsov, A.~Romeo, A.A.~Bytsenko and S.~Zerbini, 
{\em Zeta Regularization Techniques with Applications},
World Scientific, Singapore (1994).

\bibitem{eliz} E.~Elizalde, {\em Ten Physical Applications of Spectral
Zeta Functions},
Springer-Verlag, New York (1995).

\bibitem{bcvz} 
A.A.~Bytsenko, G.~Cognola, L.~Vanzo and S.~Zerbini,  
Phys.~Rep. {\bf 266}, 1 (1996),
hep-th/9505061.

\bibitem{evz}
E.~Elizalde, L.~Vanzo and S.~Zerbini,
Commun.~Math.~Phys. {\bf 194}, 613-630 (1998),
hep-th/9701060.

\bibitem{efvz}
E.~Elizalde, A.~Filippi, L.~Vanzo and S.~Zerbini,
Phys.~Rev. {\bf D 57}, 7430-7443 (1998), 
hep-th/9710171.

\bibitem{bielefeld}
A.~Filippi,
Nucl. Phys. {\bf A 642}, 222-227 (1998),
hep-th/9809050.

\bibitem{regensburg}
A.~Filippi,
{\em Chemical potential and multiplicative anomaly}, 
hep-th/9809098.

\bibitem{wod} 
M.~Wodzicki, 
{\em Non-commutative Residue, Chapter I}.
In {\em Lecture notes in Mathematics}. Yu.I. Manin, editor, volume 1289, 320.   
Springer-Verlag  (1987).

\bibitem{kas} 
C.~Kassel,  
Asterisque {\bf 177}, 199 (1989), Sem.~Bourbaki.

\bibitem{konvis} 
M.~Kontsevich and S.~Vishik, 
{\em Functional Analysis on the Eve of the 21st Century},
volume~1,  173--197,  (1993).

\bibitem{byt}
I.~Brevik, A.A.~Bytsenko, A.E.~Goncalves, F.L.~Williams,
J. Phys. {\bf A 31}, 4437-4448 (1998),
hep-th/9711159.

\bibitem{evans}
T.S. Evans,
Phys. Lett. {\bf B 457}, 127-132 (1999),
hep-th/9803184.

\bibitem{revans}
E.~Elizalde, A.~Filippi, L.~Vanzo and S.~Zerbini, 
{\em Is the multiplicative anomaly dependent on the regularization~? },  
hep-th/9804071.

\bibitem{dowker}
J.S.~Dowker,
{\em On the relevance of the multiplicative anomaly},
hep-th/9803200.

\bibitem{rdowker}
E.~Elizalde, A.~Filippi, L.~Vanzo and S.~Zerbini, 
{\em Is the multiplicative anomaly relevant~?},
hep-th/980472.

\bibitem{nuovozerbini}
E.~Elizalde, G.~Cognola and S.~Zerbini,
Nucl. Phys {\bf B 532}, 407-428 (1998),
hep-th/9804118.
 
\bibitem{mor}
V.~Moretti,
Commun. Math. Phys. {\bf 210}, 327-363 (1999),
gr-qc/9805091.

\bibitem{mcktom2}
J.J.~McKenzie-Smith and D.J.~Toms,
{\em Zero-point energies and the multiplicative anomaly},
hep-th/0005201.

\bibitem{mcktom}
J.J.~McKenzie-Smith and D.J.~Toms,
Phys. Rev. {\bf D 58}, 105001 (1998),
hep-th/9805184.

\bibitem{raybook}
R.J.~Rivers,
{\em Path Integral Methods in Quantum Field Theory},
Cambridge University Press, Cambridge (1987).

\bibitem{kap1} 
J.I.~Kapusta,  
Phys.~Rev. {\bf D 24}, 426 (1981).

\bibitem{kap2} 
J.I.~Kapusta, 
{\em Finite-Temperature Field Theory},
Cambridge University Press, Cambridge (1989).

\bibitem{habwel1} 
H.E.~Haber and H.A.~Weldon,  
Phys.~Rev.~Lett. {\bf 46}, 1497 (1981).

\bibitem{habwel2}
H.E.~Haber and H.A.~Weldon, 
J.~Math.~Phys. {\bf 23}, 1852 (1982).

\bibitem{habwel3}
H.E.~Haber and H.A.~Weldon, 
Phys.~Rev. {\bf D 25}, 502 (1982).

\bibitem{bd} 
J.~Bernstein and S.~Dodelson, 
Phys.~Rev. Lett. {\bf 66}, 683 (1991).

\bibitem{bbd} 
K.~Benson, J.~Bernstein and S.~Dodelson, 
Phys.~Rev. {\bf D 44}, 2480 (1991).

\bibitem{cogvan} 
G.~Cognola, L.~Vanzo, S.~Zerbini,
Phys.~Rev. {\bf D 47}, 4575 (1993).

\bibitem{kirtom1} 
K.~Kirsten and D.J.~Toms,  
Phys.~Lett. {\bf B 368}, 119 (1996).

\bibitem{chemical} 
A.~Filippi,  
{\em Inclusion of Chemical Potential for Scalar Fields},
hep-ph/9703323.

\end{thebibliography}

\end{document}











