Holographic model for light quarks in anisotropic hot dense QGP with external magnetic field

We present a five-dimensional twice anisotropic holographic model supported by Einstein-dilaton-three-Maxwell action describing light quarks. The first of the Maxwell fields provides finite chemical potential. The second Maxwell field serves for anisotropy, characterizing spatial anisotropy of the QGP produced in heavy-ion collisions (HIC). The third Maxwell field is related to a magnetic field that appears in HIC. The dependence of the 5-dim black hole solution and confinement/deconfinement phase diagram on this magnetic field is considered. The effect of the inverse magnetic catalyses for light quarks phase diagram is obtained. Positions of critical end points are found. We also study the behavior of the conductivity for light quarks in both isotropic and anisotropic cases and show that behaviour of the conductivity near critical points essentially depend on quark masses, meanwhile at high temperature they are similar.


Introduction
Holographic duality provides an effective approach to study nonperturbative properties of QCD [1][2][3].One of the achievements of the QCD holographic description is that holographic QCD (HQCD) models can describe QCD at all energy scales [4,5].HQCD reproduces properties of QCD obtained by others methods, namely perturbative renormalization group and lattice approach.However neither of these methods can describe QCD with non-zero chemical potential and therefore predict the phase diagram for quark-gluon plasma (QGP).Nowadays this goal is achieved by holographic methods only.There are HQCD models that describe QCD with heavy quarks [6], light quarks [7] as well as with realistic quarks contents [8].These models are based on the so-called "bottom-up" holographic approach.Holographic method can also describe anisotropic dense matter in external magnetic field.
The goal of this paper is to study the behaviour of anisotropic QCD with light quarks in strong external magnetic field.This paper is the generalization of our consideration of anisotropic QCD with heavy quarks in external magnetic field [9] on the one hand and of our consideration of anisotropic QCD with light quarks at zero magnetic field [10] to the case of non-zero external field on the other hand.
The model under consideration [11] has been built using a potential reconstruction method [1-7, 9, 10, 12-31].This model is characterized by a 5-dimensional Einstein-dilaton-Maxwell action and a 5-dimensional metric with a warp factor that tunes the blackening function and, consequently, the thermodynamics of the model.The phase transition structure of HQCD strongly depends on the warp factor choice.The simplest warp factor has an exponential form with a quadratic polynomial [6,32].It reproduces the first-order phase transition from lattice calculations (Columbia plots [33,34]) for the heavy quarks.To reproduce the first-order phase transition for the light quarks one has to use rational functions [5,7].
To describe the QGP produced in heavy ion collisions (HIC) by holographic methods one has to make more modifications of the 5-dimensional metric.These modifications are related to two types of the anisotropy that take place in HIC: i) anisotropy between the longitudinal and transversal directions and ii) anisotropy in the transversal plane related to magnetic field produced in non-central collisions of ions.It was proposed to characterize the first type of anisotropy by a parameter ν that is similar to the Lifshitz parameter of the AdS metric modification.The isotropic model corresponds to ν = 1.For the value of about ν = 4.5 the energy dependence of the entropy produced in the HIC agrees with the experimental data [35] for the energy dependence of the total multiplicity of created particles [36].The magnetic field created by non-central ion collisions [37] produces anisotropy in the plane transverse to the ion collision axis.Therefore, it is important to consider these types of anisotropy in holographic models [9,23,[38][39][40][41][42][43][44][45][46].To construct anisotropic models, one considers an additional Maxwell field to support the anisotropy in metrics [9, 10, 12, 14, 15, 18-20, 22, 23, 25, 31, 36, 47].Note that calculations in the twice anisotropic model for the light quarks are much more complicated than for the heavy quarks [9,23], since the warp factor and the kinetic gauge function f 1 have more complex form.For the heavy quark model we have found [9] the effect of the inverse magnetic catalysis (IMC) -larger absolute values of the coupling coefficient c B lead to decrease of the transition temperature.In the model for the light quarks considered in this paper we also find the IMC.This work is a generalization of results obtained in [9,10] to the case of twice anisotropic holographic light quark model.The results obtained in [10] are generalized to the case of present of an external magnetic field, and the results obtained in [9] to the case of light quarks.
We also study electrical transport properties of QGP and its relation to the thermal direct photon production in the line of previous studies of these problems for other holographic models [48][49][50][51][52][53][54][55].Photons do not interact with QGP hadronic matter and therefore provide information about various characteristics of QGP including phase structure at different time scales.It is well known that the direct photons (DP) emission rate is connected to the conductivity of QGP [56,57].In this work we generalize results of [31] concerning the conductivity and photon production to the light quarks case.We show that conductivity both in longitudinal and transverse directions essentially depends on anisotropy parameter ν.The DP flow also depends on temperature and other thermodynamic properties of QGP.Like it was in the heavy quarks case, for the light quarks the higher the magnetic field or/and chemical potential, the higher the conductivity.We will see that for ν = 1 the conductivity behavior for the heavy and light quarks is different at near-critical temperatures, but at high temperatures they both saturate to constant values.For ν = 4.5 the conductivity behavior for the heavy and light quarks is also different at near-critical temperatures as well as at high temperatures.Conductivity for the light quarks model along the collision direction monotonically decreases with temperature while transversal components have minimum and increase at large temperatures.Since we expect isotropisation [18,58,59], i.e. ν → 1 in about 1 − 5f m/c ∼ 10 −24 s, we expect essential change of conductivity at this time scale.
The paper is organized as follows.In Sect.2 the holographic model is presented.In Sect.3 the thermodynamics of the model and behaviour of the Wilson loops are described, the confinement/deconfinement phase diagram is obtained.Sect.4 describes the derivation of the direct photons and electrical conductivity properties.In Sect.5 we summarize our results and discuss future directions of investigations.

Model
In this section we present a model suitable for the light quarks' version of the confinement/deconfinement phase diagram [10] in magnetic field similar to how it was done in our work [9].

Metric and EOM
Here we take the same action as in our previous work [9].In Einstein frame it has the following form (2.1) Here we use the anzats, where the non-zero components of the electro-magnetic field and the field strengths are In action (2.1) φ = φ(z) is the scalar field, f 1 (φ), f 2 (φ) and f B (φ) are the coupling functions associated with the Maxwell fields A µ , F µν and µν correspondingly, q and q B are constants and V (φ) is the scalar field potential.Thus (2.1) is the extended version of the action used in [10,14], where we add an external magnetic field We also use the same metric anzatz [9]: Let us make some notes on the main parameters of metric (2.3).The difference between "heavy quarks" and "light quarks" cases lies in the form of the warp factor b(z).For the heavy quarks we used b(z) = e − cz 2 2 [14].To get the "light quarks" version we follow [7] and assume A(z) = − a ln(bz 2 + 1).
Here parameters ν and c B have the same meaning as in [9].The primary anisotropy parameter ν describes non-equivalence between longitudinal and transversal directions.The isotropic case corresponds to ν = 1, while ν = 4.5 reproduces the multiplicity of the charged particles production [36].
We refer to c B as to the magnetic field parameter that describes non-centrality of the heavy ion collisions.Although explicit relation of c B to the magnitude of magnetic field is not established, there are considerations suggesting that B 2 = − c B [44,60].So we have c B < 0 in all our calculations.Magnetic field parameter c B is an implicit function of the magnetic charge q B from the anzats (2.2).However, we will not need the explicit form of this relation in what follows.We usually set q = q B = 1 in numerical calculations.The metric anzats (2.3) also contains the blackening function g(z) which will be determined in the subsequent paragraphs.
Applying the stationary action principle, we get the same equations of motion (EOM) as in [9]: (2.11) Turning off the external magnetic field, i.e. putting c B = q B = f B = 0 into (2.5)-(2.11),we get the EOM from [10].Normalizing to the AdS-radius, L = 1, we get the EOM from [14].Excluding anisotropy, i.e. putting ν = 1 and f 2 = 0, we get the expressions that fully coincide with the EOM from [6,7].Thus (2.5)-(2.11)are universal anisotropic EOM, appropriate both for heavy and light quarks description, that include solution from [6,7] as an isotropic limit.We also consider the general form of the boundary conditions: (2.12) where z 0 = 0 corresponds to [7] and z 0 = z h to [14].The choice of the boundary condition for the scalar field was discussed in details in [10].

Solution
Just as it was in previous cases, to solve EOM (2.5)-(2.11)we need to determine the form of the coupling function f 1 .Let us take the same form as for the light quarks model [10].Taking into account the "light quarks" warp factor, we get (2.17) We take a = 4.046, b = 0.01613, c = 0.227 to make our solution agree with the results from [7] in the isotropic case, the results from [10] in the anisotropic case and the results from [9] for the heavy quarks version.These values are due to the mass spectrum of ρ meson with its excitations and lattice results for the phase transition temperature [7].
To obtain the blackening function we need to know the kinetic function f B describing coupling to the third Maxwell field F (B) µν .Equation (2.5) is a consequence of the remaining equations (2.6)-(2.11),which is checked directly.Subtracting (2.9) from (2.10), we get an expression for the function that allows to transform (2.7) into  and get the solution (2.20) However this solution can't be accepted as the final one and should be improved.First of all it can't be shown that our previous "light quarks" blackening function [10] serves as zero magnetic field limit for (2.20).To overcome this difficulty we simplify the expression via partial integration: where I 1 (z) is described by the corresponding expression from (2.20), but I 2 (z) was redefined as Further solving (2.5)-(2.11)system with the coupling functions (2.15) and (2.18) and the boundary condition (2.14) gives (2.25) Note, that we do not fix coupling function for the 3-rd Maxwell f B , but derive it from the EOM with intent.Fixing, for example, f B = f 1 makes system (2.5)-(2.11)not selfconsistent, therefore proper solution with the f B = f 1 condition can't be found.

Temperature and entropy
For the metric (2.3) and the light quarks warp factor temperature and entropy can be written as: On Fig. 2 temperature as a function of horizon is presented.The multivalued behavior of T (z h ) preserves starting from the chemical potential value of the so called critical end point (CEP) till another specific chemical potential value for which second horizon appears.For ν = 1 µ CEP > 0, while in cases of primary anisotropy (ν = 1.5, 3, 4.5) µ CEP = 0 for small magnetic fields.Magnetic field's increasing leads to the growth of µ CEP values and decrease in characteristic temperatures of the multivalued behavior, i.e. the temperatures of the background phase transition (corresponding to black hole -black hole collapse in AdS 5 ).Thus in the magnetic field strong enough background phase transition should occur at high chemical potentials and negligible temperatures.

Background phase transition
This is exactly the picture we observe on phase diagrams Fig. 3, obtained via free energy consideration   where we normalize the free energy to vanish at z h → ∞.The background phase transition curves become shorter with the magnetic field increasing (larger c B absolute values).Note, that two opposite tendencies can be seen: growth of the µ CEP and decrease of the maximum chemical potential value.For c B ≈ (− 0.05; − 0.03) depending on ν these tendencies come into conflict, and CEP curve reverses it's movement back towards lower values of the chemical potential (Fig. 4A).
Primary anisotropy affects the background phase transition for the light quarks QGP in much the same way as it did for the heavy quarks [9]: lowers transition temperatures and enlarges chemical potentials avaliable.But it almost doesn't influence the limit c B values at which phase transition still exists.Phase transition temperatures reach negligible values at c B ≈ − 0.1 almost regardless of primary anisotropy ν (Fig. 3A-D, Fig. 4B).It's effect on c B for which µ CEP becomes non-zero is more notable: c B < − 0.001 for ν = 1.5 and ν < − 0.005 for ν = 4.5 (Fig. 4A,C).
We present CEP for ν = 1, 1.5, 3, 4.5 in a 3-dimensional space of parameters (µ, c B , T ) on Fig. 4D.Thus, figures A, B and C are the projections of a complete picture D on the corresponding planes.Also, one can see the smooth peach coloured   surface of the first order phase transitions interpolating for intermediate values of anisotropy.The orange curve indicates boundaries of possible phase transition locations in this space of parameters.

Temporal Wilson loops
The expectation value of the temporal Wilson loop calculation does not differ from the heavy quarks case [9] untill we substitute the specific light quarks warp factor into the general expressions for the dynamical wall (DW) equations: Primary anisotropy has significant influence on the Wilson loop transformations in the magnetic field (Fig. 5).For ν = 1 the Wilson loop phase transition curve becomes less inclined and it's temperature decreases, while maximum chemical potential grows.At c B = − 0.0000985 the curve becomes horizontal and for c B < − 0.0000985 it just disappears (Fig. 5A).For ν = 4.5 larger absolute c B also leads to lesser temperatures, but chemical potential decreases as well, so the Wilson loop curve shrinks clinging to the abscissa axis µ at the end.The curve disappears for c B < − 0.0858 (Fig. 5D).For intermediate values of primary anisotropy ν = 1.5 (Fig. 5B) and ν = 3 (Fig. 5C) a mixture of both these tendences can be seen.But inverse magnetic catalysis takes place for any primary anisotropy.
Let us now consider the combination of the 1-st order phase transition and a crossover for the primary anisotropy and magnetic field fixed (while both transition curves do exist).On Fig. 6 regions of their intersection or approach are shown, and the general view of the phase diagrams can be found on the inserted mini-plots.Points of intersections, where the confinement/deconfinement main role passes from the crossover to the 1-st order phase transitions (or vice versa), are marked as (µ by 2 , T by 2 ) (Tab.1).Critical end points, where the 1-st order phase transition begins, are marked as CEP LQ .
For the primary isotropic case ν = 1 a crossover disappears very soon, so it's position relative to the 1-st order phase transition does not have time to change much.The intersection point just shifts to the larger µ and the lesser T .
For ν = 1.5 both curves are very close to each other and look cocentric.They do not intersect, so the crossover defines the phase transition for all chemical potential values.Magnetic field makes these curves closer to each other and hardly distinguishable, but the intersection appears at about c B = − 0.018 and coincides with the CEP LQ .At this moment phase transition main role clearly passes from the crossover to the 1-st order phase transition though for rather a short interval.Before it's disappearance at c B = − 0.01882 the crossover line shifts down to lesser temperatures and becomes almost horizontal like it was for ν = 1.
For ν = 3 the 1-st order phase transition and a crossover can't be distiguished properly.Magnetic field shifts the 1-st order phase transition curve under the crossover, so that they do not intersect.Under these circumstances there is no smooth transfer between the crossover and the 1-st order phase transition, but a jump to lower temperature -temperature of CEP LQ .Near the limit c B = − 0.0631 the crossover catches up with the 1-st order phase transition and a smooth transfer appears again just before the crossover vanishes.

Conductivity in anisotropic background
We use the model derived above to investigate properties of the direct photons formed in HIC.The detailed consideration of the holographic approach to electrical conductivity and direct photon emission rate of QGP is presented in [31,48,61].We briefly describe the framework for the holographic investigations of the electrical conductivity below.
In the anisotropic background1 we consider a probe Maxwell field We assume this creates no gravitational backreaction.Function f 0 (φ) is called a gauge kinetic function and is used to fit the model with lattice data for the isotropic QGP at zero chemical potential and zero magnetic field [48].
Transport coefficients may be found within the linear response theory using the Kubo relations.Namely, the electrical conductivity is related to the retarded Green's functions of electric currents.Therefore, one is interested in a plane wave solution to (4.2), which will be holographically dual to U (1) electric current J µ of the boundary theory.We take the ansatz of a plane wave propagating along the x 3 -direction.The wave amplitude is a function of the "holographic coordinate" z only.The Green's function in a low frequency limit may be found either using the prescription by Iqbal and Liu [61] or Son and Starinets [62].The equivalence of these approaches is demonstrated in [31].The components of electric conductivity tensor are then ) where z h is the event horizon position [31].

