# Multiplicative anomaly and finite charge density

## 1 INTRODUCTION

In field theory we often have to deal with functional determinants of differential operators. These, as formal products of infinite eigenvalues, are divergent objects (UV divergence) and a regularization scheme is therefore necessary. One of the most successful and powerful ones is the zeta-function regularization method . It permits us to give a meaning to the ill defined quantity $`\mathrm{ln}detA`$, where $`A`$ is a second order elliptic differential operator, through the zeta function $`\zeta (s|A)=\text{Tr}A^s`$, 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 to respect to $`s`$ at zero is well defined and the logarithm of the zeta-function regularized functional determinant will then be defined by

$`\mathrm{ln}det{\displaystyle \frac{A}{M^2}}=\zeta ^{}(0|A)\zeta (0|A)\mathrm{ln}M^2,`$ (1)   
where $`M^2`$ is a renormalization scale mass.

Sometimes, however, the differential operator takes a matrix form in the field space, as is the case with the two real components $`\varphi _i`$ of a complex scalar field. In this case we end up evaluating a quantity of the form $`\mathrm{ln}det(AB)`$, with $`A`$ and $`B`$ two commuting pseudo-differential operators. The fact is that it is not always true that the equality $`\mathrm{ln}det(AB)=\mathrm{ln}det(A)+\mathrm{ln}det(B)`$ holds. On the contrary, an additional term $`a(A,B)`$, called the multiplicative anomaly , may be present on the right hand-side and eventually have physical relevance . In this work I will introduce this quantity, compute it and analyse its physical relevance in the case of a charged scalar field at finite temperature and charge density, as well as present other possible physical systems in which it could play a role.

This work has been developed in collaboration with E. Elizalde, in Barcelona (Spain), and L. Vanzo and S. Zerbini, in Trento (Italy). My thank goes also to R. Rivers and T. Evans for stimulating discussions.

## 2 THE BOSE GAS AT FINITE CHARGE DENSITY

The relativistic complex scalar field at finite temperature in the presence of a net charge density has given rise to a certain interest during recent years .

In our recent paper  it is shown that, in a coherent regularized approach, the multiplicative anomaly, overlooked until now, could play a role in this system. I will try here to outline the results avoiding the mathematical machinery. For clarity, I will mainly restrict myself to four space-time dimensions although the system has been studied in generic $`D`$ dimensions.

The relevant quantity for my proposes is the grand canonical partition function, which, for this system, is 

$`Z_\beta (\mu )=\text{Tr}e^{\beta (H\mu Q)}={\displaystyle _{\varphi (\tau )=\varphi (\tau +\beta )}}[d\varphi _i]e^{\frac{1}{2}_0^\beta 𝑑\tau {\scriptscriptstyle d^3x\varphi _iA_{ij}\varphi _j}},`$ (2)   
where $`H`$ is the Hamiltonian of the system, $`Q`$ the charge and $`\beta `$ the inverse of the temperature. $`A_{ij}`$ is the elliptic, non-self-adjoint, matrix valued, differential operator

$`\left(\begin{array}{cc}_\tau ^2^2+m^2e^2\mu ^2& 2ie\mu _\tau \\ 2ie\mu _\tau & _\tau ^2^2+m^2e^2\mu ^2\end{array}\right).`$ (5)   
In this case computing the partition function requires taking both an algebraic determinant and a functional one. The standard procedure consists in taking the algebraic one first . Stimulated by a recent criticism , we showed the validity of this procedure .

Now, we have two possible factorizations for this algebraic determinant:

$`\mathrm{ln}Z_\beta (\mu )={\displaystyle \frac{1}{2}}\mathrm{ln}det{\displaystyle \frac{A_{ij}}{M^2}}={\displaystyle \frac{1}{2}}\mathrm{ln}det\left[{\displaystyle \frac{L_+}{M^2}}{\displaystyle \frac{L_{}}{M^2}}\right]={\displaystyle \frac{1}{2}}\mathrm{ln}det\left[{\displaystyle \frac{K_+}{M^2}}{\displaystyle \frac{K_{}}{M^2}}\right],`$ (6)   
where:

$`K_\pm =^2+m^2+\left(i_\tau \pm ie\mu \right)^2L_\pm =_\tau ^2+\left(\sqrt{^2+m^2}\pm e\mu \right)^2.`$ (7)   
I will avoid here the standard steps that lead to the computation of the logarithm of the partition function, as the reader will find them in greater detail in ref. .

