Quantum version of transport coefficients in Nambu–Jona-Lasinio model at finite temperature and strong magnetic field

We have estimated parallel and perpendicular components of electrical conductivity and shear viscosity of quark matter at finite magnetic field and temperature by using their one-loop Kubo expressions in the framework of Nambu–Jona-Lasinio (NJL) model. At finite magnetic field, a non-trivial medium dependence of those quantities can be found. Previously these NJL-profiles have been addressed in relaxation time approximation, where cyclotron motion of quarks with medium dependent mass plays the key role. With respect to the earlier estimations, the present work provides further enriched profiles via Kubo framework, where field theoretical descriptions of quark transport with medium dependent mass and (Landau) quantized energy have been identified as the key ingredients. Hence the present study can be considered as the complete quantum field theoretical description of the transport coefficients in the framework of NJL model at finite temperature and magnetic field.


I. INTRODUCTION
Production of strong magnetic fields in the early stages of relativistic heavy-ion collisions (HIC) is a longstanding topic, that is being extensively studied [1,2].The strengths of these produced fields have been estimated to be even larger than the strong-interaction scale Λ 2 QCD ≃ 0.06 GeV 2 (e.g.Pb-Pb collisions at the Large Hadron Collider estimates  ∼ 15 2  ≫ Λ 2 QCD ;   is the pion mass ∼ 0.135 GeV) [3] which subsequently indicates that these fields heavily influence various observables in the hot and dense quark matter such as the quark condensates.Many such modifications have already been studied but the interpretations of those modifications in the system are still ambiguous [4][5][6][7].One of the biggest challenges to understanding these magnetic field-induced modifications is to grasp the time dependence of the produced magnetic field in the early stages of HIC.There are several schools of thought on this topic.Some of the preliminary studies indicated that the fields weaken fast as the system expands [3,8].Then there are also studies that suggested that the induced electric currents in the expanding matter due to the produced magnetic fields can in turn produce magnetic fields again, overall changing the longevity of the early-produced fields [9][10][11][12].Very recent studies have again given emphasise on the short lifetime of the induced magnetic field [13] and further suggested that the search for any magnetic effects in the HIC would be highly challenging.The absence of CME signals in the isobar experiments [14] also concurs with this inference.In Ref. [15], the validity of Ohm's law has been argued in view of the rapidly evolving quark matter produced in HIC and the behavior of the time dependent conductivity has been discussed.All these studies give strong indications that electrical conductivity can be a really important quantity in this scenario.Hence, understanding the microscopic calculation of transport coefficients like electrical conductivity at finite magnetic fields might be considered an important topic to study.Electrical conductivity and other transport properties such as shear and bulk viscosity could have been computed unambiguously using Lattice QCD, a nonperturbative first-principles numerical method formulated in Euclidean space, if not for the crude inversion techniques to reconstruct the Minkowski spectral functions from Euclidean correlation functions.However, there are existing lattice QCD results for transport coefficients with SU(2) quenched simulations [16] and with full QCD SU(2+1) simulations [17].There are also several other recent analytical nonperturbative studies of electrical conductivity at finite magnetic fields [18][19][20][21][22][23][24][25][26][27][28][29][30].On the other hand, simulations of the field dynamics invariably involve solving relativistic and dissipative magnetohydrodynamics equations, which require other transport coefficients like shear and bulk viscosity.The magnetic field dependence of the shear viscosity was computed recently in Refs.[31][32][33][34][35][36][37][38][39][40] and of the bulk viscosity in Refs.[41][42][43][44][45].
Former formalism connect one-body kinematics with equilibrium distribution function and its deviation to many-body mechanics, appeared as thermodynamics and transport coefficients of the system.Next, the deviation is obtained in terms of relaxation time   by using the RTA based Boltzmann transport equation.On the other hand, Kubo relation defines transport coefficients as transportation probability field operators like shear stress and electro-magnetic current between two points.So, the diagrammatic quantum field theory calculation becomes the scope of this framework.It is quite fascinating that these two completely different methodologies converge into exactly same expressions of transport coefficients in absence of magnetic field, which is well established by earlier Refs.[47][48][49][50].On the other hand it develops curiosity again, when finite magnetic field extension of RTA expressions [36][37][38][39][40] and Kubo expressions [29,30,46] don't coincide.It is because of the quantum field theoretical aspect of the magnetic field, which was missing in the RTA picture, and properly incorporated in the Kubo expressions.Refs [29,46] have gone through leading order estimations of conductivity and viscosity components, where propagators at finite temperature () and magnetic field () have been considered.A rich quantum field theoretical (QFT) structure is noticed in the parallel and perpendicular components of those transport coefficients.In the present work we intend to explore the effects of that field theoretical structure on the transport coefficients within a particular system described by an effective QCD model, carrying a nontrivial  and  dependence.We have highlighted on the QFT modification of transport coefficients by comparing with their corresponding RTA expressions.These comparisons were already done in Refs.[29,30,46] for massless fermionic or bosonic system in general but the present work has extended it into a more specific and realistic system -quark matter in the light of effective QCD model framework.
Among the several existing effective QCD models, the Nambu-Jona-Lasinio (NJL) model [51,52] is well adopted for QCD phenomenology at finite temperature [53][54][55][56].At finite  and , quark condensate and constituent quark mass become ,  dependent functions via magneto-thermodynamical phase-space, built from the Landau quantization technique of thermal field theory (TFT).In the present work, we have used the NJL model of Refs.[57,58], a model that compliments the novel lattice QCD results [59][60][61][62], which first showed that strong magnetic fields have dramatic effects on the QCD phase diagram.These lattice simulations have found that though the magnitude of the light quark condensates increase with the magnetic field for low temperatures, they start to decrease for temperatures close to  pc ≃ 0.16 GeV, the region associated with the chiral symmetry restoration.Former case is connected with magnetic catalysis (MC) and latter case with inverse magnetic catalysis (IMC).These phenomena are well reviewed in Refs.[63][64][65][66][67][68].
As first attempts, Refs [28,36] have provided estimations of shear viscosity [37,38] and electrical conductivity [21,38] in presence of magnetic field by using RTA expressions within the NJL model.But the field theoretical structure rich Kubo expressions [29,46] of the transport coefficients have never been explored within the NJL model incorporating IMC.In the present work we have specifically explored this.
The paper is organized as follows.Sec.II has gone through the main results of NJL model calculation at finite magnetic field.Next in Secs.III, IV, RTA and Kubo expressions of transport coefficients at finite magnetic field are respectively addressed.Then, our Kubo or field theory based NJL estimation of transport coefficients are plotted and their additional ,  profiles with respect to their corresponding RTA or classical based estimations are discussed in Sec.V.At the end, Sec.VI has summarized the study with the new findings.

