Spectral function and dilepton rate from a strongly magnetised hot and dense medium in light of mean field model

We have calculated the electromagnetic spectral function (SF) vis-\`{a}-vis the dilepton rate (DR) by evaluating the one-loop photon polarisation tensor for a strongly magnetised hot and dense medium. The calculation is performed in the ambit of mean field models namely Nambu\textendash Jona-Lasinio (NJL) and its Polyakov loop extended version (PNJL) in the lowest Landau level approximation. These models allow for a dynamical generation of quark mass which, evidently, gets affected in the presence of a magnetised medium. We have shown that the strength of the SF gets boosted because of the presence of dynamical quarks in lieu of the current quarks. It increases as we increase the magnetic field for a given value of temperature in both NJL and PNJL models. This increment is further reflected in the DR which is enhanced as compared to the Born rate for certain values of invariant masses. We also discuss the impact of chemical potential on DR within the present scenario.


I. INTRODUCTION
Understanding Quantum Chromodynamics (QCD) in presence of magnetic background has gained lots of contemporary research interests [1]. It is important to study QCD in presence of external magnetic field not only for its relevance with the astrophysical phenomena [2][3][4][5][6][7][8] but also due to the possibility of strong magnetic field generation in non-central heavy-ion collision (HIC) [9] which sets the stage for investigation of magnetic modifications. Although the background fields produced in Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) are much smaller in comparison with the field strengths that prevailed during the cosmological electro-weak phase transition which may reach up to eB ≈ 200m 2 π [10], they are strong enough (eB ≈ m 2 π at RHIC and eB ≈ 10m 2 π at LHC) to cast significant influence on the deconfined medium of strongly interacting quarks and gluons known as quark gluon plasma (QGP) which is created in such HICs [11,12]. So far a large number of efforts has been put to understand the effects of this background magnetic field on the strongly interacting QGP. This results in a plethora of novel and interesting phenomena -finite temperature magnetic catalysis [13][14][15] and inverse magnetic catalysis [16][17][18][19][20][21][22]; chiral magnetic effect [23][24][25], chiral-and color-symmetry broken/restoration phase [26,27]; modification of dispersion properties in a magnetised hot QED medium [28][29][30], thermodynamic properties [31][32][33]; dilepton production from a hot magnetized QCD plasma [34][35][36][37][38][39][40] to name a few.
We know that dynamical properties of many particle system can be investigated through the correlation function (CF) [41][42][43]. Thus, by calculating the correlation function and its spectral representation many of the hadronic properties can be studied. These kind of properties for QCD in vacuum have been well studied [44]. Now in presence of hot and dense medium such vacuum properties are modified [45,46], because the dispersion relation of the particles moving in the medium gets modified. In presence of background magnetic field there will be further modifications which necessitates the modification of the available theoretical tools. From the temporal and spatial parts of the CFs, separately, we can extract important informations regarding the medium. As for example, the temporal CF tells us about the response of the conserved density fluctuations, whereas the spatial CF provides information regarding the mass and width of the particle in a medium. Similarly, the spectral representation of the vector current-current CF is related to the production of lepton pair, which leaves the created medium with minimum interaction. Thus, the dilepton rate (DR) is considered as one of the most reliable probes for QGP.
Being an important probe, the rate of dilepton production in presence of magnetic field has been studied using different techniques. Some of the earliest studies are done by Tuchin [34][35][36] where he used phenomenological input to estimate the DR. To obtain the DR [35,36] a semi-classical Weiszäcker-Williams method [47] is utilised. Employing formal field theoretical technique and using Ritus eigenfunction method DR from a magnetised hot and dense medium was calculated in [37]. In two other articles, where two of the present authors were involved, such formal field theoretical techniques have been used along with Schwinger method [48] to estimate the DR -in one case with imaginary time formalism [38] and real time formalism in the other [39]. The problem was revisited in the latter formalism in [40] for hot QCD matter and significant enhancement was found in the low mass region.
In all such studies the photon self energy 1 is calculated applying a formal field theoretical treatment while using quark propagators modified by the background magnetic field. In this article we would like to estimate the DR in presence of a strong magnetic field within the ambit of effective models such as Nambu-Jona-Lasinio (NJL) model [50][51][52] and its extension with Polyakov loop (PNJL) [53][54][55]. Here we also extend our previous studies [38,39] to incorporate a finite chemical potential in addition to temperature. Now, in NJL model the quark propagator will be modified by an effective mass which is dynamically generated in a magnetised hot and dense medium [56][57][58][59][60][61][62][63][64]. We shall see that this will have some impact on the DR that we wish to evaluate.
Ongoing investigations provide us a strong hint that the QGP created in the HICs is a strongly interacting one [65][66][67][68][69][70][71][72]. Such nonperturbative nature of the medium can be mimicked roughly, through PNJL model where a background gauge field is added to the usual NJL model. The properties of the magnetised medium can be affected by such nonperturbative behaviour. In one of our earlier study we have investigated this point and found that the spectral properties are influenced to quite an extent, specially the DR gets enhanced due to the nonperturbative effect through the Polyakov loop fields [49]. This makes the present study more interesting where the effect of the strong magnetic field will be further influenced by the presence of a background gauge field, specifically through PNJL model augmented with magnetic field [73][74][75][76]. We perform our calculation in the strong magnetic field limit, in which case only the lowest Landau level (LLL) remains significant. This could be relevant for the initial stages of HICs, when the magnetic field is supposed to be very strong.
The paper is organised as follows: In section II we briefly introduce the (P)NJL models. The details of the calculation of SF and DR have been given in section III. Then in section IV we show the results obtained and finally conclude in section V.

A. (P)NJL Model
Here we briefly describe the details of the effective models that we work with, namely the NJL and PNJL models. For the sake of completeness we start by recapitulating the case when the magnetic field is absent. The present work considers a two flavor NJL model. The corresponding Lagrangian is [50][51][52]77] where, m 0 =diag(m u , m d ) with m u = m d and τ 's are Pauli matrices. G S is the coupling constants for local scalar type four-quark interaction. The scalar four quark interaction term leads to the formation of chiral condensate σ = ψ ψ . The thermodynamic potential in mean field approximation is given as: where f being the energy of a quark with flavor f having constituent mass or the dynamical mass M f and Λ being a finite three momentum cut-off. The values of Λ (651 MeV), G S (10.08 GeV −2 ) and m 0 (5.5 MeV) are fitted to reproduce some of the QCD vacuum results [55]. The thermodynamic potential depends on the dynamical fermion mass M f which depends on the scalar (σ) condensates through the relation Let us now briefly discuss PNJL model [53-55, 78, 79], where, in addition to those in the NJL model we have a couple of additional mean fields in the form of the expectation value of the Polyakov Loop fields Φ and its conjugatē Φ. The Lagrangian for the two-flavor PNJL model is given by, where D µ = ∂ µ − igA µ a λ a /2, A µ a being the SU (3) background fields, λ a 's are the Gell-Mann matrices and U(Φ,Φ, T ) is the effective Polyakov Loop gauge potential. The thermodynamic potential can be obtained as The effective Polyakov Loop gauge potential is parametrized [55] as Values of different coefficients a 0 , a 1 , a 2 , a 3 , b 3 , b 4 and κ are same as those given in [79,80]. T 0 is taken as 210 MeV. The Vandermonde (VdM) determinant J(Φ,Φ) is given as [79,81] B. In presence of magnetic field The Lagrangian of NJL model in presence of magnetic field reads as [56][57][58][59][60][61][62][63][64], where / D = γ µ D µ and D µ = ∂ µ − iqA µ with q being the electric charge (q u = 2/3e and q d = −1/3e; e is the charge of a proton) and A µ being the electromagnetic gauge field. The field strength tensor is given as Now in presence of magnetic field the dispersion relation of quarks will be modified to where the magnetic field is taken to be in the z-direction. l is the Landau level (LL) and s is the spin states of the quark which we need to take care in presence of magnetic field. Note that now the particle energy depends on the flavour not only through M f but also through q f . The magnetic field will also modify the integral over the three momenta as, Thus the thermodynamic potential in presence of magnetic field [56][57][58][59][60][62][63][64] becomes Here f , l and s represent the sum over flavour, LL and the spin states, respectively. The magnetic field is applied in the z-direction and its magnitude is B, B = Bẑ. The term B 2 /2 in the thermodynamic potential arises from the background magnetic field. Note that the integral in the second term is ultraviolet divergent. This divergence will be taken care of through the three momentum cut-off (Λ) once we rewrite it in terms of Ω vac in equation 14. Also note that we can no more perform the flavour sum trivially as the expression depends on the magnitude of the quark charge. In fact one can separate out the up and down quark parts, since they couple to the magnetic field with different strength. However, although the quark condensates are different in presence of magnetic field we are still in the isospin-symmetric scenario, i.e., m u = m d = m 0 . Thus the constituent quark masses for up (M u ) and down (M d ) quarks remain the same (M u = M d = M ) [57,82]. It is given as M = m 0 − G S σ u + σ d . With this it is evident that the gap equations for up and down quarks will be the same. So we do not separate them here. Though one can treat them separately to see the isospin-symmetry breaking effects of unequal coupling strengths of up and down quarks to magnetic field [57]. Also note that in presence of magnetic field the dispersion relation in equation 9 depends on the type of flavour only through the strength of the coupling to magnetic field and, expectedly, it will be different for different flavour except for the LLL.
In the same way we can write the PNJL model in presence of magnetic field [73][74][75][76], which reads where / D now also contains background gauge field along with the magnetic field with, D µ = ∂ µ −iqA µ −igA µ a λ a /2. All the other notations carry their meanings as mentioned before. The corresponding thermodynamic potential [73][74][75][76] is Here again the diverging second term will be rewritten in terms of Ω vac in equation 15 to render it finite. We can further simplify this potential using the technique used in [56] and write it in both the conditions -(i) T = 0, µ = 0 and B = 0 and (ii) all of them being nonzero. Here we concentrate only on the second case and the expression in equation 11 can further be written in terms of Ω vac as The corresponding PNJL potential is Before we go to the calculation of vector CF in the effective models we discuss here a few of the intricacies. Some of the earliest studies using NJL model showed that a strong constant magnetic field enhances spontaneous chiral symmetry breaking through the generation of fermion dynamical mass [83,84]. This phenomenon is known as magnetic catalysis (MC) and can be understood in terms of dimensional reduction of the system. This effect of MC was further obtained by many of the effective model studies [57,58,75,82,85,86], including NJL and PNJL. In models like PNJL behaviour of MC is found for both chiral and deconfinement dynamics. It is also found by different lattice QCD simulations that the values of light quark condensates increase at temperature well below and well above the critical temperature (T c ) but it decreases near T c [16,87]. This behaviour of decreasing condensate values with the increase of magnetic field is dubbed as inverse magnetic catalysis (IMC). On the other hand IMC is counterintuitive and is probably occurring due to the dominance of sea contributions over the valence contributions of the condensate around T c [88]. Once IMC is recognised, for both chiral and deconfinement transitions [17,88,89] by LQCD calculation several attempts have been made to understand it through different effective QCD models [18-20, 60-64, 74, 76, 90, 91]. For example the IMC phenomenon is observed in Ref [60] through a magnetic field dependent coupling constant in NJL model, whereas within the same model similar behaviour is achieved in a strong magnetic field as the coupling constant becomes anisotropic and the quarks produce antiscreeing effect at LLL [61]. At this point we also want to note that in the present work we will not be discussing about the effect of IMC because we will not be focussing in temperature ranges close to T c , where IMC effects become dominant. Hence this study will only incorporate MC effects, as we will see later in our results.

A. Vector correlation function in one-loop
In this section we would like to compute the current current correlator in NJL and PNJL models in presence of a strong external background magnetic field. The two point current-current correlator C µν (Q) is related to photon self-energy Π µν (Q) as with q f is the electric charge of a given quark flavour f . The photon self energy in one loop level (figure 1) is again expressed as where Q is the external momentum, P and K = P + Q are the loop momenta. Tr c represents both color and Dirac traces whereas the f is over flavor because we have considered a two-flavor system (N f = 2) of equal current quark mass (m f = m u = m d = 5 MeV if not said otherwise). Finally the electromagnetic spectral representation is extracted from the imaginary part of the correlation function C µ µ (Q) as To write the fermionic propagator in presence of a background field we use the prescription provided by Schwinger [48]. When the background magnetic field is very large (q f B → ∞), the only relevant LL is the lowest one. This point is well described through figure 1 as in [38]. In the strong field approximation the fermion propagator gets simplified to where P is the four momentum and m f and q f are the current mass and the absolute charge of the quarks, respectively. In presence of strong magnetic field there will be dynamical mass generation and the m f will be replaced by the effective quark mass. This is incorporated here through the effective model calculation. The notations that we use in equation 19 and for the rest of the paper are as follows where and ⊥ represent parallel and perpendicular components, respectively. Now the NJL quark propagator in Hartree approximation is written as where M f is the effective mass already defined in section II; its magnetic field dependence, which arises through the gap equation, is now explicitly shown. We will elaborate the calculation for NJL model and the extension to PNJL one will be straightforward. Now we rewrite equation 19 when the fermion mass is an effective one So following equation 17 the current-current correlator can be computed as where 'sfa' indicates the strong field approximation and Tr represents the trace in the Dirac space only. The longitudinal and transverse parts are now completely separated and we can perform the Gaussian integration over the transverse momenta, which leads to .
The tensor structure S µν originates from the Dirac trace and is given by where the Lorentz indices µ and ν can only take longitudinal values.
To evaluate the spectral function, we can contract the indices µ and ν in equation 24 which leads to a further simplification as .

Magnetised hot medium
First we describe a strongly magnetised hot medium. To introduce temperature we use the imaginary time formalism technique and thus replace the p 0 integral by Matsubara sum as . (27) We should note here that the effective mass is now also dependent on temperature along with the magnetic field. One can always use the contour integration method to perform the Matsubara sum, but it becomes increasingly difficult when the number of propagators present in a diagrams increases. There is an elegant way to perform such complicated frequency sum using so called Saclay method. The essential trick, introduced by Pisarski [92], is to use propagator that are in co-ordinate representation in time but momentum representation in space: and where E k = k 2 3 + M f (T, B) 2 and n F (x) = 1/(exp(βx) + 1) is the Fermi-Dirac distribution function with β = 1/T . Using these, equation 27 can be simplified as We now perform the τ integral to get where the factors r and l can be related with the particles and antiparticles thereby representing various physical processes. The discontinuity can be computed using which leads to The general form of the delta function in equation 33 corresponds to four processes [38], depending on the values of r and l. Among them, r = 1 and l = 1 corresponds to a process where a quark and a antiquark annihilate to a virtual photon, which is the only allowed process. So, for this case, one can write from equation 33, After performing the p 3 integral, the spectral function in strong field approximation is obtained as where with the superscript s representing the strong magnetic field approximated result.

Magnetised hot and dense medium
In this section we incorporate the effect of baryon density as well in order to describe the general case of a strongly magnetised hot and dense medium. Here the effective mass further becomes dependent on chemical potential along with T and B through the gap equation. In presence of chemical potential the NJL propagator (21) gets modified to So the Matsubara sum includes chemical potential and gets modified to This replacement eventually will have effects only in the distribution function. Thus equation 35 gets modified to where n ∓ F (q s ) = n F (q s ∓ µ). We use equation 41 to calculate SF and also use it to get DR in the ambit of NJL model. It expectedly reduces to equation 35 for µ = 0.

Magnetised hot and dense medium in presence of Polyakov loop field
We now want to extend our study by including a background gauge field (chromo). This is done through the PNJL model for which the effective propagator reads as where A 4 is the background gauge field. One crucial difference with the NJL model is that now the effective mass also varies with A 4 . This dependence is denoted through the Polyakov loop (Φ). The whole procedure of calculation will remain the same, but now with A 4 the colour trace will be nontrivial. In usual NJL model the colour trace just gives a factor of N c . Whereas in PNJL model the colour trace modifies the usual Fermi-Dirac distribution functions. The distribution functions of particle and anti-particle get modified to [49,80] and respectively. f − and f + represent the distribution functions for particle and antiparticle respectively. Thus equation 35 gets modified to With Φ =Φ = 1 we get back the usual Fermi-Dirac distribution function. This essentially happens when we go to a temperature much higher than the critical temperature. For µ = 0, we have Φ =Φ which implies f − = f + . So, now we have spectral function for three different cases: (a) only magnetised hot medium denoted as ρ 1 V , (b) magnetised hot and dense medium is denoted as ρ 2 V and (c) magnetised hot and dense medium in presence of Polyakov loop field is denoted as ρ 3 V .

B. Dilepton rate
Once we have the vector spectral function (ρ V , in general) we can calculate the dilepton rate. Now in presence of magnetic field there are three possibilities: (a) only the initial quark-pair is moving in presence of the magnetic field, (b) both the initial quark-pair and the final lepton-pair are moving in the magnetised medium and (c) only the final lepton-pair gets influenced by the magnetic field. We will consider here both the cases (a) and (b). Case (a) is the most relevant scenario for a non-central heavy-ion collisions where the magnetic field is decaying very fast [93,94].
For case (a) the dilepton production rate for massless (m l = 0) leptons with N f = 2 can be written from [38] as Similarly for case (b) the dilepton production rate with N f = 2 can be written from [38] as Detailed discussions on different scenarios can be found in these articles [35,38].

A. Effective models
We have the plots for the constituent mass for different values of B in the nonzero quark mass limit for NJL model in figure 2. In the left panel 2(a) we have shown how the constituent mass increases with the increase of the magnetic field (i.e. Magnetic Catalysis (MC) effect). Depending on the value of the magnetic field generally one needs to sum up appropriate number of Landau levels. As in our main results of spectral properties, we will be using LLL approximation; in the right panel 2(b) we have also shown the validity of the LLL approximation as a function of T for a given magnetic field. The choice of LLL and its validity are discussed in numerical details at the appropriate places.

B. Spectral function in a strongly magnetized medium within the ambit of effective models
Next we move on to the main results of the present manuscript, i.e. the electromagnetic spectral properties in a strongly magnetized medium. In figure 5 the variation of the electromagnetic vector spectral function in NJL model is shown as a function of the invariant mass M of the photon for different values of magnetic field with a fixed value of temperature and chemical potential. For the left panel 5(a) we have considered vanishing chemical potential, whereas in the right panel 5(b) the effect of chemical potential on the spectral function has been shown for fixed values of both temperature and magnetic field. From fig 5(a) we learn that for given values of T and µ the strength of the spectral function increases as we increase the magnetic field. Similarly it can be seen from figure 5(b) that for a fixed value of T , with increased value of µ the strength of the spectral function have decreased. This is expected as the finite value of chemical potential adds to the medium effect. Similar plots for PNJL model is shown in figure 6 where the effect of the chemical potential is seen to be negligible for the range of temperature that we have investigated for. Also one can notice the relative small increase in the value of the electromagnetic spectral function for the PNJL model as compared to the NJL one. This increment from NJL to PNJL can be even higher depending on the range of temperature that one deals with. Now, since the spectral function calculation is done within the LLL approximation, the effective model inputs should comply with that assumption. As evident from figures 2, 3 and 4, for models the LLL approximation is not valid for all temperatures. So depending on the value of temperature that we work with, it might introduce some error. For example, in figure 5(a) for the lowest considered magnetic field (15 m 2 π ), the LLL approximation introduces a maximum error of roughly 12%. The error gets reduced to ∼ 7% for eB = 21 m 2 π . These errors for NJL model can be further reduced by choosing appropriate temperature and chemical potential; for example the error gets reduced to below 5% at T = 130 MeV, eB = 21 m 2 π and µ = 0 and it can be lowered by further reducing the temperature. The value of temperature for NJL model is so chosen that it can be compared with the PNJL one for which the maximum error from figure 6(a) is ∼ 1% for eB = 15 m 2 π and becomes even lesser with the increase of magnetic field. Thus, for the same temperature the validity of LLL approximation is better in PNJL model as compared to the NJL one. This happens because of the presence of background gauge field which mimics the QCD more closely by emulating the confinement effect (statistically).
Finally if we compare our present results (in light of effective models) with the previous results obtained in Ref. [38] considering just the current quark mass, one can immediately notice a vast enhancement in the magnitude of the electromagnetic spectral function as we had expected -it is roughly of the order of 10 4 . This can be understood in the following way. If we compare equation (26) from [38] with equation (35) in this article, then it is clear that  the major change is in terms of the mass. The effective mass is almost two orders of magnitude higher than the current quark mass that we used in the previous work. There are other changes also, for example the change in the distribution function, which is now dependent on the external magnetic field; but those are small as compared to the change in the value of quark mass.
C. Dilepton production rate in a strongly magnetized medium within the ambit of effective models After the evaluation of the electromagnetic spectral functions in view of effective models, now we present our results for the dilepton production rate. As mentioned earlier, we will be dealing with two separate cases; firstly when only the initial quarks are affected by the magnetic field, which can be applicable for a fast decaying external magnetic field or the dileptons produced at the edge of the fireball. And secondly the more general case of both the initial quarks and the produced lepton pairs being influenced by the external strong magnetic field. Below we discuss both the cases. rate also increases. In the right panel 7(b) the suppressing effect of chemical potential has been shown, which can be mapped from figure 5(b). Similar plots for PNJL are shown in figure 8. The nature of all the plots can be easily understood extrapolating from the spectral function results. As an overall comparison with the existing results we found the following. In reference [38], where this case was shown considering the current quark mass, the enhancement in the dilepton production rate (i.e. the ratio dR/dR Born > 1 ) was for a very small range of invariant mass (M < 100 MeV) which was insufficient to be detected in any of the heavy ion collisions experiment. Here for NJL model we can see that the enhancement is observed for a range of invariant mass starting at M ∼ 800 MeV and for PNJL one at M ∼ 850 MeV, which is certainly more compatible with the presently ongoing HIC experiments. The reason for such shifts in the values of invariant mass for which the DR is enhanced is the increase in effective mass which is dynamically generated in presence of magnetic field. After a certain value of invariant mass the DR falls off very fast because of LLL approximation [38]. While calculating the DR there is no more source of errors other than that already present in the spectral function. The maximum error possible in figure 7(a) is 8% while that in figure 8(a) is less than 1%. For the temperature that we worked with the error in NJL model is larger than PNJL one. PNJL model is more relevant as far as the full QCD is concerned, thus the results in PNJL model should be taken more seriously. 2. Both the initial quarks and the final lepton pairs are influenced by the magnetic field One would naively expect that the dynamical mass generation will also lead to the enhancement of DR when both initial quark pair and final lepton pair are influenced by the magnetic field. But this does not happen. In figure 9 we have shown the plots for both dimuon and dielectron in both NJL (9(a)) and PNJL (9(b)) models. First of all, in this case the DR is proportional to B 2 -where one B comes from the leptonic part and the other from the electromagnetic spectral function involving quarks (vide equation 43). Along with M f we have now another threshold in terms of the leptonic mass (m l ) and essentially the effective threshold is determined by max(M f , m l ). We also observe that now the DR is proportional to m 2 l . The mass of the lepton is expected to be affected by the magnetic field, but we are in the strong field regime (LLL approximation) and the mass correction from the magnetic field vanishes. So we use the known vacuum values of masses for leptons. Because of the presence of these multiple factors which act opposite to M f , the DR (ratio) in this case always remain below unity. The dimuon rate is ∼ 10 4 times higher than the dielectron one, since the rate is proportional to m 2 l . The discussion on error calculation done in subsection IV C 1 remains valid for these results as well.

V. SUMMARY AND CONCLUSION
In this work we have calculated the electromagnetic spectral function (SF) by evaluating the one-loop photon polarisation tensor for a strongly magnetised hot and dense medium. Because of the consideration of high enough magnetic field we use the lowest Landau level (LLL) approximation. The whole calculation is performed in the ambit of mean field models namely Nambu-Jona-Lasinio (NJL) and its Polyakov loop extended version (PNJL). These models allow for a dynamical generation of a medium (i.e. T , µ) dependent quark mass which further gets affected in the presence of a magnetised medium. This is shown here as magnetic catalysis (MC) for chiral transition in both NJL and PNJL models. The presence of magnetic field also influences the deconfinement dynamics -a catalysing effect for deconfinement transition is also observed, although the effect is milder. We have also discussed both the importance and the effects of the Vandermonde (VdM) term in the pure gauge potential in presence of magnetic field. The value of the Polyakov loop, specially the form that we work with, overshoots unity for higher values of temperature in absence of the VdM term, which should not be the case. The inclusion of VdM term further shifts the critical temperature for both chiral and deconfinement transitions to higher values.
We have argued that for the values of the magnetic field considered here the LLL approximation is not strictly valid above certain respective temperatures. Keeping the temperature within the validity region of LLL we have evaluated the electromagnetic correlation function and dilepton rate (DR) using the dynamically generated magnetic field dependent effective mass of quark in the loop where LLL approximation is used. It is found that the strength of the SF gets boosted, predominantly because of the presence of dynamical quark mass. The difference in the orders of magnitude between the current quark mass and the dynamically generated effective mass gets reflected in the vast enhancement of the SF in our newfound results, when compared with the existing ones. As we increase the strength of the magnetic field for a given value of temperature the strength of the SF increases, for both NJL and PNJL models. Being proportional to the SF this increment is further reflected in the DR which is enhanced as compared to the Born rate for certain values of invariant masses. For non-zero chemical potential the value decreases in NJL model due to increased opposing medium effect, whereas no significant change is observed in case of PNJL model for the range of temperatures that we have investigated for. Our most significant observation is the shifting of the enhanced DPR from an experimentally insignificant (M < 100 MeV, in our previous work [38]) to an experimentally favourable range of invariant mass (M ∼ 800 − 900 MeV), again mainly due to the magnetic field dependent dynamically generated effective mass.
It is important to mention here again that we have worked with the usual (P)NJL models which give rise to MC and we do not explicitly discuss the issues of IMC effects, mainly concentrating on the evaluation of the spectral function vis-a-vis dilepton rate. This is reasonable as long as we are not close to the transition temperature, because away from the T c for both T T c and T T c one gets MC [16,87,89]. It has also been established [96] that introducing a new interaction term in the Lagrangian which preserves chiral symmetry and rotations along the magnetic field direction subsequently increases the effective dynamical mass in LLL. This would account for an exciting future prospect in view of spectral properties which, as we have shown in the present context, are heavily dependent on the dynamical mass. Some of the earlier studies [97,98] have also shown that in presence of baryonic chemical potential at zero or low temperatures, strong external magnetic field usually favors the formation of spatially non-uniform inhomogeneous condensates in the form of a dual chiral density waves (DCDW). As we have dealt with finite baryon chemical potential in this study, it would be interesting to see how the DCDW affects the spectral properties, both at zero and finite temperatures. Finally, we note that the LLL approximation valid for a high enough magnetic field may be realised only in limited cases in the present HIC experiments. To generalize our claim, the DR presented in this paper should be evaluated taking all the Landau levels into account, where the restriction on the temperature range may also be avoided. Work in this line is in progress.

VI. ACKNOWLEDGEMENT
CAI would like to thank Rishi Sharma for many fruitful discussions. He would also like to thank Department of Atomic Energy, India for the financial support. He acknowledges the facilities provided by Variable Energy Cyclotron Centre, India where he was a visiting fellow when the preliminary parts of the work were completed and also his present institute, Tata Institute of Fundamental Research, India. AB is supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Govt. of Brazil, under grant CAPES.