The shear and bulk relaxation times from the general correlation functions

TL;DR

This paper derives Kubo formulas for shear and bulk relaxation times and computes the shear viscosity relaxation time in massless bb theory using a diagrammatic approach, advancing understanding of quantum field relaxation processes.

Contribution

It introduces new Kubo formulas for relaxation times and provides the first diagrammatic calculation of shear viscosity relaxation time in massless bb theory.

Findings

Kubo formulas for shear and bulk relaxation times derived

Shear viscosity relaxation time calculated in massless bb theory

Results enhance understanding of quantum field relaxation dynamics

Abstract

In this paper we present two quantum field theoretical analyses on the shear and bulk relaxation times. First, we discuss how to find Kubo formulas for the shear and the bulk relaxation times. Next, we provide results on the shear viscosity relaxation time obtained within the diagrammatic approach for the massless $\lambda\phi^4$ theory.

Figures

Equations

$${\partial }_{\mu }{T}^{\mu \nu }$$

$$T^{\mu \nu }=ϵu^{\mu }{u}^{\nu }-{\Delta }^{\mu \nu }(P+\Pi )+{\pi }^{\mu \nu }$$

$$\Pi ={\Pi }_{NS}-{\tau }_{\Pi }\dot{\Pi },{\pi }^{\mu \nu }={\pi }_{NS}^{\mu \nu }-{\tau }_{\pi }{\dot{\pi }}^{⟨\mu \nu ⟩},$$

$$0$$

$$\displaystyle k_{\alpha}\big{(}\bar{G}^{\alpha\beta,\mu\nu}(k)-g^{\beta\mu}\langle\hat{T}^{\alpha\nu}\rangle-g^{\beta\nu}\langle\hat{T}^{\alpha\mu}\rangle+g^{\alpha\beta}\langle\hat{T}^{\mu\nu}\rangle\big{)}=0.$$

$${\bar{G}}_R^{xy,xy}(\omega ,k_y)=\frac{{\omega }^2(ϵ+g_T(k_y)+i\omega A(\omega ,k_y))}{k_y^2-i\omega /D(\omega ,k_y)-{\omega }^{2}B(\omega ,k_y)}-P,$$

$$\eta =\underset{\omega \to 0}{\lim }\underset{k_y\to 0}{\lim }\frac{1}{\omega }\text{Im}{\bar{G}}_R^{xy,xy}(\omega ,k_y),\eta {\tau }_{\pi }-\frac{\kappa }{2}=-\frac{1}{2}\underset{\omega \to 0}{\lim }\underset{k_y\to 0}{\lim }{\partial }_{\omega }^2\text{Re}{\bar{G}}_R^{xy,xy}(\omega ,k_y),$$

$${\bar{G}}_L(\omega ,\mathbf{k})=\frac{{\omega }^2(ϵ+P+{\omega }^{2}Q(\omega ,\mathbf{k}))}{{\mathbf{k}}^2-{\omega }^2/Z(\omega ,\mathbf{k})+i{\omega }^{3}R(\omega ,\mathbf{k})},$$

$$\zeta =\underset{\omega \to 0}{\lim }\underset{\mathbf{k}\to 0}{\lim }\frac{1}{\omega }\text{Im}{\bar{G}}_R^{PP}(\omega ,\mathbf{k}),\zeta {\tau }_{\Pi }=-\frac{1}{2}\underset{\omega \to 0}{\lim }\underset{\mathbf{k}\to 0}{\lim }{\partial }_{\omega }^2\text{Re}{\bar{G}}_R^{PP}(\omega ,\mathbf{k}).$$

Full text

1 Introduction

Over the course of last decades relativistic viscous hydrodynamics has been shown to successfully describe and explain the behavior of the strongly interacting matter produced in heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC), see Refs. [1, 2] and references therein. In general, any viscous fluid is characterized by a set of transport coefficients. These enter hydrodynamic equations as parameters and must be obtained from the underlying microscopic theory, either in quantum field theory or kinetic theory. Many approaches have been developed to study the first-order transport coefficients [3]-[11] for various systems. There are also some studies on second-order transport coefficients [12]-[22] but their quantum-field-theoretical determination does not seem to be complete.

