\Delta self-energy at finite temperature and density and the \pi N cross-section

self-energy at finite temperature and density and the cross-section

Snigdha Ghosh Variable Energy Cyclotron Centre, 1/AF Bidhannagar, Kolkata - 700064, INDIA Homi Bhabha National Institute, Training School Complex, Anushaktinagar, Mumbai - 400085, INDIA    Sourav Sarkar Variable Energy Cyclotron Centre, 1/AF Bidhannagar, Kolkata - 700064, INDIA Homi Bhabha National Institute, Training School Complex, Anushaktinagar, Mumbai - 400085, INDIA    Sukanya Mitra Indian Institute of Technology Gandhinagar, Gandhinagar-382355, Gujarat, INDIA

The self energy of -baryon is evaluated at finite temperature and density using the real time formalism of thermal field theory. The Dyson-Schwinger equation is used to get the exact thermal propagator followed by the spectral function of . The scattering cross section obtained using explicit exchange is normalized to the experimental data in vacuum and its medium modification is implemented by means of the exact thermal propagator. A significant suppression of the peak of the cross-section is observed at higher temperature and baryon density. Effects on the mean relaxation time of nucleons and the temperature dependence of the shear viscosity of a pion nucleon gas are demonstrated.

I Introduction

The properties of hadrons under conditions of high temperature and/or baryon density have been widely studied owing to the possibility of restoration of chiral symmetry of QCD spontaneously broken by the vacuum Bazavov (). Relativistic collisions of heavy ions provide the unique opportunity to actually produce hot/dense matter under controlled conditions and provide experimental verification of such studies (see e.g. Rapp ()). The spectral changes of hadrons as a result of compression and/or heating thus play a very significant role in unraveling the dynamics of the underlying colour degrees of freedom and the vacuum structure of QCD CBMBook (). The broadening of vector spectral function in the medium observed through the invariant mass spectra of dileptons in heavy ion collisions at SPS Arnaldi () and RHIC Adare () have been interpreted (see e.g. Sarkar ()) to have non-trivial, though indirect implications on chiral symmetry restoration. In this connection the study of spectral properties of baryon resonances such as and etc. are important; though not measurable directly because of final state interactions, they provide valuable input to the evaluation of the self-energies of low mass vector mesons which can be measured in heavy ion collision experiments through electromagnetic probes. In addition to broad resonances stable hadrons such as the nucleon also develop considerable widths in the medium Ghosh2 () and this is mainly due to resonant scattering with pions involving baryon resonances.

The in-medium self-energy of the has played a key role in the understanding of dynamics of nuclear interactions particularly in the resonance region where it is excited as a resonance Ericson (). The propagation of , and thus become intimately connected. The description of pion-nucleus interaction for example, depends on the formation, propagation and decay of the in the nuclear environment Oset (); Malfliet (). Models based on this picture, known as -hole models, successfully explain many scattering phenomena in elastic as well as inelastic channels. The self-energy extracted using this approach has also found applications in transport models describing the dynamical evolution of nucleus-nucleus collisions Mosel (); Wolf ().

The spectral modification of the at finite temperature and density and the in-medium cross-section remains an interesting issue of discussion. A very useful and topical application of this study is to investigate the medium dependence of the shear viscosity of hot and dense hadronic matter which is produced towards the later stages of relativistic heavy ion collisions. In the kinetic theory approach of evaluating transport coefficients like shear and bulk viscosities, thermal conductivity etc. the cross-section enters as the dynamical input. Since pions account for most of the multiplicity in relativistic heavy ion collisions, the case of the pion gas has received some attention in recent times. Scattering amplitudes evaluated using chiral perturbation theory to lowest order have been used in Santalla (); Chen () and unitarization improved estimates were employed in Dobado () to evaluate the shear viscosity. Again, phenomenological scattering cross-section using experimental phase shifts were utilized in Prakash (); Chen (); Davesne (); Itakura () to obtain the viscous coefficients. While in Moore (); Juan () the effect of number changing processes on the bulk viscosity of a pion gas has been studied, in Fraile () unitarized chiral perturbation theory was used to demonstrate the breaking of conformal symmetry by the pion mass. It is important to point out that in all these approaches vacuum amplitudes have been used. However, cross-sections of scattering between excitations in the medium could be non-trivially affected by the presence of other particles constituting it. It follows that in addition to the usual quantum fluctuations, thermal fluctuations modify the propagation of mediating particles and consequently the invariant amplitude. The case of the cross-section in a hot pion gas was studied by Barz et al Barz () obtaining a substantial reduction of the magnitude in the region of the peak. With this motivating feature, recently Mitra1 (), the amplitudes for elastic scattering were evaluated using in-medium propagators for the exchanged and mesons incorporating loop graphs with , , and mesons in the internal lines Ghosh1 (). The cross-section so obtained showed a substantial reduction in magnitude at the peak position at higher values of temperature. When used in the collision integral of the transport equation, these medium modified scattering amplitudes were found to significantly affect the temperature dependence of the viscosities Mitra2 () and thermal conductivity Mitra3 () of a hot pion gas.

