text:
# Vector Manifestation and Violation of Vector Dominance in Hot Matter

## 1 Introduction

Spontaneous chiral symmetry breaking is one of the most important features in low energy QCD. The phenomena of the light pseudoscalar mesons (the pion and its flavor partner), which are regarded as the approximate Nambu-Goldstone bosons associated with the symmetry breaking, are well described by the symmetry property in the low energy region. This chiral symmetry is expected to be restored in hot and/or dense matter and properties of hadrons will be changed near the critical point of the chiral symmetry restoration HatsudaKunihiro ; Pisarski:95 ; Brown-Rho:96 ; Brown:2001nh ; HatsudaShiomiKuwabara ; Rapp-Wambach:00 ; Wilczek . It is important to investigate the physics near the phase transition point both theoretically and experimentally. In fact, the CERN Super Proton Synchrotron (SPS) observed an enhancement of dielectron ($`e^+e^{}`$) mass spectra below the $`\rho /\omega `$ resonance Agakishiev:1995xb . This can be explained by the dropping mass of the $`\rho `$ meson (see, e.g., Refs. Li:1995qm ; Brown-Rho:96 ; Rapp-Wambach:00 ) following the Brown-Rho scaling proposed in Ref. BR . Furthermore, the Relativistic Heavy Ion Collider (RHIC) has started to measure several physical processes in hot matter which include the dilepton energy spectra. This will further clarify the properties of vector mesons in hot matter. Therefore it is interesting to study the properties of the pion and vector mesons near the critical temperature, especially associated with the dropping mass of the vector meson.

There is no strong restriction concerning vector meson masses in the standard scenario of chiral symmetry restoration, where the pion joins with the scalar meson in the same chiral representation \[see for example, Refs. HatsudaKunihiro ; Bochkarev:1995gi ; PT:96 \]. However there is a scenario which certainly requires the dropping mass of the vector meson and supports the Brown-Rho scaling: In Ref. HYb , it was proposed that there can be another possibility for the pattern of chiral symmetry restoration, named the vector manifestation (VM). The VM was proposed as a novel manifestation of the chiral symmetry in the Wigner realization, in which the chiral symmetry is restored by the massless degenerate pion (and its flavor partners) and the longitudinal $`\rho `$ meson (and its flavor partners) as the chiral partner. In terms of the chiral representations of the low-lying mesons, there is a representation mixing in the vacuum HYb ; HYc . When we approach the critical point, there are two possibilities for the pattern of the chiral symmetry restoration: One possible pattern is the standard scenario and another is the VM. Both of them are on an equal footing with each other in terms of the chiral representations. It is worthwhile to study the physics associated with the VM as well as that with the standard scenario of the chiral symmetry restoration.

It has been shown that the VM is formulated at a large number of flavor HYb , critical temperature HSasaki  and critical density HKR  by using the effective field theory including both pions and vector mesons based on the hidden local symmetry (HLS) BKUYY ; BKYa , where a second order or weakly first order phase transition was assumed. In the VM at finite temperature and/or density, the intrinsic temperature and/or density dependences of the parameters of the HLS Lagrangian played important roles to realize the chiral symmetry restoration consistently: In the framework of the HLS the equality between the axial vector and vector current correlators at the critical point can be satisfied only if the intrinsic thermal and/or density effects are included. The intrinsic temperature and/or density dependences are nothing but the information converted from the underlying QCD through the Wilsonian matching HYa ; HSasaki . In Ref. HKRS , the predictions of the VM were made on the pion velocity and the vector and axial vector susceptibilities. It was stressed that the vector mesons in addition to the pseudoscalar mesons become relevant degrees of freedom near the critical temperature.

In this paper we shed some light on the validity of the vector dominance (VD) in hot matter. In several analyses such as the one on the dilepton spectra in hot matter carried out in Ref. Rapp-Wambach:00 , the VD is assumed to be held even in the high temperature region. There are several analyses resulting in the dropping mass consistent with the VD, as shown in Refs. GomezNicola:2002tn ; Dobado:2002xf . On the other hand, the analysis done in Ref. Pisarski  shows that the thermal vector meson mass goes up if the VD holds. Thus, it is interesting to study what the VM predicts on the VD. In the present analysis we present a new prediction of the VM in hot matter on the direct photon-$`\pi `$-$`\pi `$ coupling which measures the validity of the VD of the electromagnetic form factor of the pion. We find that the VM predicts a large violation of the VD at the critical temperature. This indicates that the assumption of the VD may need to be weakened, at least in some amounts, for consistently including the effect of the dropping mass of the vector meson.

This paper is organized as follows: In section 2, we show the HLS Lagrangian which we use in the present analysis. In section 3, we present an entire list of the hadronic thermal corrections to the two-point functions including the vector and pseudoscalar meson loop contributions in the background field gauge. We think that it is useful to summarize the entire list of the hadronic thermal corrections since only the relevant part of them was listed in Ref. HKRS . In section 4, we first give an account of the general idea of the intrinsic thermal and/or dense effects, which was not explicitly stated before. Then, following Ref. HKRS , we briefly review how to extend the Wilsonian matching to the version at non-zero temperature to incorporate the intrinsic thermal effect. In section 5, we show the explicit forms of the pole masses of the longitudinal and transverse modes of the vector meson at non-zero temperature. Including the intrinsic thermal effect through the Wilsonian matching, we show how the VM is formulated at the critical temperature. In section 6, we study the temperature dependences of the temporal and spatial pion decay constants and the pion velocity in the low temperature region. The estimation of the value of the critical temperature is also carried out in wider range of input parameters than the one used in Ref. HSasaki . In section 7, we study the validity of the VD in hot matter, and show that the VD is largely violated near critical temperature. In section 8, we give a summary and discussions. Several functions and formulas used in this paper are listed in Appendices A-D.

## 2  Hidden Local Symmetry

In this section, we show the Lagrangian based on the hidden local symmetry (HLS) which we use in the present analysis.

The HLS model is based on the $`G_{\mathrm{global}}\times H_{\mathrm{local}}`$ symmetry, where $`G=SU(N_f)_L\times SU(N_f)_R`$ is the chiral symmetry and $`H=SU(N_f)_V`$ is the HLS. The basic quantities are the HLS gauge boson and two matrix valued variables $`\xi _L(x)`$ and $`\xi _R(x)`$ which transform as

$$\xi _{L,R}(x)\xi _{L,R}^{}(x)=h(x)\xi _{L,R}(x)g_{L,R}^{},$$
(1)
where $`h(x)H_{\mathrm{local}}\text{and}g_{L,R}[\text{SU}(N_f)_{\mathrm{L},\mathrm{R}}]_{\mathrm{global}}`$. These variables are parameterized as

$$\xi _{L,R}(x)=e^{i\sigma (x)/F_\sigma }e^{i\pi (x)/F_\pi },$$
(2)
where $`\pi =\pi ^aT_a`$ denotes the pseudoscalar Nambu-Goldstone bosons associated with the spontaneous symmetry breaking of $`G_{\mathrm{global}}`$ chiral symmetry, and $`\sigma =\sigma ^aT_a`$ denotes the Nambu-Goldstone bosons associated with the spontaneous breaking of $`H_{\mathrm{local}}`$. This $`\sigma `$ is absorbed into the HLS gauge boson through the Higgs mechanism. $`F_\pi \text{and}F_\sigma `$ are the decay constants of the associated particles. The phenomenologically important parameter $`a`$ is defined as

$$a=\frac{F_{\sigma }^{}{}_{}{}^{2}}{F_{\pi }^{}{}_{}{}^{2}}.$$
(3)
The covariant derivatives of $`\xi _{L,R}`$ are given by

$`D_\mu \xi _L`$ $`=`$ $`_\mu \xi _LiV_\mu \xi _L+i\xi _L_\mu ,`$
$`D_\mu \xi _R`$ $`=`$ $`_\mu \xi _RiV_\mu \xi _R+i\xi _R_\mu ,`$ (4)
where $`V_\mu `$ is the gauge field of $`H_{\mathrm{local}}`$, and $`_\mu \text{and}_\mu `$ are the external gauge fields introduced by gauging the $`G_{\mathrm{global}}`$ symmetry.

The HLS Lagrangian with the lowest derivative terms in the chiral limit is given by BKUYY ; BKYa

$$_{(2)}=F_{\pi }^{}{}_{}{}^{2}\text{tr}\left[\widehat{\alpha }_\mu \widehat{\alpha }_{}^\mu \right]+F_{\sigma }^{}{}_{}{}^{2}\text{tr}\left[\widehat{\alpha }_\mu \widehat{\alpha }_{}^\mu \right]\frac{1}{2g^2}\text{tr}\left[V_{\mu \nu }V^{\mu \nu }\right],$$
(5)
where $`g`$ is the HLS gauge coupling, $`V_{\mu \nu }`$ is the field strength of $`V_\mu `$ and

$`\widehat{\alpha }_{}^\mu `$ $`=`$ $`{\displaystyle \frac{1}{2i}}\left[D^\mu \xi _R\xi _R^{}D^\mu \xi _L\xi _L^{}\right],`$
$`\widehat{\alpha }_{}^\mu `$ $`=`$ $`{\displaystyle \frac{1}{2i}}\left[D^\mu \xi _R\xi _R^{}+D^\mu \xi _L\xi _L^{}\right].`$ (6)
In the HLS, as first pointed in Ref. Georgi  and developed further in Refs. HYd ; Tana ; HYa ; HYc , thanks to the gauge symmetry, a systematic loop expansion can be formally performed with the vector mesons included in addition to the pseudoscalar mesons. In this chiral perturbation theory (ChPT) with HLS the vector meson mass is considered as small compared with the chiral symmetry breaking scale $`\mathrm{\Lambda }_\chi `$, by assigning $`𝒪(p)`$ to the HLS gauge coupling Georgi ; Tana :

$$g𝒪(p).$$
(7)
According to the entire list shown in Ref. Tana , there are 35 counter terms at $`𝒪(p^4)`$ for general $`N_f`$. However, only three terms are relevant when we consider two-point functions in the chiral limit:

$$_{(4)}=z_1\text{tr}\left[\widehat{𝒱}_{\mu \nu }\widehat{𝒱}^{\mu \nu }\right]+z_2\text{tr}\left[\widehat{𝒜}_{\mu \nu }\widehat{𝒜}^{\mu \nu }\right]+z_3\text{tr}\left[\widehat{𝒱}_{\mu \nu }V^{\mu \nu }\right],$$
(8)
where

$`\widehat{𝒜}_{\mu \nu }={\displaystyle \frac{1}{2}}\left[\xi _R_{\mu \nu }\xi _R^{}\xi _L_{\mu \nu }\xi _L^{}\right],`$ (9)
$`\widehat{𝒱}_{\mu \nu }={\displaystyle \frac{1}{2}}\left[\xi _R_{\mu \nu }\xi _R^{}+\xi _L_{\mu \nu }\xi _L^{}\right],`$ (10)
with $`_{\mu \nu }\text{and}_{\mu \nu }`$ being the field strengths of $`_\mu \text{and}_\mu `$.

In Ref. HYa , based on the ChPT with HLS, the Wilsonian matching was proposed, through which the parameters of the HLS Lagrangian are determined by the underlying QCD at the matching scale $`\mathrm{\Lambda }`$. (For a review of the ChPT with HLS, the Wilsonian matching and VM at zero temperature, see Ref. HYc .) The Wilsonian matching at $`T=0`$ with $`N_f=3`$, where the validity of the ChPT with HLS is essential, was shown to give several predictions in remarkable agreement with experiments HYa ; HYc . This strongly suggests that the ChPT with HLS is valid even numerically. In Ref. HSasaki , we applied the ChPT with HLS combined with the Wilsonian matching to finite temperature. There the expansion parameter is $`T/F_\pi (\mathrm{\Lambda })`$ instead of $`T/F_\pi (0)`$ used in the standard ChPT. Since $`F_\pi (\mathrm{\Lambda })>F_\pi (0)`$, we think that the present formalism can be applied in the higher temperature region than the standard ChPT.

## 3 Two-Point Functions in Background Field Gauge

In the present approach hadronic thermal effects are included by calculating loop contributions of the pseudoscalar and vector meson. In Ref. HSasaki  we used hadronic thermal corrections to the pion decay constant and vector meson mass calculated within the framework of the chiral perturbation theory with HLS in the Landau gauge. While in Ref. HKRS  hadronic thermal corrections were calculated in the background field gauge. (For the application of the background field gauge to the HLS, see, e.g., Ref. HYc .) In this section we show details of the calculation of the hadronic thermal corrections to the two-point functions in the background field gauge.

Let us consider the loop corrections to the two-point functions at non-zero temperature. As was done in Ref. HKRS  (see also section 4), we neglect the possible Lorentz symmetry violating effects caused by the intrinsic temperature dependences of the bare parameters in the present analysis. So we calculate loop corrections at non-zero temperature from the bare Lagrangian with Lorentz invariance shown in Eq. (5). For this purpose it is convenient to introduce the following functions:

$`A_0(M^2;T)`$ $``$ $`T{\displaystyle \underset{n=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \frac{d^3k}{(2\pi )^3}\frac{1}{M^2k^2}},`$ (11)
$`B_0(p_0,\overline{p};M_1,M_2;T)`$ $``$ $`T{\displaystyle \underset{n=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \frac{d^3k}{(2\pi )^3}\frac{1}{[M_1^2k^2][M_2^2(kp)^2]}},`$ (12)
$`B^{\mu \nu }(p_0,\overline{p};M_1,M_2;T)`$ $``$ $`T{\displaystyle \underset{n=\mathrm{}}{\overset{\mathrm{}}{}}}{\displaystyle \frac{d^3k}{(2\pi )^3}\frac{\left(2kp\right)^\mu \left(2kp\right)^\nu }{[M_1^2k^2][M_2^2(kp)^2]}},`$ (13)
where $`\overline{p}=|\stackrel{}{p}|`$ and the $`0`$th component of the loop momentum is taken as $`k^0=i2n\pi T`$, while that of the external momentum is taken as $`p^0=i2n^{}\pi T`$ \[$`n,n^{}`$: integer\]. Using the standard formula (see, e.g., Ref. Kap ), these functions are divided into two parts as

$`A_0(M;T)`$ $`=`$ $`A_0^{(\mathrm{vac})}(M)+\overline{A}_0(M;T),`$ (14)
$`B_0(p_0,\overline{p};M_1,M_2;T)`$ $`=`$ $`B_0^{(\mathrm{vac})}(p;M_1,M_2)+\overline{B}_0(p_0,\overline{p};M_1,M_2;T),`$ (15)
$`B^{\mu \nu }(p_0,\overline{p};M_1,M_2;T)`$ $`=`$ $`B^{(\mathrm{vac})\mu \nu }(p;M_1,M_2)+\overline{B}^{\mu \nu }(p_0,\overline{p};M_1,M_2;T),`$ (16)
where $`A_0^{(\mathrm{vac})}`$, $`B_0^{(\mathrm{vac})}`$ and $`B^{(\mathrm{vac})\mu \nu }`$ express the quantum corrections given in Eqs. (C.30)–(C.32), and $`\overline{A}_0`$, $`\overline{B}_0`$ and $`\overline{B}^{\mu \nu }`$ the hadronic thermal corrections. We summarize the explicit forms of the functions $`\overline{A}_0`$, $`\overline{B}_0`$ and $`\overline{B}^{\mu \nu }`$ in various limits relevant to the present analysis in Appendix B. Note that the $`0`$th component of the momentum $`p_0`$ in the right-hand-sides of the above expressions is analytically continued to the Minkowski variable: $`p_0`$ is understood as $`p_0+iϵ`$ ($`ϵ+0`$) for the retarded function and $`p_0iϵ`$ for the advanced function. Then, $`A_0^{(\mathrm{vac})}`$, $`B_0^{(\mathrm{vac})}`$ and $`B^{(\mathrm{vac})\mu \nu }`$ have no explicit temperature dependence, while they have intrinsic temperature dependence which is introduced through the Wilsonian matching as we will see later. In the following we write the two-point functions of $`\overline{𝒜}_\mu `$-$`\overline{𝒜}_\nu `$, $`\overline{𝒱}_\mu `$-$`\overline{𝒱}_\nu `$, $`\overline{V}_\mu `$-$`\overline{V}_\nu `$ and $`\overline{V}_\mu `$-$`\overline{𝒱}_\nu `$ as $`\mathrm{\Pi }_{}^{\mu \nu }`$, $`\mathrm{\Pi }_{}^{\mu \nu }`$, $`\mathrm{\Pi }_V^{\mu \nu }`$ and $`\mathrm{\Pi }_V^{\mu \nu }`$, respectively.

At non-zero temperature there are four independent polarization tensors, which we choose as defined in Appendix A. We decompose the two-point functions $`\mathrm{\Pi }_{}^{\mu \nu }`$, $`\mathrm{\Pi }_V^{\mu \nu }`$ $`\mathrm{\Pi }_V^{\mu \nu }`$ and $`\mathrm{\Pi }_{}^{\mu \nu }`$ as

