Thermal contribution of unstable states

TL;DR

This paper investigates the thermodynamic role of unstable states within the Lee model, analyzing phase shift derivatives and spectral functions to accurately account for their contributions.

Contribution

It provides an exact analysis of unstable state contributions to thermodynamics using the Lee model, connecting phase shifts, spectral functions, and Green functions.

Findings

The phase shift derivative correctly identifies interaction contributions.

Spectral functions and phase shifts are explicitly connected in the Lee model.

The approach clarifies how unstable states influence thermodynamic properties.

Abstract

Within the framework of the Lee model, we analyze in detail the difference between the energy derivative of the phase shift and the standard spectral function of the unstable state. The fact that the model is exactly solvable allows us to demonstrate the construction of these observables from various exact Green functions. The connection to a formula due to Krein, Friedal, and Lloyd is also examined. We also directly demonstrate how the derivative of the phase shift correctly identifies the relevant interaction contributions for consistently including an unstable state in describing the thermodynamics.

Equations

$$A_{\rho }(E)\approx A_{\rho }^{BW}(E)=\frac{{\Gamma }_{\rho \pi \pi }}{(E-m_{\rho }{)}^2+\frac{{\Gamma }_{\rho \pi \pi }^2}{4}}.$$

$$P_{\mathrm{free}}=P_{\pi }^{(0)}+P_{\rho }^{(0)},$$

$$P_{\rho }^{(0)}\to P_{\rho }=\int \frac{dE}{2\pi }{A}_{\rho }(E)\times P^{(0)}(E,T)$$

$$P_{Sch.A}=P_{\pi }^{(0)}+P_{\rho }.$$

$$\begin{array}{rl}{P}_{\mathrm{S}-\mathrm{matrix}}& =P_{\pi }^{(0)}+\int \frac{dE}{2\pi }B(E)\times P^{(0)}(E,T)\\ B(E)& =2\frac{d}{dE}{\mathcal{Q}}_{\pi \pi }(E).\end{array}$$

$$\begin{array}{rl}{\delta }_{\pi \pi }(E)& \to \pi \theta (E-m_{\rho })\\ B(E)& \to 2\pi \delta (E-m_{\rho }).\end{array}$$

$$B(E)=A_{\rho }(E)+\sum _q\Delta A_{2\pi }(E;q),$$

$$P_{\mathrm{S}-\mathrm{matrix}}=P_{\pi }^{(0)}+P_{\rho }+\Delta P_{2\pi },$$

$$\Delta P_{2\pi }=\int \frac{dE}{2\pi }\sum _q\Delta A_{2\pi }(E;q)\times P^{(0)}(E,T).$$

$$|\rho ⟩,\{|q⟩\},$$

$$\begin{array}{r}{H}_0|\rho ⟩={\Delta }_0|\rho ⟩,\\ H_0|q⟩=ϵ(q)|q⟩.\end{array}$$

$$V=\sum _q{g}_{\mathrm{eff}}C(q)\left(|\rho ⟩⟨q⟩+|q⟩⟨\rho ⟩\right),$$

$$\begin{array}{rl}V|\rho ⟩& =\sum _{q^{\prime }}{g}_{\mathrm{eff}}C(q^{\prime })|q^{\prime }⟩\\ V|q⟩& =g_{\mathrm{eff}}C(q)|\rho ⟩.\end{array}$$

$$\begin{array}{r}C(q{)}^2=\frac{4\pi q^2}{(2\pi {)}^3}\delta q\end{array}$$

$$H=H_0+V.$$

$$\begin{array}{rl}{G}_0(E)& =\frac{1}{E-H_0+iϵ}\\ G(E)& =\frac{1}{E-H+iϵ},\end{array}$$

$$\begin{array}{rl}G& =G_0+G_{0}VG\\ & =G_0+G_{0}T{G}_0,\end{array}$$

$$\begin{array}{rl}T& =V+V{G}_{0}T\\ & =V+VGV\\ & =VG{G}_0^{-1}.\end{array}$$

$$\begin{array}{rl}\hat{S}& =G_0^{\ast }{G^{\ast }}^{-1}G{G_0}^{-1}\\ & =I+\left(G_0-G_0^{\ast }\right)\times VG{G}_0^{-1}\\ & =I-2\pi i\delta (E-H_0)\times T.\end{array}$$

$$\begin{array}{rl}\mathcal{K}& =-\mathrm{Im}\left[(G-G^{\ast })-(G_0-G_0^{\ast })\right]\\ & =A-A_0\\ & =\Delta A,\end{array}$$

$$\begin{array}{rl}A(E)& =-2\mathrm{Im}G=-\mathrm{Im}(G-G^{\ast })\\ A_0(E)& =-2\mathrm{Im}{G}_0=-\mathrm{Im}(G_0-G_0^{\ast }).\end{array}$$

$$B(E)=2\frac{\partial }{\partial E}\mathcal{Q}(E).$$

$$\begin{array}{rl}{G}_{\rho }& \equiv ⟨\rho |G|\rho ⟩\\ & =⟨\rho |G_0|\rho ⟩+⟨\rho |G_{0}V{G}_0|\rho ⟩+⟨\rho |G_{0}V{G}_{0}V{G}_0|\rho ⟩+\dots ,\end{array}$$

$$\begin{array}{rl}⟨\rho |G_0|\rho ⟩& =G_{\rho }^0=\frac{1}{E-{\Delta }_0+iϵ}\\ ⟨\rho |G_{0}V{G}_0|\rho ⟩& =0,\end{array}$$

$$⟨\rho |G_{0}V{G}_{0}V{G}_0|\rho ⟩$$

$$\begin{array}{rl}{G}_{\rho }& =G_{\rho }^0+G_{\rho }^0{\Sigma }_{\rho }{G}_{\rho }^0+G_{\rho }^0{\Sigma }_{\rho }{G}_{\rho }^0{\Sigma }_{\rho }{G}_{\rho }^0+\dots \\ & =\frac{1}{(G_{\rho }^0{)}^{-1}-{\Sigma }_{\rho }},\end{array}$$

$$\begin{array}{rl}{\Sigma }_{\rho }& \equiv ⟨\rho |V{G}_{0}V|\rho ⟩\\ & =\sum _{q}C(q{)}^2{g}_{\mathrm{eff}}^2{G}_{2\pi }^0(E;q).\end{array}$$

$$\begin{array}{rl}\frac{\partial }{\partial E}{\Sigma }_{\rho }& =-1\times \sum _{q}C(q{)}^2{g}_{\mathrm{eff}}^2{G}_{2\pi }^0(E;q)G_{2\pi }^0(E;q).\end{array}$$

$$\begin{array}{rl}{G}_{2\pi }^0(E;q)& =\frac{1}{E-ϵ(q)+iϵ}\\ ⟹\frac{\partial }{\partial E}{G}_{2\pi }^0(E;q)& =-G_{2\pi }^0(E;q)\times G_{2\pi }^0(E;q)\end{array}$$

$$\begin{array}{rl}{G}_{2\pi }(E;q)& =⟨q|G|q⟩\\ & =G_{2\pi }^0(E;q)+⟨q|G_{0}T{G}_0|q⟩\\ & =G_{2\pi }^0(E;q)+G_{2\pi }^0(E;q)\times ⟨q|T|q⟩\times G_{2\pi }^0(E;q).\end{array}$$

Full text

I Introduction

Formal treatment of the interactions in a gas of particles at finite temperature is an important topic in thermal field theory Bloch et al. (2008); Campa et al. (2009); Satz (2012). In particular, a consistent description of the unstable states is imperative for understanding the hadron gas, see, e.g., Refs. Satz (2012); Florkowski (2010); Dashen et al. (1969); Venugopalan and Prakash (1992). Some questions of interest include: How are the bulk properties of the medium, such as pressure and energy density, affected by unstable particles? What are the effective ways to take these into consideration? What insights can be gained from comparing different approaches, e.g., the standard imaginary time formalism and those based on a virial expansion? Dashen et al. (1969); Osborn and Tsang (1976); Fré ; Chaichian et al. (1994); LeClair (2007); How and LeClair (2010); Liu (2013); Lo (2017) In this work we tackle these issues using an effective Hamiltonian approach. Besides an intuitive modification driven by the spectral function of the unstable particle (the resonant contribution), we shall demonstrate how the very presence of an interaction modifies the 2-body states composed of the stable particles (the nonresonant contribution). The sum of these modifications recovers a well-known result Beth and Uhlenbeck (1937); Dashen et al. (1969); Weinhold et al. (1998), according to which the derivative of the scattering phase shift can be identified as the density of state for computing the partition function.

Many thermal models have the widths of the resonances implemented but the nonresonant interactions are neglected. This can lead to misleading results when interpreting the contribution from an interaction channel Broniowski et al. (2015); Friman et al. (2015). An illustrative example is the pion-pion scattering in the $I=0$ channel, in which the famous $f_{0}(500)$ resonance, a.k.a. the $\sigma$-meson, is involved. The empirical phase shift, analyzed by the chiral perturbation theory Oller et al. (1999), reveals that there are substantial (effective) repulsive corrections coming from the exchange interactions in the t- and u-channels. Additional cancellation also comes from the $I=2$ channel. A model that corrects only for the width of the resonance is incapable of handling these effects. Some models try to remedy this by introducing extraneous repulsive forces, e.g., via an excluded volume. This, however, will generally lead to a model which contradicts the known phase shift Lo et al. (2017).

