# Finite-volume and magnetic effects on the phase structure of the three-flavor Nambu–Jona-Lasinio model

###### Abstract

In this work we analyze the finite-volume and magnetic effects on the phase structure of a generalized version of Nambu–Jona-Lasinio model with three quark flavors. By making use of mean-field approximation and Schwinger’s proper-time method in a toroidal topology with antiperiodic conditions, we investigate the gap equation solutions under the change of the size of compactified coordinates, strength of magnetic field, temperature and chemical potential. The ’t Hooft interaction contributions are also evaluated. The thermodynamic behavior is strongly affected by the combined effects of relevant variables. The findings suggest that the broken phase is disfavored due to both increasing of temperature and chemical potential, and the drop of the cubic volume of size , whereas it is stimulated with the augmentation of magnetic field. In particular, the reduction of (remarkably at fm) engenders a reduction of the constituent masses for -quarks through a crossover phase transition to the their corresponding current quark masses. On the other hand, the presence of a magnetic background generates greater values constituent quark masses, inducing smaller sizes and greater temperatures at which the constituent quark masses drop to the respective current ones.

###### pacs:

11.10.Wx, 11.30.Qc, 11.10.Kk## I Introduction

In the recent decades we have witnessed experimental and theoretical progresses in the assessment of strongly interacting matter under extreme conditions. One of the most interesting predictions of its underlying theory, Quantum Chromodynamics (QCD), is that it experiences a phase transition to a deconfined state at sufficiently high temperatures. This mentioned state composed of quarks and gluons is commonly referred to as the quark-gluon plasma (QGP) and it can now be observed in heavy-ion collisions qgpdisc ; rev-qgp . Its thermodynamic and transport properties continues to the forefront of endeavour of community to outline its rich phase structure rev-qgp ; Prino:2016cni ; Pasechnik:2016wkt .

In particular, one interesting feature on the phase structure of strongly interacting matter is the finite-volume effects. There are estimations indicating that QGP-like systems yielded in heavy-ion collisions are of the order of units or dozens of fermi Bass:1998qm ; Palhares:2009tf ; Graef:2012sh . In this sense, the influence of the size of the system on its thermodynamic behaviour and phase diagram has been widely analyzed in the literature through distinct and effective approaches, as Dyson–Schwinger equations of QCD Luecker:2009bs ; Li:2017zny ; Shi:2018swj , quark-meson model Braun:2004yk ; Braun:2005fj , extensions and generalizations of the Nambu–Jona-Lasinio (NJL) model Ferrer:1999gs ; Abreu:2006 ; Ebert0 ; Abreu:2009zz ; Abreu:2011rj ; Bhattacharyya:2012rp ; Bhattacharyya:2014uxa ; Bhattacharyya2 ; Pan:2016ecs ; Kohyama:2016fif ; Wang:2018qyq , and others Gasser:1986vb ; Damgaard:2008zs ; Fraga ; Abreu3 ; Abreu6 ; Ebert3 ; Abreu4 ; Abreu5 ; Abreu7 ; Bao1 ; PhysRevC.96.055204 ; Samanta ; Wu ; Klein:2017shl ; Shi . The main findings of these works point out that the finite volume, combined with other relevant variables, act on the thermodynamic behavior of strongly interacting matter.

In the papers cited above, important questions arise in order to have a detailed description of limited-size physical systems, like the mentioned QGP. One example is the spontaneous symmetry breaking of chiral symmetry, which occurs for systems in the bulk approximation, but is disfavored for finite volume. Thus, this context gives rise to the natural debate about what are the conditions in which an ideal bulk system seems a good approximation for confined systems. In other words, one may wonder about the range of the size of the system in which the finite-volume effects affect the phase structure and the spontaneous symmetry breaking, making modifications in the thermodynamical properties of quarkionic matter.

Specially, in a very recent paper by Wang et al. Wang:2018qyq is discussed the influence of the finite-volume effects on the chiral phase transition of quark matter at finite temperature, based on a two-flavour NJL model approach with a proper-time regularization and use of a stationary-wave condition. It is found that when the cubic volume size is larger than 500 fm, the chiral quark condensate is indistinguishable from that at , which is far greater than the estimates of the size of QGP produced at laboratory and the lattice QCD simulations. Also, they found that when the space size is less than 0.25 fm, the spontaneous symmetry breaking concept is no longer valid.