$`\mathrm{\Pi }_{}^{\mu \nu }`$ $`=`$ $`u^\mu u^\nu \mathrm{\Pi }_{}^t+(g^{\mu \nu }u^\mu u^\nu )\mathrm{\Pi }_{}^s+P_L^{\mu \nu }\mathrm{\Pi }_{}^L+P_T^{\mu \nu }\mathrm{\Pi }_{}^T,`$ (17)
$`\mathrm{\Pi }_V^{\mu \nu }`$ $`=`$ $`u^\mu u^\nu \mathrm{\Pi }_V^t+(g^{\mu \nu }u^\mu u^\nu )\mathrm{\Pi }_V^s+P_L^{\mu \nu }\mathrm{\Pi }_V^L+P_T^{\mu \nu }\mathrm{\Pi }_V^T,`$ (18)
$`\mathrm{\Pi }_V^{\mu \nu }`$ $`=`$ $`u^\mu u^\nu \mathrm{\Pi }_V^t+(g^{\mu \nu }u^\mu u^\nu )\mathrm{\Pi }_V^s+P_L^{\mu \nu }\mathrm{\Pi }_V^L+P_T^{\mu \nu }\mathrm{\Pi }_V^T,`$ (19)
$`\mathrm{\Pi }_{}^{\mu \nu }`$ $`=`$ $`u^\mu u^\nu \mathrm{\Pi }_{}^t+(g^{\mu \nu }u^\mu u^\nu )\mathrm{\Pi }_{}^s+P_L^{\mu \nu }\mathrm{\Pi }_{}^L+P_T^{\mu \nu }\mathrm{\Pi }_{}^T,`$ (20)
where $`u^\mu =(1,\stackrel{}{0})`$ and $`P_L^{\mu \nu }`$ and $`P_T^{\mu \nu }`$ are defined in Eq. (A.1). Similarly, we decompose the function $`B^{\mu \nu }`$ as

$$B^{\mu \nu }=u^\mu u^\nu B^t+(g^{\mu \nu }u^\mu u^\nu )B^s+P_L^{\mu \nu }B^L+P_T^{\mu \nu }B^T.$$
(21)
Furthermore, similarly to the division of the functions into one part for expressing the quantum correction and another for the hadronic thermal correction done in Eqs. (14)–(16), we divide the two-point functions into two parts as

$$\mathrm{\Pi }^{\mu \nu }(p_0,\overline{p};T)=\mathrm{\Pi }^{(\mathrm{vac})\mu \nu }(p_0,\overline{p})+\overline{\mathrm{\Pi }}^{\mu \nu }(p_0,\overline{p};T).$$
(22)
Accordingly each of four components of the two-point functions defined in Eqs. (17)–(20) is divided into the one part for the quantum correction and another for the hadronic thermal correction. We note that all the divergences are included in the zero temperature part $`\mathrm{\Pi }^{(\mathrm{vac})}`$.

Now, let us present the entire list of the hadronic thermal corrections to the two-point functions. (For the explicit forms of the quantum corrections $`\mathrm{\Pi }^{(\mathrm{vac})}`$, see Appendix C.)

We show the diagrams for contributions to $`\mathrm{\Pi }_{}^{\mu \nu }`$ at one-loop level in Fig. 1.

The temperature dependent parts are obtained as

$`\overline{\mathrm{\Pi }}_{}^t(p_0,\overline{p};T)`$ $`=`$ $`N_faM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,0;T)+N_f{\displaystyle \frac{a}{4}}\overline{B}^t(p_0,\overline{p};M_\rho ,0;T)`$ (23)
$`+N_f(a1)\overline{A}_0(0;T),`$
$`\overline{\mathrm{\Pi }}_{}^s(p_0,\overline{p};T)`$ $`=`$ $`N_faM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,0;T)+N_f{\displaystyle \frac{a}{4}}\overline{B}^s(p_0,\overline{p};M_\rho ,0;T)`$ (24)
$`+N_f(a1)\overline{A}_0(0;T),`$
$`\overline{\mathrm{\Pi }}_{}^L(p_0,\overline{p};T)`$ $`=`$ $`N_f{\displaystyle \frac{a}{4}}\overline{B}^L(p_0,\overline{p};M_\rho ,0;T),`$ (25)
$`\overline{\mathrm{\Pi }}_{}^T(p_0,\overline{p};T)`$ $`=`$ $`N_f{\displaystyle \frac{a}{4}}\overline{B}^T(p_0,\overline{p};M_\rho ,0;T).`$ (26)
As we will see later in section 5, we define the vector meson mass as the pole of longitudinal or transverse component of the vector meson propagator. The diagrams for contributions to $`\mathrm{\Pi }_V^{\mu \nu }`$ are shown in Fig. 2.

We obtain the hadronic thermal corrections to $`\overline{\mathrm{\Pi }}_V^{\mu \nu }`$ as

$`\overline{\mathrm{\Pi }}_V^t(p_0,\overline{p};T)`$ $`=`$ $`2N_f\overline{A}_0(M_\rho ;T)M_\rho ^2N_f\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (27)
$`+{\displaystyle \frac{a^2}{8}}N_f\overline{B}^t(p_0,\overline{p};0,0;T)+{\displaystyle \frac{9}{8}}N_f\overline{B}^t(p_0,\overline{p};M_\rho ,M_\rho ;T),`$
$`\overline{\mathrm{\Pi }}_V^s(p_0,\overline{p};T)`$ $`=`$ $`2N_f\overline{A}_0(M_\rho ;T)M_\rho ^2N_f\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (28)
$`+{\displaystyle \frac{a^2}{8}}N_f\overline{B}^s(p_0,\overline{p};0,0;T)+{\displaystyle \frac{9}{8}}N_f\overline{B}^s(p_0,\overline{p};M_\rho ,M_\rho ;T),`$
$`\overline{\mathrm{\Pi }}_V^L(p_0,\overline{p};T)`$ $`=`$ $`4p^2N_f\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (29)
$`+{\displaystyle \frac{a^2}{8}}N_f\overline{B}^L(p_0,\overline{p};0,0;T)+{\displaystyle \frac{9}{8}}N_f\overline{B}^L(p_0,\overline{p};M_\rho ,M_\rho ;T),`$
$`\overline{\mathrm{\Pi }}_V^T(p_0,\overline{p};T)`$ $`=`$ $`4p^2N_f\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (30)
$`+{\displaystyle \frac{a^2}{8}}N_f\overline{B}^T(p_0,\overline{p};0,0;T)+{\displaystyle \frac{9}{8}}N_f\overline{B}^T(p_0,\overline{p};M_\rho ,M_\rho ;T).`$
Another two-point function associated with the vector current correlator is $`\mathrm{\Pi }_V^{\mu \nu }`$. We also show the one-loop diagrams for contributions to $`\mathrm{\Pi }_V^{\mu \nu }`$ in Fig. 3.

We get the temperature dependent parts as

$`\overline{\mathrm{\Pi }}_V^t(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{a}{2}}N_f\overline{A}_0(0;T)+{\displaystyle \frac{1}{2}}N_f\overline{A}_0(M_\rho ;T)+N_fM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (31)
$`+{\displaystyle \frac{1}{8}}N_f\overline{B}^t(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{a(2a)}{8}}N_f\overline{B}^t(p_0,\overline{p};0,0;T),`$
$`\overline{\mathrm{\Pi }}_V^s(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{a}{2}}N_f\overline{A}_0(0;T)+{\displaystyle \frac{1}{2}}N_f\overline{A}_0(M_\rho ;T)+N_fM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (32)
$`+{\displaystyle \frac{1}{8}}N_f\overline{B}^s(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{a(2a)}{8}}N_f\overline{B}^s(p_0,\overline{p};0,0;T),`$
$`\overline{\mathrm{\Pi }}_V^L(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{1}{8}}N_f\overline{B}^L(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{a(2a)}{8}}N_f\overline{B}^L(p_0,\overline{p};0,0;T),`$ (33)
$`\overline{\mathrm{\Pi }}_V^T(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{1}{8}}N_f\overline{B}^T(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{a(2a)}{8}}N_f\overline{B}^T(p_0,\overline{p};0,0;T).`$ (34)
In section 7, we will define the parameter $`a`$ from the direct $`\gamma \pi \pi `$ coupling using the two-point function $`\mathrm{\Pi }_{}^{\mu \nu }`$. We show the diagrams for contributions to $`\mathrm{\Pi }_{}^{\mu \nu }`$ at one-loop level in Fig. 4.

We get the temperature dependent parts as

$`\overline{\mathrm{\Pi }}_{}^t(p_0,\overline{p};T)`$ $`=`$ $`(a1)N_f\overline{A}_0(0;T)N_fM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (35)
$`+{\displaystyle \frac{1}{8}}N_f\overline{B}^t(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{(2a)^2}{8}}N_f\overline{B}^t(p_0,\overline{p};0,0;T),`$
$`\overline{\mathrm{\Pi }}_{}^s(p_0,\overline{p};T)`$ $`=`$ $`(a1)N_f\overline{A}_0(0;T)N_fM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (36)
$`+{\displaystyle \frac{1}{8}}N_f\overline{B}^s(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{(2a)^2}{8}}N_f\overline{B}^s(p_0,\overline{p};0,0;T),`$
$`\overline{\mathrm{\Pi }}_{}^L(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{1}{8}}N_f\overline{B}^L(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{(2a)^2}{8}}N_f\overline{B}^L(p_0,\overline{p};0,0;T),`$ (37)
$`\overline{\mathrm{\Pi }}_{}^T(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{1}{8}}N_f\overline{B}^T(p_0,\overline{p};M_\rho ,M_\rho ;T)+{\displaystyle \frac{(2a)^2}{8}}N_f\overline{B}^T(p_0,\overline{p};0,0;T).`$ (38)
At the end of this section, using the relations shown in Eq. (B.18), we obtain

$`\overline{\mathrm{\Pi }}_V^t=\overline{\mathrm{\Pi }}_V^t=\overline{\mathrm{\Pi }}_{}^t`$
$`=\overline{\mathrm{\Pi }}_V^s=\overline{\mathrm{\Pi }}_V^s=\overline{\mathrm{\Pi }}_{}^s`$
$`=N_f{\displaystyle \frac{1}{4}}\left[\overline{A}_0(M_\rho ;T)+\overline{A}_0(0;T)\right]N_fM_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T).`$ (39)
Since the quantum corrections to the corresponding components are identical as we have shown in Eq. (C.43), the above relation implies that the components $`\mathrm{\Pi }_V^t`$, $`\mathrm{\Pi }_V^t`$, $`\mathrm{\Pi }_{}^t`$, $`\mathrm{\Pi }_V^s`$, $`\mathrm{\Pi }_V^s`$ and $`\mathrm{\Pi }_{}^t`$ agree:

$$\mathrm{\Pi }_V^t=\mathrm{\Pi }_V^t=\mathrm{\Pi }_{}^t=\mathrm{\Pi }_V^s=\mathrm{\Pi }_V^s=\mathrm{\Pi }_{}^s.$$
(40)
This relation is important to prove the conservation of the vector current correlator as shown in Ref. HKRS <sup>#1</sup><sup>#1</sup>#1 Actually, for the conservation of the vector current correlator, $`\mathrm{\Pi }_V^t=\mathrm{\Pi }_V^t=\mathrm{\Pi }_{}^t`$ and $`\mathrm{\Pi }_V^s=\mathrm{\Pi }_V^s=\mathrm{\Pi }_{}^s`$ are enough, and $`\mathrm{\Pi }_V^t`$ and $`\mathrm{\Pi }_V^s`$ can be generally different.

## 4 Intrinsic Thermal Effects

In this section, we first give an account of the general idea of the intrinsic temperature and/or density dependences of the bare parameters. Next, following Ref. HKRS , we briefly review how to extend the Wilsonian matching to the version at non-zero temperature in order to incorporate the intrinsic thermal effect into the bare parameters of the HLS Lagrangian. There we discuss the effect of Lorentz symmetry violation at bare level, and summarize the conditions for the bare parameters obtained in Ref. HSasaki  through the Wilsonian matching at the critical temperature.

### 4.1 General Concept of the Intrinsic Thermal and/or Density Effects

In this subsection, we give a general idea of the matching between the effective field theory (EFT) and QCD. Green’s functions have important information of system, which is generated by the following functional of a set of source fields denoted by $`J`$:

$$Z[J]=𝒟q𝒟\overline{q}𝒟Ge^{S_{\mathrm{QCD}}[J]},$$
(41)
where $`q`$ $`(\overline{q})`$ denotes quark or antiquark field, $`G`$ is gluon field and $`S_{\mathrm{QCD}}`$ represents the action expressed in terms of the quarks and gluons. The basic concept of the EFT is that the effective Lagrangian, which has the most general form constructed from the chiral symmetry, can give the same generating functional as in Eq. (41):

$$Z[J]=𝒟Ue^{S_{\mathrm{eff}}[J]},$$
(42)
where $`U`$ denotes the relevant hadronic fields such as the pion fields and $`S_{\mathrm{eff}}`$ is the action expressed in terms of these hadrons. Corresponding to each Green’s function derived from Eq. (41), we have the same Green’s function obtained from the EFT through Eq. (42). Performing the matching between these two Green’s functions, we can determine the parameters of the EFT.

In some matching schemes, the renormalized parameters of the EFT are determined from QCD. On the other hand, the matching in the Wilsonian sense is performed based on the following general idea: The bare Lagrangian of the EFT is defined at a suitable matching scale $`\mathrm{\Lambda }`$ and the generating functional derived from the bare Lagrangian leads to the same Green’s function as that derived from Eq. (41) at $`\mathrm{\Lambda }`$. Then the bare parameters of the EFT are determined through the Wilsonian matching. In other words, we obtain the bare Lagrangian of the EFT after integrating out the high energy modes, i.e., the quarks and gluons above $`\mathrm{\Lambda }`$. The information of the high energy modes is included in the parameters of the EFT. Thus when we integrate out high energy modes in hot and/or dense matter, the parameters are in general dependent on temperature and/or density. The intrinsic temperature and/or density dependences are nothing but the signature that the hadron has an internal structure constructed from the quarks and gluons. This is similar to the situation where the coupling constants among hadrons are replaced with the momentum-dependent form factor in the high energy region. Thus the intrinsic thermal and/or dense effects play more important roles in the higher temperature region, especially near the critical temperature.

### 4.2 Wilsonian Matching

The Wilsonian matching proposed in Ref. HYa  is done by matching the axial vector and vector current correlators derived from the HLS with those from the operator product expansion (OPE) in QCD at the matching scale $`\mathrm{\Lambda }`$ <sup>#2</sup><sup>#2</sup>#2 For the validity of the expansion in the HLS the matching scale $`\mathrm{\Lambda }`$ must be smaller than the chiral symmetry breaking scale $`\mathrm{\Lambda }_\chi `$. .

Here it should be noticed that there is no longer Lorentz symmetry in hot matter, and the Lorentz non-scalar operators such as $`\overline{q}\gamma _\mu D_\nu q`$ may exist in the form of the current correlators derived from the OPE HKL . This leads to a difference between the temporal and spatial bare pion decay constants. However, we neglect the contributions from these operators since they give a small correction compared with the main term $`1+\frac{\alpha _s}{\pi }`$ HSasaki . This implies that the Lorentz symmetry breaking effect in the bare pion decay constant is small, $`F_{\pi ,\mathrm{bare}}^tF_{\pi ,\mathrm{bare}}^s`$ HKRS . Thus to a good approximation we determine the pion decay constant at non-zero temperature through the matching condition at zero temperature, putting possible temperature dependences into the gluonic and quark condensates HSasaki ; HKRS :

$$\frac{F_\pi ^2(\mathrm{\Lambda };T)}{\mathrm{\Lambda }^2}=\frac{1}{8\pi ^2}\left[1+\frac{\alpha _s}{\pi }+\frac{2\pi ^2}{3}\frac{\frac{\alpha _s}{\pi }G_{\mu \nu }G^{\mu \nu }_T}{\mathrm{\Lambda }^4}+\pi ^3\frac{1408}{27}\frac{\alpha _s\overline{q}q_T^2}{\mathrm{\Lambda }^6}\right].$$
(43)
Through this condition the temperature dependences of the quark and gluonic condensates determine the intrinsic temperature dependences of the bare parameter $`F_\pi (\mathrm{\Lambda };T)`$, which is then converted into those of the on-shell parameter $`F_\pi (\mu =0;T)`$ through the Wilsonian RGEs.

Now, let us consider the Wilsonian matching near the chiral symmetry restoration point. Here we assume that the order of the chiral phase transition is second or weakly first order, and thus the quark condensate becomes zero continuously for $`TT_c`$. First, note that the Wilsonian matching condition (43) provides

$$\frac{F_\pi ^2(\mathrm{\Lambda };T_c)}{\mathrm{\Lambda }^2}=\frac{1}{8\pi ^2}\left[1+\frac{\alpha _s}{\pi }+\frac{2\pi ^2}{3}\frac{\frac{\alpha _s}{\pi }G_{\mu \nu }G^{\mu \nu }_{T_c}}{\mathrm{\Lambda }^4}\right]0,$$
(44)
which implies that the matching with QCD dictates

$$F_\pi ^2(\mathrm{\Lambda };T_c)0$$
(45)
even at the critical temperature where the on-shell pion decay constant vanishes by adding the quantum corrections through the RGE including the quadratic divergence HYa  and hadronic thermal corrections, as we will show in section 6. As was shown in Ref. HKR  for the VM in dense matter, the Lorentz non-invariant version of the VM conditions for the bare parameters are obtained by the requirement of the equality between the axial vector and vector current correlators in the HLS, which should be valid also in hot matter HKRS :

$`a_{\mathrm{bare}}^t\left({\displaystyle \frac{F_{\sigma ,\mathrm{bare}}^t}{F_{\pi ,\mathrm{bare}}^t}}\right)^2\stackrel{TT_c}{}1,a_{\mathrm{bare}}^s\left({\displaystyle \frac{F_{\sigma ,\mathrm{bare}}^s}{F_{\pi ,\mathrm{bare}}^s}}\right)^2\stackrel{TT_c}{}1,`$ (46)
$`g_{T,\mathrm{bare}}\stackrel{TT_c}{}0,g_{L,\mathrm{bare}}\stackrel{TT_c}{}0,`$ (47)
where $`a_{\mathrm{bare}}^t,a_{\mathrm{bare}}^s,g_{T,\mathrm{bare}}`$ and $`g_{L,\mathrm{bare}}`$ are the extensions of the parameters $`a_{\mathrm{bare}}`$ and $`g_{\mathrm{bare}}`$ in the bare Lagrangian with the Lorentz symmetry breaking effect included as in Appendix A of Ref. HKR .