In this study we consider a system of stable particles “$\pi$” (the pions) and unstable particle “$\rho$” (the $\rho$-mesons). 111 These notions obviously mirror the case of a thermal hadron gas with pions and $\rho(770)$’s. However, our discussion generally applies to any thermal system with unstable states. Each $\rho$ can decay into two $\pi$’s via the interaction $\rho\rightarrow\pi\pi$. Chains of interactions, e.g., $\pi\pi\rightarrow\rho\rightarrow\pi\pi$, are also included. The fact that only two types of particles are considered drastically simplifies the discussion, but the main non-trivial features are kept.

We use a Lee-Model Hamiltonian (LH) Lee (1954); Chiu et al. (1992) to describe the interactions in the $\rho\pi$ system. Similar Hamiltonians have been used to explore various areas in physics, ranging from atomic physics and quantum optics Jaynes and Cummings (1963); Babelon and Talalaev (2007); Babelon et al. (2009); Scully and Zubairy (1997) to baryon decays Liu et al. (2016). The Lee model offers a useful theoretical setup to study an interacting system Facchi et al. (2001, 2000); Facchi and Pascazio (2001, 1998); Berman and Ford (2010); Kofman et al. (1994); Giacosa (2013) and in many aspects it resembles a quantum field theory Giacosa (2012). In practice, the LH contains an unstable state $\left|\rho\right\rangle$, as well as a continuum of 2-body states $\left|\pi(\vec{q})\pi(-\vec{q})\right\rangle$ with all possible relative momenta $\vec{q}$ of the pair. The decay $\rho\rightarrow\pi\pi$ is described by mixing terms of the form: $\sim g\left(\left|\rho\right\rangle\left\langle\pi\pi\right|+\left|\pi\pi\right\rangle\left\langle\rho\right|\right)$. This gives the unstable particle $\rho$ a width, i.e., a distribution in energies (or masses) dictated by the spectral function $A_{\rho}(E)$ Giacosa (2012); Matthews and Salam (1958, 1959); Giacosa and Pagliara (2007), such that $A_{\rho}(E)\,\frac{dE}{2\pi}$ can be interpreted as the probability that the unstable $\rho$ has energy between $E$ and $E+dE$. The spectral function can be calculated from the imaginary part of the $\rho$-propagator. See Sec. III for details. For a narrow-width state this can be approximated by the Breit-Wigner (BW) formula Weisskopf and Wigner (1930a, b),

[TABLE]

Note that the width $\Gamma_{\rho\pi\pi}$ is generally energy dependent, and the mass $m_{\rho}$ can be modified by the real part of quantum loops. Thus, the decay probability is never exactly exponential. See, e.g., the theoretical treatment in Refs. Facchi et al. (2001); Fonda et al. (1978); Ersak and the experimental results in Refs. Wilkinson et al. (1997); Fischer et al. (2001); Rothe et al. (2006); and Ref. Giacosa and Pagliara (2011) for a discussion in a quantum field theory.

Based on the LH, the finite temperature properties of the system can be derived using the standard techniques of statistical mechanics. See Sec. IV for details. Here we give a synopsis of our discussion. Consider the hypothetical limit where $\rho$ and $\pi$’s are completely decoupled, i.e., $g=0$ (or $\Gamma_{\rho\pi\pi}=0$). The pressure of the system would be given by a sum of two contributions:

[TABLE]

where P${}^{(0)}_{a}=P^{(0)}(E=m_{a},T)$ denotes the pressure of an ideal gas of a species $a$, which depends on its mass ($m_{a}$) and degeneracy. The dependence on temperature $T$ is understood. Even if the interaction is switched on, Eq. (2) can still be used as an estimate of the pressure for a narrow-width $\rho$. This is the fundamental premise of the hadron resonance gas (HRG) model Hagedorn (1965); Andronic et al. (2018a): contribution of resonances to the thermodynamics is approximated by an uncorrelated gas of zero-width particles.

As a next step in improving the approximation, we take into account the width of $\rho$ via a weighed sum by $A_{\rho}(E)$:

[TABLE]

such that the total pressure is approximated as (Scheme-A):

[TABLE]

