Physics Research International

Volume 2014 (2014), Article ID 435023, 8 pages

http://dx.doi.org/10.1155/2014/435023

Research Article

Linear Sigma Model at Finite Temperature and Baryonic Chemical Potential Using the -Midpoint Technique

M. Abu-Shady

Department of Mathematics, Faculty of Science, Menoufia University, Menoufia, Shiben El-Koom, 32511, Egypt

Received 6 October 2013; Accepted 23 January 2014; Published 6 March 2014

Academic Editor: Anand Pathak

Copyright © 2014 M. Abu-Shady. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A baryonic chemical potential () is included in the linear sigma model at finite temperature. The effective mesonic potential is numerically calculated using the -midpoint rule. The meson masses are investigated as functions of the temperature () at fi…xed value of baryonic chemical potential. The pressure and energy density are investigated as functions of temperature at fi…xed value of . The obtained results are in good agreement in comparison with other techniques. We conclude that the calculated effective potential successfully predicts the meson properties and thermodynamic properties at finite baryonic chemical potential.

1. Introduction

Quantum Chromodynamics (QCD) is a fundamental theory of strong interactions since it is renormalizable. However, thermodynamics of QCD is not well understood because of its nonperturbative nature. In particular, QCD phase diagram is an essential for understanding not only natural phenomena such as compact stars and the early universe but also laboratory experiments such as relativistic heavy-ion collisions. The quantitative calculations of phase diagram from first-principle lattice QCD (LQCD) have the well-known sign problem when the baryon chemical potential is included [1].

At nonzero baryonic chemical potential, standard numerical lattice simulations do not work well. Therefore our knowledge of the QCD phase diagram at nonzero and relies exclusively on effective models. One of the effective models for describing baryon properties is the linear sigma model, which was suggested by Gell-Mann and Levy [2] to describe nucleons interacting via sigma and pion exchanges. The model is extended to describe the quark interactions by Birse and Banerjee [3]. The model is modified to provide good description of baryon properties at zero temperature as in [4–9]. At finite temperature and zero baryonic potential, the model gives a good description of meson properties and thermodynamic properties as in [10–17].

The baryonic chemical potential is studied in the linear sigma model at finite temperature as in [18–20]. Also, the isospin chemical potential is studied in the same model by using the Cornwall-Jackiw-Tomboulis formalism [21].

The aim of this work is to calculate the effective mesonic potential in the presence of the baryonic chemical potential at finite temperature using the -midpoint rule. The present paper is extended from the previous work [17] at finite baryon chemical potential. The calculated effective potential is used to obtain meson properties, the pressure, and the energy density.

This paper is organized as follows. In Section 2, the linear sigma model at zero temperature and baryonic chemical potential is briefly explained. In Section 3, the effective mesonic potential is presented at finite chemical potential. Numerical calculations and results are discussed in Section 4. Summary and conclusion are presented in Section 5.

2. The Linear Sigma Model

The interactions of quarks via the exchange of - and -meson fields is given by the Lagrangian density [3] as follows:

with

is the meson-meson interaction potential where and are the quark, sigma, and pion fields, respectively. In the mean-field approximation, the meson fields are treated as time-independent classical fields. This means that we replace the power and products of the meson fields by corresponding powers and the products of their expectation values. The meson-meson interactions in (2) lead to hidden chiral symmetry with taking on a vacuum expectation value

where  MeV is the pion decay constant. The final term in (2) is included to break the chiral symmetry explicitly. It leads to the partial conservation of axial-vector isospin current (PCAC). The parameters can be expressed in terms of , meson masses, as

3. The Linear Sigma Model at Finite Temperature and Baryonic-Chemical Potential

The purpose of this section is to include the finite temperature and chemical potential in (1); the effective potential is extended to include the chiral interacting meson fields with quarks fields at the finite temperature and baryonic chemical potential [18]. As follows:

The first term is the potential in the tree level defined in (2) and the second term is the chiral meson fields interacting with quarks at the finite temperature and baryonic chemical potential. Nonzero values of the chiral fields in the chiral broken phase dynamically generate a quark mass . The integration is taken over momentum space (for details, see [18]).

3.1. Numerical Calculations by Using the N-Midpoint Technique