We have already undertaken the task to study the shear and the bulk relaxation times consistently from first principles. The comprehensive analysis is shown in [23] and here only a brief summary is presented. Using general properties of Green functions and the gravitational Ward identity we first parametrize the stress-energy correlation functions to find their most general forms. Then Kubo formulas for the relaxation times are found in the hydrodynamic limits of the corresponding response functions. We also study shear effects in the massless scalar field theory $\lambda\phi^{4}$ and calculate the shear relaxation time within the real-time formalism.

2 Equations of viscous hydrodynamics

The behavior of a relativistic system, which is close to thermal equilibrium, can be well described by the viscous hydrodynamics, which is based on the energy-momentum conservation law

[TABLE]

where the energy-momentum tensor takes the form

[TABLE]

with $\epsilon$ being the energy density, $P$ - thermodynamic pressure, $u^{\mu}$ are the components of the flow velocity with the normalization condition $u^{\mu}u_{\mu}=1$, $\Delta^{\mu\nu}=g^{\mu\nu}-u^{\mu}u^{\nu}$ is the projection operator with $u_{\mu}\Delta^{\mu\nu}=0$, and the Minkowski metric is $g^{\mu\nu}=(1,-1,-1,-1)$. The terms $\Pi$ and $\pi^{\mu\nu}$ are the bulk viscous pressure and the shear stress tensor, respectively, which have well defined forms in the Navier-Stokes limit. Then the viscous corrections are characterized by the bulk viscosity $\zeta$ and shear viscosity $\eta$, respectively. In the second order formulation of viscous hydrodynamics, the response of medium to the thermodynamic forces is not instantaneous. The viscous corrections approach their corresponding Navier-Stokes forms within some characteristic time scales, which are the bulk and shear relaxation times, $\tau_{\Pi}$ and $\tau_{\pi}$. Consequently, the viscous corrections are subject to relaxation equations

[TABLE]

where $\Pi_{NS}$ and $\pi^{\mu\nu}_{NS}$ are the bulk pressure and the stress tensor in the Navier-Stokes approach, and we used the notation $A^{\langle\mu\nu\rangle}\equiv\Delta^{\mu\nu}_{\alpha\beta}A^{\alpha\beta}$ where $\Delta^{\mu\nu}_{\alpha\beta}\equiv(\Delta^{\mu}_{\alpha}\Delta^{\nu}_{\beta}+\Delta^{\mu}_{\beta}\Delta^{\nu}_{\alpha}-2/3\Delta^{\mu\nu}\Delta_{\alpha\beta})/2$. We will not consider the non-linear terms here. For more advanced studies on the hydrodynamic equations see, for example, [24].

In case when there are no other than energy and momentum currents occurring in the system, two hydrodynamic modes determine its dynamics. They are governed by the following dispersion relations

[TABLE]

where $\omega$ and ${\bf k}$ are the frequency and wavevector of the modes, $D_{T}=\eta/(\epsilon+P)$, $\gamma=\zeta/(\epsilon+P)$, and $v_{s}^{2}=\partial P/\partial\epsilon$ is the speed of sound squared. The dispersion relation (4) governs the propagation of the diffusion mode which occurs in the direction transverse to the the flow velocity. The sound mode is, in turn, given by the dispersion relation (2) and it is an effect of the small disturbances propagating longitudinally in the medium. Both dispersion relations are essential to determine the respective correlation functions of the stress-energy tensor.

3 Stress-energy correlation functions and Kubo formulas

Since viscous hydrodynamics is a manifestation of the linear response theory, the deviations of different observables are given in terms of corresponding equilibrium response functions. Therefore, the response functions carry dynamical information about the system. In general, these functions cannot be calculated exactly but one is able to parametrize their most general structures for the stress-energy tensor components using the following arguments. First, the real part of a correlation function is an even function of frequency and the imaginary part, which directly corresponds to the spectral function, must be an odd function of frequency. Next, since the stress-energy tensor represents at the same time the conserved currents and also generators of the space time evolution, its correlation functions must satisfy the gravitational Ward identity [25]

[TABLE]

Finally, the low-frequency and long-wavelength limits must be properly incorporated to ensure the correlation functions behave well in these limits.

Using aforementioned arguments we have parametrized the retarded Green function for the shear flow as

[TABLE]

