Magnetic Catalysis in Holographic Model with Two Types of Anisotropy for Heavy Quarks

: In our previous paper [1] we have constructed a twice anisotropic five-dimensional holographic model supported by Einstein-dilaton-three-Maxwell action that reproduced some essential features of the “heavy quarks” model. However, that model did not describe the magnetic catalysis (MC) phenomena expected from lattice results for the QGP made up from heavy quarks. In this paper we fill this gap and construct the model that improves the previous one. It keeps typical properties of the heavy quarks phase diagram, and meanwhile possesses the MC. The deformation of previous model includes the modification of the “heavy quarks” warp factor and the coupling function for the Maxwell field providing the non-trivial chemical potential.


Introduction
Quantum chromodynamics (QCD) is a theory that describes strong interactions between subatomic particles such as quarks and gluons.Complete description of the QCD phase diagram in a parameter space with temperature, chemical potential, quark masses, anisotropy, magnetic field etc. is a challenging and very important task in high energy physics.Standard methods to do calculations in QCD such as perturbation no longer work for the strongly coupled regime of this theory, while the lattice theory has problems with non-zero chemical potential calculations.Hence, to understand physics of the strongly coupled quark-gluon plasma (QGP) produced in heavy ion collisions (HIC) at RHIC and at the LHC, and future experiments, we need a non-perturbative approach [2][3][4].
The enhancement effect of the phase transition temperature under the magnetic field increasing is known as MC phenomenon, the opposite effect is called IMC.Lattice calculations show that for small chemical potential there is a substantial influence of the magnetic field on the QCD phase diagram structure.This influence essentially depends on the quark mass: for small quark mass (light quarks) IMC takes place, meanwhile for large mass (heavy quarks) MC occurs.In this context note that lattice calculations predict different types of phase transitions even for small chemical potential and zero magnetic field -we have a crossover for light quarks, and a first-order phase transition for heavy quarks.The holographic QCD models for heavy and light quarks constructed in [37,[44][45][46][47][48] reproduce these phase diagram features at small chemical potential and predict new interesting phenomena for finite chemical potential, in particular, the locations of the critical end points.In our previous papers [49], see also [29], where the light quark holographic model with non-zero magnetic field is investigated, it has been shown that IMC takes place.Our paper [1] shows that the heavy quark holographic model [45] still has IMC, not MC, that contradicts with lattice zero chemical potential calculation.This indicates that one has to modify the heavy quark holographic model [1].
In the current paper we fill this gap and construct a heavy quark model that improves the previous one [1,45,50].The main goal of the improvements is to get the MC phenomenon in holographic description of the heavy quarks' first order phase transition scenario with external magnetic field keeping typical properties of the heavy quarks phase diagram.For this purpose we can consider additional z 4 - [51][52][53] or/and z 5 -terms [54,55] into the exponent warp factor.In particular, within this holographic model we show that z 4 -term allows to produce the MC phenomenon required.
As we have emphasized in the previous papers [45,46,56], there is a reason to introduce one more parameter characterizing the QCD phase diagram -an anisotropy parameter ν.Non-central HIC produces anisotropic QGP, and the isotropisation time is estimated as 1-5 fm/c ∼ 10 −24 s [57].Anisotropic holographic models have been used to study QGP in [56,[58][59][60][61][62][63][64][65][66][67].One of the main purposes to consider anisotropic models is to describe the experimental energy dependence of total multiplicity of particles created in HIC [68].In [56] it has been shown that the choice of the primary anisotropy parameter value about ν = 4.5 reproduces the energy dependence of total multiplicity [68].Note that isotropic models could not reproduce it (for more details see [56] and references therein).In addition, it is very interesting to know how the primary (spatial) anisotropy can affect the QCD phase transition temperature.Note also that there is another type of anisotropy due to magnetic field and its effect on the QCD phase diagram is a subject of interest.
In this work we set up a twice anisotropic "heavy quarks" model.In fact, we consider 5-dim Einstein-Maxwell-dilaton action with three Maxwell fields: the first Maxwell field sets up finite non-zero chemical potential in the gauge theory, the second Maxwell field provides the primary spatial anisotropy to reproduce the multiplicity dependence on energy, and the 3-rd Maxwell field provides another anisotropy that originates from magnetic field in the gauge theory.We use an anisotropic metric as an ansatz to solve Einstein equations and the field equations self-consistently.The central question of the current investigation is the form of the warp factor able to provide the MC phenomenon within the constructed holographic model.This our consideration shows a phenomenological character of the bottom-up holographic models [1, 4, 45-47, 49, 50, 53, 54, 69-98], that is different from the top-down holographic models [99][100][101][102][103].
This paper is organized as follows.In section 2 we present a 5-dim holographic model to describe a hot dense anisotropic QCD in the magnetic field background.In section 3 we introduce an appropriate warp factor able to produce MC phenomenon in this holographic model and obtain the first order phase transition for the model parameters.In section 4 we review our main results.This work in complemented with Appendix A where we solve EOMs, Appendix B where we present expressions for the blackening function derivatives, gauge coupling functions and dilaton potential, and Appendix C where we consider the relation of our setting with the setting [53] explicitly.

Holographic Model with three Maxwell Fields
Let us take the Lagrangian in Einstein frame used in [1]: where R is Ricci scalar, ϕ is the scalar field, f 0 (ϕ), f 1 (ϕ) and f 3 (ϕ) are the coupling functions associated with stresses F 0 , F 1 and F 3 of Maxwell fields, and V (ϕ) is the scalar field potential.In this paper we considered F 0 , F 1 and F 3 as first, second and third Maxwell fields, respectively.Varying Lagrangian (2.1) over the metric we get Einstein equations of motion (EOMs): where and varying over the fields gives the fields equations Let us take the metric ansatz in the following form: and for matter fields (2.8) In (2.6) L is the AdS-radius, b(z) is the warp factor set by A(z), g(z) is the blackening function, ν is the parameter of primary anisotropy caused by non-symmetry of heavyion collision (HIC), and c B is the coefficient of secondary anisotropy related to the magnetic field F 3 .Choice of A(z) determines the heavy/light quarks description of the model.In previous works we considered A(z) = − cz2 /4 for heavy quarks [1,45,50] and A(z) = − a ln(bz 2 + 1) for light quarks [46,49].In (2.8) q 1 and q 3 are constant "charges".The explicit form of the EOM (2.3-2.5) with ansatz (2.7)-(2.8) is given in Appendix (A.12-A.18).Investigation of their self-consistency shows that there is one dependent equation in the system and all other equations are independent.Thus, system (A.12-A.18) is self-consistent and the dilaton field equation (A.12) serves as a constraint.
It is important to note that the coupling function f 0 is defined from the requirement to reproduce the Regge trajectories.Functions f 1 and f 3 are obtained from the EOM, and we find that they are different (see Appendices A and B for more details).If we take f 0 = f 1 = f 3 , then we cannot construct a solution of EOM within our ansatz.
Note that while solving equations of motions we do not actually get f 0 (ϕ), f 1 (ϕ), f 3 (ϕ) and V (ϕ) dependencies, but obtain The reason is that ϕ(z)-expression is rather complicated, so the analytic expression for the inverse function z = z(ϕ) can't be written down.There still remains the possibility to get this function via approximation, but such a result can be useful for a limited number of aspects only because of lack of accuracy.

Magnetic Catalysis for Heavy Quarks
Our goal is to generalize solution [1] to get magnetic catalysis effect on the heavy quarks version of the phase diagram.For this purpose we choose the deformation of the warp factor from [1] again.This factor has been used in [51,52] to reproduce Cornell potential and in [53] also to get magnetic catalysis effect.

Solution and thermodynamics for
Our strategy to solve the EOMs presented in Appendix A with the factor A(z) = − cz 2 /4 − pz 4 is the same as in [1] and [46].Subtracting (A.17) from (A.16) we get the expression for the third Maxwell field's coupling function and rewrite equation (A.14) as To derive the exact solutions we just need to specify the warp factor: Following [48,53,95] we take c = 4R gg /3, R gg = 1.16, p = 0.273 (this choice is dictated by the Regge spectra and lattice QCD fitting) and solve system (A.12)-(A.18)with usual boundary conditions where z 0 serves to fit the string tension behavior [46].Equation (A.13) with (C.8) gives For q 3 = 1 and c B = C the result (3.7) coincides with the expressions (2.27) and (2.31) in [53]: Density is the coefficient in A t expansion: lg(ρ/μ) and density ⇢(z h )/µ in logarithmic scale (B) for the "heavy quarks" case with (solid lines) and without (dashed lines) magnetic field, We use the following formulas of temperature and entropy: For the metric (2.6) and the warp factor (3.3) temperature and entropy can be written as: In Fig. 2 (first line) we see that magnetic "charge" q 3 a↵ects the temperature function for the fixed magnetic coe cient c B and non-zero chemical potential µ values independently from primary anisotropy.For zero chemical potential temperature obviously has no dependence on q 3 at all, see eq. (3.12).At fixed c B increasing q 3 decreases -7 - A B Figure 1.Electric potential A t (z) (A) and density ρ(z h )/µ in logarithmic scale (B) for the "heavy quarks" case with (solid lines) and without (dashed lines) magnetic field, The electric potential A t (z) and density ρ(z h )/µ in logarithmic scale are depicted in Fig. 1 in panels A and B, respectively.Equation (3.2) with (C.8) and (3.7) gives We use the following formulas of temperature and entropy: ) For the metric (2.6) and the warp factor (3.3) temperature and entropy can be written as: T q 3 = 0.5 T q 3 = 0.5 the phase transition temperature and evenually makes it monotonic (Fig. 2, second line) independently from primary anisotropy as well.All this leads to the idea used in other works to associate the magnetic "charge" with the magnetic coe cient q 3 = c B .
To understand the role of pz 4 -term in the warp factor let's consider the temperature behavior for p = 0, 0.1, 0.2, 0.273 (Fig. 3).Earlier in [1] we obtained a hypersensitivity  of the background phase transition from the magnetic field: for c B ⇠ 0.01 (for ⌫ = 1) it fully degenerated.Non-zero p makes the background phase transition in the magnetic field more stable, allowing it to achieve more realistic values (in Fig. 3 q 3 = c B = 0.2).In Fig. 4 temperature function for di↵erent chemical potential values and p = 0.273 in presented.Magnetic field amplification makes three-valued behavior of z h (T ) to become monotonic and lowers the local minimum temperature.Entropy has usual behavior (Fig. 5) and the inverse magnetic catalysis can be expected again.To make this sure let us consider a free energy function of the background phase transition from the magnetic field: for c B ⇠ 0.01 (for ⌫ = 1) it fully degenerated.Non-zero p makes the background phase transition in the magnetic field more stable, allowing it to achieve more realistic values (in Fig. 3 q 3 = c B = 0.2).In Fig. 4 temperature function for di↵erent chemical potential values and p = 0.273 in presented.Magnetic field amplification makes three-valued behavior of z h (T ) to become monotonic and lowers the local minimum temperature.Entropy has usual behavior (Fig. 5) and the inverse magnetic catalysis can be expected again.To make this sure let us consider a free energy function In Fig. 2 (first line) we see that magnetic "charge" q 3 affects the temperature function for the fixed magnetic coefficient c B and non-zero chemical potential µ values independently from primary anisotropy.For zero chemical potential temperature obviously has no dependence on q 3 at all, see eq. (3.12).At fixed c B increasing q 3 decreases the phase transition temperature and evenually makes it monotonic (Fig. 2, second line) independently from primary anisotropy as well.All this leads to the idea used in other works to associate the magnetic "charge" with the magnetic coefficient To understand the role of pz 4 -term in the warp factor let's consider the temperature behavior for p = 0, 0.1, 0.2, 0.273 (Fig. 3).Earlier in [1] we obtained a hypersensitivity of the background phase transition from the magnetic field: for c B ∼ − 0.01 (for ν = 1) it fully degenerated.Non-zero p makes the background phase transition in the magnetic field more stable, allowing it to achieve more realistic values (in Fig. 3 q In sition or Hawking-Page-like (HP) phase transition originating from the free energy is presented in Fig. 6 and confirms the inverse magnetic catalysis e↵ect.As to the pz 4term e↵ect, in Fig. 7 we can clearly see that larger p-value prevents the background phase transition degeneracy with the magnetic field amplification.Thus our choice of (3.3) as a warp factor to deform the metric did not possess the in presented.Magnetic field amplification makes three-valued behavior of z h (T ) to become monotonic and lowers the local minimum temperature.Entropy has usual behavior (Fig. 5) and the inverse magnetic catalysis can be expected again.To make this sure let us consider a free energy function MC phenomenon although could lead to good results for IMC.In the next section we consider a new warp factor to produce the proper result for MC phenomenon.

Solution for
In search for the MC phenomenon we considered a corrected version of our new warp factor.Let us remind that varying Lagrangian (2.1) over the metric (2.6) we get the same equations of motion (A.12-A.18)for di↵erent fields.Using the new warp factor MC phenomenon although could lead to good results for IMC.In the next section we consider a new warp factor to produce the proper result for MC phenomenon.

Solution for
In search for the MC phenomenon we considered a corrected version of our new warp factor.Let us remind that varying Lagrangian (2.1) over the metric (2.6) we get the same equations of motion (A.12-A.18)for di↵erent fields.Using the new warp factor Note that for non-zero chemical potential we should integrate to second horizon z h 2 and for µ = 0 we have The result for the first order (background) black hole-black hole (BB) phase transition or Hawking-Page-like (HP) phase transition originating from the free energy is presented in Fig. 6 and confirms the inverse magnetic catalysis effect.As to the pz 4term effect, in Fig. 7 we can clearly see that larger p-value prevents the background phase transition degeneracy with the magnetic field amplification.
Thus our choice of (3.3) as a warp factor to deform the metric did not possess the MC phenomenon although could lead to good results for IMC.In the next section we consider a new warp factor to produce the proper result for MC phenomenon.

Solution for
In search for the MC phenomenon we considered a corrected version of our new warp factor.Let us remind that varying Lagrangian (2.1) over the metric (2.6) we get the same equations of motion (A.12-A.18)for different fields.Using the new warp factor where c = 4R gg /3, R gg = 1.16, p = 0.273, one can solve system of EOMs (A.12-A.18)with the same boundary conditions (3.4-3.6).
To possess the linear Regge trajectories for the meson mass spectra in our model, in comparison with [53], we considered the kinetic function (C.8) and the new warp factor (3.16).Then, at T = µ = B = 0 one can produce linear mass spectrum m 2 n = 4R gg n, in such a way that the parameter R gg can be fitted by Regge spectra of meson, such as J/ψ.Note also that the parameters R gg and p can be fixed for the zero magnetic field with R gg = 1.16 and p = 0.273 [48,95].

Blackening function
For the corrected factor (3.16) the EOM on the gauge field A t (z) is the same as before and has the same solution (3.7).Therefore, the equation (3.2) with (C.8), (3.7) and the corrected warp factor (3.16) gives The behavior of the blackening function g in terms of the holographic coordinate z for different values of the magnetic coefficient c B and different primary anisotropy background, i.e. ν = 1 (green lines) and ν = 4.5 (blue lines), normalized to the horizon z h = 1, is depicted in Fig. 8.A.Here blackening function is monotonic, and larger values of magnetic coefficient c B correspond to lower g values both in isotropic and anisotropic background.But comparing isotropic and anisotropic backgrounds we see that at 0 < |c B | < 0.3 the blackening function gets lower values for ν = 4.5 than for ν = 1, while at 0.3 < |c B | < 1 it is not sensible to changes in primary anisotropies.ϕ The dilaton field (z) for di↵erent c B , q 3 = 5 (A) and di↵erent q 3 , c B = 0.5 (B); ⌫ = 1 (green lines) and ⌫ = 4.5 (blue lines); R gg = 1.16, p = 0.273, z 0 = 0.001;

Scalar field
The scalar, i.e. dilaton field can be obtained from equation (A.15)2 with the boundary condition (3.6) imposed 2 This is a generalization of the method [104,105] of reconstructing the dilaton potential, see also [106][107][108].
-14 -A B The effect of chemical potential on the blackening function for different primary anisotropies of the background is demonstrated in Fig. 8.B.Larger µ decreases the blackening function value in both isotropic and anisotropic background cases.But, for the fixed chemical potential the blackening function value is smaller in the background with larger primary anisotropy.Note also that for large chemical potential values one has to deal with the second black hole horizon.

Scalar field
The scalar, i.e. dilaton field ϕ can be derived from equation (A.15) 2 with the boundary condition (3.6) imposed νξ . (3.20) Expanding the integrand of the dilaton field we have ϕ(z) Therefore, the dilaton field has no divergency at z 0 = 0 on the primary isotropic background ν = 1, while on anisotropic background a logarithmic divergency exists.It is important to note that we generalize the boundary condition for the dilaton field as ϕ(z 0 ) = 0 [46,96], where z 0 can be some function of z h .The fact is that the scalar field boundary conditions can affect the temperature dependence of the string tension, i.e. the coefficient of the linear term of the Cornell potential.The string tension should decrease with the temperature growth and become zero at the confinement/deconfinement phase transition [109][110][111].To preserve this feature and also avoid divergences in anisotropic backgrounds we considered the dilaton boundary condition as ϕ(z 0 ) = 0.For special cases one can consider z 0 = 0 [48] and z 0 = z h [45].
Fig. 9 shows that the scalar field is a monotonically increasing function of the holographic coordinate z both in primary isotropic and anisotropic cases, i.e. for ν = 1 and ν = 4.5, but larger primary anisotropy shifts the dilaton curve up to larger ϕ-values.Larger absolute value of the magnetic coefficient c B and larger magnetic charge q 3 make ϕ(z) to grow faster (Fig. 9.A and B respectively).

Coupling function f 3
The function f 3 that describes the coupling between the third Maxwell field F 3 and the dilaton field ϕ is still calculated by the expression (3.1).The detailed formula obtained by substituting the blackening function, it's derivative and the corrected warp factor (3.16) is presented in Appendix B. For zero magnetic coefficient c B = 0 the coupling function f 3 obviously equals zero.It's behavior depending on the holographic direction z is plotted in Fig. 10.We see that for zero chemical potential µ = 0 it is positive and preserves the NEC both in isotropic ν = 1 and primary anisotropic ν = 4.5   According to the boundary condition g(z h ) = 0 expression (3.1) is simplified as At the first horizon (the one with smaller z h value, that really matters) the blackening function is decreasing, so g 0 (z h ) < 0 and temperature T (z h ) = |g 0 (z h )/4⇡| is positive.
If we also take c B < 0, their product is positive, all the other multipliers in   According to the boundary condition g(z h ) = 0 expression (3.1) is simplified as At the first horizon (the one with smaller z h value, that really matters) the blackening function is decreasing, so g 0 (z h ) < 0 and temperature T (z h ) = |g 0 (z h )/4⇡| is positive.
If we also take c B < 0, their product is positive, all the other multipliers in backgrounds (Fig. 10.A).However in the isotropic background we see the decreasing monotonic behavior of f 3 (z), while in the anisotropic background it demonstrates a multivalued behavior with a local minimum and a local maximum.Note also that for larger magnetic coefficient (larger c B absolute values) f 3 is positive not everywhere beyond the fixed horizon z h , therefore we need to choose appropriate parameters in our theory and in particular the correct value for the second horizon to have positive value for coupling function f 3 (see the Fig. 10.B for the second horizon z h 2 = 0.55 for ν = 1 and z h 2 = 0.45 for ν = 4.5).
According to the boundary condition g(z h ) = 0 expression (3.1) is simplified as At the first horizon (the one with smaller z h value, that really matters) the blackening function is decreasing, so g ′ (z h ) < 0 and temperature T (z h ) = |g ′ (z h )/4π| is positive.
If we also take c B < 0, their product is positive, all the other multipliers in (3.21) are positive as well, therefore f 3 (z h ) > 0 for any negative c B in the z interval we need for 0 < z < z h min , where z h min is not the fixed horizon, but the horizon with smaller value, i.e.T (z h min ) = 0.For zero chemical potential z h min = z h and for µ = 1, for example, z h min = z h 2 (Fig. 8.B).In Fig. 11 the third coupling function f 3 in terms of dilaton field ϕ is displayed.It demonstrates a nonmonotonic behavior, quite sensible to the magnetic field presence -larger c B absolute value leads to larger f 3 .Neither fixed horizon (Fig. 11.A), nor chemical potential (Fig. 11.B) have no significant effect.

Coupling function f 1
We also need to check the NEC for the function f 1 in our model.It describes coupling between the second Maxwell field F 1 and the dilaton field ϕ: The exact formula obtained by substituting the blackening function, its derivative and the new warp factor is presented in Appendix B. It is clearly seen that the coupling function f 1 is zero for ν = 1, as according to the holographic dictionary the second Maxwell field F 1 serves for primary anisotropy of the background.Fig. 12 shows the coupling function f 1 in terms of holographic coordinate z.It gets positive value for zero chemical potential (Fig. 12.A, light blue), so the NEC is fulfilled in our model.For perceptible chemical potential value µ = 0.5 (Fig. 12.A, dark blue) the coupling function f 1 can be not positive for z < z h in magnetic field strong enough.But the second horizon shift under the first (fixed) one must be taken in account again.For z < z h 2 the coupling function f 1 stays positive similar to the coupling function f 3 , thus NEC is respected (Fig. 12.B, the second horizon z h min = z h 2 = 0.39).For f 1 it can be proven algebraically, like it was done for it in [1] and for f 3 in previous section, that it doesn't break the NEC either.According to the boundary condition g(z h ) = 0 expression for f 1 can be simplified as At the first horizon z h = z h min the blackening function is decreasing, so g ′ (z h ) < 0 and temperature T (z h ) = |g ′ (z h )/4π| is positive.As f 1 has sense for ν > 1 only, all the multipliers in (3.23) are positive, so their product is positive, therefore f 1 (z h ) > 0 for any z h from the interval 0 < z h < z h min , where z h min is the actual horizon, i.e.T (z h min ) = 0.
For the completeness of the results presented the second coupling function f 1 (ϕ) is also plotted (Fig. 13).It turns out to be insensitive to the magnetic field for small ϕ values and almost insensitive to the fixed horizon value even for large ϕ (Fig. 13.A).Besides, the magnetic field shifts f 1 up to larger values, while the chemical potential, on the contrary, downs to smaller values (Fig. 13.B).
In Fig. 11 the third coupling function f 3 in terms of dilaton field is displayed.It demonstrates a nonmonotonic behavior, quite sensible to the magnetic field presence -larger c B absolute value leads to larger f 3 .Neither fixed horizon (Fig. 11.A), nor chemical potential (Fig. 11.B) have no significant e↵ect.

Coupling function f 1
We also need to check the NEC for the function f 1 in our model.It describes coupling between the second Maxwell field F 1 and the dilaton field : The exact formula obtained by substituting the blackening function, its derivative and the new warp factor is presented in Appendix B. It is clearly seen that the coupling function f 1 is zero for ⌫ = 1, as according to the holographic dictionary the second Maxwell field F 1 serves for primary anisotropy of the background.Fig. 12 shows the coupling function f 1 in terms of holographic coordinate z.It gets positive value for zero chemical potential (Fig. 12.A, light blue), so the NEC is fulfilled in our model.For perceptible chemical potential value µ = 0.5 (Fig. 12.A, dark blue)   the coupling function f 1 can be not positive for z < z h in magnetic field strong enough.But the second horizon shift under the first (fixed) one must be taken in account again.For z < z h 2 the coupling function f 1 stays positive similar to the coupling function f 3 , thus NEC is respected (Fig. 12.B, the second horizon z h min = z h 2 = 0.39).For f 1 it can be proven algebraically, like it was done for it in [1] and for f 3 in previous section, that it doesn't break the NEC either.According to the boundary condition g(z h ) = 0 expression for f 1 can be simplified as At the first horizon z h = z h min the blackening function is decreasing, so g 0 (z h ) < 0 and temperature T (z h ) = |g 0 (z h )/4⇡| is positive.As f 1 has sense for ⌫ > 1 only, all the multipliers in (3.23) are positive, so their product is positive, therefore f 1 (z h ) > 0 for any z h from the interval 0 < z h < z h min , where z h min is the actual horizon, i.e.T (z h min ) = 0.For the completeness of the results presented the second coupling function f 1 ( ) is also plotted (Fig. 13).It turns out to be insensitive to the magnetic field for small values and almost insensitive to the fixed horizon value even for large (Fig. 13.A).Besides, the magnetic field shifts f 1 up to larger values, while the chemical potential, on the contrary, downs to smaller values (Fig. 13.B).

Scalar potential
The equation of motion for the scalar potential V as a function of holographic coordinate z can be obtained from the equation (A.18): (3.24) The exact formula obtained by substituting the blackening function, its first and second derivatives and the new corrected warp factor is presented in Appendix B. The scalar potential V as a function of dilaton field is described by equations (3.20) and (3.24) and is depicted in Fig. 14.Within the interval 0  z  z h , that we are interested in, V ( ) is stable and has no crucial dependence either on the horizon nor on the magnetic field.More focus in the scale allows to see that the magnetic field increases the absolute value of the scalar potential, while z h growth decreases it.Primary anisotropy (Fig. V

Scalar potential
The equation of motion for the scalar potential V as a function of holographic coordinate z can be obtained from the equation (A.18): The exact formula obtained by substituting the blackening function, its first and second derivatives and the new corrected warp factor is presented in Appendix B. The scalar potential V as a function of dilaton field ϕ is described by equations (3.20) and (3.24) and is depicted in Fig. 14.Within the interval 0 ≤ z ≤ z h , that we are interested in, V (ϕ) is stable and has no crucial dependence either on the horizon nor on the magnetic field.More focus in the scale allows to see that the magnetic field increases the absolute value of the scalar potential, while z h growth decreases it.Primary anisotropy (Fig. 14.B for ν = 4.5) shifts the scalar potential to larger absolute values (Fig. 14.A for ν = 1), causes a constant region to appear at small ϕ and then makes V (ϕ) to decrease faster.

Temperature and entropy
Using the metric (2.6) and the warp factor b(z) = e − cz 2 /2 − 2(p−c B q 3 )z 4 (3.16) one can obtain the temperature and entropy from equations (3.12) and (3.13) respectively: The behavior of the temperature T as a function of the horizon radius z h for different values of the magnetic coefficient c B (1-st line) and the magnetic "charge" q 3 (2-nd line) in backgrounds with different primary anisotropy ν = 1 (A) and ν = 4.5 (B) and zero chemical potential µ = 0 is shown in Fig. 15.The system temperature is obviously sensible to the magnetic field parameters in such a way that the temperature to both magnetic field parameters, and the Hawking-Page phase transition temperature T HP increases.Primary anisotropy makes the temperature local minimum to decrease again and therefore the Hawking-Page phase transition T HP requires less energy.To investigate the process precisely we need to plot a free energy for di↵erent chemical potential values.T q 3 = 0 q 3 = 1 T q 3 = 0 q 3 = 1  T

Free energy and magnetic catalysis
To investigate the first order phase transition we start from the free energy F as a function of temperature T consideration for di↵erent values of the magnetic field parameters c B and q 3 at zero chemical potential.Regardless of the primary background anisotropy (⌫ = 1 in Fig. 17 minimum grows with the magnetic coefficient c B absolute value (Fig. 15, 1-st line) and the magnetic "charge" q 3 (Fig. 15, 2-nd line).Naively, it means that the Hawking-Page phase transition temperature T HP value increases with the magnetic field in all the cases considered.Note that to investigate the magnetic field effect on the phase transition temperature below, we consider the magnetic coefficient c B , although the magnetic F q 3 = 0.75 q 3 = 1.5  "charge" q 3 acts on the transition temperature in a similar way.As to the primary anisotropy, it makes the temperature local minimum to decrease both for various c B and q 3 .Therefore, one can expect the Hawking-Page phase transition temperature T HP to increase during the isotropisation process.But to check this point and the BB phase transition explicitly, we need to calculate the free energy of the system.Now, we consider the non-zero chemical potential as we need to include matter to investigate the realistic QGP with high baryonic density.In Fig. 16 the system temperature T in terms of the horizon z h at the chemical potential µ = 0.3 for different values of c B (1-st line) and q 3 (2-nd line) is plotted.The temperature is still very sensible to both magnetic field parameters, and the Hawking-Page phase transition temperature T HP increases.Primary anisotropy makes the temperature local minimum to decrease again and therefore the Hawking-Page phase transition T HP requires less energy.To investigate the process precisely we need to plot a free energy for different chemical potential values.

Free energy and magnetic catalysis
To investigate the first order phase transition we start from the free energy F as a function of temperature T consideration for different values of the magnetic field parameters c B and q 3 at zero chemical potential.Regardless of the primary background anisotropy (ν = 1 in Fig. 17.A and ν = 4.5 in Fig. 17.B) the free energy is a multivalued function of temperature and at µ = 0 has two branches.One of them is positive but asymptotically decreasing to zero, the other goes down into the region of negative values sharply, almost vertically.The negative part of this branch describes the large stable black hole, positive free energy values correspond to the small unstable black hole and at F = 0 thermal gas is formed.Temperature of the phase transition from a small black hole to thermal gas known as Hawking-Page phase transition is so called Hawking-Page temperature, F (T HP ) = 0.According to the holographic dictionary this process corresponds to the first order phase transition in the dual 4-dim gauge theory.To reveal the magnetic field effect on the free energy behavior and, consequently, on the phase transition temperature T HP , we increase the absolute values of the magnetic parameters c B (Fig. 17, 1-st line) and q 3 (Fig. 17, 2-nd line).In both cases T HP grows, thus the magnetic catalysis phenomenon takes place for any primary anisotropy, but in the background with higher primary anisotropy phase transition requires lower temperature, like it was in previous works [1,[37][38][39][40]49].
At a non-zero chemical potential there appears the free energy third (almost horizontal) branch, and a closed structure called a "swallow-tail" appears (Fig. 18, µ = 0.3).Its self-intersection point describes the first order phase transition, and its temperature is marked as BB temperature T BB .All other tendencies preserve: magnetic field parameters c B and q 3 cause the phase transition growth while the primary anisotropy ν lowers it.This allows us to expect the MC-effect to be general on the first order phase diagram.All these results are summarized in Fig. 19, where the phase transition temperature dependence on the magnetic coefficient c B for µ = 0 (A) and µ = 0.3 (B) in the background with different primary anisotropy is depicted; the results for various chemical potentials in primary isotropic case are confronted in Fig. 19.C.Hawking-Page temperature, F (T HP ) = 0.According to the holographic dictionary this process corresponds to the first order phase transition in the dual 4-dim gauge theory.
To reveal the magnetic field e↵ect on the free energy behavior and, consequently, on the phase transition temperature T HP , we increase the absolute values of the magnetic parameters c B (Fig. 17, 1-st line) and q 3 (Fig. 17, 2-nd line).In both cases T HP grows, thus the magnetic catalysis phenomenon takes place for any primary anisotropy, - To obtain the full picture of the first order phase transition the free energy and consequently the temperature dependence on the chemical potential should be considered.At µ = 0 the temperature has minimum T min = T (z min ) (Fig. 20.A).Black hole solutions with z h > z min (small black holes) are unstable, and the essence of the phase transition lies in the collapse of such a black hole into a stable black hole with the same temperature, but z h < z min (large black holes).There is no black hole solution for T < T min , therefore Hawking-Page phase transition occurs at T HP ≥ T min .At the very point T min a large stable black hole transforms into the thermal gas.This point is clearly seen on the free energy plot as F (T min ) = 0 (Fig. 20

.B).
For µ > 0 the function T (z h ) also has a local maximum T max = T (z max ), where z max > z min (Fig. 20.A), and this leads to the "swallow-tail" appearance in the free energy plot again (Fig. 20.B).Temperature of the phase transition, i.e. the collapse from a small unstable black hole to a large stable black hole, can be determined as the temperature of the self-intersection point in the "swallow-tail" base of the F (T )curve.As the chemical potential increases both T (z min ) and T (z max ) decrease so that the difference between them reduces and eventually disappears at some critical value of the chemical potential.For example, in the magnetic field with c B = − 0.5 for ν = 1 and ν = 4.5 it happens at µ CEP HQ ≈ 0.94 and µ CEP HQ ≈ 0.81 correspondingly (Fig. 20.A).This process is reflected by the "swallow-tail" decrease on the free energy plot (Fig. 20.B).For µ > µ CEP HQ the black hole temperature becomes a monotonic function of the horizon and its free energy becomes smooth.Note that additional details on this subject can be found in papers on previous considerations [1,45,46].

Phase diagrams
Phase diagram, i.e. the diagram of the confinement/deconfinement phase transition displays the dependence of the phase transition temperature on the chemical potential.This phase diagram consists of two different phase transition lines, i.e. first order phase transition line and phase transition line for temporal Wilson loops on the (µ,T )plane.In this research we considered the first order phase transition only.In the "heavy quarks" version of our holographic model the phase transition line stretches over the interval µ ∈ [0, µ max ], and the temperature in this interval drops (Fig. 21).The rightmost point of the curve with coordinates µ max , T (µ max ) is called a critical end but in the background with higher primary anisotropy phase transition requires lower temperature, like it was in previous works [1,[37][38][39][40]49].
At a non-zero chemical potential there appears the free energy third (almost horizontal) branch, and a closed structure called a "swallow-tail" appears (Fig. 18, µ = 0.3).Its self-intersection point describes the transition from the confinement to the deconfinement phase, and it's temperature is marked as BB temperature T BB .All other tendencies preserve: magnetic field parameters c B and q 3 cause the phase transition growth while the primary anisotropy ⌫ lowers it.This allows us to expect the MC-e↵ect to be general on the first order phase diagram.All these results are summarized in Fig. 19, where the phase transition temperature dependence on the magnetic coe cient c B for µ = 0 (A) and µ = 0.3 (B) in the background with di↵erent primary anisotropy is depicted; the results for various chemical potentials in primary isotropic case are confronted in Fig. 19.C.
To obtain the full picture of the first order phase transition the free energy and consequently the temperature dependence on the chemical potential should be considered.At µ = 0 the temperature has minimum T min = T (z min ) (Fig. 20.A).Black hole solutions with z h > z min (small black holes) are unstable, and the essence of the phase transition lies in the collapse of such a black hole into a stable black hole with the same temperature, but z h < z min (large black holes).There is no black hole solution for T < T min , therefore Hawking-Page phase transition occurs at T HP T min .At the very point T min a large stable black hole transforms into the thermal gas.This point is clearly seen on the free energy plot as F (T min ) = 0 (Fig. 20.B).
For µ > 0 the function T (z h ) also has a local maximum T max = T (z max ), where z max > z min (Fig. 20.A), and this leads to the "swallow-tail" appearance in the free energy plot again (Fig. 20.B).Temperature of the phase transition, i.e. the collapse from a small unstable black hole to a large stable black hole, can be determined as point (CEP HQ for heavy quarks) and marks the free (not attached to the axis) end of the phase transition line (dots on the right end of the T (µ)-curves).
First of all, we really see the MC effect in the background with any amount of anisotropy.For ν = 1 CEP chemical potential µ CEP HQ grows by increasing |c B | for the region 0 < |c B | < 0.5, after that µ CEP HQ decreases by increasing |c B | and gets the temperature of the self-intersection point in the "swallow-tail" base of the F (T )curve.As the chemical potential increases both T (z min ) and T (z max ) decrease so that the di↵erence between them reduces and eventually disappears at some critical value of the chemical potential.For example, in the magnetic field with c B = 0.5 for ).We presented critical end points for different c B with ν = 1 and ν = 4.5 in Table 1.
In both cases the T HP/BB rise slows down with increasing magnetic field, before the complete disappearance of phase transition lines it stops and even turns back.We observed that for the small region the IMC appears at near-limit c B values, namely |c B | > 5 for ν = 1 and |c B | > 8 for ν = 4.5 (light green and blue curves, respectively).But it is obvious that in the primary isotropic background the phase transition degenerates at a magnetic field weaker than in the primary anisotropic one.Besides for ν = 4.5 phase transition temperature is lower and drops faster as µ increases (Fig. 21.C).We

NOT FOR DISTRIBUTION JHEP_144P_0723 v1
⌫ = 1 and ⌫ = 4.5 it happens at µ CEP HQ ⇡ 0.94 and µ CEP HQ ⇡ 0.81 correspondingly (Fig. 20.A).This process is reflected by the "swallow-tail" decrease on the free energy plot (Fig. 20.B).For µ > µ CEP HQ the black hole temperature becomes a monotonic function of the horizon and its free energy becomes smooth.Note that additional details on this subject can be found in papers on previous considerations [1,45,46].-27 -

A B
⌫ = 1 and ⌫ = 4.5 it happens at µ CEP HQ ⇡ 0.94 and µ CEP HQ ⇡ 0.81 correspondingly (Fig. 20.A).This process is reflected by the "swallow-tail" decrease on the free energy plot (Fig. 20.B).For µ > µ CEP HQ the black hole temperature becomes a monotonic function of the horizon and its free energy becomes smooth.Note that additional details on this subject can be found in papers on previous considerations [1,45,46].

Conclusion and Discussion
In this research we studied the influence of the magnetic field on the first order phase transition temperature.For this purpose we used the "bottom-up" approach and chose 5-dim Einstein-dilaton-Maxwell holographic model with three Maxwell fields.In previous paper [1] only the IMC phenomenon was obtained.In this research we look for a warp factor, that serves for a deformation of metric, providing the MC phenomenon for the heavy quark model [1].As our holographic model is phenomenological, there is no systematic way to construct it.We have to use the trial and error method to get results compatible with existing experiments and others theoretical methods.
Let us summarize our main results.
• A new 5-dim exact analytical solution for anisotropic holographic model of QGP is presented.One of its important features is inclusion of two types of anisotropy, caused by the spatial anisotropy, called primary, and the external magnetic field.Their influences on the background physical properties such as background phase transition are investigated.
• Choosing the warp factor that deforms the metric is a key point of the current research.The warp factor b(z) = e 2A(z) = e − cz 2 /2 − 2pz 4 leads to IMC in the sense of decreasing of the critical temperature with increasing the magnetic field B.
• The MC effect is achieved for the warp factor b(z) = e − cz 2 /2 − 2(p−c B q 3 )z 4 for c B not to be large (see discussion of the Figure 21).This takes place for zero chemical potential, i.e. for Hawking-Page-like (HP) phase transition, and for non-zero chemical potential, i.e. black hole-black hole first order.In both cases we get MC in the dual 4-dim gauge theory.
• The effect of primary anisotropy on black hole-black hole (BB) and Hawking-Page-like (HP) phase transition is investigated.It is found that anisotropy decreases both T HP and T BB for all values of magnetic field.
• The phase diagram, i.e. the dependence of the phase transition temperature on chemical potential, is built for different magnetic field magnitudes and different primary anisotropy values within the model constructed.
• Complete disappearance of the phase transition lines for the primary isotropic background, ν = 1, occurs at weaker magnetic field (|c B | ∼ 6) in comparison to the anisotropic one, ν = 4.5 (|c B | ∼ 10).
• Even for the near-limit chemical potential values the NEC to be preserved, as the consideration has sense and can be performed within the physical interval between the boundary and the second horizon in this model.
• It is expected that for the presented model different quantities such as baryon density, entanglement entropy, electrical conductivity should have jump in the vicinity of the first order phase transition.This jump should strongly depend on the model parameters (anisotropy, magnetic field, chemical potential etc.) similar to previously considered models [49,83,97].
• In [1,45,50,112] the temporal and spatial Wilson loops were considered on the background for heavy quarks with two types of anisotropy with the warp factor b(z) = e − cz 2 /2 .It would be interesting to study phase transition on this background with the new corrected warp factor.Also one can investigate energy loss and jet quenching on this background similar to [113].This will allow to obtain full confinement/deconfinement phase transition structure, that is determined by the Wilson loop and first order phase transition interplay.
In this paper we do not make calculations of the chiral condensate ⟨ ψψ⟩.We would like to emphasize that to calculate ⟨ ψψ⟩ one has to consider a new action, i.e. chiral action including a few new fields χ, and solve corresponding equations of motion.We do not perform these calculations in the present paper, this will be the subject of the future investigations.Similar calculations have been performed in [54,114,115] in different models.The chiral condensate ⟨ ψψ⟩ has been calculated for light quark models in [114,115] for zero magnetic field.The chiral condensate ⟨ ψψ⟩ has been calculated in [54] for non zero B for the heavy quark holographic model with IMC in the sense of decreasing of the critical temperature with increasing B (i.e.without extra term z 4 in the warp factor).It has been found that ⟨ ψψ⟩ increasing with increasing B that could be called MC in the sense of the condensate value.and inserting the equation of g(z) (3.17After some algebra, one can obtain + µ 2 (c B (−1 + q 3 ) + 2R gg ) Ĩ2 (z h ) Ĩ1 (z h ) −e 1 2 z 2 (4Rgg+6(p−cB q 3 )z 2 −4c B +c B q 3) µ 2 (c B (−1 + q 3 ) + 2R gg ) × e C Comparison with [53] In this section we intend to compare the geometry of [1] with the metric and Lagrangian introduced in [53]: electro-magnetic ansatz F : A µ = A t (z), 0, 0, A 3 (x 2 ), 0 , µ = 0, 1, 2, 3, 4. (C.4) Components F 00 ∼ A t (z) and F 23 ∼ A 3 (x 2 ) of the united electro-magnetic field F µν from [53] formally correspond to electric Maxwell field F 0 µν and first magnetic Maxwell field F 1 µν from [1].Magnetic component F 23 acts along the x 1 -direction in (C.1), but it is real magnetic field, not an effective source of primary anisotropy ν, as [53] describes an isotropic model with magnetic field.Magnetic field (analogous to F µν from [1]) has non-zero component F 23 , acts along the x 1 -direction and influences g 11 making it differ from g 22 and g 33 .Therefore we can say that for models [53] and [1]  We also see that In [53] the coupling function f 0 and the warp factor are: where d(z) = − (R gg /3) z 2 −pz 4 , R gg = 1.16, p = 0.273, c 1 = − 1.Therefore, via simple comparison with [1] we have c ←→ 4R gg /3 = 4 • 1.16/3 = 1.54 (6).In addition for the coupling function f 0 we considered in this research.We'll leave c B as a model parameter to preserve an opportunity to fit magnetic field back reaction on the metric (2.6) and fix AdS-radius to be L = 1 in all numerical calculations.

Fig. 4
temperature function for different chemical potential values and p = 0.273
14.B for ⌫ = 4.5) shifts the scalar potential to larger absolute