Massive $2$-form field and holographic ferromagnetic phase transition

In this paper, we investigate in some detail the holographic ferromagnetic phase transition in an AdS${_4}$ black brane background by introducing a massive 2-form field coupled to the Maxwell field strength in the bulk. In the two probe limits, one is to neglect the back reaction of the 2-form field to the background geometry and to the Maxwell field, and the other to neglect the back reaction of both the Maxwell field and the 2-form field, we find that the spontaneous magnetization and the ferromagnetic phase transition always happen when the temperature gets low enough with similar critical behavior. We calculate the DC resistivity in a semi-analytical method in the second probe limit and find it behaves as the colossal magnetic resistance effect in some materials. In the case with the first probe limit, we obtain the off-shell free energy of the holographic model near the critical temperature and compare with the Ising-like model. We also study the back reaction effect and find that the phase transition is always second order. In addition, we find an analytical Reissner-Norstr\"om-like black brane solution in the Einstein-Maxwell-2-form field theory with a negative cosmological constant.

So far, most attentions about the duality application to condensed matter physics have been focused on the electronic properties of materials. In condensed matter physics, magnetism also plays an important role in materials including high temperature superconductors and heavy fermion metals in many strongly correlated electronic systems. The gauge/gravity duality provides an approach and perspective to understand these challenging propblems. Though there exist a few works involving magnetism in holographic superconductor models, magnetism does not pay the central role. Yet the scarcity of models on magnetism is due to various technical challenges in holographic context. In a previous work [23], we first proposed a new example of the application of the AdS/CFT correspondence by real-izing the paramagnetism/ferromagnetism phase transition in a dyonic Reissner-Nordatröm-AdS black brane. This model also was extended to realize the paramagnetism/antiferromangnetism phase transition by introducing two real antisymmetric tensor fields with interaction between them and they are coupled to the background gauge field strength [24]. In Ref. [25], by combining the holographic p-wave superconductor model [26] and the holographic ferromagnetism model, we studied the coexistence and competition of ferromagnetism and p-wave superconductivity. It is found that the results depend on the self-interaction of magnetic moment of the complex vector field and which phase (superconductivity or ferromagnetism) appears first. We noted that the ferromagnetic superconductivity has been also discussed in Ref. [27] by introducing two SU (2) gauge fields in the bulk. In that model superconductivity and ferromagnetism happens simultaneously and the magnetic susceptibility is finite at the ferromagnetic critical temperature. However, the ferromagnetic p-wave superconductor in heavy fermion systems such as UCoGe, URhGe and UGe 2 , shows the two critical temperatures are different in general [28] and magnetic susceptibility diverges at ferromagnetic critical temperature [29].
The reasons of using an antisymmetric tensor field to model ferromagnetic phase transition in holographic setup are as follows. On the viewpoint of symmetry breaking, ferromagnetic phase transition breaks the time reversal symmetry spontaneously in low temperature (if spatial dimension is more than 2, it also breaks spatial rotation symmetry), which in general is not associated with other symmetry breakings such as U(1), SU(2) and so on. Thus the dual operator does not carry U(1) or SU(2) charge, which implies that we need to use a real field corresponding to such an operator. From the point of view of covariance, magnetic field is the spatial component of a SO(1,3) tensor F µν , so magnetic moment should be also a component of a tensor. Furthermore, according to the origin of magnetic moment, magnetism of material is controlled by the intrinsic magnetism of electrons which is obtained from the Lagrangian for electrons, where A µ is the gauge potential, F µν is the gauge strength field and j µ is current density which is defined as Then the magnetic moment density − → N is the spatial component of an antisymmetric tensor field, Here i, j, k = 1, 2, 3. Thus we see that an effective field to describe magnetic moment in the boundary field in a covariant manner needs an antisymmetric tensor and its spatial component corresponds to the magnetic moment. These considerations can be generalized into 2+1 dimensions, where magnetic moment should still be regarded as the spatial components of an antisymmetric tensor field. Then according to the spirit of AdS/CFT, we can use an antisymmetric tensor field in the bulk to dual such a tensor operator.
Recently, it has been found in Ref. [30] that the original model in [23] contains a ghost and a new model has been proposed, there the exterior differential is used to replace the covariant derivative in the kinetic term of the antisymmetric tensor field. With such a simple modification it has been shown that there does not exist any ghost and causality violation does not appear in the new model, while the key results in the original model are kept qualitatively. While in Ref. [30] we paid our main attention to the health of the model and showed that the spontaneous magnetization can happen when the temperature gets low enough, in this paper, we are going to study this model in some detail with some directions. Before going on, we mention here that in this paper we assume the physical phase is homogeneous and will not consider inhomogeneous phases. This assumption is just for simplifying our discussions in technology. The possibility whether the inhomogeneous phase could appear spontaneously in this model under the homogeneous boundary conditions is left to study in the future.
This paper is organized as follows. In section II, we will describe the holographic model. We will give our ansatz for matter field and derive equations of motion in section III. In this section, we will show that there exists an analytical black hole solution different from the AdS Reissner-Nordtrsöm solution in this model and discuss associated properties of the solution. In sections IV and V, we will investigate the paramagnetism/ferromagnetism phase transition in two different probe limits. One is to neglect the back reaction of the 2-form field to the black brane geometry and to the Maxwell field, and the other to neglect the back reaction of both the Maxwell field and the form field. In the former case, we will also calculate the off-shell free energy of the holographic model near the phase transition point by using the Sturm-Liouville eigenvalue method and obtain a free energy form like the Ginzburg-Landau one, and make a comparison with the Ising like university class model. In the latter case, we calculate the DC resistivity in the ferromagnetic phase, which shows the behavior of the colossal magnetic resistance effect in condensed matter physics. We study the full back reaction effect in section VI by solving the full equations of motion of the model and find that the phase transition is aways second order. The summary and some discussions are included in section VII.

