Holographic DC Conductivity for Backreacted Nonlinear Electrodynamics with Momentum Dissipation

We consider a holographic model with the charge current dual to a general nonlinear electrodynamics (NLED) field. Taking into account the backreaction of the NLED field on the geometry and introducing axionic scalars to generate momentum dissipation, we obtain expressions for DC conductivities with a finite magnetic field. The properties of the in-plane resistance are examined in several NLED models. For Maxwell-Chern-Simons electrodynamics, negative magneto-resistance and Mott-like behavior could appear in some parameter space region. Depending on the sign of the parameters, we expect the NLED models to mimic some type of weak or strong interactions between electrons. In the latter case, negative magneto-resistance and Mott-like behavior can be realized at low temperatures. Moreover, the Mott insulator to metal transition induced by a magnetic field is also observed at low temperatures.


I. INTRODUCTION
Gauge/gravity duality [1-3] has provided powerful tools for exploring the behavior of strongly coupled quantum phases of matter, and some remarkable progresses have been made [4][5][6][7][8]. Conductivity is an important transport quantity in condensed matter, and the gauge/gravity duality provides a framework to compute it for strongly interacting field theories.
Studying the behavior of the conductivity in the presence of external magnetic fields can help us to better understand the transport properties of materials. For normal metals, the resistance is a monotonically increasing function of the magnetic field [9], which appears as positive magneto-resistance. However, negative magneto-resistance has been observed in several experiments [10][11][12]. On the other hand, the behavior of negative magneto-resistance was found in strongly coupled holographic chiral anomalous systems [13][14][15][16]. In [17], it showed that negative magneto-resistance could also arise in nonanomalous relativistic fluids due to the distinctive gradient expansion. Note that the transport phenomena in the presence of Weyl corrections have also been discussed in [18][19][20]. Recently, the magnetotransport of a strongly interacting system in 2 + 1 dimensions was examined in a holographic Dirac-Born-Infeld model in [21,22]. Negative magnetoresistance was found for a family of dynoic solutions in [22]. The DC conductivity in the probe DBI case with the vanishing magnetic field was also discussed in [23].
Mott insulators can be parent materials of high T c cuprate superconductors. A Mott insulator has an insulating ground state driven by Coulomb repulsion. Mott-like behavior is that strong interactions between electrons would prevent the charge carriers to efficiently transport charges. Constructing a holographic model describing Mott insulators is still a challenging task. In [24][25][26][27], dynamically generating a Mott gap has been proposed in holographic models by considering fermions with dipole coupling. A holographic construction of the large-N Bose-Hubbard model was presented in [28], and the model admitted Mott insulator ground states in the limit of large Coulomb repulsion. Some other holographic models dual to Mott insulators include [29][30][31]. Recently, a holographic model using a particular type of NLED, namely iDBI, was proposed in [32] to mimic interactions between electrons by self-interactions of the NLED field. It showed that Mott-like behavior appeared for large enough self-interaction strength.
In this paper, we extend the analysis of the magnetotransport in a holographic Dirac-Born-Infeld model in [22] to a general NLED model. As in [22], our analysis is performed in a full backreacted fashion. To break translational symmetry, we follow the method in [33] to add axionic scalars, which depend on the spatial directions linearly.
The rest of this article is organized as follows. In section II, we set up our holographic model. The expressions for the DC conductivities with a finite magnetic field are obtained in section III. Some limiting cases, including high temperature limit, are then discussed.
In section IV, the dependence of the in-plane resistance on the temperature, the charge density and the magnetic field are investigated for Maxwell, Maxwell-Chern-Simons, Born-Infeld, square and logarithmic electrodynamics. In section V, we summarize our results and conclude with a brief discussion.

II. HOLOGRAPHIC SETUP
Consider a 4-dimensional model of gravity coupled to a nonlinear electromagnetic field A a and two axions ψ I with action given by where Λ = − 3 l 2 , and we take 16πG = 1 for simplicity. In the action (1), we assume that the generic NLED Lagrangian is L (s, p), where we build two independent nontrivial scalars using F ab = ∂ a A b − ∂ b A a and none of its derivatives: where, for later convenience, we define θ ≡ L (0,1) (0, 0). Note that we set the AdS radius l = 1 hereafter.
Varying the action (1) with respect to g ab , A a , and ψ I , we find that the equations of motion are where T ab is the energy-momentum tensor: and we define To construct a black brane solution with asymptotic AdS spacetime, we take the following ansatz for the metric, the NLED field and the axions ψ 1 = αx, and ψ 2 = αy, where h denotes the magnitude of the magnetic field. The axions are responsible for the breaking the translational invariance and generating momentum dissipation. The equations of motion then take the form: It can be shown that eqns. (8) and (10) guarantee that eqn. (9) always holds. Eqn. (10) leads to where ρ is a constant. One has f (r h ) = 0 at the horizon r = r h , and the Hawking temperature of the black brane is given by Hence at r = r h , eqn. (8) reduces to where

III. DC CONDUCTIVITY
Via gauge/gravity duality, the black brane solution (7) describes an equilibrium state at finite temperature T , which is given by eqn. (13). The NLED field is a U(1) gauge field and dual to a conserved current J µ in the boundary theory. In this section, we calculate the DC conductivities for J µ using the method developed in [34,35].

A. Derivation of DC Conductivity
To calculate the DC conductivities, we consider the perturbations of the form: where i = x, y and I = 1, 2. The fields a i (r) do not appear explicitly in the NLED Lagrangian L (s, p). Thus, the conjugate momentum of the field a i (r) with respect to rfoliation is radially independent: where the conjugate current is Similarly, the conjugate momentum of the field A t (r) is also a constant flux where one has We can then compute the expectation value of J t for the boundary theory by Using eqns. (6), (7), (11), (15), and (19), we find that, at the linearized order, which means that ρ can be interpreted as the charge density in the dual field theory. Evaluating eqn. (19) at r = r h gives which lead to To express J i in terms of E i , we first consider the constraints of regularity on the metric and fields around the horizon [35]: We then consider the tx and ty component of the perturbed Einstein's equations: Using the regularity conditions (25), eqns. (26) reduce to where s h and p h are given by eqns. (14). Solving eqns. (27) for h ti (r h ) in terms of E i and using the regularity conditions (25) to evaluate eqns. (24) at r = r h , one can relate the currents J i to the electric fields E i via where the DC conductivities σ ij are given by To express σ ij in terms of ρ, h and T , one needs to solve eqns. (13) and (22) for r h and A ′ t (r h ) in terms of ρ, h and T and plug the r h and A ′ t (r h ) expressions into eqns. (29). Therefore, σ ij are in general functions of the temperature T , the charge density ρ, the magnetic field h and the strength of momentum dissipation α. Notice that the conductivities are left invariant under the separate scaling symmetries given by for constant λ. The resistivity matrix is the inverse of the conductivity matrix:

B. Various Limiting Cases
In section IV, we will use eqns. (29) to discuss the properties of the DC conductivities in some NLED models. Before focusing on a specific model, we now consider some limiting cases of the general formulae for σ ij or R ij .

Weak and Strong Dissipation Limits
When α = 0, the system will restore Lorentz invariance. In a Lorentz invariant theory, it showed [37] that the DC conductivities in the presence of a magnetic field were As a check, we find that, in the weak dissipation limit with α 2 ≪ 1, the DC conductivities in eqns. (29) become which are consistent with eqns. (32).
In the strong dissipation limit with α 2 ≫ 1, we find that the DC conductivities become It is noteworthy that eqns. (34) agree with the results in [36], where the DC conductivities were computed for a probe NLED field. In fact, when α 2 ≫ 1, the geometry is almost determined by the contributions from the axionic sector, and hence the NLED field can be approximated as a probe one.

Vanishing Magnetic Field and Charge Density
For the h = 0 case, the DC conductivities reduce to For the iDBI Lagrangian, our results reduce to eqn. (3.1) in [32].
At zero charge density ρ = 0, the DC conductivities become These DC conductivities are in general non-zero and can be interpreted as incoherent contributions [38], known as the charge conjugation symmetric contribution σ ccs . There is another contribution from explicit charge density relaxed by some momentum dissipation, σ diss , which depends on the charge density ρ. Our results show that, for a general NLED model, the DC conductivities usually depend on σ diss and σ ccs in a nontrivial way.

High Temperature Limit
Finally, we consider the high temperature limit T ≫ √ h, √ ρ, α . In this limit, eqn. The blue region indicates insulating behavior.
Green Region indicates negative magneto-resistance.
FIG. 1: Parameter space in terms of θ and h/ρ in the high temperature limit. Note that θ ≡ The resistance then reduces to which only depends on θ = L (0,1) (0, 0) and is independent of the nonlinear effects of the NLED field. This is understood as the nonlinear terms being suppressed by the high temperature. One can define a metal and an insulator for dR xx /dT > 0 and dR xx /dT < 0, respectively. Eqn. (39) shows that, for any NLED model in the high temperature, a metalinsulator transition (MIT) occurs when the term (1 + θ 2 ) h 2 + 2θρh − 1−θ 2 1+θ 2 ρ 2 changes the sign. In FIG. 1(a), we plot the parameter space for a metal and an insulator with respect to FIG. 1(b), where we find • Green Region: In this region, one has that ∂R xx /∂ |h| < 0. To describe how the electrical resistance responds to an externally-applied magnetic field, one can define magneto-resistance as density.
• Yellow Region: In this region, one has that ∂R xx /∂ |ρ| > 0. This is Mott-like behavior, which can be explained by the electronic traffic jam: strong enough e-e interactions prevent the available mobile charge carriers to efficiently transport charges. In particular, when h = 0, eqn. (39) gives that ∂R xx /∂ |ρ| > 0 as long as θ 2 > 1.
If the NLED Lagrangian L (s, p) is CP invariant, one has θ = 0. In this case, R xx becomes which gives that the system displays metallic behavior for |h/ρ| < 1 and insulating behavior for |h/ρ| > 1. Moreover, one always has that ∂R xx /∂ |h| > 0 and ∂R xx /∂ |ρ| < 0. Therefore, there is no negative magneto-resistance or Mott-like behavior for CP invariant NLED models in the high temperature limit. Note that, in [22], eqn. (41) was also obtained for the high temperature limit of the DBI model.

IV. EXAMPLES
In this section, we will use eqns. (13), (22) and (29) to study the dependence of the in-plane resistance R xx on the temperature T , the charge density ρ and the magnetic field h in Maxwell, Maxwell-Chern-Simons, Born-Infeld, square and logarithmic electrodynamics.
The behavior of R xx in the high temperature limit has already been discussed in section III.
So we will focus on the behavior of R xx around T = 0 in this section.

A. Maxwell Electrodynamics
To study the effects of the nonlinear and θ terms on R xx , we first consider Maxwell electrodynamics, in which L (s, p) = s. At T = 0, the resistance R xx is given by which is plotted against ρ/α 2 and h/α 2 in FIG. 2(a). At T /α = 1, we also plot R xx against ρ/α 2 and h/α 2 in FIG. 2(b). Both figures show the saddle surfaces, which imply that We now discuss the dependence of R xx on ρ and h at T = 0. The resistance R xx can be expressed in terms of r h , ρ and h: At zero temperature, the resistance R xx becomes In FIG. 4(a), we plot R xx versus ρ/α 2 and h 2 /α 2 with θ = 3. Similar to Maxwell electrodynamics, we have a saddle surface. However, the valley in FIG. 4(a) is at ρ/α 2 + 3h 2 /α 2 = 0, instead of ρ/α 2 = 0. This twist of the valley would result in the appearance of negative magneto-resistance and Mott-like behavior.
When h = 0, we find that which means that there is no Mott-like behavior for h = 0 if θ 2 ≤ 1. In FIG. 5, we plot the parameter space for ∂ |ρ| R xx > 0 in terms of ρ/α 2 and h/α 2 with θ = 1 and θ = 3. In the yellow region, one has ∂ |ρ| R xx > 0. As expected, the line h = 0 is in the yellow region for θ = 3 while it is not for θ = 1. In FIG. 4(b), we plot R xx versus ρ/α 2 for various values of h/α 2 . One has h/α 2 = 0 for the blue line, and it has a minimum at ρ/α 2 = 0. As one increases the magnitude of the charge density, the value of R xx first increases, then reaches a maximum, and then decreases monotonically. For h/α 2 = 0, the behavior of R xx is different when one moves along the positive and negative ρ/α 2 directions. Along the positive ρ/α 2 direction, the behavior of R xx is similar to the h/α 2 = 0 case. However, if one increases the magnitude of the charge density along the negative ρ/α 2 direction, the value of R xx first decreases until reaching a minimum, then increases until reaching a maximum, and then decreases monotonically.