Model tuning
To tune the model parameters we use the results of lattice calculations from [63].These results are known for the isotropic QGP (ν = 1) with a number of flavours N f = 3 at zero chemical potential (µ = 0) and without magnetic field (c B = 0).Thus, we fit the behaviour of electrical conductivity to the lattice data in this simplest case using the gauge kinetic function f 0 (φ).Then it is possible to calculate the conductivity in anisotropic case and for non-zero parameters.We take the gauge kinetic function in the following form where C em = 2e 2 /3 is electromagnetic constant for N f = 3 and φ is the dilaton field.
Numerical calculation for the isotropic conductivity and the gauge function are shown on Fig. 7.Here the critical temperature is T c = 0.155 GeV.It was calculated as a temperature of the BB phase transition for a given set of parameters.The conductivity increases with temperature growth and saturates some value close to the N = 4 SYM plasma conductivity.

Isotropic case
The conductivity tensor in the isotropic (ν = 1) case with zero magnetic field reduces to the single value σ = σ ii (i = 1, 2, 3), so the QGP provides equal resistance to direct photons in all directions.The dependence of the conductivity to the temperature ratio σ ii /T on the temperature for different values of chemical potential is shown in Fig. 8A.The non-zero magnetic field kills the degeneracy between two transverse directions (Fig. 8B).The plots presented in Fig. 8 demonstrate that the ratios σ ii /T grow with the normalized temperature till some temperatures around 2T c .The precise location of the maximum depends on the chemical potential and the magnetic field values.At high temperatures the σ/T ratio decreases to some constant values near the conductivity of N = 4 SYM plasma.Plots on Fig. 8A and 8B demonstrate the growth of conductivity with the magnetic field and/or chemical potential growing for the same temperature values.Also, one can see that in the directions orthogonal to the magnetic field the conductivity is lower than in a parallel one (brown curves are lower than the green ones of the same thickness).Therefore, QGP is more opaque along the external magnetic fields.