Based on our experience with the pion gas and keeping in view the upcoming CBM experiment at FAIR it is natural to ask how the presence of a finite baryon density in addition to temperature is likely to affect the shear viscosity. To study such effects one has to include nucleons and consequently the cross-section becomes the principal dynamical factor in its evaluation. Incorporating the in-medium cross-section calculated using the modified self-energy a more reliable estimate of the shear viscosity, in particular it’s dependence on temperature and baryon density can be obtained. When used as input in the viscous hydrodynamic equations a more realistic scenario of space time evolution of the later stages of heavy ion collisions is likely to be achieved.

In this work we obtain the self-energy at finite temperature and baryon density evaluating several one-loop diagrams with , , and in the internal lines using standard thermal field theoretic methods. The in-medium propagator of the is then used in the scattering amplitudes to obtain the cross-section. This is utilized to evaluate the relaxation times in a hadronic gas mixture of pions and nucleons. Finally, the temperature and density dependence of the shear viscosity is obtained.

Medium modifications of the resonance has been studied mostly in nuclear matter Oset (); Wehrberger (); Korpa1 (); Kim (); Xia (). Whereas a many body expansion in terms of particle-hole excitations has been used in Oset () to evaluate the self-energy in nuclear matter, in Wehrberger () and Kim () its decay in the medium is investigated using the quantum hadrodynamics (QHD) model. A self-consistent treatment of pions and ’s in nuclear matter at zero and finite temperature may be found in Xia () and Korpa2 () respectively. In Korpa1 (), the self-energy due to loop in nuclear matter is obtained in a relativistic approach. More recently, modification of the spectral function at finite temperature and density due to resonant scattering off thermal pions has been obtained in vanHees (). In addition to these theoretical studies, properties of the have been studied experimentally using invariant mass analyses of pairs Hjort (); Mishra ().

In the next section we recall some basic features of the real-time version of thermal field theory. Then we evaluate the self-energy and discuss numerical results. In the subsequent section we evaluate the amplitudes of scattering leading to the cross-section in the medium. This is followed by a section on the shear viscosity of a gas and finally by summary and discussions. Some details of the calculation is provided in the Appendix.

Ii The self-energy in the medium

ii.1 The in-medium propagators in the real time formalism

In the real time formalism of thermal field theory, all two-point functions including the self-energy take the form of matrices Mallik (); Bellac (). But each of the matrices may be diagonalized, when it is given essentially by a single analytic function which determines completely the dynamics of the corresponding two-point function. This function being given by any one, say the 11-component of the matrix, we need to evaluate only this component of the self-energy matrix. In the following we specify only the 11-component of the thermal propagator for the particles involved in the one-loop graphs that are considered in this work.

The 11-component of a free thermal propagator matrix for a particle consists of its vacuum propagator and a term depending on the on-shell distribution function of like particles in the medium through which it propagates. The form of the latter term depends only on whether the particle in question is a boson or a fermion.

The 11-component of the thermal pion propagator is given by





and for and 0 for .

The thermal part remaining same, the propagator with the familiar polarization sum follows as