Assuming, as in the precedent literature, the validity of the identity $`\mathrm{ln}det(AB)=\mathrm{ln}det(A)+\mathrm{ln}det(B)`$ and disregarding the anomaly, we obtain, for the $`K_\pm `$ factorization,

$`\mathrm{ln}Z_\beta (K_+,K_{})`$ $`=`$ $`{\displaystyle \frac{\beta V}{32\pi ^2}}\left[m^4(\mathrm{ln}\frac{m^2}{M^2}3/2)\right]`$ (8)   
$``$ $`V{\displaystyle \frac{d^3k}{(2\pi )^3}\mathrm{ln}(1e^{\beta (\sqrt{k^2+m^2}e\mu )})}V{\displaystyle \frac{d^3k}{(2\pi )^3}\mathrm{ln}(1e^{\beta (\sqrt{k^2+m^2}+e\mu )})},`$   
where the expected contributions for vacuum, particles and antiparticles are manifest. From now on I will represent the thermal contributions as $`S(\beta ,\mu )`$.

Similar manipulations can be done for the other factorization $`L_\pm `$. In this case though, the chemical potential does not appear in the sum over Matsubara frequencies, but remains with the momentum integral and therefore the term linear in $`\beta `$ will be chemical potential dependent:

$`\mathrm{ln}Z_\beta (L_+,L_{})`$ $`=`$ $`{\displaystyle \frac{\beta V}{32\pi ^2}}\left[m^4(\mathrm{ln}\frac{m^2}{M^2}3/2)\right]+{\displaystyle \frac{\beta V}{8\pi ^2}}\left({\displaystyle \frac{e^4\mu ^4}{3}}e^2\mu ^2m^2\right)+S(\beta ,\mu ).`$ (9)   
In this system the importance of the multiplicative anomaly is therefore manifest. Despite having $`\mathrm{ln}(K_{}K_+)=\mathrm{ln}(L_{}L_+)`$, these two options give two different results for a zeta-function regularized partition function if the multiplicative anomaly is disregarded.

## 3 THE MULTIPLICATIVE ANOMALY

The multiplicative anomaly  is defined as

$`a_D(A,B)=\mathrm{ln}det(AB)\mathrm{ln}det(A)\mathrm{ln}det(B)`$ (10)   
where the determinants of the two elliptic operators are defined by means of the zeta-function method. I recall that $`D`$ are the space-time dimensions. In principle, it could be computed directly as difference of the involved quantities. In reality, actual calculations are very complicated even for simpler operators. We can fortunately resort to Wodzicki’s results for a remarkably neat recipe.

For any classical pseudo-differential operator $`A`$ there exists a complete symbol $`A(x,k)=e^{ikx}Ae^{ikx}`$. This admits an asymptotic expansion for $`|k|\mathrm{}`$,

$`A(x,k){\displaystyle \underset{j=0}{}}A_{aj}(x,k),`$ (11)   
where the coefficients (their number is infinite) fulfil the homogeneity property $`A_{aj}(x,tk)=t^{aj}A_{aj}(x,k)`$, for $`t>0`$. The number $`a`$ is called the order of $`A`$. Now, Wodzicki  proved that for two invertible, self-adjoint, elliptic, commuting, pseudodifferential operators on a smooth compact manifold without boundaries $`M_D`$:

$`a(A,B)={\displaystyle \frac{\text{res}\left[(\mathrm{ln}(A^bB^a))^2\right]}{2ab(a+b)}}=a(B,A),`$ (12)   
where $`a>0`$ and $`b>0`$ are the orders of $`A`$ and $`B`$, respectively. Here the quantity $`\text{res}(A)`$ is the Wodzicki non-commutative residue. It can be computed easily using the homogeneous component $`A_D(x,k)`$ of order $`D`$ of the complete symbol,

$`\text{res}(A)={\displaystyle _{M_D}}{\displaystyle \frac{dx}{(2\pi )^D}}{\displaystyle _{|k|=1}}A_D(x,k)𝑑k.`$ (13)   
All this can be applied to our operators. As an example:

$`A(x,k)_{K_\pm }=\left[\mathrm{ln}(k^2+m^2e^2\mu ^2+i2e\mu k_\tau )\mathrm{ln}\left(k^2+m^2e^2\mu ^2i2e\mu k_\tau \right)\right]^2.`$ (14)   
Simply expanding (14) and performing the above integration (13) we obtain the non-commutative residue. Remembering (12) and that the order of our operators is 2, we have the related multiplicative anomaly as

$`a_4(K_+,K_{})={\displaystyle \frac{\beta V}{8\pi ^2}}\left[e^2\mu ^2(m^2{\displaystyle \frac{e^2\mu ^2}{3}})\right].`$ (15)   
The same can be done for $`L_\pm `$, obtaining another expression for $`a_4(L_+,L_{})`$.

Finally, including this two results in (8) and (9) respectively, we obtain

$`\mathrm{ln}Z_\beta (K_+,K_{})=\mathrm{ln}Z_\beta (L_+,L_{})`$ $`=`$ $`{\displaystyle \frac{\beta V}{32\pi ^2}}\left[m^4(\mathrm{ln}\frac{m^2}{M^2}3/2)\right]+S(\beta ,\mu )`$ (16)   
$``$ $`{\displaystyle \frac{\beta V}{16\pi ^2}}\left[e^2\mu ^2(m^2{\displaystyle \frac{e^2\mu ^2}{3}})\right],`$   
and the logarithm of the partition function turns to be the same for the two different approaches. Although consistent now, our result is remarkably different from the one in the literature where the multiplicative anomaly was disregarded. The physical relevance of this additional term will be discussed in the next section.

More generally, this term can be easily computed for any space-time dimension $`D`$ and turns out to be always vanishing for odd $`D`$ . It has also been computed for the self-interacting field , but there many difficulties arise when dealing with the regularized determinants of the complicated operators involved .

## 4 PHYSICAL RELEVANCE

To investigate the physical relevance of the multiplicative anomaly the crucial quantity is the effective potential in presence of external sources, which can be expressed as a function of the charge density $`\rho =\frac{1}{\beta V}\frac{\mathrm{ln}Z_\beta (\mu ,J_i)}{\mu }=\frac{<Q>}{V}`$ and the mean field $`x^2=\mathrm{\Phi }^2`$ as 

$`F(\beta ,\rho ,x)={\displaystyle \frac{1}{\beta V}}\mathrm{ln}Z_\beta (\mu )+{\displaystyle \frac{\mu }{\beta V}}{\displaystyle \frac{\mathrm{ln}Z_\beta (\mu )}{\mu }}+{\displaystyle \frac{1}{2}}(m^2+e^2\mu ^2)x^2,`$ (17)   
$`\rho ={\displaystyle \frac{1}{\beta V}}{\displaystyle \frac{\mathrm{ln}Z_\beta (\mu )}{\mu }}+e^2\mu x^2,`$ (18)   
where the later is an implicit expression for the chemical potential as a function of $`\rho `$.

The physical states correspond to the minima of the effective potential, located in $`\frac{F}{x}=x(m^2e^2\mu ^2)=0`$. We find therefore: 1) an unbroken phase, $`x=0`$, $`e\mu <m`$; 2) a symmetry breaking solution, $`x0`$, $`e\mu =\pm m`$, giving the relativistic Bose-Einstein condensation. For our system, explicitly, the unbroken and broken phase are respectively

$`_\beta =minF`$ $`=`$ $`_V{\displaystyle \frac{1}{\beta V}}S(\beta ,\mu )+\mu \rho +{\displaystyle \frac{1}{16\pi ^2}}\left[e^2\mu ^2(m^2{\displaystyle \frac{e^2\mu ^2}{3}})\right],`$ (19)   
$`\rho `$ $`=`$ $`{\displaystyle \frac{1}{\beta V}}{\displaystyle \frac{S(\beta ,\mu )}{\mu }}{\displaystyle \frac{e}{8\pi ^2}}\left[e\mu (m^2{\displaystyle \frac{2e^2\mu ^2}{3}})\right],`$ (20)   
where $`_V`$ is the vacuum contribution, and

