Advances in High Energy Physics

Volume 2017, Article ID 4135329, 7 pages

https://doi.org/10.1155/2017/4135329

## Chiral Phase Transition in Linear Sigma Model with Nonextensive Statistical Mechanics

Key Laboratory of Quark & Lepton Physics (MOE) and Institute of Particle Physics, Central China Normal University, Wuhan 430079, China

Correspondence should be addressed to Ben-Wei Zhang; nc.ude.uncc.liam@gnahzwb

Received 9 September 2017; Accepted 7 November 2017; Published 29 November 2017

Academic Editor: Edward Sarkisyan-Grinbaum

Copyright © 2017 Ke-Ming Shen et al. 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. The publication of this article was funded by SCOAP^{3}.

#### Abstract

From the nonextensive statistical mechanics, we investigate the chiral phase transition at finite temperature and baryon chemical potential in the framework of the linear sigma model. The corresponding nonextensive distribution, based on Tsallis’ statistics, is characterized by a dimensionless nonextensive parameter, , and the results in the usual Boltzmann-Gibbs case are recovered when . The thermodynamics of the linear sigma model and its corresponding phase diagram are analysed. At high temperature region, the critical temperature is shown to decrease with increasing from the phase diagram in the plane. However, larger values of cause the rise of at low temperature but high chemical potential. Moreover, it is found that different from zero corresponds to a first-order phase transition while to a crossover one. The critical endpoint (CEP) carries higher chemical potential but lower temperature with increasing due to the nonextensive effects.

#### 1. Introduction

Quantum Chromodynamics (QCD) is a basic theory of describing the strong interactions among quarks and gluons, the fundamental constituents of matter. More and more attentions have already been attracted to the QCD phase transition both theoretically and experimentally [1–8]. Though experimental studies and lattice Monte-Carlo simulations have made it possible to research on the phase diagram quantitatively, there still remains uncertainty at high baryon density region [9]. Consequently, the phase transition is also a vital topic in high-energy physics, where the thermal vacuum created by heavy-ion collisions differs from the one at zero temperature and chemical potential [7]. In order to study certain essential features of it, the linear sigma model has been proposed to illuminate the restoration of chiral symmetry and its spontaneous breaking [8].

Near the phase transition boundary, one should be cautious when using the Boltzmann-Gibbs (BG) statistics for the appearance of critical fluctuations due to a large correlation length [10]. It is of interest to investigate the phase transition in the formalism beyond conventional BG statistical mechanics. Recently nonextensive statistics firstly proposed in [11] has attracted a lot of attention and discussions [12]. In this formalism, instead of the exponential function, a generalized -exponential function is defined as [11, 12]where the parameter is called the nonextensive parameter, which accounts for all possible factors violating assumptions of the usual BG case. Its inverse function is also listed [11, 12]Both of them return to the usual exponential and logarithm function with .

The purpose of this paper is to clarify the nonextensive effects on physical quantities of the chiral phase transition in the generalized linear sigma model. We focus on the situations where both of temperature and chemical potential are not vanished, which then indicates the influence of the Tsallis distribution on the whole phase diagram in the plane. Whereas the nonextensive parameter is still a phenomenological parameter [12], not only the case of but also , in this work, is computed. For comparisons, we have presented discussions on this issue for finite temperature but vanishing chemical potential [13, 14]. Given the consistency of nonextensive generalizations with the initial BG approaches, we also list the results of which were investigated [15]. We close our researches with the comparisons to the nonextensive Nambu Jona-Lasinio model (-NJL model) [16, 17] of the critical endpoint (CEP), whose location is still the hot topic for experiments as well as its theoretical researches [3, 5].

This paper is organized as follows. In Section 2 we introduce the theoretical framework, where the nonextensive -linear sigma model is stated. Their consequences for various thermodynamic quantities with different nonextensive parameters, , are explored in Section 3; more detailed discussions on the results are also contained. Section 4 is our brief summary and outlook.

#### 2. Theoretical Framework

Within the linear sigma model, the chiral effective Lagrangian with quark degrees of freedom reads [15, 18, 19]where stands for the spin- two flavors’ light quark fields and the scalar field and the pion field together form a chiral field , with its potentialConsidering the obvious symmetry breaking term , is invariant under chiral transformations. The chiral symmetry is spontaneously broken in the vacuum with the expectation values: and , where MeV is the pion decay constant. By the partially conserved axial vector current (PCAC) relation [15], the quantity with the constant , where MeV is the pion mass. The coupling constant is fixed as by , where MeV is the sigma mass. Another coupling constant is usually determined by the requirement of the constituent quark mass in vacuum, , which is about of the nucleon mass, leading to . [15]

In order to investigate the temperature and the chemical potential dependence in this model, let us consider a system of both quarks and antiquarks in the thermodynamical equilibrium. Here quark chemical potential . And the grand partition function goes likewhere with being the volume of the system.

Thus the grand canonical potential can be obtainedwith the (anti)quark contribution beingwhere is the internal degrees of freedom of quarks and is the valence (anti)quark energy, with the mass of constituent (anti)quark defined asHere the first divergent term of is absorbed in the coupling constant in the results which comes from the negative energy states of the Dirac sea.

It is inadequate to apply naively the BG statistics in such a system; critical fluctuations of energy and particle numbers will appear as well as a large correlation. In order to investigate the phase transition of systems departing from the classical thermal equilibrium, the nonextensive statistics [11] are introduced. The so-called Tsallis entropy and density matrix are given, respectively, as and , where is the Boltzmann constant (set to 1 for simplicity next), describes the degree of nonextensivity, and is the corresponding generalized partition function.

Recently these generalized statistics have been of great interest theoretically [20–22] and widely applied in many fields [23–26]. In the following, we investigate the linear sigma model within the nonextensive statistics. Firstly we rewrite the partition function of (5) aswhere the -exponential is seen in (1). Considering the -thermodynamics [27], we haveBefore carrying it out, we give out the generalized identities with respect to -sums and integralsand the generalized sum over , in our assumptions,where is used and the index appears because of the dualityIntegrating over and dropping terms that are independent of and , we finally obtain

In our calculations within the mean-field approximation, we follow [15, 28] where the expectation value of the pion field is set to zero; . By solving the gap equationthe value for constituent (anti)quark mass can be determined, which will be also affected by different nonextensive parameters, . Here we have replaced and in the exponent by their expectation values in the mean-field approximation.

With such a -thermal effective potential, we then explore the nonextensive effects on the physical quantities, as well as the phase transition, in the linear sigma model. The numerical results will be shown in the next section.

#### 3. Results and Discussion

In virtue of the fact that there still exist fierce controversies over the possible interpretations of the nonextensive parameter , we shall discuss the nonextensive effects in the -linear sigma model for both the and case. Meanwhile, we give out the result of as the baseline for better understanding.

It is worthy to note that the value of nonextensive parameter cannot be much smaller than since, in the expression of (14), the corresponding generalized exponentialwhere the part of the base should be larger than : ; namely,Thus some upper limitation of energy of the integral in (7) should be given in case of divergence when . On the other hand, too much smaller values of are not necessary to be computed physically during our investigation on the nonextensive effects on the phase transition. Therefore, here we just list the results of and as baseline.

We start our discussions with presenting in Figure 1 the resulting thermodynamical potential as a function of , the constituent (anti)quark mass. Different evidently results in a large change of the thermodynamical potential which shows that the effects caused by nonextensivity are quite strong whether the quark chemical potential vanishes or not.