We now consider the fermionic propagators whose thermal (matrix) parts are different from the bosonic ones discussed above. The 11-components of the nucleon and propagators are respectively given by




where, is given in Appendix-A.

The complete propagator is given by the Dyson equation in terms of the free fermion propagator and self-energy ,


where each is a 22 matrix in the thermal indices. They can be diagonalized to get the respective analytic functions, denoted by bar, so that


The self-energy function can be obtained from any single component of the self energy matrix as discussed in Appendix-A. It is related to, say, the 11-component by


where for and for .

ii.2 The self-energy

Figure 1: Feynman Diagrams for -Self Energy

We begin by writing down the expressions for the one-loop self-energy graphs for the in vacuum. For the four cases shown in fig. 1 they can be generally expressed as


where contains terms coming from the two vertices and the spin factors appearing in the propagators for the internal lines. They can be read off from the expressions for given in Appendix-B and appear as




We now proceed to write down the corresponding expressions in the medium. As discussed above, we need to evaluate only the 11-component which is obtained by replacing the vacuum propagators by the 11-component of thermal propagators given in the previous section. The self-energy in the medium is thus given by


Expanding and we obtain in addition to , two terms which are linear in the thermal distribution function and the fourth term non-linear in the distribution function which is purely imaginary. Performing the integral and using eqns. (8) we obtain the imaginary and real parts of the self-energy function as


where and .

Each of the terms in the imaginary part can be related to scattering and decay of the -baryon. The delta functions in the four terms define the kinematic domains where these processes occur. The regions where these are non-vanishing correspond to branch cuts in the complex plane. The first and third terms in (LABEL:im1) are non-vanishing for . This is the usual unitary cut already present in vacuum. The second and fourth terms which are non-zero for defines the Landau cut and is purely a medium effect. Confining ourselves to the kinematic region and , the first and fourth terms only contribute. The first corresponds to absorption of the due to decay into a baryon-meson pair such as , etc. and is thus weighted by a thermal factor , indicating Bose enhancement of the meson and Pauli blocking of the baryon in the process plus the usual thermal factors for the initial state in the time reversed process where the is produced. The fourth term is due to absorption of the due to scattering from a meson producing a baryon in the final state and vice versa as is evident from the thermal weight factor .

To take into account the finite width of unstable particles in the loop graphs the self-energy is folded with their spectral functions. As a consequence the sharp thresholds of the branch cuts get smeared. For unstable mesons () we use Gonzalez (),


with and . Here so that


where . This is obtained using with .

For baryons () with non-trivial decay width in the loops we use


with and . In this case for which the decay formula is given by eq. (22) in Section III.

ii.3 Numerical results

Let us begin with the results for the self-energy in the medium. We show numerical results for the spin-averaged real and imaginary parts of the function given by Wehrberger (); Kim ()


where denote Rarita-Schwinger spinors. The factors given in eqs. (13) then go over to . In fig. 2 we plot the Landau cut contribution to the imaginary part coming from the different loop graphs at MeV for MeV and MeV in panels (a) and (b) respectively. The corresponding results for MeV are shown in panels (c) and (d). As discussed in the last section, it follows that for the Landau cut extends up to the difference of masses of the baryon and the meson in the loop graph in the stable limit. These contributions at lower values of are a result of scattering processes in the thermal medium and are absent in vacuum. Comparing fig. 2(a) with (b) (and (c) with (d)) we see that the and loops start contributing to the imaginary part only at larger baryon densities.

Figure 2: The Landau-cut contribution to the imaginary part of the self energy function for different and

The unitary cut contributions to the imaginary part are much larger in magnitude than the ones coming from the Landau cut. The thresholds lie at higher energies, for stable particles starting from the sum of the nominal masses of the particles in the loop graph. Shown in fig. 3 are the unitary cut contributions from the four loops at MeV. We observe the sequential opening up of heavier decay channels. These are the same as in vacuum but are now weighted by the Pauli blocking and Bose enhancement factors in the final state. The form factor suppresses the usual monotonous rise of these contributions at higher .

Figure 3: The unitary cut contribution to the imaginary part of the self-energy at =100 MeV.

