Axionic charged black branes with arbitrary scalar nonminimal coupling

In this paper, we construct four-dimensional charged black branes of a nonminimally coupled and self-interacting scalar field. In addition to the scalar and Maxwell fields, the model involves two axionic fields homogeneously distributed along the two-dimensional planar base manifold providing in turn a simple mechanism of momentum dissipation. Interestingly enough, the horizon of the solution can be located at two different positions depending on the sign of the parameter associated to the axionic field, and in both cases there exists a wide range of values of the nonminimal coupling parameter yielding to physical acceptable solutions. For a negative parameter that sustains the axionic fields, the allowed nonminimal coupling parameters take discrete values and the solution turns out to be extremal since its has zero temperature. A complete analysis of the thermodynamic features of the solutions is also carried out. Finally, thanks to the mechanism of momentum dissipation, the holographic DC conductivities of the solutions are computed in terms of the black hole horizon data, and we analyze the effects of the nonminimal coupling parameter on these conductivities. For example, in the purely electric case, we notice that as long as the nonminimal coupling parameter takes the discrete values associated to the extremal solution, the DC conductivity vanishes identically reproducing in turn an insulator behavior. In the non extremal case, we point out the existence of a particular value of the nonminimal coupling parameter (which is greater than the conformal one in four dimensions) yielding an infinite conductivity; this is due to the fact that the translation invariance is restored at this point. Finally, in the dyonic case, we show that the conductivity matrix for the extremal solution has a Hall effect-like behavior.


