Critical Line Back-Bending Induced either by Finite Nc Corrections or by a Repulsive Vector Channel

We analyze the two flavor version of the Nambu-Jona-Lasinio model with a repulsive vector coupling , at finite temperature and quark chemical potential, in the strong scalar coupling   V G   S G regime. Considering GV = 0, we review how finite Nc effects are introduced by means of the Optimized Perturbation Theory (OPT) which adds a


Introduction
Although most of the results obtained up to now seem to support the quantum chromodynamics (QCD) critical point (CP), an interesting observation against its existence comes from the numerical simulations of QCD at imaginary chemical potential by de Forcrand and Philipsen [1][2][3] which shows that the region of quark masses (m c ) where the transition is presumably of the first order (for quark masses smaller than the physical ones), tends to shrink for small positive values of the chemical potential as shown in the upper panel of Figure 1.Conversely, according to models supporting the critical point, the first order region should expand when the chemical potential increases, so that the physical quark mass point hits the critical line at some finite value of the temperature and chemical potential as shown in the bottom panel of Figure 1.A possible explanation for the disagreement between the "exotic" scenario (Figure 1, upper panel) and the "standard" scenario (Figure 1, bottom panel) has been given in [4,5] where it was suggested that a strong (repulsive) vector coupling may account for the initial shrinkage of the first order region, that would then start expanding again at larger values of the chemical potential leading to the back-bending of the critical surface and the recovery to the CP at the physical quark mass values.As a result, two critical points should appear for a given range of (small) quark masses, as argued by Bowman and Kapusta [6] who investigated the Linear Sigma Model (LSM) including thermal fluctuations and considered small values for the pion mass.A pictorial view of this peculiar situation is given by Figure 2 which illustrates a possible back-bending scenario in the two-flavor case.In the more traditional T − μ plane, these two critical points are located at the end of two first order transi-  tion lines where one of them represents the usual line which starts at zero temperature and chemical potential of the order of the constituent quark mass while the other is an unusual line which starts at zero chemical potential and high temperature [6,7].
In [4] the author has considered the three flavor Nambu-Jona-Lasinio model (NJL) at large-N c with an explicit repulsive vector interaction, with coupling G V , in order to produce a back-bending that would conciliate the lattice results obtained by de Forcrand and Philipsen with most model predictions.It is well known that within the NJL this type of interaction weakens the first order transition line [8] in opposition to the scalar coupling which tends to favor the appearance of first order phase transitions [9].The explicit presence of a vector term was decisive in order to produce the back-bending scenario within the large-N c application of [4].It was explained that the net effect produced by a repulsive vector channel is to add a term like −G V  q 2 (here ρ q represents the quark number density) to the pressure and as a result the size of the first order covers a smaller range of temperatures as compared to the G V = 0 case.
At the same time, the value of the coexistence chemical potential for a given temperature occurs at a higher value when G V ≠ 0 and, as a consequence, the critical end point happens at smaller temperatures to be higher chemical potentials than in the case of vanishing G V .Although such a vector term is known to be important at high densities in theories such as the Walecka model for nuclear matter, its consideration is more delicate within a non renormalizable model such as the NJL where usually
However, fixing G V poses and additional problem since this quantity should be fixed using the ρ meson mass which, in general, happens to be higher than the maximum energy scale set by Λ.Then, G V is usually considered to be a free parameter whose estimated value ranges between 0.25 G S and 0.5G S [11,12].
Alternatively, when going beyond the large-N c (or mean field) level one may induce quantum (loop) corrections which mimic the physical effects caused by a classical (tree) term such as G V .This is precisely what has been observed in an application of the nonperturbative Optimized Perturbation Theory (OPT) method to the two flavor NJL model with vanishing G V [13].The OPT results for phase diagram for this model show that 1 c N corrections induced by this approximation reproduce the same qualitative features obtained by considering the model at large-N c with an explicit repulsive vector channel.The reason is that the OPT two loop contributions add a term like to the pressure.The relationship between the OPT, at G V = 0, and the large-N c approximation, at G V ≠ 0, has been recently investigated in great detail in the framework of the abelian NJL at finite densities and zero temperature in [14].In the context of the eventual back-bending behavior of the critical line in the μ − m c plane, the OPT has also been previously employed with success in [7].There, the strategy was to use very high values for G S in order to obtain a T − μ phase diagram dominated by first order chiral transitions only.
Then, the OPT with its term was used at different quark mass values showing that, in this case, two critical points emerge at low m c due to the weakening of the first order line at intermediate μ values leading to the back-bending behavior observed in without the need to explicitly include a vector channel in the lagrangian density.
Note that a repulsive vector type of coupling was not explicitly considered in the two flavor LSM application performed by Bowman and Kapusta which, on the other hand, was carried out beyond the mean field level through the consideration of thermal fluctuations.
In the present work, we extend the comparison between the OPT (at G V = 0) and the large-N c approximation (at G V ≠ 0) to the non abelian NJL model at finite temperature and density in the strong coupling and small quark mass regime showing that, as expected, both methods agree from the qualitative point of view leading to a back-bending which would be completely missed by a standard large-N c evaluation.
Our results also emphasize the importance played by 1 c N terms which are easily taken into account by the OPT so that this method may be viewed as a robust alternative to investigate nonperturbative effects related to the chiral transition of strongly interacting matter.The work is organized as follows.In the next section, we perform a large-N c application to the two flavor NJL version in the strong coupling regime for the G V ≠ 0 case.In Section 3 we review the OPT results, at G V = 0, which were originally obtained in [7].We then compare, in Section 4, the analytical and numerical results obtained with the two different model approximations.Our conclusions are presented in Section 5.

The NJL in the Strong Coupling Regime
The standard version of the two flavor Nambu-Jona-Lasinio model lagrangian density L with a repulsive vector channel reads [10,15] where  (a sum over flavors and color degrees of freedom is implicit) represents a flavor isodoublet (u and d type of quarks) N c -plet quark fields while  are isospin Pauli matrices.As emphasized in [16], the introduction of a repulsive vector interaction term of the form 1) is also allowed by the chiral symmetry.Such a term can become important at finite densities, generating a saturation mechanism depending on the vector coupling strength that provides better matter stability [10,16].Here, we will show that this term also influences the phase diagram, especially at the low temperature and high density region.A standard parametrization for this model is Λ = 587.9MeV, G S Λ 2 = 2.44, and m c = 5.6 MeV so that, with these inputs, one obtains f π = 93 MeV, m π = 135 MeV, and M = 400 MeV at T = 0, and μ = 0 [10].
However, as will be shown in the next subsection, in order to simulate the back-bending behavior in the present model we will keep Λ = 587.9MeV considering 2 4 S G   while varying the current quark mass from m c = 0 to m c = m phys = 5.6 MeV.As discussed in [7] this nonstandard choice does not affect very much the predicted values for observables such as f π , m π and  .
On the other hand the quark effective mass, which is directly proportional to G S , assumes very high values (around 800 MeV > Λ).However, this is not a problem for our present purpose of simulating the back-bending in a qualitative way (see [7] for a more complete discussion).

Thermodynamical Potential at Large-N c with Finite G V Contributions
The large-N c (or MFA) evaluation of the thermodynamical potential within this model is standard and yields [10]   , , , where the dressed "free gas" term is given by the scalar density is given by and the quark number density is given by The integrals appearing in the above equations are defined by and , where Here, the divergent contributions corresponding to the first term on the right hand side of Equations ( 6) and ( 7) are regulated by Λ.The effective quark mass, M, and the effective chemical potential, µ, are obtained from solving the following coupled self consistent equations and Let us now review the main effects of the repulsive vector interaction in the phase diagram by considering the standard parametrization.Figure 3 shows the situation for G S Λ 2 = 2.44, Λ = 587.9MeV and m c = 0 (chiral limit) at G V /G S = 0, 0.1, and 0.4.As expected the G V term has little effect at low chemical potential values where the second order chiral transition dominates since m c = 0.
One also notices that, as G V increases, the first order chiral transition line weakens so that the tricritical point occurs at smaller temperatures.It is also clear that for a given the coexistence chemical potential takes place at higher values with increasing G V .This scenario is also observed away from the chiral limit except that now the second order transition is replaced by a cross over region and the tricritical point turns into a critical end point.This type of phase diagram where a CP naturally appears at the physical quark mass point is predicted by most model approximations.Now, in order to simulate the lattice results by de Forcrand and Philipsen we first need to obtain a first order phase transition at μ = 0 for low m c values so that at vanishing densities our results would be consistent with the lattice results furnished by the Columbia plot  0 c m   [17].Within the NJL this is easily achieved by increasing G S [9] so that, at G V = 0 and m c = 0, the whole phase diagram is dominated by a first order phase transition.plane showing that, for G V = 0, the first order line recedes from the temperature axis as m c increases without reproducing the back-bending scenario which would conciliate the de Forcrand and Philipsen results with model predictions.Now, if we turn on the vector interaction, still at m c = 0, the weakening of the first order transition happens in a different way so that two segments of first order chiral transitions appear.One of them is the usual one which starts at T = 0 while the other is an unusual first order line which starts at μ = 0 as Figure 4 shows.So, at the expense of considering a strong G S and a finite G V we have managed to induce the appearance of two tricritical points at vanishing m c which, as will be shown in Section 3, leads to the back-bending scenario.