When we use the conditions for the parameters $`a^{t,s}`$ in Eq. (46) and the above result that the Lorentz symmetry violation between the bare pion decay constants $`F_{\pi ,\mathrm{bare}}^{t,s}`$ is small, we can easily show that the Lorentz symmetry breaking effect between the temporal and spatial bare sigma decay constants is also small, $`F_{\sigma ,\mathrm{bare}}^tF_{\sigma ,\mathrm{bare}}^s`$ HKRS . While we cannot determine the ratio $`g_{L,\mathrm{bare}}/g_{T,\mathrm{bare}}`$ through the Wilsonian matching since the transverse mode of vector meson decouples near the critical temperature. <sup>#3</sup><sup>#3</sup>#3 In Ref. HKR , the analysis including the Lorentz non-invariance at bare HLS theory was carried out. Due to the equality between axial vector and vector current correlators, $`(g_{L,\mathrm{bare}},g_{T,\mathrm{bare}})(0,0)`$ is satisfied when we approach the critical point. This implies that at the bare level the longitudinal mode becomes the real NG boson and couples to the vector current correlator, while the transverse mode decouples. Furthermore $`g_L0`$ is a fixed point for the RGE Sasaki:2003qj . Thus in any energy scale the transverse mode decouples from the vector current correlator.  However this implies that the transverse mode is irrelevant for the quantities studied in this paper. Therefore in the present analysis, we set $`g_{L,\mathrm{bare}}=g_{T,\mathrm{bare}}`$ for simplicity and use the Lorentz invariant Lagrangian at bare level. In the low temperature region, the intrinsic temperature dependences are negligible, so that we also use the Lorentz invariant Lagrangian at bare level as in the analysis by the ordinary chiral Lagrangian in Ref. GL .

At the critical temperature, the axial vector and vector current correlators derived in the OPE agree with each other for any value of $`Q^2`$. Thus we require that these current correlators in the HLS are equal at the critical temperature for any value of $`Q^2\text{around}\mathrm{\Lambda }^2`$. As we discussed above, we start from the Lorentz invariant bare Lagrangian even in hot matter, and then the axial vector current correlator $`G_A^{(\mathrm{HLS})}`$ and the vector current correlator $`G_V^{(\mathrm{HLS})}`$ are expressed by the same forms as those at zero temperature with the bare parameters having the intrinsic temperature dependences:

$`G_A^{(\mathrm{HLS})}(Q^2;T)`$ $`=`$ $`{\displaystyle \frac{F_\pi ^2(\mathrm{\Lambda };T)}{Q^2}}2z_2(\mathrm{\Lambda };T),`$
$`G_V^{(\mathrm{HLS})}(Q^2;T)`$ $`=`$ $`{\displaystyle \frac{F_\sigma ^2(\mathrm{\Lambda };T)[12g^2(\mathrm{\Lambda };T)z_3(\mathrm{\Lambda };T)]}{M_{\rho }^{}{}_{}{}^{2}(\mathrm{\Lambda };T)+Q^2}}2z_1(\mathrm{\Lambda };T).`$ (48)
By taking account of the fact $`F_\pi ^2(\mathrm{\Lambda };T_c)0`$ derived from the Wilsonian matching condition given in Eq. (44), the requirement $`G_A^{(\mathrm{HLS})}=G_V^{(\mathrm{HLS})}`$ is satisfied only if the following conditions are met HSasaki :

$`g(\mathrm{\Lambda };T)`$ $`\stackrel{TT_c}{}`$ $`0,`$ (49)
$`a(\mathrm{\Lambda };T)=F_\sigma ^2(\mathrm{\Lambda };T)/F_\pi ^2(\mathrm{\Lambda };T)`$ $`\stackrel{TT_c}{}`$ $`1,`$ (50)
$`z_1(\mathrm{\Lambda };T)z_2(\mathrm{\Lambda };T)`$ $`\stackrel{TT_c}{}`$ $`0.`$ (51)
Note that the intrinsic thermal effects act on the parameters in such a way that they become the values of Eqs. (49)-(51).

Through the Wilsonian matching at non-zero temperature mentioned above, the parameters appearing in the hadronic thermal corrections calculated in section 3 have the intrinsic temperature dependences: $`F_\pi `$, $`a`$ and $`g`$ appearing there should be regarded as

$`F_\pi `$ $``$ $`F_\pi (\mu =0;T),`$
$`a`$ $``$ $`a\left(\mu =M_\rho (T);T\right),`$
$`g`$ $``$ $`g\left(\mu =M_\rho (T);T\right),`$ (52)
where $`M_\rho `$ is determined from the on-shell condition:

$$M_\rho ^2M_\rho ^2(T)=a(\mu =M_\rho ;T)g^2(\mu =M_\rho ;T)F_\pi ^2(\mu =M_\rho ;T).$$
(53)
From the RGEs for $`g\text{and}a`$ in Eqs.(C.58) and (C.57), we find that $`(g,a)=(0,1)`$ is the fixed point. Therefore, Eqs. (49) and (50) imply that $`g`$ and $`a`$ at the on-shell of the vector meson take the same values:

$`a(\mu =M_\rho (T);T)\stackrel{TT_c}{}1,`$
$`g(\mu =M_\rho (T);T)\stackrel{TT_c}{}0,`$ (54)
where the parametric vector meson mass $`M_\rho (T)`$ is determined from the condition (53). The above conditions with Eq. (53) imply that $`M_\rho (T)`$ also vanishes:

$$M_\rho (T)\stackrel{TT_c}{}0.$$
(55)
## 5 Vector Meson Mass

In Ref. HSasaki  we have shown that the vector manifestation (VM) in hot matter can actually be formulated by using the hadronic thermal correction to the vector meson pole mass calculated in the Landau gauge HS . In Ref. HKRS  the background field gauge was used to calculate the hadronic thermal corrections. There the vector meson pole mass was not explicitly studied. In this section we first define the vector meson pole masses of both the longitudinal and transverse modes at non-zero temperature from the vector current correlator in the background field gauge and show the explicit forms of the hadronic thermal corrections from the vector and pseudoscalar meson loop. Then, including the intrinsic temperature dependences of the parameters near the critical temperature determined in the previous section, we show that the vector meson mass vanishes at the critical temperature.

Let us define pole masses of longitudinal and transverse modes of the vector meson from the poles of longitudinal and transverse components of the vector current correlator in the rest frame HKRS :

$$G_V^{\mu \nu }=P_L^{\mu \nu }G_V^L+P_T^{\mu \nu }G_V^T,$$
(56)
where

$`G_V^L`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Pi }_V^s\left(\mathrm{\Pi }_V^L+2\mathrm{\Pi }_V^L\right)}{\mathrm{\Pi }_V^s\mathrm{\Pi }_V^L}}+\mathrm{\Pi }_{}^L,`$
$`G_V^T`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Pi }_V^s\left(\mathrm{\Pi }_V^T+2\mathrm{\Pi }_V^T\right)}{\mathrm{\Pi }_V^s\mathrm{\Pi }_V^T}}+\mathrm{\Pi }_{}^T.`$ (57)
Then, the pole masses are obtained as the solutions of the following on-shell conditions:

$`0=\text{Re}\left[\mathrm{\Pi }_V^s(p_0=m_\rho ^L(T),0;T)\mathrm{\Pi }_V^L(p_0=m_\rho ^L(T),0;T)\right],`$
$`0=\text{Re}\left[\mathrm{\Pi }_V^s(p_0=m_\rho ^T(T),0;T)\mathrm{\Pi }_V^T(p_0=m_\rho ^T(T),0;T)\right],`$ (58)
where $`m_\rho ^L(T)`$ and $`m_\rho ^T(T)`$ denote the pole masses of the longitudinal and transverse modes, respectively. As we have calculated in section 3, $`\mathrm{\Pi }_V^s(p_0,\overline{p};T)`$, $`\mathrm{\Pi }_V^L(p_0,\overline{p};T)`$ and $`\mathrm{\Pi }_V^T(p_0,\overline{p};T)`$ in the HLS at one-loop level are expressed as

$`\mathrm{\Pi }_V^s(p_0,\overline{p};T)`$ $`=`$ $`F_\sigma ^2(M_\rho )+\stackrel{~}{\mathrm{\Pi }}_V^S(p^2)+\overline{\mathrm{\Pi }}_V^s(p_0,\overline{p};T),`$
$`\mathrm{\Pi }_V^L(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{p^2}{g^2(M_\rho )}}+\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2)+\overline{\mathrm{\Pi }}_V^L(p_0,\overline{p};T),`$
$`\mathrm{\Pi }_V^T(p_0,\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{p^2}{g^2(M_\rho )}}+\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2)+\overline{\mathrm{\Pi }}_V^T(p_0,\overline{p};T),`$ (59)
where the explicit forms of the finite renormalization effects $`\stackrel{~}{\mathrm{\Pi }}_V^S(p^2)`$ and $`\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2)`$ are given in Eqs. (C.60) and (C.61), and those of the hadronic thermal effects $`\overline{\mathrm{\Pi }}_V^s(p_0,\overline{p};T)`$, $`\overline{\mathrm{\Pi }}_V^L(p_0,\overline{p};T)`$ and $`\overline{\mathrm{\Pi }}_V^T(p_0,\overline{p};T)`$ are given in Eqs. (39), (29) and (30). Substituting Eq. (59) into Eq. (58), we obtain

$`\left[m_\rho ^L(T)\right]^2`$ $`=`$ $`M_\rho ^2+g^2(M_\rho )[\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^S(p_0^2)+\text{Re}\overline{\mathrm{\Pi }}_V^s(p_0,0;T)`$
$`\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^T(p_0^2)\text{Re}\overline{\mathrm{\Pi }}_V^L(p_0,0;T)]_{p_0=m_\rho ^L(T)},`$
$`\left[m_\rho ^T(T)\right]^2`$ $`=`$ $`M_\rho ^2+g^2(M_\rho )[\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^S(p_0^2)+\text{Re}\overline{\mathrm{\Pi }}_V^s(p_0,0;T)`$ (60)
$`\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^T(p_0^2)\text{Re}\overline{\mathrm{\Pi }}_V^T(p_0,0;T)]_{p_0=m_\rho ^T(T)}.`$
We can replace $`m_\rho ^L(T)`$ and $`m_\rho ^T(T)`$ with $`M_\rho `$ in the hadronic thermal effects as well as in the finite renormalization effect, since the difference is of higher order. Then, noting that $`\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^S(p^2=M_\rho ^2)=\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2=M_\rho ^2)=0`$ as shown in Eq. (C.48), we obtain

$`\left[m_\rho ^L(T)\right]^2`$ $`=`$ $`M_\rho ^2g^2(M_\rho )\text{Re}\left[\overline{\mathrm{\Pi }}_V^L(M_\rho ;0;T)\overline{\mathrm{\Pi }}_V^s(M_\rho ;0;T)\right],`$ (61)
$`\left[m_\rho ^T(T)\right]^2`$ $`=`$ $`M_\rho ^2g^2(M_\rho )\text{Re}\left[\overline{\mathrm{\Pi }}_V^T(M_\rho ;0;T)\overline{\mathrm{\Pi }}_V^s(M_\rho ;0;T)\right].`$ (62)
Let us compute the explicit form of the pole mass using the expression of $`\overline{\mathrm{\Pi }}_V^s`$, $`\overline{\mathrm{\Pi }}_V^L`$ and $`\overline{\mathrm{\Pi }}_V^T`$ calculated in section 3. Here we note that Eqs. (B.28) and (B.29) imply that $`\overline{B}^L\overline{B}^s`$ agrees with $`\overline{B}^T\overline{B}^s`$ in the rest frame. Then in the rest frame $`\overline{\mathrm{\Pi }}_V^L\overline{\mathrm{\Pi }}_V^s`$ agrees with $`\overline{\mathrm{\Pi }}_V^T\overline{\mathrm{\Pi }}_V^s`$. Thus the longitudinal pole mass is the same as the transverse one:

$$m_\rho ^L(T)=m_\rho ^T(T)m_\rho (T).$$
(63)
By using the low momentum limits of the functions shown in Eqs. (B.27) and (B.28), $`\overline{\mathrm{\Pi }}_V^L\overline{\mathrm{\Pi }}_V^s`$ in the rest frame is expressed as

$`\overline{\mathrm{\Pi }}_V^L(p_0,\overline{p}=0;T)\overline{\mathrm{\Pi }}_V^s(p_0,\overline{p}=0;T)`$ (64)
$`=`$ $`N_f\left[{\displaystyle \frac{a^2}{24}}\stackrel{~}{G}_2(p_0;T)\stackrel{~}{J}_1^2(M_\rho ;T)+\left({\displaystyle \frac{M_\rho ^2}{4}}p_0^2\right)\stackrel{~}{F}_3^2(p_0;M_\rho ;T)+{\displaystyle \frac{3}{8}}\stackrel{~}{F}_3^4(p_0;M_\rho ;T)\right],`$
where the functions $`\stackrel{~}{I}_n`$, $`\stackrel{~}{J}_m^n`$, $`\stackrel{~}{F}_m^n`$ and $`\stackrel{~}{G}_m^n`$ are defined in Appendix D. Substituting the above expression into Eq. (61) and using the relation

$`{\displaystyle \frac{1}{3M_\rho ^2}}\stackrel{~}{F}_3^4(M_\rho ;M_\rho ;T)`$ $`=`$ $`{\displaystyle \frac{1}{4}}\stackrel{~}{F}_3^2(M_\rho ;M_\rho ;T){\displaystyle \frac{1}{3M_\rho ^2}}\stackrel{~}{J}_1^2(M_\rho ;T),`$ (65)
we find that the vector meson pole mass is expressed as

$`m_\rho ^2(T)`$ $`=`$ $`M_{\rho }^{}{}_{}{}^{2}+N_fg^2\left[{\displaystyle \frac{a^2}{12}}\stackrel{~}{G}_2(M_\rho ;T)+{\displaystyle \frac{4}{5}}\stackrel{~}{J}_1^2(M_\rho ;T)+{\displaystyle \frac{33}{16}}M_\rho ^2\stackrel{~}{F}_3^2(M_\rho ;M_\rho ;T)\right].`$ (66)
The contribution in this expression agrees with the result in Ref. HSasaki  which is derived from the hadronic thermal correction calculated in the Landau gauge in Ref. HS  by taking the on-shell condition (61). <sup>#4</sup><sup>#4</sup>#4The functions used in this paper are related to the ones used in Ref. HSasaki  as $`\stackrel{~}{J}_1^2=\frac{1}{2\pi ^2}J_1^2`$, $`\stackrel{~}{G}_2=\frac{1}{2\pi ^2}\overline{G}_2`$, and so on.

Before going to the analysis near the critical temperature, we study the behavior of the pole mass in the low temperature region, $`TM_\rho `$. In this region the functions $`F_m^n`$ and $`J_m^n`$ are suppressed by $`e^{M_\rho /T}`$, and thus give negligible contributions. Since $`G_2\frac{\pi ^4}{15}\frac{T^4}{M_{\rho }^{}{}_{}{}^{2}}`$, the vector meson pole mass increases as $`T^4`$ in the low temperature region, dominated by the $`\pi `$-loop effects:

$$m_{\rho }^{}{}_{}{}^{2}(T)M_{\rho }^{}{}_{}{}^{2}+\frac{N_f\pi ^2a}{360F_{\pi }^{}{}_{}{}^{2}}T^4\text{for}TM_\rho .$$
(67)
Note that the lack of $`T^2`$-term in the above expression is consistent with the result by the current algebra analysis Dey-Eletsky-Ioffe .

Now, let us study the vector meson pole mass near the critical temperature. As shown in section 4, the intrinsic temperature dependences of the parameters of the HLS Lagrangian determined through the Wilsonian matching imply that the parametric vector meson mass $`M_\rho `$ vanishes at the critical temperature. Then, near the critical temperature we should take $`M_\rho T`$ in Eq. (66) instead of $`TM_\rho `$ which was taken to reach the expression in Eq. (67) in the low temperature region. Thus, by noting that

$`\stackrel{~}{G}_2(M_\rho ;T)`$ $`\stackrel{M_\rho 0}{}`$ $`\stackrel{~}{I}_2(T),`$
$`\stackrel{~}{J}_1^2(M_\rho ;T)`$ $`\stackrel{M_\rho 0}{}`$ $`\stackrel{~}{I}_2(T),`$
$`M_{\rho }^{}{}_{}{}^{2}\stackrel{~}{F}_3^2(M_\rho ;M_\rho ;T)`$ $`\stackrel{M_\rho 0}{}`$ $`0,`$ (68)
the pole mass of the vector meson at $`TT_c`$ becomes

$`m_{\rho }^{}{}_{}{}^{2}(T)`$ $`=`$ $`M_{\rho }^{}{}_{}{}^{2}+N_fg^2{\displaystyle \frac{15a^2}{12}}\stackrel{~}{I}_2(T)`$ (69)
$`=`$ $`M_{\rho }^{}{}_{}{}^{2}+N_fg^2{\displaystyle \frac{15a^2}{144}}T^2.`$
In the vicinity of $`a1`$ the hadronic thermal effect gives a positive correction, and then the vector meson pole mass is actually larger than the parametric mass $`M_\rho `$. However, the intrinsic temperature dependences of the parameters obtained in section 4 lead to $`g0`$ and $`M_\rho 0`$ for $`TT_c`$. Then, from Eq. (69) we conclude that the pole mass of the vector meson $`m_\rho `$ also vanishes at the critical temperature:

$$m_\rho (T)0\text{for}TT_c.$$
(70)
This implies that the VM is realized at the critical temperature, which is consistent with the picture shown in Refs. BR ; Brown-Rho:96 ; Brown:2001jh ; Brown:2001nh .

## 6  Pion Parameters and the Critical Temperature

In this section we first study the temperature dependences of the temporal and spatial pion decay constants and pion velocity in the low temperature region, and briefly review how the two decay constants vanish at the critical temperature following Ref. HKRS . We also estimate the value of the critical temperature in a wider range of input parameters than the one used in Ref. HSasaki .

### 6.1 Pion Decay Constant and Pion Velocity