We now consider the thermal component of the real part of the self-energy consisting of principal value integrals. As seen in fig. 4 the magnitudes are quite small compared to the imaginary parts and are not expected to contribute to a thermal shift in the pole position of the .

Figure 4: The real part of the self energy function for different and

Having obtained the imaginary and real parts of the self-energy we now plot the spectral function which is the imaginary part of the complete propagator (7). We do so in two parts. The low region depicting the Landau cut contribution shown in panel (a) of fig. 5 is purely a thermal contribution. The high region in the vicinity of the bare mass consisting of the contributions from the unitary cuts is shown in panel (b) of fig. 5. The spectral density is seen to have a significant dependence on the temperature and chemical potential of the medium. However, comparing the curves at and MeV, both for MeV, the effect of baryon density only shows up at higher temperatures and has its origin in the thermal distribution functions for the baryons. In general, we find a gradual suppression of the peak with increasing temperature and density owing to the larger imaginary parts in the denominator of the in-medium propagator. As seen from eq. (LABEL:im1), the increase in the imaginary part comes from two factors: (i) the Bose enhancement factor for the pions and rho mesons in the final state in the first term which is the unitary cut contribution and (ii) the Landau cut contribution coming from scattering of the mesons off the propagating as given by the fourth term. The second and third terms do not contribute because of kinematic reasons.

Figure 5: The spectral function of the for (a) up to 1 GeV and (b) from 1 to 1.8 GeV

In fig. 6(a) we have made a comparison of the spectral function obtained in our approach with that of vanHees () for two sets of values of temperature and nucleon density representative of conditions likely to be achieved in relativistic heavy ion collisions at RHIC and in the CBM experiment. A reasonable agreement is observed between the present work and that of vanHees () as depicted by the continuous lines and symbols respectively. In addition to differences in the Lagrangian and associated parameters used in the two approaches the disparities in the spectral functions in the two cases could arise due to contributions coming from higher order effects introduced through dressed nucleon and pion propagators in the self-energy considered in vanHees () wherein vertex corrections were included through Migdal parameters in the pion propagator.

Figure 6: The spectral function as a function of : (a) in comparison with vanHees () and (b) for different form factors. The symbols correspond to the results of Ref. vanHees ()

To take into account the finite size of the vertices a phenomenological hadronic form factor has been introduced. The details are provided in Appendix-B. The numerical results presented here correspond to a monopole type form factor, denoted by Form Factor I in fig. 6(b) with MeV which produces a good fit to the phase shift data and cross section. We also plot the spectral function using Form Factor I with MeV and find a small reduction at the peak though the cross section remains largely unchanged. This is because in the spectral function the square of the form factor appears multiplicatively but in the cross-section its effect is largely canceled at the -channel pole position as seen from eq. (23) in the next section. For a comparison we also plot the spectral function using an exponential form factor denoted by Form Factor II with GeV (as used in Korpa1 () with GeV). No appreciable difference is found with the one with 600 MeV. All the plots in this figure correspond to MeV, MeV and MeV. The symbols denote the results of vanHees ().

At this point a few comments on the hadronic form factors at the vertices are in order. It is well known that in local field theory involving spin-3/2 fields the redundant degrees of freedom associated with unphysical spin-1/2 fields are eliminated through the Rarita-Schwinger constraints Rarita (). The interacting case is more complex and suffers due to participation of the spurious spin-1/2 components. A coupling consistent with gauge invariance of the Rarita-Schwinger field was constructed in Pascalutsa () preserving the correct number of degrees of freedom. However, unphysical behavior in the computed tree-level cross-section results if the reaction is cut off by standard hadronic form factors Vrancx (). In this work we have used the conventional vertex, which does suffer from a small presence of spin-1/2 components both in vacuum and in nuclear matter at saturation density Korpa1 (). As described above, the cut-off in the form factor used here was obtained by fitting the cross-section. The spectral functions evaluated using this (conventional) vertex were found Korpa3 () to differ slightly with the ones calculated using the consistent coupling discussed above if the same form factor is applied in both the cases. It was further shown that this difference could be eliminated if the additional momentum factor stemming from the higher derivative nature of the consistent interaction Pascalutsa () was compensated either by an additional form factor term or by adjusting the cut-off values of the original form factor.