Thermodynamical Potential at Vanishing G V with Finite N c Contributions
The basic idea of the OPT method is to deform the original lagrangian density by adding a quadratic term like (1 )     to the original lagrangian density as well as by multiplying all coupling constants by δ [7].The new parameter δ is just a bookkeeping label and η represents an arbitrary mass parameter.Perturbative calculations are then performed in powers of the dummy parameter δ which is formally treated as small and set to the original value, δ = 1, at the end 1 .
Therefore, the fermionic propagator is dressed by η which may also be viewed as an infrared regulator in the case of massless theories.After a physical quantity, such as the thermodynamical potential (Ω), is evaluated to the k-order and δ set to the unity only a residual η depend-ence remains.Then, optimal nonperturbative results can be obtained by requiring that Ω(k) (η) be evaluated where it is less sensitive to variations of the arbitrary mass parameter.This requirement translates into the criterion known as the Principle of Minimal Sensitivity (PMS) [18] In general, the solution to this equation implies in self consistent relations generating a nonperturbative coupling dependence.In most cases nonperturbative 1 c N corrections appear already at the first nontrivial order while the large-N c (or MFA) results can be recovered at any time simply by considering N c → ∞.Finally, note that the OPT has the same spirit as the Hartree and the Hartree-Fock approximation in which one also adds and subtracts a mass term.However, within these two traditional approximations the topology of the dressing is fixed from the start: direct (tadpole) terms for Hartree and direct plus exchange terms for Hartree-Fock.On the other hand, within the OPT, the dressed mass term    acquires characteristics which change order by order progressively incorporating direct, exchange, vertex corrections, etc, effects.The differences between these three different methods have been recently discussed in [14].To implement the OPT within the NJL model at G V = 0 one follows the prescription used in [13] to write Then, the order-δ thermodynamical potential can be written as (see [13] for technical details) where now all the integrals defining the quantities ω FG , ρ s , and ρ q are redefined as . Then, for each pair of   , T  values the optimum mass parameter,  , can be obtained by solving the PMS equation given by [13] Note that when N c → ∞ the PMS optimization procedure sets  = −2G S ρ s exactly reproducing the large-N c result (for the standard NJL model) with no ρ q depend-ence.By comparing the OPT thermodynamical potential given by Equation ( 13) with its large-N c counterpart given by Equation ( 2) one notices that the OPT induces a finite N c correction of the form as discussed in [4], the G V term gives a net contribution of the form G V ρ q 2 .Then, one could expect that the two different model approximations given by the OPT (at G V = 0) and the large-N c approximation (at G V ≠ 0) lead to the same qualitative picture of the phase diagram.A recent comparison performed with one flavor at vanishing temperature in [14] showed that this is indeed the case.
Here, we are now in position to extend that comparison to the more realistic two flavor case at finite temperatures.ing of the critical line so that the CP will be recovered at m c = m phys even if initially (at low values of μ) the line bends in such a way which is reminiscent of the "exotic" scenario displayed by the right panel of Figure 1.The physical nature and even the critical exponents of the two different critical points which occur at small m c have been discussed in great detail in [7].From our results it is clear that as m c → m phys the unusual first order line disappears and only the usual "liquid-gas" type of first order line survives in accordance with most model predictions.

