Fermionic response in nonlinear arcsin electrodynamics

We consider certain blackhole solution in non-linear arcsin electrodynamics coupled with gravity and axions. We have studied the behaviour of the fermionic operators in the dual (2+1)-dimensional theory. We consider holographic spectral function for both the backreacted solutions and probe limit over the range of physical parameters. We find that with a variation of the charge density the system changes from Fermi liquid to non-Fermi liquid and the transition point depends on the temperature.


Introduction
Strongly coupled systems arose in condensed matter turns out to be suitable arena [1][2][3][4][5][6][7][8] for application of holographic techniques of gauge/Gravity duality [9][10][11][12]. These techniques can address strongly coupled non-gravitational systems in the set up of a gravity theory with a black hole background. In particular, analysis of fermionic excitations [13][14][15] exhibits scaling behaviour of non-Fermi liquids in this set up. Further studies regarding scaling exponents and dimension of dual theories appear in [16]. Introduction of dipole coupling on the gravity side leads to dynamically generated gap [17,18]. Charged Lifshitz black branes were considered in [19,20]. Effect of doping parameters on fermionic excitations in the holographic set up was studied in [21,22] and gives rise to transition between Fermi liquid to non-Fermi liquid phase. Similar studies of fermions in top down approach appear in [23][24][25][26][27][28].
All these studies involve Maxwell electrodynamics in the gravity setup. Even though these linear electrodynamic models are successful in explaining various features in the context of condensed matter systems, it fails to capture the behaviour of some of the systems, one such is cuprates high T c superconductors in the strange metal phase where the resistivity depends linearly on temperature and so called the anomaa e-mail: nishalrai10@gmail.com lous scaling behaviour of cuprates. Perhaps that triggers further explorations of the non-linear electrodynamic models. A string inspired Dirac-Born-Infeld (DBI) model was analysed in [29] to obtain DC conductivity along the method purposed in [30,31]. This was further generalised in [32] and [33] leading to linear scaling of resistivity with a temperature. More general non-linear electrodynamics was considered [34,35].
Considering the success of non-linear electrodynamic models it is natural to study the Fermi surface and excitations around it. On that score, [36] has studied the fermionic behaviour for string inspired DBI model. In the present work, we will consider another non-linear model of electrodynamics known as arcsin electrodynamics, which was introduced in [37]. An attractive feature of this model is the finite electric field of a pointlike charge at the origin and static electric energy of a particle is also finite [38,39]. The dual boundary theory of this model has been analysed in several works. By introducing a scalar field in the gravity theory, it was shown [39], that the boundary theory may undergo condensation, that corresponds to spontaneous breaking of a U (1) symmetry and thus it admits a superconducting phase [39]. This model was further extended by introducing a neutral scalar field and axions, where the backreaction on the metric was also considered. This extended model exhibits metal/insulator transition with a variation of certain parameter of the theory, showing additional phase structure. Thus among the possible non-linear electrodynamics, it has been established that this particular theory has a wide spectrum of phases. In particular, as has been shown in [35], this model shows a linear temperature dependence of resistivity in a probe limit and thus it can capture the properties of direct conductivity of cuprates in the strange metal phase. However, analysis of Fermi surface and excitations around it has not been done yet and the aim of the present work is to fill in this gap. For this model, we have studied the nature of fermionic behaviour over the range of several parameters. We find there is a transition/crossover from Fermi liquid to non-Fermi liquid phase occurring due to variation of charge density of the system. We have also studied the fermionic excitations in this probe limit as well, which may give some insight on the underlying mechanism for such anomalous behaviour of cuprates.
The plan of the article is as follows. In the next Sect. 2, we will briefly introduce the nonlinear arcsin electrodynamics model and its blackhole solution. We will also discuss the probe limit and its background solution. In Sect. 3, we introduce Fermions and compute the Green's function of the dual operator. In Sect. 4, we present our results where we numerically solve the Dirac equation to obtain the Green's function and study its nature with the variation of different parameters. In Sect. 5 we conclude with a discussion.