Anisotropic case
In this subsection we study the electrical conductivity of QGP for the anisotropy parameter value ν = 4.5.The results of numerical calculations of σ 33 /T (conductivity along the external magnetic field) for different values of magnetic field and chemical potential are presented in Fig. 9.
We see that at the ratio σ 33 /T decreases with temperature starting from T = T c .The minimum position depends on the chemical potential and the external magnetic field value.The σ 33 /T ratio after the minimum is greater for zero chemical potential and for lower magnetic fields.Stronger magnetic field and larger chemical potential increase the conductivity at the fixed temperature.One can see that at high temperatures the curves in Fig. 9 are getting closer to each other and almost coincide.Thus, anisotropic QGP with ν = 4.5 along the external magnetic field direction is the most opaque at critical temperatures, becomes more transparent around T ≈ 2T c for µ = 0, Fig. 9A, and T ≈ 3.5T c for µ = 1, Fig. 9B, and after these values becomes more opaque with increasing temperature again.
Thermodynamic properties of σ 22 /T are similar to those of σ 33 /T and are shown in Fig. 10.We see that near the phase transition point σ 22 /T decreases to a minimal value, and then starts to increase linearly at high temperatures.Plots in Fig. 10A,B show that higher the chemical potential is, higher the ratio σ 22 /T becomes.However, the growth rate at high temperatures is larger for µ = 0 case.The curves for   strong and weak magnetic fields are almost indistinguishable, but the conductivity in stronger magnetic fields is greater than in weak ones.For zero magnetic field there is no difference between σ 22 and σ 33 , it can be seen in non-zero magnetic fields only.The value of σ 33 /T is always larger than σ 22 /T .The comparison of these two conductivity components can be found on Fig. 11.
The results of similar calculations for σ 11 /T are shown in Fig. 12.We see that σ 11 /T monotonically decreases to zero with temperature growth, while magnetic field and chemical potential push it down.Or, in other words, QGP along the collision line near the critical temperature is almost opaque and becomes transparent at high temperatures.
On Fig. 13 one can see the change of the conductivity behavior with the primary anisotropy parameter increasing.The green curves denote the isotropic case ν = 1 and for the anisotropic cases the thicker line corresponds to the larger ν = 1.As anisotropy increases, qualitative changes occur in all directions, where as usually two transverse and a longitudinal directions differ crucially.At low energies conductivity in longitudinal direction increases at low temperatures and starts decreasing after approximately T = 2.5T c .At high values of anisotropy parameter ν longitudinal conductivity decreases monotonically.On the contrary, transverse components increase at high temperatures.