Conclusions
We have considered the two flavor NJL model in the strong scalar coupling regime (G S Λ 2  4) in order to compare two distinct model approximations.The first is the traditional large-N c approximation which was applied by explicitly considering a finite repulsive vector interaction (proportional to G V ) which was introduced at the classical (tree) level.The second is the alternative OPT method which was applied to the standard version of the NJL model .

 
0 V G  Our first step towards the simulation of the backbending behavior was to tune G S at G V = 0 so that the large-N c approximation predicts that the first order transition line, which usually starts at T = 0, will touch the T axis at μ = 0 for very small m c values.
For mass values closer to the physical ones, this approximation recovers the expected cross over behavior at small μ with the appearance of a single critical point at intermediate chemical potentials.Next, we have shown that by considering finite values for G V at c phys this single first order transition line splits into two lines in the T − μ plane.One of them is similar to the usual "liquid-gas" line which starts at T = 0 and ends at intermediate temperature values.The other one, which has a more "chiral" behavior according to the analysis of m m  [7], is located at the high temperature region and disappears as m c approaches m phys from below.In this way, we were able to induce the back-bending behavior for two flavors in a manner analogous to the one adopted in [4] for the three flavor case.The OPT results for this strong coupling and small regime obtained in [7] were then reviewed so that a numerical comparison could be performed.At the first non trivial order, this approximation includes one and two loop terms which would belong to the and N type of expansion.In particular, the two loop terms generate a negative contribution to the pressure given by   2 S q f c , where ρ q represents the fermionic density.This term is similar to the net −G V ρ q 2 contribution considered at large-N c G N N   [4].Within the OPT the 1 c N suppressed vector term competes with its scalar counterpart, ρ s , weakening the first order line at intermediate values of μ and enhancing the appearance of two critical points in the T − μ plane for m c values which are smaller than the physical ones.Finally by scanning the values of m c , we have mapped the T − μ phase diagram into the μ − m c plane observing that, for strong couplings, the large-N c approximation, at finite G V , and the OPT, at vanishing G V , predict that the first order transition region shrinks for low values of μ as observed in the lattice simulations of [1].But then, at intermediate chemical potentials, the vector terms −G V ρ q 2 (large-N c ) and (OPT) change the first order phase transition region into a cross over region.Finally, at higher chemical potentials, the first order transition region reappears and then expands as μ is increased.So, our results suggest that even if an initial shrinkage of the first order region is confirmed by lattice simulations, it does not necessary rule out the existence of the CP which is expected to occur at intermediate chemical potentials for physical quark masses.In this case, a back-bending will be observed on the μ − m c plane outlining the importance of a repulsive vector contribution in agreement with N N [4].Our comparison allows us to conclude that similar qualitative results will be obtained either by explicitly considering such a contribution at the classical level as in the large-N c case or by radiatively generating it by going beyond the mean field level.Within the NJL model, the advantage of the second procedure, which can easily be implemented within the OPT, is that it does not require the fixing of G V which is a drawback of the first procedure.The results obtained in the present application support those obtained in [14], for the simpler abelian NJL model at vanishing temperature, showing the robustness of the OPT method.Finally, note that at these non standard high coupling values our results are to be taken only as qualitative predictions.

Figure 1 .
Figure 1.Upper panel: The "exotic" scenario with the shrinkage of the first order transition region which excludes the presence of a CP at the physical point.Bottom panel: The "standard" scenario with the expansion of the first order transition region which provides the appearance of a CP at the physical point.

Figure 2 .
Figure 2. Schematic view of the back-bending scenario in the µ  m c plane.The shaded area represents the region of first order phase transitions; the dashed line represents second order phase transitions while the light part of the figure represents the cross over region.Two critical points (c 1 and c 2 ) appear at m c << m phys and the expected CP is present at m c = m phys .
the integrals are regulated by a momentum cut-off, Λ.Within this model, G S and Λ are usually fixed to reproduce the pion mass  

Figure 4 Figure 3 .
Figure 3.The large-N c phase diagram for a standard parametrization ( = 587.9MeV and G S  2 = 2.44) in the chiral limit (m c = 0).The continuous lines represent first order phase transition and the dashed lines represent second order phase transitions.The solid symbols represent the tricritical points for the ratios G V /G S = 0, 0.1 and 0.4.Note that G V weakens the first order phase transition and, for a given T, shifts the coexistence chemical potential to higher values.

Figure 4 .Figure 5 .
Figure 4.The large-N c phase diagram in the strong coupling regime ( = 587.9MeV and G S  2  4).Upper panel: the chiral limit (m c = 0) where the continuous lines represent first order phase transitions and the dashed lines represent second order phase transitions.The solid symbols represent the tricritical points for the ratios G V /G S = 0, 0.1 and 0.3.Bottom panel: different current mass values (m c /m phys = 0, 0.5, 1.0, 1.5, 2.0 and 2.5) have been used with fixed G V /G S = 0.5.For m c ≠ 0 the dashed lines represent the cross over region.Both panels show that G V weakens the first order phase transition but for small m c and high G S its presence leads to the emergence of two (tri) critical points.Note how m c has a greater influence over the high-T region while G V affects the intermediate to low-T region.

7 .Figure 6 .
Figure 6.The large-N c (at G V /G S = 0.5) and the OPT (at G V = 0) phase diagrams in the strong coupling regime ( = 587.9MeV and G S  2 = 4) for weakly broken chiral symmetry (m c = 0.1 MeV).The figure shows the usual CP (c 2 ) and the exotic high-T CP (c 1 ) at the end of the first order transitions for the OPT (continuous line) and large-N c (dotdashed).For comparison the figure also shows the large-N c result (at G V = 0 and m c = 0) with the dashed line representing the first order transition.

Figure 7 .
Figure 7.The µ − m c plane at strong (scalar) coupling ( = 587.9MeV, G S  2 = 4) showing the back-bending.The shaded area represents the region of first order phase transitions; the dashed dark line represents second order phase transitions while the light part of the figure represents the cross over region.Left panel: The large-N c result at G V /G S = 0.5.Right panel: The OPT result at G V = 0.In both cases one observes two CP at m c << m phys but only the one at high µ survives at the physical point (m c = m phys ).