II. THE MODEL
In this paper, the model we are considering is just general relativity with a negative cosmological constant Λ = −3/L 2 , a U(1) field A µ and a massive 2-form field M µν in 4-dimension space-time. The ghost free action reads [30] where L is the AdS radius which is set to be unity and 2κ 2 = 16πG is associated to the gravitational constant in the bulk, which will be set to be unity in the following. g is the determinant of the bulk metric g µν . dM is the exterior differential of 2-form field M µν . m 2 = 0 is the squared mass of 2-form field M µν and must be greater than zero, which will be explained shortly. λ and J are two real model parameters with J < 0 in order the magnetization to happen spontaneously. λ 2 characterizes the back reaction of the 2-form field M µν to the background geometry and to the Maxwell field strenghth 1 . V (M µν ) is a nonlinear potential of the 2-form field. It describes the self-interaction of the polarization tensor, which should be expanded as the even power of M µν . In this model, we take the following form, Here * is the Hodge-star operator. The choice of nonlinear potential is not unique. We choose this form just for simplicity. As shown in Ref. [30], this potential shows a global minimum at some nonzero value of ρ. By varying action (3), we can get the equations of motion for the matter fields and gravitational field as The energy-momentum tensor T µν reads In the AdS/CFT correspondence, a hairy black hole with appropriate boundary conditions can be explained as a condensed phase of the dual field theory, while a black hole without hair is dual to a normal phase. Clearly when M µν vanishes, the model admits the AdS Reissner-Nordström (RN) black brane solution, which corresponds 1 Note that here λ can be also understood as the measurement of the coupling between the tensor field Mµν and the Maxwell field strength Fµν by rescaling the tensor field and the parameter J.
to the normal phase in the dual field theory. When M µν appears, we will see that there exists an analytical black brane solution differing from the AdS-RN solution. The new solution corresponds to the normal phase rather than the usual AdS RN solution. When we lower the temperature, the system exhibits an instability which triggers to break time reversal symmetry spontaneously as well as spatial rotation symmetry since the condensate will pick out one direction as special (if spatial dimension is more than 2) and the paramagnetism/ferromagnetism phase transition happens.

III. ANSATZ AND TRIVIAL SOLUTION
As we will discuss the full solution to the action (3) including the back reaction to the spacetime geometry, we start with the following ansatz for the metric and take the self-consistent ansatz for polarization field and U(1) field as with some real functions f (r), a(r), φ(r), p(r),and ρ(r). The bulk field B is a constant magnetic field, which can be regarded as external magnetic field in the dual boundary theory. We will denote the position of the horizon as r h and the conformal boundary will be at r → ∞. Since we would like to study a dual theory with finite chemical potential or charge density accompanied by a U(1) symmetry, we turn on A t in the bulk. Suppose the black brane horizon is at r h with f (r h ) = 0. According to gauge/gravity duality, the Hawking temperature of the black brane is identified with the temperature of boundary thermal state, which is given by and the thermal entropy S is given by the Bekenstein-Hawking entropy of the black brane where A denotes the area of the horizon and V 2 = dxdy. Put the ansatz (7) and (8) into equations (5), the in-dependent equations of motion read where a prime denotes the derivative with respect to r. The first three equations are for the polarization field and Maxwell field, and the last two are the two independent components of gravitational field equations. In fact, there are three nonzero components of gravitational field equations, but only two of them are independent due to the Bianchi identity. We are interested in the black brane configurations which have a regular event horizon located at r h . Therefore, in addition to f (r h ) = 0, one must require φ(r h ) = 0 in order for g µν A µ A ν being finite at the horizon. We require the regularity conditions at the horizon r = r h , which means that all functions will have finite values at r h and admit a series expansion in terms of (r − r h ). Then, at the horizon, we have the following relations, The black brane solution should be asymptotically AdS. Thus one has the following asymptotic solutions near the AdS boundary, where δ = √ 1 + 4m 2 and the dots stand for the higher order terms of 1/r, ρ ± and f 0 are all constants. The Breitenlohner-Freedman (BF) bound requires m 2 > − 1 4 according to the asymptotic solution of ρ. According to the AdS/CFT dictionary, up to a normalization, the coefficients µ and σ are directly related to the chemical potential and charge density in the dual system, respectively. The magnetic moment density in the dual theory is defined by the integration [23,30] Note that there is an additional restriction on m 2 when we treat the magnetic field B as the source of ρ. In order to keep ρ damping, we have to impose boundary condition for ρ such that ρ + = 0 and a restriction on the parameter such that m 2 > 0 [30], otherwise the integration (14) will diverge or the equation for ρ can not be linearized near the AdS boundary. These two restrictions can be understood in another way: when B = 0, the constant ρ + should be viewed as the source of the corresponding operator in the boundary field theory, according to AdS/CFT duality. In order the symmetry to be broken spontaneously, one has to turn off the source term.
In addition, in order to make ρ condense in the case without external magnetic field when the temperature is low enough, the model parameters should violate the BF bound for ρ near the horizon, whose geometry is an AdS 2 for an extremal black brane.
To see how this requirement restricts the parameters, we first consider the solution for (11) in the case of ρ(r) = B = 0 and a(r) = V (M ) = 0. The equations (9) become, The equations admit an analytical solution when m 2 = 0, i.e., Then the temperature of the black brane is Compared with the usual planar AdS RN black brane solution [31], we see that µ 2 in the metric of usual AdS RN solution is just replaced byμ in the new solution. Namely, the new solution has the same form as the planar AdS RN solution with the same Maxwell field φ(r), but also with a nontrivial profile of p(r) for the massive 2form field. Then we can easily calculate the free energy Ω and the charge density σ of the system as 2 , We see that the properties of this black hole solution depends on the value of 1 − λ 2 /4m 2 . When 1 − λ 2 /4m 2 = 0, the solution is just the planar AdS Schwarzschild black brane with nonzero U(1) gauge field and polarization field. To investigate the physical properties when 1−λ 2 /4m 2 = 0, it is useful to compute the partial derivative of charge density with respect to chemical potential, It is easy to see that Eq. (19) is positive when 1 − λ 2 /4m 2 > 0. However, when 1−λ 2 /4m 2 < 0, the value of (19) can become negative, which leads the dual boundary system to be chemical instability. To understand it, one can image that a box is submerged into the environment with fixed temperature T and chemical potential µ in equilibrium state. The energy and charge can exchange between the interior of the box and environment through the wall (see Fig. 1). Now because of thermal fluctuation, the chemical potential in the box then is µ − δµ < µ, which leads to some charge coming into the box. However, because Eqs. (19) is negative, adding the charge of the system will decrease the chemical potential of the system. Then we see that the chemical potential will decrease again, which will lead that more charge come into the box. So the charge inside the box will be more and more and the system is unstable. From Eq. (19), we can find the dual boundary system is chemical stable region only when temperature and chemical potential satisfy 3 , If Eq. (19) is positive, we see that adding charge will increase the chemical potential, so after some time, the system can be in equilibrium again. In fact, in the case of 1 − λ 2 /4m 2 > 0, the solution shares all the properties of usual AdS RN black brane except we have to replace µ 2 by µ. But in the case of 1 − λ 2 /4m 2 < 0, the properties are very different from those of usual AdS RN black brane. For example, the temperature gets increased when we increase the chemical potential while 2 The total charge density σ can be also calculated from the second equation in (5), which gives the same value as the one in (18). This also confirms that the chemical potential of the system is given by µ. 3 It is worth noting that this condition is only valid in grand canonical ensemble. In other ensembles, the stable condition is different in general. fixing the horizon radius and there does not exist zero temperature entropy.
In this paper, we only focus on the situation of 1 − λ 2 /4m 2 > 0, which can give a chemical stable dual boundary system and zero temperature black hole solution whose IR geometry has the form of asymptotic AdS 2 geometry. With this, we consider the possibility of spontaneous symmetry breaking of this system in low temperature. This can be analyzed by solving the following equation in the background (16), The condition of instability for ρ in some low temperature is that the BF bound in AdS 4 is retained but the BF bound in AdS 2 is violated. Near the horizon for an extremal black brane, the geometry is asymptotic AdS 2 .
In this region, the asymptotic solution for equation (15) is, So the conditions for spontaneous condensate are, One can immediately see that these inequalities have solution only when J < 0. Under the restriction of (23), there is a critical temperature, lower than which the nonzero ρ begins to appear. Then the full coupled equations of motion do not admit an analytical solution. Therefore, we have to solve them numerically. We will use shooting method to solve equations (11). In order to find the solutions for all the five functions F = {ρ, φ, p, a, f }, where we can directly solve φ(r) = ∞ r h dr(m 2 − 4Jρ 2 /r 4 )p and then put it into the equation of φ and obtain an equation for p, we must impose suitable boundary conditions at both the horizon r = r h and conformal boundary r → ∞. There are two kinds of scaling symmetries which are useful when we perform numerical computations: and Under the above two scaling symmetries, we finally have four independent parameters {a(r h ), ρ(r h ), p(r h )} and horizon radius r h at hand. In this paper, we will fix the boundary chemical potential to be unitary. When these four parameters are given, we can integrate the equations out of the horizon to obtain the whole solutions. When we perform the numerical computations, we can first set {r h = 1, a(r h ) = 0}. Then for a given ρ(r h ), we can adjust the value of p(r h ) to match the boundary condition ρ + = 0 at r → ∞. After solving the coupled differential equations, one should use the second scaling symmetry (26) to satisfy the asymptotic conditions a(∞) = 0. Then we use the first transformation (24) to fix the chemical potential for each solution the same. By this method, we get one-parameter solution ρ(r; T ) for T < T c . And then we can get the behavior of magnetic moment density N with respect to temperature T . When the parameters satisfy the inequalities (23), the solution of ρ = 0 will appear even in the case of B = 0. Because under the time reversal transformation, B → −B, in order to make action be invariant, we have to have ρ → −ρ. So when B = 0 but ρ = 0 in the source free case, the time reversal symmetry is broken spontaneously.
In order to see the main properties of the model, let us first study the model in probe limit for simplicity. Here we may consider two kinds of probe limit. The first one is to take the model parameter λ → 0 as we did in Ref. [23]. This case is just to neglect the back reaction of the massive 2-form field to the black brane geometry and the Maxwell field, but considers the effect of the Maxwell field to the background geometry. This probe limit corresponds to the case in such materials that their electromagnetic response properties are very weak compared with the external field and have little effects on the their transport properties. In other words, in this probe limit, the influence of external field on the materials is considered, but we neglect the back reaction of electromagnetic response to the external field and structures of materials such as crystal structure or energy band. The other is to neglect all back reaction of matter fields including the Maxwell field to the background geometry. In this probe limit, the interaction between the electromagnetic response and external field is taken into account so that we can study how spontaneous magnetization influences the electric transport, but they both have little influence on the structures of materials. We will study two cases in the following sections separately.