II. NJL MODEL IN PRESENCE OF MAGNETIC FIELD
For our purpose we have chosen here the isospin-symmetric two-flavor NJL model, whose Lagrangian density in presence of an electromagnetic (EM) field (  ) is given by where  depicts  and  quark fields iso-doublet, each being an   (= 3)−plet,   being the number of colors.m = diag(, ) is the quark-mass matrix within exact isospin symmetry and hence in the rest of the paper we will work with the notation ) is the charge matrix in the flavor space and  is the coupling constant of the NJL model.Solution of this model in the mean-field approximation corresponds to the leading-order approximation in the 1/  expansion.Since the NJL model is unrenormalizable, in presence of an external magnetic field isolating the divergences from the vacuum structures comprised of various Landau levels requires extra care.Based on well known results of [69][70][71], and recently explored in Ref. [72], the magnetic field independent regularization (MFIR) present itself to be a satisfactory method to study NJL model in a constant and external magnetic field.The advantage of this regularization scheme is the complete separation of the magnetic field contributions from the vacuum or thermal contributions and subsequent regularization of the vacuum term using standard procedures.On the other hand, most of the non-MFIR based regularization schemes try to remove the divergences at each individual Landau levels by using sharp cutoff regulator functions depending on that particular Landau level.Sharply cutting off the divergences at each Landau levels generate abrupt transitions between Landau levels which can lead to unphysical results [73,74], as e.g.oscillations in the chiral quark condensate.Although use of smoother regulator functions improve the situation and also help us identify the unphysical oscillations from the possible physical ones [72].
In the mean field approximation, the gap equation for the constituent quark mass  at finite temperature  and in the presence of a magnetic field  is given by where ⟨ ψ    ⟩ is the quark condensate of flavor  .In presence of an external magnetic field , ⟨ ψ    ⟩ can be written as a sum of three contributions [57,70,71]: with where Γ(  ) is the Euler gamma function,   =  2 /(2|  |).Within MFIR regularization scheme, the vacuum term ⟨ ψ    ⟩ vac given in Eq.( 4) is ultraviolet divergent and we use a three-dimensional cutoff Λ.The parametrization of the model is determined treating the coupling , the current quark mass  and the cutoff Λ as free parameters, which are fixed by fitting the vacuum values of the pion mass   , pion decay constant   and quark condensate ⟨ ψ    ⟩.In addition, in Eq. ( 6)  represents Landau levels, with   = 2 −  ,0 being the spin degeneracy factor and  0 (  , ) is the Fermi-Dirac distribution function: where As mentioned earlier, in the present work, our aim is to include the effects of the inverse magnetic catalysis (IMC) phenomenon on the quasi particle effective/constituent quark mass .Usual NJL model with a fixed coupling constant , in the quasi-particle approximation is unable to describe IMC [59,60].One effective way to resolve this issue is to impose that the coupling constant  of the model is − and −dependent [57,58,75].Using  (, ) a precise description of the lattice results for the  and  quark condensates has been obtained within the NJL model with the parametrization [58]: where (), (),   () and () depend only on the magnitude of  and their values for selected values of  are given in Table 11 The expression given in Eq. ( 9) has been adopted in Ref. [58] for mere convenience since it is well adapted to fit LQCD results.All following numerical results refer to the parametrization adopted in Ref. [58].The described model with  (, ) was used to study the effects of a magnetic field on neutral pion mass in Refs.[76,77], which subsequently produced results that agree with corresponding lattice QCD results [78].For further evidences and inputs about the IMC phenomenon readers can look into refs.[61, 66-68, 79, 80].