Since the Lorentz symmetry is violated at non-zero temperature, there exist two pion decay constants $`f_\pi ^t`$ and $`f_\pi ^s`$ PT:96 . These $`f_\pi ^t`$ and $`f_\pi ^s`$ are associated with the temporal and spatial components of the axial vector current, respectively. In the present analysis, they are expressed as HKRS

$`f_\pi ^t(\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Pi }_{}^t(E,\overline{p};T)}{\stackrel{~}{F}(\overline{p};T)}}={\displaystyle \frac{F_\pi ^2+\overline{\mathrm{\Pi }}_{}^t(\overline{p},\overline{p};T)}{\stackrel{~}{F}(\overline{p};T)}},`$
$`f_\pi ^s(\overline{p};T)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Pi }_{}^s(E,\overline{p};T)}{\stackrel{~}{F}(\overline{p};T)}}={\displaystyle \frac{F_\pi ^2+\overline{\mathrm{\Pi }}_{}^s(\overline{p},\overline{p};T)}{\stackrel{~}{F}(\overline{p};T)}},`$ (71)
where $`E`$ denotes the pion energy and we replaced $`E`$ by $`\overline{p}`$ in the hadronic thermal corrections $`\overline{\mathrm{\Pi }}_{}^t(E,\overline{p};T)`$ and $`\overline{\mathrm{\Pi }}_{}^s(E,\overline{p};T)`$ since the difference is of higher order. $`\stackrel{~}{F}(\overline{p};T)`$ is the wave function renormalization constant of the $`\pi `$ field given by MOW ; HKRS

$`\stackrel{~}{F}^2(\overline{p};T)`$ $`=`$ $`\text{Re}\mathrm{\Pi }_{}^t(E,\overline{p};T)`$ (72)
$`=`$ $`F_\pi ^2+\text{Re}\overline{\mathrm{\Pi }}_{}^t(\overline{p},\overline{p};T).`$
By using the above decay constants, the pion velocity at one-loop level is expressed as PT:96 ; HKRS  <sup>#5</sup><sup>#5</sup>#5 This form of the pion velocity looks slightly different from the one used in Ref. HKRS . However, this is actually equivalent to the one in Ref. HKRS  at one-loop order, and more convenient to study the temperature dependence of the pion velocity in the low temperature region.

$$v_\pi ^2(\overline{p};T)=1+\frac{\stackrel{~}{F}(\overline{p};T)}{F_\pi ^2}\left[\text{Re}f_\pi ^s(\overline{p};T)\text{Re}f_\pi ^t(\overline{p};T)\right].$$
(73)
In Ref. HSasaki , following Ref. Bochkarev:1995gi , we defined the pion decay constant $`f_\pi `$ as the pole residue of the axial vector current correlator. It is related to the above $`f_\pi ^t`$ and $`f_\pi ^s`$ at one-loop order as

$$f_\pi ^2(\overline{p};T)=f_\pi ^t(\overline{p};T)f_\pi ^s(\overline{p};T).$$
(74)
Furthermore we have

$$\stackrel{~}{F}(\overline{p};T)=\text{Re}f_\pi ^t(\overline{p};T).$$
(75)
Substituting Eqs. (B.5) and (B.6) into Eqs. (23) and (24), we obtain $`\overline{\mathrm{\Pi }}_{}^t`$ and $`\overline{\mathrm{\Pi }}_{}^s`$ for the on-shell pion as

$`\overline{\mathrm{\Pi }}_{}^t(p_0=\overline{p}+iϵ,\overline{p};T)=N_f(a1)\stackrel{~}{I}_2(T)`$
$`+{\displaystyle \frac{N_f}{2}}a\left[{\displaystyle \frac{1}{2}}\overline{B}^t(\overline{p}+iϵ,\overline{p};M_\rho ,0;T)2M_\rho ^2\overline{B}_0(\overline{p}+iϵ,\overline{p};M_\rho ,0;T)\right],`$ (76)
$`\overline{\mathrm{\Pi }}_{}^s(p_0=\overline{p}+iϵ,\overline{p};T)=N_f(a1)\stackrel{~}{I}_2(T)`$
$`+{\displaystyle \frac{N_f}{2}}a\left[{\displaystyle \frac{1}{2}}\overline{B}^s(\overline{p}+iϵ,\overline{p};M_\rho ,0;T)2M_\rho ^2\overline{B}_0(\overline{p}+iϵ,\overline{p};M_\rho ,0;T)\right],`$ (77)
where we put $`ϵ+0`$ to make the analytic continuation for the frequency $`p_0=i2\pi nT`$ to the Minkowski variable. We show the explicit forms of the functions $`\overline{B}^t`$, $`\overline{B}^s`$ and $`\overline{B}_0`$ in Eqs. (LABEL:B0\_A\_on-shell), (LABEL:Bt\_A\_on-shell) and (LABEL:Bs\_A\_on-shell).

In general the pion velocity in medium does not agree with the value at $`T=0`$ due to the interaction with the heat bath. Below $`T_c`$, since the temporal component does not agree with the spatial one due to the thermal vector meson effect, $`\overline{\mathrm{\Pi }}_{}^t\overline{\mathrm{\Pi }}_{}^s`$, the pion velocity $`v_\pi (\overline{p};T)`$ is not the speed of light. As we will see below, in the framework of HLS the pion velocity receives a change from the $`\rho `$-loop effect for $`0<T<T_c`$. When we take the low temperature limit ($`TM_\rho `$) and the soft pion limit ($`\overline{p}M_\rho `$ and $`\overline{p}T`$), the real parts of $`\overline{\mathrm{\Pi }}_{}^t`$ and $`\overline{\mathrm{\Pi }}_{}^s`$ are approximated as

$`\text{Re}\overline{\mathrm{\Pi }}_{}^t(p_0=\overline{p}+iϵ,\overline{p};T)`$ $``$ $`N_f\stackrel{~}{I}_2(T)+N_f{\displaystyle \frac{a}{M_\rho ^2}}\stackrel{~}{I}_4(T)N_fa\sqrt{{\displaystyle \frac{M_\rho }{8\pi ^3}}}e^{M_\rho /T}T^{3/2},`$ (78)
$`\text{Re}\overline{\mathrm{\Pi }}_{}^s(p_0=\overline{p}+iϵ,\overline{p};T)`$ $``$ $`N_f\stackrel{~}{I}_2(T)N_f{\displaystyle \frac{a}{M_\rho ^2}}\stackrel{~}{I}_4(T)+N_fa\sqrt{{\displaystyle \frac{M_\rho }{8\pi ^3}}}e^{M_\rho /T}T^{3/2}.`$ (79)
By using Eqs. (78) and (79) and neglecting the terms proportional to the suppression factor $`e^{M_\rho /T}`$, the pion velocity is expressed as

$`v_\pi ^2(\overline{p};T)`$ $``$ $`1N_f{\displaystyle \frac{2a}{F_\pi ^2M_\rho ^2}}\stackrel{~}{I}_4(T)`$ (80)
$`=`$ $`1{\displaystyle \frac{N_f}{15}}a\pi ^2{\displaystyle \frac{T^4}{F_\pi ^2M_\rho ^2}}<1.`$
This shows that the pion velocity is smaller than the speed of light already at one-loop level in the HLS due to the $`\rho `$-loop effect. This is different from the result obtained in the ordinary ChPT including only the pion at one-loop \[see for example, Ref. PT:96  and references therein\]. Generally, the longitudinal $`\rho `$ contribution to $`\overline{\mathrm{\Pi }}_{}^t`$ expressed by $`\overline{B}^t`$ in Eq. (76) differs from that to $`\overline{\mathrm{\Pi }}_{}^s`$ by $`\overline{B}^s`$ in Eq. (77), which implies that, below the critical temperature, there always exists a deviation of the pion velocity from the speed of light due to the longitudinal $`\rho `$-loop effect:

$$v_\pi ^2(\overline{p};T)<1\text{for}0<T<T_c.$$
(81)
Next, we study the temporal and spatial pion decay constants in hot matter defined by Eq. (71). The imaginary parts of $`\overline{\mathrm{\Pi }}_{}^t`$ and $`\overline{\mathrm{\Pi }}_{}^s`$ in the low temperature region are given by

$`\text{Im}\overline{\mathrm{\Pi }}_{}^t(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ $`\stackrel{\overline{p}T}{}`$ $`{\displaystyle \frac{N_f}{4}}a\text{Im}\overline{B}^t(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ (82)
$`\stackrel{TM_\rho }{}`$ $`{\displaystyle \frac{N_f}{8\pi }}aM_\rho ^2e^{M_\rho /T},`$
$`\text{Im}\overline{\mathrm{\Pi }}_{}^s(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ $`\stackrel{\overline{p}T}{}`$ $`{\displaystyle \frac{N_f}{4}}a\text{Im}\overline{B}^s(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ (83)
$`\stackrel{TM_\rho }{}`$ $`{\displaystyle \frac{N_f}{8\pi }}aM_\rho ^2e^{M_\rho /T}.`$
Thus the contributions from the imaginary parts $`\text{Im}\overline{\mathrm{\Pi }}_{}^{t,s}`$ are small because of the suppression factor $`e^{M_\rho /T}`$. From Eqs. (78) and (79) with $`\stackrel{~}{I}_2(T)=T^2/12`$ and $`\stackrel{~}{I}_4(T)=\pi ^2T^4/30`$, we obtain the following results for the real parts of $`f_\pi ^t`$ and $`f_\pi ^s`$:

$`\left[\text{Re}f_\pi ^t\right]\stackrel{~}{F}`$ $``$ $`F_\pi ^2N_f\left[{\displaystyle \frac{T^2}{12}}{\displaystyle \frac{a\pi ^2}{30M_\rho ^2}}T^4+a\sqrt{{\displaystyle \frac{M_\rho }{8\pi ^3}}}e^{M_\rho /T}T^{3/2}\right],`$
$`\left[\text{Re}f_\pi ^s\right]\stackrel{~}{F}`$ $``$ $`F_\pi ^2N_f\left[{\displaystyle \frac{T^2}{12}}+{\displaystyle \frac{a\pi ^2}{30M_\rho ^2}}T^4a\sqrt{{\displaystyle \frac{M_\rho }{8\pi ^3}}}e^{M_\rho /T}T^{3/2}\right].`$ (84)
We note that $`\stackrel{~}{F}=\text{Re}f_\pi ^t`$ as shown in Eq. (75). Then, neglecting the terms suppressed by $`e^{M_\rho /T}`$, we obtain the difference between $`(f_\pi ^t)^2`$ and $`f_\pi ^tf_\pi ^s`$ as

$$(f_\pi ^t)^2f_\pi ^tf_\pi ^s\frac{N_f}{15}a\pi ^2\frac{T^4}{M_\rho ^2}>0.$$
(85)
This implies that the contribution from the $`\rho `$-loop (the second and third terms in the brackets) generates a difference between the temporal and spatial pion decay constants in the low temperature region. Similarly, at general temperature below $`T_c`$, there exists a difference between $`f_\pi ^t\stackrel{~}{F}`$ and $`f_\pi ^s\stackrel{~}{F}`$ due to the $`\rho `$-loop effects: $`\left[\text{Re}f_\pi ^t\right]\stackrel{~}{F}>\left[\text{Re}f_\pi ^s\right]\stackrel{~}{F}`$. Since we can always take $`\stackrel{~}{F}`$ to be positive, we find that $`\text{Re}f_\pi ^t`$ is larger than $`\text{Re}f_\pi ^s`$:

$$\text{Re}f_\pi ^t(\overline{p};T)>\text{Re}f_\pi ^s(\overline{p};T)\text{for}0<T<T_c.$$
(86)
This result is consistent with Eq. (81) because $`v_\pi ^21=(\stackrel{~}{F}/F_\pi ^2)[\text{Re}f_\pi ^s\text{Re}f_\pi ^t]`$ by definition. The difference between $`(f_\pi ^t)^2`$ and $`f_\pi ^tf_\pi ^s`$ shown in Eq. (85) is tiny, and the hadronic thermal corrections to them are dominated by the first term ($`T^2/12`$) in the bracket in Eq. (84). Then, in the very low temperature region, the above expressions are further reduced to

$$f_\pi ^2=\left(f_\pi ^t\right)\left(f_\pi ^s\right)\left(f_\pi ^t\right)^2F_{\pi }^{}{}_{}{}^{2}\frac{N_f}{12}T^2,$$
(87)
which is consistent with the result obtained in Ref. GL .

Now we study how $`f_\pi ^t`$ and $`f_\pi ^s`$ behave near $`T_c`$. As we have shown in Eq. (44), the Wilsonian matching determines the bare pion decay constant in terms of the parameters appearing in the OPE. Furthermore, as we discussed around Eq. (54), at the critical temperature the intrinsic thermal effects lead to $`(g,a)(0,1)`$ which is a fixed point of the coupled RGEs. Then the RGE for $`F_\pi `$ becomes

$$\mu \frac{dF_{\pi }^{}{}_{}{}^{2}}{d\mu }=\frac{N_f}{(4\pi )^2}\mu ^2.$$
(88)
This RGE is easily solved and the relation between $`F_\pi (\mathrm{\Lambda };T_c)\text{and}F_\pi (0;T_c)`$ is given by

$$F_\pi ^2(0;T_c)=F_\pi ^2(\mathrm{\Lambda };T_c)\frac{N_f}{2(4\pi )^2}\mathrm{\Lambda }^2.$$
(89)
Finally we obtain the pion decay constants as follows:

$`f_\pi ^t\stackrel{~}{F}=F_\pi ^2(0;T)+\overline{\mathrm{\Pi }}_{}^t(\overline{p},\overline{p};T),`$
$`f_\pi ^s\stackrel{~}{F}=F_\pi ^2(0;T)+\overline{\mathrm{\Pi }}_{}^s(\overline{p},\overline{p};T),`$ (90)
where the second terms are the hadronic thermal effects. In the VM limit $`(g,a)(0,1)`$ the temperature dependent parts become

$$\overline{\mathrm{\Pi }}_{}^t(\overline{p},\overline{p};T)\stackrel{TT_c}{}\frac{N_f}{24}T_{c}^{}{}_{}{}^{2},\overline{\mathrm{\Pi }}_{}^s(\overline{p},\overline{p};T)\stackrel{TT_c}{}\frac{N_f}{24}T_{c}^{}{}_{}{}^{2}.$$
(91)
From Eqs. (74), (90) and (91), the order parameter $`f_\pi `$ becomes

$$f_\pi ^2(\overline{p};T)\stackrel{TT_c}{}f_\pi ^2(\overline{p};T_c)=F_\pi ^2(0;T_c)\frac{N_f}{24}T_c^2.$$
(92)
Since the order parameter $`f_\pi `$ vanishes at the critical temperature, this implies

$$F_\pi ^2(0;T)\stackrel{TT_c}{}F_\pi ^2(0;T_c)=\frac{N_f}{24}T_c^2.$$
(93)
Thus we obtain

$$(f_\pi ^t)^2\stackrel{TT_c}{}0,f_\pi ^tf_\pi ^s\stackrel{TT_c}{}0.$$
(94)
From Eq. (86) or equivalently Eq. (81), the above results imply that the temporal and spatial pion decay constants vanish simultaneously at the critical temperature HKRS :

$$f_\pi ^t(T_c)=f_\pi ^s(T_c)=0.$$
(95)
Comparing Eq. (92) with the expression in the low temperature region in Eq. (87) where the vector meson is decoupled, we find that the coefficient of $`T_{c}^{}{}_{}{}^{2}`$ is different by the factor $`\frac{1}{2}`$. This is the contribution from the $`\sigma `$-loop (longitudinal $`\rho `$-loop). In the low temperature region the $`\pi `$-loop effects give the dominant contributions to $`f_\pi (T)`$ and the $`\rho `$-loop effects are negligible. However by the vector manifestation the $`\rho `$-loop contributions are also incorporated into $`f_\pi (T)`$ near the critical temperature.

### 6.2 Critical Temperature

In this subsection we estimate the value of the critical temperature $`T_c`$ where the order parameter $`f_\pi ^2`$ vanishes.

We first determine $`T_c`$ by naively extending the expression (87) to the higher temperature region to get $`T_c^{(\mathrm{hadron})}=180\text{MeV}`$ for $`N_f=3`$. However this naive extension is inconsistent with the chiral restoration in QCD since the axial vector and vector current correlators do not agree with each other at that temperature. As is stressed in Ref. HSasaki , the disagreement between two correlators is cured by including the intrinsic thermal effect. As can be seen from Eq. (43), the intrinsic temperature dependence of the parameter $`F_\pi `$ is determined from $`\frac{\alpha _s}{\pi }G_{\mu \nu }G^{\mu \nu }_T`$ and $`\overline{q}q_T`$, and gives only a small contribution compared with the main term $`1+\frac{\alpha _s}{\pi }`$. However it is important that the parameters in the hadronic corrections have the intrinsic temperature dependences as $`(a,g)(1,0)`$ for $`TT_c`$, which carry the information of QCD. Actually the inclusion of the intrinsic thermal effects provides the formula (92) for the pion decay constant at the critical temperature, in which the second term has an extra factor of $`1/2`$ compared with the one in Eq. (87).

In Ref. HSasaki  we determined the critical temperature from $`f_\pi (T_c)=0`$ and estimated the value which is dependent on the matching scale $`\mathrm{\Lambda }`$: From Eq. (93) we obtain

$$T_c=\sqrt{\frac{24}{N_f}}F_\pi (0;T_c).$$
(96)
Using Eqs.(44) and (89) we get HSasaki

$$\frac{T_c}{\mathrm{\Lambda }}=\sqrt{\frac{3}{N_f\pi ^2}}\left[1+\frac{\alpha _s}{\pi }+\frac{2\pi ^2}{3}\frac{\frac{\alpha _s}{\pi }G_{\mu \nu }G^{\mu \nu }_{T_c}}{\mathrm{\Lambda }^4}\frac{N_f}{4}\right]^{\frac{1}{2}}.$$
(97)
We would like to stress that the critical temperature is expressed in terms of the parameters appearing in the OPE. We evaluate the critical temperature for $`N_f=3`$. The gluonic condensate at $`T_c`$ is about half of the value at $`T=0`$ Miller ; Brown:2001nh  and we use $`\frac{\alpha _s}{\pi }G_{\mu \nu }G^{\mu \nu }_{T_c}=0.006\text{GeV}^4`$. We show the predicted values of $`T_c`$ for several choices of $`\mathrm{\Lambda }_{\mathrm{QCD}}`$ and $`\mathrm{\Lambda }`$ in Table 1.

Note that the Wilsonian matching describes the experimental results very well for $`\mathrm{\Lambda }_{QCD}=0.4\text{GeV}`$ and $`\mathrm{\Lambda }=1.1\text{GeV}`$ at $`T=0`$ HYc . At non-zero temperature, however, the matching scale $`\mathrm{\Lambda }`$ may be dependent on temperature. The smallest $`\mathrm{\Lambda }`$ in Table 1 is determined by requiring $`(2\pi ^2/3)\frac{\alpha _s}{\pi }G_{\mu \nu }G^{\mu \nu }/\mathrm{\Lambda }^4<0.1`$. Here we note that the critical temperature may be changed by including the higher order corrections.

## 7 Parameter $`a`$ and Violation of Vector Dominance

In this section we study the validity of vector dominance (VD) of electromagnetic form factor of the pion in hot matter. In Ref. Harada:2001rf  it has been shown that VD is accidentally satisfied in $`N_f=3`$ QCD at zero temperature and zero density, and that it is largely violated in large $`N_f`$ QCD when the VM occurs. At non-zero temperature there exists the hadronic thermal correction to the parameters. Thus it is nontrivial whether or not the VD is realized in hot matter, especially near the critical temperature. Here we will show that the intrinsic temperature dependences of the parameters of the HLS Lagrangian play essential roles, and then the VD is largely violated near the critical temperature.

We first study the direct $`\gamma \pi \pi `$ interaction at zero temperature. We expand the Lagrangian (5) in terms of the $`\pi `$ field with taking the unitary gauge of the HLS ($`\sigma =0`$) to obtain

$`_{(2)}`$ $`=`$ $`\text{tr}\left[_\mu \pi ^\mu \pi \right]+ag^2F_\pi ^2\text{tr}\left[\rho _\mu \rho ^\mu \right]+2i\left({\displaystyle \frac{1}{2}}ag\right)\text{tr}\left[\rho ^\mu [_\mu \pi ,\pi ]\right]`$ (98)
$`2\left(agF_\pi ^2\right)\text{tr}\left[\rho _\mu 𝒱^\mu \right]+2i\left(1{\displaystyle \frac{a}{2}}\right)\text{tr}\left[𝒱^\mu [_\mu \pi ,\pi ]\right]+\mathrm{},`$
where the vector meson field $`\rho _\mu `$ is introduced by

$$V_\mu =g\rho _\mu ,$$
(99)
and vector external gauge field $`𝒱_\mu `$ is defined as

$$𝒱_\mu \frac{1}{2}\left(_\mu +_\mu \right).$$
(100)
At the leading order of the derivative expansion in the HLS, the form of the direct $`\gamma \pi \pi `$ interaction is easily read from Eq. (98) as

$$\mathrm{\Gamma }_{\gamma \pi \pi (\mathrm{tree})}^\mu =e(qk)^\mu (1\frac{a}{2}),$$
(101)
where $`e`$ is the electromagnetic coupling constant and $`q`$ and $`k`$ denote outgoing momenta of the pions. This shows that, for the parameter choice $`a=2`$, the direct $`\gamma \pi \pi `$ coupling vanishes (the VD of the electromagnetic form factor of the pion).

At the next order there exist quantum corrections. In the background field gauge adopted in the present analysis the background fields $`\overline{𝒜}_\mu `$ and $`\overline{𝒱}_\mu `$ include the photon field $`A_\mu `$ and the background pion field $`\overline{\pi }`$ as

$`\overline{𝒜}_\mu `$ $`=`$ $`{\displaystyle \frac{1}{F_\pi (0)}}_\mu \overline{\pi }+{\displaystyle \frac{ie}{F_\pi (0)}}A_\mu [Q,\overline{\pi }]+\mathrm{},`$
$`\overline{𝒱}_\mu `$ $`=`$ $`eQA_\mu {\displaystyle \frac{i}{2F_\pi ^2(0)}}[_\mu \overline{\pi },\overline{\pi }]+\mathrm{},`$ (102)
where $`F_\pi (0)`$ is the on-shell pion decay constant (residue of the pion pole) in order to identify the field $`\overline{\pi }`$ with the on-shell pion field, and the charge matrix $`Q`$ is given by

$$Q=\left(\begin{array}{ccc}2/3& & \\ & 1/3& \\ & & 1/3\end{array}\right).$$
(103)
The direct $`\gamma \pi \pi `$ interaction including the next order correction is determined from the two-point functions of $`\overline{𝒱}_\mu `$-$`\overline{𝒱}_\nu `$ and $`\overline{𝒜}_\mu `$-$`\overline{𝒜}_\nu `$ and three-point function of $`\overline{𝒱}_\mu `$-$`\overline{𝒜}_\alpha `$-$`\overline{𝒜}_\beta `$. We can easily show that contribution from the three-point function vanishes in the low energy limit as follows: Let $`\mathrm{\Gamma }_{\overline{𝒱}\overline{𝒜}\overline{𝒜}}^{\mu \alpha \beta }`$ denote the $`\overline{𝒱}_\mu `$-$`\overline{𝒜}_\alpha `$-$`\overline{𝒜}_\beta `$ three-point function. Then the direct $`\gamma \pi \pi `$ coupling derived from this three-point function is proportional to $`q_\alpha k_\beta \mathrm{\Gamma }_{\overline{𝒱}\overline{𝒜}\overline{𝒜}}^{\mu \alpha \beta }`$. Since the legs $`\alpha `$ and $`\beta `$ of $`\mathrm{\Gamma }_{\overline{𝒱}\overline{𝒜}\overline{𝒜}}^{\mu \alpha \beta }`$ are carried by $`q`$ or $`k`$, $`q_\alpha k_\beta \mathrm{\Gamma }_{\overline{𝒱}\overline{𝒜}\overline{𝒜}}^{\mu \alpha \beta }`$ is generally proportional to two of $`q^2`$, $`k^2`$ and $`qk`$. Since the loop integral does not generate any massless poles, this implies that $`q_\alpha k_\beta \mathrm{\Gamma }_{\overline{𝒱}\overline{𝒜}\overline{𝒜}}^{\mu \alpha \beta }`$ vanishes in the low energy limit $`q^2=k^2=qk=0`$. Thus, the direct $`\gamma \pi \pi `$ interaction in the low energy limit can be read from the two-point functions of $`\overline{𝒱}_\mu `$-$`\overline{𝒱}_\nu `$ and $`\overline{𝒜}_\mu `$-$`\overline{𝒜}_\nu `$ as

$$\mathrm{\Gamma }_{\gamma \pi \pi }^\mu =e\frac{1}{F_\pi ^2(0)}\left[q_\nu \mathrm{\Pi }_{}^{\mu \nu }(q)k_\nu \mathrm{\Pi }_{}^{\mu \nu }(k)\frac{1}{2}(qk)_\nu \mathrm{\Pi }_{}^{\mu \nu }(p)\right],$$
(104)
where $`p_\nu =(q+k)_\nu `$ is the photon momentum. Substituting the decomposition of the two-point function given in Eq. (C.41) and taking the low-energy limit $`q^2=k^2=p^2=0`$, we obtain

$$\mathrm{\Gamma }_{\gamma \pi \pi }^\mu =e(qk)^\mu \left[1\frac{1}{2}\frac{\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2=0)}{F_\pi ^2(0)}\right],$$
(105)
where we used $`\mathrm{\Pi }_{}^{(\mathrm{vac})S}(q^2=0)=F_\pi ^2(0)`$. Comparing the above expression with the one in Eq. (101), we define the parameter $`a(0)`$ at one-loop level as

$$a(0)=\frac{\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2=0)}{F_\pi ^2(0)}.$$
(106)
We note that, in Ref. HYb , $`a(0)`$ is defined by the ratio $`F_\sigma ^2(M_\rho )/F_\pi ^2(0)`$ by neglecting the finite renormalization effect which depends on the details of the renormalization condition. While the above $`a(0)`$ in Eq. (106) defined from the direct $`\gamma \pi \pi `$ interaction is equivalent to the one used in Ref. HKRS . In the present renormalization condition (C.45), $`\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2=0)`$ is given by

$$\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2=0)=F_\sigma ^2(M_\rho )+\frac{N_f}{(4\pi )^2}M_\rho ^2\left(2\sqrt{3}\mathrm{tan}^1\sqrt{3}\right).$$
(107)
By adding this, the parameter $`a(0)`$ is expressed as

$$a(0)=\frac{F_\sigma ^2(M_\rho )}{F_\pi ^2(0)}+\frac{N_f}{(4\pi )^2}\frac{M_\rho ^2}{F_\pi ^2(0)}\left(2\sqrt{3}\mathrm{tan}^1\sqrt{3}\right).$$
(108)
Using $`M_\rho =771.1\text{MeV}`$, $`F_\pi (\mu =0)=86.4\text{MeV}`$ estimated in the chiral limit Gas:84 ; Gas:85a ; Gas:85b  <sup>#6</sup><sup>#6</sup>#6 In Ref. HYc , it was assumed that the scale dependent $`F_\pi (\mu )`$ agrees with the scale-independent parameter $`F_\pi `$ in the ChPT at $`\mu =0`$ to obtain $`F_\pi (\mu =0)=86.4\pm 9.7`$MeV. Here, we simply use this value to get a rough estimate of the value of the parameter $`a(0)`$ as done in Ref. HYc .  and $`F_\sigma ^2(M_\rho )/F_\pi ^2(0)=2.03`$ obtained through the Wilsonian matching for $`\mathrm{\Lambda }_{\mathrm{QCD}}=400\text{MeV}`$ and the matching scale $`\mathrm{\Lambda }=1.1\text{GeV}`$ in Ref. HYc , we estimate the value of $`a(0)`$ at zero temperature as

$$a(0)2.31.$$
(109)
This implies that the VD is well satisfied at $`T=0`$ even though the value of the parameter $`a`$ at the scale $`M_\rho `$ is close to one Harada:2001rf .

Now, let us study the direct $`\gamma \pi \pi `$ coupling in hot matter. In general, the electric mode and the magnetic mode of the photon couple to the pions differently in hot matter, so that there are two parameters as extensions of the parameter $`a`$. Similarly to the one obtained at $`T=0`$ in Eq. (104), in the low energy limit the direct $`\gamma \pi \pi `$ interaction derived from $`\overline{𝒜}_\mu `$-$`\overline{𝒜}_\nu `$ and $`\overline{𝒱}_\mu `$-$`\overline{𝒱}_\nu `$ two-point functions is expressed as

$`\mathrm{\Gamma }_{\gamma \pi \pi }^\mu (p;q,k)`$
$`={\displaystyle \frac{1}{\stackrel{~}{F}(\overline{q};T)\stackrel{~}{F}(\overline{k};T)}}[q_\nu \mathrm{\Pi }_{}^{\mu \nu }(q_0,\overline{q};T)k_\nu \mathrm{\Pi }_{}^{\mu \nu }(k_0,\overline{k};T)`$
$`{\displaystyle \frac{1}{2}}(qk)_\nu \mathrm{\Pi }_{}^{\mu \nu }(p_0,\overline{p};T)],`$ (110)
where $`q`$ and $`k`$ denote outgoing momenta of the pions and $`p=(q+k)`$ is the photon momentum. $`\stackrel{~}{F}`$ is the wave function renormalization of the background $`\overline{\pi }`$ field given in Eq. (72). Note that each pion is on its mass shell, so that $`q_0=v_\pi (\overline{q})\overline{q}`$ and $`k_0=v_\pi (\overline{k})\overline{k}`$. To define extensions of the parameter $`a`$, we consider the soft limit of the photon: $`p_00`$ and $`\overline{p}0`$<sup>#7</sup><sup>#7</sup>#7Note that we take the $`p_00`$ limit first and then take the $`\overline{p}0`$ limit for definiteness. We think that this is natural since the form factor in the vacuum is defined for space-like momentum of the photon.  Then the pion momenta become

$$q_0=k_0,\overline{q}=\overline{k}.$$
(111)
Note that while only two components $`\mathrm{\Pi }_{}^t`$ and $`\mathrm{\Pi }_{}^s`$ appear in $`q_\nu \mathrm{\Pi }_{}^{\mu \nu }`$ or $`k_\nu \mathrm{\Pi }_{}^{\mu \nu }`$, $`(qk)_\nu \mathrm{\Pi }_{}^{\mu \nu }`$ includes all four components $`\mathrm{\Pi }_{}^t`$, $`\mathrm{\Pi }_{}^s`$, $`\mathrm{\Pi }_{}^L`$ and $`\mathrm{\Pi }_{}^T`$. Since the tree part of $`\mathrm{\Pi }_{}^L`$ and $`\mathrm{\Pi }_{}^T`$ is $`2z_2p^2`$ which vanishes at $`p^2=0`$, it is natural to use only $`\mathrm{\Pi }_{}^t`$ and $`\mathrm{\Pi }_{}^s`$ to define the extensions of the parameter $`a`$. By including these two parts only, the temporal and the spatial components of $`\mathrm{\Gamma }_{\gamma \pi \pi }^\mu `$ are given by

$`\mathrm{\Gamma }_{\gamma \pi \pi }^0(0;q,q)`$ $`=`$ $`{\displaystyle \frac{2q_0}{\stackrel{~}{F}^2(\overline{q};T)}}\left[\mathrm{\Pi }_{}^t(q_0,\overline{q};T){\displaystyle \frac{1}{2}}\mathrm{\Pi }_{}^t(0,0;T)\right],`$
$`\mathrm{\Gamma }_{\gamma \pi \pi }^i(0;q,q)`$ $`=`$ $`{\displaystyle \frac{2q_i}{\stackrel{~}{F}^2(\overline{q};T)}}\left[\mathrm{\Pi }_{}^s(q_0,\overline{q};T){\displaystyle \frac{1}{2}}\mathrm{\Pi }_{}^s(0,0;T)\right].`$ (112)
Thus we define $`a^t(T)`$ and $`a^s(T)`$ as

$`a^t(\overline{q};T)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Pi }_{}^t(0,0;T)}{\mathrm{\Pi }_{}^t(q_0,\overline{q};T)}},`$
$`a^s(\overline{q};T)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Pi }_{}^s(0,0;T)}{\mathrm{\Pi }_{}^s(q_0,\overline{q};T)}}.`$ (113)
Here we should stress again that the pion momentum $`q_\mu `$ is on mass-shell: $`q_0=v_\pi (\overline{q})\overline{q}`$.

In the HLS at one-loop level the above $`a^t(\overline{q};T)`$ and $`a^s(\overline{q};T)`$ are expressed as

$`a^t(\overline{q};T)`$ $`=`$ $`a(0)\left[1+{\displaystyle \frac{\overline{\mathrm{\Pi }}_{}^t(0,0;T)a(0)\overline{\mathrm{\Pi }}_{}^t(\overline{q},\overline{q};T)}{a(0)F_\pi ^2(0;T)}}\right],`$ (114)
$`a^s(\overline{q};T)`$ $`=`$ $`a(0)\left[1+{\displaystyle \frac{\overline{\mathrm{\Pi }}_{}^s(0,0;T)a(0)\overline{\mathrm{\Pi }}_{}^s(\overline{q},\overline{q};T)}{a(0)F_\pi ^2(0;T)}}\right],`$ (115)
where $`a(0)`$ is defined in Eq. (106). From Eq. (39) together with Eq. (B.13) we obtain $`\overline{\mathrm{\Pi }}_{}^t`$ and $`\overline{\mathrm{\Pi }}_{}^s`$ in Eqs. (114) and (115) as

$`\overline{\mathrm{\Pi }}_{}^t(0,0;T)`$ $`=`$ $`\overline{\mathrm{\Pi }}_{}^s(0,0;T)`$ (116)
$`=`$ $`{\displaystyle \frac{N_f}{4}}\left[a^2\stackrel{~}{I}_2(T)\stackrel{~}{J}_1^2(M_\rho ;T)+2\stackrel{~}{J}_1^0(M_\rho ;T)\right].`$
Before going to the analysis near the critical temperature, let us study the temperature dependence of the parameters $`a^t(\overline{q};T)`$ and $`a^s(\overline{q};T)`$ in the low temperature region. At low temperature $`TM_\rho `$ the functions $`\stackrel{~}{J}_1^2(M_\rho ;T)`$ and $`\stackrel{~}{J}_1^0(M_\rho ;T)`$ are suppressed by $`e^{M_\rho /T}`$. Then the dominant contribution is given by $`\stackrel{~}{I}_2(T)=T^2/12`$. Combining this with Eq. (78) and (79), we obtain

$$a^ta^sa(0)\left[1+\frac{N_f}{12}\left(1\frac{a^2}{4a(0)}\right)\frac{T^2}{F_\pi ^2(0;T)}\right],$$
(117)
where $`a`$ is the parameter renormalized at the scale $`\mu =M_\rho `$, while $`a(0)`$ is defined in Eq. (106). We think that the intrinsic temperature dependences are small in the low temperature region, so that we use the values of parameters at $`T=0`$ to estimate the temperature dependent correction to the above parameters. By using $`F_\pi (0)=86.4\text{MeV}`$, $`a(0)2.31`$ given in Eq. (109) and $`a(M_\rho )=1.38`$ obtained through the Wilsonian matching for $`(\mathrm{\Lambda }_{\mathrm{QCD}},\mathrm{\Lambda })=(0.4,\mathrm{\hspace{0.17em}1.1})\text{GeV}`$ and $`N_f=3`$ HYc , $`a^t`$ and $`a^s`$ in Eq. (117) are evaluated as