Introduction
In the last decades, the ideas underlying the anti de-Sitter/conformal field theory (AdS/CFT) correspondence have been applied to get a better understanding of phenomena that occur in the condensed matter physics as the quantum Hall effect, the superconductivity or the superfluidity, see e.g. [1,2]. As a significative example, we can mention the case of charged hairy black holes with a planar horizon that may be relevant to describe the behavior of unconventional superconductors [3]. In this scenario, the nonzero condensate behavior of the unconventional superconductors is mimicked by the existence of a hair at low temperature that must disappear in the high temperature regime [3]. Nevertheless, finding black holes with such features is a highly nontrivial problem that is rendered even more difficult by the various no-hair theorems with scalar fields existing in the current literature, see e.g. [4]. Fortunately, the precursor works of Refs. [5,6] have established that scalar fields nonminimally coupled seem to be an excellent laboratory in order to escape the standard scalar no-hair theorems. Indeed, as shown independently in Refs. [5,6], conformal scalar field nonminimally coupled to Einstein gravity can support black hole configuration with a nontrivial scalar field. These black hole solutions have been dubbed BBMB solution in the current literature. However, the BBMB solution suffers from some pathology essentially because of the divergence of the scalar field at the event horizon. This inconvenient makes its physical interpretation and the problem of its stability a subject of debate. A way of circumventing this pathology consists in introducing a cosmological constant whose effect is to precisely push the singularity behind the horizon, and as a direct consequence, the scalar field becomes well-defined at the event horizon [7]. It is important to stress that for the BBMB solution or its extensions with cosmological constant dubbed as the MTZ solution [7,8], the parameter ξ that couples nonminimally the scalar field to the curvature is always the conformal one in four dimensions, namely ξ = 1 6 , and the horizon topology of the black hole solutions is either spherical or hyperbolic depending on the sign of the cosmological constant. There also exist examples of black hole solutions with the conformal coupling ξ = 1 6 but with a potential term that breaks the conformal invariance of the matter source, see e.g. [9,10]. Nevertheless, charged black holes with planar horizon topology for a scalar field nonminimally coupled to Einstein gravity with or without a cosmological constant are not known. 1 It seems that extra matter source is needed in order to sustain planar charged black hole with a nonminimal scalar field. This intuition is based on the works done in [14,15] where a planar version of the MTZ solution was rendered possible thanks to the introduction of two 3-forms that were originated from two Kalb-Raimond potentials. Interestingly enough, this construction was also extended for arbitrary nonminimal coupling in [16]. Very recently, it has also been shown that (charged) planar AdS black holes can arise as solutions of General Relativity with a source composed by a conformal scalar field together with two axionic fields depending linearly on the coordinates of the planar base manifold [17]. The existence of these planar black holes is mainly due to the presence of the axionic scalar fields which, in addition of inducing an extra scale, allow the planar solutions to develop an event horizon. More precisely, as proved in [17], the black hole mass is related to the parameter associated to the axionic fields, and hence these axionic (charged) black branes can be interpreted as extremals in the sense that all their Noetherian charges are fixed in term term of the axionic intensity parameter. Many other interesting features are inherent to the presence of axionic fields for planar solutions. Among others, axionic fields depending linearly on the coordinates of the planar base manifold provide a very simple mechanism of momentum dissipation [18]. From an holographic point of view, this feature has a certain interest since, as established in Refs. [19,20], the computation of the DC conductivities can be uniquely expressed in term of the black hole horizon data. Mainly because of these results, the study of axionic black hole configurations in different contexts has considerably grow up the last time, see e.g. [21][22][23][24][25][26][27][28].
In the present, we plan to extend the work done in the conformal situation [17] to the case of a four-dimensional scalar field with an arbitrary nonminimal coupling with two axionic fields. The model also involves a parameter b that enters in the scalar potential and in the function that minimally couples the scalar and the axionic fields, see below (2.3a-2.3b). The range of this extra parameter will be fixed by some reality conditions as well as by demanding the solutions to have positive entropy. Asymptotically AdS planar 1 In higher dimensions D > 4, it has been shown that locally AdS black hole solutions with planar base manifold can emerge for scalar field nonminimally coupled to Lovelock gravity [11][12][13]. dyonic black hole solutions will be presented with axionic fields homogeneously distributed along the orthogonal planar coordinates of the base manifold for a priori any positive value of the nonminimal coupling parameter ξ . However, the positive entropy condition will considerably reduce the range of the permissible values of the nonminimal coupling parameter. As in the conformal case [17], the solutions only contain an integration constant denoted by ω. Interestingly enough, the location of the event horizon can be at two different positions depending on the sign of the parameter ω. Moreover, for ω < 0, the range of permissible values of the nonminimal coupling parameter is discrete, and the solution is shown to be extremal since its has zero temperature. Finally, the full DC conductivities associated to the charged black brane solutions will be computed following the recipes given in [19,20]. One of our motivations is precisely to identify the impact on the conductivities of the nonminimal coupling parameter.
The plan of the paper is organized as follows. In the next section, we present the model which consists of the Einstein gravity action with a negative cosmological constant with a source given by a self-interacting and nonminimally scalar field coupled to two axionic fields. In Sect. 3, the asymptotically AdS planar dyonic solutions are presented. In Sect. 4, a detailed analysis of the thermodynamic properties of the solutions through the Hamiltonian method is provided allowing to identify correctly the mass of the dyonic solutions. In the next section, following the perturbative method presented in [19,20], the full DC conductivities of the (non) extremal solutions are computed and the effects of the nonminimal coupling parameter are analyzed. The last section is devoted to our conclusions.

