text:
# On the thermal sunset diagram for scalar field theories

## I Introduction

The so-called sunset diagram, which is depicted in Fig.1, is one of two-loop self-energy diagrams. The sunset diagram appears in theories with 4-point vertices such as $`𝒪(4)`$ linear sigma model, $`\varphi ^4`$ theory, and so on.

At finite temperature, two-loop self-energy diagrams have remarkable features that are not seen at zero temperature. Consider the discontinuity of the self-energy. At zero temperature, the discontinuity of one-loop diagrams is due to two particle intermediate states and that of two-loop diagrams comes from three particle states in addition to two particle ones. In general, in higher loop diagrams new processes appear. However, they contribute to the discontinuity only at higher energies. On the other hand, at finite temperature, even though the number of particles participating in the process at a given order of loops is the same as that at zero temperature, new processes appear even at low energy. This is because at finite temperature some of particles participating in the process can be particles in the heat bath. Furthermore, as will be shown in this paper, there exist processes which are possible at arbitrary energy. Accordingly, the discontinuity of two-loop self-energies is non-vanishing in all the energy region. Therefore, in some cases at finite temperature, to extend calculations to two-loop order has a meaning more than just making more precision.

One of such cases is the spectral function for the sigma meson in the linear sigma model. The spectral function of the sigma meson at finite temperature was studied by Chiku and Hatsuda at one-loop level . The sigma meson at zero temperature has a large width due to the strong coupling with two pions. However, they found that at finite temperature the spectral function near $`\sigma \pi \pi `$ threshold is enhanced as a typical signal of chiral phase transition. This is because, as the temperature increases, the mass of $`\sigma `$ decreases while that of $`\pi `$ increases due to partial restoration of chiral symmetry and, accordingly, the phase space available for the $`\sigma \pi \pi `$ decay is squeezed to zero. At finite temperature, however, there exist processes, collision with or absorption of a thermal particle in the heat bath, which contributes to the discontinuity at arbitrary energy for $`\sigma `$, as discussed above. If we include them, the structure of the spectral function might be significantly modified. Their effects on the spectral function cannot be taken into account till we extend the calculation including two-loop self-energies.

In this paper, we take a step forward in the calculation of two-loop self-energy diagrams, i.e. we evaluate the thermal sunset diagram for scalar field theories. After giving a brief review of the real time formalism \[2-14\] in the second section, which we use for the calculation of diagrams, we first examine the bubble diagram, a one-loop diagram which appears as a subdiagram of the sunset diagram, at finite temperature in the third section. Then, we discuss the structure of the discontinuous part of the thermal sunset diagram in the fourth section and explain how we calculate the sunset diagram in the fifth section. Using our result for this diagram, we study the effects of the diagram on the spectral function of $`\sigma `$ at finite temperature in $`𝒪(4)`$ linear sigma model in the sixth section. Finally we summarize the paper in the seventh section.

## II Brief review of the real time formalism

For the calculations of thermal Feynman diagrams we adopt the real time formalism throughout this paper. We briefly review the formalism in this section.

In the real time formalism propagators are given by $`2\times 2`$ matrices. The 4-components of the free propagator of a scalar particle with mass $`m`$ are given by

$`i\mathrm{\Delta }_{11}^F(k;m)`$ $`=`$ $`{\displaystyle \frac{i}{k^2m^2+i\eta }}+n(k_0)2\pi \delta (k^2m^2),`$ (II.1)
$`i\mathrm{\Delta }_{12}^F(k;m)`$ $`=`$ $`e^{\sigma k_0}[n(k_0)+\theta (k_0)]2\pi \delta (k^2m^2),`$ (II.2)
$`i\mathrm{\Delta }_{21}^F(k;m)`$ $`=`$ $`e^{\sigma k_0}[n(k_0)+\theta (k_0)]2\pi \delta (k^2m^2),`$ (II.3)
$`i\mathrm{\Delta }_{22}^F(k;m)`$ $`=`$ $`{\displaystyle \frac{i}{k^2m^2i\eta }}+n(k_0)2\pi \delta (k^2m^2),`$ (II.4)
where $`n(k_0)`$ is the Bose-Einstein distribution function at temperature $`T1/\beta `$:

$$n(k_0)=\frac{1}{\mathrm{exp}(\beta |k_0|)1}.$$
(II.5)
Off-diagonal elements depend on a free parameter $`\sigma `$. In this paper we make the symmetrical choice $`\sigma =\beta /2`$, leading to

$`i\mathrm{\Delta }_{12}^F(k;m)`$ $`=`$ $`i\mathrm{\Delta }_{21}^F(k;m)=e^{\beta |k_0|/2}n(k_0)2\pi \delta (k^2m^2).`$ (II.6)
We denote the matrix of the full propagator by

$`\stackrel{~}{\mathrm{\Delta }}(k;m)=\left(\begin{array}{cc}\mathrm{\Delta }_{11}(k;m)& \mathrm{\Delta }_{12}(k;m)\\ \mathrm{\Delta }_{21}(k;m)& \mathrm{\Delta }_{22}(k;m)\end{array}\right).`$ (II.9)
Each component is the solution of the Schwinger-Dyson equation:

$`\mathrm{\Delta }_{ab}(k;m)=\mathrm{\Delta }_{ab}^F(k;m)+\mathrm{\Delta }_{ac}^F(k;m)\mathrm{\Pi }_{cd}(k)\mathrm{\Delta }_{db}(k;m),`$ (II.10)
where $`\mathrm{\Pi }_{cd}(k)`$ is the self-energy.

The following matrix

$$U(k)=\left(\begin{array}{cc}\sqrt{1+n(k_0)}& \sqrt{n(k_0)}\\ \sqrt{n(k_0)}& \sqrt{1+n(k_0)}\end{array}\right),$$
(II.11)
‘diagonalizes’ the free propagator:

$$U(k)^1\left(\begin{array}{cc}\mathrm{\Delta }_{11}^F(k;m)& \mathrm{\Delta }_{12}^F(k;m)\\ \mathrm{\Delta }_{21}^F(k;m)& \mathrm{\Delta }_{22}^F(k;m)\end{array}\right)U(k)^1=\left(\begin{array}{cc}\mathrm{\Delta }_0^F(k;m)& 0\\ 0& \mathrm{\Delta }_0^F(k;m)^{}\end{array}\right),$$
(II.12)
where $`\mathrm{\Delta }_0^F(k;m)`$ is the Feynman propagator at $`T=0`$

$$\mathrm{\Delta }_0^F(k;m)=\frac{1}{k^2m^2+i\eta }.$$
(II.13)
It is known that $`U(k)`$ also diagonalizes the full propagator  as well as $`\mathrm{\Pi }_{ab}(k)`$

$$\mathrm{\Pi }_{ab}(k)=(U(k)^1)_{ac}\left(\begin{array}{cc}\overline{\mathrm{\Pi }}(k)& 0\\ 0& \overline{\mathrm{\Pi }}(k)^{}\end{array}\right)_{cd}(U(k)^1)_{db}.$$
(II.14)
This matrix equation gives relations between matrix elements:

$`\mathrm{Re}\overline{\mathrm{\Pi }}(k)`$ $`=`$ $`\mathrm{Re}\mathrm{\Pi }_{11}(k),`$ (II.15)
$`\mathrm{Im}\overline{\mathrm{\Pi }}(k)`$ $`=`$ $`\mathrm{tanh}(\beta |k_0|/2)\mathrm{Im}\mathrm{\Pi }_{11}(k),`$ (II.16)
$`\mathrm{Im}\overline{\mathrm{\Pi }}(k)`$ $`=`$ $`i\mathrm{sinh}(\beta |k_0|/2)\mathrm{\Pi }_{12}(k).`$ (II.17)
Let us next find the expression of the spectral function. $`\mathrm{\Delta }_{11}`$ has the following spectral representation:

$$\mathrm{\Delta }_{11}(k;m)=d\omega ^2\rho (\omega ,𝐤)\mathrm{\Delta }_{11}^F(k_0,\omega ),$$
(II.18)
where we denote the argument of $`\mathrm{\Delta }_{11}^F`$ by $`k_0`$ and $`\omega \sqrt{𝐤^2+m^2}`$ instead of $`k`$ and $`m`$. One can prove that, if $`\mathrm{\Delta }_{11}`$ is given in Eq.(II.18), the Green function

$$\mathrm{\Delta }(t,𝐱)=\frac{i}{(2\pi )^4}\mathrm{d}^4ke^{ikx}\mathrm{\Delta }_{11}(k;m)$$
(II.19)
obeys the Kubo-Martin-Schwinger (KMS) condition . From Eq.(II.18) we obtain

$$\rho (k)=\frac{1}{\pi }\mathrm{tanh}(\beta |k_0|/2)\mathrm{Im}\mathrm{\Delta }_{11}(k;m).$$
(II.20)
Our next task is thus to find the expression of $`\mathrm{Im}\mathrm{\Delta }_{11}(k;m)`$. Denoting the diagonalized full propagator by