A. Spontaneous magnetization and susceptibility
Let us first investigate the spontaneous magnetization in the limit of λ → 0. In this limit, we neglect the back reaction of polarization field to the gauge field and background geometry. The background geometry and the Maxwell field can be taken as where we have set the horizon radius r h = 1. In this case, the equations of the polarization field read Note that the expression of p(r) can be solved directly. We put it into the equation of ρ(r), and get This equation shows that ρ will be spontaneously condensed below a critical temperature only when m 2 > 0. When m 2 > 0, near the critical temperature where ρ is very small and therefore the term ρ 2 can be neglected in (29), thus increasing the chemical potential 4 will decrease the effective mass at the horizon (note that J < 0), which leads that ρ can be condensed spontaneously below a critical temperature. 5 Furthermore, the restrictions on the parameters in the case of λ → 0 are, In what follows, since different parameters will give similar results, we will choose m 2 = −J = 1/8 as a typical example when we perform numerical computations.
To study the spontaneous magnetization, we need set B = 0. Since near the critical temperature, ρ is very small, the nonlinear term of ρ can be neglected in (29).To find the critical temperature, we can solve the linearized equation of ρ with initial condition at the horizon r = 1 and boundary condition ρ + = 0 at r → ∞. Without loss the generality, we take the initial value of ρ at horizon to be unity, and treat T as the shooting parameter to match the source free condition. There will be many solutions for shooting parameters T , we choose the highest one as the critical temperature T c . We find the critical temperature T c /µ 1.78 for the case of given model parameters.
When the temperature is lower than the critical temperature T c , in order to examine whether the polarization field ρ can make spontaneous magnetization when the external magnetic field B = 0. We plot the value of ρ + versus ρ(r h ) at the horizon in the top panel of Fig. 2. Each curve in this plot corresponds to different temperature. This plot shows a typical example in the high and low temperature cases which correspond to the red and the black lines, respectively, while the case with the critical temperature is shown by the blue curve. In the case of high temperature, we see that the curve has no intersecting point with horizontal axis except a trivial point at the origin, which corresponds to a trivial solution with ρ = 0 and it describes the paramagnetic phase. On the other hand, when the temperature is low enough, we find that there exists a nontrivial solution which locates at ρ(r h ) = 0. This solution breaks the time reversal symmetry and the system gets into a ferromagnetic phase. As a result, we see that the model indeed can give rise to a paramagnetism/ferromagnetism phase transition in the case without external magnetic field.
Furthermore, when the temperature is lower than the critical temperature T c , we have to solve the equation (29) with the condition of ρ + = 0 by shooting method. The spontaneous magnetization N is defined by the integration (14) with a = 0. The spontaneous magnetic moment with respect to temperature is shown in the bottom panel of Fig. 2. In addition, since the magnetic moment of the polarization field obtains an expectation value, the time reversal symmetry is broken spontaneously.
By fitting this curve near the critical temperature, we find that there is a square root behavior for the magnetic moment versus temperature, which is a typical behavior for a second order phase transition, specifically, for m 2 = −J = 1/8, we have This gives the critical exponent 1/2, the same as the one from the mean field theory. Except for the magnetic moments which is one of the characteristic properties of ferromagnetic material, another remarkable one is the behavior of susceptibility density of the material in the external magnetic field. The static susceptibility density is defined by When we turn on the external magnetic field B, the function ρ is nonzero in any temperature. In order to compute the susceptibility density, we need to shoot for the boundary conditions with one parameter ρ(r h ) for equation (29) under the given external magnetic B and temperature T . From the definition of χ, which involves only the behavior of B → 0, so the computation can be simplified in the following way.
In the case of T > T c , considering the result ρ = 0 when B = 0, we except that ρ is in the same order as B, then ρ 2 is as the same order as B 2 . So in the case with weak magnetic field, we can neglect nonlinear terms in equation (29). In that case, we just need to solve the "linearized" equation Because ρ ∝ B, we can set B = 1, Solving the equation (34) with the boundary condition ρ + = 0, the susceptibility density then can be computed by In Fig. 3, we show the magnetic susceptibility as a function of temperature and find it satisfies the Curie-Weiss law of ferromagnetism in the region of T → T + c . Concretely, for the chosen model parameters, we have So we can conclude that the dual system is in a paramagnetic phase in high temperature and ferromagnetic phase in low temperature. The model can describe a paramagnetism/ferromagnetism phase transition. We will see later that the exponents in Eqs. (32) and (36) are exact and can be obtained by analytical methods.