Conclusion
In this paper we have considered a twice anisotropic holographic model for the light quarks [11].It is characterized by the Einstein-dilaton-three-Maxwell action and the 5-dimensional metric with the warp factor that has been considered in isotropic case for the light quarks [7].Our model has two different types of anisotropy mentioned above: the anisotropy related to the parameter ν (responsible for reproduction of the experimental energy dependence of multiplicity), and the anisotropy related to the external magnetic field (described by the third Maxwell field).The twice anisotropic model for the light quarks is much more complicated than the twice anisotropic model for heavy quarks [9,23], since the deformation factor and the kinetic gauge function f 1 have more compound form.
We have investigated the dependence of the phase transition structure on the anisotropy, characterized by parameter ν and magnetic fields.The phase transition structure is set up by the mutual arrangement of two types of the phase transition lines -the 1-st order phase transition, originating from the metric background, and the crossover, originating from the temporal Wilson loop behavior.We have got the following results concerning the phase transition structure.
• For ν = 1 the 1-st order phase transition in the light quarks model lasts from (µ CEP , T CEP ) to (µ max , 0).Magnetic field c B = 0 shortens the 1-st order phase transition, increasing µ CEP and decreasing µ max .In the heavy quarks model [9] the 1-st order phase transition lied from (0, T max ) to (µ CEP HQ , T CEP HQ ) almost horizontally.Magnetic field shortened this line just shifting µ CEP HQ to zero.For the light quarks model the 1-st order phase transition line disappears at c B about 7 times larger than it was for the heavy quarks model.Inverse magnetic catalysis is found for both models.
• For ν > 1 primary anisotropy decreases temperature of the 1-st order phase transition and increases µ max .It almost doesn't influence the c B value that the 1-st order phase transition line disappears at, like it was for the heavy quarks model [9].
• For ν = 1 magnetic field decreases the Wilson loop curve inclination until it becomes horizontal.This leads to the inverse magnetic catalysis effect for µ near zero and the direct magnetic catalysis for larger chemical potential values.
In the heavy quarks model [9] the Wilson loop curve underwent the inverse magnetic catalysis only and shrinked tending to vertical T -axis.
• Primary anisotropy ν > 1 decreases the Wilson loop phase transition temperature and the region of the direct magnetic catalysis.It also increases the absolute value of c B at which the Wilson loop survives.For ν = 4.5 the Wilson loop curve shrinks tending to horizontal µ-axis.For the heavy quarks model the larger anisotropy was the weaker it decreased the Wilson loop temperature.
• For ν = 1 in the light quarks model the confinement/deconfinement phase transition is determined by the crossover from µ = 0 to the intersection with the 1-st order phase transition µ by 2 , where it picks up the main role.On the opposite, for the heavy quarks model [9] the 1-st order phase transition was decisive from µ = 0 to the intersection with the crossover µ by 2 .
• Magnetic field c B = 0 for ν = 1 shifts the intersection point (µ by 2 , T by 2 ) to the larger chemical potentials and the lower temperatures, but doesn't actually change the general scheme for the light quarks model.For the heavy quarks model [9] the 1-st order phase transition moved to above the Wilson loop curve and lost it's influence on the confinement/deconfinement process.
• Anisotropy ν > 1 makes the 1-st order phase transition to move under the Wilson loop curve.Therefore the crossover determines the confinement/deconfinement process from µ = 0 to µ CEP , where the 1-st order phase transition begins and takes the main role with a jump.For the heavy quarks model [9] the 1-st order phase transition also moved under the Wilson loop curve and determined the confinement/deconfinement process from µ = 0 to µ CEP HQ , where the 1-st order phase transition ends and gives the main contribution to the crossover with a jump.
We have also studied the DC conductivity of the QGP in the deconfined phase under different conditions.The conductivity has different values and thermodynamic properties in the HIC direction and in the orthogonal plane presented, see Fig. 8-12.
• For ν = 1 the conductivity σ ii /T near the critical temperature grows with the temperature up to ≈ 2T c and then starts to decrease to a constant value near the conductivity of N = 4 SYM plasma along all the directions; along the orthogonal to the magnetic field directions is lower than in a parallel one; grows with the chemical potential and the external magnetic field values.
• For ν = 4.5 the conductivity σ 33 /T along the magnetic field direction has a minimum and increases at high temperatures (Fig. 9); σ 22 /T along the transversal direction also has a minimum followed by the increment at high temperatures (Fig. 10); σ 33 and σ 22 can be distinguished at the non-zero magnetic fields, and σ 33 is always greater than σ 22 (Fig. 11); σ 11 /T along the longitudinal direction monotonically decreases (Fig. 12).
It is instructive to compare the conductivity for the light and heavy quarks.Let us remind that for the heavy quark model we had [31] the following picture.
• For ν = 1 at temperatures higher than the critical value the ratios σ ii /T increase to the constant values depending on µ and c B ; increasing µ or magnetic field increases the values of these constant (saturation) values; the conductivity σ 22 and σ 33 are almost indistinguishable for small values of the external magnetic field.
• For ν = 4.5 the conductivity σ 22 /T and σ 33 /T along the transversal directions increase with the temperature growth and saturate a certain value σ 22 /T ≈ σ 33 /T depending on the µ and magnetic field values; σ 11 /T along the longitudinal direction approaches zero at high temperatures.
Thus, we see that for ν = 1 the conductivity behavior for the heavy and light quarks is different at near-critical temperatures, but at high temperatures they saturate some constant values and these values increase with the chemical potential and magnetic field growth.The saturation constants are equivalent at zero µ and c B .Isotropic QGP for the light quarks as well as for the heavy quarks is almost transparent in all directions near the critical temperature and become opaque at high temperature.The main difference between them may be seen at intermediate temperatures.
For ν = 4.5 the conductivity behavior for the heavy and light quarks is different at near-critical temperature as well as for high temperatures.In particular, there is no saturation of σ 33 and σ 22 at high temperature for the light quarks, meanwhile these components of the conductivity for the heavy quarks reach a constant value at high temperature.The longitudinal component of the conductivity for the heavy quarks decreases to zero as it does for the light quarks.
There are some obvious open questions regarding our model that are worth exploring.Let us mention some of them.It would be interesting to study the quark potential for the light quark model considered here like it was done for the heavy quark model [64].
It would also be interesting to study the entanglement entropy (EE) behavior for our model.Note that EE was studied in the [14] model and was not studied either in the presence of an external magnetic field for the heavy quark model [9], or for the light quark model in the absence of an external magnetic field [10].
In the spirit of the jet quenching studying [65,66] for the simplest [14] model, it would be interesting to investigate the behavior of the jet quenching for the model considered in this article.Of particular interest is the consideration of the external magnetic field influence on the results.
In this paper we determine the MC/IMC as the increase/decrease of the first order phase transition temperature with increasing the magnetic field value.In lattice calculations the MC and IMC are determined as increase/decrease of the chiral phase transition temperature with increasing the value of magnetic field [67][68][69][70].Note that the chiral phase transition coincides with the first order phase transition in the holographic model for heavy quarks [21], see also studies of the chiral condensate in holography in [71][72][73][74][75] and references therein.The study of the quark condensate and the chiral phase transition is certainly of interest for the structure of the phase transitions for light quarks.

Figure 7 .
Figure 7. Gauge kinetic function f 0 as a function of the horizon z h (A) and relation of electrical conductivity to temperature σ/T as a function of the normalised temperature T /T c (B): numerical calculations in our model for ν = 1, c B = 0 and µ = 0 (the blue curve), lattice calculations from [63] (the blue dots with error bars) and the conductivity of N = 4 SYM plasma (the red dash-dotted line).

Figure 11 .
Figure 11.The dependence of σ 22 /T (red curves) and σ 33 /T (blue curves) on the normalised temperature T /T c in anisotropic case ν = 4.5 for the magnetic field parameter c B = − 0.07 with zero chemical potential µ = 0 (A) and large chemical potential µ = 1 (B); critical temperature T c = 0.121 GeV.