$`_\beta `$ $`=`$ $`_V{\displaystyle \frac{1}{\beta V}}S(\beta ,e\mu =m)+{\displaystyle \frac{m}{e}}\rho +{\displaystyle \frac{1}{8\pi ^2}}{\displaystyle \frac{m^4}{3}},`$ (21)   
$`\rho `$ $`=`$ $`{\displaystyle \frac{1}{\beta V}}{\displaystyle \frac{S(\beta ,\mu )}{\mu }}|_{e\mu =m}{\displaystyle \frac{e}{8\pi ^2}}{\displaystyle \frac{m^3}{3}}+emx^2.`$ (22)   
It is possible to see, under a detailed generic $`D`$ analysis, that these expressions, as they are, give some inconsistencies in the broken phase . We have to remember, though, that we worked until now with regularized but “unrenormalized” charge density. Since it appears in the partition function multiplied by $`\mu `$, any ambiguity in it will correspond to an uncertainty in the free energy density of the kind $`\mu K`$. This $`K`$ has to be fixed following physical requirements. A very reasonable one is that the symmetry is unbroken at $`T=0`$, $`\rho =0`$. For $`D=4`$, $`K`$ will be $`K=\frac{em^3}{24\pi ^2}`$. For $`D=4`$ only, this choice also removes the multiplicative anomaly contribution to the charge density, so that the anomaly does not alter the broken phase in any aspect and we get

$`_\beta `$ $`=`$ $`_V{\displaystyle \frac{1}{\beta V}}S(\beta ,m)+{\displaystyle \frac{m}{e}}\rho ^R,`$ (23)   
$`ex^2`$ $`=`$ $`{\displaystyle \frac{1}{m}}\left(\rho ^R+{\displaystyle \frac{1}{\beta V}}{\displaystyle \frac{S(\beta ,\mu )}{\mu }}|_{e\mu =m}\right).`$ (24)   
Also the critical temperature ($`x=0`$, $`e\mu =m`$) remains unchanged with this renormalization. Is different the unbroken phase, where the anomalous term remains:

$`_\beta `$ $`=`$ $`_V{\displaystyle \frac{1}{\beta V}}S(\beta ,\mu )+\mu \rho ^R{\displaystyle \frac{\mu em^3}{12\pi ^2}}+{\displaystyle \frac{1}{8\pi ^2}}\left(e^2\mu ^2m^2{\displaystyle \frac{1}{3}}e^4\mu ^4\right),`$ (25)   
$`\rho ^R`$ $`=`$ $`{\displaystyle \frac{1}{\beta V}}{\displaystyle \frac{S(\beta ,\mu )}{\mu }}{\displaystyle \frac{1}{8\pi ^2}}\left(e^2\mu m^2{\displaystyle \frac{2}{3}}e^4\mu ^3\right)+{\displaystyle \frac{em^3}{24\pi ^2}}.`$ (26)   
We should now observe that the anomalous contribution to the free energy is non leading at ultra relativistic temperatures $`T>m`$, since the thermal terms go as $`T^4`$. On the other hand, it does not even contribute to the low temperature limit, corresponding to the broken phase, so that it could give relevant corrections only in a intermediate range $`Tm`$. Notice, finally, that the anomalous term is vanishing as $`e0`$, and the correct expression of the free energy density for the uncharged boson gas is recovered.

## 5 GENERAL CONSIDERATIONS

This analysed is just one of the many possible physical systems were the multiplicative anomaly could play a role . The first question is therefore if this additional terms will always have physical relevance . This will in general depend on the system. As an example, for the above case the anomaly is vanishing for any odd dimension. In other cases  the anomaly could be simply non physical, as it could be reabsorbed in the renormalization procedure. Work is currently in progress on other systems, including fermionic ones. It is not difficult to see that for a single free fermionic field the anomaly is always vanishing, too. On the other side, it could play a role for neutrino mixing, in relation with recent results regarding inequivalent representations of the vacuum  and, in general, any time when there is a possible mixing or rotation in the functional space of the fields . This needs further investigation due to the deep connection between the multiplicative anomaly and the functional measure, which goes to the roots of the definition of the functional integral itself.

The other relevant question is, of course, if the anomaly is regularization dependent . Zeta function regularization is just one example of a wider class of regularizations called: “generalized proper-time regularizations” , for which we showed the anomaly to be present . This topic is also under further investigations and created a vivid debate lately. Here too, the answers are probably to be found in a proper and consistent definition of the ill-defined functional determinant itself, where this regularization approach has, up to now, proved to be rigorous and coherent .

Note added in proof: After my talk, a work  by McKenzie-Smith and Toms appeared. There, the relevance of considering the multiplicative anomaly within a functional integral approach is recognized although they do not agree on its physical relevance for the relativistic charged bosonic field.