This scheme is employed in many versions of the HRG models, see e.g., Refs. Satz (2012); Florkowski (2010); Andronic et al. (2009); Alba et al. (2014); Torrieri et al. (2005); Vovchenko et al. (2018). See also the K-matrix-based approach Doring and Koch (2007); Dash et al. (2018). However, Eq. (4) is not yet complete. According to the S-matrix formulation of statistical mechanics by Dashen et al. Dashen et al. (1969) (see also the discussion by Weinhold et al. Weinhold et al. (1998)), the correct result of the pressure at arbitrary $g$ (or $\Gamma_{\rho\pi\pi}$) is given by (S-matrix scheme):

[TABLE]

where ${\mathcal{Q}}_{\pi\pi}(E)$ is the phase shift for the scattering process $\pi\pi\rightarrow\pi\pi$. Here we summarize some key features of the S-matrix scheme:

• (i)

Observe that there is no explicit $\rho$ contribution in Eq. (5): The pressure is determined based on the scattering information of the asymptotic (stable) states alone. In fact, it is not compulsory to introduce the $\rho$ state as an explicit degree of freedom. Its presence is encoded in the phase shift. This point will be made clear by direct model calculations.

• (ii)

Eqs. (4) and (5) reduce to the free gas result (2) in the limit of $g\rightarrow 0$ ( or $\Gamma_{\rho\pi\pi}\rightarrow 0$ ). 222The $g=0$ limits of $A_{\rho}(E)$ and $B(E)$ are not well-defined. For the latter, we have

[TABLE]

• (iii)

Generally, $B(E)\neq A_{\rho}(E)$, and Eqs. (4) and  (5) are thus different. Systems which show substantial deviation are plenty: In addition to the case of $f_{0}(500)$ mentioned, nonresonant contribution is found to be important in the study of $\kappa(700)$ Broniowski et al. (2015); Friman et al. (2015), the $N^{*}$ and $\Delta$ resonances Lo et al. (2018); Andronic et al. (2018b), the $S=-1$ hyperons Fernandez-Ramirez et al. (2018), etc. It is also the case for the recently discovered $X,$ $Y,$ $Z$ states Lebed et al. (2017); Chen et al. (2016); Esposito et al. (2015). As shown in a recent work Ortega et al. (2018), the state $X(3872)$ makes only a small contribution to the thermodynamics due to nonresonant effects. For what concerns other states, future studies based on Eq. (5) are needed.

In this work we verify Eq. (5), instead of Eq. (4), gives the correct description of the thermodynamics of an interacting system. This point had been raised in previous works, see e.g., Ref. Weinhold et al. (1998), but the actual adoption of the scheme remains limited Broniowski et al. (2003); Andronic et al. (2018b). We hope that a more detailed account of the different spectral functions can raise the awareness of the issue in the community and further promote the use of the correct formula. We also establish their formal relations to the resolvent and the density of states. This gives an interesting perspective in describing the thermodynamics of an interacting system.

Using an LH, we derive the mismatch analytically, the result takes the form

[TABLE]

where the second term in the R.H.S. describes the the modification of the spectral function of the 2-body $|\pi\pi\rangle$ state. Such a term is present even in the absence of the resonance, e.g., taking the large $m_{\rho}$ limit. (See Sec. III) The correct expression of the pressure can be decomposed as:

[TABLE]

where

[TABLE]

As we shall see, the last term is in general not negligible and can even be dominant at low temperatures.

The paper is organized as follows: In Sec. II the details of the Lee model are presented. Then, in Sec. III, the spectral functions for both $\rho$ and $\pi\pi$, and the phase shift are introduced. Here the important Eq. (7) is derived. In Sec. IV the thermodynamic properties of the system are determined analytically, with special focus on the pressure with its various contributions. A numerical example shows that $\Delta P_{2\pi}$ can be sizable and in general should not be neglected. Finally, discussions and conclusions are given in Sec. V.

II The Lee Model

The Lee model Lee (1954) describing the $\rho\leftrightarrow\pi\pi$ system can be formulated as follows. 333 In the following we measure energy with respect to $2m_{\pi}$, and the nonrelativistic dispersion relation $\epsilon(q)={q^{2}}/{m_{\pi}}$ is used. We also choose to present our model in a discretized form. This prepares for the later numerical treatment of solving the system on a momentum grid. Hall et al. (2013). It is easy to go to the continuum by taking $\sum_{q}C(q)^{2}\,(\cdots)\rightarrow\int\frac{d^{3}q}{(2\pi)^{3}}\,(\cdots)$.

Introducing the basis states in the center-of-mass (CM) frame:

[TABLE]

where $q$ is the momentum label for the two-pion state $|\pi(\vec{q})\pi(-\vec{q})\rangle$. The Hamiltonian of the system can be represented as an $(1+N_{q})\times(1+N_{q})$-matrix. The non-interacting Hamiltonian $H_{0}$ is a diagonal matrix with