$$a^ta^sa(0)\left[1+0.066\left(\frac{T}{50\text{MeV}}\right)^2\right].$$
(118)
This implies that the parameters $`a^t`$ and $`a^s`$ increase with temperature in the low temperature region. However, since the correction is small, we conclude that the vector dominance is well satisfied in the low temperature region.

At higher temperature the intrinsic thermal effects become more important. As we have shown in section 4, the parameters $`(g,a)`$ approach $`(0,1)`$ for $`TT_c`$ by the intrinsic temperature dependences, and then the parametric vector meson mass $`M_\rho `$ vanishes. As we have shown in section 6, near the critical temperature $`\mathrm{\Pi }_{}^t`$ and $`\mathrm{\Pi }_{}^s`$ in Eqs. (114) and (115) approach the following expressions:

$`\overline{\mathrm{\Pi }}_{}^t(\overline{q},\overline{q};T)\stackrel{TT_c}{}{\displaystyle \frac{N_f}{24}}T^2,`$
$`\overline{\mathrm{\Pi }}_{}^s(\overline{q},\overline{q};T)\stackrel{TT_c}{}{\displaystyle \frac{N_f}{24}}T^2.`$ (119)
On the other hand, the functions $`\mathrm{\Pi }_{}^t`$ and $`\mathrm{\Pi }_{}^s`$ in Eq. (116) at the limit of $`M_\rho /T0`$ and $`a1`$ become

$$\overline{\mathrm{\Pi }}_{}^t(0,0;T)=\overline{\mathrm{\Pi }}_{}^s(0,0;T)\frac{N_f}{2}\stackrel{~}{I}_2(T)=\frac{N_f}{24}T^2.$$
(120)
Furthermore, from Eq. (108), the parameter $`a(0)`$ approaches $`1`$ for $`M_\rho 0`$ and $`F_\sigma ^2(M_\rho )/F_\pi ^2(0)1`$:

$$a(0)1.$$
(121)
From the above limits in Eqs. (119), (120) and (121), the numerators of $`a^t(\overline{q};T)`$ and $`a^s(\overline{q};T)`$ in Eqs. (114) and (115) behave as

$`\overline{\mathrm{\Pi }}_{}^t(0,0;T)a(0)\overline{\mathrm{\Pi }}_{}^t(\overline{q},\overline{q};T)0,`$
$`\overline{\mathrm{\Pi }}_{}^s(0,0;T)a(0)\overline{\mathrm{\Pi }}_{}^s(\overline{q},\overline{q};T)0.`$ (122)
Thus we obtain

$$a^t(\overline{q};T),a^s(\overline{q};T)\stackrel{TT_c}{}1.$$
(123)
This implies that the vector dominance is largely violated near the critical temperature.

## 8 Summary and Discussions

In this paper we first showed the detailed calculations of the two-point functions in the background field gauge. Then, we showed how to extend the Wilsonian matching to the version at non-zero temperature. Throughout this paper, we assume that the chiral phase transition is of second or weakly first order so that the current correlators agree with each other at the critical temperature. Thus the Wilsonian matching leads to the following conditions:

$`g(\mathrm{\Lambda };T)\stackrel{TT_c}{}0,a(\mathrm{\Lambda };T)\stackrel{TT_c}{}1.`$ (124)
It is important that these conditions for the bare parameters are realized by the intrinsic thermal effect only, and further $`(g,a)=(0,1)`$ is a fixed point of RGEs. In order to be consistent with the chiral symmetry restoration in QCD, the axial vector current correlator must agree with the vector current correlator in the limit $`TT_c`$. This agreement is satisfied when we incorporate not only the hadronic corrections but also intrinsic effects into the physical quantities.

In section 5, we studied the vector meson pole mass including the hadronic thermal corrections calculated in the background field gauge. In the low temperature region $`TM_\rho `$, the vector meson mass increases as $`T^4`$ dominated by the pion loop. When we take the limit $`TT_c`$, it is crucial that the parameter $`M_\rho `$ goes to zero by the intrinsic thermal effect. Then we showed that the pole mass of the vector meson vanishes:

$`m_\rho (T)\stackrel{TT_c}{}0.`$
This implies that the vector manifestation is realized in hot matter as was first shown in Ref. HSasaki  where we used the hadronic thermal corrections calculated in the Landau gauge HS . We should stress that the conditions in Eq. (124) realized by the intrinsic thermal effects, which we call “the VM conditions in hot matter”, are essential for the VM to take place in hot matter. One might think that the VM is same as Georgi’s vector realization Georgi , in which the order parameter $`f_\pi `$ is non-zero although the chiral symmetry is unbroken. However in the VM the order parameter certainly becomes zero by the quadratic divergence and chiral symmetry is restored by massless vector meson. Therefore the VM is consistent with the Ward-Takahashi identity since it is the Wigner realization Yamawaki .

In section 6, we studied the temperature dependence of the pion velocity together with those of the temporal and spatial pion decay constants in the low temperature region. In hot matter, the pion velocity $`v_\pi `$ is not the speed of light because of the lack of Lorentz invariance. Actually we confirmed that

$`v_\pi ^2(\overline{p};T)1{\displaystyle \frac{N_f}{15}}{\displaystyle \frac{a\pi ^2}{F_\pi ^2M_\rho ^2}}T^4<1\text{for}TM_\rho .`$
We briefly reviewed how the temporal and spatial pion decay constants vanish at the critical point simultaneously following Ref. HKRS . Then we estimated the value of the critical temperature. From the order parameter $`f_\pi `$ including the intrinsic thermal effect as well as the hadronic correction, we obtained the following values for several choices of the matching scale ($`\mathrm{\Lambda }=0.8\text{-}\mathrm{\hspace{0.17em}1.1}`$ GeV) and the scale of QCD ($`\mathrm{\Lambda }_{\mathrm{QCD}}=0.30\text{-}0.45`$ GeV):

$`T_c=200\text{-}\mathrm{\hspace{0.17em}250}\text{MeV}.`$
We should note that the above values may be changed when we adopt a different way to estimate the matrix elements of operators in the OPE side, e.g., an estimation with the dilute pion gas approximation HKL  or that by the lattice QCD calculation lattice . Even when we choose one way to estimate the matrix elements in the OPE, some temperature effects are supposed to be left out due to the truncation of the OPE to neglect higher order operators, inclusion of which will cause a small change of the above values of the critical temperature. The important point is that as a result of the Wilsonian matching, $`T_c`$ is obtained as in Eq. (97) in terms of the quark and gluonic condensates, not hadronic degrees of freedom. It is expected that the value of $`T_c`$ may become smaller than that obtained in this paper by including the higher order corrections.

Finally in section 7, we presented a new prediction associated with the validity of vector dominance (VD) in hot matter. In the HLS model at zero temperature, the Wilsonian matching predicts $`a2`$ HYa ; HYc  which guarantees the VD of the electromagnetic form factor of the pion. Even at non-zero temperature, this is valid as long as we consider the thermal effects in the low temperature region, where the intrinsic temperature dependences are negligible. We showed that, as a consequence of including the intrinsic effect, the VD is largely violated at the critical temperature:

$`a^t(\overline{p};T)\stackrel{TT_c}{}1,a^s(\overline{p};T)\stackrel{TT_c}{}1.`$
In general, full temperature dependences include both hadronic and intrinsic thermal effects. Then there exist the violations of VD and universality of the $`\rho `$-coupling at generic temperature, although at low temperature the VD and universality are approximately satisfied.

In several analyses such as the one on the dilepton spectra in hot matter done in Ref Rapp-Wambach:00 , the VD is assumed to be held even in the high temperature region. Our result indicates that the assumption of the VD may need to be weakened, at least by some amount, for consistently including the effect of the dropping mass of the vector meson into the analysis.

Several comments are in order:

The unitarity of the scalar channel of $`\pi `$-$`\pi `$ scattering will be broken around 400-500 MeV in the present model although the existence of the $`\rho `$ meson helps a little, as was studied in Ref. Harada:1995dc . The unitarity of the vector sector, on the other hand, is not violated due to the existence of the vector mesons. Then, we do not consider the scalar sector in the present analysis, and we restrict ourselves to the quantities related to the vector and axial vector sector. In Ref. HYc , it was discussed that in the limit approaching the critical point, the scalar meson may decouple from the theory since the scalar meson belongs to the different chiral representation from both the pion and the longitudinal vector meson.

In section 5, we studied the pole mass of the vector meson at the critical temperature. It is also important to study the screening mass of the vector meson which is determined from the vector current correlator as done in the lattice calculation. It is not necessary that the correlation function of vector current be divergent at the critical temperature since $`\overline{q}\gamma ^\mu q`$ is not the order parameter of the spontaneous chiral symmetry breaking. The vector susceptibility $`\chi _V`$ is related to the correlation length of the vector current in the static limit and is in fact finite at the critical temperature as studied in Ref. HKRS , which is consistent with the result by the lattice QCD calculation Gottlieb:ac ; Allton:2002zi . This may imply that the screening mass of the vector meson is finite at the critical temperature.

In the present analysis we neglected the Lorentz symmetry violating effects at bare level since they are small as was shown in Ref. HKRS  (see also section 4). It is interesting to include such small corrections at bare level and study their effects to the physical quantities (see e.g., Refs. Sasaki:2003qj ; Harada:2003at ).

It is important to study which universality class the VM belongs to. The analysis of critical exponents gives the answer of this problem. We will study this in a future publication.

In this paper, we studied the physical quantities including both hadronic and intrinsic thermal effects at the limit $`TT_c`$ as well as in the low temperature region, $`TM_\rho `$. It is important to study the behavior away from $`T_c`$ since the properties of vector mesons in the temperature region from 100 MeV to 200 MeV are needed for the analysis of RHIC data. This will be done in future works.

## Acknowledgment

