Thermal contribution of unstable states

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.


Introduction
Formal treatment of the interactions in a gas of particles at finite temperature is an important topic in thermal field theory [1][2][3]. In particular, a consistent description of the unstable states is imperative for understanding the hadron gas, see, e.g., Refs. [3][4][5][6]. 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? [5,[7][8][9][10][11][12][13] 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 [5,14,15], according to which the Electronic supplementary material The online version of this article (https://doi.org/10.1140/epjc/s10052-019-6844-y) contains supplementary material, which is available to authorized users. a e-mail: pmlo@gsi.de 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 [16,17]. 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 σ -meson, is involved. The empirical phase shift, analyzed by the chiral perturbation theory [18], 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 [19].
In this study we consider a system of stable particles "π " (the pions) and unstable particle "ρ" (the ρ-mesons). 1 Each ρ can decay into two π 's via the interaction ρ → ππ. Chains of interactions, e.g., ππ → ρ → ππ, 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) [20,21] to describe the interactions in the ρπ system. Similar Hamiltonians have been used to explore various areas in physics, ranging from atomic physics and quantum optics [22][23][24][25] to baryon decays [26]. The Lee model offers a useful theoretical setup to study an interacting system [27][28][29][30][31][32][33] and in many aspects it resembles a quantum field theory [34]. In practice, the LH contains an unstable state |ρ , as well as a continuum of 2-body states |π(q)π(−q) with all possible relative momenta q of the pair. The decay ρ → ππ is described by mixing terms of the form: ∼ g (|ρ ππ| + |ππ ρ|). This gives the unstable particle ρ a width, i.e., a distribution in energies (or masses) dictated by the spectral function A ρ (E) [34][35][36][37] can be interpreted as the probability that the unstable ρ has energy between E and E + d E. The spectral function can be calculated from the imaginary part of the ρ-propagator. See Sect. 3 for details. For a narrow-width state this can be approximated by the Breit-Wigner (BW) formula [38,39], Note that the width ρππ is generally energy dependent, and the mass m ρ 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. [27,40] and the experimental results in Refs. [41][42][43]; and Ref. [44] 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 Sect. 4 for details. Here we give a synopsis of our discussion. Consider the hypothetical limit where ρ and π 's are completely decoupled, i.e., g = 0 (or ρππ = 0). The pressure of the system would be given by a sum of two contributions: 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 ρ. This is the fundamental premise of the hadron resonance gas (HRG) model [45,46]: 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 ρ via a weighed sum by A ρ (E): such that the total pressure is approximated as (Scheme-A): This scheme is employed in many versions of the HRG models, see e.g., Refs. [3,4,[47][48][49][50]. See also the K-matrix-based approach [51,52]. However, Eq. (4) is not yet complete. According to the S-matrix formulation of statistical mechanics by Dashen et al. [5] (see also the discussion by Weinhold et al. [15]), the correct result of the pressure at arbitrary g (or ρππ ) is given by (S-matrix scheme): where Q ππ (E) is the phase shift for the scattering process ππ → ππ. Here we summarize some key features of the S-matrix scheme: (i) Observe that there is no explicit ρ 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 ρ 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) Equations (4) and (5) reduce to the free gas result (2) in the limit of g → 0 ( or ρππ → 0 ). 2 For the latter, we have (iii) Generally, B(E) = A ρ (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 κ(700) [16,17], the N * and Δ resonances [53,54], the S = −1 hyperons [55], etc. It is also the case for the recently discovered X, Y, Z states [56][57][58]. As shown in a recent work [59], 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. [15], but the actual adoption of the scheme remains limited [54,60]. 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 where the second term in the R.H.S. describes the the modification of the spectral function of the 2-body |ππ state. Such a term is present even in the absence of the resonance, e.g., taking the large m ρ limit (see Sect. 3). The correct expression of the pressure can be decomposed as: where 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 Sect. 2 the details of the Lee model are presented. Then, in Sect. 3, the spectral functions for both ρ and ππ, and the phase shift are introduced. Here the important Eq. (7) is derived. In Sect. 4 the thermodynamic properties of the system are determined analytically, with special focus on the pressure with its various contributions. A numerical example shows that ΔP 2π can be sizable and in general should not be neglected. Finally, discussions and conclusions are given in Sect. 5.

The Lee model
The Lee model [20] describing the ρ ↔ ππ system can be formulated as follows. 3 Introducing the basis states in the center-of-mass (CM) frame: where q is the momentum label for the two-pion state The Hamiltonian of the system can be represented as an 3 In the following we measure energy with respect to 2m π , and the nonrelativistic dispersion relation (q) = q 2 /m π 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. [61]. It is easy to go to the continuum by taking The interaction V describes the coupling of ρ with the |q states such that We use to implement the spherical degeneracy. Note that δq = 2π L in the finite volume formulation, L being the size of the box. The coupling g eff is generally q-dependent. The full Hamiltonian then reads With the Hamiltonian defined, we can construct the resolvent operators which can be understood again as (1 + N q ) × (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: and 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 [13,62] Since ρ is not an asymptotic state, the actual scattering matrix S is extracted from the (lower-right) N q × N q block ofŜ. In addition, we introduce an operator K , due to Krein, Friedal and Lloyd (KFL) [63,64], defined as the difference of the spectral functions where we have identified the spectral function operator These operators are deeply connected with the scattering phase shift Q and the effective spectral function B. The latter is defined as There are multiple ways to extract the phase shift Q from the resolvents (or Green functions). (See Eqs. (36), (48)). These and the explicit relations among the various observables will be demonstrated in the context of the Lee model in Sect. 3.

Phase shift and effective spectral function
In the Lee model various theoretical quantities, e.g. the propagator and the self-energy of ρ, 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 K and the derivative of the phase shift.

The ρ-propagator
The ρ-propagator G ρ can be computed from ρ|G(E)|ρ , i.e. the first diagonal entry of the matrix G(E), as where each term can be directly worked out: and It is clear that the propagator G ρ can be re-summed to all orders via where the self-energy of ρ can be explicitly computed by A relation that will prove useful later is the energyderivative of ρ : (29) which is easily seen by noting

The 2-pion propagator
Now we turn to the 2-pion states. The propagator G 2π (E; q) can be explicitly worked out from the corresponding diagonal entries of G(E): We now show that the diagonal T-matrix is directly related to the full propagator G ρ . 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: The first term is 0 since V is off-diagonal. The second term gives and hence This relates the amplitude of ππ scatterings to the ρpropagator. Indeed, in the simple setting of the Lee model, all the physical information concerning the unstable state ρ can be extracted from the diagonal T-matrix q|T |q . Finally, the full propagator G 2π (E; q) can be obtained in closed form as

Effective spectral function
We are now ready to examine the expressions of the phase shift Q, the effective spectral function B and the operator K in the context of the Lee model. The phase shift can most simply be extracted from G ρ via 4 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 4 One can obtain the phase shift directly from the S-matrix. See Eq. (48).
using the relation previously obtained we get Comparing with the expression of the KFL operator K in Eq. (21), we obtain Relation (40) constitutes the main result of this work. It demonstrates how the operator B extracts the physical content, including the contribution from the unstable state ρ, of the system. From the first line of Eq. (40), we see that B includes the contribution from the full spectral function A ρ and the 2-pion nonresonant interaction q ΔA 2π (E; q). The second line offers an alternative, but equivalent interpretation: B includes the contribution from the bare-ρ, together with the interaction contribution contained in tr K . The latter includes contributions from the change in the energy spectra of both the ρ and the 2-pion states.
To understand relation (40) better, we consider the interesting limit of vanishing coupling g eff → 0. At this limit, the KFL operator K vanishes by definition. However the phase shift derivative operator B would give that is, it becomes a Dirac-delta function for the bare ρ state. It follows that the phase shift Q would becomes a step function Q → π ×θ(E −Δ 0 ). This is an intuitive limit for describing a ρ that decouples from the pions: ρ 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 Δ 0 >> E. In this case, the resonant structure is suppressed, and the nonresonant term dominates. One can show that as 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 Im( ρ ) are subleading relative to ∂ ∂ E Im( ρ ), as the latter is of O(q 2l−1 ). This should be distinguished from the residual effect of the resonance width at threshold, which is of 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 [65].

Numerical results
As a numerical exercise, we solve the Lee Model on a momentum grid following the method of Ref. [61]. The Hamiltonians are constructed directly as an (1 + N q ) × (1 + N q ) matrix. The various Green functions are computed by matrix inversions. The method is very robust, and the procedure is as follows: and invert (iii) Extract the quantities of interest. For example, the propagators are obtained from 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 1). (46) In this exercise, we have chosen an appropriate P-wave coupling to describe the physical ρ-meson: with parameters g = 23.5 GeV −3/2 , Δ 0 = 0.64 GeV, and = 0.4 GeV. The form factor renders the real part of the self-energy finite. 5 It can be motivated in a quantum field theory with non-local interactions or dressing of vertex [66][67][68][69]. For what concerns us here it parametrizes the finite-size effects of the hadrons. See Ref. [70] for a rigorous treatment within the QFT framework.
The numerical results for 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 ≈ 800 fm; ε = 0.01 GeV), in two ways: one uses G ρ via Eq. (36) (grey squares), the other uses the scattering matrix S via Eq. (48) (black circles). Both results agree quite well with the continuum limit, although they appear to have different convergence property. 6 A key feature to note is the apparent shift of the strength of B, compared to A ρ , towards lower energies. This effect originates from the nonresonant scattering term q ΔA 2π (E; q), and is needed, in addition to A ρ , for a complete description of the interacting system.

Phase shift from S-matrix
In Eq. (36) we extracted the phase shift Q from the ρpropagator. The same information is also available from the 2-pion sector via the S-matrix. In Ref. [13], the following recipe has been proposed to extract the phase shift from an N q × N q S-matrix where we have introduced thet operator, defined aŝ  28) and (34) we see and indeed it is straightforward to verify the same result of the phase shift as in Eq. (36): The approximation in the last line of Eq. (48) 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 tr ln S = tr ln(I + i ×t) The approximation in Eq. (48) becomes exact ift satisfies the following property Note the effective replacement of I with 1 in the above. Inspecting the interaction term of the Lee Model in Eq. (13), we see that the matrix element oft takes the product form It follows that and hence tr (t 2 ) = (trt) 2 and so on for higher power. This demonstrates that criteria (53) is satisfied by a class oft 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 oft essentially vanishes. In fact, all eigenvalues oft are zero, except one, which is trt. This is another way to understand why the approximation made in Eq. (48) is justified. The full implication of this result is not yet completely clear, and will be explored in a future work.

Thermodynamics
The change in the density of state due to interactions, as revealed by the KFL K operator or the B function, is the key input for the S-matrix formulation of statistical mechanics [5,6]. 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 where ξ is the fugacity, related to the particle chemical potential via ξ = e μ/T . Z N is the N-body partition function. The corresponding expansion for the logarithm of Z reads and one can work out The interacting part of the partition function satisfies 7 The latter can be re-expressed via where we have integrated out the CM motion in the total energy of the 2-body system 7 For simplicity we neglect corrections due to quantum statistics and focus on the μ = 0 case. with m tot being the total mass. E is the energy of the relative motion where m red is the reduced mass. We have made use of the form of the nonrelativistic dispersion (see e.g., Ref. [13] for a relativistic formulation) to integrate out the CM momentum P and obtain the thermal wavelength λ T Note that the bare ρ state is not counted in the trace of Z (0) 2 . The remaining integral in Eq. (61) 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 ρ based on the input from the Lee model. The first thing to notice is that the L.H.S. of Eq. (61) can be directly computed from the eigenvalues (λ n ) of the Hamiltonian: 8 Tr 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. [5] dictates the choice of D(E) → B(E) = 2 ∂ ∂ E Q(E). This means that we compute the thermodynamic pressure via ΔP B includes the contribution of the unstable state |ρ and nonresonant ππ interaction: From Eq. (40) 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) → trK , and in this case the contribution of A 0 ρ needs to be added separately as Here ΔP K contains the contribution from ΔA ρ and ΔA 2π . Note that as g → 0, ΔP K → 0; while ΔP B → P ρ . The two results in Eqs. (66) and (68) are equivalent, i.e.
which is just restating the relation (40). The various partial pressures are shown in Fig. 2. Due to the contribution from the nonresonant ππ state, the pressure based on B is substantially larger than the one based on A ρ 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 -derivable approach [15]. The B-functionbased 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 [74], the same density of states would be obtained as long as the phase shifts agree. And when an empirical phase shift Q(E) is used, the function B(E) becomes model independent, while the splitting into A ρ (E) and q ΔA 2π (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.

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 (ΔA 2π ), in addition to, and independent of, the effect coming from the width of the unstable state ΔA ρ .
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 [75,76]: 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 [77,[79][80][81][82] that such a time delay, instead of the inverse width 1/ (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 Fonda et al. [40], the standard exponential decay law is valid only in the limited case of an energy-independent Breit-Wigner distribution. Rescattering effect, apparently related to ΔA 2π , would lead to non-exponential behavior [83]. A clearer theoretical understanding of B and A ρ 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 [13,21,55,84,85], 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 [60,[86][87][88]. We defer this more challenging problem to future research. ence number: 41977). FG acknowledges financial support from the Polish National Science Centre (NCN) through the OPUS project no. 2015/17/B/ST2/01625.

Data Availability Statement
This manuscript has associated data in a data repository. [Authors' comment: The phase shift data analyzed are extracted from Refs. [89,90]. Results of our model calculation are available in the supplementary datafile.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .