An NJL model analysis of a magnetised nonextensive QCD medium

We investigate the effect of a background magnetic field when applied to a nonextensive QCD medium in a $2+1$ flavour Nambu\textemdash Jona-Lasinio model. The effect of a constant as well as an $eB$-dependent coupling is considered. For the constant coupling, the well-known magnetic catalysis effect is observed with reduced strength in both the condensates and the transition temperatures compared to the standard extensive medium. In the case of a field-dependent coupling, we observe a competition between the nonextensive parameter, $q$ and $eB$. With sufficiently high $q$-values, the phenomenon of inverse magnetic catalysis \textemdash\, a well-known trait in effective models featuring $eB$-dependent coupling, appears to be eliminated within the range of magnetic fields examined in our present study.


I. INTRODUCTION
The creation of a strong magnetic field in heavy-ion collisions (HIC) is now an established fact.Its strength in non-central collisions can be huge as compared to any known terrestrial or extraterrestrial magnetic field-generating phenomenon.Depending on the collisional energy the strength of the magnetic field can vary between ∼ 1 − 15 m 2 π from Relativistic Heavy Ion Collider (RHIC) to Large Hadron Collider (LHC) [1].
From the later half of the first decade in this century, the effect of the magnetic field on the properties of the QCD medium have been explored extensively using both theoretical and experimental techniques [2,3].Numerous novel phenomena have been proposed and studied in that context, some of the most studied ones are magnetic catalysis (MC) [4], inverse magnetic catalysis (IMC) [5,6], chiral magnetic effect [7], chiral separation effect [8,9] etc.Those interested can delve into an engaging review article that examines a spectrum of potential novel non-dissipative phenomena in the presence of magnetic fields and beyond [10].
On the other hand, in the pursuit of understanding the data from HIC experiments both in RHIC and LHC one statistical distribution, recently, has gained a lot of attention [11][12][13][14][15].This is widely known as the Tsallis distribution after Constantino Tsallis who in his famous paper [16] investigated the possibility of a nonextensive form of entropy involving a parameter q.The extensive nature of the entropy is obtained back in the limit q → 1.As compared to the exponentials in traditional Boltzmann-Gibbs (BG) extensive statistics, the nonextensive statistics lead to power laws.Its application is found in almost every branch of physics [17].
Its applicability to describe the fireball created in HICs stems from the following facts.The fireball is finite, rapidly evolving in nature and also experiences strong intrinsic fluctuations and long-range correlations.These characteristics render the system non-uniform and global equilibrium cannot be established.As a consequence, some quantities develop power law tailed behaviour and the nonextensive statistics becomes applicable.
There are multiple versions of the Tsallis distribution available in the literature [16,[18][19][20].Theoretical issues concerning the thermodynamic consistency while implementing Tsallis distribution in the case of relativistic high energy distribution [21,22] have been clarified in Ref. [23,24].A recent review on the topic concerning nuclear and particle physics can be found in Ref. [25].
There have also been multiple developments in applying nonextensive statistics in QCD-based effective model describing both the hadronic and quark degrees of freedoms [26][27][28][29][30][31].These studies, among other things, discussed scalar and vector meson fields, properties of π and σ mesons, chiral symmetry restoration, effect of background gauge field (Polyakov loop field), characteristics of magnetisation, thermodynamic properties of QCD medium etc.
Existing literatures on utilising the effective models for a nonextensive QCD medium do not incorporate the background magnetic field [27][28][29]31].On the other hand, there have been ample efforts to understand the effects of MC, IMC and magnetic properties like magnetisation, magnetic susceptibility, etc., of the QCD medium in different effective models [32][33][34][35][36][37][38][39], albeit using extensive statistics.Thus, there is currently a gap in the existing literature regarding effective models for a nonextensive QCD medium that is also magnetized.
Keeping this in mind, in the present work, we shall present a nonextensive version of the Nambu-Jona-Lasinio (NJL) model describing a magnetised QCD medium.Specifically, chiral symmetry restoration while intertwined with nonextensiveness and magnetic field will occupy the central stage of our investigation.Such a study, besides addressing the existing theoretical gap, could have relevance for understanding magnetized QCD media with strong intrinsic fluctuations, long-range interactions, etc., potentially encountered in heavy-ion collision experiments.We start our investigation for the nonextensive statistics with a zero magnetic field (eB) for 2 + 1 flavours to put things into perspective.Then we move to the case of nonzero strength of eB.To gauge the effect of the nonextensive statistics we execute the analysis for BG statistics in parallel.We perform the analysis in presence of eB considering different possible scenarios.
First, we check results for the standard NJL model where the coupling constant is independent of eB.Our result shows that at a specific magnetic field strength, the impact of increasing q resembles that of the zero magnetic field case.It leads to a decreasing trend in both the condensates and the transition temperatures.We find out that the model produces the MC effect in Tsallis statistics as well.However, it results in reduced strength in both the condensates and the transition temperatures compared to the BG statistics.
Then, we take our analysis to the case of a magnetic field dependent coupling constant.We choose to work with a coupling constant whose dependence on eB has been determined by fitting to the available lattice QCD (LQCD) data, which produces the IMC effect [33].We discover that there is a competition between the two parameters q and eB.This is one of the most important findings of the present analysis.Depending on their strengths the observed IMC effect in BG statistics can be eliminated.
Our findings also suggest that with eB-dependent G S the trend of the transition temperature is non-monotonic.For smaller values of q, such as q = 1.05, it remains constant initially and then starts decreasing with increasing magnetic field.With further increment in q, it exhibits a non-monotonic trend as a function of eB: initially it increases and then starts falling.To render the present analysis complete, we also explore the basic form of the coupling constant inspired by Ref. [40], in which the strong coupling α s decreases with increasing eB and put it in the appendix.
The remainder of the manuscript is organised as follows: First, we provide a brief review on the formalism for 2 + 1 flavour NJL model in section II.We start with the standard NJL model in BG statistics in subsection II A and then extend it to the nonextensive statistics in subsection II B. In section III, we thoroughly discuss our findings and compare them with the known results for detailed analysis and comparison.Finally, in section IV, we summarise the conclusion.