We would like to thank Doctor Youngman Kim, Professor Mannque Rho and Professor Koichi Yamawaki for useful discussions and comments. We are grateful to Professor Gerry Brown for critical reading of the manuscript. MH would like to thank Professor Dong-Pil Min for his hospitality during the stay in Seoul National University where part of this work was done. The work of MH was supported in part by the Brain Pool program (#012-1-44) provided by the Korean Federation of Science and Technology Societies. This work is supported in part by the 21st Century COE Program of Nagoya University provided by Japan Society for the Promotion of Science (15COEG01).

Appendices

## Appendix A Polarization Tensors at Finite Temperature

At non-zero temperature there exist four independent polarization tensors, $`u^\mu u^\nu ,(g^{\mu \nu }u^\mu u^\nu ),P_L^{\mu \nu }`$ and $`P_T^{\mu \nu }`$. The rest frame of medium is shown by $`u^\mu =(1,\stackrel{}{0})`$. $`P_L^{\mu \nu }`$ and $`P_T^{\mu \nu }`$ are given by Kap

$`P_{T\mu \nu }`$ $`=`$ $`g_{\mu i}\left(\delta _{ij}{\displaystyle \frac{\stackrel{}{p_i}\stackrel{}{p_j}}{\overline{p}^2}}\right)g_{j\nu },`$
$`P_{L\mu \nu }`$ $`=`$ $`\left(g_{\mu \nu }{\displaystyle \frac{p_\mu p_\nu }{p^2}}\right)P_{T\mu \nu },`$ (A.1)
where we define $`p^\mu =(p_0,\stackrel{}{p})`$ and $`\overline{p}=|\stackrel{}{p}|`$. They have the following properties:

$`P_{L\mu }^\mu =1,P_{T\mu }^\mu =2,`$
$`P_{L\mu \alpha }P_L^{\alpha \nu }=P_{L\mu }^{}{}_{}{}^{\nu },P_{T\mu \alpha }P_T^{\alpha \nu }=P_{T\mu }^{}{}_{}{}^{\nu },`$
$`P_{L\mu \alpha }P_T^{\alpha \nu }=0.`$ (A.2)
Let us decompose a tensor $`\mathrm{\Pi }^{\mu \nu }(p_0,\overline{p})`$ into

$$\mathrm{\Pi }^{\mu \nu }=u^\mu u^\nu \mathrm{\Pi }^t+(g^{\mu \nu }u^\mu u^\nu )\mathrm{\Pi }^s+P_L^{\mu \nu }\mathrm{\Pi }^L+P_T^{\mu \nu }\mathrm{\Pi }^T.$$
(A.3)
Each component is obtained as

$`\mathrm{\Pi }^t`$ $`=`$ $`\mathrm{\Pi }^{00}+{\displaystyle \frac{\stackrel{}{p}_i}{p_0}}\mathrm{\Pi }^{i0},`$
$`\mathrm{\Pi }^s`$ $`=`$ $`{\displaystyle \frac{\stackrel{}{p}_i}{\overline{p}}}\mathrm{\Pi }^{ij}{\displaystyle \frac{\stackrel{}{p}_j}{\overline{p}}}{\displaystyle \frac{p_0\stackrel{}{p}_i}{\overline{p}^2}}\mathrm{\Pi }^{i0},`$
$`\mathrm{\Pi }^L\mathrm{\Pi }^s`$ $`=`$ $`{\displaystyle \frac{\stackrel{}{p}_i}{\overline{p}}}\mathrm{\Pi }^{ij}{\displaystyle \frac{\stackrel{}{p}_j}{\overline{p}}}+{\displaystyle \frac{\stackrel{}{p}_i}{p_0}}\mathrm{\Pi }^{i0},`$
$`\mathrm{\Pi }^T\mathrm{\Pi }^s`$ $`=`$ $`{\displaystyle \frac{1}{2}}P_{T\mu \nu }\mathrm{\Pi }^{\mu \nu }.`$ (A.4)
## Appendix B Loop Integrals at Finite Temperature

In this appendix we list the explicit forms of the functions appearing in the hadronic thermal corrections, $`\overline{A}_0`$, $`\overline{B}_0`$ and $`\overline{B}^{\mu \nu }`$ \[see Eqs. (11)–(16) for definitions\] in various limits relevant to the present analysis.

The functions $`\overline{A}_0`$, which are independent of external momentum $`p_\mu `$, is given by

$`\overline{A}_0(M_\rho ;T)`$ $`=`$ $`\stackrel{~}{J}_1^2(M_\rho ;T),`$ (B.5)
$`\overline{A}_0(0;T)`$ $`=`$ $`\stackrel{~}{I}_2(T),`$ (B.6)
where $`\stackrel{~}{J}_1^2`$ and $`\stackrel{~}{I}_2`$ are defined in Eqs. (D.62) and (D.63).

We list the relevant limits of the functions $`\overline{B}_0(p_0,\overline{p};M_\rho ,0;T)`$ and $`\overline{B}^{\mu \nu }(p_0,\overline{p};M_\rho ,0;T)`$ which appear in the two-point function of $`\overline{𝒜}_\mu `$. When the pion momentum is taken as its on-shell, the function $`\overline{B}_0`$ becomes

$`\overline{B}_0(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$
$`=`$ $`{\displaystyle }{\displaystyle \frac{d^3k}{(2\pi )^3}}[{\displaystyle \frac{1}{2\omega _\rho }}{\displaystyle \frac{1}{e^{\omega _\rho /T}1}}\{{\displaystyle \frac{1}{(\omega _\rho \overline{p})^2(\omega _\pi ^p)^22iϵ(\omega _\rho \overline{p})}}+{\displaystyle \frac{1}{(\omega _\rho +\overline{p})^2(\omega _\pi ^p)^2+2iϵ(\omega _\rho +\overline{p})}}\}`$
$`+{\displaystyle \frac{1}{2\omega _\pi ^p}}{\displaystyle \frac{1}{e^{\omega _\pi ^p/T}1}}\{{\displaystyle \frac{1}{(\omega _\pi ^p\overline{p})^2\omega _\rho ^22iϵ(\omega _\pi ^p\overline{p})}}+{\displaystyle \frac{1}{(\omega _\pi ^p+\overline{p})^2\omega _\rho ^2+2iϵ(\omega _\pi ^p+\overline{p})}}\}],`$
where we put $`ϵ+0`$ to make the analytic continuation of the frequency $`p_0=i2\pi nT`$ to the Minkowski variable, and we defined

$$\omega _\rho =\sqrt{|\stackrel{}{k}|^2+M_\rho ^2},\omega _\pi ^p=|\stackrel{}{k}\stackrel{}{p}|.$$
(B.8)
Two components $`\overline{B}^t`$ and $`\overline{B}^s`$ of $`\overline{B}^{\mu \nu }`$ in the same limit are given by

$`\overline{B}^t(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$
$`=`$ $`{\displaystyle }{\displaystyle \frac{d^3k}{(2\pi )^3}}[{\displaystyle \frac{1}{2\omega _\rho }}{\displaystyle \frac{1}{e^{\omega _\rho /T}1}}\{{\displaystyle \frac{(2\omega _\rho \overline{p})^2}{(\omega _\rho \overline{p})^2(\omega _\pi ^p)^22iϵ(\omega _\rho \overline{p})}}+{\displaystyle \frac{(2\omega _\rho +\overline{p})^2}{(\omega _\rho +\overline{p})^2(\omega _\pi ^p)^2+2iϵ(\omega _\rho +\overline{p})}}`$
$`{\displaystyle \frac{\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})}{\overline{p}}}({\displaystyle \frac{2\omega _\rho \overline{p}}{(\omega _\rho \overline{p})^2(\omega _\pi ^p)^22iϵ(\omega _\rho \overline{p})}}{\displaystyle \frac{2\omega _\rho +\overline{p}}{(\omega _\rho +\overline{p})^2(\omega _\pi ^p)^2+2iϵ(\omega _\rho +\overline{p})}})\}`$
$`+{\displaystyle \frac{1}{2\omega _\pi ^p}}{\displaystyle \frac{1}{e^{\omega _\pi ^p/T}1}}\{{\displaystyle \frac{(2\omega _\pi ^p\overline{p})^2}{(\omega _\pi ^p\overline{p})^2\omega _\rho ^22iϵ(\omega _\pi ^p\overline{p})}}+{\displaystyle \frac{(2\omega _\pi ^p+\overline{p})^2}{(\omega _\pi ^p+\overline{p})^2\omega _\rho ^2+2iϵ(\omega _\pi ^p+\overline{p})}}`$
$`+{\displaystyle \frac{\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})}{\overline{p}}}({\displaystyle \frac{2\omega _\pi ^p\overline{p}}{(\omega _\pi ^p\overline{p})^2\omega _\rho ^22iϵ(\omega _\pi ^p\overline{p})}}{\displaystyle \frac{2\omega _\pi ^p+\overline{p}}{(\omega _\pi ^p+\overline{p})^2\omega _\rho ^2+2iϵ(\omega _\pi ^p+\overline{p})}})\}],`$
$`\overline{B}^s(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$
$`=`$ $`{\displaystyle }{\displaystyle \frac{d^3k}{(2\pi )^3}}[{\displaystyle \frac{1}{2\omega _\rho }}{\displaystyle \frac{1}{e^{\omega _\rho /T}1}}`$
$`\times \{{\displaystyle \frac{(2\stackrel{}{k}\stackrel{}{p}\overline{p}^2)^2}{\overline{p}^2}}({\displaystyle \frac{1}{(\omega _\rho \overline{p})^2(\omega _\pi ^p)^22iϵ(\omega _\rho \overline{p})}}+{\displaystyle \frac{1}{(\omega _\rho +\overline{p})^2(\omega _\pi ^p)^2+2iϵ(\omega _\rho +\overline{p})}})`$
$`{\displaystyle \frac{\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})}{\overline{p}}}({\displaystyle \frac{2\omega _\rho \overline{p}}{(\omega _\rho \overline{p})^2(\omega _\pi ^p)^22iϵ(\omega _\rho \overline{p})}}{\displaystyle \frac{2\omega _\rho +\overline{p}}{(\omega _\rho +\overline{p})^2(\omega _\pi ^p)^2+2iϵ(\omega _\rho +\overline{p})}})\}`$
$`+{\displaystyle \frac{1}{2\omega _\pi ^p}}{\displaystyle \frac{1}{e^{\omega _\pi ^p/T}1}}`$
$`\times \{{\displaystyle \frac{(2\stackrel{}{k}\stackrel{}{p}\overline{p}^2)^2}{\overline{p}^2}}({\displaystyle \frac{1}{(\omega _\pi ^p\overline{p})^2\omega _\rho ^22iϵ(\omega _\pi ^p\overline{p})}}+{\displaystyle \frac{1}{(\omega _\pi ^p+\overline{p})^2\omega _\rho ^2+2iϵ(\omega _\pi ^p+\overline{p})}})`$
$`+{\displaystyle \frac{\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})}{\overline{p}}}({\displaystyle \frac{2\omega _\pi ^p\overline{p}}{(\omega _\pi ^p\overline{p})^2\omega _\rho ^22iϵ(\omega _\pi ^p\overline{p})}}{\displaystyle \frac{2\omega _\pi ^p+\overline{p}}{(\omega _\pi ^p+\overline{p})^2\omega _\rho ^2+2iϵ(\omega _\pi ^p+\overline{p})}})\}].`$
We further take the $`M_\rho 0`$ limits of the above expressions. The limit of $`\overline{B}_0(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ includes the infrared logarithmic divergence $`\mathrm{ln}M_\rho ^2`$. However, it appears multiplied by $`M_\rho ^2`$ in $`\overline{\mathrm{\Pi }}_{}^t`$ and $`\overline{\mathrm{\Pi }}_{}^s`$, and the product $`M_\rho ^2\overline{B}_0(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ vanishes at $`M_\rho 0`$ limit:

$$\underset{M_\rho 0}{lim}M_\rho ^2\overline{B}_0(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)=0.$$
(B.11)
As for $`\overline{B}^{t,s}(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$, we obtain

$`\underset{M_\rho 0}{lim}\overline{B}^t(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ $`=`$ $`2\stackrel{~}{I}_2(T),`$
$`\underset{M_\rho 0}{lim}\overline{B}^s(p_0=\overline{p}+iϵ,\overline{p};M_\rho ,0;T)`$ $`=`$ $`2\stackrel{~}{I}_2(T).`$ (B.12)
In the static limit ($`p_00`$), the functions $`M_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,0;T)`$ and $`\overline{B}^t(p_0,\overline{p};M_\rho ,0;T)\overline{B}^L(p_0,\overline{p};M_\rho ,0;T)`$ become

$`\underset{p_00}{lim}M_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,0;T)`$
$`=M_\rho ^2{\displaystyle \frac{d^3k}{(2\pi )^3}\frac{1}{\omega _\rho ^2(\omega _\pi ^p)^2}\left[\frac{1}{\omega _\rho }\frac{1}{e^{\omega _\rho /T}1}\frac{1}{\omega _\pi ^p}\frac{1}{e^{\omega _\pi ^p/T}1}\right]},`$ (B.13)
$`\underset{p_00}{lim}\left[\overline{B}^t(p_0,\overline{p};M_\rho ,0;T)\overline{B}^L(p_0,\overline{p};M_\rho ,0;T)\right]`$
$`={\displaystyle \frac{d^3k}{(2\pi )^3}\frac{4}{\omega _\rho ^2(\omega _\pi ^p)^2}\left[\frac{\omega _\rho }{e^{\omega _\rho /T}1}\frac{\omega _\pi ^p}{e^{\omega _\pi ^p/T}1}\right]}.`$ (B.14)
Taking the low momentum limits of these expressions, we obtain

$`\underset{\overline{p}0}{lim}M_\rho ^2\overline{B}_0(p_0=0,\overline{p};M_\rho ,0;T)=\stackrel{~}{I}_2(T)\stackrel{~}{J}_1^2(M_\rho ;T),`$ (B.15)
$`\underset{\overline{p}0}{lim}\left[\overline{B}^t(p_0=0,\overline{p};M_\rho ,0;T)\overline{B}^L(p_0=0,\overline{p};M_\rho ,0;T)\right]`$
$`={\displaystyle \frac{4}{M_\rho ^2}}\left[\stackrel{~}{J}_1^2(M_\rho ;T)\stackrel{~}{I}_4(T)\right].`$ (B.16)
Moreover, in the VM limit ($`M_\rho 0`$), these are reduced to

$`\underset{\overline{p}0}{lim}M_\rho ^2\overline{B}_0(p_0=0,\overline{p};M_\rho ,0;T)\stackrel{M_\rho 0}{}0,`$
$`\underset{\overline{p}0}{lim}\left[\overline{B}^t(p_0=0,\overline{p};M_\rho ,0;T)\overline{B}^L(p_0=0,\overline{p};M_\rho ,0;T)\right]\stackrel{M_\rho 0}{}2\stackrel{~}{I}_2(T).`$ (B.17)
We consider the functions $`\overline{B}_0`$ and $`\overline{B}^{\mu \nu }`$ appearing in the two-point functions of $`\overline{𝒱}_\mu `$ and $`\overline{V}_\mu `$. $`\overline{B}^t`$ and $`\overline{B}^s`$ are calculated as

$`\overline{B}^t(p_0,\overline{p};0,0;T)=\overline{B}^s(p_0,\overline{p};0,0;T)=2\overline{A}_0(0;T)=2\stackrel{~}{I}_2(T),`$
$`\overline{B}^t(p_0,\overline{p};M_\rho ,M_\rho ;T)=\overline{B}^s(p_0,\overline{p};M_\rho ,M_\rho ;T)=2\overline{A}_0(M_\rho ;T)=2\stackrel{~}{J}_1^2(M_\rho ;T).`$
(B.18)
The function $`\overline{B}_0`$ is expressed as

$`\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$ (B.19)
$`=`$ $`{\displaystyle }{\displaystyle \frac{d^3k}{(2\pi )^3}}[{\displaystyle \frac{1}{2\omega _\rho }}{\displaystyle \frac{1}{e^{\omega _\rho /T}1}}\{{\displaystyle \frac{1}{(\omega _\rho p_0)^2(\omega _\rho ^p)^2}}+{\displaystyle \frac{1}{(\omega _\rho +p_0)^2(\omega _\rho ^p)^2}}\}`$
$`+{\displaystyle \frac{1}{2\omega _\rho ^p}}{\displaystyle \frac{1}{e^{\omega _\rho ^p/T}1}}\{{\displaystyle \frac{1}{(\omega _\rho ^pp_0)^2\omega _\rho ^2}}+{\displaystyle \frac{1}{(\omega _\rho ^p+p_0)^2\omega _\rho ^2}}\}],`$
where we define

$$\omega _\rho ^p=\sqrt{|\stackrel{}{k}\stackrel{}{p}|^2+M_\rho ^2}.$$
(B.20)
In the static limit $`(p_00)`$, the functions $`\overline{B}_0`$ and $`\overline{B}^L`$ are expressed as

$`\underset{p_00}{lim}M_\rho ^2\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)`$
$`=M_\rho ^2{\displaystyle \frac{d^3k}{(2\pi )^3}\frac{1}{\omega _\rho ^2(\omega _\rho ^p)^2}\left[\frac{1}{\omega _\rho }\frac{1}{e^{\omega _\rho /T}1}\frac{1}{\omega _\rho ^p}\frac{1}{e^{\omega _\rho ^p/T}1}\right]},`$ (B.21)
$`\underset{p_00}{lim}\overline{B}^L(p_0,\overline{p};M_\rho ,M_\rho ;T)`$
$`={\displaystyle }{\displaystyle \frac{d^3k}{(2\pi )^3}}{\displaystyle \frac{1}{\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})}}[{\displaystyle \frac{1}{\omega _\rho }}{\displaystyle \frac{1}{e^{\omega _\rho /T}1}}\{4\omega _\rho ^2\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})\}`$
$`{\displaystyle \frac{1}{\omega _\rho ^p}}{\displaystyle \frac{1}{e^{\omega _\rho ^p/T}1}}\{4(\omega _\rho ^p)^2+\stackrel{}{p}(2\stackrel{}{k}\stackrel{}{p})\}].`$ (B.22)
Taking the low momentum limit $`(\overline{p}0)`$, these expressions become

$`\underset{\overline{p}0}{lim}M_\rho ^2\overline{B}_0(p_0=0,\overline{p};M_\rho ,M_\rho ;T)={\displaystyle \frac{1}{2}}\left[\stackrel{~}{J}_1^0(M_\rho ;T)\stackrel{~}{J}_1^2(M_\rho ;T)\right],`$ (B.23)
$`\underset{\overline{p}0}{lim}\overline{B}^L(p_0=0,\overline{p};M_\rho ,M_\rho ;T)=2\left[M_\rho ^2\stackrel{~}{J}_1^0(M_\rho ;T)+2\stackrel{~}{J}_1^2(M_\rho ;T)\right].`$ (B.24)
In the VM limit $`(M_\rho 0)`$, these are reduced to

$`\underset{\overline{p}0}{lim}M_\rho ^2\overline{B}_0(p_0=0,\overline{p};M_\rho ,M_\rho ;T)\stackrel{M_\rho 0}{}0,`$ (B.25)
$`\underset{\overline{p}0}{lim}\overline{B}^L(p_0=0,\overline{p};M_\rho ,M_\rho ;T)\stackrel{M_\rho 0}{}4\stackrel{~}{I}_2(T).`$ (B.26)
On the other hand, the functions $`\overline{B}_0,\overline{B}^L\overline{B}^s`$ and $`\overline{B}^T\overline{B}^s`$ in the low momentum limit $`(\overline{p}0)`$ are given by

$`\underset{\overline{p}0}{lim}\overline{B}_0(p_0,\overline{p};M_\rho ,M_\rho ;T)={\displaystyle \frac{1}{2}}\stackrel{~}{F}_3^2(p_0;M_\rho ;T),`$ (B.27)
$`\underset{\overline{p}0}{lim}\left[\overline{B}^L(p_0,\overline{p};M_\rho ,M_\rho ;T)\overline{B}^s(p_0,\overline{p};M_\rho ,M_\rho ;T)\right]={\displaystyle \frac{2}{3}}\stackrel{~}{F}_3^4(p_0;M_\rho ;T),`$ (B.28)
$`\underset{\overline{p}0}{lim}\left[\overline{B}^T(p_0,\overline{p};M_\rho ,M_\rho ;T)\overline{B}^s(p_0,\overline{p};M_\rho ,M_\rho ;T)\right]={\displaystyle \frac{2}{3}}\stackrel{~}{F}_3^4(p_0;M_\rho ;T),`$ (B.29)
where the function $`\stackrel{~}{F}_3^n`$ is defined in Appendix D.

## Appendix C  Quantum Corrections

In this appendix we briefly summarize the quantum corrections to the two-point functions.

Let us define the functions $`B_0^{(\mathrm{vac})}`$, $`B^{(\mathrm{vac})\mu \nu }`$ and $`A_0^{(\mathrm{vac})}`$ by the following integrals HYc :

$`A_0^{(\mathrm{vac})}(M)`$ $`=`$ $`{\displaystyle \frac{d^nk}{i(2\pi )^4}\frac{1}{M^2k^2}},`$ (C.30)
$`B_0^{(\mathrm{vac})}(p;M_1,M_2)`$ $`=`$ $`{\displaystyle \frac{d^nk}{i(2\pi )^4}\frac{1}{[M_1^2k^2][M_2^2(kp)^2]}},`$ (C.31)
$`B^{(\mathrm{vac})\mu \nu }(p;M_1,M_2)`$ $`=`$ $`{\displaystyle \frac{d^nk}{i(2\pi )^4}\frac{(2kp)^\mu (2kp)^\nu }{[M_1^2k^2][M_2^2(kp)^2]}}.`$ (C.32)
In the present analysis it is important to include the quadratic divergences to obtain the RGEs in the Wilsonian sense. Here, following Refs. HY:conformal ; HYa ; HYc , we adopt the dimensional regularization and identify the quadratic divergences with the presence of poles of ultraviolet origin at $`n=2`$ Veltman . This can be done by the following replacement in the Feynman integrals:

$$\frac{d^nk}{i(2\pi )^n}\frac{1}{k^2}\frac{\mathrm{\Lambda }^2}{(4\pi )^2},\frac{d^nk}{i(2\pi )^n}\frac{k_\mu k_\nu }{\left[k^2\right]^2}\frac{\mathrm{\Lambda }^2}{2(4\pi )^2}g_{\mu \nu }.$$
(C.33)
On the other hand, the logarithmic divergence is identified with the pole at $`n=4`$:

$$\frac{1}{\overline{ϵ}}+1\mathrm{ln}\mathrm{\Lambda }^2,$$
(C.34)
where

$$\frac{1}{\overline{ϵ}}\frac{2}{4n}\gamma _E+\mathrm{ln}(4\pi ),$$
(C.35)
with $`\gamma _E`$ being the Euler constant.

By using the replacements in Eqs. (C.33) and (C.34), the integrals in Eqs. (C.30), (C.31) and (C.32) are evaluated as

$`A_0^{(\mathrm{vac})}(M)={\displaystyle \frac{\mathrm{\Lambda }^2}{(4\pi )^2}}{\displaystyle \frac{M^2}{(4\pi )^2}}\mathrm{ln}{\displaystyle \frac{\mathrm{\Lambda }^2}{M^2}},`$ (C.36)
$`B_0^{(\mathrm{vac})}(p^2;M_1,M_2)={\displaystyle \frac{1}{(4\pi )^2}}\left[\mathrm{ln}\mathrm{\Lambda }^21F_0(p^2;M_1,M_2)\right],`$ (C.37)
$`B^{(\mathrm{vac})\mu \nu }(p;M_1,M_2)`$
$`=g^{\mu \nu }{\displaystyle \frac{1}{(4\pi )^2}}\left[2\mathrm{\Lambda }^2M_1^2\mathrm{ln}{\displaystyle \frac{\mathrm{\Lambda }^2}{M_1^2}}M_2^2\mathrm{ln}{\displaystyle \frac{\mathrm{\Lambda }^2}{M_2^2}}(M_1^2M_2^2)F_A(p^2;M_1,M_2)\right]`$
$`\left(g^{\mu \nu }p^2p^\mu p^\nu \right){\displaystyle \frac{1}{(4\pi )^2}}\left[{\displaystyle \frac{1}{3}}\mathrm{ln}\mathrm{\Lambda }^2F_0(p^2;M_1,M_2)+4F_3(p^2;M_1,M_2)\right],`$ (C.38)
where $`F_0`$, $`F_A`$ and $`F_3`$ are defined by

$`F_0(s;M_1,M_2)`$ $`=`$ $`{\displaystyle _0^1}𝑑x\mathrm{ln}\left[(1x)M_1^2+xM_2^2x(1x)s\right],`$
$`F_A(s;M_1,M_2)`$ $`=`$ $`{\displaystyle _0^1}𝑑x(12x)\mathrm{ln}\left[(1x)M_1^2+xM_2^2x(1x)s\right],`$
$`F_3(s;M_1,M_2)`$ $`=`$ $`{\displaystyle _0^1}𝑑xx(1x)\mathrm{ln}\left[(1x)M_1^2+xM_2^2x(1x)s\right].`$ (C.39)
We consider the quantum corrections denoted by $`\mathrm{\Pi }^{(\mathrm{vac})\mu \nu }`$ \[see Eq. (22)\]. At tree level the two-point functions of $`\overline{𝒜}_\mu `$, $`\overline{𝒱}_\mu `$ and $`\overline{V}_\mu `$ are given by