Bosonic part
In this section, we revisit the model of non-linear arcsin electrodynamics which is interacting with gravity introduced in [35]. The action is given by where ψ (I ) are the two axions and φ is the scalar field whose potential is given by V (φ). S 0 is the action for non-linear electrodynamics where F = 1 4 F μν F μν . Z 1 (φ) and Z 2 (φ) are the coupling constant which depends on the scalar φ. One may observe that the Eq. (2.2) reduces to the usual Maxwell action at the limit of vanishing gauge field A μ → 0 with the choice of Z 2 = 1 and Z 1 = 1.
The equation of motion for the above action turns out to be as follows. From the Einstein part of the action we have where θ μν can be expressed in terms of energy momentum tensor of gauge fields T μν as μν = T μν − 1 2 g μν T λ λ and in the second term on right hand side sum over I is implied.
The equation for scalar and axion fields φ and ψ I are given by The equation for gauge field turns out to be The quantity θ μν , which is related to the energy-momentum tensor of the gauge field is given by The expression of the current J μ in the dual field theory is which is to be evaluated at the boundary. In order to get the black hole solution, we will consider the following ansatz for the metric, gauge field, scalar field and axions where h is the magnetic field which lies on x − y plane. The scalar field φ only depends on radial direction r . The axions ψ I is set up in x and y direction with the magnitude of momentum dissipation k p which breaks the translational symmetry at the boundary. Since we are only interested on the electrically charged solution of the black hole, we are setting h to be zero. The equations are quite involved, so we will consider the particular family of the solution by following the ansatz given in [32,35] (2.9) z 1 and V 0 are constants where we consider z 1 = 1 for further simplification. After taking these assumptions, the above equation of motion can be solved exactly. The solution for the metric coefficients are given in [35], where C(r ) is solved exactly and is given by (2.10) Substituting C(r ) in the above equation of motion for scalar, the equation for the metric coefficient D(r ) becomes Solving the above first order differential equation and the equation of motion for the gauge field the coefficients D(r ) of the metric and the gauge field are given by where ρ is the charge density of the system and M mass of the blackhole. M is fixed by D(r ) which vanishes at the horizon (D(r = r h )=0). The temperature of the blackhole is given by where we have defined = 1 √ 1−F 2 . In the expression of temperature, F is evaluated at horizon r = r h .

Probe geometry
As has been pointed out that strong momentum dissipation limit may be considered as a probe limit [35], where one can ignore the backreaction on the metric. In this probe limit we will consider hyperscaling violating geometry as the background, which may reproduce linear temperature dependence of resistivity [34]. The hyperscaling geometry appears as a solution [34] with the coupling given by V (φ) ∼ −V 0 e ηφ and Y (φ) ∼ e αφ , where α and η are two constants. The ansatz for the metric and other fields for hyperscaling violating geometry are given by (2.14) where κ, k p and L are constants. z and θ are the Lifshitz scaling and hyperscaling violation exponents respectively. The solution with given ansatz has been studied in [35], where the metric coefficient f (r ) and other parameters are Here, r h is the radius of horizon satisfying f (r h ) = 0. The gauge field can be easily obtained by solving the equation of motion given in [26]. Here, we will only give the expression for the gauge field for this background which turns out to be where we have definedC = (Z 1 Z 2 r 2(1−θ) ) 2 with Z 1 ∼ e γ φ and Z 2 ∼ e δφ . The temperature is given by We would like to comment on the parameters z and θ before proceeding. z = 1 and θ = 0 corresponds to AdS background. To have a well-defined geometry and resolvable singularity, the range of these two parameters is restricted. This restriction arises by considering the Gubser criterion in conjunction with the null energy condition In addition to this, in order to have linear temperature dependence of resistivity of cuprates one need to satisfy (2.20) provided λ ≥ 0, for simplicity we will be taking λ = 0. Even though this behaviour arises in high temperature limit, one can choose the appropriate value of the parameters to retain the behaviour of cuprate at sufficiently low temperature.