On the other hand, the importance of a magnetic background on the phase diagram of strongly interacting matter has also been a subject of great interest. The motivation comes from both phenomenology of heavy ion collisions and compact stars, in which a strong magnetic background is produced Kharzeev ; Skokov:2009qp ; Chernodub:2010qx ; Ayala1 ; Tobias ; Heber ; MAO ; Ayala2 ; Mamo:2015dea ; Pagura ; Magdy ; Ayala0 . Taking into account the case of heavy-ion collisions at RHIC and LHC, a notable dependence of phase diagram with the strength of magnetic field is estimated in the hadronic scale, i.e. ( MeV is the pion mass). In this scenario, intriguing distinct phenomena appearing for different ranges of temperatures have been proposed recently by use of effective approaches: for sufficiently small magnetic field, magnetic catalysis (enhancement of broken phase) occurs at smaller temperatures, while the inverse magnetic catalysis effect (stimulation of the restoration of chiral symmetry) takes place at higher temperatures Tobias ; MAO ; Mamo:2015dea ; Pagura ; Magdy ; Ayala0 . It should be also noted that inverse magnetic catalysis can also occur even at , in the region of the phase diagram involving chemical potential versus for intermediate/high values of field strength.

Consequently, a natural discussion emerges about the combined effects on the phase structure of a thermal gas of quark matter (which might be described as a first approximation by the NJL model) restricted to a reservoir and under the effect of a magnetic background.

Hence, in view of these recent theoretical and experimental advances on the physics of strongly interacting matter, we believe that there is still room to other contributions. To this end, the main interest of this paper is to extend the analysis performed in Ref. Wang:2018qyq . In the present work we will address the questions raised above under the framework of a generalized version of NJL model. We investigate the finite-size effects on the phase structure of three-flavor NJL model with ’t Hooft interaction, which engenders a six-fermion interaction term, without and with the presence of a magnetic background. We make use of mean-field approximation and Schwinger’s proper-time method in a toroidal topology with antiperiodic conditions, manifested by the utilization of generalized Matsubara prescription for the imaginary time and spatial coordinate compactifications. To deal with the nonrenormalizable nature of the NJL model, the ultraviolet cut-off regularization procedure is employed. The behavior of constituent quark masses and , engendered by the solutions of gap equations, are studied under the change of the size of compactified coordinates, temperature and chemical potential . The ’t Hooft interaction contributions are also evaluated.

We organize the paper as follows. In Section II, we calculate the -dependent effective potential and gap equations obtained from the three-flavor NJL model in the mean-field approximation, using Schwinger’s proper-time method generalized Matsubara prescription. The phase structure of the system and the behavior of constituent quark masses and are shown and analyzed in Section III, without and with the presence of a magnetic background. Finally, Section IV presents some concluding remarks.

## Ii The Model

### ii.1 Lagrangian and thermodynamics

Let us start by introducing the -flavor version of the NJL model, whose most commonly used Lagrangian density reads NJL ; NJL1 ; Vogl ; Klevansky ; Hatsuda ; Buballa ,

(1) |

where is the quark field carrying flavors; denotes the mass term,

(2) |

with being the corresponding mass matrix; represents the four-fermion-interaction term,

(3) |

with being the respective coupling constant and the generators of in flavor space; and designates the so-called ’t Hooft interaction,

(4) |

with being the corresponding coupling constant.

We remark some features of the model introduced above. The version of NJL model presented here is generalized to any number of flavors, but we are particularly interested in the case of . This means that the quark field has its components related to the lightest quark flavors: . This allows to identify the generators of in flavor space as follows: are the Gell-Mann matrices and . Also, the corresponding mass matrix becomes . Besides, henceforth the isospin symmetry on the Lagrangian level is assumed, i.e. , while -flavor symmetry is explicitly broken, generating .

From the perspective of -symmetry in Lagrangian terms, we see that is a symmetric four-point interaction. On the other hand, the ’t Hooft interaction is a determinant in flavor space, which means that it is a flavor-mixing -point interaction, involving an incoming and an outgoing quark of each flavor. Thus, for three flavors it engenders a six-fermion-interaction term, and is -symmetric, but breaks the symmetry which was left unbroken by . Thus, express the anomaly, which appears in the QCD context from the gluon sector, in terms of a tree-level interaction in the present pure quark model. As discussed in Refs. Klevansky ; Hatsuda ; Buballa , the ’t Hooft interaction is phenomenologically relevant to yield the mass splitting of and mesons. Moreover, it can be shown that in the chiral limit, , the acquires a finite value of mass due to , whereas the other pseudoscalar mesons stay massless.

The focus of the present analysis concerns the lowest-order estimate of the phase structure of this model. In this sense, the calculations for obtaining relevant quantities will be performed within the mean-field (Hartree) approximation. To this end, the non-vanishing quark condensates,

(5) |

(), are assumed to be the only allowed expectation values which are bilinear in the quark fields. In this approximation, the interaction terms in are linearized in the presence of , that is,