A. Extensive statistics
In this section, we briefly revisit the NJL model in the standard extensive scenario.The 2+1 flavour NJL Lagrangian in presence of a magnetic field can be written as [32,35,[41][42][43][44][45][46][47], where and q being the electric charge (q u = 2/3e, q d = −1/3e and q s = −1/3e); e is the charge of a proton, F µν = ∂ µ A ν − ∂ ν A µ is the field strength tensor and A µ being the electromagnetic gauge field.The Lagrangians L 1 and L 2 are [48][49][50] The Lagrangian L 1 is usually called the symmetric Lagrangian as it respects the full global symmetry of QCD, SU On the other hand, L 2 breaks the U (1) A symmetry and is responsible for the axial anomaly found in QCD.It is known as the 't Hooft determinant term or simply the determinant term.Before obtaining the thermodynamic potential, one must note a few important things.Because of the presence of the magnetic field there are some important modifications in the calculation.First, the dispersion relation of the quarks will be modified as, where M f is the effective quark mass for a quark of flavour f .The magnetic field is in the z direction ( ⃗ B = B ẑ), which results in discretising the momentum in the transverse plane, here in the x − y plane.The transverse component of the momentum now depends on both the Landau level (LL), l and the spin, s.It is important to note how in the magnetised version the energy of the individual flavour is dependent not only on the respective effective masses but also on the respective charges.Another important change, which is somewhat expected from Eq. 4, is in the integral over the three momenta, The integral is now performed only over the continuous component, p z .Also, note the change in the lower limit of the integral as compared to the zero magnetic field case, where it is 0. With all these knowledge in hand, we can get the thermodynamic potential potential for the NJL model in presence of the magnetic field as, where, σ f = ⟨ ψf ψ f ⟩ with f = u, d and s 1 .The first two terms are the vacuum contributions from the condensates of different flavours.The second term arises from the 't Hooft determinant term and is the explicit manifestation of preservation of the axial anomaly.Ω f mev is a term which we call magnetic field entangled vacuum term.It can be separated into a purely magnetic field dependent term and the standard vacuum term, which we show in a moment.Ω f med is the term in presence of both the magnetic field and the temperature and can be termed as the medium term.These two terms are given as, respectively.f , l and s stand for the summation over flavour, Landau levels and spin states of the quarks, respectively.Following the technique in Ref. [41], Ω f mev can be further decomposed into the standard vacuum term (Ω f vac ) and a purely magnetic field dependent term ( where ) is the Riemann-Hurwitz zeta function.The detail of obtaining the expression for Ω mag can be found in Ref. [41].
Thus, the final expression for the thermodynamic potential in NJL model becomes, We take the model parameters as: Number of light flavours = 2, number of colours = 3, m u = m d = 0.0055 GeV, m s = 0.1407 GeV, three momentum cut-off Λ = 0.6023 GeV, the scalar coupling constant3 , G S = 3.67/Λ2 , the coupling constant for the 't Hooft determinant term, G D = 12.36/Λ 5 , pion mass m π = 0.135 GeV and the pion decay constant, f π = 0.0924 GeV (0.0879 GeV, in the chiral limit), kaon mass m k = 0.4977 GeV, eta mesons, m ′ η = 0.9578 GeV and m η = 0.5148 GeV.All these parameters are taken from the Ref. [51].
The magnetic field dependent coupling constant is fitted as [33], where, Because of the 't Hooft determinant term there is a cross-flavour coupling which is called "flavour mixing" [50].As a result the effective masses are given by for up, down and strange quarks, respectively.

B. Nonextensive statistics
We know that the entropy is considered to be an extensive property, i.e., its value depends on the system size.Constantino Tsallis proposed a nonextensive statistics [16], which provides a generalisation of the traditional BG entropy.This generalisation is made through a parameter q, which, in the limit q → 1, gives back the BG entropy.
As compared to the exponentials in traditional extensive statistics, the nonextensive statistics lead to power laws.This is the most fundamentally different manifestation of the entropy being a nonextensive quantity.
Mathematically speaking, in nonextensive statistics or the so-called Tsallis statistics, the usual exponential are replaced by an equivalent q-exponential, where, Using these expressions one can get the generalised versions of the standard statistical distributions.One can derive the Tsallis distribution starting from the Boltzmann equation [52].It is straightforward to check that the above equation reduces to the standard exponential in the limit q → 1 giving back the standard distributions.Any values of q other than unity gives a measure of the "nonextensiveness" of the system.In the same fashion the inverse function, i.e., the equivalent to standard logarithm can as well be defined, ln q (x) ≡ It can be again trivially checked that this reduces to the standard logarithmic function of the variable x in the limit q → 1.Now, we will introduce the NJL model for a nonextensive statistics.In that case, only Ω f med will be modified from the expression of the full thermodynamic potential given in Eq. 11.Ω f med is given in Eq. 8.It is evident that the logarithms in that equation will be modified to incorporate the nonextensive parameter q.This will be done using Eqs.15 and 16.
For the analysis in the nonextensive scenarios, we assume that there will not be any modifications to the model parameters due to the nonextensive nature of the system.So, the treatment of the nonextensive parameter q will be in the same footing as the temperature and chemical potential in the model.The model parameters fitted at zero temperature and chemical potential are used to evolve the system at nonzero temperature and chemical potential.The guideline for q remains the same.We extend this argument further in the presence of the magnetic field and the treatment of the magnetised medium remains analogous.
The modification in Eq. 8 can be summarised as, In the present analysis, we stick to the zero chemical potential scenario. 4In that case, ensuring e q (x) to be always a real number the condition 1 + (1 − q)x 1/(1−q) ≥ 0 has to be satisfied.Considering q > 1 always satisfies the condition.On the other hand, q-values > 1 are also phenomenologically well motivated [12,24], which we elaborate in the result section.

A. In the absence of a magnetic field
To explore the properties of a nonextensive magnetised medium we begin by revisiting the scenario in the absence of an external field.To meet that purpose we employ a 2 + 1-flavour NJL model.We are more interested to explore the QCD phase diagram.Thus, chiral dynamics will hold the central attention while we utilise the NJL model.In the upper panel of Fig. 1, we show the plots of condensates for up/down (left panel) and strange (right panel) quarks as a function of T for different values of the nonextensive parameter q.We observe that the chiral condensates decreases with increasing nonextensivity around the transition region.The transition also tends to be sharper for both the light and heavy quarks.Eventually, we will be more focussed on the fate of the light quarks.The plots of the strange quark condensates are displayed to keep track of the heaviest constituent of the system.To obtain the plots for the corresponding effective masses is straightforward (Eq.13) and they reflect the same decreasing behaviour as the quark condensates with increasing q values as shown in the lower panel of Fig. 1.Such weakening of chiral condensates and effective masses when the medium departs from its extensive nature agree with existing results.[27,29].
For the choice of q > 1, in principle, one can take any higher values but in practice, it is never far away from unity, particularly in the case of HICs.In such experiments, typical values of q obtained from the fits to the transverse momentum distribution for identified charged particles lies within the range of 1.1 − 1.2 [12,24].Thus, the maximum deviation of q should not exceed 20%.Keeping this in mind we have taken three different values of q, 1.05, 1.10 and1.15 to illustrate the effect of "nonextensiveness"; with the highest value remaining well within the expected limit.
In Table I, we quote the transition temperatures (T CO ) as a function of q.They are obtained by tracking the inflection points of the condensates.The quoted numbers are only for the light quarks (same for both up and down quarks).To have a pictorial idea of the decrease in chiral transition temperature, we show a plot of the T CO as a function of q with the T CO being scaled with T CO (q = 1) in Fig. 2.There is roughly 10 percentage decrease in the transition temperature for the highest value of q that we consider here.B. In the presence of a magnetic field In the pursuit of exploring the magnetised nonextensive medium, we proceed chronologically, i.e., we follow the steps how the NJL model in presence of the magnetic field has been developed over the years.Thus, we start with the constant G S .In that case, the NJL model is known to produce the effect of MC throughout the temperature range [53,54].It is interesting to check how the magnetised nonextensive medium behaves with a constant coupling.
After the discovery of the IMC effect by LQCD [5], effective models were tweaked to reproduce that effect.One of the ways in the NJL model is to make the coupling constant dependent on eB [33] or on both eB and T [32,35].In our present study, we choose to work with the parametrisation given in [33].The parametric form solely depends on eB and there is no involvement of temperature.Since we will discuss (I)MC effects in detail, we should define them appropriately and leave no scope for confusion.The MC or IMC is always defined by looking at the behaviour of the condensates -its increasing value signifies catalysis and decreasing value implies to inverse catalysis.The trend of T CO plays no role there.It is a coincidence that for the physical parameters IMC effect is accompanied with a decreasing T CO .In fact, one can conjure up scenarios in which the condensate values enhance (thus MC) but the T CO still decreases [55].
Apart from the NJL model [32,33,35,46], there are other effective models which have been utilised to study QCD in presence of a magnetic field.example, linear sigma model (LSM) with quarks is used to explore the IMC effect in Ref. [34], whereas its Polyakov loop extended version (PLSM) and hadron resonance gas (HRG) models are utilised to explore thermodynamic quantities like magnetisation etc in Ref. [36].PLSM model has also been used further in the presence of magnetic field in Refs.[37][38][39].

Constant coupling
With the constant coupling, we begin the discussion by showing the plots of uniquely (10 m 2 π ) magnetised condensates for different values of the q parameter in Fig. 3.Because of their different coupling strength to the magnetic field the up and down quark condensates are now different.They are shown in the first and second panel.The overall effect of increasing q value remains the same as it is in zero eB case (Fig. 1), i.e., the condensate decreases with increasing q for a given value of T .In the third panel, we describe the condensate average with the same criteria.It is important as it will be used to define the chiral crossover temperature5 as well as the the fate of (I)MC will be decided by looking at its behaviour.On the other hand, to track the behaviour of the individual quark condensate is important, as they respond differently in the presence of a magnetic field.The values of the T CO for different values of q are given in Table II.As expected from the behaviour of the condensates, T CO decreases as a function of q.To understand the percentage decrement, we provide with a scaled T CO plot in Fig. 4. The reduction in the crossover temperature can be as high as 10% for the highest used value of q. q 1 1.05 1.10 1.15 T CO (MeV) 182 177 172 167 If the extent of the nonextensivity of the medium is fixed and the magnetic field is varied then the effects on the condensates are shown in Fig. 5.For the sake of comparison, the scenarios for the extensive case are also shown with the dashed lines with increasing strength of the magnetic field.The known MC effect is evident there.The MC effect also prevails in the nonextensive medium with two noticeable differences.First, at a given eB value, the strength of the condensate for q = 1.1 is smaller than the extensive case, particularly in the vicinity of the transition region.In fact, it is within this transition region that the differences between the extensive and nonextensive media are most conspicuous, especially when dealing with lower magnetic field values.Secondly, the fall of the curves are comparatively flattened in a nonextensive medium.These affect the fate of the chiral dynamics of the system which we discuss next.With MC prevailing for both extensive and nonextensive cases, the weakening of the chiral condensates and the corresponding transition temperatures are evident in Fig. 5.This is what reflects in the Table III, where we provide T CO 's for both q → 1 and q = 1.10 with different values of eB.Though the T CO for a given value of eB is smaller for a nonextensive medium as compared to the extensive one, the percentage change is always higher for the former for all values of eB considered here.It is depicted in Fig. 6, where the T CO 's are scaled with the respective values at zero eB.

Magnetic field dependent coupling constant
For the field dependent coupling constant, we first show the condensate plots for the light quarks along with the average for the known extensive scenario.It is important to note that we are not using the modified definition of the condensates utilised in Ref. [33] (Eq.14) which are taken from LQCD [5].Such definitions are particularly used in effective models to have results on the same footing and thus comparable with that in LQCD.
In the present study, we do not have such obligation as we are not comparing with any existing data.Rather, we are interested to investigate the fate of the (I)MC effects for a nonextensive medium in the model itself.Thus for us, knowing the condensates and their average will be adequate to decide on the fate of the (I)MC effects and the trend of T CO .To determine the nature of the effect (whether MC or IMC) the condensate averages will be sufficient as argued in the third paragraph of Sec.III B.
Before going to discuss the results a few comments are necessary.In Ref. [40], it was shown that the coupling constant decreases with the increasing magnetic field strength as α s (eB) = 1/ b ln(|eB|/Λ 2 QCD ) , where α s is the strong interaction coupling constant and b = (11N c − 2N f )/6π.If we consider G S to be analogous to α s , then it should as well decrease with eB.Thus, this form can be used as an ansatz to model the running coupling in NJL model.Since, we plan to make this study an extensive one, we explore this possibility as well and put it in the appendix IV.In Fig. 7, the condensates in an extensive medium with a magnetic field dependent coupling constant are plotted.It is evident from the plots that IMC effect, expectedly, prevails around the transition region.This is similar to what is found in the existing literatures [32][33][34][35][36][37][38][39]46].Though the status of IMC will be finally decided by the trend of the condensate average, it is also interesting to observe the individual flavour.What we notice is that the up and down quarks are not affected by the IMC effect with equal magnitude.This comes as no surprise considering the different coupling strengths of up and down quarks with the magnetic field.In fact, if one calculates the crossover temperature for individual flavour, they will be different.Though, as already mentioned, the crossover temperature that we quote will always be obtained using the condensate average.It will be interesting to study how these scenarios are affected in a nonextensive medium, which we discuss next.In Fig. 8, we show the same plots as those in Fig. 7 but for a nonextensive medium with q = 1.10.We observe some interesting changes in the behaviour of the individual as well as the average condensates.It is clear that the IMC effect for up quark is spoiled for all the strengths of the magnetic field considered here.For the down quark, the effect is present almost for all values of eB but with a reduced strength.Thus, two light quarks, expectedly, behave differently in presence of a magnetic field.As for the condensate average, which we use to decide on the types of the effect in the medium, the IMC effect is ruined by the "nonextensiveness" of the medium except for the highest strengths of eB considered here.This is elaborated with an inset plot zooming in the crossing regions.It is obvious that there is an interplay between the strengths of q and eB, which eventually determines the types of catalysing effects (MC or IMC).We elaborate it with the next two figures.In Fig. 9, we plot the condensate averages for different values of eB with two different values of q.These two figures help us understand the interplay between eB and q on determining the types of the catalysing effects in the medium.For q = 1.05, the IMC effect reappears starting from eB = 30 m 2 π as shown in the figure in the left panel.Thus, decreasing q from 1.10 to 1.05 makes the IMC effect favourable for lower values of eB.On the other hand, if we increase q to 1.15 the IMC effect completely disappears as evident from the figure in the right panel.For both the panels, the transition regions are zoomed in for the inset plots to make the crossings of the condensates conspicuous.Now, we make a comparison of the T CO 's for different values of q including the extensive case with zero magnetic field in Table .IV.For q → 1, the decreasing nature of the crossover temperature is captured.As we crank up the nonextensive nature of the medium by increasing q, T CO starts behaving differently.For q = 1.05,Initially it does not change and then it decreases with increasing eB.With further increase in q the trend becomes non-monotonic.For both q = 1.10 and 1.15, T CO first increases and then decreases.With higher q values the increase is greater.
All these points are well illustrated in the Fig. 10, where we can also read out the percentage change in the crossover temperatures and their trend following the liens which are solely intended to guide the eyes.We also observe that decreasing T CO is not always accompanied with IMC.For example, at q = 1.15 the IMC effect is absent but we still observe decreasing T CO for the higher values of eB.Thus, simultaneous occurrence of the two is a mere coincidence as has been discussed in detail in Ref. [55].There it is shown that with higher pion masses there is MC effect but with decreasing T CO .
It is also noteworthy that such non-monotonic behaviors arise in the context of an effective model due to the interplay between eB and q, with q ≥ 1.10.Such q-values render the medium non-extensive.On the other hand, all the available LQCD studies [6,56,57] demonstrate a monotonic decreasing behaviour of the crossover temperature as a function of the magnetic field, which are obviously based on extensive thermodynamics.Finally, in Fig. 11 we show the plots for the condensates for different values of q with a given eB.This is analogous to Fig. 3 but for a varying coupling constant.The condensates along with their averages still display the decreasing features below the T CO 's with increasing q values just like Fig. 3 with a crossing just above the transition temperatures.This crossing appears because of the running coupling constant.
In Table .V, we quote the T CO values for different q at eB = 20 m 2 π .It is not affected by the crossings of the condensates and decreases with increasing q as in the case with a constant coupling.In Fig. 12, the scaled transition temperature is shown as a function of q to elucidate the percentage change.