The purpose of this subsection is to calculate the effective mesonic potential, the and -meson masses, the pressure, and the energy density using the -midpoint technique. We rewrite (5) as follows:

Equation (6) is written in the dimensionless form as follows:

with

where

The is the dimensionless form of . Substituting into (7), we obtain

and we write dimensionless form of as follows:

Using midpoint rule (for details, see [22, 23]), we obtain the approximated integral as follows:

where

In [12, 17], the authors applied the second derivation of the effective potential respect to and to obtain the effective meson masses. The first derivative of the effective potential is given by

where

Then, we obtain the effective sigma mass as follows:

where

We obtain the derivation of , and with respect to as follows:

Similarly, we obtain the effective pion mass as follows:

where

We obtain the derivation of , and with respect to as follows:

In [24], the local pressure is given in the nondimensional form

Therefore, the energy density is defined as in [25] as follows:

4. Discussion of Results

In (6), the integration is solved by the -midpoint algorithm as in [17]. The index is taken to get a good accuracy for the numerical integration. In the present work, the parameters of the model such as  MeV,  MeV,  MeV, and the coupling constant at zero temperature are used as the initial parameters at the finite temperature. The calculated effective potential in (12) is used to study the effect of finite temperature on the effective pion and sigma masses, the pressure, and the energy density in the presence of the baryonic chemical potential. In Figure 1, the sigma and pion masses are plotted as functions of temperature in the absence of baryonic chemical potential () and the presence of the baryonic chemical potential  MeV. We note that the sigma mass decreases with increasing temperature while pion mass increases with increasing temperature. The curves crossed at the critical temperature are equal to 191 MeV which is in agreement with the previous work [17].

Figure 1: Sigma and pion masses are plotted as functions of temperature in the absence and presence of baryonic chemical potential.

At finite baryonic chemical potential, it is more difficult to predict the phase transition using lattice QCD. Phenomenology models such as the quark sigma model and the NJL model are used to describe the phase transition. Chen et al. [26] pointed out that the prediction of the chiral phase transition at finite baryonic chemical potential is still under considerable speculation. In present work, the baryonic chemical potential taken is equal to  MeV. The qualitative features of sigma and pion masses are not changed by involving the baryonic potential while the critical temperature shifted from 191 MeV to 99 MeV by including the baryonic potential in the model. In comparison with Scavenius et al. [18], they found that the phase transition is first-order at coupling constant and the phase transition changed to crossover by decreasing coupling constant. Therefore, the phase transition and critical point depend on the parameters of the models. In comparison with Hansen et al. [27], they found the phase transition is crossover when they took the baryonic chemical potential  MeV in the Nambu-Jona-Lasinio (NJL) model and PNJL model. Hence, our results are in agreement with the results of [27].

In Figure 2, the sigma and pion masses are plotted as functions of temperature at coupling constant , , and  MeV. We note that the behavior is not changed with increasing coupling constant . The critical temperature decreases with increasing coupling constant .

Figure 2: Sigma and pion masses are plotted as functions of temperature for two values of coupling constant .

In Figure 3, we examine the effect of sigma mass on the behavior of the effective sigma and pion masses and the critical temperature in the presence of baryonic chemical potential. We note that the behavior is unchanged with strong increase in the sigma mass while the critical temperature shifted to higher value from  MeV to  MeV corresponding to  MeV and  MeV, respectively. The similar effect is obtained in [17] in the absence of baryonic chemical potential. Also, we obtained the similar effect using the logarithmic potential where the same technique is used [28].

Figure 3: Sigma and pion masses are plotted as functions temperature for two values of sigma mass.

Next, we need to examine the behavior of the pressure at finite temperature. In Figure 4, the pressure is plotted as a function of temperature in the absence and presence of the chemical potential. We note that the curve increases continuously with increasing temperature. The values of the pressure at lower-values of temperature are not sensitive in comparison with the values of the pressure at the higher values of temperature. Hence the effect of largest values of temperature has more effect on the pressure. By increasing baryon chemical potential up to  MeV, we note the pressure shifted to higher values. This means the effect of the baryonic chemical potential has more effect on the pressure at higher values of temperature. Berger and Christov [29] found that the pressure increases with increasing temperature in hot medium using the NJL model in the mean-field approximation. Bowman and Kapusta [25] studied the pressure as a function of baryon chemical potential for different values of temperatures. They found that the pressure increases with increasing chemical potential. Hence, the present results of the pressure are in agreement with [25, 29].