Near Horizon limit
The nature of the fermionic excitation is closely related to the near horizon structure of the metric. In this regard, we will analyse near horizon structure of the metric given in Eqs. (2.8) and (2.14). To study the fermionic spectral function we will consider two different backgrounds; First, we will consider fully backreacted background which is given in (2.8) and second we will consider the hyperscaling violating geometry which arises in the probe limit. We will begin with the discussion of the former with extremal case (T = 0).

Backreacted geometry
The full metric for fully backreacted geometry is given in (2.8), we would like to rewrite this metric in the following form vanishes, similarly we will define a point r * such that the derivative of f (r ) at r * vanishes ( f (r * ) = 0). However, in the limit of zero temperature, these two points coincide with each other i.e. r * = r h . First, we begin the near horizon expansion of the metric in this limit of zero temperature. The metric coefficient f (r ) at the near horizon (r → r h ) develops a double zero which is given by where all the quantities in L 2 are evaluated at (r = r h ). A natural scaling factor L 2 of the metric occurs in the near horizon limit. We have used Eq. (2.11) to find D(r ) and L 2 .
Following [16] we will consider the scaling limit (2.23) In this limit the metric (2.8) and the gauge field (2.12) are given by where e = √ −2F L 2 2 evaluated at horizon (r = r h ). We can see from the above Eq. (2.24), that the structure of the metric (2.8) at near horizon becomes Ad S 2 × R 2 , with L 2 being the radius of curvature of Ad S 2 .
To generalise this in the case of finite temperature, where r * = r h one considers the following equation in addition to (2.23) [16,22] the metric under the additional equation turns out to be blackhole in Ad S 2 × R 2 given by where the temperature is defined as T = 1 2πζ 0 . It may be interesting to check the asymptotic limit of this geometry which turns out to be Ad S 4 given by where we have defined b 1 as with c 2 1 = −1 + 1 + ρ 4 .

Probe geometry
In the hyperscaling violating geometry, one can consider the similar near horizon expansion (r → r h ). Before proceeding, we will consider the following general coordinate transformation r = u 1 z−θ , under which the metric (2.14) turns out to be In the extremal limit, the metric at near horizon turns out to be Ad S 2 × R 2 with the metric and gauge field given by where e = −2F and a 1 = F (u) are both evaluated at u * . 1 To deduce the above near horizon metric, we have chosen the following scaling limits In the finite temperature, one can consider the additional scaling limit given in (2.25) and follow the similar procedure. On doing so, the metric turns out to be a blackhole in Ad S 2 ×R 2 , with metric and the gauge field given by where L 2 and e are given in (2.30).

Green's function
In order to probe the system, we will consider Dirac action in fully backreacted geometry and hyperscaling violating geometry as well. The action is given by where the first term corresponds to the kinetic term, the second term is the mass of the spin 1/2 particle and the last term corresponds to the coupling of Dirac field with the electrodynamic field, where the strength of the coupling is given by the dipole coupling p. The covariant derivative and Gamma matrices are given by,