III. RELAXATION TIME APPROXIMATED (RTA) EXPRESSIONS
Our main aim is to apply the NJL model for calculating transport coefficients, based on quantum field theoretical calculation at finite temperature and magnetic field.To realize QFT contribution in transport coefficients for quark matter within NJL model, we will first revisit their expressions based on kinetic theory, as we intend to compare between them.First we will revisit the expressions of transport coefficients without Landau quantization, which we have broadly mentioned as RTA expressions throughout this paper.Refs.[28,36] have already gone through the NJL model estimations of the RTA expressions of transport coefficients, whose framework will be quickly and briefly revisited in the next subsections for the sake of completeness.We want to emphasize here again that the present article is aimed to zoom into the quantum aspects of transportation and the transition from classical to quantum estimations, within the framework of NJL model.

A. RTA expressions of Electrical Conductivity
The conductivity tensor    can be realized as a proportional connector between current density   and electric field   via macroscopic Ohm's law   =      .In the microscopic picture of dissipation, we can assume that the equilibrium distribution function of quark at zero quark chemical potential,  0 = 1    +1 undergoes a small deviation   .Therefore, one can express (dissipative) current density as [21,38,39,81] where 4  = 2 × 2 ×   are respectively spin, particle-anti-particle, and color degeneracy factors of medium constituent quark with electric charge   and energy  = { ì  2 +  2 } 1/2 .Using relaxation time approximation (RTA) of the Boltzmann transport equation (BTE) one can know the form of   [21,38,39,81].
To do the calculations, electromagnetic field geometry should be fixed first.Here, we fix our external magnetic field direction along the z-axis and electric field will be considered first along the x or y direction for calculating the perpendicular conductivity    or    and then subsequently along the z-axis for calculating the parallel conductivity   .So for the first case,  is along z direction and  is along x direction.Using this condition of electromagnetic field geometry in RTA based BTE, we will get [38,39] where   is relaxation time and    =     is inverse of cyclotron frequency.The second term in the RHS of Eq. ( 12) is related to the Hall conductivity and since the sign of    for particles and anti-particles will be opposite, so the net value of the Hall component will be zero at zero quark chemical potential.At non-zero quark chemical potential, due to an imbalance of particle and anti-particle densities, one can get non-zero net Hall conductivity [40].But that is not our matter of interest here as we restrict ourselves to the zero (net) quark density zone.Hence, from the first term of Eq. ( 12), we will get the perpendicular (   ) component of conductivity tensor [21,38,39,81] Similar to    RTA , one can get