(6) |

which means that terms quadratic in the fluctuations will be neglected. Moreover, terms in channels without condensate or nondiagonal in flavor space, like , and , are are excluded. Therefore, in this context the NJL Lagrangian density can be rewritten as

(7) |

where we have introduced as the constituent quark mass matrix, which is diagonal in flavor space: . Notice that the constant terms in have been neglected, since they give trivial contributions.

At this point we can present the thermodynamic potential density at temperature and quark chemical potential , which is defined by

(8) | |||||

where is the grand canonical partition function, , the Hamiltonian density (the Euclidean version of Lagrangian density ) and the functional trace over all states of the system (spin, flavor, color and momentum). Thus, the integration over fermion fields generates the mean-field thermodynamic potential in the form

(9) | |||||

where is the free Fermi-gas contribution,

(10) |

Here the sum over denotes the sum over the fermionic Matsubara frequencies, .

Then, the minimization of the thermodynamic potential in Eq. (10) with respect to quark condensates allows us to obtain the following gap equations,

(11) |

The physical solutions of Eq. (11) are determined from the stationary points of the thermodynamic potential, which lead to the standard expression for the quark condensates,

(12) |

where and

(13) |

We have assumed for simplicity the quark chemical potentials as . Eq. (11) contains a non-flavor mixing term proportional to the coupling constant (coming from the four-point interaction contribution) and a flavor-mixing term proportional to the coupling constant (coming from the six-point interaction). Hence, the role of the ’t Hooft term is to engender flavor-mixing contributions in the constituent masses.

### ii.2 Generalized Matsubara prescription and proper-time formalism

To take into account finite-size effects on the phase structure of the model, we denote the Euclidean coordinate vectors by , where and , with being the length of the compactified spatial dimensions. As a consequence, the Feynman rules explicited in the arguments of the sum-integral mixing in Eq. (12) must be replaced according to the generalized Matsubara prescription livro ; PR2014 ; Emerson , i.e.,

(14) |

such that

(15) |

where Due to the KMS conditions livro ; PR2014 ; Bellac ; Kapusta , which determines the thermal Matsubara frequencies according to the field statistics, the fermionic nature of the system under study imposes an antiperiodic condition in the imaginary-time coordinate.

However, in the case of spatial compactified coordinates there are no restrictions with respect to the periodicity. So, the choice of the boundary conditions of the fermion fields in the spatial directions results in a relevant issue, especially because it affects the meson masses and other relevant observables, as claimed in studies mentioned in Introduction. Current lattice QCD simulations usually employ periodic quark boundary conditions Klein:2017shl , thanks to its tendency to reduce the finite volume effects, favoring the obtention of the thermodynamic limit. Nevertheless, as has been supported by a large number of studies that make use of effective QCD models Gasser:1986vb ; Klein:2017shl ; Shi:2018swj ; Wang:2018qyq , the fields should take the same boundary condition in their spatial and temporal directions. One important aftermath of this selection is the symmetry in the formalism between the spatial and temporal directions, through the exchange , causing -independent coupling constants and . So, we can fit the model parameters at zero temperature and infinite volume and hold them in the finite temperature and volume analysis. Thus, following this assumption, we also adopt antiperiodic boundary conditions for the compactified spatial coordinates: in Eq. (15) assumes the value .

In the present work the thermodynamic potential and the gap equations will be treated within the Schwinger proper-time method Schwinger ; Farina . Accordingly, the kernel of the propagator in Eq. (12) can be rewritten as

(16) |

where is the so-called proper time. Therefore, the application of Eq. (16) and the generalized Matsubara prescription (15) into (12) yields the quark chiral condensates expressed as

(17) | |||||

In addition, to perform the necessary manipulations in a relatively simple and more tractable way, it is convenient to employ the Jacobi theta functions Bellman . Noticing the very useful property of the second Jacobi theta function, ,

(18) | |||||

then the -dependent chiral quark condensate in Eq. (17) can be expressed in the following way:

(19) | |||||

It is worthy noticing that in order to regularize the proper-time approach, we introduce a ultraviolet cutoff , namely,

(20) |

We conclude this section by remarking that the bulk form of the system can be studied straightforwardly. To do this, the continuum limit is performed by reverting the Matsubara prescription in Eq. (14), which yields in Eq. (19) Gaussian integrals in momentum space. So, the expression for the chiral quark condensate acquires the form

Hence, we have obtained above the thermodynamic potential and gap equations whose solutions produce the -dependent constituent quark masses. Other thermodynamic quantities such as pressure, entropy and others can be obtained by similar procedures to those described above.