$`U(k)^1\stackrel{~}{\mathrm{\Delta }}(k;m)U(k)^1`$ $`=`$ $`\left(\begin{array}{cc}\mathrm{\Delta }(k;m)& 0\\ 0& \mathrm{\Delta }(k;m)^{}\end{array}\right),`$ (II.23)
we obtain

$`\stackrel{~}{\mathrm{\Delta }}`$ $`=`$ $`U(k)\left(\begin{array}{cc}\mathrm{\Delta }(k;m)& 0\\ 0& \mathrm{\Delta }(k;m)^{}\end{array}\right)U(k)`$ (II.26)
$`=`$ $`\left(\begin{array}{cc}[1+n(k_0)]\mathrm{\Delta }(k;m)n(k_0)\mathrm{\Delta }(k;m)^{}& \sqrt{n(k_0)[1+n(k_0)]}[\mathrm{\Delta }(k;m)\mathrm{\Delta }(k;m)^{}]\\ \sqrt{n(k_0)[1+n(k_0)]}[\mathrm{\Delta }(k;m)\mathrm{\Delta }(k;m)^{}]& [1+n(k_0)]\mathrm{\Delta }(k;m)^{}+n(k_0)\mathrm{\Delta }(k;m)\end{array}\right).`$ (II.29)
Taking the $`(1,1)`$-components in the both sides of the above equation yields

$$\mathrm{Im}\mathrm{\Delta }_{11}(k;m)=\mathrm{coth}(\beta |k_0|/2)\mathrm{Im}\mathrm{\Delta }(k;m).$$
(II.31)
The expression for $`\mathrm{\Delta }(k;m)`$ can be found in the following manner. In Eq.(II.10), iteratively using Eq.(II.10) itself in the right hand side and diagonalizing, we obtain

$`U(k)^1\stackrel{~}{\mathrm{\Delta }}(k;m)U(k)^1=\left(\begin{array}{cc}\mathrm{\Delta }_0^F(k;m)& 0\\ 0& \mathrm{\Delta }_0^F(k;m)^{}\end{array}\right)`$ (II.34)
$`+\left(\begin{array}{cc}\mathrm{\Delta }_0^F(k;m)& 0\\ 0& \mathrm{\Delta }_0^F(k;m)^{}\end{array}\right)\left(\begin{array}{cc}\overline{\mathrm{\Pi }}(k)& 0\\ 0& \overline{\mathrm{\Pi }}(k)^{}\end{array}\right)\left(\begin{array}{cc}\mathrm{\Delta }_0^F(k;m)& 0\\ 0& \mathrm{\Delta }_0^F(k;m)^{}\end{array}\right)+\mathrm{}.`$ (II.41)
The (1,1) component of this equation gives

$`\mathrm{\Delta }(k;m)`$ $`=`$ $`(U(k)^1\stackrel{~}{\mathrm{\Delta }}(k;m)U(k)^1)_{11}`$ (II.43)
$`=`$ $`\mathrm{\Delta }_0^F(k;m)+\mathrm{\Delta }_0^F(k;m)\overline{\mathrm{\Pi }}(k)\mathrm{\Delta }_0^F(k;m)+\mathrm{}`$ (II.44)
$`=`$ $`{\displaystyle \frac{1}{k^2m^2\overline{\mathrm{\Pi }}(k)}}.`$ (II.45)
Therefore, from Eqs.(II.20), (II.31) and (II.45), we obtain the desired expression of the spectral function:

$$\rho (k)=\frac{1}{\pi }\frac{\mathrm{Im}\overline{\mathrm{\Pi }}(k)}{[k^2m^2\mathrm{Re}\overline{\mathrm{\Pi }}(k)]^2+[\mathrm{Im}\overline{\mathrm{\Pi }}(k)]^2}.$$
(II.46)
Thus the spectral function can be written only with one component of the self-energies while all the components enter $`\mathrm{\Delta }_{11}(k;m)`$ (see Eq.(II.10)).

## III Bubble diagram

In this section we examine the “bubble” diagram, which is a one-loop diagram shown in Fig.2. The bubble diagram appears as a subdiagram of the sunset diagram and therefore we need the expression of the former in the calculation of the latter.

Among the 4-components of the bubble diagram we need only the (1,1) component since the sunset diagram does not have any internal vertices. The (1,1) component of the bubble diagram for scalar particles is given by

$$_{\mathrm{bub}}(p;m_1,m_2)_{11}=\frac{\mathrm{d}^4k}{(2\pi )^4}i\mathrm{\Delta }_{11}^F(p+k;m_1)i\mathrm{\Delta }_{11}^F(k;m_2).$$
(III.47)
Eq.(III.47) can be expressed as the sum of terms with different numbers of Bose-Einstein factors, $`n`$:

$$_{\mathrm{bub}}(p;m_1,m_2)_{11}=I^{(2)}(p^2;m_1,m_2)+(F^{(2)}(p;m_1,m_2)+(12))+F^{(3)}(p;m_1,m_2).$$
(III.48)
Each term in the right hand side is respectively given by

$`I^{(2)}(p^2;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}^4k}{(2\pi )^4}\frac{i}{(p+k)^2m_1^2+i\eta }\frac{i}{k^2m_2^2+i\eta }},`$ (III.49)
$`F^{(2)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}^4k}{(2\pi )^4}\frac{i}{(p+k)^2m_1^2+i\eta }n(k_0)2\pi \delta (k^2m_2^2)},`$ (III.50)
$`F^{(3)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}^4k}{(2\pi )^4}n(p_0+k_0)2\pi \delta ((p+k)^2m_1^2)n(k_0)2\pi \delta (k^2m_2^2)}.`$ (III.51)
Let us carry out the integration in the above equations. The $`T`$-independent part $`I^{(2)}`$ is given in textbooks . In $`d`$-dimension it is given by

$`I^{(2)}(p^2;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{i}{16\pi ^2}}\left({\displaystyle \frac{1}{\overline{ϵ}}}+x_+\mathrm{ln}{\displaystyle \frac{m_1^2}{\kappa ^2}}x_{}\mathrm{ln}{\displaystyle \frac{m_2^2}{\kappa ^2}}2I\right),`$ (III.52)
where $`\kappa `$ is the renormalization point and $`\frac{1}{\overline{ϵ}}`$ and $`x_\pm `$ are

$$\frac{1}{\overline{ϵ}}\frac{1}{ϵ}\gamma +\mathrm{ln4}\pi (ϵ=(4d)/2,\gamma :\mathrm{Euler}\mathrm{constant}),$$
(III.53)
$`x_\pm \pm {\displaystyle \frac{1}{2}}+{\displaystyle \frac{m_1^2m_2^2}{2p^2}}.`$ (III.54)
$`I`$ is given by

$`I`$ $`=`$ $`\{\begin{array}{cc}\sqrt{C}\left[\mathrm{ln}\frac{(x_+\sqrt{C})(x_{}+\sqrt{C})}{(x_{}\sqrt{C})(x_++\sqrt{C})}+i2\pi \right]\hfill & \mathrm{for}p^2>(m_1+m_2)^2,\hfill \\ 2\sqrt{D}\left[\mathrm{arctan}\frac{x_+}{\sqrt{D}}\mathrm{arctan}\frac{x_{}}{\sqrt{D}}\right]\hfill & \mathrm{for}(m_1m_2)^2<p^2<(m_1+m_2)^2,\hfill \\ \sqrt{C}\mathrm{ln}\frac{(\sqrt{C}x_+)(\sqrt{C}+x_{})}{(\sqrt{C}x_{})(\sqrt{C}+x_+)}\hfill & \mathrm{for}p^2<(m_1m_2)^2,\hfill \end{array}`$ (III.58)
where $`C`$ and $`D`$ are

$`C`$ $``$ $`{\displaystyle \frac{1}{4}}\left[1{\displaystyle \frac{(m_1+m_2)^2}{p^2}}\right]\left[1{\displaystyle \frac{(m_1m_2)^2}{p^2}}\right],DC.`$ (III.60)
We see from Eqs.(III.49) and (LABEL:I) that the discontinuous part, $`\mathrm{Im}iI^{(2)}`$, is non-vanishing for $`p^2>(m_1+m_2)^2`$.

The divergent part in Eq.(III.49) is renormalized by applying a counter term to the Lagrangian. The finite part is determined in such a way that the resultant self-energy satisfies a proper normalization condition.

Let us next turn to the $`T`$-dependent part, $`F^{(2)}`$ and $`F^{(3)}`$. We first discuss their discontinuous parts. As shown in Appendix A the discontinuous part of $`F^{(2)}(p;m_1,m_2)`$ is given analytically as follows:

1. For $`p^2>(m_1+m_2)^2`$ or $`0<p^2<(m_1m_2)^2`$,

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{16\pi |𝐩|}}{\displaystyle \frac{1}{\beta }}\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta \omega _+}}{1\mathrm{e}^{\beta \omega _{}}}}\right|.`$ (III.61)
2. For $`(m_1m_2)^2<p^2<(m_1+m_2)^2`$,

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`0.`$ (III.62)
3. For $`p^2<0`$,

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{16\pi |𝐩|}}{\displaystyle \frac{1}{\beta }}\mathrm{ln}\left|(1\mathrm{e}^{\beta \omega _+})(1\mathrm{e}^{\beta \omega _{}})\right|.`$ (III.63)

Here,

$$\omega _\pm =\frac{1}{2}\left|\sqrt{\left(1+\frac{m_2^2m_1^2}{p^2}\right)^2}|p_0|\pm \sqrt{\left[1\frac{(m_2+m_1)^2}{p^2}\right]\left[1\frac{(m_1m_2)^2}{p^2}\right]}|𝐩|\right|.$$
(III.64)
$`F^{(3)}`$ has only discontinuous part. Its calculation can be done in a way similar to that for $`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$. We show the final results:

1. For $`p^2>(m_1+m_2)^2`$,

$`F^{(3)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|\beta }}{\displaystyle \frac{1}{\mathrm{e}^{\beta |p_0|}1}}\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta \omega _+}}{1\mathrm{e}^{\beta \omega _{}}}}{\displaystyle \frac{\mathrm{e}^{\beta (|p_0|\omega _{})}1}{\mathrm{e}^{\beta (|p_0|\omega _+)}1}}\right|.`$ (III.65)
2. For $`(m_1m_2)^2<p^2<(m_1+m_2)^2`$,

$`F^{(3)}(p;m_1,m_2)`$ $`=`$ $`0.`$ (III.66)
3. For $`0<p^2<(m_1m_2)^2`$,

$`F^{(3)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|\beta }}{\displaystyle \frac{1}{\mathrm{e}^{\beta |p_0|}1}}`$ (III.68)
$`\times \left[\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta \omega _+}}{1\mathrm{e}^{\beta \omega _{}}}}\right|\mathrm{e}^{\beta |p_0|}\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta (|p_0|+\omega _+)}}{1\mathrm{e}^{\beta (|p_0|+\omega _{})}}}\right|\right].`$
4. For $`p^2<0`$,

$`F^{(3)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|\beta }}\{{\displaystyle \frac{1}{e^{\beta |p_0|}1}}[\mathrm{ln}|1\mathrm{e}^{\beta \omega _+}|+\mathrm{e}^{\beta |p_0|}\mathrm{ln}|1\mathrm{e}^{\beta (\omega _+|p_0|)}|]`$ (III.70)
$`+{\displaystyle \frac{1}{e^{\beta |p_0|}1}}[\mathrm{ln}|1\mathrm{e}^{\beta \omega _{}}|+\mathrm{e}^{\beta |p_0|}\mathrm{ln}|1\mathrm{e}^{\beta (\omega _{}+|p_0|)}|]\}.`$

In the above equations $`\omega _\pm `$ are given by Eq.(III.64) with $`m_1`$ and $`m_2`$ replaced by $`\mathrm{max}\{m_1,m_2\}`$ and $`\mathrm{min}\{m_1,m_2\}`$, respectively, since $`F^{(3)}(p;m_1,m_2)`$ is symmetric with respect to $`m_1`$ and $`m_2`$ by definition.

It should be noted that the $`T`$-dependent part has two cuts in the complex $`p^2`$ plane: one starts from $`p^2=(m_1+m_2)^2`$ to the right along the real axis and the other from $`p^2=(m_1m_2)^2`$ to the left.

Let us next discuss the continuous part of $`F^{(2)}`$, which is

$`\mathrm{Re}iF^{(2)}(p;m_1,m_2)={\displaystyle \frac{\mathrm{d}^4k}{(2\pi )^4}𝒫\frac{1}{(p+k)^2m_1^2}n(k_0)2\pi \delta (k^2m_2^2)},`$ (III.72)
where $`𝒫`$ stands for the prescription of Cauchy’s principal value. After integration over $`k_0`$ and angle, we can express Eq.(III.72) as a one-dimensional integral:

$`\mathrm{Re}iF^{(2)}(p;m_1,m_2)`$ (III.73)
$`=`$ $`{\displaystyle \frac{1}{16\pi ^2|𝐩|}}{\displaystyle _{m_2}^{\mathrm{}}}d\omega n(\omega )`$ (III.74)
$`\times `$ $`𝒫\mathrm{ln}\left|{\displaystyle \frac{(p^2+2p_0\omega 2|𝐩|\sqrt{\omega ^2m_2^2}+m_2^2m_1^2)(p^22p_0\omega 2|𝐩|\sqrt{\omega ^2m_2^2}+m_2^2m_1^2)}{(p^2+2p_0\omega +2|𝐩|\sqrt{\omega ^2m_2^2}+m_2^2m_1^2)(p^22p_0\omega +2|𝐩|\sqrt{\omega ^2m_2^2}+m_2^2m_1^2)}}\right|,`$ (III.75)
whose integration will be carried out numerically.

## IV Structure of the discontinuous part of the sunset diagram

Before proceeding to the calculation of the thermal sunset diagram, we study the structure of the discontinuous part of the diagram. In ref., Weldon analyzed the discontinuous part of the bubble diagram. We will generalize it for the sunset diagram. For this purpose it is convenient to use Eq.(II.17) , namely

$`\mathrm{Im}i\overline{}_{\mathrm{sun}}(k;m_1,m_2,m_3)=i\mathrm{sinh}(\beta |k_0|/2)i_{\mathrm{sun}}(k;m_1,m_2,m_3)_{12}.`$ (IV.77)
Here $`\overline{}_{\mathrm{sun}}(k;m_1,m_2,m_3)`$ is the $`(1,1)`$ component of the diagonalized self-energy matrix for the sunset diagram. We note that the self-energy which actually enters spectral functions is diagonalized one (see Eq.(II.46)). $`_{\mathrm{sun}}(k;m_1,m_2,m_3)_{12}`$ is the $`(1,2)`$ component of the sunset diagram, which is given by the following integral:

$`i_{\mathrm{sun}}(k;m_1,m_2,m_3)_{12}`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}^4p}{(2\pi )^4}i\mathrm{\Delta }_{12}^F(p;m_1)i_{\mathrm{bub}}(kp;m_2,m_3)_{12}}.`$ (IV.78)
For $`i\mathrm{\Delta }_{12}^F(p;m_1)`$ it is convenient to use another form:

$`i\mathrm{\Delta }_{12}^F(p;m_1)`$ $`=`$ $`e^{\beta p_0/2}f(p_0)ϵ(p_0)2\pi \delta (p^2m_1^2),`$ (IV.79)
$`f(p_0)`$ $`=`$ $`{\displaystyle \frac{1}{e^{\beta p_0}1}}.`$ (IV.80)
$`i_{\mathrm{bub}}(p;m_1,m_2)_{12}`$ denotes the $`(1,2)`$ component of the bubble diagram:

$`i_{\mathrm{bub}}(kp;m_2,m_3)_{12}`$ $`=`$ $`i{\displaystyle \frac{\mathrm{d}^4q}{(2\pi )^4}i\mathrm{\Delta }_{12}^F(q;m_2)i\mathrm{\Delta }_{12}^F(kpq;m_3)}.`$ (IV.81)
One can reduce this equation to

$`i_{\mathrm{bub}}(p;m_1,m_2)_{12}`$ (IV.82)
$`=`$ $`ie^{\beta (k_0p_0)/2}f(k_0p_0)`$ (IV.85)
$`\times {\displaystyle }{\displaystyle \frac{\mathrm{d}^4q}{(2\pi )^2}}{\displaystyle \frac{1}{4E_2E_3}}\{(1+n_2+n_3)[\delta (k_0p_0E_2E_3)\delta (k_0p_0E_2+E_3)]`$
$`(n_2n_3)[\delta (k_0p_0E_2+E_3)\delta (k_0p_0+E_2E_3)]\},`$
where $`E_2=\sqrt{|𝐪|^2+m_2^2}`$ and $`E_3=\sqrt{|𝐤𝐩𝐪|^2+m_3^2}`$. We also define $`E_1=\sqrt{|𝐩|^2+m_1^2}`$ for later use. $`n_i`$ is the Bose-Einstein factor defined by $`n_i=n(E_i)`$. Using Eq.(IV.77) and Eq.(IV.78) in which Eq.(IV.85) is substituted yields

$`\mathrm{Im}i\overline{}_{\mathrm{sun}}(k;m_1,m_2,m_3)`$ (IV.86)
$`=`$ $`\pi ϵ(k_0){\displaystyle \frac{\mathrm{d}^3p}{(2\pi )^3}\frac{\mathrm{d}^3q}{(2\pi )^3}\frac{1}{8E_1E_2E_3}}`$ (IV.91)
$`\times \{((1+n_1)(1+n_2)(1+n_3)n_1n_2n_3)\delta (k_0E_1E_2E_3)`$
$`+(n_1n_2n_3(1+n_1)(1+n_2)(1+n_3))\delta (k_0+E_1+E_2+E_3)`$
$`+[(n_2n_3(1+n_1)n_1(1+n_2)(1+n_3))\delta (k_0E_1+E_2+E_3)`$
$`+(n_1(1+n_2)(1+n_3)n_2n_3(1+n_1))\delta (k_0+E_1E_2E_3)+(\mathrm{permutations})]\}.`$
In this equation “permutations” stands for the terms obtained by permuting the particle labels of the third and the forth terms.