where $g_{T}(k_{y})=P+g_{\pi\pi}(k_{y})$ with $g_{\pi\pi}(k_{y})$ being a function of $k_{y}$ of the order at least $O(k_{y}^{2})$ in the small $k_{y}$ limit. The functions $A$, $B$, and $D$ have the form $D(\omega,k_{y})=D_{R}(\omega,k_{y})-i\omega D_{I}(\omega,k_{y})$, where $D_{R}(\omega,k_{y})$ and $D_{I}(\omega,k_{y})$ are real-valued even functions of $\omega$ and $k_{y}$. The real parts $D_{R}$ and $B_{R}$ must have a non-zero limit as $\omega\rightarrow 0$ and $k_{y}\rightarrow 0$. All other parts of $A$, $B$, and $D$ must have finite limits as $\omega\rightarrow 0$ and $k_{y}\rightarrow 0$. By matching the denominator of the Eq. (7) to the dispersion relation (4) one is able to express $D_{R}$, $D_{I}$, $B_{R}$, etc in the hydrodynamic limits via the transport coefficients and obtain the following Kubo relations

[TABLE]

where $\kappa$ is the static susceptibility. Accordingly, both formulas have to be used to find the shear relaxation time in the leading order. The coefficient $\kappa$ was shown in [26, 27] to be $\kappa=\mathcal{O}(\lambda^{0}T^{2})$ in the weak coupling limit. On the other hand, both $\eta$ and $\tau_{\pi}$ are related to the mean free path which depends inversely on the coupling constant. Hence, $\kappa$ may be omitted in further analysis.

To examine the bulk viscosity and its relaxation time one needs to consider the response function to the longitudinal fluctuations. Employing the analogous procedure as in case of shear flow we obtain the most general form of the longitudinal fluctuation correlation function in the form

[TABLE]

where the functions $Q$, $R$, and $Z$ are constrained in a similar way as the functions $A$, $B$, and $D$ occurring in Eq. (7). When hydrodynamic limits are applied one obtains some Kubo formulas employing $\zeta$, $\eta$, $\tau_{\Pi}$, $\tau_{\pi}$, and $\kappa$. However, with the help of Eqs. (8) we can extract simple Kubo formulas involving only bulk viscosity and its relaxation time. These are then given in terms of the pressure-pressure correlation function

[TABLE]

These formulas constitute the way of computation of the bulk relaxation time.

4 Shear relaxation time in scalar field theory

It is elastic scatterings that constitutes the shear effects and therefore to evaluate the shear relaxation time within the scalar massless field theory $\lambda\phi^{4}$ one requires expanding the retarded Green function to the leading order. For the comprehensive calculation see [23].

Following Kubo relations one obtains integral equations for the imaginary and real parts of the retarded Green functions, which then are solved numerically and the shear relaxation time is extracted. It has been shown that both the shear viscosity and the relaxation time are related to the inverse of the thermal width, as expected. In the left panel of Fig. 1 the relaxation time is plotted as a function of the thermal mass over temperature and $m_{\textrm{th}}/T$ which should be identified with $\sqrt{\lambda}$ with $\lambda$ being the coupling constant. In the left panel of Fig. 1 we present the ratio $\langle\epsilon+P\rangle\tau_{\pi}/\eta$ as a function of $m_{\textrm{th}}/T$ showing that it changes between 6.11 up to 6.55 in the range shown. The results presented here are in agreement with those obtained within kinetic theory [12, 15].

5 Conclusions

In this article we presented a set of Kubo formulas which enable us to compute the shear and the bulk relaxation times. The main focus was to obtain the Kubo formula to compute the bulk relaxation time and the formula found is consistent with the one obtained via the metric perturbation analysis [22] but it is different from the one obtained from the projection operator method [21].

We also computed the shear relaxation time and the ratio $\langle\epsilon+P\rangle\tau_{\pi}/\eta$ within scalar massless theory and their values are consistent with previous studies. In general the ratio seems to vary between 5 up to 7 across different theories.

Acknowledgments

A. Czajka acknowledges the support from the program Mobility Plus of the Polish Ministry of Science and Higher Education. S.J. is supported in part by the Natural Sciences and Engineering Research Council of Canada. Discussions with C. Gale and G.D. Moore are very much appreciated.