In the next section, we will discuss the thermodynamic behavior of the present model.

### ii.3 Presence of a magnetic background

Considering now the system in the presence of an uniform external magnetic field, we redefine the Lagrangian density in Eq. (1) with the derivative in kinetic term replaced by the modified covariant derivative: , where denotes the four-potential associated to the magnetic background, and is the quark charge electric matrix, , with . We choose for convenience the Landau gauge, , where is the intensity of the external magnetic field in the direction . Therefore, the field operators are given by the set of normalized eigenfunctions of the Landau basis, which means that the energy eingenvalues associated to the solutions of equation of motion read,

(22) |

where is the so-called cyclotron frequency; is the spin variable; and , denotes the so-called Landau levels.

A fundamental consequence of presence of magnetic background is the modification of Feynman rules: the four-momentum integrals suffer the dimensional reduction, i.e.

(23) |

When combined with finite-temperature and finite-size effects, the Euclidean coordinate vectors must be rewritten with and , being the size of the compactified spatial dimension. Then, taking the Matsubara prescription discussed in Eq. (14) the integration over momenta space is modified accordingly,

(24) |

Hence, the thermodynamic potential and gap equations given by Eqs. (10) and (11) must be analyzed keeping in mind the thermodynamic relations and chiral condensates with the modified integration over momenta as in Eq. (24). This means that the expressions for in Eq. (19) should be replaced by

(25) | |||||

which can be rewritten, after performing the proper time regularization procedure and the sum of geometrical series and spin polarization, as

(26) | |||||

In the next section, we will discuss the -dependence of thermodynamic quantities introduced above.

## Iii Phase structure

We devote this section to the analysis of the phase structure of the system, focusing on how it behaves with the change of the relevant parameters of the model and, in special, the influence of the boundaries on the behavior of constituent quark masses and , which are solutions of expressions given by Eq. (11); explicitly, they are

(27) |

with We simplify the present study by fixing , which means that the system consists in a -quark gas constrained in a cubic box.

### iii.1 Absence of magnetic field

Once the regularization procedure is chosen, one must determine the complete set of input parameters that provides a satisfactory description of hadron properties at zero temperature and density. In this sense, we use the model parameters reported in Ref. Kohyama:2016fif , which have been estimated by the fitting in light of the following observables: MeV, MeV, MeV and MeV. Our choice of the set is the following:

(28) |

The contribution of the ’t Hooft interaction will be always taken into account unless explicitly stated otherwise, as in the case of Fig. 2, where the parameters have been refitted in order to obtain the same dressed masses at zero temperature and density that correctly reproduce the relevant observables. Moreover, since we work in the limit of isospin symmetry for current and quark masses, in absence of magnetic field the solutions of Eq. (27) satisfy . Hence, throughout this subsection we will show and discuss the results keeping this fact in mind.

We start, for completeness, by studying the behavior of the constituent quark masses under a change of parameters but without the presence of boundaries, which is basically the scenario described in Refs. Buballa ; Franceses . In Fig. 1 we have plotted the values of (top panel) and (bottom panel) that are solutions of the gap equations in Eq. (27) as functions of the temperature , in the bulk limit and taking different values of the chemical potential . It can be seen that the set of parameters given by Eq. (28) engenders the constituent quark masses and at vanishing temperature. In the range of smaller temperatures, there is no relevant modification. This phenomenon remains up to a certain temperature, where the masses start to decrease with the augmentation of . Therefore the broken phase is inhibited and a crossover transition takes place. Above a given temperature value, the dressed quark masses approach the magnitudes of the corresponding current quark masses, i.e. . It is also worth mentioning that the constituent mass for the -quark has a smoother drop than the one for the -quark, falling to the magnitude at larger . Another feature is that at higher temperatures the system tends faster toward the symmetric chiral phase as the chemical potential increases.

We now investigate in detail the influence of boundaries on the phase structure of the system, considering the -quark gas in a region delimited by a cubic box. In Fig. 2 we have plotted the values of that are solutions of the gap equations in Eq. (27) as function of the inverse length , at vanishing temperature and chemical potential. For the sake of comparison with the existing literature (e.g., Refs. Buballa ; Franceses ; Wang:2018qyq ), we have taken into account the situations with and without the ’t Hooft interaction term (i.e., or , respectively). Note that in the case of a vanishing , the parameters in Eq (28) have been readjusted keeping the values of and fixed and varying and to obtain the same dressed masses estimated at zero temperature and density for the situation with . Then, in this context we have worked with the set,

(29) |