[TABLE]

The interaction $V$ describes the coupling of $\rho$ with the $|q\rangle$ states

[TABLE]

such that

[TABLE]

We use

[TABLE]

to implement the spherical degeneracy. Note that $\delta q=\frac{2\pi}{L}$ in the finite volume formulation, $L$ being the size of the box. The coupling $g_{\rm eff}$ is generally $q$-dependent. The full Hamiltonian then reads

[TABLE]

With the Hamiltonian defined, we can construct the resolvent operators

[TABLE]

which can be understood again as $(1+N_{q})\times(1+N_{q})$ matrices. In the remainder of this paper we will suppress the $E$ dependence unless there is a chance for confusion.

The well-known relations from the Lippmann-Schwinger equation can also be directly realized:

[TABLE]

and

[TABLE]

In this paper, we investigate the inclusion of an unstable state in the description of thermodynamics using the S-matrix formulation. The key operators of interest in this scheme is the scattering operator Lo (2017); Taylor (2012)

[TABLE]

Since $\rho$ is not an asymptotic state, the actual scattering matrix $S$ is extracted from the (lower-right) $N_{q}\times N_{q}$ block of $\hat{S}$. In addition, we introduce an operator $\mathcal{K}$, due to Krein, Friedal and Lloyd (KFL) Cvitanović et al. (2016); Texier (2016), defined as the difference of the spectral functions

[TABLE]

where we have identified the spectral function operator

[TABLE]

These operators are deeply connected with the scattering phase shift $\mathcal{Q}$ and the effective spectral function $B$. The latter is defined as

[TABLE]

There are multiple ways to extract the phase shift $\mathcal{Q}$ from the resolvents (or Green functions). (See Eqs. (35), (47)). These and the explicit relations among the various observables will be demonstrated in the context of the Lee model in Sec. III.

III Phase shift and effective spectral function

In the Lee model various theoretical quantities, e.g. the propagator and the self-energy of $\rho$, can be analytically computed. It is a useful exercise to revisit these formulas as it helps to build an understanding the physical content of the KFL operator $\mathcal{K}$ and the derivative of the phase shift.

III.1 The $\rho$-propagator

The $\rho$-propagator $G_{\rho}$ can be computed from $\langle\rho|G(E)|\rho\rangle$, i.e. the first diagonal entry of the matrix $G(E)$, as

[TABLE]

where each term can be directly worked out:

[TABLE]

and

[TABLE]

It is clear that the propagator $G_{\rho}$ can be re-summed to all orders via

[TABLE]

where the self-energy of $\rho$ can be explicitly computed by

[TABLE]

A relation that will prove useful later is the energy-derivative of $\Sigma_{\rho}$:

[TABLE]

which is easily seen by noting

[TABLE]

III.2 The 2-pion propagator

Now we turn to the 2-pion states. The propagator $G_{2\pi}(E;q)$ can be explicitly worked out from the corresponding diagonal entries of $G(E)$:

[TABLE]

We now show that the diagonal T-matrix is directly related to the full propagator $G_{\rho}$. This is an important relation, as it dictates how the properties of the unstable state can be inferred from the scattering of the stable particles. To see that we employ the following expression of the T-matrix:

[TABLE]

The first term is [math] since $V$ is off-diagonal. The second term gives

[TABLE]

and hence

[TABLE]

This relates the amplitude of $\pi\pi$ scatterings to the $\rho$-propagator. Indeed, in the simple setting of the Lee model, all the physical information concerning the unstable state $\rho$ can be extracted from the diagonal T-matrix $\langle q|T|q\rangle$.

Finally, the full propagator $G_{2\pi}(E;q)$ can be obtained in closed form as

[TABLE]

III.3 Effective spectral function

We are now ready to examine the expressions of the phase shift $\mathcal{Q}$, the effective spectral function $B$ and the operator $\mathcal{K}$ in the context of the Lee model.

The phase shift can most simply be extracted from $G_{\rho}$ via 444One can obtain the phase shift directly from the S-matrix. See Eq. (47).

[TABLE]

The effective spectral function $B$ is known to be related to the interacting part of the density of state. In the Lee model, it is

[TABLE]

using the relation previously obtained

[TABLE]

we get

[TABLE]

Comparing with the expression of the KFL operator $\mathcal{K}$ in Eq. (20), we obtain

[TABLE]