2)
(e μ ) a is the vielbein, ω μν is the spin connection and q is the charge of the spin 1/2 particle. The equation of motion for the above action is given by, To eliminate the spin connection in the Dirac Eq. (3.3), we will consider the following ansatz ψ = (−gg rr ) 1 4 e −iωt+ik i x i λ(r ). For further simplification, we will now introduce the projection operator [24,26] 1 u * =u h in the case of extremal limit (T = 0). and rewrite λ(r ) as λ ± = ± λ(r ), wherek is the unit vector along the spatial direction. Our choices of Gamma matrices are We can write λ ± in terms of two component spinor given . In addition, we use the spatial rotational symmetry to set k i = k for i = x and zero for i = y. With the above choices together with Gamma matrices the Dirac equation reduces to Before discussing solution of this Dirac Eq. (3.5) let us consider the near horizon limit of it. We will consider only the backreacted geometry for this purpose. For the extremal case (T = 0), using (2.24), Dirac equation simplifies into The solution to this equation is given by λ ± = a E + ζ −ν + bE − ζ ν , where E ± are the real eigen vector of U and the exponent is the eigen value of U given by For the finite temperature, calculation remains similar except one has to consider the near horizon metric and gauge field for the finite temperature given in Eq. (2.26). In this case, the eigenvalue ν is given by ν plays a central role in determining the nature of the Fermi surface. The lifetime of the fermionic excitations around the Fermi surface is given by [22] τ = ω −ν at zero temperature and τ ∼ T −2ν at finite temperature. As the lifetime of the excitation depends on ν, the metallic behaviour of the system is controlled by ν.
In the region of phase diagram with ν ≥ 1, we have a normal metal phase, where the excitation of the quasiparticles is governed by Landau Fermi liquid. However, in the region where 1/2 < ν < 1 we get stable quasiparticles, whose lifetime scales differently from that of Landau Fermi liquid theory. In the region with ν < 1/2, we have a short lived quasiparticles which are the characteristic of strange metal phase, this region is called non-Fermi liquid region. In between the regime of stable and non-stable quasiparticle, there is a transition point where ν = 1/2, this region corresponds to a marginal Fermi liquids regime. The conformal dimension of the dual operator in the IR CFT will be = ν + 1 2 . There is a range of momentum k, known as oscillatory region [24,26], for which ν becomes complex thus leading to complex dimension of the dual operator.
Coming back to the solution of the Dirac Eq. (3.5) for the full metric, we need a boundary condition. For that purpose we define η± = λ 1± λ 2± , which yield the flow equation given by where v ± = g rr g x x (ω + q A t ) ± p∂ r A t . (3.10) To solve the above Eq. (3.9) numerically, we impose the infalling boundary condition given by The boundary retarded Green's function is given by In the background with hyperscaling violating geometry in the UV limit, the boundary Green's function in the regime of linear response has been derived in [40] and is given by Comparing (3.12) and (3.13), one may notice that for the massless case, the two Green's function matches. We can see from the Eq. (3.9) that the two diagonal part of the retarded Green function are related to each other by flipping the sign of k i.e. G 11 (ω, k) = G 22 (ω, −k). So it is sufficient to evaluate the Green's function for only one from the two, hence we will be only considering G 22 for our calculation.
In the next section, we will be studying the behaviour of the spectral function given by A(ω, k) = I m(G 22 ) with respect to variation in charge density ρ, charge of the particle q and the momentum dissipation term k p in zero and finite temperature for fully backreacted background. In case of the hyperscaling violating geometry, we will study the fermionic spectrum for a specific choice of background parameter which leads to linear temperature response of resistivity.

Result
In this section, we consider the behaviour of the Fermi surface associated with the operators dual to the fermionic modes. As mentioned earlier, we will limit ourselves to G 22 . Considering G 22 means that we only consider η − , which is sufficient, as the behaviour of the other Fermion for (η + ) will be similar i.e. if we find a Fermi surface for η − at k F than the Fermi surface for η + will be in −k F . For the sake of simplicity, we will be considering the massless case with m = 0. We will also set the value of the scalar potential to be V 0 = 6 and Pauli's coupling p = 0. To begin, we numerically solve (3.9) where we impose the infalling boundary condition given in Eq. (3.11), and obtain the spectral function A(ω, k) as function of ω and k. In order to locate the Fermi surface, we will look for the poles of the Green's function in case of the zero temperature and in case of finite temperature we will follow [24,28], where the position of Fermi surface is given by sharp peak in the plot of a spectral function around k = k F at ω = 0, which has a small enough width. In addition, if we plot the spectral function vs. ω at k = k F that should show a peak around ω= 0.