We can observe that the bulk approach appears as a good approximation in the range of greater values of (up to fm). The dressed masses are affected by the presence of boundaries as the size of the system decreases, acquiring smaller values smoothly. In this region of lower , we see that the presence of the ’t Hooft interaction, which engenders a system of coupled gap equations shown in Eq. (27), keeps the essential features unchanged when the parameter set is properly chosen to reproduce the masses in the bulk. Only in the case of the contributions coming from yield a different decreasing more pronounced than the situation of a vanishing . Also, we notice that below a given value of , the dressed quark masses tend to the corresponding current quark masses. In particular, this happens at - fm for -quarks. So, analogously to the behavior described in the previous figures, has a smoother falloff to the current quark mass than . Thus, the dependence of the phase structure of the system on the inverse length is qualitatively similar to the one found for the temperature, which is expected, due to the equivalent nature between and , manifested in the generalized Matsubara prescription in Eq (15).

The plot in Fig. 3 is the same as in Fig. 2 for the case , but with the solutions for and obtained at vanishing temperature and chemical potential taking the chiral limit for the quarks, i.e., setting . Once more, as in the previous situation, in the chiral limit we have refitted the parameters to give a good prescription of the hadronic observables at zero temperature and density and in the bulk. Therefore, keeping the value of and fixed and varying and to obtain the same dressed masses estimated at , we have employed the set

(30) |

Again, the restoration of chiral symmetry is induced at smaller sizes of the system, lost in the bulk. In this context, the results suggest the existence of a critical volume at which the constituent mass for -quark reaches zero value via a second-order phase transition. On the other hand, goes down to the value of softly.

Let us concentrate at this point on the evaluation of the combined effects of boundaries, finite temperature and finite chemical potential on the system. They can be better described from Fig. 4, in which are plotted the values of effective masses that are solutions of the gap equation in Eq. (27) as function of inverse length for different values of , keeping the temperature fixed at 150 MeV. As suggested in previous figures, in the bulk ( or ), and suffer a decrease as the chemical potential grows, and diminishes even more with the drop of and tend to the respective current quark masses. Also, at smaller the dressed masses associated to different values of converge to the values of corresponding current quark masses.

This analysis is completed with Figs. 5 and 6, where are plotted the values of and as functions of the temperature , taking different values of chemical potential and size , respectively. When compared with Fig. 1, the plots in Fig. 5 manifest the smaller values of constituent quark masses reached in the case of finite volume. Fig. 6 shows the values of the size of cubic box in which and stand with the corresponding values of current quarks masses in all range of temperature. Hence, our findings suggest that the presence of boundaries disfavors the maintenance of long-range correlations, inducing the inhibition of the broken phase.

### iii.2 Magnetic background effects

Now we investigate the combined effects of finite volume and magnetic background on the phase structure of the thermal gas of quark matter, taking into account that here we have only one spatial compactified coordinate. We remember that in gap equations (27) of present situation we use a different expression for the condensates, i.e. Eq. (26) rather than (19), in which a dimensional reduction of four-momentum integrals in Feynman rules took place. The parameter set used is shown in Eq. (28), which gives the correct hadronic observables at .

Furthermore, noticing that the coupling of the quarks to the electromagnetic field depends on the quark flavor (see discussion in subsection II.3), then in principle we cannot adopt the solutions of Eq. (27) satisfying , as done previously. That being so, we solve the system of three gap equations in (27) to assess the differences between and engendered by the coupling to magnetic background.

In Fig. 7 is plotted the values of constituent quark masses as functions of inverse of thickness for different values of , keeping temperature and chemical potential fixed. For higher values of the considered range of magnetic field strength, we remarked that the approximate isospin symmetry for the constituent masses of the and quarks is broken, with being bigger than due to the larger coupling of quark to magnetic background. Also, we see that in the bulk ( or ), the presence of an external magnetic field engenders the augmentation of . In other words, the raise of field strength enhances the broken phase. This feature is expected and called in literature as magnetic catalysis effect. This phenomenon remains in all considered range of . On the other hand, as previously the reduction of engenders a reduction of constituent masses, but in a more slightly way than the situation without magnetic field. Also, in the range of smaller values of the dressed masses associated to different values of converge to the values of corresponding current quark masses. The interesting point is that the growth of field strength induces smaller values for at which the system stand with the values . Thus, the combination of finite size and magnetic effects generates a competition between then, since the former inhibits the broken phase whereas the last one yields its stimulation.

We complement this analysis with Figs. 8 and 9, where are plotted the values of as functions of the temperature , taking different values of field strength and size , respectively, which can be compared with Figs. 5 and 6. Again, they manifest the smaller values of constituent quark masses reached in the case of finite size and greater temperatures. Since the magnetic field increases the constituent quark masses, the critical size at which stand with the corresponding values of current quarks masses in all range of temperature should be even smaller than the situation without magnetic field (0.5-1.0 fm).

