## Brazilian Journal of Physics

##
*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007

#### http://dx.doi.org/10.1590/S0103-97332007000100012

**Linear sigma model at finite baryonic density and symmetry breakings **

**F. L. Braghin**

Instituto de Física, Universidade de São Paulo; C.P. 66.318; CEP 05315-970; São Paulo, Brazil

**ABSTRACT**

The linear sigma model at finite baryonic density with a massive vector field is investigated considering that all the bosonic fields develop non zero expected classical values, eventually associated with condensates and corresponding to dynamical symmetry breakings which might occur in the QCD phase diagram. A modified equation for the classical vector field is proposed with its respective solution. Some *in medium* properties of the model (mainly masses) are investigated within reasonable prescriptions. In particular the behavior of the *in medium* pion and sigma masses and a particular way of calculating *in medium* coupling to baryons is investigated. A symmetry radius for finite baryonic densities is proposed and calculated in different ways in terms of the other variables of the model and these different ways of calculating it agree quite well. However, assuming that the pion and sigma masses go to zero close to the restoration of chiral symmetry a too high value for the critical density is obtained r_{c} @ 4.3 r_{0}.

** Keywords:** Linear sigma model; Finite baryonic density; QCD phase diagram

**I. INTRODUCTION**

Matter at high energy densities has been continuously investigated to provide a deepeer understanding of strong interactions. Experimental (relativistic and high energy) heavy ion collisions (r.h.i.c. and h.e.h.i.c.) provide very important collection of data to construct this knowledge. Asymptotic freedom is a key property in this program. With different approaches, hadronic models based on quantum chromodynamics (QCD) at finite energy density have been extensively investigated and different condensates are usually expected to appear at different energy densities. In the vacuum, the lightest strong interacting particles are known to respect, approximatedly at least, chiral symmetry *SU _{L}*(2) ×

*SU*(2) which is expected to be spontaneously broken down to

_{R}*SU*(2). This is expected to give rise to a scalar (quark-antiquark) condensate, which might have relevant effects in experimental conditions and in realistic calculations, it rearranges the theory [1-3]. As a consequence, the hadronic properties which depend on this condensate (which can be the order parameter or proportionally to it [4]) are expected to vary with energy density. In particular the behavior of the rho vector meson has been considered as a possible signature of the chiral symmetry restoration via dilepton emission [5] although recent experimental analysis of the rho vector meson spectral function shows to be incompatible with energy shift [7].

In this work the Linear Sigma Model (LSM) [8] at finite baryonic density, r_{B}, is investigated with a massive classical vector field. It is based in [9,10]. All the mesons in the model are considered to develop classical counterparts. There is a renewed interest in the pseudoscalar condensation indeed [9,11,12]. It provides a way of enhancing CP violation at finite density [13]. The exact field equations and the stability equation are truncated such as analytical solutions are obtained by considering particular prescriptions for the stability condition. The numerical solutions have a self consistency although the so-called "full self consistency" is only achieved in a level of approximation and for the interactions considered in a model, not in the complete self consistency of the exact realistic quantum theoretical many body problem which still is too difficult to obtain. The truncations in the effective action, done in the next section, are based in the following considerations: (1) the effective potential of spin zero bosons keeps the same form of that at the tree level calculation (i.e., quantum fluctuations basically rearrange the tree level model), (2) each component of the system, i.e. baryons/ spin zero bosons/ spin one fields, have nearly independent stability conditions. Hopefully this assumption might go along with the observation of different slope parameters and temperature freeze-out [14] for each of the hadrons emerging from relativistic heavy ion collisions - and eventually different contributions for the corresponding Hydrodynamics. The corresponding (dynamical) equation for each of the fields is satisfied. The complete numerical investigation of the results will be presented elsewhere [9,15]. Several properties of in medium hadrons are investigated, namely scalar and pseudoscalar meson masses and couplings, and their relation to the behavior of the (chiral) symmetry is worked out within a particular prescription which provides results in agreement with the usually expected behavior. A symmetry radius is defined for the investigation of the symmetry behavior and its estimation (and dependence with the baryonic density) is done in several ways which yield very close result. In spite of being a quite simple model, without several degrees of freedom which should be relevant at high energy densities, results show to be consistent. Assuming that pion and sigma masses go to zero close to the density in which chiral symmetry would be restored a too high critical density is obtained. This can either signal that the model is too simple for describing Physics at too high densities or that their masses should not be expected to be zero close to the phase transition.

**II. THE LINEAR SIGMA MODEL AT FINITE r _{B}**

The Lagrangian density of linear sigma sodel (LSM) with baryons, *N _{i}*(

**x**), sigma and pions, (s,p), covariantly coupled to a vector field,

*V*

_{µ}, is given by [8]:

where the covariant derivatives are given in [8] and they will not be completely considered in this communication. The other terms and parameters are standard [8,9]. The introduction of a chemical potential, with an extra term d = - *g*_{0} µ_{chem}* N*, is nearly equivalent to a shift of the classical temporal component of the vector field *V*_{0} coupled to the nucleons. This field however is a dynamical degree of freedom (d.o.f.) and will be treated as such. Since the condensates (such as º < s > ) depend on the density, so do most of the hadronic masses. Part of the baryon masses are considered to come from the the coupling to the scalar mesonic field and part from an explicit mass term for the baryons in the Lagrangian: *M*^{*} = *M ± g _{S}* .

The spin zero fields will be treated in the framework of the variational Gaussian approach with a truncation [3,9]. With the truncation of the effective potential of sigma and pion it can be written keeping the same form of the tree level effective potential. The equations for expected values of the sigma and pion are found accordingly, and shifts in these classical parts are considered due to the rearrangement brought by quantum fluctuations, such as ® .

The total energy density is written in terms of the four variational parameters for the fields s,, plus baryonic densities and vector field variables [9]. To investigate the behavior the temporal component of (classical) vector field, the total energy density is varied with respect to *V*_{0}, which is not quantized, instead of using its Euler-Lagrange equation. In this approach, *V*_{0} is found either by writing a reasonable (or exact) expression for r_{B} as a function of *V*_{0}, as it is given below, or *V*_{0} is treated like a variational parameter such that one can determine a parametric function r_{B} = r_{B} [*V*_{0}] which satisfy the equations of mouvement and stability. For this derivation, _{V} was kept constant. This second procedure yields a sort of variational equation for the corresponding parametric dependence. The corresponding (variational) equation for a constant background field component *V*_{0} can be given by:

where the Euler-Lagrange equation for *V*_{0} can be recovered by neglecting the derivative term above.