Figure 4: The pressure is plotted as a function of temperature for two values of and  MeV.

Next, we need to examine the energy density in the presence of the baryon chemical potential. In Figure 5, the energy density is plotted as a function of temperature at and  MeV; we note that the energy density increases smoothly with increasing temperature at and  MeV. The energy density is not sensitive to the baryonic chemical potential up to  MeV. Therefore, the effect of the baryon chemical potential appears at higher values of temperature. An increase in the baryonic chemical potential leads to increasing energy density as in Figure 6. Berger and Christov [29] studied the energy density in the NJL model and found that the energy density increases with increasing temperature up to  MeV at finite density. Bowman and Kapusta [25] studied the energy density in the linear sigma model. They found that energy density increases as a continuous curve with increasing temperature for a vacuum pion mass of 200 MeV. This indicated that the phase transition is crossover. Therefore, the present result refers to the same conclusion that the phase transition is crossover at pion mass equals 140 MeV. In the present quark model, we included the fermion sector term which is ignored in other works such as in [10, 12, 15]. So we need to examine the effect of this factor on the energy density. In Figure 6, the energy density is plotted as a function of temperature at the coupling constant and , and  MeV. We note that the energy density of quarks decreases with increasing the coupling content due to the bound state is more tight.

Figure 5: The energy density is plotted as a function of temperature for two values of and  MeV.

Figure 6: The energy density is plotted as function of temperature for two values of coupling constant at  MeV.

5. Summary and Conclusion

In this work, we investigated the effect of the baryon chemical potential on the meson properties, the pressure, and the energy density at finite temperature. The effective mesonic potential is calculated by using the N-midpoint method at finite of chemical potential. We found that the behavior of the meson properties is in good agreement in comparison with other models. The pressure and energy density are examined as functions of temperature at finite temperature. A comparison is presented with other works. In addition, we included the fermion sector in the linear sigma model. This sector is ignored in many previous studies such as in [10, 12, 15].

We conclude that the present technique successfully predicts the behavior of the meson properties, the critical point temperature, and the energy density in comparison with other works at finite chemical potential and avoid the difficulty that previously mentioned in [17].

Appendix

In this appendix, we write the basic steps of the -midpoint rule to calculate the integral . We divide the interval into subintervals of length

and the midpoint of the th interval is

The the midpoint rectangle is the rectangle of height over the subinterval . This rectangle signed area

Thus, is equal to the sum of the signed areas of these rectangles

by increasing the parameter , we get a good accuracy for the numerical integration.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

M. Fukugita, “Lattice quantum chromodynamics with dynamical quarks,” Nuclear Physics B, vol. 4, pp. 105–129, 1988.

M. Gell-Mann and M. Levy, “The axial vector current in beta decay,” Il Nuovo Cimento, vol. 16, no. 4, pp. 705–726, 1960.

M. C. Birse and M. K. Banerjee, “Chiral model for nucleon and delta,” Physical Review D, vol. 31, no. 1, pp. 118–127, 1985.

W. Broniowski and B. Golli, “Approximating chiral quark models with linear sigma-models,” Nuclear Physics A, vol. 714, no. 3-4, pp. 575–588, 2003.

M. Abu-Shady, “Effect of the A-term on the nucleon properties in the extended linear sigma model,” International Journal of Theoretical Physics, vol. 48, no. 4, pp. 1110–1121, 2009.

M. Abu-Shady, “Effect of coherent-pair approximation on nucleon properties in the extended linear sigma model,” Acta Physica Polonica B, vol. 40, no. 8, pp. 2225–2238, 2009.

M. Abu-Shady, “Effect of logarithmic mesonic potential on nucleon properties,” Modern Physics Letters A, vol. 24, no. 20, pp. 1617–1629, 2009.

M. Abu-Shady and M. Rashdan, “Effect of a logarithmic mesonic potential on nucleon properties in the coherent-pair approximation,” Physical Review C, vol. 81, no. 1, Article ID 015203, 8 pages, 2010.