C. Born-Infeld Electrodynamics
Born-Infeld electrodynamics is described by the Lagrangian density where the coupling parameter a is related to the string tension α ′ as a = (2πα ′ ) 2 . When |a| ≪ 1, we can recover the Maxwell Lagrangian. For the a > 0 case, the properties of the resistance R xx were discussed in [22]. It showed that the behavior of R xx obtained in [22] was quite similar to that in the Maxwell case, which has been investigated in section IV A. On the other hand, the a < 0 case turns out more interesting. For the case with vanishing magnetic filed, the properties of σ xx have been analyzed in depth in [32], in which it showed that the conductivity could decrease with increasing charge density for large enough selfinteraction strength [1] . Here, we extend the analysis to the non-vanishing magnetic field case. We can solve eqn. (11) for A ′ t (r): which shows that there is a singularity at r = r s ≡ 1/4 −a (h 2 + ρ 2 ) for a < 0. To have a physical solution, we need to hide the singularity behind the horizon: r s < r h , which could put an upper bound on h 2 + ρ 2 .
We first discuss behavior of R xx at T = 0. At zero temperature, the condition r s < r h gives that there is an upper bound on h 2 + ρ 2 for a < − 1 6 : We plot R xx versus ρ/α 2 and h 2 /α 2 at T = 0 with a = −0.4 and a = −1 in FIGs. 6(b) and 6(c), respectively. It is noteworthy that, in FIGs. 6(b) and 6(c), the domains of  . 7(a). We find that R xx increases monotonically as one increases the magnitude of the charge density with the magnetic field fixed, which shows that Mott-like behavior occurs in all of the allowed parameter range.
We also display the resistance R xx for a = −1 in FIG. 7(c). When h = 0, R xx increases monotonically with increasing the magnitude of the charge density. For a small but nonvanishing h/α 2 , e.g. h/α 2 = 0.05 and 0.1, the non-monotonic behavior at large values of ρ/α 2 appears. As the value of ρ/α 2 increases, R xx increases first and then decreases after reaching a maximum. However for a larger value of h/α 2 , e.g. h/α 2 = 0.15 and 0.19, we find that R xx decreases monotonically with increasing the magnitude of the charge density, and hence Mott-like behavior disappears. In summary, Mott-like behavior always occurs in When α/ √ ρ = 1, the RHS of the above equation is negative, which explains why the curves in FIG. 8(b) could not go to zero temperature. Since Born-Infeld electrodynamics is CP invariant, one has insulating behavior for h/ρ > 1 and metallic behavior for h/ρ < 1 in the high temperature limit, which is clearly shown in FIGs. 8(b) and 8(d). However, at low temperatures, the temperature dependence of R xx in the a < 0 case is quite different from those in the a > 0 and Maxwell cases. For α/ √ ρ = 1, FIG. 8(b) displays insulating behavior at low temperatures. As one increases the temperature, the system would start to exhibit metallic behavior. If one keeps increasing the temperature, the system would stay metallic behavior for h/ρ < 1, but it would return to insulating behavior for h/ρ < 1. In the α/ √ ρ = 10 case, according to FIG. 8(d), one has insulating behavior for h/ρ < 1 and metallic behavior for h/ρ > 1 at low temperatures. Therefore, increasing the magnitude of the magnetic field would induce a transition or crossover from insulating to metallic behavior at low temperatures and that from metallic to insulating behavior at high temperatures.
We find that the system does not exhibit negative magneto-resistance or Mott-like behavior at high temperatures. If one has negative magneto-resistance or Mott-like behavior at low temperatures, they would disappear at a high enough temperature. We plot R xx as a function of h/ρ for various values of T / √ ρ with α/ √ ρ = 10 in FIG. 9(a). With a fixed value of the charge density, FIG. 9(a) shows that R xx decreases with increasing h for T / √ ρ = 0.25, 0.5 and 1, and R xx increases with increasing h for T / √ ρ = 2 and 3. Similarly, FIG. 9(b) shows that, with a fixed value of the magnetic field, the system displays Mott-like behavior for T / √ h = 0.25, 0.5 and 1, and R xx decreases with increasing ρ for T / √ h = 2 and 3.