Backreacted background
In this subsection, we have studied the behaviour of the fermionic spectrum with variation of the parameters for the fully backreacted background.
From (2.13) fixing the value of the temperature T , we get a relation between the charge density ρ and the momentum dissipation term k p and on choosing the value of one determines the value of the other. Here, the value of k p is determined by choosing the value of ρ and vice versa.
First, we begin with the extremal case T = 0, where we vary the charge density ρ and look for the poles of the Green's function for two cases with charge q = 2 and q = 1. We have plotted k F versus the inverse of charge density (1/ρ) for q = 2 in Fig. 1a and q = 1 in Fig. 2a. Blue dots correspond to Fermi momenta (k F ), which corresponds to poles of the Green's function. In the same figure, the shaded portion which is tapering toward the right represents an oscil-    Fig. 2a, with the decrease of charge density of the system Fermi surface, moves toward the oscillatory region and enters around ρ = 1.333. In order to study the change in nature of the fermionic excitations further, we have plotted ν versus inverse of charge density (1/ρ) of the system in Fig 1b for q = 1 and Fig 2b for  q = 2. As we can observe that on decreasing the charge density (ρ), ν starts decreasing and vanishes around ρ = 0.8849. Initially, at sufficiently large charge density, the system was in a Fermi regime, as one start decreasing the charge density the system starts to move from the regime of Fermi liquid enters the marginal Fermi regime at ρ = 1.1628 where ν = 1/2 and reaches the non-Fermi regime. Further decrease of ρ leads to vanishing of ν signalling non existence of the Fermi surface, where Fermi momenta are inside the oscillatory region. Thus with the decrease of charge density ρ, there is a transition from Fermi to the non-Fermi regime. This transition/crossover is taking place at zero temperature and may be related to a quantum phase transition. In addition, we can see from the comparison of the Figs. 1b and 2b, that with the increase of charge of the particle q the point of transition from Fermi to non-Fermi get shifted to the right of the parameter space.
In order to make the change of behaviour with momentum dissipation more explicit, we have plotted ν vs.
√ ρ/k p given in Fig. 3a  ν = 1/2, the upper plane above ν = 1/2 is the Fermi regime and lower plane correspond to the non-Fermi regime. There is a critical value of the momentum dissipation at which ν = 1/2. As we move toward the upper plane, the quasiparticles become stable and finally enters the region of Landau Fermi liquid which is given by the region with and above the green dashed line in Fig. 3a, b. Comparing Fig. 3a, b we see for higher charge, q = 2 the transition is occurring at a higher value of k p . The transition is driven by the changing the magnitude k p of axionic scalars and can be thought of as a disorder driven transition [32].
After discussing the zero temperature case we will move to the finite temperature. We have plotted spectral function A(ω = 0, k) vs k for four different temperatures starting from T / √ ρ = 0.035. We have set the value of charge density ρ = 2 and q = 2. We can see that with the increase in temperature the height of peak starts decreasing and becomes broader and move towards the left. With the increase in temperature, the density gets spread out due to thermal excitation leading to a broader peak (Fig. 4).
In order to see the change in the nature of fermionic excitations over the range of momentum dissipation with variation of temperature, we have given a density plot in T / √ ρ vs. √ ρ/k p plane in Fig. 5a, b, with a colour gradient representing values of ν, for two different values of charge q = 1 and q = 2. The value of ν increases from left to right, with darker and lighter shades corresponding to lower and higher values of ν respectively. The white region on the left is the oscillatory region and the green dashed line is the curve for ν = 1/2, where on the left of the curve we have strange metal phase and on the right, we have stable quasiparticle phase. The density plot extends for the higher value of ν > 1 which is given by the yellow dashed line representing a normal metal phase. We can see from the figures that, with an increase in temperature the critical value of momentum dissipation k p for a transition from the strange metal to the normal metal phase decreases. Even though for a given temperature the transition occurs on the specific value of momentum dissipation for different values of charge of the particle q, the qualitative behaviour does not depend on this variable q. This is similar to the normal metal to strange metal transition found in [22] where the transition occurred due to tuning of doping instead of a variation of the magnitude of axionic scalar.
In order to see the similar changes with charge density, in Fig. 6a we have a density plot in T / √ ρ vs. ρ plane, where we can see the critical transition value of ρ from strange √ ρ/k p plane for q = 1.     metal to normal metal phase increases with the increase in temperature. It may be worth mentioning that, as the system moved from the Fermi to the non-Fermi region, the height of the peak of the spectral function decreased exponentially. This can be observed from Fig. 6b, where we have plotted height of the spectral function (A(ω = 0, k = k p )) vs. ρ for q = 1 at T /ρ = 0.025, with the blue dots corresponding to the Fermi momenta. In order to see the role played by the arcsin parameter Z 2 in the fermionic behaviour, in addition to Z 2 = 1 we have also considered Z 2 = 0.25 and 0.5, where we have chosen T = 0, q = 2 and z 1 = 1. In Fig. 7a we have plotted 1/ρ vs. ν for different values of Z 2 , with blue, brown and purple curve showing Z 2 = 1, 0.5 and 0.25 respectively, with red and green dashed line corresponding to ν = 1/2 and ν = 1 line respectively. We can see from this Fig. 7a that with an increase of Z 2 the transition point from Fermi to Non-Fermi liquid regime changes, for higher value of Z 2 the transition occurs at higher value of charge density ρ. From the Fig.  7a or more explicitly from Table 1, one may also observe that for a given value of ρ there is a transition from Fermi to Non-Fermi liquid with an increase of arcsin parameter Z 2 . We have also plotted Fermi momenta k F vs. inverse of charge density 1/ρ in Fig. 7b, with blue, brown and purple curve representing Z 2 = 1, 0.5 and 0.25 respectively, with black dashed line showing ρ = 0.91. From this Fig. 7b we can see that for a given value of ρ the Fermi momenta k F decreases with the increase of arcsin parameter Z 2 as given in Table 1. A similar kind of behaviour has been observed in [36], where the Fermi momenta increases with the decrease of BI parameter and the system which is already in the non-Fermi liquid regime becomes "more non Fermi".