𝑦 𝑦
RTA by considering electric field along y direction instead of x direction and it can be easily found that both have same expressions.Both components can be called perpendicular components in general.
Next for the parallel component of conductivity tensor (  ), electric and magnetic field both are applied along the z-axis, so the Lorentz force will not work along the z direction.Hence   comes out to be This parallel component remains the same as we get isotropic expressions in the absence of the magnetic field.We will see later that this magnetic field independent behavior of parallel conductivity in classical case will be modified in the quantum picture.

B. Introducing Landau quantization in electrical conductivity: QM expressions
After getting the RTA expressions of conductivity, here, we will add their quantum aspects.We know that Schrodinger's equation in presence of a magnetic field converts into a quantum harmonic oscillator-type problem with quantized energy, which is known as Landau quantization.For relativistic spin 0 and spin 1/2 particles, one has to respectively solve Klein-Gordon and Dirac's equations in a magnetic field to find out the quantized energy with Landau levels .Here, we will straightaway impose the Landau quantization directly.For the rest of this manuscript, we have mentioned the Landau quantized RTA expressions as the QM (quantum mechanical) expressions.The main changes will be within the phase space integration 2 ∫  3  (2  ) 3 , which will be changed to =   .Based on these modifications, Eqs. ( 13) and ( 14) will now be expressed as :

C. RTA expressions of Shear viscosity
After electrical conductivity, let us also quickly revisit the semi-classical expression of shear viscosity, whose detailed calculations can be found in Refs.[35][36][37][38][39]81].The main motivation of this approach is to build a connection between macroscopic fluid description and microscopic kinetic theory description.Similar to electrical conductivity, which is realized as a proportionality constant between current density and electric field, shear viscosity also becomes the proportionality constant between viscous stress tensor and velocity gradient.In absence of magnetic field, the general form of Newton-Stoke relation is given by    = U   between viscous pressure    and velocity gradient Now in presence of a magnetic field, one can construct five independent velocity gradient tensors, and hence, we get five shear viscosity components - 0,1,2,3,4 .So, the magnetized medium modified Newton-Stoke relation can be written as where velocity gradient U    will be constructed by unit vector of magnetic field along with the fluid velocity.There are two different possible sets of five independent trace less tensors, prescribed in Ref. [81] and Refs.[42,43].Readers can find the detailed structure of U    in those Refs.[42,43,81].In the microscopic picture of viscous dissipation, we can assume that due to velocity gradients U    , the equilibrium distribution function  0 undergoes a small deviation    along the directions of different velocity gradients U    .Therefore, one can express the (dissipative) viscous pressure as [21,38,39,81] Using RTA based BTE, one can know the form of    [35-39, 82, 83] and get the final microscopic expressions of   's.If we analyze the earlier existing references, we can find that Refs.[35][36][37][38][39] have adopted former set of tensors and Ref. [38,82,83] have adopted latter set of tensors.Now five shear viscosity components, obtained from two different tensors are inter-connected [38,42,43], so one can proceed to use any one of the sets tensors.Here, we will use the tensors, prescribed by Refs.[42,43] and RTA expressions of  0,1,2,3,4 , whose detailed derivation can be found in Ref. [38].Among them,  3,4 is the Hall viscosity, which will be zero for vanishing quark chemical potential, just like the Hall conductivity.In the present work, we will focus on parallel ( ∥ ) and perpendicular ( ⊥ ) components of shear viscosity (considering the magnetic field is along the z-axis), which are related with  0,2 as [38,84,85] where RTA expression of  2 is [35][36][37][38][39]81]:

IV. KUBO EXPRESSIONS
In Sec.III, we have discussed about RTA based expressions of transport coefficients like shear viscosity and electrical conductivity.In this section we will discuss an alternative methodology, where one can find a quantum field theoretical structure of these transport coefficients.Owing to the Kubo relation, one can represent different transport coefficients as zero momentum limit of thermal correlator for corresponding field operator.These correlators can be expressed as one-loop self energy diagrams by considering the kinetic/free part of Lagrangian density.With further inclusion of thermal width Γ  (which can also be considered as the inverse of relaxation time   in RTA approach) in propagator, one can get a non-divergent expression of transport coefficients.In absence of magnetic field, this procedure yields one loop Kubo expressions of different transport coefficients, which are exactly same as their RTA expressions.But in presence of magnetic field the RTA and Kubo expressions differ from each other [29,30,46].In the present article, we have estimated those Kubo expressions within the NJL model.In the present section we discuss those expressions for electrical conductivity and shear viscosity.For more details on the following formalism one can look into Refs.[29,30,46].

A. Kubo expressions of Electrical conductivity
In the Kubo formalism, the electrical conductivity tensor   is calculated from the long-wavelength limit of the in-medium (EM) spectral function   as Here, the spectral function   is related to the the Fourier transform of the vector current-current correlator given by in which ⟨...⟩ denotes the ensemble average, T  is the time ordering with respect to symmetric Schwinger-Keyldish contour  in the complex time plane as used in the real time formalism (RTF) of finite temperature field theory, and, 11 refers to the fact that the two points (in time ordering) are on the real horizontal segment of the contour .Now, using vector current   = −   of 2-flavor NJL quark matter for charge  and Wick's theorem in Eq. ( 23), we get where, Φ  (, );  ∈ {, } is the gauge dependent phase factor and   11 ( ) is the 11-component of the momentum space thermo-magnetic quark propagator in RTF of flavor  , explicitly given by [86,87] in which  () = Θ()  0 () + Θ(−)  0 (−),  0 () denotes the Fermi-Dirac thermal distribution function already defined in Eq. ( 7),  is the Landau level index,    = √︃  2 + 2   and D  ( ) is with the convention  −1 () =  1 −1 () = 0 for the Laguerre polynomials   .While writing the propagator in presence of magnetic field in ẑ direction, we took 27) into Eq.( 26) followed by substituting in Eq. ( 23), we get [29,30,46] where, Now substituting Eq. ( 29) into (22), we obtain the conductivity tensor as Having obtained the conductivity tensor, it is easy to extract the parallel and perpendicular components of conductivity using where, the projectors P ∥ ,⊥  are given by in which,   = 1 2         ,   is the fluid four-velocity,   =    ,ext −    ,ext is the field strength tensor, antisymmetric in Lorentz indices,   =        and Ξ  = Δ  +     .In the local rest frame (LRF) of the fluid,   LRF = (1, 0, 0, 0) and   LRF = (0, 0, 0, 1) points along the direction of external magnetic field.Substituting Eq. ( 31) into Eq.( 32) and performing the  2  ⊥ integral analytically we finally arrive at where, It is to be noted that  −1 −1 = 0 which follows from the convention  −1 () =  1 −1 () = 0 used in Eq. ( 28).Few comments on using the thermo-magnetic quark propagator of Eq. ( 26) and ( 27) while calculating the in-medium spectral function   in Eq. ( 23) are in order here.We note that, the local current   () = −()  () appearing in Eq. ( 23) consists of Heisenberg quark field () corresponding to the full interacting Hamiltonian of the model that includes the effects of effective 'strong' interaction as well as of static (due to external magnetic field) and dynamical EM interaction at finite temperature and/or density.This in turn implies that, the quark propagator  11 (, ) = T  ()() 11 appearing in the expression of the spectral function   in Eq. ( 25) should be the complete 'dressed' propagator of the quarks carrying the effects of all the interactions.Now, the quark propagator in Eq. ( 27) which is used to calculate the spectral function, does contain both the effects of effective 'strong' as well as static EM interaction due to external magnetic field.In particular, the temperature and magnetic field dependent constituent quark mass  =  (, ) entering in the expression of the propagator, captures the effect of 'strong' interaction; whereas, the use of Schwinger proper-time formalism along with RTF takes into account the effects of static EM interaction due to external magnetic field and finite temperature to all orders.However, in this work, we have neglected the effect of dynamical EM interaction at finite temperature in the quark propagator which is of the order  em  2 [88] where  em = 1/137 is the fine structure constant.Clearly, the effect of dynamical EM interaction is sub-leading as compared to the effect of 'strong' interaction as well as to the external magnetic field scale  ≫  em  2 for the temperature and magnetic field ranges considered in this work.
From Eq. ( 35) one can clearly see the difference between Kubo and RTA expressions.In Sec.III B, we have explicitly imposed Landau quantization in RTA expressions (Eqs.( 15), ( 16)) but in the Kubo expression, Landau quantization is implicitly accounted for, originating from the thermo-magnetic quark propagator.

B. Kubo expressions for Shear viscosity
Similar to the previous subsection, the viscosity components can analogously calculated from the Kubo formalism.The fourth rank viscosity tensor V   is calculated from the long-wavelength limit of the in-medium spectral function    of the energy-momentum tensor (EMT) as Here, the spectral function    is related to the the Fourier transform of the two-point EMT-EMT correlator given by Now, for 2-flavor NJL quark matter in MFA, the EMT is given by, Substituting Eq. ( 40) into Eq.( 39), we get [29,30,46] where, in which, Finally substituting Eq. ( 41) into (38), we obtain the viscosity tensor as Having obtained the viscosity tensor, it is easy to extract the parallel and perpendicular components of shear viscosity using where, the projectors P ∥ ,⊥   are given by Finally substituting Eq. ( 46) into Eq.( 47) and performing the  2  ⊥ integral analytically we finally arrive at where, in which, At the end, let us put our working formulae in Table (II).
Table II.Equation numbers of  ⊥, ∥ and  ⊥, ∥ , which will be the working formulae in the result section.

V. RESULTS
In this section we will discuss the numerical results of Kubo expressions within NJL model and compare them with previously explored RTA and QM expressions to reveal their additional contributions.Since NJL model at finite temperature and magnetic field provides us magneto-thermodynamic phase-space information of QCD from non-perturbative to perturbative domain, our plan is to identify them.For this purpose we have to take the massless quark results, which are close to pQCD results, as a reference to identify the non-perturbative effects by looking at the differences between massless and NJL model results.As a result, we will be careful in our graphical representations to reveal the different changes like classical to quantum, or pQCD to non-pQCD transformations.
Let us start with electrical conductivity curves, shown in Fig. 1 where RTA, QM, and Kubo results are denoted by dotted, solid, and dashed curves.RTA curves are built from Eqs. ( 13) and ( 14) by using NJL based quark mass  (, ), obtained from Eq. ( 2).If we take the massless limit, parallel conductivity becomes proportional to  2 .So by choosing appropriate normalized quantities, we can get  and  independent values [38,39], e.g.: Horizontal dotted line in the left-upper panel of Fig. 1 indicate this  and  independent nature of mass less RTA curves.Similar to Stephan-Boltzmann (SB) lines in the plots of lattice quantum chromodynamics (LQCD) based thermodynamics [59,60], we may treat these horizontal curves as a reference line of transport coefficients.One can notice that when we consider  (, ), we get suppressed values of  ∥    2 concerning its massless values.The same trend is noticed in LQCD thermodynamics [59,60].This suppression can be realized as a non-pQCD effect, observed in both thermodynamics and transport coefficients of QCD medium.Usually, it is the LQCD calculations [59,60], effective QCD models like NJL model [57,58] or LQCD mapping quasi-particle model [39], which attempt to map non-pQCD contributions at finite temperature and magnetic field.
Next, we go for the QM and Kubo curves of parallel components by using Eq. ( 16) and Eq.(35), where Landau quantization is incorporated.We draw QM and Kubo curves separately to establish that Kubo curves cover all the field theoretical quantum effects, whereas QM curves cover a partial aspect of it.So introducing Landau quantization by hand in RTA expressions, through which we get the QM curves, might not be recommended for getting the full quantum effect in transport coefficients.Here, our final focal interest will be the Kubo curves.We present the QM curves as a bridge between the RTA and Kubo curves.At low  in the left-upper panel of Fig. 1, we notice that massless QM and Kubo curves deviate from their horizontal RTA curve.In Ref. [59], a similar pattern has observed for massless RTA and QM thermodynamical quantities like entropy density.Landau quantization plays an important role for this deviation from the horizontal line in both cases -the entropy density graph in Ref. [59] and the plot of  ∥    2 in the left-upper panel of Fig. 1.When we used  (, ) in the QM and Kubo expressions of  ∥ , we again get suppressed results in low  concerning their massless results.Right-upper panel of Fig. 1 shows that the enhancement of quantum estimations of  ∥ compared to its classical values will increase with the magnetic field.Classical and quantum curves will tend to merge at high  and low  domains, which one can understand as the classical domain.On the other hand, quantum estimations of  ∥ remain larger than classical estimations in the low  and high  domain, which can be understood as the quantum domain.This phenomena is quite interesting and could possibly connect to the magneto-resistance in the domain of condensed matter physics [89,90], but a systematic and comparative study in future would be required to better comment on this topic.
Next, we explore the perpendicular component of the electrical conductivity.From Eq. ( 13), one can identify the effective relaxation time [38,39] with an approximate value inverse cyclotron frequency of massless quark matter [38,39] :  Due to this term,  ⊥ /(   2 ) curve of massless matter does not remain horizontal along  and  axes like  ∥ /(   2 ).Replacing the zero mass by NJL based  (, ), we will get further suppression due to non-pQCD effect.Red and black dotted lines for massless and  (, ) in left-lower panel of Fig. 1 display this fact.Next, QM (solid) and Kubo (dashed) curves of  ⊥ /(   2 ) are generated by using Eqs.( 15) and ( 35) respectively.They are different in massless case [46] as well as for NJL based  (, ).The main reason is the transformation to an effective relaxation time expression: In both the cases of RTA and QM, we consider classical concept of cyclotron motion using   , while in Kubo case, a transition between two energy levels, separated by unit Landau level difference [46] has come into the picture to describe the perpendicular conductivity components.Similar to parallel conductivity, perpendicular conductivity increases with respect to its classical values as we increase the magnetic field.One can observe this in the right-lower panel of Fig. 1.Reader may notice that the right panel graphs, plotted against -axis, carry few data points because of unavailability of LQCD data of .However, it does not hinder guessing the approximated trend of the curves.We would also like to mention at this point that generating the QM or Kubo results at  → 0 is an impossible task (as we have to consider summing over an infinite number of Landau levels).But the coincidence of all three curves at  = 0 is expected and alternatively checked.Next, let us come to the other transport coefficient considered in our present study -shear viscosity, whose parallel and perpendicular components are our matter of interest.In absence of a magnetic field, both components are the same and their massless limit follows  4 dependence as conductivity follows  2 dependence.So we choose the dimensionless ratio /(   4 ), which will remain horizontal along  axis for  = 0 case.At finite , RTA expressions of parallel and perpendicular components of the shear viscosity are given in Eqs (20), whose massless expressions can be simplified to [38]   where and same   's as given in Eq. ( 60).According to above simplified expressions, normalized values of  ∥ /(   4 ),  ⊥ /(   4 ) of massless quark matter will increase with  and decrease with .When we use constituent quark mass  (, ) from the NJL model, we will get suppressed values of RTA curves (solid line) concerning their massless curves (red solid line).The suppression represents the non-pQCD effect in the parallel and perpendicular components of shear viscosity.When we go for Kubo expressions, given in Eqs. ( 50)-( 52), we will get enhanced values of shear viscosity compared to their RTA values as we have noticed for electrical conductivity case.However, these enhanced Kubo values of transport coefficients are expected in the quantum domain only i.e. at high  and low  domain.
At the end of this result section, we want to highlight again that the new ingredient of the present work is the Kubo estimation of NJL matter, whose RTA and QM estimations are already addressed in earlier Ref. [28] for electrical conductivity and Ref. [36] for shear viscosity.For zooming in on the Kubo contribution, the RTA and QM curves are presented here for comparison.