Iii The -N Cross Section

Figure 7: The phase shift compared to data.

Having studied the spectral modification of the in the medium we are now in a position to investigate how these changes affect the cross-section. We aim to set up a dynamical framework wherein medium effects can be implemented using thermal field theoretic methods and which at the same time is normalized to the experimental data in vacuum. We consider the interaction (46) and first check with the phase shift data Koch () defining with the partial wave amplitude given by . The decay width which follows from the imaginary part of eq. (42) is


where the c.m. momentum . As seen in fig. 7, a reasonable agreement is obtained using MeV and MeV.

Figure 8: The -N elastic scattering cross section with medium effects.

Next we evaluate the matrix elements for elastic scattering in the isospin basis in which we replace the free vacuum propagator by an effective one containing the vacuum self energy due to the loop diagrams mentioned above. Averaging over isospin, the squared invariant amplitude for the process is given by


where , and are given by


in which


where is defined in (14).

The cross-section given by turns out to be in good agreement with the isospin averaged total elastic cross-section given in Prakash () (obtained using phase shift and inelasticity data from Bareyre () and Koch ())up to about 1.5 GeV as seen from the solid curve in fig. 8. It is to be noted at this point that we have considered only exchange in the evaluation of elastic scattering amplitude with the aim of fixing the parameters (see Appendix B) and thus obtaining a baseline for estimating the effect of the modified propagator on the cross-section. For a more general treatment it is necessary to consider the exchange of as well as nearby resonances like the Roper(1440), etc. in the evaluation of the scattering amplitudes. In such a case, however, it could be quite challenging to obtain a satisfactory agreement with the experimental data especially in the region beyond the peak.

Having thus normalized the framework with the experimental data we now turn on the medium effects. We replace the vacuum self energy in the above expressions by the in-medium ones evaluated in the real-time formalism described above. A significant suppression of the peak with increasing temperature is obtained owing to the increase in the imaginary part due to reasons explained earlier. The small upward shift at higher baryon densities comes from the small positive contribution of the real part of the self-energy. As seen in fig. 4, there are substantial cancellations between the contributions from various loops depending essentially on the attractive or repulsive nature of the effective interactions considered.

Iv The shear viscosity of a gas

Transport coefficients can be obtained in (a) the kinetic theory approach using the transport equation and (b) the diagrammatic approach using Kubo formulae which relates them to retarded two-point functions. The latter was used in Sabya () to obtain the shear viscosity of a pion gas. However, for our present purpose which is to highlight the effect of the in-medium cross-section on the shear viscosity, the kinetic theory approach is more suited (see e.g Mitra1 (); Mitra2 ()). The transport equation describing the evolution of the phase space density of pions and nucleons in a hadronic gas mixture slightly away from local equilibrium is given by


where is the particle velocity and is the collision integral. The distribution function for such a system assumes a form where the local equilibrium distribution function is given by , the plus and minus signs correspond to nucleons and pions respectively and is the deviation function. , and denote the local temperature, fluid velocity and chemical potentials. Assuming all except the particle to be in equilibrium the collision integral simplifies to Gavin ()


where is the relaxation time which characterizes the rate of change of the distribution function due to interaction with the species in the medium. For binary elastic collisions it is Prakash ()




where , and the function if represents a fermion (boson). The dynamical input which goes into the determination of the distribution function appears in the interaction rate .

To extract the shear viscosity we turn to the energy-momentum tensor. For small gradients of the local fluid velocity the shear dissipative part is well-known to be Gavin ()


Now, in terms of the distribution function the correction to the component of the stress-energy tensor is given by Itakura ()


From (29) and (30) we get to lowest order


where in the last line we retained only the (traceless) part appropriate for shear viscosity. Putting this in (34) and equating with (33) we obtain the shear viscosity of the pion-nucleon mixture,

Figure 9: The mean relaxation time of nucleons (a) and shear viscosity (b) of the gas as a function of for two values of . The legend ’Vacuum’ indicates the use of cross-section evaluated in vacuum in eq. (32).