The stability condition for the ground state, with binding energy *E*_{0}/*A* = /r_{B} < 0, can be written as > 0, where r_{0} is the stability density. The expressions for the energy density and its derivative with respect to r_{B} is separated into three parts such that each component of the hadronic matter satisfies the stability equation above separatedly. With this prescription the solutions for the variational equations of each of the components satisfy the respective stability equation. The reliability of this factorization is not mathematically proven although some arguments for being reasonable are given in the Introduction and in [9]. The resulting equations (prescriptions) are the following:

Where _{V} = *g _{V} V*

_{0}r

_{B}- . The complete set of solutions for the equations will be investigated elsewhere and compared to the exact numerical solutions. The variation of the sigma and pion masses arise from the corresponding classical fields. In the vacuum can go to zero as long as

^{2}® 0 satisfying the Goldstone theorem when

*c*® 0 in the Lagrangian.

**A. Densities, coupling constants and masses**

The baryon fields, which depend on the bosonic fields through the Dirac equation coupled to the mesons, are quantized in terms of creation and annihilation operators. The baryonic degrees of freedom sum up into the densities: baryonic (r_{B}), scalar (r_{S}) and pseudo-scalar (r_{ps}) densities. These quantitites will not be explicitely evaluated here although they are partially used below [9,10]. The energy density due to the fermions (antifermions) () and the density of baryons (antibaryons) () can be written, in the leading order, in terms of (1) their momenta for each kind of baryons (*i*), up to the last occupied level with momentum *k _{F}*, and of (2) the classical vector field as [10,16]:

In these expressions *E*_{±} = *g V*_{0} ± are the eigenvalues of the corresponding Dirac equation [10,16]. These expressions still correspond to an approximation and show deviations from the Fermi liquid picture.

The scalar and pseudoscalar densities, which appear in the equations of s, , can be expanded in terms of the scalar and pseudoscalar condensates, for example, as:

where the coefficients are obtained from the expressions calculated with the solution of the corresponding Dirac equation. They are such that r_{S} = r_{ps} = 0 when = _{vac} and = 0, with higher order terms are indicated by *o*(||) and *o*(). By substituting these expressions into the variational equations for the respective condensates [9] it is found that the terms proportional to and can yield contributions to the *in medium* masses of sigma and pions. The terms proportional to ^{3} and ^{3} can produce (effective) contributions for the coupling constant, i.e., l ® l ^{*} = l ± . These two corrections for the in medium effective coupling constant may also be different from each other, eventually leading to different interactions of the pion and sigma in the baryonic medium.

The second equation of prescriptions (3) can produce the same solution of the equation (2). The third of prescriptions (3) can define a symmetry radius in the medium:

In this expression is a constant to be determined from the parameters of the model. In the vacuum: ^{2} = *v*^{2} = as discussed above. Modifications in the equation (3-(*iii*)) will produce different dependences on the baryonic density. A more general symmetry radius, corresponding to particular modifications of the corresponding differential equation due to diverse couplings for example, might be written as: where , , *c, *g are constants to be related to the parameters of the model.

One way of calculating is found by assuming that this symmetry radius is valid over a range of baryonic densities. This is a crude approximation because, heavier hadrons as well as quark and gluon d.o.f. are expected to be relevant for high energy densities. In the high density when chiral symmetry should be restored: = ® 0. This critical density is written as r_{c} = *u* r_{0},. Thus in this point: = . At the saturation density (r_{B}) the expression for can be written as:

Four values are considered: (i) *u* = 2, (ii) u = 3, (iii) *u* = 3.5 and (iv) *u* = 4. Considering the branch of solutions for which ^{2} + ^{2} < ^{2}, at r_{B} = r_{0}, it follows respectively:

Other solutions are not presented. Since the squared value ^{2} is a scalar which appear very often in the expressions, it may be that , i.e., a pion classical field could be responsible for modifications in the pion decay constant and consequently measurable, even if competing with other effects. The topological Skyrme model can provide some argument in favor of such interpretation for a classical pion field inside hadrons [17]. These expressions may be therefore useful for relating descriptions of different ranges of the matter phase diagram.

The values obtained for from estimates (8) are respectively given by:

Another way of estimating is shown by considering the meson masses in the medium. With the expressions for meson masses in terms of the classical fields and *v* [8,9], the symmetry radius can be written as: , where at a given density r_{B}. In these expressions the coupling l was also kept constant (and positive) and *c* = 0, in the Lagrangian term. Two possible behaviors are obtained in this picture for the restoration of chiral symmetry: the sum of these masses may decrease or increase depending on the sign of . For r_{0} = 0.15fm^{-3} and @ - 0.15fm^{-} the above expression yields approximated values ()^{2} (r_{0}) @ ( 1 ±0.53) (r_{B} = 0). If one considers that the pion and sigma masses disappear close to the chiral symmetry restoration point (i.e., if (µ^{*})^{2} ® 0), with the values above we obtain that r_{c} @ 4.3 r_{0}. Seemingly it is a too high baryonic density and the reasons are quite apparent. Firstly, as emphasized above, the present work only takes into account the light sector of hadrons and it does not consider quark and gluon degrees of freedom. Furthermore, it is a controversial subject whether pion and sigma masses (two point Green's functions) should be expected to be so close to zero (as it was assumed to obtain such high value for the critical density) close to (and at the) deconfinement critical point.

These ways of calculating provide crude (but interesting and curious) estimations. The corresponding *in medium* hadron properties are qualitatively in agreement with other estimations [5]. The inclusion of other relevant d.o.f. will be presented elsewhere as well as a corresponding calculation at finite temperature.

**B. Summary and Conclusions**

In this work some aspects of the Linear Sigma Model at finite baryonic density were investigated with a massive classical vector field, based in [9,10]. All the mesons in the model were considered to develop classical counterparts. In part this is due to independent new investigations on pseudoscalar condensates which have shown a renewed interest in the pseudoscalar condensation indeed [9,10,12]. The exact field equations and the stability equation were truncated for obtaining analytical solutions which capture the expected behavior of the system. These solutions have a self consistency although the so-called "full self consistency" is only achieved in a level of approximation and for the interactions considered in a model, not in the complete self consistency of the exact realistic quantum theoretical many body problem which still is too difficult to obtain. The truncations in the effective action, done in the next section, are based in the following considerations: (1) the effective potential of spin zero bosons keeps the same form of that at the tree level calculation (i.e., quantum fluctuations basically rearrange the tree level model), (2) each component of the system, i.e. baryons/ spin zero bosons/ spin one fields, have nearly independent stability conditions. Hopefully this assumption might go along with the observation of different slope parameters and temperature freeze-out for each of the hadrons emerging from relativistic heavy ion collisions - and eventually different contributions for the corresponding hydrodynamic. The corresponding (dynamical) equation for each of the fields are satisfied. The complete numerical investigation of the results will be presented elsewhere [9,15]. Several properties of in medium hadrons were investigated, namely scalar and pseudoscalar meson masses and couplings, and their relation to the behavior of the (chiral) symmetry is worked out within a particular prescription which provides results in agreement with the expected behavior. A symmetry radius was defined for the investigation of the symmetry properties and its estimation (and dependence with the baryonic density) is done in several ways with fair agreement, in spite of being a quite simple model, without several degrees of freedom which should be relevant at high energy densities. Related aspects to matter-antimatter asymmetry in relativistic heavy ion collisions and in the Early Universe will be discussed and investigated elsewhere [16].

**Acknowledgements**

This work was partially supported by FAPESP, Brazil. F.L.B. thanks brief discussions with L.McLerran, R. Rapp, D. Zschiesche, H. Stoecker, R. Pisarski, M. Munhoz and G. Krein.

[1] S. Weinberg, *The Quantum Theory of Fields -II*, Cambridge University Press, (1996). [ Links ]

[2] F.L. Braghin, Doctoral thesis, (1996) unpublished. [ Links ]

[3] P.M. Stevenson, B. Alles, and R. Tarrach, Phys. Rev. D **35** 2407 (1987). [ Links ]F.L. Braghin, Phys. Rev. D **57**, 3548 (1998); [ Links ]Phys. Rev. D **64**, 125001 (2001). [ Links ]F.L. Braghin and F.S. Navarra, Phys. Lett. B **508**, 243 (2001). [ Links ]

[4] H. Umezawa, Advanced Field Theory, AIP Press, (1995). [ Links ]J. Stern, Nucl. Phys B-Proc. Suppl. **64**, 232 (1998). [ Links ]

[5] T. Hatsuda and T. Kunihiro, nucl-th/0112027. [ Links ]G. Chanfray, M. Ericson, and P.A.M. Guichon, Phys. Rev. C **63**, 055202 (2001). [ Links ]D. Roder, J. Ruppert, and D.H. Rischke, Phys. Rev. D **68** 016003 (2003). [ Links ]J. Wambach, Prog. in Part. and Nuc. Phys. **50** 615 (2003). [ Links ]R. Rapp and J. Wambach, Adv. Nucl. Phys. **25**, 1 (2000). [ Links ]D. Jido, T. Hatsuda, and T. Kunihiro, Phys. Rev. D **63**, 011091 (2001). [ Links ]J. Delorme, G. Chanfray, and M. Ericson, Nuc. Phys. A **603**, 239 (1996). [ Links ]L.S. Celenza, H.S. Wang, and C.M. Shakin, Phys. Rev. C **63**, 025209 (2001). [ Links ]G.E. Brown and M. Rho, Phys. Rep. **269**, 333 (1996). [ Links ]N. Kaiser, S. Fritsch, and W. Weise, Nuc. Phys. A **697**, 255 (2002). [ Links ]P.K. Panda and F.L. Braghin, Phys. Rev. C **66**, 055207 (2002). [ Links ]

[6] K. Suzuki *et al*, Phys. Rev. Lett. **92**, 072302 (2004); [ Links ]K. Yokokawa, T. Hatsuda, A. Hayashigaki, and T. Kunihiro, hep-ph/0204163. [ Links ]

[7] Damjanovic for the NA60 Collaboration, arXiv: nucl-ex/0609026 (2006). [ Links ]

[8] S. Gasiorowicz and D.A. Geffen Rev. of Mod. Phys. **41**, 531 (1969). [ Links ]

[9] F.L. Braghin, pré-print IF-1570/2001; hep-ph/0206107; [ Links ]IF-2006. F.L. Braghin, in *HADRON Physics 2002*, ed. by C. Vasconcellos *et al* World Scientific (2002). [ Links ]

[10] F.L. Braghin, Int. Journ. of Mod. Phys. D **13**, 1267 (2004). F.L. Braghin, accepted to be published in Int. Journ. of Mod. Phys. D. [ Links ]

[11] H. Shen, S. Tamenaga, H. Toki, nucl-th/0409057. [ Links ]

[12] A.B. Migdal *et al*, Phys. Rep. **192**, 179 (1990). [ Links ]

[13] See, for example, Y. Mori, K. Saito Phys. Lett. B **552**, 21 (2003). [ Links ]

[14] For example in: F. Grassi, Braz. J. Phys. **35**, 52 (2005) [ Links ]and references therein. HADRON-RANP Joint Meeting, RJ, Brazil, 2004, ed. by M. Bracco *et al*, AIP Conference Proceedings 739, APS (2004). [ Links ]

[15] F.L. Braghin, pre-print IF-USP (2005). [ Links ]

[16] F.L. Braghin, works presented in: ISMD 2006, Paraty, RJ, Brazil, September 2006; [ Links ]Physics in Collisions, Buzios, RJ, Brazil, July 2006 -e-Proceedings: arXiv:hep-ph/0609242; [ Links ]18 th International Conference on Few Body, Santos, SP Brazil, August 2006. [ Links ]

[17] For example: T.H.R. Skyrme, Proc. Roy. Soc. Ser. A **260**, 127 (1961). [ Links ]G.S. Adkins, in *Chiral Solitons* ed. by G.S. Adkins, Add. Wesley, 1988. [ Links ]F.L. Braghin, I.P. Cavalcante, Phys. Rev. C **67** 065207 (2003) and references therein. [ Links ]

Received on 29 September, 2006