Let us now consider the physical content of Eq.(LABEL:imsunset2). The first term in Eq.(LABEL:imsunset2) may be interpreted as the probability for the decay $`\mathrm{\Phi }\varphi _1\varphi _2\varphi _3`$ with the statistical weight $`(1+n_1)(1+n_2)(1+n_3)`$ for stimulated emission minus the probability for the creation $`\varphi _1\varphi _2\varphi _3\mathrm{\Phi }`$ with the weight $`n_1n_2n_3`$ for absorption. The second term is the anti-particle counter part of the first term. The third term represents the probability for $`\mathrm{\Phi }\overline{\varphi }_2\overline{\varphi }_3\varphi _1`$ with the weight $`n_2n_3(1+n_1)`$ minus that for $`\varphi _1\mathrm{\Phi }\overline{\varphi }_2\overline{\varphi }_3`$ with the weight $`n_1(1+n_2)(1+n_3)`$. Here $`\overline{\varphi }_i`$ stands for the anti-particle of $`\varphi _i`$. The forth term is the anti-particle counter part of the third term. It represents the probability for $`\mathrm{\Phi }\overline{\varphi }_1\varphi _2\varphi _3`$ with the weight $`n_1(1+n_2)(1+n_3)`$ minus that for $`\varphi _2\varphi _3\mathrm{\Phi }\overline{\varphi }_1`$ with the weight $`n_2n_3(1+n_1)`$. All processes are shown in Fig.3.

We next find the region of $`k^2`$ where the physical processes contained in Eq.(LABEL:imsunset2) are possible, which is equivalent to looking for the condition under which the integral over q in Eq.(LABEL:imsunset2) survives. For the first and the second terms in Eq.(LABEL:imsunset2) to be non-vanishing, $`k^2`$ must satisfies the condition $`k^2>(m_1+M_{23})^2`$, where $`M_{23}`$ is the invariant mass of $`\varphi _2`$ and $`\varphi _3`$. Therefore, the processes in Fig.3 (a), (b), (c) and (d) are possible for $`k^2>(m_1+m_2+m_3)^2`$ since $`M_{23}>m_2+m_3`$. The third and the forth terms survive when $`k^2<(m_1M_{23})^2`$. The processes in Fig.3 (e), (f), (g) and (h), therefore, take place at arbitrary $`k^2`$. This is reasonable since the processes in Fig.3 (g) and (h) are scattering ones and since those in Fig.3 (e) and (f) can be also regarded as scattering by interpreting the incoming $`\mathrm{\Phi }`$ in (e) as an outgoing anti-$`\mathrm{\Phi }`$ and by doing the outgoing $`\mathrm{\Phi }`$ in (f) as an incoming anti-$`\mathrm{\Phi }`$. Accordingly, the discontinuous part of the sunset diagram is non-vanishing for arbitrary $`k^2`$, which is a remarkable feature of the thermal self-energy at and beyond two-loop order.

## V Calculation of the thermal sunset diagram

In this section, we explain how to calculate the thermal sunset diagram by reducing it to an expression written in terms of the bubble diagram previously obtained.

The (1,1) component of the sunset diagram shown in Fig.1 is given by

$$_{\mathrm{sun}}(k;m_1,m_2,m_3)_{11}=\frac{\mathrm{d}^4p}{(2\pi )^4}i\mathrm{\Delta }_{11}^F(p;m_1)\frac{\mathrm{d}^4q}{(2\pi )^4}i\mathrm{\Delta }_{11}^F(q;m_2)i\mathrm{\Delta }_{11}^F(kpq;m_3).$$
(V.93)
We decompose Eq.(V.93) into terms without and with Bose-Einstein factors, which we denote by $`_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$ and $`_{\mathrm{sun}}^{\mathrm{th}}(k;m_1,m_2,m_3)_{11}`$, respectively:

$$_{\mathrm{sun}}(k;m_1,m_2,m_3)_{11}=_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}+_{\mathrm{sun}}^{\mathrm{th}}(k;m_1,m_2,m_3)_{11}.$$
(V.94)
They are given by

$`_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}^4p}{(2\pi )^4}\frac{i}{p^2m_1^2+i\eta }I^{(2)}((kp)^2;m_2,m_3)},`$ (V.95)
$`_{\mathrm{sun}}^{\mathrm{th}}(k;m_1,m_2,m_3)_{11}`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}^4p}{(2\pi )^4}\frac{i}{p^2m_1^2+i\eta }}`$ (V.97)
$`\times [(F^{(2)}(kp;m_2,m_3)+(23))+F^{(3)}(kp;m_2,m_3)]`$
$`+`$ $`{\displaystyle }{\displaystyle \frac{\mathrm{d}^4p}{(2\pi )^4}}n(p_0)2\pi \delta (p^2m_1^2)[I^{(2)}((kp)^2;m_2,m_3)`$ (V.99)
$`+(F^{(2)}(kp;m_2,m_3)+(23))+F^{(3)}(kp;m_2,m_3)].`$
Here we have expressed the second integral in Eq.(V.93) in terms of $`I^{(2)}`$, $`F^{(2)}`$ and $`F^{(3)}`$. The $`T`$-independent part, Eq.(V.95), has a subdivergence coming from the nested bubble diagram and a two-loop overall divergence. On the other hand, $`T`$-dependent part, Eq.(LABEL:sunset1), has only a subdivergence, which can be removed by carrying out renormalization at one-loop level. Hereafter, $`I^{(2)}`$ expresses that with the divergence removed. We calculate Eqs.(V.95) and (LABEL:sunset1) separately.

### V.1 $`T`$-independent part

The $`T`$-independent part of sunset type diagrams has been calculated by several authors \[19-22\] so far. In this paper we calculate its finite part using a dispersion relation.

We first find the discontinuous part of $`_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$ and then compute the finite continuous part using the obtained discontinuous part, via the twice-subtracted dispersion relation.

In order to calculate the discontinuous part of $`_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$, we substitute a dispersion relation for $`iI^{(2)}((kp)^2;m_2,m_3)`$:

$`iI^{(2)}((kp)^2;m_2,m_3)`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}{\displaystyle _{(m_2+m_3)^2}^{\mathrm{}}}dM^2{\displaystyle \frac{\mathrm{Im}iI^{(2)}(M^2;m_2,m_3)}{M^2(kp)^2i\eta }},`$ (V.101)
into Eq.(V.95). Then, we obtain

$`_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}{\displaystyle _{(m_2+m_3)^2}^{\mathrm{}}}dM^2\mathrm{Im}iI^{(2)}(M^2;m_2,m_3)`$ (V.103)
$`\times {\displaystyle }{\displaystyle \frac{\mathrm{d}^4p}{(2\pi )^4}}{\displaystyle \frac{i}{p^2m_1^2+i\eta }}{\displaystyle \frac{i}{(kp)^2M^2+i\eta }}`$
$`=`$ $`{\displaystyle \frac{1}{\pi }}{\displaystyle _{(m_2+m_3)^2}^{\mathrm{}}}dM^2\mathrm{Im}iI^{(2)}(M^2;m_2,m_3)I^{(2)}(k^2;m_1,M).`$ (V.104)
We take the discontinuous part of this equation:

$`\mathrm{Im}i_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}{\displaystyle _{(m_2+m_3)^2}^{\mathrm{}}}dM^2\mathrm{Im}iI^{(2)}(M^2;m_2,m_3)\mathrm{Im}iI^{(2)}(k^2;m_1,M).`$ (V.105)
Using the expression for $`\mathrm{Im}iI^{(2)}(p^2;m_,m_2)`$ we can easily evaluate Eq.(LABEL:sunsetvac-im) numerically.

Next we turn to the continuous part. The continuous part of Eq.(V.95) is divergent. Therefore one needs corresponding counter terms in the Lagrangian. The second and third diagrams of Fig.4 appear at two-loop order, which cancel the divergences in the bare sunset diagram, the first diagram of Fig.4. In this paper, we do not explicitly go through the renormalization procedure but concentrate on the finite part.

In order to calculate the finite part, we use the twice-subtracted dispersion relation:

$`\mathrm{Re}i\stackrel{~}{}_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}`$ (V.107)
$``$ $`\mathrm{Re}\left\{i_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}i_{\mathrm{sun}}^{\mathrm{vac}}(0;m_1,m_2,m_3)_{11}k^2\left[{\displaystyle \frac{}{k^2}}i_{\mathrm{sun}}^{\mathrm{vac}}(k^2;m_1,m_2,m_3)_{11}\right]_{k^2=0}\right\}`$ (V.108)
$`=`$ $`{\displaystyle \frac{k^4}{\pi }}{\displaystyle _{(m_1+m_2+m_3)^2}^{\mathrm{}}}dM^2𝒫{\displaystyle \frac{\mathrm{Im}i_{sunset}^{vac}(M^2;m_1,m_2,m_3)}{M^4(M^2k^2)}}.`$ (V.109)
The second and the third terms in the left hand side are divergent and have to be renormalized but the right hand side is finite and remain unchanged by renormalization. Using the results for $`\mathrm{Im}i_{\mathrm{sun}}^{\mathrm{vac}}(M^2;m_1,m_2,m_3)`$ we can easily evaluate Eq.(V.109) numerically. After $`\mathrm{Re}i\stackrel{~}{}_{\mathrm{sun}}^{\mathrm{vac}}(m_{phys}^2;m_1,m_2,m_3)_{11}`$ is subtracted, Eq.(V.109) coincides with that in modified minimal subtraction scheme up to an irrelevant overall factor.

### V.2 $`T`$-dependent part

We rewrite the integrand of the first term in the $`T`$-dependent part Eq.(LABEL:sunset1) so that it becomes the same form as that in the second term, i.e. (a delta function)$`\times `$(some functions). This can be done by recombining the two factors in the integrand of $`F^{(2)}(kp;m_2,m_3)+(23)`$ and $`F^{(3)}(kp;m_2,m_3)`$, $`T`$-independent part and $`T`$-dependent part of $`i\mathrm{\Delta }_{11}^F`$. The result is as follows:

$`_{\mathrm{sun}}^{\mathrm{th}}(k;m_1,m_2,m_3)_{11}`$ (V.110)
$`=`$ $`{\displaystyle }{\displaystyle \frac{\mathrm{d}^4p}{(2\pi )^4}}\{n(p_0)2\pi \delta (p^2m_1^2)[I^{(2)}(kp;m_2,m_3)+(F^{(2)}(kp;m_2,m_3)+(23))`$ (V.114)
$`+F^{(3)}(kp;m_2,m_3)]`$
$`+n(p_0)2\pi \delta (p^2m_2^2)\left[I^{(2)}(kp;m_1,m_3)+F^{(2)}(kp;m_1,m_3)\right]`$
$`+n(p_0)2\pi \delta (p^2m_3^2)I^{(2)}(kp;m_1,m_2)\}.`$
When we put $`k=(k_0,\mathrm{𝟎})`$, Eq.(V.114) is reduced to

$`_{\mathrm{sun}}^{\mathrm{th}}(k_0,\mathrm{𝟎};m_1,m_2,m_3)_{11}`$ (V.115)
$`=`$ $`{\displaystyle \frac{1}{4\pi ^2}}{\displaystyle \underset{\tau =\pm }{}}\{{\displaystyle _{m_1}^{\mathrm{}}}\mathrm{d}\omega n(\omega )\sqrt{\omega ^2m_1^2}[I^{(2)}(k_0+\tau \omega ,\sqrt{\omega ^2m_1^2};m_2,m_3)`$ (V.121)
$`+(F^{(2)}(k_0+\tau \omega ,\sqrt{\omega ^2m_1^2};m_2,m_3)+(23))`$
$`+F^{(3)}(k_0+\tau \omega ,\sqrt{\omega ^2m_1^2};m_2,m_3)]`$
$`+{\displaystyle _{m_2}^{\mathrm{}}}\mathrm{d}\omega n(\omega )\sqrt{\omega ^2m_2^2}[I^{(2)}(k_0+\tau \omega ,\sqrt{\omega ^2m_2^2};m_1,m_3)`$
$`+F^{(2)}(k_0+\tau \omega ,\sqrt{\omega ^2m_2^2};m_1,m_3)]`$
$`+{\displaystyle _{m_3}^{\mathrm{}}}\mathrm{d}\omega n(\omega )\sqrt{\omega ^2m_3^2}I^{(2)}(k_0+\tau \omega ,\sqrt{\omega ^2m_3^2};m_1,m_2)\},`$
Using the expressions of $`I^{(2)}`$, $`F^{(2)}`$ and $`F^{(3)}`$ obtained in the previous section, we can evaluate the continuous and discontinuous parts of Eq.(V.121) numerically.

## VI Contribution of the thermal sunset diagram to the spectral function of the sigma meson at finite temperature

The purpose of this section is to see how two-loop diagrams affect observables by evaluating the contribution of the thermal sunset diagram to $`\sigma `$ spectral function at finite temperature in the $`𝒪(4)`$ linear sigma model.

It is known that naive perturbation theory breaks down at $`T0`$ and that resummation of higher orders is necessary . We adopt here a resummation technique called optimized perturbation theory (OPT) . We first briefly review the procedure of OPT applied to $`𝒪(4)`$ linear sigma model. The original linear sigma model Lagrangian is as follows:

$$=\frac{1}{2}\left[(_\mu \varphi _i)^2\mu ^2\varphi _i^2\right]\frac{\lambda }{4!}(\varphi _i^2)^2+h\varphi _0+\mathrm{counter}\mathrm{terms},$$
(VI.122)
where $`\varphi _i=(\sigma ,\stackrel{}{\pi })`$ and $`h\varphi _0`$ being the explicitly symmetry breaking term. For the renormalized couplings $`\mu ^2`$, $`\lambda `$ and $`h`$ and the renormalization point $`\kappa `$ we use the values determined in : $`\mu ^2=(283\mathrm{MeV})^2`$, $`\lambda =73.0`$, $`h=(123\mathrm{MeV})^3`$, $`\kappa =255\mathrm{MeV}`$.

In OPT one adds and subtracts a new mass term with the mass $`m`$ to the Lagrangian. Thus, we have

$``$ $`=`$ $`{\displaystyle \frac{1}{2}}[(_\mu \varphi _i)^2m^2\varphi _i^2]+{\displaystyle \frac{1}{2}}\chi \varphi _i^2{\displaystyle \frac{\lambda }{4!}}(\varphi _i^2)^2+h\varphi _0+\text{(counter term)},`$ (VI.123)
where $`\chi m^2\mu ^2`$. The idea of OPT is reorganization of perturbation theory: one treats the added one as a tree-level mass term while the subtracted one as perturbation.

When the spontaneous symmetry breaking takes place, tree level masses of $`\pi `$ and $`\sigma `$ read, respectively,

$$m_{0\pi }^2=m^2+\frac{\lambda }{6}\xi ^2\text{}m_{0\sigma }^2=m^2+\frac{\lambda }{2}\xi ^2,$$
(VI.124)
where $`\xi `$ is the vacuum expectation value of $`\sigma `$ and determined by the stationary condition for the thermal effective action $`V(\xi ,T,m^2)`$ :

$$\frac{V(\xi ,T,m^2)}{\xi }=0.$$
(VI.125)
Note that the derivative with respect to $`\xi `$ does not act on $`m^2`$.

If Green’s functions are calculated in all orders in OPT, they should not depend on the arbitrary mass, $`m`$. However, if one truncates perturbation series at a certain order they depend on it. One can determine this arbitrary parameter so that the correction terms are as small as possible. We adopt the following condition :

$$\mathrm{\Pi }_\pi (k^2=m_{0\pi }^2)+\mathrm{\Pi }_\pi (k=0;T)=0,$$
(VI.126)
where the first and second terms are respectively $`T`$-independent part and $`T`$-dependent part of the one-loop self-energy of $`\pi `$.

Let us now turn to the discussion on the spectral function of $`\sigma `$ defined by Eq.(II.46):

$$\rho _\sigma (k_0,𝐤)=\frac{1}{\pi }\frac{\mathrm{Im}\overline{\mathrm{\Pi }}_\sigma (k_0,𝐤)}{(k^2m_{0\sigma }^2\mathrm{Re}\overline{\mathrm{\Pi }}_\sigma (k_0,𝐤))^2+(\mathrm{Im}\overline{\mathrm{\Pi }}_\sigma (k_0,𝐤))^2}.$$
(VI.127)
Here $`\overline{\mathrm{\Pi }}_\sigma (k_0,𝐤)`$ is the self-energy of $`\sigma `$. Its real and imaginary parts are related with (1,1) component, $`\mathrm{\Pi }_\sigma ^{11}`$, via Eqs.(II.15) and (II.16).

As was already mentioned, the spectral function at one-loop order was studied by Chiku and Hatsuda . We want to see how the spectral function at one-loop order is modified by adding the thermal sunset diagram. Thus, we take

$$\mathrm{\Pi }_\sigma ^{11}(k_0,𝐤)=\mathrm{\Pi }_\sigma ^{11}(k_0,𝐤)_{1\mathrm{loop}}+\mathrm{\Pi }_\sigma ^{11}(k_0,𝐤)_{\mathrm{sun}}.$$
(VI.128)
The first term is the renormalized one-loop self-energy calculated in . The second term is the renormalized sunset diagram depicted in Fig.5 and given by

$`\mathrm{\Pi }_\sigma ^{11}(k_0,𝐤)_{\mathrm{sun}}`$ $`=`$ $`{\displaystyle \frac{\lambda ^2}{6}}i_{\mathrm{sun}}(k_0,𝐤;m_{0\pi },m_{0\pi },m_{0\sigma })_{11}.`$ (VI.129)
We show $`\rho _\sigma (k_0,𝐤=\mathrm{𝟎})`$ at $`T=200\mathrm{MeV}`$ and $`T=145\mathrm{MeV}`$ in Fig.6.

At $`T=200\mathrm{MeV}`$, the spectral function at one-loop order consists of a $`\delta `$-function peak for $`\sigma `$ and a continuum. By including the sunset diagram $`\sigma `$ acquires a width of the order of $`10\mathrm{MeV}`$. At lower temperature, $`T=145\mathrm{MeV}`$, an enhancement of the spectrum near the threshold is observed at one-loop order. When we include the sunset diagram, this feature is not lost.

Let us discuss the above results. At high temperature ($`T=200\mathrm{MeV}`$), the mass of $`\sigma `$ is smaller than twice the pion mass and the decay $`\sigma \pi \pi `$ is forbidden. As a result, within the one-loop calculation $`\sigma `$ has zero width. However, at finite temperature $`\sigma `$ can interact with thermal particles in heat bath and change into other states. Among such processes, those which are taken into account by including the sunset diagram are represented by Fig.3 with $`(\mathrm{\Phi },\varphi _1,\varphi _2,\varphi _3)`$ assigned to, for example, $`(\sigma ,\sigma ,\pi ,\pi )`$. The processes which correspond to (a) and (b) in Fig.3 are possible for $`k_0>2m_{0\pi }+m_{0\sigma }`$. This affect the spectrum at high energy. (c) and (d) in Fig.3 drop off for positive $`k_0`$. (e) and (f) affect the spectrum at low energy since they are possible for $`k_0<m_{0\sigma }2m_{0\pi }`$. The processes which correspond to (g) and (h) in Fig.3 are shown in Fig.7. They are the most important and give a finite width to $`\sigma `$ since they are allowed at arbitrary positive $`k_0`$. However, we observe that their effects at lower temperature ($`T=145\mathrm{MeV}`$) are small. The reason is traced back to Eq.(LABEL:imsunset2). The term representing the probabilities for (g) and (h) in Fig.3 is the fourth term. That integral at lower temperature is suppressed due to the statistical weight since at lower temperature the masses of $`\sigma `$ and $`\pi `$ are large.

Finally we note that, as a consequence of the non-vanishing discontinuous part of the sunset diagram, the spectral function is also non-vanishing in the all range of $`k_0`$.

## VII Summary and perspective

We have studied the sunset diagram for scalar field theories at finite temperature in the real time formalism. We have explained how we can reduce it to an expression written in terms of one-loop self-energy integrals, which can be easily evaluated numerically. We have also discussed what physical processes are contained in the discontinuous part of the diagram. We have found that there exist processes which occur only at finite temperature and some of them are allowed at arbitrary energy. As a result, the discontinuous part of the sunset diagram is non-vanishing in all the energy region, which is a remarkable feature at finite temperature and manifests itself at two-loop order.

As an application of the result, we have demonstrated how the spectral function of $`\sigma `$ at finite $`T`$ at one-loop order is modified when we include the thermal sunset diagram. At high temperature, where $`\sigma \pi \pi `$ is forbidden, $`\sigma `$ acquires a finite width of the order of $`10\mathrm{MeV}`$ due to collisions with thermal particles in the heat bath while $`\sigma `$ does not have a width at one-loop order. At lower temperature the spectrum is almost unchanged.

Finally we comment on the effect of other two-loop diagrams on the spectral function. We have seen that the threshold enhancement, which was first found in the one-loop calculation, is retained if we include the sunset diagram. However, in the present calculation the effect of the thermal width of $`\pi `$ in the $`\sigma \pi \pi `$ decay is not included. This is taken into account by including the diagrams such as shown in Fig.8, in which an internal $`\pi `$ changes into $`\sigma `$ absorbing a thermal $`\pi `$ in the heat bath.

In the one-loop calculation, $`\pi `$ has a width, $`\mathrm{\Gamma }_\pi =50\mathrm{MeV}`$, at the temperature, $`T=145\mathrm{MeV}`$, at which the threshold enhancement is observed for $`\sigma `$. If we include this effect as a constant complex mass shift for $`\pi `$ in the one-loop self-energy for $`\sigma `$, we expect $`\sigma `$ to acquire a width twice as that for $`\pi `$, i.e. $`\mathrm{\Gamma }_\sigma =2\mathrm{\Gamma }_\pi `$ . This implies that when the two-loop diagrams such as Fig.8 are included the spectral function for $`\sigma `$ would be significantly modified with the width of about $`100\mathrm{MeV}`$. The calculations of those diagrams are now in progress .

###### Acknowledgements.

The authors would like to thank M. Ohtani for useful discussion. This work was partially supported by Grants-in-Aid of the Japanese Ministry of Education, Science, Sports, Culture and Technology (No.06572).

## Appendix A Calculation of $`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ and $`F^{(3)}(p;m_1,m_2)`$