Probe geometry
In this background, we will follow a strategy similar to the one mentioned in the previous section of finite temperature. We will set p = 0 and study the massless fermionic spectrum in finite temperature with q = 2, where the temperature is controlled by the position of the horizon (r h ). In this background, we will move our UV to r → 0 by considering the transformation r → 1/r . In this regards, we will consider the condition given in (2.19) z > 2 and θ < 2. We will set λ = 0, α = 1 and z = 2.2 which in turn determines the value of other parameters for cuprates like scaling given in (2.20).
We have plotted spectral function A(ω = 0, k) vs k for ρ = 1 in Fig. 8a where we find a very sharp peak at k F = −6.305, which is confirmed as the existence of Fermi surface from A(ω, k F ) vs k plot given in Fig. 8b where we again see a very sharp peak at ω = 0 with k = k F . Our analysis of probe limit at finite temperature shows existence of distinct Fermi surface. On comparison with the result of fully backreacted geometry in finite temperature, it can be observed that the peak is much sharper for the probe limit, which may be due to fact that the hyperscaling violating geometry leads to sharper peak. On the other hand, comparing the result with [40] where they have considered a linear Maxwell electrodynamics with the same background of hyperscaling violating geometry as used in our probe limit, we find that the present analysis leads to much sharper peak.
. ω for q = 2, ρ = 1. In this case, one can attribute sharpness of the peak to the non-linear electrodynamics considered in this work.