IV. CONCLUSION
We have presented an NJL model analysis of a magnetised QCD medium in a nonextensive ambience.There exist study discussing chiral symmetry restoration and other properties like meson masses in such scenarios.We add to that ongoing effort by including the magnetic field.
To achieve our goal, we utilised a 2 + 1 flavour NJL model.We start the analysis by revisiting the zero magnetic field case and found that the strength of the chiral condensate decreases with increasing nonextensive parameter q.With increasing q the chiral transition temperature decreases as well.These findings are in accord with what we already know and put our study into a perspective.
While coming to the focal point of the present work, we include the magnetic field and investigate the chiral transition for a q-value greater than unity.We also analysed and discussed the results for q → 1 wherever deemed necessary.The analysis is mainly divided into two parts -with i) a constant coupling (G 0 S ) and ii) a magnetic field dependent coupling constant (G S (eB)).
For a constant coupling at a specific magnetic field strength, the impact of increasing q resembles that of the zero magnetic field case, leading to a decreasing trend in both the condensates and the transition temperatures.The well-known magnetic catalysis effect persists with increasing eB for a q-value greater than unity.However, it results in reduced strength in both the condensates and the transition temperatures compared to the extensive scenario.We also discovered that the percentage increase in the transition temperatures is consistently higher in a nonextensive scenario.
For a magnetic field dependent coupling the observations are more interesting.Depending on the strength of q and eB the observed IMC effect in BG statistics can be destroyed.For q = 1.10, the IMC effect, which is nothing but the decreasing condensates with increasing strength of eB, is missing, except for the highest strength of magnetic field (40m 2 π ) that we examined in this study.There is a competition between the two parameters q and eB which becomes evident when we analysed for other two values of q = 1.05 and 1.15.For q = 1.05, the IMC effect can be observed for eB = 30m 2 π , a field value lower than that required for q = 1.10.On the other hand, at q = 1.15, the IMC effect is totally eliminated for all the strengths of eB explored here.
With eB-dependent G S the behaviour of the transition temperature, argued to be independent of the IMC effect, turns out to be little subtle.It also depends on both q and eB.For q = 1.05, it remains constant initially and then starts decreasing with increasing magnetic field.At even higher values of q, it exhibits a non-monotonic trend as a function of eB -initially increasing before subsequently starts falling.The initial increment is higher for higher values of q.
We also explored the condensates using G S (eB) for different values of q with a given eB.There is a crossing of the condensates at a temperature slightly higher than the transition temperature as compared to no-crossing with constant G S .Also, in this case, the percentage decrease of the transition temperature is slightly lower as compared to the constant G S .
However, we must keep in mind that the temperature in Tsallis statistics can have a more generalised meaning due to its non-extensive nature.In this regard, using a more generalised statistics [58,59] could be useful.Such a statistics is characterised by a unique pair of scaling exponents and the Tsallis statistics along with the Boltzmann-Gibbs statistics become special cases.It has been recently utilised in Ref. [60] to successfully reproduce the particle ratios at two different energies in HICs.There exist other works [61,62] along with a review on the topic [63].Thus, it will be interesting to check how much of the present conclusion gets affected in such a generalised scenario and demands a separate investigation.
Having discussed these interesting findings within a mean-field approach, it becomes intriguing to explore the effects of fluctuations by transitioning to beyond mean-field approximations.In the pursuit of robust predictions, it will be interesting to test some of the present observations using beyond mean field approaches.For example, the weakening of the chiral condensates with increasing q parameter, the non-monotonic behaviour of the T CO as a function eB for q-values greater than unity etc.Such exercises are certainly beyond the scope of the present study and will be reported elsewhere.
As an immediate extension to the present analysis of the chiral symmetry for a nonextensive magnetised medium it would be really interesting to include the background gauge field to bring out the effect of statistical confinement.This will be particularly useful in the investigation of transport properties [64,65] with effective degrees of freedom following the technique outlined in [66].Such investigations could bear implications for HIC related observations.It also remains to be seen whether, by employing the techniques outlined in Ref. [67], the findings of the present effective model can be linked to certain phenomenological consequences.We plan to report on these in the near future.field case.These are in conformity with the observation in Ref [33].Thus, though this simple ansatz can provide a hint of IMC effect and decreasing T CO it fails to be in accord with the established evidence [5,6].It is useful as a hint and needs to be tweaked to be successfully used in an effective model [32,33,35].
As we extend this analysis for a nonextensive medium (q = 1.10), no important change is observed and the discussion in the above paragraph remains more or less valid.The only change that takes place is the reduction in the gap of T CO 's between the zero eB and highest value of eB considered.Thus we can safely conclude that the magnetic field dependent coupling constant with a simple ansatz does not behave much differently in a nonextensive medium.

35 T 15 FIG. 1
FIG. 1: (B = 0): Plot of the condensates (upper panel) and effective masses (lower panel) for different values of the q parameter.Left panel: light quarks.Right panel: strange quark.

15 FIG. 3 :
FIG.3:(B ̸ = 0, G 0 S ): Plot of the condensates for different values of q at eB = 10 m 2 π .First two panels show up and down quark condensates, respectively.The last panel is for condensate average.

m π 2 FIG. 8 :
FIG. 8: (B ̸ = 0, GS(eB)): Plot of the condensates and their average for different values of eB in a nonextensive medium.Panels are designated in the same way as in Fig. 7.

15 FIG. 11 :
FIG.11:(B ̸ = 0, GS(eB)): Plot of the condensates and their averages for different values of the q parameter.Left to right panels represent up quark, down quark and average condensates, respectively.

m π 2 FIG. 13 :
FIG. 13: (G 0 S /GS(eB)): Plot of the condensate averages for different values of eB both in the extensive (left panel) and nonextensive (right panel) scenarios.

TABLE II
S ): Scaled TCO as a function of q for eB = 10 m 2 π .

TABLE III :
(G 0 S ): Chiral transition temperature for different values of eB and q.