In this appendix, we derive Eqs.(III.61)-(III.63) and Eqs.(III.65)-(III.68).

From Eq.(III.50) we obtain

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{\mathrm{d}^4k}{(2\pi )^4}2\pi \delta ((p+k)^2m_1^2)n(k_0)2\pi \delta (k^2m_2^2)}`$ (A.130)
$`=`$ $`{\displaystyle \frac{1}{4\pi }}{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑k_0{\displaystyle _0^{\mathrm{}}}d|𝐤||𝐤|^2n(k_0)\delta (k^2m_2^2)`$ (A.132)
$`\times {\displaystyle _1^1}d\mathrm{cos}\theta \delta (p^2+2p_0k_02|𝐩||𝐤|\mathrm{cos}\theta +k^2m_1^2).`$
$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ receives contribution from $`k`$ which satisfies

$$\frac{|p^2+2p_0k_0+m_2^2m_1^2|}{2|𝐩||𝐤|}<1$$
(A.133)
and $`k_0^2=𝐤^2+m_2^2`$. We make a square of Eq.(A.133) and obtain

$`4p_0k_0(p^2+m_2^2m_1^2)<(p^2+m_2^2m_1^2)^24p_0^2m_2^24p^2|𝐤|^2.`$ (A.134)
There are three cases when Eq.(A.134) holds:

1. LHS$`>0`$ and RHS$`>0`$, (RHS)$`{}_{}{}^{2}`$(LHS)$`{}_{}{}^{2}>0`$.
2. LHS$`<0`$ and RHS$`<0`$, (RHS)$`{}_{}{}^{2}`$(LHS)$`{}_{}{}^{2}<0`$.
3. LHS$`<0`$ and RHS$`>0`$.

First, we calculate (RHS)$`{}_{}{}^{2}`$(LHS)<sup>2</sup>.

$`(\mathrm{RHS})^2(\mathrm{LHS})^2`$ (A.135)
$`=`$ $`16p^4|𝐤|^4+\{8p^2[(p^2+m_2^2m_1^2)^2+4p_0^2m_2^2]16(p^2+m_2^2m_1^2)^2p_0^2]|𝐤|^2`$ (A.137)
$`+\left[(p^2+m_2^2m_1^2)^2+4p_0^2m_2^2\right]^216(p^2+m_2^2m_1^2)^2p_0^2m_2^2`$
$`=`$ $`16p^4|𝐤|^4+8\left[(p^22p_0^2)(p^2+m_2^2m_1^2)^2+4p^2p_0^2m_2^2\right]|𝐤|^2`$ (A.139)
$`+\left[(p^2+m_2^2m_1^2)^24p_0^2m_2^2\right]^2`$
Whether there exists $`|𝐤|^2`$ which satisfies $`(\mathrm{RHS})^2(\mathrm{LHS})^2<0`$ depends on the sign of the following expression:

$`D`$ $`=`$ $`16\left[(p^22p_0^2)(p^2+m_2^2m_1^2)^2+4p^2p_0^2m_2^2\right]^216p^4\left[(p^2+m_2^2m_1^2)^24p_0^2m_2^2\right]^2`$ (A.140)
$`=`$ $`64p_0^2|𝐩|^2(p^2+m_2^2m_1^2)^2\left[p^2(m_1+m_2)^2\right]\left[p^2(m_1m_2)^2\right].`$ (A.141)
Thus, for $`(m_1m_2)^2<p^2<(m_1+m_2)^2`$

$$(\mathrm{RHS})^2(\mathrm{LHS})^2>0,$$

for $`p^2<(m_1m_2)^2`$ or $`(m_1+m_2)^2<p^2`$

$`|𝐤|^2<|𝐤|_{}^2\mathrm{or}|𝐤|^2>|𝐤|_+^2(\mathrm{RHS})^2(\mathrm{LHS})^2>0,`$
$`|𝐤|_{}^2<|𝐤|^2<|𝐤|_+^2(\mathrm{RHS})^2(\mathrm{LHS})^2<0,`$
where $`|𝐤|_\pm ^2`$ are given by

$`|𝐤|_\pm ^2`$ $`=`$ $`{\displaystyle \frac{1}{4}}\left\{\sqrt{\left(1+{\displaystyle \frac{m_2^2m_1^2}{p^2}}\right)^2}|𝐩|\pm \sqrt{\left[1{\displaystyle \frac{(m_2+m_1)^2}{p^2}}\right]\left[1{\displaystyle \frac{(m_1m_2)^2}{p^2}}\right]}|p_0|\right\}^2.`$ (A.143)
$`\mathrm{?}`$
Secondly, for $`p^2>m_1^2m_2^2`$

$$p_0k_00\mathrm{LHS}0,$$

for $`p^2<m_1^2m_2^2`$

$$p_0k_00\mathrm{LHS}0.$$

Thirdly, for $`p^2>0`$

$$\mathrm{RHS}<0$$

for $`p^2<0`$

$$|𝐤|^2\frac{(p^2+m_2^2m_1^2)^2+4p_0^2m_2^2}{4p^2}|𝐤|_0^2\mathrm{RHS}0$$

Therefore, the conditions for the above three cases are respectively given by

1. $`p^2<0`$, $`p_0k_0(p^2+m_2^2m_1^2)>0`$ and $`|𝐤|^2>|𝐤|_+^2`$.
2. $`p^2>(m_1+m_2)^2`$, $`p_0k_0(p^2+m_2^2m_1^2)<0`$ and $`|𝐤|_{}^2<|𝐤|^2<|𝐤|_+^2`$,

$`0<p^2<(m_1m_2)^2`$ and $`p_0k_0(p^2+m_2^2m_1^2)<0`$ and $`|𝐤|_{}^2<|𝐤|^2<|𝐤|_+^2`$,

$`p^2<0`$, $`p_0k_0(p^2+m_2^2m_1^2)<0`$ and $`|𝐤|_{}^2<|𝐤|^2<|𝐤|_0^2`$.
3. $`p^2<0`$, $`p_0k_0(p^2+m_2^2m_1^2)<0`$ and $`|𝐤|^2>|𝐤|_0^2`$.

The case C can be combined with the third of the case B as

1. $`p^2<0`$, $`p_0k_0(p^2+m_2^2m_1^2)<0`$ and $`|𝐤|^2>|𝐤|_{}^2`$.

1. For $`p^2>(m_1+m_2)^2`$ or $`0<p^2<(m_1m_2)^2`$,

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}{\displaystyle _{|𝐤|_{}}^{|𝐤|_+}}d|𝐤|{\displaystyle \frac{|𝐤|}{\omega _k}}n(\omega _k)`$ (A.144)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}{\displaystyle _\omega _{}^{\omega _+}}d\omega n(\omega )`$ (A.145)
$`=`$ $`{\displaystyle \frac{1}{16\pi |𝐩|}}{\displaystyle \frac{1}{\beta }}\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta \omega _+}}{1\mathrm{e}^{\beta \omega _{}}}}\right|.`$ (A.146)
2. For $`(m_1m_2)^2<p^2<(m_1+m_2)^2`$,

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`0.`$ (A.147)
3. For $`p^2<0`$,

$`\mathrm{Im}iF^{(2)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}\left({\displaystyle _{|𝐤|_{}}^{\mathrm{}}}+{\displaystyle _{|𝐤|_+}^{\mathrm{}}}\right)d|𝐤|{\displaystyle \frac{|𝐤|}{\omega _k}}n(\omega _k)`$ (A.148)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}\left({\displaystyle _\omega _{}^{\mathrm{}}}+{\displaystyle _{\omega _+}^{\mathrm{}}}\right)\mathrm{d}\omega n(\omega )`$ (A.149)
$`=`$ $`{\displaystyle \frac{1}{16\pi |𝐩|}}{\displaystyle \frac{1}{\beta }}\mathrm{ln}\left|(1\mathrm{e}^{\beta \omega _+})(1\mathrm{e}^{\beta \omega _{}})\right|.`$ (A.150)