D. Square Electrodynamics
Consider a Born-Infeld like Lagrangian which gives the same result for σ ij /R ij as Born-Infeld electrodynamics in the h = 0 case.
We now study the dependence of R xx on ρ and h at T = 0. Since the behavior of R xx in the a > 0 case is similar to that in Maxwell electrodynamics, we focus on the a < 0 case.
In FIG. 10(a), we plot R xx as a function of ρ/α 2 and h 2 /α 2 with a = −0.4, which displays negative magneto-resistance at fixed charge density for all of the allowed parameter range.

E. Logarithmic Electrodynamics
Finally, we consider logarithmic electrodynamics, whose Lagrangian is described by We display the dependence of R xx on ρ/α 2 and h/α 2 in the the a = −1 case at T = 0 in FIG.   11(a), which shows that ∂R xx /∂ |h| < 0. However , FIG 11(b) shows that ∂R xx /∂ |ρ| > 0 for small values of h/α 2 , and ∂R xx /∂ |ρ| < 0 for large enough values of h/α 2 , which means a strong enough magnetic field would destroy Mott-like behavior.

V. DISCUSSION AND CONCLUSION
In this paper, we used gauge/gravity duality to investigate the properties of the DC conductivities with a finite magnetic field of a strongly correlated system in 2+1 dimensions.
The charge current in the boundary field theory is dual to a NLED field in bulk. In our holographic setup, we considered the backreaction effects of the NLED field on the geometry and introduced axionic scalars to generate momentum dissipation. We then presented the expressions for the DC conductivities for a general NLED field. Specifically, one can use eqns. (13), (22) and (29) to express the DC conductivities in terms of the temperature T , the charge density ρ and the magnetic field h of the dual field theory. In the second part of our paper, we discussed the properties of the in-plane resistance R xx in some interesting NLED models, where there appeared Mott-like behavior or negative magneto-resistance in some cases. In Table I, we summarize the results for the ρ and h dependences of R xx at zero temperature in the NLED models discussed above. (55) For the NLED models discussed in section IV, we have L (2,0) (0, 0) = a. For a > 0, the See FIG. 2(a). ∂ |h| R xx > 0. See FIG. 2(a).  See FIG. 6(a). See FIG. 6(a). See FIG. 7(a). ∂ |h| R xx < 0. See FIG. 7(b).

for small values of
h/α 2 and ρ/α 2 . See FIG. 7(c). See FIG. 7(d).  non-linearities correction tends to reduce the strength of the repulsive force between two electrons. However for a < 0, the correction tends to increase the strength of the force.
So it is natural to expect that a negative a may correspond to strong interactions between electrons in our holographic model, which could lead to Mott-like behavior [1] .
It is interesting to note that Born-Infeld electrodynamics with a = −1 could describe Mott [1] In [32], the role of a has been discussed for iDBI model.
insulator to metal transition (IMT) induced by a magnetic field. In fact, FIG. 8(d) shows that, at low temperatures, the system has insulating/metallic behavior for a weak/strong magnetic field. On the other hand, FIG. 7(c) displays that, for a weak magnetic field at low temperatures, the system usually has Mott-like behavior and hence is a Mott-Insulator.
When the magnetic field grows strong enough, Mott-like behavior disappears, and meanwhile, the system exhibits metallic behavior. A magnetic field-induced IMT for a Mott system, namely a bilayer ruthenate, Ti-doped Ca 3 Ru 2 O 7 , was presented in [39]. Our analysis for IMT is rather qualitative, and it deserves future more detailed studies.
For a negative enough a, we found that, at low temperatures, the resistance R xx always decreases with increasing magnetic field, which appears as negative magneto-resistance. On the other hand, the behavior of R xx in the a > 0 NLED models is similar to that in Maxwell electrodynamics, in which one always have positive-resistance. It seems that, decreasing a from a positive value to a negative one, which corresponds to increasing the strength of the interactions between electrons, would lead the magneto-resistance to change from positive to negative. It showed in [40] that the magneto-resistance could change from positive to negative by gradually introducing artificial disorder through Ga + ion irradiation to pristine graphene. To relate our results to the experiments, we need to better understand how a is related to external control parameters.
In this paper, we found the σ ij /R ij expressions for a general NLED field. Our analyses for the properties of the resistance in NLED models are preliminary. One can use these expressions to find or construct a NLED model to realize some interesting experimental results, such as the scaling relationship between applied magnetic field and temperature observed in the magneto-resistance of the pnictide superconductor.