Discussion
We have studied the fermionic behaviour of the nonlinear arcsin electrodynamics model given in [35]. We have evaluated the spectral function to determine the existence of a Fermi surface in zero and finite temperature. In case of zero temperature, we find that for a given value of charge density ρ the particle with higher charge q has a higher chance of forming a Fermi surface. For a given charge of the particle q, the system with large charge density ρ is in the region of normal metallic phase and on decreasing the charge density the system starts to move toward the marginal Fermi liquid regime where ν = 1/2. On the further reduction of the charge density ρ results in the transition/crossover toward the strange metal phase. The system shows a similar kind of behaviour with the variation of the momentum dissipation k p , where there is a transition/crossover from normal metal to strange metal phase, which may be considered as a disorder driven transition [32]. Next, we investigated the finite temperature case, where we see for a given charge q with the increase in temperature the height of the spectral function A(ω, k) starts decreasing and becomes broader and broader, showing that as the temperature rises the density gets spread out due to the thermal excitation. We have also studied the nature of the Fermi surface with the variation of momentum dissipation k p and charge density of the system ρ for finite temperature case. We see that with an increase in temperature the critical value of momentum dissipation for the transition decreases. In case of variation with the charge density ρ, with increase in temperature the critical value of charge density for the transition increases. For a given value of momentum dissipation and charge density with the increase in temperature, there is a transition from normal metal to strange metal phase, which is similar to the case in cuprates high T c superconductors where the transition depends on the doping.
In the present model as already mentioned in the introduction, we have considered a particular phase of the dual boundary theory. Comparing these results with the results of [39], we see that the boundary charge density ρ plays an important role in both the phases of the dual boundary theory. In the case of [39], ρ determines the critical phase transition temperature of the superconductor and also in the present case ρ determines the critical transition temperature from the regime of Fermi to Non-Fermi. In both cases, the transition temperature increases with the increase of charge density ρ. In particular, we have observed that in [39] for a given temperature as β increases the critical value of charge density increases similarly in the present model if we keep the temperature fixed at T = 0 and increase the analogous parameter Z 2 the critical value of charge density also increases. We have also studied the fermionic behaviour with the variation of the arcsin parameter Z 2 , where we find that with the increase of the arcsin parameter there is a transition from Fermi to non-Fermi regime. We have also analysed the behaviour of Fermi momenta k F with the variation of the arcsin parameter Z 2 , where we find that for a given boundary charge density ρ, Fermi momenta decreases with the increase of arcsin parameter. Comparing this result with another model of nonlinear electrodynamics studied in [36], where the nonlinearity was controlled by BI parameter, we see that in both the cases with the increase of nonlinear parameter the Fermi momenta decreases and the system which is already in the non-Fermi regime moves toward "more non-Fermi" regime.
We have also considered the probe limit of this nonlinear system with hyperscaling violating geometry as the back-ground. It was shown in [33,35] that with the particular choices of the parameters in this background, the system exhibits a linear temperature dependence of resistivity as in the case of cuprates. We have chosen those particular parameters and studied the fermionic spectrum in a finite temperature, where we find a very sharp peak in the spectral function indicating the existence of the Fermi surface. Comparing these results with [40], where linear Maxwell electrodynamics was considered in the same background as our probe background, we find that the height of the peak in the our case is much higher and sharper. This feature can be attributed to non-linear electrodynamics.
In this work, we have chosen the neutral scalar field φ to vanish in the case of backreacted geometry and consider a simpler solution. Considering more general configurations with non-zero neutral scalar φ along with non-trivial choices of the coupling Z 1 may lead to further interesting solution and it would be interesting to study the fermionic behaviour for such backgrounds. It has been shown in [35] that by turning on the magnetic field h may give rise to metal-insulator transition with a variation of h. Analysis of fermionic excitations in presence of magnetic field may lead to further insight on such transitions. In [39], the author have analysed the superconducting phase of this nonlinear electrodynamic model in finite temperature. It may be interesting to study the fermionic behaviour in this superconducting phase at finite temperature as well as its zero temperature limit which may correspond to a domain wall solutions along the line of [41][42][43]. As we have mentioned the dual boundary theory associated with this model admits a rich phase structure and it would be interesting to incorporate all of these phases in a single gravity theory. Another natural extension of the present work is to consider U (1) × U (1) gauge fields in this setup and introduce doping along the line of [22]. We would like to report on some of these in future.