M. Abu-Shady, “Chiral logarithmic quark model of N and ∆ with an A-term in the meanfield approximation,” International Journal of Modern Physics A, vol. 26, no. 2, p. 235, 2011.

N. Petropoulos, “Linear sigma model and chiral symmetry at finite temperature,” Journal of Physics G, vol. 25, no. 11, pp. 2225–2241, 1999.

G. Amelino-Camelia and S.-Y. Pi, “Self-consistent improvement of the finite-temperature effective potential,” Physical Review D, vol. 47, no. 6, pp. 2356–2362, 1993.

Y. Nemoto, K. Naito, and M. Oka, “Effective potential of the O(N) linear sigma-model at finite temperature,” European Physical Journal A, vol. 9, no. 2, pp. 245–259, 2000.

M. Abu-Shady, “Meson properties at finite temperature in the linear sigma model,” International Journal of Theoretical Physics, vol. 49, no. 10, pp. 2425–2436, 2010.

M. Abu-Shady and H. Mansour, “Quantized linear sigma model at finite temperature and nucleon properties,” Physical Review C, vol. 85, Article ID 055204, 2012.

J. T. Lenaghan and D. H. Rischke, “The O(N) model at finite temperature : renormalization of the gap equations in Hartree and large-N approximations,” Journal of Physics G, vol. 26, no. 4, p. 431, 2000.

H. C. G. Caldas, A. L. Mota, and M. C. Nemes, “The chiral fermion meson model at finite temperature,” Physical Review D, vol. 63, no. 5, Article ID 056011, 2001.

M. Abu-Shady, “N-midpoint rule for calculating the effective mesonic potential at finite temperature,” International Journal of Theoretical Physics, vol. 52, pp. 1165–1174, 2013.

O. Scavenius, A. Mocsy, I. N. Mishustin, and D. H. Rischke, “Chiral phase transition within effective models with constituent quarks,” Physical Review C, vol. 64, no. 4, Article ID 045202, 2001.

T. Brauner, “Spontaneous symmetry breaking in the linear sigma model at finite chemical potential,” Physical Review D, vol. 74, no. 8, Article ID 085010, 13 pages, 2006.

N. Bilic and H. Nikolic, “Chiral-symmetry restoration in the linear sigma model at nonzero temperature and baryon density,” European Physical Journal C, vol. 6, no. 3, pp. 513–523, 1999.

H. Mao, N. Petropoulos, S. Shu, and W.-Q. Zhao, “The linear sigma model at a finite isospin chemical potential,” Journal of Physics G, vol. 32, no. 11, pp. 2187–2198, 2006.

D. Jogerman, “Investigation of a modified mid-point quadrature formula,” Mathematics of Computation, vol. 20, pp. 79–89, 1966.

F. Stetter, “On a generalization of the midpoint rule,” Mathematics of Computation, vol. 22, pp. 661–663, 1968.

M. Nahrgang and M. Bleiher, “Fluid dynamics with a critical point,” Acta Physica Polonica B, vol. 2, p. 405, 2009.

E. S. Bowman and J. I. Kapusta, “Critical points in the linear sigma model with quarks,” Physical Review C, vol. 79, no. 1, Article ID 015202, 7 pages, 2009.

J.-W. Chen, H. Kohyama, and U. Raha, “Model analysis of thermal UV-cutoff effects on the chiral critical surface at finite temperature and chemical potential,” Physical Review D, vol. 83, no. 9, Article ID 094014, 2011.

H. Hansen, W. M. Alberico, A. Beraudo, A. Molinari, M. Nardi, and C. Ratti, “Mesonic correlation functions at finite temperature and density in the Nambu-Jona-Lasinio model with a Polyakov loop,” Physical Review D, vol. 75, no. 6, Article ID 065004, 2007.

M. Abu-shady, “A new technique for the calculation of effective mesonic potential at finite temperature in the logarithmic quark-sigma model,” Journal of Fractional Calculus and Applications, vol. 3(s), p. 6, 2012.

J. Berger and C. V. Christov, “Phase transition and the nucleon as a soliton in the Nambu-Jona-Lasinio model at finite temperature and density,” Nuclear Physics A, vol. 609, no. 4, pp. 537–561, 1996.