$`\mathrm{\Pi }_{}^{(\mathrm{tree})\mu \nu }(p)`$ $`=`$ $`F_{\pi ,\mathrm{bare}}^2g^{\mu \nu }+2z_{2,\mathrm{bare}}(p^2g^{\mu \nu }p^\mu p^\nu ),`$
$`\mathrm{\Pi }_{}^{(\mathrm{tree})\mu \nu }(p)`$ $`=`$ $`F_{\sigma ,\mathrm{bare}}^2g^{\mu \nu }+2z_{1,\mathrm{bare}}(p^2g^{\mu \nu }p^\mu p^\nu ),`$
$`\mathrm{\Pi }_V^{(\mathrm{tree})\mu \nu }(p)`$ $`=`$ $`F_{\sigma ,\mathrm{bare}}^2g^{\mu \nu }{\displaystyle \frac{1}{g_{\mathrm{bare}}^2}}(p^2g^{\mu \nu }p^\mu p^\nu ),`$
$`\mathrm{\Pi }_V^{(\mathrm{tree})\mu \nu }(p)`$ $`=`$ $`F_{\sigma ,\mathrm{bare}}^2g^{\mu \nu }+z_{3,\mathrm{bare}}(p^2g^{\mu \nu }p^\mu p^\nu ).`$ (C.40)
Thus the one-loop contributions to $`\mathrm{\Pi }_{}^{(\mathrm{vac})\mu \nu }`$ give the quantum corrections to $`F_\pi ^2`$ and $`z_2`$. Similarly, each of the one-loop contributions to $`\mathrm{\Pi }_{}^{(\mathrm{vac})\mu \nu }`$, $`\mathrm{\Pi }_V^{(\mathrm{vac})\mu \nu }`$ and $`\mathrm{\Pi }_V^{(\mathrm{vac})\mu \nu }`$ includes the quantum corrections to two parameters shown above. For distinguishing the quantum corrections to two parameters included in the two-point function, it is convenient to decompose each two-point function as

$$\mathrm{\Pi }^{(\mathrm{vac})\mu \nu }(p)=\mathrm{\Pi }^{(\mathrm{vac})S}(p)g^{\mu \nu }+\mathrm{\Pi }^{(\mathrm{vac})LT}(p)(g^{\mu \nu }p^2p^\mu p^\nu ).$$
(C.41)
It should be noticed that, since we use the Lagrangian with Lorentz invariance, the form of the quantum corrections is Lorentz invariant. Then, the relation between four components given in Eqs. (17)-(20) and two components shown above are given by

$`\mathrm{\Pi }^{(\mathrm{vac})t}=\mathrm{\Pi }^{(\mathrm{vac})s}=\mathrm{\Pi }^{(\mathrm{vac})S},`$
$`\mathrm{\Pi }^{(\mathrm{vac})L}=\mathrm{\Pi }^{(\mathrm{vac})T}=\mathrm{\Pi }^{(\mathrm{vac})LT}.`$ (C.42)
Using the decomposition in Eq. (C.41), we identify $`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})S}(p^2)`$ with the quantum correction to $`F_\pi ^2`$, $`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})LT}(p^2)`$ with that to $`z_2`$, and so on. It should be noticed that the following relation is satisfied HYa ; HYc :

$$\mathrm{\Pi }_V^{(\mathrm{vac})S}(p^2)=\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2)=\mathrm{\Pi }_V^{(\mathrm{vac})S}(p^2).$$
(C.43)
Then the quantum correction to $`F_\sigma ^2`$ can be extracted from any of $`\mathrm{\Pi }_{}^{\mu \nu }`$, $`\mathrm{\Pi }_V^{\mu \nu }`$ and $`\mathrm{\Pi }_V^{\mu \nu }`$.

We note that in Refs. HYa ; HYc  the finite corrections of $`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})\mu \nu }`$ are neglected by assuming that they are small. In this paper we include these finite contributions in addition to the divergent corrections. As in Refs. HS ; HSasaki , we adopt the on-shell renormalization condition. They are expressed as <sup>#8</sup><sup>#8</sup>#8 In the framework of the ChPT with HLS, the renormalization point $`\mu `$ should be taken as $`\mu M_\rho `$ since the vector meson decouples at $`\mu =M_\rho `$. Below the scale $`M_\rho `$ the parameter $`F_\pi `$, which is expressed as $`F_\pi ^{(\pi )}(\mu )`$ in Refs. HYa ; HYc , runs by the quantum correction from the pion loop alone. Then $`F_\pi ^2(\mu =0)`$ in the right-hand-side of the renormalization condition Eq. (C.44) is defined by $`[F_\pi ^{(\pi )}(\mu =0)]^2`$. Due to the presence of the quadratic divergence $`F_\pi (M_\rho )`$ is not connected smoothly to $`F_\pi ^{(\pi )}(M_\rho )`$. The relation between them is expressed as HYa ; HYc  $`[F_\pi ^{(\pi )}(M_\rho )]^2=F_\pi ^2(M_\rho )+[N_f/(4\pi )^2][a(M_\rho )/2]M_\rho ^2`$, where $`a(\mu )`$ for $`\mu >M_\rho `$ is defined as $`a(\mu )F_\sigma ^2(\mu )/F_\pi ^2(\mu )`$. By using this relation, the renormalization condition (C.44) determines the condition for $`F_\pi ^2(M_\rho )`$. Strictly speaking, we should use $`F_\pi (M_\rho )`$ in the calculations in the present analysis. However, the difference between $`F_\pi (0)`$ and $`F_\pi (M_\rho )`$ inside the loop correction coming from the finite renormalization effect is of higher order, and we use $`F_\pi (0)`$ inside one-loop corrections in this paper. :

$`\text{Re}\left[\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2=0)\right]`$ $`=`$ $`F_\pi ^2(\mu =0),`$ (C.44)
$`\text{Re}\left[\mathrm{\Pi }_V^{(\mathrm{vac})S}(p^2=M_\rho ^2)\right]`$ $`=`$ $`F_\sigma ^2(\mu =M_\rho ),`$ (C.45)
$`\text{Re}\left[\mathrm{\Pi }_V^{(\mathrm{vac})LT}(p^2=M_\rho ^2)\right]`$ $`=`$ $`{\displaystyle \frac{1}{g^2(\mu =M_\rho )}},`$ (C.46)
where $`\mu `$ denotes the renormalization point. Using the above renormalization conditions, we can write the two-point functions as

$`\mathrm{\Pi }_{}^{(\mathrm{vac})S}(p^2)`$ $`=`$ $`F_\pi ^2(0)+\stackrel{~}{\mathrm{\Pi }}_{}^S(p^2),`$
$`\mathrm{\Pi }_V^{(\mathrm{vac})S}(p^2)`$ $`=`$ $`F_\sigma ^2(M_\rho )+\stackrel{~}{\mathrm{\Pi }}_V^S(p^2),`$
$`\mathrm{\Pi }_V^{(\mathrm{vac})LT}(p^2)`$ $`=`$ $`{\displaystyle \frac{1}{g^2(M_\rho )}}+\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2),`$ (C.47)
where $`\stackrel{~}{\mathrm{\Pi }}_{}^S(p^2)`$, $`\stackrel{~}{\mathrm{\Pi }}_V^S(p^2)`$ and $`\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2)`$ denote the finite renormalization effects. Their explicit forms are listed below. From the above renormalization conditions they satisfy

$$\stackrel{~}{\mathrm{\Pi }}_{}^S(p^2=0)=\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^S(p^2=M_\rho ^2)=\text{Re}\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2=M_\rho ^2)=0.$$
(C.48)
Thus in the present renormalization scheme, all the quantum effects for the on-shell parameters at leading order are included through the renormalization group equations.

In the following, using the above functions, we summarize the quantum corrections to two components $`\mathrm{\Pi }^{(\mathrm{vac})S}`$ and $`\mathrm{\Pi }^{(\mathrm{vac})LT}`$ of the two-point functions $`\mathrm{\Pi }_{}^{(\mathrm{vac})\mu \nu },\mathrm{\Pi }_{}^{(\mathrm{vac})\mu \nu },\mathrm{\Pi }_V^{(\mathrm{vac})\mu \nu }`$ and $`\mathrm{\Pi }_V^{(\mathrm{vac})\mu \nu }`$ which are defined in Eq. (C.41). For $`\overline{𝒜}_\mu `$-$`\overline{𝒜}_\nu `$ two-point function, we obtain

$`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})S}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}[{\displaystyle \frac{2a}{2}}\mathrm{\Lambda }^2+{\displaystyle \frac{3}{4}}aM_\rho ^2\mathrm{ln}\mathrm{\Lambda }^2+{\displaystyle \frac{1}{4}}aM_\rho ^2\mathrm{ln}M_\rho ^2`$
$`aM_\rho ^2\{1+F_0(p^2;M_\rho ,0)+{\displaystyle \frac{1}{4}}F_A(p^2;M_\rho ,0)\}],`$
$`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})LT}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}{\displaystyle \frac{a}{4}}\left[{\displaystyle \frac{1}{3}}\mathrm{ln}\mathrm{\Lambda }^2F_0(p^2;M_\rho ,0)+4F_3(p^2;M_\rho ,0)\right].`$ (C.49)
Corrections to $`\overline{𝒱}_\mu `$-$`\overline{𝒱}_\nu `$ two-point function are given by

$`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})S}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}\left[{\displaystyle \frac{a^2+1}{4}}\mathrm{\Lambda }^2+{\displaystyle \frac{3}{4}}M_\rho ^2\mathrm{ln}M_\rho ^2+M_\rho ^2\left\{{\displaystyle \frac{1}{4}}\mathrm{ln}M_\rho ^21F_0(p^2;M_\rho ,M_\rho )\right\}\right],`$
$`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})LT}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}{\displaystyle \frac{1}{8}}[{\displaystyle \frac{a^24a+5}{3}}\mathrm{ln}\mathrm{\Lambda }^2F_0(p^2;M_\rho ,M_\rho )`$ (C.50)
$`+4F_3(p^2;M_\rho ,M_\rho )4(2a)^2\mathrm{ln}M_\rho ^2].`$
As for $`\overline{V}_\mu `$-$`\overline{V}_\nu `$ we have

$`\mathrm{\Pi }_{V\text{(1-loop)}}^{(\mathrm{vac})S}(p^2)`$ $`=`$ $`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})S}(p^2),`$
$`\mathrm{\Pi }_{V\text{(1-loop)}}^{(\mathrm{vac})LT}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}[{\displaystyle \frac{a^287}{24}}\mathrm{ln}\mathrm{\Lambda }^2+4+{\displaystyle \frac{23}{8}}F_0(p^2;M_\rho ,M_\rho )`$ (C.51)
$`+{\displaystyle \frac{9}{2}}F_3(p^2;M_\rho ,M_\rho ){\displaystyle \frac{a^2}{8}}\{F_0(p^2;0,0)4F_3(p^2;0,0)\}].`$
Finally, corrections to $`\overline{V}_\mu `$-$`\overline{𝒱}_\nu `$ two-point function are expressed as

$`\mathrm{\Pi }_{V\text{(1-loop)}}^{(\mathrm{vac})S}(p^2)`$ $`=`$ $`\mathrm{\Pi }_{\text{(1-loop)}}^{(\mathrm{vac})S}(p^2),`$
$`\mathrm{\Pi }_{V\text{(1-loop)}}^{(\mathrm{vac})LT}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}{\displaystyle \frac{1}{8}}[{\displaystyle \frac{1+2aa^2}{3}}\mathrm{ln}\mathrm{\Lambda }^2a(2a)\mathrm{ln}M_\rho ^2`$ (C.52)
$`F_0(p^2;M_\rho ,M_\rho )+4F_3(p^2;M_\rho ,M_\rho )].`$
The renormalization conditions in Eqs. (C.44)- (C.46) lead to the following relations among the bare and renormalized parameters:

$`F_{\pi ,\mathrm{bare}}^2{\displaystyle \frac{N_f}{4(4\pi )^2}}\left[2(2a)\mathrm{\Lambda }^2+3aM_\rho ^2\mathrm{ln}\mathrm{\Lambda }^2\right]`$
$`=F_\pi ^2(0){\displaystyle \frac{N_f}{4(4\pi )^2}}aM_\rho ^2\left[3\mathrm{ln}M_\rho ^2+{\displaystyle \frac{1}{2}}\right],`$ (C.53)
$`F_{\sigma ,\mathrm{bare}}^2{\displaystyle \frac{N_f}{4(4\pi )^2}}\left[(a^2+1)\mathrm{\Lambda }^2+3M_\rho ^2\mathrm{ln}\mathrm{\Lambda }^2\right]`$
$`=F_\sigma ^2(M_\rho ){\displaystyle \frac{N_f}{4(4\pi )^2}}M_\rho ^2\left[3\mathrm{ln}M_\rho ^24(1\sqrt{3}\mathrm{tan}^1\sqrt{3})\right],`$ (C.54)
$`{\displaystyle \frac{1}{g_{\mathrm{bare}}^2}}{\displaystyle \frac{N_f}{(4\pi )^2}}{\displaystyle \frac{87a^2}{24}}\mathrm{ln}\mathrm{\Lambda }^2`$
$`={\displaystyle \frac{1}{g^2(M_\rho )}}{\displaystyle \frac{N_f}{8(4\pi )^2}}\left[{\displaystyle \frac{87a^2}{3}}\mathrm{ln}M_\rho ^2{\displaystyle \frac{1475a^2}{9}}+41\sqrt{3}\mathrm{tan}^1\sqrt{3}\right].`$ (C.55)
From these relations, we obtain the RGEs for the parameters $`F_\pi ,g`$ and $`a`$ as HYa

$`\mu {\displaystyle \frac{dF_{\pi }^{}{}_{}{}^{2}}{d\mu }}`$ $`=`$ $`{\displaystyle \frac{N_f}{2(4\pi )^2}}\left[3a^2g^2F_{\pi }^{}{}_{}{}^{2}+2(2a)\mu ^2\right],`$ (C.56)
$`\mu {\displaystyle \frac{da}{d\mu }}`$ $`=`$ $`{\displaystyle \frac{N_f}{2(4\pi )^2}}(a1)\left[3a(a+1)g^2(3a1){\displaystyle \frac{\mu ^2}{F_{\pi }^{}{}_{}{}^{2}}}\right],`$ (C.57)
$`\mu {\displaystyle \frac{dg^2}{d\mu }}`$ $`=`$ $`{\displaystyle \frac{N_f}{2(4\pi )^2}}{\displaystyle \frac{87a^2}{6}}g^4.`$ (C.58)
The finite renormalization effects of the two-point functions are expressed as

$`\stackrel{~}{\mathrm{\Pi }}_{}^S(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}aM_\rho ^2\left[\left(1{\displaystyle \frac{M_\rho ^2}{4p^2}}\right)\left\{1\left(1{\displaystyle \frac{M_\rho ^2}{p^2}}\right)\mathrm{ln}\left(1{\displaystyle \frac{p^2}{M_\rho ^2}}\right)\right\}{\displaystyle \frac{1}{8}}\right],`$ (C.59)
$`\stackrel{~}{\mathrm{\Pi }}_V^S(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{(4\pi )^2}}M_\rho ^2\left[\sqrt{3}\mathrm{tan}^1\sqrt{3}+2\sqrt{{\displaystyle \frac{4M_\rho ^2p^2}{p^2}}}\mathrm{tan}^1\sqrt{{\displaystyle \frac{p^2}{4M_\rho ^2p^2}}}\right],`$ (C.60)
$`\stackrel{~}{\mathrm{\Pi }}_V^{LT}(p^2)`$ $`=`$ $`{\displaystyle \frac{N_f}{8(4\pi )^2}}[{\displaystyle \frac{a^2}{3}}\mathrm{ln}\left({\displaystyle \frac{p^2}{M_\rho ^2}}\right)24(1{\displaystyle \frac{M_\rho ^2}{p^2}})+41\sqrt{3}\mathrm{tan}^1\sqrt{3}`$ (C.61)
$`{\displaystyle \frac{2(12M_\rho ^2+29p^2)}{p^2}}\sqrt{{\displaystyle \frac{4M_\rho ^2p^2}{p^2}}}\mathrm{tan}^1\sqrt{{\displaystyle \frac{p^2}{4M_\rho ^2p^2}}}].`$
## Appendix D Functions

In this appendix, we list the integral forms of the functions which appear in the expressions of the physical quantities and the several limits of the loop integrals shown in Appendix B. The functions $`\stackrel{~}{I}_n(T)`$ and $`\stackrel{~}{J}_m^n(M;T)`$ ($`n`$, $`m`$: integers) are given by

$`\stackrel{~}{I}_n(T)={\displaystyle \frac{d^3k}{(2\pi )^3}\frac{|\stackrel{}{k}|^{n3}}{e^{k/T}1}}={\displaystyle \frac{1}{2\pi ^2}}\widehat{I}_nT^n,`$
$`\widehat{I}_2={\displaystyle \frac{\pi ^2}{6}},\widehat{I}_4={\displaystyle \frac{\pi ^4}{15}},`$ (D.62)
$`\stackrel{~}{J}_m^n(M;T)={\displaystyle \frac{d^3k}{(2\pi )^3}\frac{1}{e^{\omega /T}1}\frac{|\stackrel{}{k}|^{n2}}{\omega ^m}},`$ (D.63)
where

$$\omega =\sqrt{k^2+M^2}.$$
(D.64)
The functions $`\stackrel{~}{F}_3^n(p_0;M;T)`$ and $`\stackrel{~}{G}_n(p_0;T)`$, which appear in the vector meson pole mass in section 5, are defined as

$`\stackrel{~}{F}_3^n(p_0;M;T)={\displaystyle \frac{d^3k}{(2\pi )^3}\frac{1}{e^{\omega /T}1}\frac{4|\stackrel{}{k}|^{n2}}{\omega (4\omega ^2p_{0}^{}{}_{}{}^{2})}},`$
$`\stackrel{~}{G}_n(p_0;T)={\displaystyle \frac{d^3k}{(2\pi )^3}\frac{|\stackrel{}{k}|^{n3}}{e^{k/T}1}\frac{4|\stackrel{}{k}|^2}{4|\stackrel{}{k}|^2p_{0}^{}{}_{}{}^{2}}}.`$ (D.65)