B. Holographic Ginzburg-Landau formulation
In action (3), two parameters appear in the Lagrangian. One is the mass of the polarization field, which corresponds to the conformal dimension of the dual operator by the standard dictionary. The other is the selfinteraction coupling constant J. So a natural question is what the meaning of this constant is in the dual boundary theory? In this subsection, we will try to give a quantitative interpretation. Since we are working in the probe limit, the geometry and external Maxwell field are fixed, this leads to a simplification to compute the partition function of the bulk theory. This is the reason that we consider this kind of probe limit here. Our method is to compute the effective grand thermodynamic potential in both sides and to equate them (which are equivalent to compute the partition functions in both sides and to equate them). Then we can "read off" the meaning of J in the dual boundary theory.
We assume that the the spontaneous magnetization in the boundary theory can be described by so-called Ising universality class model at least in the vicinity of critical temperature, whose Hamiltonian is, where r, r are the positions of lattices, J (r − r ) > 0 is the exchange integration of lattices r and r , s(r) the z-component of spin at lattice r, K is a positive constant which corresponds to the deviation from the Ising model. When K → ∞, Hamiltonian (37) reduces to that of the usual Ising model.
In this paper, the dual boundary is a (2+1)dimensional spacetime, i.e., a film system. Then the spin (or magnetic moment) in fact is just a pseudoscalar. So we have − → s = s. To compare with our holographic model, it is convenient to change Hamiltonian (37) into the one in continuous limit, where l is the lattice spacing and J = r J (r), R 2 J = r r 2 J (r). The summations for J and R 2 J are only in one crystal lattice. We see that the Ising like Hamiltonian is a λφ 4 theory if we use s = s √ J R J . From Appendix A, we can find that the grand thermodynamic potential in the mean field approximation in the high temperature limit reads, with Here we have used s cl to represent the classical value of s. We see that there is a critical temperature T c = 4πm 2 s /λ s γ E . When T > T c , the thermodynamic equilibrium phase corresponds to s cl = 0, while T < T c , the thermodynamic equilibrium phase has s cl = 0. Now let us consider the gravity side of the holographic model to compute the grand thermodynamic potential. It is convenient to make a coordinate transformation by z = r h /r. As p can be solved directly, then we put it into the equation of ρ(r), and get As we will care about the behavior of T → T c , the value of ρ will be a small quantity near the transition point.
In this case, we can make a Taylor's expansion on the nonlinear term of ρ in Eq. (41) as, Neglecting the high order terms, Eq. (41) can be rewritten as Up to the order of ρ 4 , the part of polarization field in action (3) can be written as, which is a function of T and B, but a functional of ρ. The asymptotic solution for (43) is The source free condition is ρ + = 0 as z → 0 + . Under this, the grand thermodynamic potential or free energy in grand canonical ensemble Ω is, According to thermodynamic relationship, It seems that the magnetic moment should be, However, in our previous papers, we defined the magnetic moment as, There is a difference factor 1/2 between (49) and (48). Now we show that the expression (48) is not true. The reason is as follows. The relation (47) is on-shell, while the equation (46) is off-shell. In order to obtain the differential relation of Ω(T, B; ρ) with respect to (T, B, ρ), we need to use the Euler homogenous function theorem. We should first note that under the scaling transformation z → kz, (t, x, y) → k(t, x, y), we have Ω(kT, k 2 B; k 2 ρ) = kΩ(T, B; ρ), which gives, Submitting (46) into (50) and considering the on-shell condition δΩ/δρ = 0, we find that, This is just the definition (49) .
The key step for computing the grand thermodynamic potential is to structure the Sturm-Liouville problem 6 , which is the following ODE: with the boundary conditions: The weight function ω(z) can be an arbitrary positive continuous function in the region of [0, z h ]. From a practical point of view, we choose weight function such that the values of λ n will not influence the asymptotic behaviors of equation (52) when z → 0 + . There are many choices for weight function. Here we choose ω(r) = z k with an integer k > 2.
Once note that the asymptotic solution when r → ∞ for equation (52) is, one can find that the boundary condition (b) corresponds to ρ + = 0. Let L 2 ([0, z h ], ω(z), dz) be the Hilbert space of square integrable functions on [0, z h ], i.e., with the inner product and D be the subspace of L 2 ([0, z h ], ω(z), dz) that satisfies the boundary conditions of (a) and (b), i.e., Then we can prove that P is the self-adjoint operator on D, i.e., According to the properties of SL problem, the solutions of (52) form a function basis on D with which one can expand any functions belonging to D, i.e., and with c n = ρ n , h . Note that all the parameters in equation (52) depend on temperature, so do the eigenvalues λ n and eigenfunctions ρ n . Indeed, the minimal eigenvalue, i.e., the first eigenvalue λ 1 is the function of temperature. There is a critical temperature T c , at which we have λ 1 = 0. We can find that when T > T c , λ 1 > 0. In order to show this, we take m 2 = 1/8 as an example to plot ρ + with respect to λ. The boundary condition (b) is equivalent to ρ + (λ n ) = 0. Fig. 4 shows the value of y(λ) = arctan[10ρ + (λ)] with respect to λ 1/3 (it is just for convenience). The x-coordinate values of zero points (y = 0) give the eigenvalue λ n . We can see that all the eigenvalues are positive when T > T c and λ 1 < 0 when T < T c . For any temperature, the equation has infinite eigenvalues and the smallest eigenvalue exists, which shows the system is stable.
Let us now turn our attention to the free energy (46). For convenience, we will use scaling transformation to set z h = 1 in the process of computation, and then transform into the case of fixing chemical potential in the final results.
Let ρ(r) = ρ(r) + B/m 2 be any function configuration belonging to D, in which ρ(r) dose't need to be the solution of EoM (43). We can use the eigenfunction ρ n to expand ρ(r) and magnetic moment as, where c n and N n are coefficients, defined as, Here we have assumed that {ρ n } is an unit base. Then the variational principle of Ω(T, B; ρ) underlying the equations of motion, or finding a solution of EoM (43), is equivalent to minimizing Ω(T, B, c n ) with respect to c n 's. Let us consider the case of spontaneous magnetization, i.e., the case with B = 0. In this case, we have with ρ = ∞ n=1 c n ρ n . Using the orthogonal relationship, we have, If T > T c , then λ n > 0. Because of J f < 0, we can find that Ω(T, c n ) ≥ 0. The minimization of Ω(T, c n ) = 0 is achieved only when c n = 0, i.e., ρ = 0. So the nonzero solution appears only when λ 1 < 0, i.e., T < T c . This is just what we have obtained in the pervious section. When T → T − c , we can set λ 1 = a 0 (T /T c −1) with a 0 > 0 and assume that the off-shell solution is dominated by the first term in (60) only, i.e., |c 1 | c n for n ≥ 2 in (61). As a result, we have, with a 1 = 1 4 1 0 ρ 4 1 z 8 dz| T =Tc > 0 and, Put (65) into (64), we can obtain, We can see that this is just the Ginzburg-Landau (GL) theory of ferromagnetic model. Now the grand thermodynamic potentials of the holographic model in the gravity side and in the boundary theory side are in hand, we can use the AdS/CFT duality to relate them, which gives, Here we should mention that the dual theory in the boundary being an Ising-like model is an assumption. With this equality, At the critical temperature, comparing (66) with (39), we can obtain, Note that in materials, the spontaneous magnetization is proportional to the expectation value of z-component of spin, i.e., N ∝ s cl . Following the definition of J f in Eqs. (43) and λ s in Eqs. (40), we see that if J = 0, then Eq. (68) leads to K = 0. Therefore we can interpret the parameter J in action (3) as the deviation from the standard Ising model in the boundary theory.
Here some remarks are in order. First, the equality Eq. (67) is exact according to the AdS/CFT correspondence. However, the expressions for the grand thermodynamic potentials in both sides are just some approximations. In gravity side, because we used the probe limit which neglects the back reaction of polarization field to the background geometry and to external field and we only computed the grand thermodynamic potential in the classical level which neglects the quantum correction. In the boundary side, we used mean field expansion and only took the tree level of the quantum fluctuation into account. So the relationship (68) is only approximately valid and our explanation for J is only a qualitative description. Second, the Ising-like model (37) can only describe the local moment system, which is usually suitable for insulator. For metal magnetic materials, we need to use different model. So what the meaning of J is in these materials needs to be considered in the future.
With the grand thermodynamic potential in Eq. (66), we can obtain the expression of magnetic moment in the ferromagnetic phase as, This confirmed the critical behavior obtained in the numerical calculations and the critical exponent 1/2 is an exact result. We can compute all the coefficients appearing in (69) and compare them with the numerical ones. We can also compute the grand thermodynamic potential when magnetic field is nonzero and get the magnetic susceptibility and hysteresis loop. The details will be given in Appendix B.