The shear viscosity is thus made up of contributions from the pion and nucleon components which are coupled through the momentum dependent relaxation times and . The relative importance at a given value of and is a consequence of interplay between the phase space factors as well as scattering cross-sections. In order to focus on the effect of the in-medium cross-section we take the and cross-sections in vacuum.

The mean relaxation time of the specie is defined in terms of the thermal average of the momentum dependent inverse relaxation time . With


the mean relaxation time is given by . We plot in fig. 9(a) the mean relaxation time of nucleons as a function of for two values of with and without medium effects. The features of the numerical results can be understood by realizing that the relaxation time for binary collision approximately goes as so that . Thus, the increase of density of the species, in this case pions, with plays the dominant role and accounts for the decreasing nature of the curves. It also follows that the relative increase of nucleon density for larger results in a relative decrease of the relaxation time. Again, since the pion density is considerably more than the nucleon density and is also much larger than , the nucleon relaxation time is expected to be dominated by the pion component. Since is smaller compared to vacuum, as seen in fig. 8, the relaxation time is larger in the medium. Similar arguments hold also for the pion relaxation time. Its magnitude is decided by the (vacuum) cross-section which being much larger, overshadows the medium dependence of the cross-section. We thus do not show it separately. These features as well as the results with the vacuum cross-section are quite in agreement with Prakash ().

In fig. 9(b) the shear viscosity is plotted as a function of . As discussed above the behavior of the pion and nucleon components and their relative magnitude decides that of the mixture. For lower nucleon densities the pion component dominates the viscosity of the mixture. As the nucleon density increases there is a substantial increase in the nucleon component and a decrease in the pion component, the sum being more for MeV compared to 200 MeV. This feature is irrespective of the cross-section (vacuum or medium) and is due to the interplay between the relative abundance of the species and the magnitude of their interaction cross-section. The fact that the shear viscosity in the medium is more than that in vacuum can be understood in terms of the relaxation times which are larger basically due to the lower cross-section in the medium. The increase in magnitude of the medium effect with temperature and nucleon density can be attributed to the phase space factors.

V Summary and Outlook

In this work we have studied the spectral modification of the baryon in the medium. The self-energy was evaluated from one-loop graphs comprising of , , and using the real-time formulation of thermal field theory. In addition to the contributions from decay processes occurring above thresholds which arise due to the usual thermally weighted unitary cut in the complex plane there are significant contributions coming from the Landau type discontinuities in the low region stemming from scattering processes leading to the absorption of in the medium. The cross-section is then evaluated with the effective propagator of the leading to a suppression at finite temperature and density with no significant shift of the peak position. This is expected to have non-trivial consequences on the mean free path and should consequently affect the thermalization rate of pions Barz () and nucleons produced in heavy ion collisions. We finally make an estimate of the shear viscosity of a gas of pions and nucleons using the kinetic theory approach and observe an enhancement corresponding to the in-medium cross-section which increases with temperature and nucleon density. When used as input in the hydrodynamic equations this is expected to have an observable consequence of the space-time evolution of the latter stages of heavy ion collisions.

Vi Appendix

vi.1 Diagonalising the fermion propagator

The procedure of diagonalization of the thermal matrices appearing in the real time formulation of thermal field theory only concerns the bosonic or fermionic nature of the field Mallik (). The spin sum which appears in the numerator of the propagator can thus be factored out. The free thermal fermion propagator matrix is written as where the matrix has components


The terms and containing thermal factors are given by



The matrix can be diagonalized as




It can be shown Mallik () that the complete thermal propagator matrix is also diagonalized by . From the Dyson equation (6) it then follows that the self energy matrix is diagonalized by resulting in eq. (7). The diagonal element is given by any one of the four components of as given in eq. (8).

vi.2 The self-energy in vacuum : Lagrangian and parameters

The self energy in vacuum for the one loop diagrams shown in fig. 1 are given by


where and are the scalar and vector propagators in vacuum. The ones for the spin 1/2 and 3/2 fermions are given by and respectively. The vertex factors come from the well-known interactions Krehl ()


where Krehl (); Korpa1 (), , , , and . In the above