Relation (39) constitutes the main result of this work. It demonstrates how the operator $B$ extracts the physical content, including the contribution from the unstable state $\rho$, of the system. From the first line of Eq. (39), we see that $B$ includes the contribution from the full spectral function $A_{\rho}$ and the 2-pion nonresonant interaction $\sum_{q}\Delta A_{2\pi}(E;q)$. The second line offers an alternative, but equivalent interpretation: $B$ includes the contribution from the bare-$\rho$, together with the interaction contribution contained in ${\rm tr}\,\mathcal{K}$. The latter includes contributions from the change in the energy spectra of both the $\rho$ and the 2-pion states.

To understand relation (39) better, we consider the interesting limit of vanishing coupling $g_{\rm eff}\rightarrow 0$. At this limit, the KFL operator $\mathcal{K}$ vanishes by definition. However the phase shift derivative operator $B$ would give

[TABLE]

that is, it becomes a Dirac-delta function for the bare $\rho$ state. It follows that the phase shift $\mathcal{Q}$ would becomes a step function $\mathcal{Q}\rightarrow\pi\times\theta(E-\Delta_{0})$. This is an intuitive limit for describing a $\rho$ that decouples from the pions: $\rho$ becomes a stable state, its width ceases to exist and the state should be included as particles in the asymptotic state. These are automatically implemented when the effective spectral function $B$ is used.

Another interesting limit is that of large bare resonant mass $\Delta_{0}>>E$. In this case, the resonant structure is suppressed, and the nonresonant term dominates. One can show that as $E\rightarrow 0$

[TABLE]

where $a_{l}$ is the scattering length of the channel and $l$ is the relative orbital angular momentum between the pions. Note that terms that are proportional to ${\rm Im}(\Sigma_{\rho})$ are subleading relative to $\frac{\partial}{\partial E}{\rm Im}(\Sigma_{\rho})$, as the latter is of $\mathcal{O}(q^{2l-1})$. This should be distinguished from the residual effect of the resonance width at threshold, which is of $\mathcal{O}(q^{2l+1})$. Thus, even an energy dependent Breit-Wigner model can not capture the effect of this term. In addition, the scattering lengths are well constrained by the chiral perturbation theory, and indeed the stated form of $B$ was derived Gerber and Leutwyler (1989).

III.4 Numerical results

As a numerical exercise, we solve the Lee Model on a momentum grid following the method of Ref. Hall et al. (2013). The Hamiltonians are constructed directly as an $(1+N_{q})\times(1+N_{q})$ matrix. The various Green functions are computed by matrix inversions. The method is very robust, and the procedure is as follows:

• (i)

Construct the matrices $H$ and $H_{0}$.

• (ii)

For each $E$, construct the matrices

[TABLE]

and invert

[TABLE]

• (iii)

Extract the quantities of interest. For example, the propagators are obtained from

[TABLE]

where $q_{i}$ is the discrete momentum of the $i$-th grid.

• (iv)

The spectral functions $A$’s can be obtained from the propagators by simply taking the imaginary part. To calculate the $B$ function, one possible method is

[TABLE]

In this exercise, we have chosen an appropriate P-wave coupling to describe the physical $\rho$-meson:

[TABLE]

with parameters $g=23.5\,{\rm GeV}^{-3/2}$, $\Delta_{0}=0.64$ GeV, and $\Lambda=0.4$ GeV. The form factor renders the real part of the self-energy finite. 555The imaginary part is finite even without a regulator. We checked that its value is only slightly modified by the form factor. It can be motivated in a quantum field theory with non-local interactions or dressing of vertex Terning (1991); Burdanov et al. (1996); Faessler et al. (2003); Giacosa et al. (2005). For what concerns us here it parametrizes the finite-size effects of the hadrons. See Ref. Herrmann et al. (1993) for a rigorous treatment within the QFT framework.

The numerical results for $\mathcal{Q}$ and various spectral functions are shown in Fig. 1. We have computed the phase shift on the momentum grid, ($1+N_{q}=800$; $L\approx 800$ fm; $\epsilon=0.01$ GeV), in two ways: one uses $G_{\rho}$ via Eq. (35) (grey squares), the other uses the scattering matrix $S$ via Eq. (47) (black circles). Both results agree quite well with the continuum limit, although they appear to have different convergence property. 666 Within this model one can also study how the phase shift and the spectral functions depend on the lattice size. This may help in understanding phase shift extraction in LQCD Hall et al. (2013); Briceno et al. (2018).

A key feature to note is the apparent shift of the strength of $B$, compared to $A_{\rho}$, towards lower energies. This effect originates from the nonresonant scattering term $\sum_{q}\Delta A_{2\pi}(E;q)$, and is needed, in addition to $A_{\rho}$, for a complete description of the interacting system.

III.5 Phase shift from S-matrix