A. Spontaneous magnetization and susceptibility
In the probe limit by neglecting the back reaction of all matter fields including the Maxwell field, the background geometry is just the planar AdS Schwarzschild black brane with, In the background, the equations for Maxwell field read, We can directly obtain φ(r) = ∞ r h dr(m 2 − 4Jρ 2 /r 4 )p with φ(r h ) = 0 and put it into the equation of φ , then the equation for p can be obtained. So in this probe limit, we need only solve the three equations of matter fields numerically, Note that in this case, there does not exist the AdS 2 geometry, the near horizon geometry of an extremal black brane with vanishing temperature. To have the spontaneous magnetization when the temperature is lowered, the restrictions (23) need to be reconsidered. To find the restriction about the parameters, let us consider Eqs. (72) in the high temperature region where ρ vanishes and the solutions for φ and p are the same expressions shown in Eqs. (16). We can read off the effective mass square of ρ at the horizon as, (73) Because of J < 0, the temperature term contributes a negative term into the effective mass square, which is divergent when T → 0. Thus we see that in the grand canonical ensemble, the instability always appears provided that the temperature is low enough. As a result, in this case, we need not the restrictions in (23). But, the parameters have to satisfy the condition λ 2 < 4m 2 . To see this, let us recall the charge density in Eq. (18). In this probe limit, we have, Then we reach, We see that it is positive only when λ 2 < 4m 2 . Now we consider the spontaneous magnetization in this probe limit. First, let us compute the critical temperature T c when B = 0. Similar to the case in the first kind of probe limit, the polarization field ρ is a small quantity near the critical temperature, we can neglect the nonlinear terms of ρ in the equations of φ and p. Then we can get φ(r) = µ(1 − r h /r) and p(r) = µr h /r 2 m 2 which are identical with expressions in Eqs. (16). The equation of ρ is just Eq. (21). At the horizon, the initial conditions are, When we perform the numerical computation, we can first fix the horizon radius r h = 1 so the temperature is also fixed. By adjusting the chemical potential, we shoot the boundary condition ρ + = 0. In the end, we can use the scaling transformations (24) and (26) to transform our results into the case in grand canonical ensemble where the chemical potential is fixed. As a typical example, we also choose parameters as m 2 = −J = 1/8 and λ = 1/2. The critical temperature is T c /µ 1.7871. Similarly, we can plot the relationship between ρ + and shooting parameter µ, in order to examine whether ρ gets spontaneous condensation when T < T c . We find that the solution of source free always appears, which results in the spontaneous magnetization of the system, and breaks the time reversal symmetry in low temperatures.
When the temperature is lower than the critical one T c , we have to solve Eq. (72) to get the solution of the order parameter ρ, and then compute the value of magnetic moment N , which is also defined by Eq. (14) with a = 0. In the top panel in Fig. 5, we plot the value of magnetic moment N as a function of temperature. We see that when the temperature is lower than T c , the non-trivial solution ρ = 0 and spontaneous magnetic moment appears. It corresponds a time reversal symmetry breaking spontaneously. In addition, let us stress here that if the boundary spatial dimension is three, the spatial rotational symmetry is also broken spontaneously, since a nonvanishing magnetic moment chooses a direction as a special. The numerical results show that this phase transition is a second order one with the behavior N ∝ 1 − T /T c near the critical temperature. The result is still consistent with the mean field theory description of the paramagnetism/ferromagnetism phase transition.
Next we calculate the static magnetic susceptibility in this probe limit, defined by Eq. (33). Based on the previous analysis, the magnetic susceptibility is still obtained by solving Eq. (34), and the only difference is the form of f (r). Thus we can also set the magnetic field B