VI. SUMMARY
Present work is aimed to highlight quantum field theoretical contribution at finite temperature and magnetic field in the transport coefficients of quark matter within the framework of NJL model.Earlier Refs.[28,36] have calculated the transport coefficients like shear viscosity and electrical conductivity of quark matter within the NJL model in the framework of relaxation time approximation (RTA).With temperature and magnetic field dependent constituent mass those estimations might be considered classical or semi-classical approaches.Present work provides their quantum field theoretical version using their Kubo expressions, obtained in Refs.[29,46].According to the Kubo relations, one can realize the transportation from one point to another as the propagation probability at finite temperature and magnetic field of relevant field operators like electrical current   and viscous stress tensor   .On the other hand, RTA framework describe a classical picture of transport phenomena in terms of cyclotron motion of quarks due to Lorentz force.Hence in the present work we have explored this transition from classical or RTA to quantum or Kubo picture of transport coefficients in view of the NJL model.Additional temperature and magnetic field profile in Kubo estimation is the primary content of the present work, which we have subsequently compared with corresponding RTA estimations, addressed earlier in Refs.[28,36].We also present another transition here, from massless quarks to constituent quark masses within NJL model.In low temperature and high magnetic field domain, we notice an enhancement of transport coefficients due to the transition from classical to quantum picture.On the other hand we also observe that their values got reduced during the transition from the massless case to NJL matter.Cumulative effects from both the transitions finally give us the complete field theoretical non-pQCD estimations, which remain little higher than their semi-classical values, addressed in Refs.[28,36].These differences vanish in the classical domain, i.e. high temperature and low magnetic field domain.At this point we want to point out that recently some studies have been directed towards exploring the extended version of RTA [91,92] which incorporates quantum corrections to the RTA dissipative current considered in the present study.This puts forward a good future avenue which can possibly bridge the gap between the semi-classical and the field theoretical pictures observed in the present study.
Finally we emphasize again that the present microscopic calculation is a complete quantum field theoretical estimation of transport coefficients within the framework of NJL model at finite temperature and magnetic field.The enhanced values of the conductivity at high magnetic field domain incorporating the full quantum effect may be a good signal to show that it helps slowing down the rapid decay of the magnetic field [4] produced in heavy ion collision experiments.But as discussed in the introduction, further studies are required to firmly conclude that.

Table I .
(9)ues of the fitting parameters in Eq.(9).Units are in appropriate powers of GeV.