In Eq. (35) we extracted the phase shift $\mathcal{Q}$ from the $\rho$-propagator. The same information is also available from the 2-pion sector via the S-matrix. In Ref. Lo (2017), the following recipe has been proposed to extract the phase shift from an $N_{q}\times N_{q}$ S-matrix

[TABLE]

where we have introduced the $\hat{t}$ operator, defined as

[TABLE]

From Eqs. (27) and (33) we see

[TABLE]

and indeed it is straightforward to verify the same result of the phase shift as in Eq. (35):

[TABLE]

The approximation in the last line of Eq. (47) was shown to be valid for some simple cases such as s-channel-only interaction or structureless scattering. We shall now show that the approximation is exact in the Lee model.

Consider the expansion of the logarithm of the S-matrix via the Mercator series

[TABLE]

The approximation in Eq. (47) becomes exact if $\hat{t}$ satisfies the following property

[TABLE]

giving

[TABLE]

Note the effective replacement of $I$ with $1$ in the above.

Inspecting the interaction term of the Lee Model in Eq. (12), we see that the matrix element of $\hat{t}$ takes the product form

[TABLE]

It follows that

[TABLE]

and hence ${\rm tr}\,({\hat{t}}^{2})=({\rm tr}\,\hat{t})^{2}$ and so on for higher power. This demonstrates that criteria (52) is satisfied by a class of $\hat{t}$ which involves separable potentials, i.e. the ones that are from a Kronecker (direct) product.

Furthermore, we add that from Cayley Hamilton theorem, the determinant of $\hat{t}$ essentially vanishes. In fact, all eigenvalues of $\hat{t}$ are zero, except one, which is ${\rm tr}\,\hat{t}$. This is another way to understand why the approximation made in Eq. (47) is justified. The full implication of this result is not yet completely clear, and will be explored in a future work.

IV Thermodynamics

The change in the density of state due to interactions, as revealed by the KFL $\mathcal{K}$ operator or the $B$ function, is the key input for the S-matrix formulation of statistical mechanics Dashen et al. (1969); Venugopalan and Prakash (1992). The approach is based on the method of cluster expansions, and for the second virial coefficient the result is exact. We retrace a few basic steps in relating the scattering phase shift to the thermal partition function.

Our starting point is the cluster expansion of the grand partition function

[TABLE]

where $\xi$ is the fugacity, related to the particle chemical potential via $\xi=e^{\mu/T}$. $Z_{N}$ is the N-body partition function. The corresponding expansion for the logarithm of $Z$ reads

[TABLE]

and one can work out

[TABLE]

The interacting part of the partition function satisfies 777For simplicity we neglect corrections due to quantum statistics and focus on the $\mu=0$ case.

[TABLE]

The latter can be re-expressed via

[TABLE]

where we have integrated out the CM motion in the total energy of the 2-body system

[TABLE]

with $m_{\rm tot}$ being the total mass. $E$ is the energy of the relative motion

[TABLE]

where $m_{\rm red}$ is the reduced mass. We have made use of the form of the nonrelativistic dispersion (see e.g., Ref. Lo (2017) for a relativistic formulation) to integrate out the CM momentum $\vec{P}$ and obtain the thermal wavelength $\lambda_{T}$

[TABLE]

Note that the bare $\rho$ state is not counted in the trace of $Z_{2}^{(0)}$. The remaining integral in Eq. (60) requires the input of $D(E)$, which is the change in the density of state due to interactions.

Now we are ready to examine the thermal contribution of an unstable state $\rho$ based on the input from the Lee model. The first thing to notice is that the L.H.S. of Eq. (60) can be directly computed from the eigenvalues ($\lambda_{n}$) of the Hamiltonian: 888In fact, computations involving the difference between the fully interacting system and the free case, e.g., $G(z)-G_{0}(z)$, are numerically more stable Balian and Bloch (1970). This also opens up the possibility of solving the system with other numerical techniques, such as the use of a harmonic oscillator basis in the expansion Kruppa (1998).

[TABLE]

The result is shown as points in Fig. 2. This in turns indicates what is the correct $D(E)$ to use. The S-matrix formulation of statistical mechanics by Dashen et al. Dashen et al. (1969) dictates the choice of $D(E)\rightarrow B(E)=2\frac{\partial}{\partial E}\mathcal{Q}(E)$. This means that we compute the thermodynamic pressure via

[TABLE]

$\Delta P_{B}$ includes the contribution of the unstable state $|\rho\rangle$ and nonresonant $\pi\pi$ interaction:

[TABLE]

From Eq. (39) and the discussion it is clear that one can incorporate the same physical content of the thermal medium with a different choice of $D(E)$, though with a different interpretation. For example, one can choose instead $D(E)\rightarrow{\rm tr}\mathcal{K}$, and in this case the contribution of $A^{0}_{\rho}$ needs to be added separately as