B. DC conductivity in the ferromagnetic phase
The electric transport is also an important property in the materials involving spontaneous magnetization. Now let us study how the DC conductivity is influenced by spontaneous magnetization in this model. In order to simplify our computation in technology, we will work in the probe limit by neglecting back reactions of all the matter fields. This limit can give out the main features near the critical temperature. However, in the case of near zero temperature, we have to consider the model with full back reaction. We will consider this in the future.
To compute the conductivity, we have to consider some perturbations for gauge field with harmonically time varying electric field. Due to the planar symmetry at the boundary, the conductivity is isotropic. Thus for simplicity, we just compute the conductivity along the x-direction. According to the dictionary of AdS/CFT, we consider the perturbation δA x = a x (r)e −iωt . In the probe limit, this perturbation will also lead to the perturbations of polarization field in the first order of . As a result, we have to consider the perturbations for all the components of gauge field and polarization field. However, if we only care the conductivity in the low frequency limit, i.e., T ω → 0, the problem can be simplified. In the low frequency limit, we only need turn on the three perturbations, and corresponding equations for the three perturbations in the low frequency limit read [ with p and ρ determined by Eqs. (72). Here O(ω) is the terms with order of ω which can be neglected when ω → 0. In general, the term ω 2 /r 2 f (r) can not be neglected since f (r) is zero at the horizon, which leads to that the limit of ω → 0 is ambiguous. However, at the horizon, if we impose the ingoing conditions for C rx , C ty and a x , with r * = dr/(r 2 f ), we find the system has a welldefined limit when ω → 0 if T = 0. At the AdS boundary with the source free condition, we have the following asymptotic solutions, Then the gauge/gravity duality implies that electric current J x = a x− and the DC conductivity is given by, As a holographic application of the membrane paradigm of black holes, we can direct obtain the DC conductivity from Eqs. (78) using the method proposed by Iqbal and Liu in [34]. In fact, the transport coefficients in the dual field theory can be obtained from the horizon geometry of the dual gravity in the low frequency limit. Applying this into U (1) gauge field, this conclusion implies that the DC conductivity is given by the coefficient of the gauge field kinetic term evaluated at the horizon. To see this, we assume that T > 0 and ω → 0, then we can neglect all the terms of ω in Eqs. (78). We first note that, Eq. (78c) shows that this quantity is conserved along the direction r. So at the horizon, using Eqs. (78b) and (78c), we have, (83) Combining Eqs. (78a) and (78b) and considering the fact that C ty is regular at the horizon, we have, at r → r + h . Thus we have from Eqs. (84) and (83) that Now let us take the ingoing condition for a x at the horizon, which tells us that, finally we get, .
Here p 0 and ρ 0 are the initial values of p(r) and ρ(r) at the horizon, which can be computed from Eqs. (72). In the low frequency limit, Eqs. (78) imply that the electric field is constant, i.e., lim r=r h a x (r) = a x+ . So we have seen that we can obtain the DC conductivity near the horizon as, In the top panel of Fig. 6, we plot the numerical result on the DC resistivity 1/σ as a function of temperature. We see that the resistivity shows a metallic behavior when the temperature is below the Curie temperature.
With the appearance of ferromagnetism, DC resistivity decreases when the sample gets cooling, which shows in many interesting phenomena in condensed matter physics, especially in a class of manganese oxides which are widespread because of the discovery of colossal magnetoresistance (CMR) and receive a lot of interest both in theory and experiment [36,37]. Note that this effect has a completely different physical origin from the "giant" magnetoresistance observed in layered and clustered compounds. In recent twenty years, CMR is among the main topics of study within the area of strongly correlated electron systems and its popularity is reaching the level comparable to that of the high-temperature superconducting cuprates. In the bottom panel of Fig. 6, we show the experimental data from a typical CMR material La 1−x Sr x MnO 3 as an example. We see that our model gives a very similar behavior as the latter with x ≥ 0.175. Of course, we should mention here that there still exist some differences between our model result and experimental data on CMR. In general, when T > T c , the material shows a semiconductor or insulator behavior and the DC resistivity increases with cooling the sample. In our model, however, the DC resistivity is a constant when T > T c . So this model only gives partial property of CMR when T < T c . But this is an exciting and enlightening result, because it implies that this model can lead to a possibility to build a holographic CMR model and to investigate this typical and important strong correlated electrons system in the AdS/CFT setup. We are going to investigate this issue in the future.

VI. PHASE TRANSITION WITH BACK REACTION
In previous sections, we have studied the spontaneous magnetization in two kinds of probe limit. However, probe limit may lead to some information lost. For example, in some holographic superconductor models, it will lead to the appearance of the first or zeroth order phase transition when the strength of back reaction gets beyond some values (see Refs. [26,38,39], for example). In addition, when the temperature is low enough, the probe limit may lose its validness. To have a complete phase diagram for the holographic model, we need to go beyond the probe approximation and to include the back reaction.
A. On-shell free energy The model admits various solutions, in order to determine which phase is thermodynamically favored, we should calculate the free energy of the system for both normal phase and condensed phase. In gauge/gravity duality the grand potential Ω of the boundary thermal state is identified with temperature T times the on-shell bulk action with Euclidean signature. Since we are considering a stationary problem, the Euclidean action is related to the Minkowski one by a minus sign as where L m = −F µν F µν + λ 2 L 2 and g is the determinant of the metric. We first show that, when evaluated on a solution, this action reduces to a simple surface term at the AdS boundary. From the symmetries of the solution (5) and (6), the yy component of the stress energy tensor only has a contribution from the terms proportional to the metric. Thus, the gravitational filed equations imply that The Euclidean action is then The surface term on the horizon vanishes since f (r h ) = 0. Thus we can get the Euclidean action that contains the surface term at r ∞ as, As the first item of Eq. (92) diverges when r → ∞ and must be regulated. This counter terms we need to regulate the action are the standard ones (see for example [25]): where g ∞ is the induced metric on the boundary r = r ∞ and K = g µν ∞ ∇ µ n ν is the trace of the extrinsic curvature (n µ is the outward pointing unit normal vector to the boundary). The summation S on-shell Euclidean = S E + S c.t. is now finite in the limit r → ∞. The regularized action becomes, after considering the asymptotical forms in (13),

B. Tensor hairy solutions and phase transition
We are interesting in the black brane solutions with nontrivial space-space component ρ of the tensor field. For this, we have to adopt the numerical method to find such solutions. Without loss of generality, the location of r h can be fixed to be unity in our numerical calculation. We are then left with two independent parameters {ρ(r h ), p(r h )}. By choosing p(r h ) as the shooting parameter to match the source free condition at r → ∞, i.e., ρ + = 0, we finally have a one-parameter family of solutions labeled by the value of ρ(r h ). Other coefficients can be expressed in terms of those parameters. After solving the set of equations, we can calculate the spontaneous magnetization N and free energy density.
With a fixed m 2 , we scan a wide range of J and λ within the limitation Eq.(23) in 3-dimension plane in order to trace out the evolution of critical temperature T c versus these parameters. Fig. 7 plots the critical temperature T c as a function of J and λ in the cases with m = 1/4, 1/8. 1/16 and 1/160, respectively. We can see from Fig. 8 that the critical temperature is weakly dependent of the parameter λ. Moreover, the smaller the value of m 2 , the larger the phase transition temperature T c when the same value of J and λ are considered. For each value of λ, the analytical black brane solution (16) always exists. However, for sufficiently low temperature, we always find additional solutions with non-vanishing ρ which are thermodynamically favored. That is to say, for each value of λ we take, there is a phase transition occurring at a certain temperature T c where the black brane developing a new "tensor hair" with nontrival ρ becomes thermodynamically favored. In the dual field theory side, it means that magnetic moment acquires a vacuum expectation value breaking the time reversal symmetry spontaneously.
The top panel of Fig. 9 presents the condensate as a function of temperature for λ = 1/8, 1/2 and 7/10, from  which one can see that N rises continuously from zero at T c . For small λ, the curve is similar to the case with the probe limit (compare the case of λ = 1/8 in Fig. 9 and the plot in the bottom panel of Fig. 2), We can see from the figure that when λ gets larger, the condensate increases. However, near the critical temperature, the square root behavior still holds as The free energy difference of the condensed phase and the normal phase is expressed by ∆Ω, which is plotted in the bottom panel of Fig. 9. It is obvious that below the critical temperature T c , the state with non-vanishing magnetic "tensor hair" is indeed thermodynamically favored over the normal phase because ∆Ω always be less than zero. Moreover, our numerical calculation indicates that the order of the phase transition is only second order, no matter how we increase the strength of the back reaction in parameters allowed region, i.e., λ 2 /4m 2 < 1.
The first order or even the zero order phase transition does not appear in this model.

VII. SUMMARY AND DISCUSSION
In this paper we have presented a holographic model to realize the paramagnetism/ferromagnetism phase transition in AdS black brane background by introducing a massive 2-form field in the bulk. This 2-form field couples to the background Maxwell field strength and carries self interaction.
The model admits a new analytical black brane solution with a non-trivial time-space component of the tensor field. The properties of the black brane solution depend on the value of λ 2 /4m 2 . When λ 2 /4m 2 < 1, this new solution is very similar to the planar AdS RN black hole. But if λ 2 /4m 2 > 1, the black brane solution is chemically unstable in grand canonical ensemble unless the temperature and chemical potential satisfy some additional conditions. In that case, there is not corresponding extreme black hole, i.e., the horizon radius is zero when temperature goes to zero. A very special case is that when λ 2 /4m 2 = 1, the space-time geometry is just the planar AdS Schwarzschild geometry but both the U(1) field and 2-form field nonvanish.
For our goal in this paper we focused on the case with λ 2 /4m 2 < 1 so that an asymptotic AdS 2 geometry near the horizon emerges when the temperature tends to zero. By this emergent AdS 2 geometry, we obtained the condition in which the spontaneously symmetry breaking can happen. If the parameters satisfy the condition, time reversal symmetry will be broken spontaneously and the paramagnetism/ferromagnetism phase transition can happen when the temperature is lower than a critical value.
In order to understand the properties of this holographic ferromagnetic phase transition, we investigated the paramagnetism/ferromagnetism phase transition in two kinds of probe limit and in the case with full back reaction, respectively.
In the case of the first kind of probe limit where the model parameters λ → 0. This probe limit neglects the back reactions of the 2-form field to the Maxwell field and background geometry and makes it simple to study the behavior of spontaneous magnetization. In this probe limit, we computed the critical exponents by both numerical and analytical approaches, which are agreement with the mean field results. In addition, we obtained a Ginzburg-Landau-like free energy near the critical temperature for the holographic model and argued that it is agreement with results from the Ising-like Hamiltonian in film system of the condensed matter physics and related the model parameter J to the parameters in the Ising-like model.
In the second kind of probe limit, we neglected all back reaction of matter fields to the background geometry but considered the interaction between the tensor field and the Maxwell field. In this probe limit, we are able to study the influence of spontaneous magnetization on the electric transport properties in a relatively simple way, where the background geometry is fixed. We found that the critical exponents are the same as ones in the first kind of probe limit. We also computed the DC resistivity in this probe limit. It was found that the DC resistivity is suppressed by spontaneous magnetization and shows a metallic behavior. This is very similar to the strong correlated phenomenon named CMR effect found in the some manganites.
Next we considered the case with full back reaction and solved the full equations of motion numerically. It was found that the free energy difference between the condense phase and normal phase, ∆Ω, is zero at critical temperature and always negative when T < T c . The phase transition is always a second order one as one increases the strength of the back reaction.
Main calculations are made in 3 + 1 dimensions in this paper, but it can be easily extended into the higher dimension case. In the latter case, the space rotation symmetry will be broken spontaneously when the phase transition happens. In addition, we can also generalize this model to the case with the Lifshitz symmetry in the bulk and study the influences of the Lifshitz dynamical exponent z on the condensate both in the probe limit and in the case with the back reaction. In addition, in all the calculations in this paper, we assumed that the solution in the bulk or the phase at the boundary is homogeneous, which is a strong assumption. In fact, inhomogeneous solution may exist even in a model whose Lagrangian has translation symmetry and the boundary conditions are homogeneous (see Ref. [17], for example). In many materials, the inhomogeneous phase can appear spontaneously in a chemical homogeneous materials. Therefore it is of great interest to study whether the inhomogeneity could appear spontaneously in this model. Although we focused on the ferromagnetic phase transition only in this paper, which has been understood well in condensed matter physics, it offers a framework rather than only a specific model, which can be regarded as basic starting point to understand more complicated phenomenon involving spontaneous magnetization. As a result, there are various prospects to study and we expect more exciting results could be reported in the future. the coefficients in (69) and compare with the numerical results in the previous sections. Let us first compute N 1 and a 1 . For this we have to first find the eigenfunction ρ 1 , which is the solution of, at T = T c with the conditions, For convenience, here we do not assume that {ρ n } form an unit base. Thus we have, here C 1 is the normalization coefficient and In order to compute a 0 , we need to solve equation (B1) in the limit T → T − c . To clarify that the results are independent of the specific form of weight function, we choose k = 3, 4 as two examples. Then we fit the relation λ 1 = a 0 (T /T c − 1) to find a 0 . Figure 10 shows that λ 1 and T /T c − 1 indeed satisfies a linear relation very well. Numerical results show that N 1 2.0412, a 1 0.7276,a 0 1.1309 for k = 4 and N 1 1.8370, a 1 0.4772,a 0 0.9155 for k = 3. We have, with a 2 4.966 for k = 4 and a 2 4.964 for k = 3. We see that different weight functions give different values for N 1 , a 1 and a 0 , but the same value of magnetic moment N (up to a numerical error).
The value of a 0 can also be obtained directly by solving ODE (52). In the region near the critical temperature, we assum λ 1 = a 0 (T /T c −1). Note that all quantities in (52) are the functions of temperature, so take derivative with respect to T and evaluate at T = T c , we get, Here ρ 1 is the eigenfunction of (B1). Now treat ρ T = dρ1 dT as an unknown function to be solved, then the task to find a 0 becomes to solve a non-homogenous eigenvalue problem, At the AdS boundary, ρ T has the same asymptotic behavior as (53), thus we can impose the boundary conditions as |ρ T (1)| < ∞, ρ T + = 0.
We find that ρ T ∈ D. We then use the basis {ρ n } to expand ρ T , i.e., Here C n are the modules of ρ n . Using the fact λ 1 = 0 at T = T c and It is very useful to find its equivalent form in the case with fixed r h = 1, since it is convenient when we perform numerical computation. If we fix r h = 1, the shooting parameter is chemical potential µ. The relation between temperature in grand canonical ensemble and chemical potential is given by, (1 − T /T c ).
(B15) We can see that it is determined by the equation (B1), but independent of the weight function! The expression (B14) depends on the weight function, because it depends on C 1 . In order to check the formula (B14), let us compute the values of a 0 in the cases of k = 4 and k = 3. The results gives 1.1286 and 0.9140, respectively. We see that, up to numerical errors, they are the same as what we have obtained by fitting the curve in figure 10. In addition, the expression (B15) gives a 2 = 4.9560, which is very close to what we obtained in the numerical calculation.

Susceptibility and hysteresis loop
When B = 0, the susceptibility for T > T c is defined as, In the case with T > T c and B → 0, we can neglect the non-linear term, i.e., setting J f = 0. The solution of equation (43)  (B18) Multiplying a factor ρ n /C n and integrating the above equation from 0 to 1, we can obtain c n = Bγ n λ n , with γ n = 1 0 4Jµ 2 z 4 C n m 6 ρ n dz.
Thus we can get the magnetic moment density as and the magnetic susceptibility When T → T + c , we have λ 1 = a 0 (T /T c −1) → 0 + . Thus χ is dominated by the first term in the summation of (B21) and its inverse can be expressed as In the case of m 2 = −J = 1/8, we have λ 2 χ −1 /µ c 4.0520(T /T c − 1), which is very close to our numerical result λ 2 χ −1 /µ c 4.0499(T /T c − 1) given in the numerical calculation. Next let us move to the case with B = 0. In this case from (43) we have, ρ n J f ρ 3 z 8 dz = 0, n = 1, 2, · · · .
(B25) For convenience, we assume that {ρ n } is an unit base, i.e., C n =1. The equations (B25) are equivalent to (43) if we take all the terms in (60) into account. In the case of T → T − c , assume that the first term in (60) dominates only, i.e., |c 1 | c n for n ≥ 2 in (61), we have Taking n = 1 in (B25), we have, c 1 λ 1 +Bγ 1 −c 3 1 J f 1 0 ρ 4 1 z 8 dz = c 1 λ 1 +Bγ 1 −4c 3 1 J f a 1 = 0. (B27) For a given temperature T → T c , we can combine (B26) and (B27) to obtain a relation between external magnetic field B and magnetic moment density N . Figure 11 shows the results with T = 1.05T c , T = 0.9T c and T = T c , respectively, in the case of m 2 = −J = 1/8. We see that it is very similar to what we have obtained in the previous work [23].
Finally, let us notice that in GL theory, the equation for magnetic moment density is, with two positive coefficients A 1 and A 2 . However, it is easy to see that the equation for N in our model is different from the usual form (B28) from the GL theory, which can be obtained by combining (B26) and (B27) to eliminate c 1 . Namely, alhought our model gives the similar results near the critical temperature in GL theory, there exist some differences between the holographic model and the GL theory even in the region near the critical temperature.