Model, field equations and black brane solutions
We consider a four-dimensional Einstein-Maxwell model with a negative cosmological constant, a nonminimally coupled and self-interacting scalar field φ together with two axionic fields ψ i for i = 1, 2 whose action is given by Here the Maxwell field strength is F μν = ∂ μ A ν − ∂ ν A μ , the nonminimal coupling parameter is denoted by ξ while the potential U b and the coupling ε b will be given below. For latter convenience, we shall opt for the following parametrization of the nonminimal coupling parameter ξ which maps the region ξ > 0 to n ∈ R\[0, 1]. 2 Note that the particular value ξ = 1 4 which corresponds to the limit n → ∞ will be treated separately. The potential U b (φ) and the coupling ε b (φ) associated to the axionic fields depend on a positive constant denoted by b whose range will be conditioned by some reality conditions as shown below It is also important to stress that in order to ensure a mechanism of momentum dissipation as well as a correct sign for the axionic kinetic term, the minimal coupling function ε b (φ) must be strictly positive everywhere. Before proceeding, we would like to stress that the minimal case ξ = 0 or equivalently n = 1 which has been fully treated in [18] is clearly excluded from our analysis because of the form of the expressions U b (φ) and ε b (φ). Also, we may note that for the conformal situation ξ = 1 6 or equivalently n = 3, the potential vanishes identically and the axionic minimal coupling ε b (φ) = 1, and this situation corresponds to the case already studied in [17].
The field equations obtained by varying the action (2.1) with respect to the metric, the Maxwell field A μ , the scalar field φ and the axionic fields ψ i read respectively where the different energy-momentum tensors are given by 2 A simple calculation shows that the range of values ξ < 0 is not compatible with the reality condition as defined by Eq. (2.6).
Dyonic black brane solutions of the field equations for = − 3 and with axionic fields homogeneously distributed along the two-dimensional planar base manifold can be found, and are given by where the intensity of the axionic fields ω is tied to the electric and magnetic charges by the relation We will now be as exhaustive as possible to enumerate the properties of these dyonic black brane solutions. For a nonminimal coupling parameter ξ as parameterized by (2.2), we have to distinguish between the case n > 1 corresponding to ξ ∈]0, 1 4 [, and the case n < 0, i.e. ξ ∈] 1 4 , ∞[. For n > 1, in order to deal with a real constant ω as defined by (2.6), the parameter b of the model must belong to the interval while for n < 0, one must have b ∈ ]b 0 , ∞[. Also, the constraint (2.6) allows both sign for the intensity of the axionic fields ω. This remark has an interesting consequence concerning the location of the horizon r + since for ω > 0, while for ω < 0, the event horizon is located at . In this latter case, because of the form of the metric function (2.5), the temperature of the solution vanishes identically and the black brane configuration can be interpreted as an extremal solution. It remains to ensure that the scalar field (2.5) is well-defined in the region outside the horizon, namely for r ∈ [r + , ∞[. For ω > 0, the scalar field is regular everywhere in the region r > 0. On the other hand, for a negative axionic intensity parameter ω < 0, the reality condition on the scalar field on the region r ∈ [r + = − ω √ 12κ , ∞[ restricts the parameter n to Sign of ω horizon r + Permissible values of n or equivalently of ξ and range of b be an odd negative integer n = 1 − 2k with k ∈ N\{0} or equivalently the nonminimal coupling parameters is forced to take the discrete values given by ξ = k 2(2k−1) . All these details are summarized in the Table 1. It is also interesting to stress that for ω < 0, even if the scalar field vanishes at the horizon, the expressions of the potential (2.3a) and the coupling (2.3b) remain finite once evaluated on the solution at the horizon, i.e.
This remark will be of importance for the finiteness of the DC conductivities of the extremal solutions. Finally, we would like to make a comment concerning the neutral configuration q e = q m = 0. In this case, the reality condition on ω given by Eq. (2.6) is replaced by a constraint on b given by To conclude this section, we briefly report the solution for the special coupling ξ = 1 4 which was excluded from the previous study. In contrast with the other couplings, the black brane solution can not be charged, and the neutral solution is given by for a potential U and a coupling ε free of any couplings that read One can easily see that for ω > 0, the metric solution admits a single root located at r + ∈] ω 2 3κ , ∞[.

Thermodynamics of the solutions by means of the Hamiltonian method
We now study the thermodynamic properties of the nonextremal solution (2.5) with ω > 0 whose event horizon is located at with n ∈ R\[0, 1] or equivalently 1 4 , ∞[. In order to achieve this task, we will proceed using the Euclidean approach where the time is imaginary and periodic with period β = T −1 . Here, T stands for the temperature which is fixed by requiring regularity at the horizon. The temperature of the non-extremal solution (2.5) reads The Euclidean action I Euc is related to the Gibbs free energy F by where M is the mass, S the entropy and μ i are the extra potentials with their corresponding charges Q i [29].
In order to construct a well-defined Euclidean action, we will consider a mini superspace where the metric is static and given by The radial coordinate r ranges from the horizon to infinity, i.e. r ≥ r + , and the planar coordinates both are assumed to belong to a compact set, that is x ∈ x and y ∈ y with dx dy = x y . We also assume a specific ansatz for the scalar, axionic and electromagnetic , and ψ 2 = ψ 2 (y). In doing so, the reduced Euclidean action takes the form where B is a boundary term that will be properly fixed below. Here p is the conjugate momentum of A τ , p = − r 2 N (r ) A τ (r ) , and the reduced Hamiltonian H is given by The boundary term B is fixed by requiring that the reduced action has an extremum, that is δ I Euc = 0, yielding to where this expression has to be evaluated at the limit R → ∞. The field equations obtained from varying the Euclidean action (which do not depend on the boundary term B) imply that N is a constant, and this latter can be chosen without any loss of generality to be N (r ) = 1. On the other hand, Gauss law implies that p = cst = q e . For the axionic fields, we note that their contribution strongly depends on the integral r r + ε b (φ(r )) dr, which in our case can be computed yielding to The final variation of δ B at the infinity is finite and given by where η is defined as The emergence of this extra contribution proportional to η which is effective only for n = 2 and for any nonzero negative integer can be explained from the fact that the "binomial" expression appearing between the large brackets in (3.4) will contribute with a linear term in the radial coordinate only for integer values n ∈ N less or equal to three. 3 Working in the grand canonical ensemble, where β and all the potentials are fixed, the boundary term at the infinity can be integrated as For the contribution at the horizon of the Euclidean action, we require the following variations For reason that will become clear, see below Eq. (3.10b), we definẽ and we finally get where A + = x y r 2 + is the horizon area, e = q e r + , m = q m r + and 1 and 2 correspond to the axionic potentials, defined as ω times the horizon term of (3.4), i.e., Therefore, the boundary term at the horizon is Plugging (3.6) and (3.9) we obtain our boundary term, B = B(∞) − B(r + ). The relation between the boundary term and the Gibbs free energy allows us to identify the mass M, the entropy S, the electric charge Q e , the magnetic charge Q m and the axionic charges Q 1 and Q 2 , It is reassuring to check that in the conformal case n = 3, these thermodynamic quantities reduce to those obtained in [17]. On the other hand, even if the first law is satisfied by means of the constraint (2.6), the signs of the entropy and the mass must be analyzed carefully. Firstly, the positivity of the entropy imposes a lower bound for the parameter b; this is because the constantG as defined by (3.7) must be strictly positive, and this leads to for n > 1, and b < 1 4 b 0 , for n < 0.
Referring to the Table 1, one can see that the reality condition imposed by (2.6) for ω > 0 make possible to deal with solutions with positive entropy only if n > 1 or equivalently ξ ∈]0, 1 4 [ and the range of the parameter b must then be reduced to b ∈] 1 4 b 0 , b 0 [. For the other values of the nonminimal coupling parameter n < 0 or equivalently ξ > 1 4 , the entropy of the solutions turns out to be always negative because of the reality condition (2.6). Concerning the mass, it is interesting to note that the mass is always positive for η = 0 as defined by (3.5) or for the conformal case n = 3. On the other hand, for n = 2, the positivity of the mass requires b > 1 2 b 0 , and hence our static solution with n = 2 can have positive entropy and mass by demanding the parameter b to belong to the set b ∈] 1 2 b 0 , b 0 [. Finally, for the remaining values for which η = 1, namely n = −1, −2, . . ., the solutions will always have negative mass and entropy.
For the extremal solution corresponding to ω < 0 and ξ taking discrete values (cf. Table 1) with a scalar field vanishing at the horizon φ(r + ) = 0, it is not safe to consider the Euclidean approach since the inverse of the temperature is infinite. Nevertheless, one can compute the entropy of the extremal solution by means of the Wald's formula [30] yielding S extremal = 4πκ x y r 2 + . (3.11) It is interesting to note that this result matches with the non extremal expression (3.10b) with (3.7) by taking φ(r + ) = 0.
In addition, as it can be seen from Eq. (3.11), the entropy of the extremal solution is always positive and given by the Hawking formula, S = Area 4G , after restoring correctly the value of κ in term of the Newton gravitational constant, i.e. κ = 1 16π G . On the other hand, the mass of the extremal solution is given by and it is easy to see that for the permissible discrete values of the nonminimal coupling parameters n = 1 − 2k with k ∈ N\{0} this expression is always positive for b > b 0 (cf. Table 1) even if ω < 0.

Holographic DC conductivities
In Refs. [19,20], it was established that the full DC conductivities for black holes enjoying a momentum dissipation can be computed in terms of the black hole horizon data. The main idea is to construct some conserved currents independent of the holographic radial coordinate r . To that end, we follow the prescription as described in Refs. [19,20], and we first turn on the following relevant perturbations on the black brane solution (2.5) δg t x = r 2 h t x (r ), δg r x = r 2 h r x (r ), δg ty = r 2 h ty (r ), δg r y = r 2 h ty (r ) where E x , E y are two constants. Consequently, the two perturbed Maxwell equations read f a x + f a x + q e h t x + q m ( f h r y + f h r y ) = 0, f a y + f a y + q e h ty − q m ( f h r x + f h r x ) = 0, which allow to define the following two conserved currents, According to the AdS/CFT correspondence, the holographic DC conductivities are determined by the conserved currents in the asymptotic boundary, and since the expressions obtained above (4.1a) and (4.1b) are independent of the radial coordinate, we can evaluate them at the horizon r + . The next step is to impose boundary conditions on the linearized perturbations at the black hole horizon. In order to achieve this task, it is convenient to use the Eddington-Finkelstein coordinates (v, r ) such that v = t + dr f (r ) . In this case, the gauge field will be well-defined by demanding Table 2 Signs of the entropy and mass of the dyonic solutions whose reality conditions are fixed by Table 1 while we will require the axionic fields to be constant near the horizon, and for the metric perturbations we will need to impose All the previous conditions can be substituted into the r x and r y components of the linearized Einstein equations, obtaining a system of equations for h t x (r + ) and h ty (r + ), whose solutions are given by Since the DC conductivities depend on the location of the horizon, we will first consider the non extremal case ω > 0 where the horizon is at r + = 3ω √ 12κ while the extremal situation ω < 0 and r + = − ω √ 12κ will be treated at the end of the section. Hence, the final step is to insert (4.2-4.4b) into (4.1a), (4.1b) to obtain the conductivities through where we have defined The conductivity matrix is antisymmetric in accordance with the invariance under the SO(2) symmetry. The purely electric DC conductivity σ DC is given by , where we have explicitly used the constraint (2.6). We are now in position to analyze the effects of the parameter b and the nonminimal coupling parameter n on the electrical DC conductivity. We also recall that the parameter b is subjected to some reality conditions (cf. Table 1) corresponding to the mathematical range of b but also to physical constraints that ensure the positivity of the mass and entropy (cf. Table 2) that we will refer as the physical range of b. Firstly, it is straightforward to prove that there exists a value of the nonminimal coupling parameter n 0 ≈ 3.4681 corresponding to ξ 0 ≈ 0.1779 (2.2) such that for each n > n 0 , there exists a precise value of the parameter b denoted by b 1 with b 0 4 < b 1 < b 0 (cf. Table 2) such that the DC conductivity is strictly positive for b ∈]b 1 , b 0 [ and σ DC becomes infinite at (n, b 1 ). In other words, this means that for any coupling greater than n 0 , one can always choose a parameter b ∈]b 1 , b 0 [ of the theory that yields to welldefined physical solutions, namely solutions with positive mass and entropy and having a positive conductivity. Additionally, for the choice b = b 1 , the DC conductivity becomes infinite at the point (n, b 1 ) with n > n 0 . This is due to the fact that the minimal coupling function ε b (φ) evaluated at the horizon vanishes at this point. This is not surprising since in order to ensure the mechanism of momentum dissipation yielding to finite conductivity, the coupling function ε b (φ) must not vanish.  Fig. 1 The conductivity σ DC in term of the parameter b for nonminimal couplings 0 < ξ < 1 4 Now, in order to visualize the influence on the electric DC conductivity of the parameters b and n, we will plot the graphics of σ DC in function of b and n. In Fig. 1, we plot the conductivity vs the parameter b for two distinct values of the nonminimal coupling parameter, namely n = 2 < n 0 and n = 4 > n 0 . The graphic given by Fig. 1a for n = 2 is in fact representative of all the cases 1 < n ≤ n 0 . One can see that for the physical range of b, Table 2), the conductivity is strictly positive and finite. 4 In fact, the positivity of the conductivity is always ensured even for the full mathematical range of b, i.e. b ∈]0, b 0 [. On the other hand, the graphic represented in Fig. 1b for n = 4 will be similar for any value n > n 0 . In this case, as mentioned before, one can see the existence of a a vertical asymptote at b = b 1 (for n = 4 we have b 1 ≈ 0.66534) where the conductivity becomes infinite. Nevertheless, in contrast with the previous case, the conductivity is positive only for b ≥ b 1 , and hence part of the physical range, namely b ∈] 1 4 b 0 , b 1 [ yields to negative conductivity. We also include a plot of the non extremal situation for n < 0 or equivalently ξ > 1 4 , see Fig. 2. The mathematical range of b is located at the right of the blue line while its physical range is at the left of the dotted line. One can see that for a mathematically well-defined solution, the electric conductivity is always positive even if it has neg- 4 It is worth mentioning that the values yielding negative conductivities are those for which the minimal coupling function ε b (φ) evaluated at the horizon is negative. Hence, in order to have a kinetic axionic term with the right sign to avoid ghost problems, the values associated to negative conductivities must be ruled out. ative mass and entropy. We can also appreciate the influence of the nonminimal coupling parameter on σ DC by drawing the graphics of this latter in function of n, see Fig. 3. In order to achieve this task, we must be careful with the election of b since its range of permissible values depends explicitly on n, see Tables 1 and 2 where we have defined b 0 in (2.7). As explained before, one can see that for n ≤ n 0 , the electric conductivity is positive while for n > n 0 , one has σ DC < 0. Finally, we would like to mention that as in the conformal case [17], the DC conductivity is temperature-independent.
To be more precise, on one hand we have seen that the temperature depends on the axionic parameter ω. On the other hand, in order to have solutions that respect the AdS symmetry at the boundary, the axionic charge must be related with the electric and magnetic charges through the relation (2.6). It is important to stress that this particular relation is due to the fact that the fall off of the nominally coupled scalar modifies the asymptotic behavior of our solutions. Moreover this relation makes our thermodynamic analysis integrable. However, as one uses the relation (2.6) in the conductivity matrix, we notice that ω-dependence disappears as it can be appreciated in the expression (4.7).
The remaining DC conductivities namely the thermoelectric and the heat conductivities can also be obtained by turning on a more general time-dependent perturbation 5 δ A x = t (−E + ζ a(r )) + a x (r ), where ζ can be shown to parameterize a time dependent source for the heat current. As a consequence, in addition to the conserved current J x , the following quantity is also conserved along the radial coordinate [19]. Hence, the thermoelectric conductivities α andᾱ as well as the heat conductivityκ can be computed yielding where β as usual is the inverse of the temperature. Finally, for the extremal solutions, namely ω < 0 (cf. Tables 1, 2), the analysis is completely analogue. This is due to the fact that, even if the scalar field vanishes at the event horizon, the minimal coupling functionε b remains finite at the horizon (2.8). As a direct consequence, the expressions derived previously remain valid with the difference that in the extremal case, the horizon is located at r + = − ω √ 12κ . Interestingly enough, the conductivity matrix becomes an antisymmetric matrix given by It is somehow appealing that for the extremal solution, the conductivity matrix has a Hall effect-like behavior with a Hall conductivity that looks like that B/ρ instead of ρ/B, (here ρ is the density of charge and B the orthogonal magnetic field). Indeed, in the AdS/CFT dictionary, (3+1)-dimensional AdS dyonic black holes are conjectured to be dual to a (2 + 1) CFT. In this picture, the electric bulk gauge field does not have a counterpart in the dual field theory but instead it fixes the electric charge density ρ to be proportional to the electric charge of the black hole, i.e. ρ ∝ q e . On the other hand, the magnetic bulk gauge field is in correspondence with an external magnetic field in the CFT side with a field strength B ∝ q m . Hence, the Hall conductivity in our case (4.9) is proportional to the ratio between the magnetic field and the (a) σ DC vs n for b = 0.6 (b) σ DC vs n for b = 0.95 Fig. 3 Influence of the nonminimal coupling parameter n on the electric conductivity electric charge density, i.e. σ xy ∝ B ρ . As a final comment, one can notice that in the purely electric case q m = 0, the DC conductivity for the extremal solution vanishes identically reproducing in turn an insulator behavior.

Conclusion
Here, we have considered a self-interacting scalar field nonminimally coupled to the four-dimensional Einstein gravity with a negative cosmological constant. The matter source is also supplemented by the Maxwell action with two axionic fields minimally coupled to the scalar field. Our model is specified from the very beginning by two parameters that are the nonminimal coupling parameter denoted ξ or equivalently n (2.2) and the constant b that enters in the minimal coupling as well as in the potential. For this model, we have obtained dyonic planar black holes with axionic fields depending linearly on the coordinates of the planar base manifold. We have noticed that these charged solutions depend on a unique integration constant denoted by ω and the horizon can be located at two different positions depending on the sign of ω. Surprisingly, for ω < 0, the temperature of the solution vanishes identically and hence one can interpret the solution as an extremal black brane. We have also shown that some reality conditions (cf. Table 1) supplemented by the requirement of having solutions with positive entropy and mass restrict considerably the permissible values of the nonminimal coupling parameter and of the parameter b, see Table 2. For a positive ω > 0, the set of physically acceptable values of the nonminimal coupling parameter is given by ξ ∈]0, 1 4 [ while for the extremal solution corresponding to ω < 0, only discrete values of the nonminimal coupling parameter given by ξ = k 2(2k−1) with k ∈ N\{0} yield to solutions with positive mass and entropy. These restrictions on the nonminimal coupling parameter are to be expected since, even for purely scalar field nonminimally coupled to Einstein gravity, black hole configurations have been shown to be ruled out for ξ < 0 and ξ ≥ 1 2 , see Ref. [31]. In the last part of this work, we have taken advantage of the momentum dissipation ensured by the axionic fields to compute the different conductivities by means of the recipes given in Refs. [19,20]. Many interesting results can be highlighted from the study of the holographic DC conductivities inherent to these dyonic solutions. For the non extremal solutions, we have shown that for n ≤ n 0 ≈ 3.4681 or equivalently ξ ≤ ξ 0 ≈ 0.1779, the dyonic solutions always enjoy a positive conductivity for any mathematically permissible value of b. On the other hand, for n > n 0 , the positive conductivity condition restricts the interval of b to be ]b 1 , b 0 [ with b 1 > 1 4 b 0 . In other words, this means that the physical solutions (in the sense of having positive mass and entropy) for n > n 0 with b ∈] 1 4 b 0 , b 1 [ will have a negative conductiv-ity. Also, we have shown that for n > n 0 there always exist a value of the parameter b denoted by b 1 which restores the translation invariance, in the sense that σ DC (n, b 1 ) → ∞. Finally, for the extremal solution, we have shown that the diagonal elements of the conductivity matrix precisely vanish and its off-diagonal elements are similar to those inherent to the Hall effect. An interesting extension of our model will be to consider an additional k-essence term for the axionic part of the action and to analyze the effects on the conductivities of the nonminimal coupling parameter ξ conjugated with the k-essence parameter, see Ref. [24]. In the same lines a natural generalization of these solutions would be the extension to higher dimensional scenarios following the lines of [32].
Finally, it will be very interesting to explore more deeply some of the properties of our solutions such as the extremality, the perfect conductivity or the Hall effect-like behavior of the extremal solutions. With this respect, in Ref. [33], it was shown that the Reissner-Nordstrom at the extremal limit experiences a sort of Meissner effect in the sense that the magnetic flux lines are expelled. Hence, a work to be done consists precisely in investigating the extremal solutions found here can exhibit a kind of Meissner effect.