where

$$\omega _\pm =\frac{1}{2}\left|\sqrt{\left(1+\frac{m_2^2m_1^2}{p^2}\right)^2}|p_0|\pm \sqrt{\left[1\frac{(m_2+m_1)^2}{p^2}\right]\left[1\frac{(m_1m_2)^2}{p^2}\right]}|𝐩|\right|$$
(A.151)
Similarly,

1. For $`p^2>(m_1+m_2)^2`$,

$`\mathrm{Im}iF^{(3)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}{\displaystyle _{|𝐤|_{}}^{|𝐤|_+}}d|𝐤|{\displaystyle \frac{|𝐤|}{\omega _k}}n(\omega _k)n(|p_0|\omega _k)`$ (A.152)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}{\displaystyle _\omega _{}^{\omega _+}}d\omega n(\omega )n(|p_0|\omega )`$ (A.153)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|\beta }}{\displaystyle \frac{1}{\mathrm{e}^{\beta |p_0|}1}}\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta \omega _+}}{1\mathrm{e}^{\beta \omega _{}}}}{\displaystyle \frac{\mathrm{e}^{\beta (|p_0|\omega _{})}1}{\mathrm{e}^{\beta (|p_0|\omega _+)}1}}\right|.`$ (A.154)
2. For $`(m_1m_2)^2<p^2<(m_1+m_2)^2`$,

$`\mathrm{Im}iF^{(3)}(p;m_1,m_2)`$ $`=`$ $`0.`$ (A.155)
3. For $`0<p^2<(m_1m_2)^2`$,

$`\mathrm{Im}iF^{(3)}(p;m_1,m_2)`$ $`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}{\displaystyle _{|𝐤|_{}}^{|𝐤|_+}}d|𝐤|{\displaystyle \frac{|𝐤|}{\omega _k}}n(\omega _k)n(|p_0|+\omega _k)`$ (A.156)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}{\displaystyle _\omega _{}^{\omega _+}}d\omega n(\omega )n(|p_0|+\omega )`$ (A.157)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|\beta }}{\displaystyle \frac{1}{\mathrm{e}^{\beta |p_0|}1}}`$ (A.159)
$`\times \left[\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta \omega _+}}{1\mathrm{e}^{\beta \omega _{}}}}\right|\mathrm{e}^{\beta |p_0|}\mathrm{ln}\left|{\displaystyle \frac{1\mathrm{e}^{\beta (|p_0|+\omega _+)}}{1\mathrm{e}^{\beta (|p_0|+\omega _{})}}}\right|\right].`$
4. For $`p^2<0`$,

$`\mathrm{Im}iF^{(3)}(p;m_1,m_2)`$ (A.160)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}\left({\displaystyle _{|𝐤|_{}}^{\mathrm{}}}d|𝐤|{\displaystyle \frac{|𝐤|}{\omega _k}}n(\omega _k)n(|p_0|+\omega _k)+{\displaystyle _{|𝐤|_+}^{\mathrm{}}}d|𝐤|{\displaystyle \frac{|𝐤|}{\omega _k}}n(\omega _k)n(|p_0|\omega _k)\right)`$ (A.161)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|}}\left({\displaystyle _\omega _{}^{\mathrm{}}}d\omega n(\omega )n(|p_0|+\omega )+{\displaystyle _{\omega _+}^{\mathrm{}}}d\omega n(\omega )n(|p_0|\omega )\right)`$ (A.162)
$`=`$ $`{\displaystyle \frac{1}{8\pi |𝐩|\beta }}\{{\displaystyle \frac{1}{e^{\beta |p_0|}1}}[\mathrm{ln}|1\mathrm{e}^{\beta \omega _+}|+\mathrm{e}^{\beta |p_0|}\mathrm{ln}|1\mathrm{e}^{\beta (\omega _+|p_0|)}|]`$ (A.164)
$`+{\displaystyle \frac{1}{e^{\beta |p_0|}1}}[\mathrm{ln}|1\mathrm{e}^{\beta \omega _{}}|+\mathrm{e}^{\beta |p_0|}\mathrm{ln}|1\mathrm{e}^{\beta (\omega _{}+|p_0|)}|]\}.`$

## Appendix B Physical meaning of $`|𝐤|_\pm `$

In this appendix we explain the physical meaning of $`|𝐤|_\pm `$ defined by Eq.(A.143).

Let us first consider the case $`p^2>(m_1+m_2)^2`$ or $`0<p^2<(m_1m_2)^2`$. For simplicity we suppose $`p_0>0`$. Since internal particles are on shell, in the center-of-mass frame we have

$$p^2=\left(\sqrt{m_1^2+|𝐤|^2}\pm \sqrt{m_2^2+|𝐤|^2}\right)^2,$$
(B.165)
where $`+`$ is for $`p^2>(m_1+m_2)^2`$ and $``$ is for $`0<p^2<(m_1m_2)^2`$. In either case