[TABLE]

Here $\Delta P_{\mathcal{K}}$ contains the contribution from $\Delta A_{\rho}$ and $\Delta A_{2\pi}$. Note that as $g\rightarrow 0$, $\Delta P_{\mathcal{K}}\rightarrow 0$; while $\Delta P_{B}\rightarrow P^{(0)}_{\rho}$. The two results in Eqs. (65) and (67) are equivalent, i.e.

[TABLE]

which is just restating the relation (39).

The various partial pressures are shown in Fig. 2. Due to the contribution from the nonresonant $\pi\pi$ state, the pressure based on $B$ is substantially larger than the one based on $A_{\rho}$ alone. We stress that only the former one gives a consistent description of the thermodynamics, as can be verified by the direct construction of the partition function from the eigenvalues of the Hamiltonian (black circles). Therefore Eq. (5) should be used instead of Eq. (4).

These observations are in accord with the previous analysis based on a $\Phi$-derivable approach Weinhold et al. (1998). The $B$-function-based description requires only the input from the scattering of asymptotic states. This underlines an important concept in the formulation: In computing the density of states it is not mandatory to introduce the unstable state as an explicit degree of freedom. For approaches that use only stable states as degrees of freedom, such as an effective field theory where resonances are dynamically generated Kaplan et al. (1996), the same density of states would be obtained as long as the phase shifts agree. And when an empirical phase shift $\mathcal{Q}(E)$ is used, the function $B(E)$ becomes model independent, while the splitting into $A_{\rho}(E)$ and $\sum_{q}\Delta A_{2\pi}(E;q)$ is model dependent. The Lee model studied here provides a clear picture of such a splitting, and demonstrates how an unstable state should be included in the description of the thermodynamics.

V Conclusion

In the context of the Lee model we have clarified the relation between the energy derivative of the phase shift and the spectral functions of the degrees of freedom composing the system. We have also illustrated how these quantities enter the thermal description of the system via the S-matrix formulation of the statistical mechanics. This consolidates our understanding of the connection between this and the standard approach based on thermal Green functions. In particular, we have shown that the thermodynamic trace requires the inclusion of the nonresonant contribution ($\Delta A_{2\pi}$), in addition to, and independent of, the effect coming from the width of the unstable state $\Delta A_{\rho}$.

Besides acting as an effective density of state, an alternative interpretation of the energy derivative of the phase shift is the concept of time delay Danielewicz and Pratt (1996); Kelkar and Nowakowski (2008): particles spend longer or shorter in the interaction region due to the attractive or repulsive nature of the interaction. In the contexts of transport models and resonance identification, it was argued Bass et al. (1998); Leupold (2001, 2003); Kelkar et al. (2003, 2004) that such a time delay, instead of the inverse width $1/\Gamma(E)$, should be used to measure the life-time of a resonance. A related problem is the study of the survival probability of an unstable state. According to Ersak Ersak ; Fonda et al. (1978), the standard exponential decay law is valid only in the limited case of an energy-independent Breit-Wigner distribution. Re-scattering effect, apparently related to $\Delta A_{2\pi}$, would lead to non-exponential behavior Bohm and Sato (2005). A clearer theoretical understanding of $B$ and $A_{\rho}$ could provide further insights into these topics.

So far we have restricted our discussion to Fock space up to two body. It will be extremely interesting to extend the scheme to include multi-channel and multi-body scatterings Chiu et al. (1992); Kaminski et al. (1997); Mai and Dring (2017); Lo (2017); Fernandez-Ramirez et al. (2018), and understand how these interactions would influence thermodynamic quantities. This can be a useful framework to analyze the observables in Heavy Ion Collision experiments, such as hadron yields and the momentum distributions of light hadrons Abelev et al. (2013); Sollfrank et al. (1990); Broniowski et al. (2003); Lo (2018). We defer this more challenging problem to future research.

acknowledgments

PML thanks Eric Swanson for stimulating discussions. He also acknowledges fruitful discussions with Hans Feldmeier, Bengt Friman, and Piotr Bozek. The authors are also grateful for the constructive conversations with Wojciech Broniowski, Wojciech Florkowski, and Stanislaw Mrowczynski. PML was partly supported by the Polish National Science Center (NCN), under Maestro Grant No. DEC-2013/10/A/ST2/00106 and by the Short Term Scientific Mission (STSM) program under COST Action CA15213 (reference number: 41977). FG acknowledges financial support from the Polish National Science Centre (NCN) through the OPUS project no. 2015/17/B/ST2/01625.