We conclude this section by discussing some issues. We notice that the relevance of the regularization procedure and parametrization in effective models can not be underestimated. The findings obtained above are naturally dependent on the set of parameters considered as input in Eq. (27). Modifications in these choices will obviously alter the magnitudes of the constituent quark masses and the ranges of which engender changes on their behavior, making them drop faster or slower and accordingly reshape the curves plotted in previous figures. Thus, since in the present approach we have chosen the set of input parameters that provides a satisfactory description of hadron properties at (according to Ref. Kohyama:2016fif ), the comparison with existing literature only makes sense with those works that have proceeded in the same way.

Notwithstanding the points raised above, we stress the main results of this work as follows. In the end, we see that the phase structure of the system is strongly affected by the combined variation of relevant variables, depending on the competition among their respective effects. At smaller temperatures, with the system in the broken phase, the bulk approach seems a good approximation in the range of greater values of the size of the system , since the constituent masses do not suffer modification. Nevertheless, keeping fixed, the reduction of from a given value engenders a reduction of via a crossover phase transition to their corresponding current quark masses. Yet at the same temperature and at the range , the broken phase is disfavored due to both increasing of chemical potential and the drop of . On the other hand, when a magnetic background is present the broken phase is stimulated, acquires greater values and the increase of field strength induces even smaller values for to reach symmetric phase.

Finally, we should point out the combined effect of an external magnetic field and the choice of boundary conditions in this model. The boundary effect in the context of antiperiodic condition has been studied in detail and summarized in previous paragraph. Then, we deserve some comments to the contribution of the compactified dimension length in the situation of periodic boundary conditions. After the compactification of the -coordinate in the energy eigenvalues in Eq. (22), it can be seen that for the periodic condition, we have the prescription , and no energy gap is yielded between the vacuum and the lowest Landau level at zero momentum in the chiral limit. This makes possible the vacuum condensation of fermion pairs in the infrared range of spectrum (see Ref. Ferrer:1999gs for a detailed discussion). In this sense, the presence of magnetic field catalyzes the dynamical breaking of chiral symmetry. Another perspective to see this is looking at the chiral condensate in Eq. (26) and gap equations in (27), in which the Jacobi -function must be replaced by the in the scenario of periodic condition. Thanks to the different behaviors of the -functions, the consequence of this replacement is relevant: while in the antiperiodic case the magnetic catalysis is inhibited, with the length playing a role similar to the inverse of temperature , in the periodic situation the boundary has the opposite effect: symmetry breaking is enhanced, and accordingly the constituent quark masses acquire higher values as diminishes. Thus, there is not a critical value of the size in which the symmetry is restored, and the combined effect of magnetic background and boundary conditions in periodic case strengthens the broken phase.

We must emphasize, however, that in our understanding it seems more appropriate to adopt the same boundary condition in both spatial and temporal directions, as remarked in subsection II.2 and supported by a large number of studies using effective models Gasser:1986vb ; Klein:2017shl ; Shi:2018swj ; Wang:2018qyq .

## Iv Concluding Remarks

In this work we have analyzed the finite volume and magnetic effects on the phase structure of a generalized version of the Nambu–Jona-Lasinio model with three quark flavors. By making use of mean-field approximation and the Schwinger proper time method in a toroidal topology with antiperiodic conditions, we have investigated the gap equation solutions under the change of the size of compactified coordinates, temperature, chemical potential and strength of external magnetic field. The thermodynamic behavior is strongly affected by the combined effects of relevant variables, depending on the range of their change.

Looking at the influence of a finite volume, we can remark the main results of this work: at smaller temperatures, with the system in the broken phase, the bulk approach seems a good approximation in the range of greater values of the volume of the system , since the constituent masses are not modified. Nevertheless, while keeping fixed, there is a range in which the reduction of (remarkably at fm) engenders a reduction of the constituent masses for quarks through a crossover phase transition to the their corresponding current quark masses. After all, the broken phase is disfavored due to both increasing of temperature and chemical potential, and the drop of . However, in the region in which the broken phase is inhibited, the current quark masses are not affected by the change of variables.

Looking at the magnetic background influence, only phenomena associated to the magnetic catalysis take place: an enhancement of broken phase occurs in all considered ranges of temperature and size. When a magnetic background is present, constituent quark masses acquire greater values and the increase of field strength induces smaller values for and greater values of at which the tend to the respective values of current quark masses. For higher values of magnetic field strength, we remarked that the approximate isospin symmetry for the constituent masses of the and quarks is no longer valid, with being bigger than due to the larger coupling of quark to magnetic background.