$$|𝐤|^2=\frac{[p^2(m_1+m_2)^2][p^2(m_1m_2)^2]}{4p^2}.$$
(B.166)
If we boost the system to the positive z-direction, the external and internal momenta become

$`\{\begin{array}{cc}p_0^{}=p_0\mathrm{cosh}\theta ,\hfill & \\ 𝐩_{}=0,p_z^{}=p_0\mathrm{sinh}\theta \hfill & (\theta >0),\hfill \end{array}`$ (B.169)
$`\{\begin{array}{c}k_0^{}=k_0\mathrm{cosh}\theta +k_z\mathrm{sinh}\theta ,\hfill \\ 𝐤_{}^{}{}_{}{}^{}=𝐤_{},k_z^{}=k_z\mathrm{cosh}\theta +k_0\mathrm{sinh}\theta .\hfill \end{array}`$ (B.172)
Since $`|𝐤|<k_z<|𝐤|`$,

$$|𝐤^{}|_{}^2|𝐤^{}|^2|𝐤^{}|_+^2,$$
(B.173)
where

$`|𝐤^{}|_\pm ^2`$ $`=`$ $`(\pm |𝐤|\mathrm{cosh}\theta +|k_0|\mathrm{sinh}\theta )^2`$ (B.174)
$`=`$ $`\left(\pm |𝐤|{\displaystyle \frac{p_0^{}}{p_0}}+|k_0|{\displaystyle \frac{p_z^{}}{p_0}}\right)^2`$ (B.175)
$`=`$ $`{\displaystyle \frac{1}{4}}\left\{\pm \sqrt{\left[1{\displaystyle \frac{(m_1+m_2)^2}{p^2}}\right]\left[1{\displaystyle \frac{(m_1m_2)^2}{p^2}}\right]}p_0^{}+\left|1+{\displaystyle \frac{m_2^2m_1^2}{p^2}}\right||𝐩^{}|\right\}^2.`$ (B.176)
This is nothing but Eq.(A.143). Therefore, for $`p^2>(m_1+m_2)^2`$ or $`0<p^2<(m_1m_2)^2`$, $`|𝐤|_+`$ and $`|𝐤|_{}`$ are respectively the maximum and minimum values of $`|𝐤|`$ such that the internal particles are on-shell.

Let us next consider the case $`p^2<0`$. In the Breit frame, $`p=(0,𝐩)`$, with $`𝐩`$ in the positive z-direction we write the momentum of the particle with mass $`m_2`$ by $`k=(k_0,𝐤)`$ where $`k_0`$ can be positive or negative. Then, the momentum of the particle with mass $`p+k=(k_0,𝐩+𝐤)`$. Since internal particles are on shell, we have

$$k_0^2=m_2^2+k_z^2+|𝐤_{}|^2=m_1^2+(p_z+k_z)^2+|𝐤_{}|^2.$$
(B.177)
From this equation we obtain

$`k_z`$ $`=`$ $`{\displaystyle \frac{p_z^2+m_1^2m_2^2}{2p_z}}={\displaystyle \frac{p^2+m_1^2m_2^2}{2\sqrt{p^2}}},`$ (B.178)
$`k_0^2`$ $`=`$ $`{\displaystyle \frac{[p^2+(m_1+m_2)^2][p^2+(m_1m_2)^2]}{4p^2}}^2+|𝐤_{}|^2,`$ (B.179)
When boosting the frame to the z-direction, we suppose the momenta become

$`p^{}=(p_0^{},𝐩^{}),k^{}=(k_0^{},𝐤^{}),p^{}+k^{}=(p_0^{}+k_0^{},𝐩^{}+𝐤^{}).`$ (B.180)
Then, $`p^{}`$ and $`k^{}`$ are given by

$`\{\begin{array}{c}p_0^{}=|𝐩|\mathrm{sinh}\theta ,\hfill \\ |𝐩^{}|=|𝐩|\mathrm{cosh}\theta ,\hfill \end{array}`$ (B.183)
$`\{\begin{array}{c}k_0^{}=k_0\mathrm{cosh}\theta +k_z\mathrm{sinh}\theta ,\hfill \\ 𝐤_{}^{}{}_{}{}^{}=𝐤_{},k_z^{}=k_z\mathrm{cosh}\theta +k_0\mathrm{sinh}\theta .\hfill \end{array}`$ (B.186)
Using Eqs.(B.179), (B.183) and (B.186), we obtain

$$|𝐤^{}|^2=(k_z\mathrm{cosh}\theta +k_0\mathrm{sinh}\theta )^2+|𝐤_{}|^2.$$
(B.187)
$`|𝐤^{}|^2`$ has a minimum when $`|𝐤_{}|^2=0`$ but not a maximum. Depending on the sign of $`k_0`$, the minimum value is given by

$`|𝐤^{}|_\pm ^2`$ $`=`$ $`\left(|k_z|{\displaystyle \frac{|𝐩^{}|}{\sqrt{p^2}}}\pm \sqrt{k_z^2+m_2^2}{\displaystyle \frac{p_0^{}}{\sqrt{p^2}}}\right)^2`$ (B.188)
$`=`$ $`{\displaystyle \frac{1}{4}}\left\{\left|1{\displaystyle \frac{m_1^2m_2^2}{p^2}}\right||𝐩^{}|\pm \sqrt{\left[1+{\displaystyle \frac{(m_1+m_2)^2}{p^2}}\right]\left[1+{\displaystyle \frac{(m_1m_2)^2}{p^2}}\right]}p_0^{}\right\}^2.`$ (B.189)
Eq.(B.189) is nothing but Eq.(A.143). Therefore, in the case of $`p^2<0`$, both of $`|𝐤|_+`$ and $`|𝐤|_{}`$ are the minimum values of $`|𝐤|`$ such that the internal particles can be on-shell.

References

1. S. Chiku and T. Hatsuda, Phys. Rev. D57, R6 (1998); Phys. Rev. D58, 76001 (1998)
2. J. Schwinger, J. Math. Phys. 2, 407 (1961).
3. L.V. Keldysh, JETP(Sov. Phys.) 20, 1018 (1964).
4. R.A. Craig, J. Math. Phys. 9, 605 (1968).
5. G.-Z. Zhou, Z.-B. Su, B.-L Hao and Lu Yu, Phys. Rev. B22, 3385 (1980).
6. Y. Takahashi and H. Umezawa, Coll. Phen. 2, 55 (1975).
7. G.W. Semenoff and Y. Takahashi, Nucl. Phys. B220\[FS8\], 196 (1983).
8. l.P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, New York 1969); D.C. Langreth, in 1975 NATO ADI on Linear and Non-Linear Electron Transport in Solids, J.T. Devreese and E. van Boem eds. (Plenum, New York 1976)
9. R. Mills, Propagators for many-particle systems (Gordon and Breach, New york 1969)
10. E.M. Lifshitz and L.P. Pitaevskii, Course of theoretical physics, vol.10: Physical kinetics (Pergamon, New York 1981).
11. H. Umezawa, H. Matsumoto and M. Tachiki, Thermofield Dynamics and Condensed States (North-Holland, Amsterdam, 1982).
12. A.J. Niemi and G.W. Semenoff, Ann. Phys. 152, 105 (1984).
13. A.J. Niemi and G.W. Semenoff, Nucl. Phys. B230, 181 (1984).
14. M. Lebellac, Thermal Field Theory, Cambridge University Press 1996.
15. G.W. Semenoff and H. Umezawa, Nucl. Phys. B220, 196 (1983)
16. See, for example, P. Ramond, Field Theory, A Modern Primer (Benjamin/Cummings, New York, 1981).
17. Y. Fujimoto, M. Morikawa and M. Sasaki, Phys. Rev. D33, 590 (1986).
18. H.A. Weldon, Phys. Rev. D28, 2007 (1983).
19. P. Post and J.B. Tausk, Mod. Phys. Lett. A11, 2115 (1996).
20. J. Gasser and M.E. Sainio, Eur. Phys. J. C6, 297 (1999).
21. J.O. Andersen and E. Braaten, Phys. Rev. D62, 45004 (2000); Phys. Rev. D51, 6990 (1995).
22. S. Groote, J.G. Koelner and A.A. Pivovarov, Phys. Lett. B443, 269 (1998); Nucl. Phys. B542, 515 (1999); Eur. Phys. J. C11, 279 (1999).
23. S. Weinberg, Phys. Rev. D9, 3357 (1974);
L. Dolan and R. Jackiw, Phys. Rev. D9, 3320 (1974).
24. D.A. Kirzhnits and A.D. Linde, Ann. Phys. (N.Y.)101, 195 (1976).
25. T. Nishikawa, O. Morimatsu, Y. Hidaka and M. Ohtani, in preparation.
26. T. Nishikawa, O. Morimatsu and Y. Hidaka, in preparation.