Finally, we emphasize that the results outlined above can give us insights about the finite-volume and magnetic effects on the critical behavior of quark matter produced in heavy-ion collisions. In particular, estimations have been done on the range of size of the system at which the bulk approximation looks a reasonable one. Further studies are necessary in order to verify the efficacy of the proposed framework.

###### Acknowledgements.

This paper is dedicated to the memory of Edson B. S. Corrêa. The authors would like to thank the Brazilian funding agencies CNPq, CAPES. L.M.A. would like to thank the Brazilian funding agencies CNPq (contracts 308088/2017-4 and 400546/2016-7) and FAPESB (contract INT0007/2016) for financial support. E.B.S.C. would like to thank E. Cavalcanti for helpful discussions.## References

- (1) J. Adams et al. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005).
- (2) P. Braun-Munzinger, V. Koch, T. Schafer, and J. Stachel, Phys. Rep. 621, 76 (2016).
- (3) F. Prino and R. Rapp, J. Phys. G 43, no. 9, 093002 (2016).
- (4) R. Pasechnik and M. Šumbera, Universe 3, no. 1, 7 (2017).
- (5) S. A. Bass et al., Prog. Part. Nucl. Phys. 42, 313 (1999).
- (6) L. F. Palhares, E. S. Fraga and T. Kodama, J. Phys. G 38, 085101 (2011).
- (7) G. Graf, M. Bleicher and Q. Li, Phys. Rev. C 85, 044901 (2012).
- (8) J. Luecker, C. S. Fischer and R. Williams, Phys. Rev. D 81, 094005 (2010).
- (9) B. L. Li, Z. F. Cui, B. W. Zhou, S. An, L. P. Zhang and H. S. Zong, Nucl. Phys. B 938, 298 (2019).
- (10) C. Shi, W. Jia, A. Sun, L. Zhang and H. Zong, Chin. Phys. C 42, no. 2, 023101 (2018). doi:10.1088/1674-1137/42/2/023101
- (11) J. Braun, B. Klein and H.-J. Pirner, Phys. Rev. D 71, 014032 (2005) doi:10.1103/PhysRevD.71.014032 [hep-ph/0408116].
- (12) J. Braun, B. Klein, H.-J. Pirner and A. H. Rezaeian, Phys. Rev. D 73, 074010 (2006) doi:10.1103/PhysRevD.73.074010 [hep-ph/0512274].
- (13) E. J. Ferrer, V. P. Gusynin and V. de la Incera, Phys. Lett. B 455, 217 (1999) doi:10.1016/S0370-2693(99)00470-0 [hep-ph/9901446].
- (14) L. M. Abreu, M. Gomes and A. J. da Silva, Phys. Lett. B 642, 551 (2006).
- (15) D. Ebert, K. G. Klimenko, A. V. Tyukov and V. Ch. Zhukovsky, Phys. Rev. D 78, 045008 (2008).
- (16) L. M. Abreu, A. P. C. Malbouisson, J. M. C. Malbouisson and A. E. Santana, Nucl. Phys. B 819, 127 (2009).
- (17) L. M. Abreu, A. P. C. Malbouisson and J. M. C. Malbouisson, Phys. Rev. D 83, 025001 (2011).
- (18) A. Bhattacharyya, P. Deb, S. K. Ghosh, R. Ray and S. Sur, Phys. Rev. D 87, no. 5, 054009 (2013).
- (19) A. Bhattacharyya, R. Ray and S. Sur, Phys. Rev. D 91, no. 5, 051501 (2015).
- (20) A. Bhattacharyya, R. Ray, S. Samanta and S. Sur, Phys. Rev. C 91, 041901 (2015).
- (21) H. Kohyama, D. Kimura and T. Inagaki, Nucl. Phys. B 906, 524 (2016)
- (22) Z. Pan, Z. F. Cui, C. H. Chang and H. S. Zong, Int. J. Mod. Phys. A 32, no. 13, 1750067 (2017).
- (23) Q. Wang, Y. Xiq and H. Zong, Mod. Phys. Lett. A 33, no. 39, 1850232 (2018)
- (24) J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987).
- (25) P. H. Damgaard and H. Fukaya, JHEP 0901, 052 (2009)
- (26) E. S. Fraga and L. F. Palhares, Phys. Rev. D 86, 016008 (2012).
- (27) L. M. Abreu, C. A. Linhares, A. P. C. Malbouisson, J. M. C. Malbouisson, Phys. Rev. D 88, 107701 (2013).
- (28) L. M. Abreu, A. P. C. Malbouisson, J. M. C. Malbouisson, E. S. Nery and R. Rodrigues da Silva, Nucl. Phys. B 881, 327-342 (2014).
- (29) D. Ebert, T. G. Khunjua, K. G. Klimenko and V. C. Zhukovsky, Phys. Rev. D 91, 105024 (2015).
- (30) L. M. Abreu, E. S. Nery and A. P. C. Malbouisson, Phys. Rev. D 91, 087701 (2015).
- (31) L. M. Abreu and E. S. Nery, Int. J. Mod. Phys. A 31, 1650128 (2016).
- (32) L. M. Abreu, A. P. C. Malbouisson and E. S. Nery, Mod. Phys. Lett. A 31, 1650121 (2016).
- (33) S.S. Bao and H. Shen, Phys. Rev. C 93, 025807 (2016).
- (34) L. M. Abreu, and E. S. Nery, Phys. Rev. C 96, 055204 (2017).
- (35) S. Samanta, S. Ghosh and B. Mohanty, J. Phys. G: Nucl. Part. Phys. 45, 075101 (2018).
- (36) X.H. Wu and H. Shen, Phys. Rev. C 96, 025802 (2017).
- (37) B. Klein, Phys. Rept. 707-708, 1 (2017).
- (38) C. Shi, Y. Xia, W. Jia et al., Sci. China Phys. Mech. Astron. 61 082021 (2018).
- (39) D.E Kharzeev, L.D. Mclerran and H.J. Warringa, ArXiv:0711.0950 [hep-ph], 2007.
- (40) V. Skokov, A. Y. Illarionov and V. Toneev, Int. J. Mod. Phys. A 24, 5925 (2009).
- (41) M. N. Chernodub, Phys. Rev. D 82, 085011 (2010).
- (42) A. Ayala, M. Loewe, J.C. Rojas and C. Villavicencio, Phys. Rev. D 86, 076006 (2012).
- (43) M. Ferreira, P. Costa, O. Lourenco, T. Frederico, and C. Providencia, Phys. Rev. D 89, 116011 (2014).
- (44) A. Haber, F. Preis and A. Schmitt, Phys. Rev. D 90, 125036 (2014).
- (45) S. Mao, Phys. Lett. B 758, 195 (2016).
- (46) A. Ayala, P. Mercado and C. Villavicencio, Phys. Rev. C 95, 014904 (2017).
- (47) K. A. Mamo, JHEP 1505, 121 (2015)
- (48) V. P. Pagura, D. Gomez Dumm, S. Noguera, and N. N. Scoccola Phys. Rev. D 95, 034013 (2017).
- (49) N. Magdy, M. Csanad and R. A. Lacey, J. Phys. G: Nucl. Part. Phys. 44, 025101 (2017).
- (50) A. Ayala, C. A. Dominguez, S. Hernandez-Ortiz, L. A. Hernandez, M. Loewe, D. M. Paret, and R. Zamora, Phys. Rev. D 98, 031501(R) (2018).
- (51) Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961).
- (52) Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 (1961).
- (53) U. Vogl and W. Weise, Prog. Part. Nucl. Phys. 27, 195 (1991).
- (54) S. P. Klevansky, Rev. Mod. Phys. 27, 195 (1991).
- (55) T. Hatsuda and T. Kunikiro, Phys. Rep. 247, 221 (1994).
- (56) M. Buballa, Phys. Rep. 407, 205 (2005).
- (57) F.C. Khanna, A.P.C. Malbouisson, J.M.C. Malbouisson, and A.E. Santana, Thermal Quantum Field Theory: Algebraic Aspects and Applications, World Scientific, Singapore (2009).
- (58) F.C. Khanna, A.P.C. Malbouisson, J.M.C. Malbouisson, and A.E. Santana,, Phys. Rep. 539, 135 (2014).
- (59) E. B. S. Corrêa, C. A. Linhares, A. P. C. Malbouisson, J. M. C. Malbouisson, and A. E. Santana, Eur. Phys. J. C, 77, 261 (2017).
- (60) M. Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge, UK, 1996).
- (61) J.I. Kapusta and C. Gale, Finite-Temperature Field Theory: Principles and Applications (Cambridge University Press, Cambridge, UK, 2006).
- (62) J. Schwinger, Phys. Rev. 82, 914 (1951).
- (63) F.A. Barone, H. Boschi-Filho, C. Farina, Am. J. Phys. 71, 483 (2003).
- (64) R. Bellman, A Brief Introduction to Theta Functions, Holt, Rinehart and Winston, Inc., New York (1961).
- (65) F. Gastineau, R. Nebauer and J. Aichelin, Phys. Rev. C 65, 045204 (2002).