Transport phenomena and Weyl correction in effective holographic theory of momentum dissipation

TL;DR

This paper develops a higher derivative holographic model with axionic fields and Weyl corrections to study momentum dissipation, revealing metal-insulator transitions and particle-vortex duality in transport properties.

Contribution

It introduces a novel higher derivative effective holographic theory incorporating Weyl corrections and axionic fields, analyzing momentum dissipation and transport phenomena.

Findings

Observed metal-insulator transitions at zero temperature.

Derived exact particle-vortex duality for specific dissipation strengths.

Calculated DC and optical conductivities in the model.

Abstract

We construct a higher derivative theory involving an axionic field and the Weyl tensor in four dimensional spacetime. Up to the first order of the coupling parameters, the charged black brane solution with momentum dissipation in a perturbative manner is constructed. Metal-insulator transitions are implemented when varying the system parameters at zero temperature. Also, we study the transports including DC conductivity and optical conductivity at zero charge density. We observe the exact particle-vortex duality for some specific momentum dissipation strength.

Figures

Equations

$$\displaystyle S_{0}=\int\mathrm{d}^{4}x\sqrt{-g}\left(R+\frac{6}{L^{2}}-\bar{\Phi}\right)\,,\$$

$$S_A=\int {\mathrm{d}}^{4}x\sqrt{-g}\left(-\frac{L^2}{8g_F^2}{F}_{\mu \nu }{X}^{\mu \nu \rho \sigma }{F}_{\rho \sigma }\right),$$

$$\displaystyle X_{\mu\nu}^{\ \ \rho\sigma}=I_{\mu\nu}^{\ \ \rho\sigma}-4\gamma_{1,0}L^{2}\bar{\Phi}I_{\mu\nu}^{\ \ \rho\sigma}-8\gamma_{0,1}L^{2}C_{\mu\nu}^{\ \ \rho\sigma}-8\gamma_{1,1}L^{4}\bar{\Phi}C_{\mu\nu}^{\ \ \rho\sigma}\,,\$$

$$\bar{\Phi }\equiv Tr[\Phi ]\equiv {\Phi }_{\mu }^{\mu },{\Phi }_{\nu }^{\mu }=\frac{1}{2}\sum _{I=x,y}{\partial }^{\mu }{\varphi }_I{\partial }_{\nu }{\varphi }_I.$$

$$\displaystyle\nabla_{\mu}\Big{[}\nabla^{\mu}\phi_{I}\Big{(}1-\gamma_{1,0}L^{4}F^{2}-\gamma_{1,1}L^{6}C^{\mu\nu\rho\sigma}F_{\mu\nu}F_{\rho\sigma}\Big{)}\Big{]}=0\,,$$

$${\nabla }_{\nu }(X^{\mu \nu \rho \sigma }{F}_{\rho \sigma })=0,$$

$$\displaystyle R_{\mu\nu}-\frac{1}{2}Rg_{\mu\nu}-\frac{3}{L^{2}}g_{\mu\nu}-\frac{L^{2}}{2}(1-4\gamma_{1,0}L^{2}\bar{\Phi})\Big{(}F_{\mu\rho}F_{\nu}^{\ \rho}-\frac{1}{4}g_{\mu\nu}F_{\rho\sigma}F^{\rho\sigma}\Big{)}$$

$$\displaystyle-\frac{L^{2}}{2}(1-\gamma_{1,0}L^{4}F^{2}-\gamma_{1,1}L^{6}C^{\mu\nu\rho\sigma}F_{\mu\nu}F_{\rho\sigma})\Big{(}\sum_{I=x,y}\partial_{\mu}\phi_{I}\partial_{\nu}\phi_{I}\Big{)}+\frac{1}{2}g_{\mu\nu}\bar{\Phi}$$

$$-L^4(\gamma +{\gamma }_{1,1}{L}^2\bar{\Phi })(G_{1\mu \nu }+G_{2\mu \nu }+G_{3\mu \nu })=0,$$

$$G_{1\mu \nu }$$

$$G_{2\mu \nu }$$

$$\displaystyle+2R_{\nu\alpha}F_{\mu}^{\ \beta}F^{\alpha}_{\ \beta}+2R_{\alpha\beta}F^{\alpha}_{\ \mu}F^{\beta}_{\ \nu}+2R_{\alpha\mu}F^{\alpha\beta}F_{\nu\beta}\,,\$$

$$G_{3\mu \nu }$$

$$\displaystyle S_{B}=\int\mathrm{d}^{4}x\sqrt{-g}\Big{(}-\frac{L^{2}}{8\hat{g}_{F}}G_{\mu\nu}\widehat{X}^{\mu\nu\rho\sigma}G_{\rho\sigma}\Big{)}\,,$$

$$\displaystyle\widehat{X}_{\mu\nu}^{\ \ \rho\sigma}=-\frac{1}{4}\varepsilon_{\mu\nu}^{\ \ \alpha\beta}(X^{-1})_{\alpha\beta}^{\ \ \gamma\lambda}\varepsilon_{\gamma\lambda}^{\ \ \rho\sigma}\,,\$$

$$\frac{1}{2}(X^{-1}{)}_{\mu \nu }^{\rho \sigma }{X}_{\rho \sigma }^{\alpha \beta }\equiv I_{\mu \nu }^{\alpha \beta },$$

$$(X^{-1}{)}_{\mu \nu }^{\rho \sigma }=\frac{1}{1-4{\gamma }_{1,0}{L}^2\bar{\Phi }}{I}_{\mu \nu }^{\rho \sigma }.$$

$${\hat{X}}_{\mu \nu }^{\rho \sigma }=(X^{-1}{)}_{\mu \nu }^{\rho \sigma }=\frac{1}{1-4{\gamma }_{1,0}{L}^2\bar{\Phi }}{I}_{\mu \nu }^{\rho \sigma }.$$

$$\displaystyle\mathrm{d}s^{2}=-\frac{r^{2}}{L^{2}}f(r)\mathrm{d}t^{2}+\frac{L^{2}}{r^{2}f(r)}\mathrm{d}r^{2}+\frac{r^{2}}{L^{2}}g(r)(\mathrm{d}x^{2}+\mathrm{d}y^{2})\,,\$$

$$A=A_t(r)\mathrm{d}t,$$

$$\displaystyle f(r)=f_{0}(r)+\gamma_{1,0}Y_{1,0}(r)+\gamma Y(r)+\gamma_{1,1}Y_{1,1}(r)\,,\$$

$$\displaystyle g(r)=1+\gamma G(r)+\gamma_{1,1}G_{1,1}(r)\$$

$$A_t(r)=A_{t0}(r)+{\gamma }_{1,0}{H}_{1,0}(r)+\gamma H(r)+{\gamma }_{1,1}{H}_{1,1}(r),$$

$$\displaystyle f_{0}(r)=1-\frac{M}{r^{3}}+\frac{q^{2}}{r^{4}}-\frac{\alpha^{2}L^{4}}{2r^{2}}\,,~{}~{}~{}~{}~{}~{}A_{t0}(r)=\mu-\frac{2q}{rL^{2}}\,,\$$

$$\displaystyle Y_{1,0}(r)=-\frac{4\alpha^{2}q^{2}L^{4}}{3r^{6}}\,,~{}~{}~{}~{}~{}~{}H_{1,0}(r)=\frac{8\alpha^{2}qL^{2}}{3r^{3}}\,,\$$

$$Y(r)=\frac{c_0{q}^2}{r^5}-\frac{c_{0}M}{2r^4}+\frac{c_1{\alpha }^2{L}^4}{2r^2}-\frac{c_0}{r}+\frac{20M{q}^2}{9r^7}-\frac{104q^4}{45r^8}+\frac{10{\alpha }^2{q}^2{L}^4}{9r^6}-\frac{32q^2}{9r^4},$$

$$G(r)=-\frac{c_0}{r}+c_1+\frac{4q^2}{9r^4},$$

$$\displaystyle H(r)=-\frac{c_{0}q}{r^{2}L^{2}}-\frac{4Mq}{r^{4}L^{2}}+\frac{296q^{3}}{45r^{5}L^{2}}-\frac{8\alpha^{2}qL^{2}}{9r^{3}}\,,\$$

$$Y_{1,1}(r)=\frac{d_{0}M}{2r^4}-\frac{d_0{q}^2}{r^5}+\frac{{\alpha }^2{d}_1{L}^4}{2r^2}+\frac{d_0}{r}+\frac{64{\alpha }^{2}M{q}^2{L}^4}{45r^9}-\frac{496{\alpha }^2{q}^4{L}^4}{315r^{10}}+\frac{28{\alpha }^4{q}^2{L}^8}{45r^8}-\frac{32{\alpha }^2{q}^2{L}^4}{45r^6},$$

$$G_{1,1}(r)=\frac{d_0}{r}+d_1+\frac{8{\alpha }^2{q}^2{L}^4}{45r^6},$$

$$H_{1,1}(r)=\frac{d_{0}q}{r^2{L}^2}-\frac{8{\alpha }^{2}Mq{L}^2}{3r^6}+\frac{208{\alpha }^2{q}^3{L}^2}{45r^7}-\frac{8{\alpha }^{4}q{L}^6}{15r^5}.$$

$$\displaystyle r\to r+\frac{1}{2}\gamma c_{0}-\frac{1}{2}d_{0}\gamma_{1,1}\,,\$$

$$(x,y)\to (x,y)\left(-\frac{1}{2}{d}_1{\gamma }_{1,1}-\frac{\gamma c_1}{2}+1\right),$$

$$\alpha \to \alpha \left(\frac{1}{2}{d}_1{\gamma }_{1,1}+\frac{\gamma c_1}{2}+1\right),$$

$$q$$

$$\displaystyle+\gamma_{1,1}\left(\frac{\alpha^{4}\mu L^{10}}{5r_{h}^{3}}+\frac{11\alpha^{2}\mu^{3}L^{10}}{90r_{h}^{3}}-\frac{2\alpha^{2}\mu L^{6}}{3r_{h}}\right)\,,\$$

$$M$$

$$+{\gamma }_{1,1}\left(\frac{8{\alpha }^4{\mu }^2{L}^{12}}{45r_h^3}+\frac{71{\alpha }^2{\mu }^4{L}^{12}}{630r_h^3}-\frac{22{\alpha }^2{\mu }^2{L}^8}{45r_h}\right).$$

$$r\to r_{h}r,(t,\vec{x})\to \frac{L^2}{r_h}(t,\vec{x}),A_t\to \frac{r_h}{L^2}{A}_t,M\to M{r}_h^2,Q\to Q{r}_h^2,\alpha \to \frac{r_h}{L^2}\alpha .$$

$$T=-\frac{2{\alpha }^2+{\mu }^2-12}{16\pi }-{\gamma }_{1,0}\frac{{\alpha }^2{\mu }^2}{12\pi }+\gamma \frac{{\mu }^2\left({\mu }^2-60\right)}{720\pi }+{\gamma }_{1,1}\frac{{\alpha }^2{\mu }^2\left(8{\alpha }^2+3{\mu }^2-84\right)}{360\pi }.$$

$$f(u)=(1-u)p(u),$$

$$p(u)=-\frac{1}{4}{\mu }^2{u}^3-\frac{{\alpha }^2{u}^2}{2}+u^2+u+1-\frac{1}{3}{\gamma }_{1,0}{\alpha }^2{\mu }^2{u}^3\left(u^2+u-1\right)$$

$$+\gamma \frac{1}{180}{\mu }^2{u}^3\left(2{\mu }^2\left(13u^4-14\right)+50{\alpha }^2\left(u^3-1\right)+\left({\mu }^2-100\right)\left(u^3+u^2+u\right)+240\right)$$

Full text

I Introduction

The quantum critical (QC) system has long been a central and challenging subject in condensed matter physics Sachdev:2000qpt . It is believed to account for the most interesting phenomena, such as the strange metal and pseudo-gap phase, in strongly correlated quantum materials. The QC system is associated with a QC phase transition and a QC phase. Since the QC system is strongly correlated, the conventional perturbative tools in traditional field theory, unfortunately, lose their power. We need to develop novel non-perturbative techniques and methods.

The AdS/CFT correspondence Maldacena:1997re ; Gubser:1998bc ; Witten:1998qj ; Aharony:1999ti , mapping a strongly coupled quantum field theory to a weakly coupled gravitational theory in the large N limit, provides a powerful tool to the study of QC physics and has led to great progress. Especially, the metal-insulator transition (MIT), a special example of the QC phase transition, has been widely studied in the holographic framework; for instance see Donos:2012js ; Donos:2013eha ; Donos:2014uba ; Ling:2015ghh ; Ling:2015dma ; Ling:2015epa ; Ling:2015exa ; Ling:2016wyr ; Ling:2016dck ; Baggioli:2014roa ; Baggioli:2016oqk ; Baggioli:2016oju ; Donos:2014oha ; Kiritsis:2015oxa and the references therein. To implement an MIT in a holographic framework, the key point is to deform the infrared (IR) geometry to a new fixed point by the introduction of momentum dissipation Donos:2012js ; Donos:2013eha .

Holographic QC phase at zero density has also been intensely explored in Myers:2010pk ; Ritz:2008kh ; WitczakKrempa:2012gn ; WitczakKrempa:2013ht ; Witczak-Krempa:2013nua ; Witczak-Krempa:2013aea ; Katz:2014rla ; Sachdev:2011wg ; Hartnoll:2016apf ; Bai:2013tfa ; Myers:2016wsu ; Lucas:2017dqa . By studying transport phenomena, in particular the optical conductivity, from a probe Maxwell field coupled to the Weyl tensor $C_{\mu\nu\rho\sigma}$ on top of the Schwarzschild-AdS (SS-AdS) black brane background Myers:2010pk ; Sachdev:2011wg ; Hartnoll:2016apf ; Ritz:2008kh ; WitczakKrempa:2012gn ; WitczakKrempa:2013ht ; Witczak-Krempa:2013nua ; Witczak-Krempa:2013aea ; Katz:2014rla , one observed a non-trivial frequency dependent conductivity attributed to the introduction of the Weyl tensor. It exhibits a peak, which resembles the particle response and we refer to this as the Damle-Sachdev (DS) peak Damle:1997rxu , or a dip, which is similar to the behavior of the vortex response, and is analogous to the one in the superfluid-insulator quantum critical point (QCP) 111Those kinds of peak and dip features have also been observed in probe branes and DBI action Chen:2017dsy and just higher terms in $F^{2}$ with $F$ being the Maxwell field strength Baggioli:2016oju . Myers:2010pk ; Sachdev:2011wg ; Hartnoll:2016apf .

But the peak is not the standard Drude peak and the DC conductivity has a bound, which cannot approach zero. When higher derivative (HD) terms are introduced, an arbitrarily sharp Drude-like peak can be observed at low frequency in the optical conductivity and the bound of conductivity is violated such that a zero DC conductivity can be obtained at a specific parameter 222We would like to point out that this bound in the conductivity is formalized in “almost” general theories in Grozdanov:2015qia ; Ikeda:2016rqh . But in more generic theories Baggioli:2016oqk ; Gouteraux:2016wxj ; Baggioli:2016pia , this bound is also violated. Witczak-Krempa:2013aea . Another step forward is the construction of a neutral scalar hair black brane by coupling the Weyl tensor with a neutral scalar field, which provides a framework to describe the QC phase and a transition away from QCP Myers:2016wsu ; Lucas:2017dqa .

In this paper, we shall construct a higher derivative theory including the four derivative terms, a simple summation of the Weyl tensor as well as a term from the trace of axions coupling with the gauge field, and a six derivative term, a mixed term of the product of the Weyl tensor and the axionic field coupling with the gauge field, and we obtain a charged black brane solution in a perturbative manner. By using a perturbative method, some charged black brane solutions from higher derivative gravity theory have been constructed; for instance see Ling:2016dck ; Myers:2009ij ; Liu:2008kt ; Cai:2011uh ; Dey:2015poa ; Dey:2015ytd and the references therein. Especially, in Ling:2016dck , it is the first time that an MIT is realized in the framework of higher derivative gravity. Along the line of Ling:2016dck , we shall study the MIT physics of our present model. Also, we explore the QC phase of this model at zero charge density.

We organize this paper as follows. In Sect. II, we construct the higher derivative model coupling axionic field and Weyl tensor with the gauge field. Then the perturbative black brane solution is obtained in Sect. III. In Sect. IV, we calculate the DC conductivity at finite charge density and study the MIT at zero temperature. The conductivity at zero charge density is explored in Sect. V. A brief discussion is presented in Sect. VI. The constraint on the coupling parameters is obtained in Appendix A.

II Holographic model

We construct a higher derivative holographic effective theory including metric, axions and gauge field as follows:

[TABLE]

where

[TABLE]

A pair of spatial linear dependent axionic fields, $\phi_{I}=\alpha x_{I}$ with $I=x,y$ and $\alpha$ being a constant, are introduced in the above action, which are responsible for dissipating the momentum of the dual boundary field. $L$ is the radius of the AdS spacetimes. $g_{F}$ and $\gamma_{m,n}$ with $m,n=0,1$ are the dimensionless coupling parameters. In what follows, we shall set $g_{F}=1$. $\Phi^{\mu}_{\ \nu}$ is the second order derivative term with respect to axions. The first term in the tensor $X$ gives the standard Maxwell term. $I_{\mu\nu}^{\ \ \rho\sigma}$ is an identity matrix defined as $I_{\mu\nu}^{\ \ \rho\sigma}=\delta_{\mu}^{\ \rho}\delta_{\nu}^{\ \sigma}-\delta_{\mu}^{\ \sigma}\delta_{\nu}^{\ \rho}$. The second term can be classified as a four derivative term, which is the term with $n=0,m=1$ in Gouteraux:2016wxj ; Baggioli:2016pia (Eq. (2.13) in Gouteraux:2016wxj ). The third term is also a four derivative one, constructed by the Weyl tensor, which has been well studied in Ling:2016dck ; Wu:2016jjd . For consistency with the current literature Myers:2010pk ; Wu:2016jjd ; Fu:2017oqa , we denote $\gamma_{0,1}=\gamma$ in what follows. The last term is a $6$ derivative term constructed by axions and the Weyl tensor. More higher derivative terms can be constructed in terms of axions, the Weyl tensor and the gauge field, which we leave for future study. It is easy to see that the new tensor $X$ possesses the same symmetry as $X_{\mu\nu\rho\sigma}=X_{[\mu\nu][\rho\sigma]}=X_{\rho\sigma\mu\nu}$ , like in Myers:2010pk ; Wu:2016jjd ; Fu:2017oqa .

The equations of motion (EOMs) can be straightforwardly derived from the above action (1),

[TABLE]

where

[TABLE]

Following Ref. Myers:2010pk , we can construct the electromagnetic (EM) dual theory of (1) with (2), which is

[TABLE]

where $\hat{g}_{F}^{2}\equiv 1/g_{F}^{2}$ and $G_{\mu\nu}\equiv\partial_{\mu}B_{\nu}-\partial_{\nu}B_{\mu}$. The tensor $\widehat{X}$ is defined by

[TABLE]

where $\varepsilon_{\mu\nu\rho\sigma}$ is a volume element. The tensor $\widehat{X}$ possesses the same symmetry as $X$, i.e., $\widehat{X}_{\mu\nu\rho\sigma}=\widehat{X}_{[\mu\nu][\rho\sigma]}=\widehat{X}_{\rho\sigma\mu\nu}$.

When $X_{\mu\nu}^{\ \ \rho\sigma}=I_{\mu\nu}^{\ \ \rho\sigma}$, the modified Maxwell theory (1b) is reduced the standard Maxwell one. In this case, one can easily deduce that $X^{-1}=X$ and so $\widehat{X}_{\mu\nu}^{\ \ \rho\sigma}=I_{\mu\nu}^{\ \ \rho\sigma}$ from Eqs. 7 and 6. Hence, the actions (1b) and (5) are identical, which demonstrates that the standard Maxwell theory is self-dual.

It has been shown in Myers:2010pk that, when the higher derivative term $\gamma$ is introduced, the EM self-duality is violated. Here, we demonstrate that even if only the $\gamma_{1,0}$ coupling term is introduced, the EM self-duality is also violated. We first evaluate the inverse of $X$ in terms of (7), which is

[TABLE]

Immediately, from Eq. (6), we find

[TABLE]

Since $\widehat{X}\neq X$, the EM self-duality is violated.

III Black brane solution

Since the EOM (3) are a set of third order differential equations with high nonlinearity, it has been hard to solve it analytically or even numerically so far. So following the strategy in Ling:2016dck (also see Myers:2009ij ; Liu:2008kt ; Cai:2011uh ; Dey:2015poa ; Dey:2015ytd ), we shall construct analytical solutions up to the first order of those coupling parameters. 333When the Weyl terms are turned off, i.e., $\gamma=0$ and $\gamma_{1,1}=0$, the black brane can be worked out analytically Gouteraux:2016wxj ; Baggioli:2016pia . We shall make a qualitative comparison on the DC conductivity between Gouteraux:2016wxj ; Baggioli:2016pia and our present results in Sec. IV.2. To this end, we take the following ansatz:

[TABLE]

where the UV boundary is at $r\rightarrow\infty$. Note that, when we take the following ansatz of $\phi_{I}$: $\phi_{I}=\alpha x_{I}$, Eq. (3a) satisfied automatically. So we only need to expand the functions $f(r)$, $g(r)$ and $A_{t}(r)$ in powers of $\gamma_{0,1}$, $\gamma$ and $\gamma_{1,1}$ up to the first order as

[TABLE]

where $f_{0}(r)$ and $A_{t0}(r)$ are the zeroth order solutions, which have been worked out in Andrade:2013gsa , while $Y_{i,j}(r)$, $G_{i,j}(r)$ and $H_{i,j}(r)$ are the first order solutions of $\gamma_{i,j}$. Note that we do not include the correction from $\gamma_{0,1}$ into the function of $g(r)$, so that we can make a direct comparison with the analytical solution in Gouteraux:2016wxj ; Baggioli:2016pia .

By directly solving Eq. (3) to the zeroth and first order of the coupling parameters, we can determine these functions:

[TABLE]

$(\mu,q,M,c_{0},c_{1},d_{0},d_{1})$ are seven integration constants, which are not independent from one another. Below, we shall derive the relations among them.

First, we can make the coordinate transformations

[TABLE]

and a redefinition of the axionic charge $\alpha$

[TABLE]

such that the integration constants $(c_{0},c_{1},d_{0},d_{1})$ can be eliminated. Using the conditions that $f$ and $A_{t}$ vanish at the horizon $r=r_{h}$, we obtain the relations for $(\mu,q,M)$:

[TABLE]

It is convenient to work with dimensionless quantities. So we make the following rescaling:

[TABLE]

Under this rescaling, we can set $L=1$ and $r_{h}=1$. Then the dimensionless temperature can be given by

[TABLE]

Note that all the above quantities $q$, $M$ and $T$ have been expanded to the first order of the coupling parameters $(\gamma_{1,0},\gamma,\gamma_{1,1})$. This black brane is characterized by two parameters, i.e., the temperature $T/\mu$ and the strength of the momentum dissipation $\alpha/\mu$. $\mu$ is interpreted as the chemical potential of the dual field and can be treated as the unit for the grand canonical system. For later convenience, we denote $\bar{T}\equiv T/\mu$ and $\bar{\alpha}\equiv\alpha/T$.

In addition, for the convenience of calculation, we shall work with the coordinate $u=1/r$. Then, in terms of $\mu$, we reexpress $f(u)$, $g(u)$ and $A_{t}(u)$ as follows:

[TABLE]

IV DC conductivity at finite density

IV.1 The derivation of the DC conductivity

In this section, we follow the procedure in Donos:2014uba ; Donos:2014cya ; Blake:2014yla ; Ling:2016dck to calculate the DC conductivity. To this end, we turn on the following consistent perturbations

[TABLE]

Then, one can define a radial conserved current in the bulk as

[TABLE]

Up to the first order of the coupling parameters, this conserved current can be evaluated as

[TABLE]

We have defined $Q=J^{t}$ in the above equation. It is the conserved electric charge density. Once $J^{x}$ is at hand, the DC conductivity can be evaluated in terms of Ohm’s law

[TABLE]

Since $J^{x}$ is a radial conserved quantity, the DC conductivity can be evaluated at the horizon $u=1$. Fist, we extract the value of $h_{tx}$ at the horizon from the $t$, $x$ component of the Einstein equation, which reads

[TABLE]

Notice that the above equation has taken value at $u=1$. In addition, we also need to add a regular boundary condition of $a_{x}$ at the horizon, which is

[TABLE]

Collecting Eqs. (21)-(24), we can obtain the DC conductivity:

[TABLE]

When $\gamma_{1,0}=0$ and $\gamma_{1,1}=0$, the result (25) reduces to Eq.(37) in Ling:2016dck . To compare with our present results, involving more coupling terms, we would like to present a brief review Ling:2016dck :

• •

There is a relation,

[TABLE]

which can be seen to hold when $\bar{\alpha}$ is fixed. It can be viewed as a special particle-vortex duality in Burgess:2000kj ; Murugan:2016zal .

• •

A metal-insulator transition (MIT) happens at zero temperature for a given nonzero $\gamma$ when we change the axionic charge $\bar{\alpha}$.

• •

There is a mirror symmetry at zero temperature 444This mirror symmetry also may be applicable at finite temperature.

[TABLE]

Next, we shall analyze the behavior of the DC conductivity and explore the MIT. Before proceeding, we introduce the definition of metallic phase and insulating phase adopted in many holography references Baggioli:2016rdj ; Donos:2012js ; Donos:2013eha ; Donos:2014uba ; Ling:2015ghh ; Ling:2015dma ; Ling:2015epa ; Ling:2015exa ; Ling:2016wyr ; Ling:2016dck ; Baggioli:2014roa ; Baggioli:2016oqk ; Baggioli:2016oju ; Donos:2014oha ; Kiritsis:2015oxa :

• •

Metallic phase: $\partial_{T}\sigma_{0}<0$.

• •

Insulating phase: $\partial_{T}\sigma_{0}>0$.

• •

Critical point (line): $\partial_{T}\sigma_{0}=0$.

IV.2 DC conductivity without Weyl term

In Appendix A, we analyze the causality and instabilities of the vector modes at zero density. When we only consider the $\gamma_{1,0}$ term, the analysis and the requirement of the positive DC conductivity indicate $-3/40\leq\gamma_{1,0}\leq 1/40$. But it is hard to analyze the causality and instabilities of the vector modes at finite density even if we have an analytical perturbative black brane solution. We shall leave this problem for future study. Here, we only approximately impose a further constraint from the requirement of the positive DC conductivity at finite density.

Figure 1 show the DC conductivity $\sigma_{0}$ as a function of the temperature $\bar{T}$ with different $\gamma_{1,0}$ and $\bar{\alpha}$. We find that, when $\gamma_{1,0}=0.025$, $\sigma_{0}$ is negative for small $\bar{\alpha}$ and low temperature $\bar{T}$. Further detailed exploration indicates that the positive definiteness of the DC conductivity constrains $\gamma_{1,0}$ in the range

[TABLE]

Also, we, respectively, show the DC conductivity as a function of $\bar{\alpha}$ for $\gamma_{1,0}$ belonging to the range (28) at zero temperature and finite temperature in Fig. 2. Figures 1 and 2 show that our result is qualitatively the same as that found in Gouteraux:2016wxj :

• •

At zero temperature, the DC conductivity monotonously decreases in terms of $\bar{\alpha}$.

• •

At finite temperature, the DC conductivity is qualitatively similar to that at zero temperature when $\gamma_{1,0}>0$. Meanwhile for $-3/40\leq\gamma_{1,0}<0$, the DC conductivity no longer monotonously decreases but has a minimum at some finite value of $\bar{\alpha}$.

• •

When $\bar{\alpha}$ is fixed, the DC conductivity monotonously decreases in terms of $\bar{T}$ for $\gamma_{1,0}>0$, which demonstrates a metal phase. When the sign of $\gamma_{1,0}$ changes, an opposite behavior is found, which is an insulator phase.

Therefore, our system up to the first order of the coupling parameters captures the main properties as shown in Gouteraux:2016wxj .

Finally, we present some comments on comparing with the $4$ derivative Weyl term studied in Ling:2016dck .

• •

Different from that for the four derivative Weyl term, no MIT happens for a given nonzero $\gamma_{1,0}$ when changing $\bar{\alpha}$ (see Fig. 3). But the mirror symmetry on $\frac{\partial\sigma_{0}}{\partial\bar{T}}(\bar{\alpha})$ (27) at zero temperature holds when the sign of $\gamma_{1,0}$ changes.

• •

Equation (26) holds when the sign of $\gamma_{1,0}$ changes and $\bar{\alpha}$ is fixed.

IV.3 DC conductivity from four derivative theory

When only the four derivative Weyl term $\gamma$ is involved, an MIT occurs at zero temperature by varying the axionic charge $\bar{\alpha}$. In particular, the quantum critical line is independent of the coupling parameter $\gamma$ Ling:2016dck .

In this section, we consider the mixed effect on DC conductivity in the four derivative theory including both $\gamma_{1,0}$ and $\gamma$ terms. The main properties are summarized as follows:

• •

Equations (26) and (27) hold for fixed $\bar{\alpha}$ and changing the signs of $\gamma$ and $\gamma_{1,0}$ (Fig. 4 and left plot in Fig. 5).

• •

For positive (negative) small $\gamma_{1,0}$, a MIT can be observed for negative (positive) $\gamma$ (see right plot in FIG.5). But different from the case only involving the four derivative term in Ling:2016dck , the quantum critical line is dependent on $\gamma$ (Fig.6). It provides a new platform of QCP such that we can study the holographic entanglement entropy and the butterfly effect close to QCP as in Ling:2016wyr ; Ling:2016ibq ; Ling:2016dck . We shall explore them in our present model in the future.

Before proceeding, we present some comments on the phase diagram for the MIT from four derivative theory at zero temperature (Fig. 6). For $\gamma_{1,0}<0$ and $\gamma>0$, with the increase of the strength of momentum dissipation, there is a phase transition from metallic phase to insulating one. This phenomenon is consistent with that of the usual charged particle excitations. On the other hand, for $\gamma_{1,0}>0$ and $\gamma<0$, we find that with the increase of the strength of momentum dissipation, the phase transition is opposite, i.e., the stronger momentum dissipates, the more insulating is the material. A better description of this phenomenon is provided by considering the excitations of vortices. Just as described Myers:2010pk , the EM duality of the bulk theory, which is related by changing the sign of $\gamma$, corresponds to the particle-vortex duality in the dual holographic CFT. Figure 6 shows such a duality; when we change the sign of $\gamma$, there is a duality between metallic and insulating phase. In fact, the phenomena can be easily concluded from Eq. (26). Finally, we would like to mention two corresponding examples. One is the transition observed in Myers:2010pk from the Drude-like peak at low frequency optical conductivity, which is interpreted as the charged particle excitations, to the dip, which resembles the excitations of vortices. Another one is the observation in Wu:2016jjd that the momentum dissipation drives the Drude-like peak into the dip of the low frequency optical conductivity for $\gamma>0$. Meanwhile for $\gamma<0$, the opposite scenario appears. When the sign of $\gamma$ changes, an approximate duality in optical conductivity is also observed for fixed strength of momentum dissipation. This duality is also observed in the next section.

IV.4 DC conductivity from six derivative theory

Now, we turn to a study of the effect of the six derivative term. For simplicity, we turn off the four derivative terms, i.e., we set $\gamma_{1,0}=0$ and $\gamma=0$. Figure. 7 exhibits the DC conductivity $\sigma_{0}$ as a function of the temperature $\bar{T}$ for some representative $\bar{\alpha}$ and $\gamma_{1,1}$. The left plot in Fig. 8 shows $\partial_{\bar{T}}\sigma_{0}$ as a function of $\bar{\alpha}$ at zero temperature, while right plot shows the phase diagram in the $(\gamma_{1,1},\bar{\alpha})$ plane for the MIT at zero temperature. We find that the properties of DC transport from six derivative theory is very similar to that from four derivative theory only involving the Weyl term Ling:2016dck ; here the mirror symmetries (26) and (27) hold for fixed $\bar{\alpha}$ and changing the signs of $\gamma_{1,1}$ in the phase diagram over $(\gamma_{1,1},\bar{\alpha})$ plane for the MIT at zero temperature. One difference is that the quantum critical line is shifted to $\bar{\alpha}\simeq 0.9$. 555The quantum critical line from four derivative theory only involving the Weyl term is located at $\bar{\alpha}\simeq 0.82$. It is also interesting to explore the DC conductivity at finite density from the six derivative theory only involving Weyl terms and compare the results with present results, including the mixed effect of both axions and Weyl tensor. We leave this problem for future study.

V Transports at zero density

In this section, we study the transports at zero density. In this case, the black brane solution reduces to the neutral one Andrade:2013gsa

[TABLE]

Note that we have parameterized this black brane solution by one scaling-invariant quantity $\hat{\alpha}=\alpha/4\pi T$ with $T=p(1)/4\pi$. Based on this neutral geometry background, we shall study the transport starting from four derivative and six derivative theory, respectively.

V.1 Four derivative theory

In this section, we study the properties of the conductivity in four derivative theory and see how the new higher derivative coupling term $\gamma_{1,0}$ affects them. Figure. 9 shows the optical conductivity $\sigma(\hat{\omega})$ as a function of $\hat{\omega}$ with representative $\gamma_{1,0}$, $\gamma$ and $\hat{\alpha}$. Comparing Fig. 9 with Fig. 1 in our previous work Wu:2016jjd , we observe that, for the system with positive (negative) $\gamma$ and $\gamma_{1,0}$, the transition from peak (dip) to dip (peak) appears to go easier with the increase of $\hat{\alpha}$.

As have been revealed in Wu:2016jjd , the particle-vortex duality is recovers with the change of $\gamma\rightarrow-\gamma$ for a specific value of $\hat{\alpha}=2/\sqrt{3}$. Now we want to explore if this phenomenon is generic when a new higher derivative coupling term $\gamma_{1,0}$ is taken into account. Figure. 10 shows the DC conductivity $\sigma_{0}$ as a function of $\hat{\alpha}$ for the representative $\gamma$ and $\gamma_{1,0}$. We find that, for a given $\gamma_{1,0}$, all the lines of $\sigma_{0}(\hat{\alpha})$ with different $\gamma$ intersect at one point $\hat{\alpha}=2/\sqrt{3}$, which is similar to that found for only the Weyl term $\gamma$ being involved. It indicates that $\sigma_{0}(\hat{\alpha})$ is independent of $\gamma$ for $\hat{\alpha}=2/\sqrt{3}$, which can also be deduced from the expression for DC conductivity (45). But we note that the value of $\sigma_{0}(\hat{\alpha}=2/\sqrt{3},\gamma)$ is not equal to unity. Also, the relation $\sigma_{0}(\hat{\alpha}=2/\sqrt{3},\gamma)=\frac{1}{\sigma_{0}(\hat{\alpha}=2/\sqrt{3},-\gamma)}$ does not hold. It indicates the exact duality of the DC conductivity only with the Weyl term for $\hat{\alpha}=2/\sqrt{3}$ is violated when the $\gamma_{1,0}$ term is taken into account. Furthermore, we study the optical conductivities of both the original EM theory and its dual theory for the specific value of $\hat{\alpha}=2/\sqrt{3}$, shown in Fig. 11, and we find that the exact particle-vortex duality is indeed violated when $\gamma\rightarrow-\gamma$ and $\gamma_{1,0}\rightarrow-\gamma_{1,0}$. It is easy to check that if we fix $\gamma_{1,0}$, the particle-vortex duality is also violated when $\gamma\rightarrow-\gamma$.

V.2 Six derivative theory

Now, we turn to a study of the case in six derivative theory. Figure 12 shows the optical conductivity with $\gamma_{1,1}$ being turned on. We observe that for positive $\gamma_{1,1}$ and small $\hat{\alpha}$, a small peak is displayed in the low frequency region. With the increase of $\hat{\alpha}$, the small peak starts to develop into a dip (left plot in Fig. 12). Meanwhile for negative $\gamma_{1,1}$, an opposite scenario is found (right plot in Fig. 12). The phenomenon is similar to that with the $\gamma$ term.

Also, we note that, for the specific value of $\hat{\alpha}=2/\sqrt{3}$, the DC conductivity $\sigma_{0}=1$ and is independent of $\gamma_{1,1}$ (see Fig. 23), which is similar to that with only the Weyl term Wu:2016jjd . Furthermore, we study the particle-vortex duality of this case, which is shown in Fig. 13. It is obvious that for small $\gamma_{1,1}$, the particle-vortex duality approximately holds. Meanwhile, for the specific value of $\hat{\alpha}=2/\sqrt{3}$, the duality exactly holds. Though here we do not work out the analytical understanding on the particle-vortex duality for the specific value of $\hat{\alpha}=2/\sqrt{3}$, it seems to originate from the Weyl term. The additional $\gamma_{1,0}$ term violates this exact duality. Further, we examine the duality from another six derivative term with $X_{\mu\nu}^{\ \ \rho\sigma}=-4\gamma_{1}C^{2}I_{\mu\nu}^{\ \ \rho\sigma}$, of which the original theory has been studied in our previous work Fu:2017oqa . Again, the particle-vortex duality exactly holds for $\hat{\alpha}=2/\sqrt{3}$ when $\gamma_{1}\rightarrow-\gamma_{1}$ (see Fig. 14). In future, we will further test the robustness of this phenomena by exploring that with the higher order terms of the Weyl coupling.

VI Discussions

In this work, we extend our previous work Ling:2016dck ; Wu:2016jjd to constructing a higher derivative theory including the coupling among the axionic field, the Weyl tensor and the gauge field. To be more specific, we construct four derivative terms, a simple summation of the Weyl term $C_{\mu\nu\rho\sigma}$ coupling with the gauge field, as well as a term from the trace of axions coupling with the gauge field, and a six derivative term, a mixed term by the product of Weyl tensor and the axionic field, coupling with the gauge field.

Following the strategy in Ling:2016dck , we construct the charged black brane solution with momentum dissipation in a perturbative manner up to the first order of the coupling parameters. We study the QCP from four and six derivative theory, respectively. For four derivative theory, because of the introduction of $\gamma_{1,0}$, the quantum critical line is independent of $\gamma$, which is different from the case only involving the $4$ derivative term in Wu:2016jjd . It provides a new platform of QCP such that we can study holographic entanglement entropy and butterfly effect close to QCP, which may inspire new insight. For six derivative theory, the quantum critical line is independent of the coupling parameter $\gamma_{1,1}$, which is similar that in Ling:2016dck .

Also, we study the transport phenomena including DC conductivity and optical conductivity at zero charge density, which is away from the QC phase. For four derivative theory, the momentum dissipation makes the transition from peak (dip) to dip (peak) easier, comparing with that in our previous work Wu:2016jjd . In addition, we find that for the specific value of $\hat{\alpha}=2/\sqrt{3}$, the exact particle-vortex duality, holding for only the $\gamma$ term, survives Wu:2016jjd and is violated when the $\gamma_{1,0}$ term is turned on. For the six derivative theory, the particle-vortex duality exactly holds for $\hat{\alpha}=2/\sqrt{3}$. Meanwhile the effect of the momentum dissipation on the transition between the gap and the dip is similar to that in four derivative theory.

It is definitely a novelty and an interesting matter to compute the optical conductivity at finite chemical potential $\mu$. However, even if we have obtained the perturbative black brane solution to the first order of $\gamma$ in Sect. III, we still need to solve the linear perturbative differential equations beyond the second order to obtain the optical conductivity. It is a hard task and so we shall leave it for the future. In addition, this simple model including the mixed terms between the Weyl tensor and the axions can be straightforwardly generalized to include the charge complex scalar field such that we can study the superconducting phase. It is also interesting and valuable to further explore the transport of our present model at full momentum and energy spaces, which certainly will reveal more information of the systems. This work deserves further study and we plan to publish our results in the near future.

Acknowledgements.

We in particular thank Zhenhua Zhou for the very helpful discussion on the calculation of DC conductivity at finite density. We are grateful to Peng Liu and Wei-Jia Li for helpful discussions. We also thank the anonymous referee for his/her very valuable suggestions, greatly improving our manuscript. This work is supported by the Natural Science Foundation of China under Grant nos. 11775036, 11305018 and the Natural Science Foundation of Liaoning Province under Grant no. 201602013.

Appendix A Bounds on the coupling

In this appendix, we explore the constraints on the coupling parameters. We mainly examine the causality of the dual boundary theory, the instabilities of the vector modes and the positive definiteness of the DC conductivity at zero charge density. We also discuss the constraint from the requirement that the graviton mass is real, i.e., $m_{g}^{2}>0$.

A.1 Bounds on the coupling at zero charge density

To examine the causality of the dual boundary theory and the instabilities of the vector modes, we decompose the perturbations of gauge field in the Fourier space as $A_{\mu}(t,x,y,u)\sim\mathrm{e}^{i\bf{q}\cdot\bf{x}}A_{\mu}(u,\bf{q}),$ with $\textbf{q}\cdot\textbf{x}=-\omega t+q^{x}x+q^{y}y$, and write down the EOMs as follows:

[TABLE]

where the prime denotes the derivative with respect to $u$ and the dimensionless frequency and momentum $\hat{\omega}\equiv\frac{\omega}{4\pi T}=\frac{\omega}{\mathfrak{p}}\,,\hat{q}\equiv\frac{q}{4\pi T}=\frac{q}{\mathfrak{p}},$ with $\mathfrak{p}\equiv p(1)=4\pi T$, are introduced. Due to the rotational symmetry in $xy$-plane, we have set $\textbf{q}^{\mu}=(\omega,q,0)$. Also we choose the gauge as $A_{u}(u,\textbf{q})=0$. At the same time, a tensor $X_{\mu\nu}^{\ \ \rho\sigma}$ defined as $X_{A}^{\ B}=\{X_{1}(u),X_{2}(u),X_{3}(u),X_{4}(u),X_{5}(u),X_{6}(u)\},$ with $A,B\in\{tx,ty,tu,xy,xu,yu\},$ has been introduced to simplify the expression of the perturbative EOMs. Since the background is rotationally symmetric in the $xy$-plane, we have $X_{1}(u)=X_{2}(u)$ and $X_{5}(u)=X_{6}(u)$. Combining Eqs. (30) and (31), one has a decoupled EOM for $A_{t}(u,\hat{\bf{q}})$, which is

[TABLE]

By making a transform as $A_{\mu}\rightarrow B_{\mu}$ and $X_{i}\rightarrow\widehat{X}_{i}$, we can obtain the EOMs of the dual EM theory from the above equations. Note that from Eq. (6), it is easy to deduce that $\widehat{X}_{A}^{\ B}$ is also diagonal with $\widehat{X}_{i}=1/X_{i}$.

Since $A_{x}$ can be expressed by $A_{t}$ in terms of Eq. (30), there are only two independent vector modes, $A_{t}$ and $A_{y}$, which correspond to EOMs (34) and (33). They can be formulated in Schrödinger form as

[TABLE]

Notice that we have made a coordinate transformation, $dz/du=\mathfrak{p}/f$, and a separation of variable, $A_{i}(u)=G_{i}(u)\psi_{i}(u)$, where $A_{\bar{t}}(u):=A^{\prime}_{t}(u)$ and $i=\bar{t},y$. For later convenience, we decompose the effective potential $V_{i}(u)$ into both a momentum dependent part and an independent one,

[TABLE]

where Witczak-Krempa:2013aea

[TABLE]

Before proceeding, we present the main ingredients constraining the coupling parameters as follows.

• •

If $V_{i}(u)$ satisfies,

[TABLE]

the modes meet the requirements of both causality and the stability of the dual boundary theory Buchel:2009tt ; Brigante:2008gz ; Myers:2007we .

• •

When $V_{i}(u)$ violates the lower bound, the modes may be instable. We need further analyze the zero energy bound state of the potential.

• •

An additional condition is the requirement of positive definiteness of the real part of the conductivity, especially the DC conductivity.

Next, we analyze the constraint on the coupling parameters.

A.1.1 Four derivative theory

When only the coupling parameter $\gamma_{1,0}$ survives, some related discussions have been explored in Gouteraux:2016wxj . But here one only discusses the Schr$\ddot{o}$dinger potential of the perturbation $A_{x}$. Here, we shall present a more detailed discussion in our present framework.

In terms of the expression of the DC conductivity Ritz:2008kh ; Myers:2010pk

[TABLE]

we can explicitly write it down when only $\gamma_{1,0}$ survives,

[TABLE]

Figure 15 shows $\sigma_{0}$ as a function of $\hat{\alpha}$ for sample values of $\gamma_{1,0}$. We see that, for $\gamma_{1,0}\leq 1/24$, $\sigma_{0}$ is positive for all values of $\hat{\alpha}$. Meanwhile, for $\gamma_{1,0}>1/24$, it vanishes for some finite $\hat{\alpha}$. This can also be seen from the following: when $\hat{\alpha}\rightarrow+\infty$, $\sigma_{0}=1-24\gamma_{1,0}$. Therefore, a non-negative $\sigma_{0}$ gives a constraint on $\gamma_{1,0}$ as $\gamma_{1,0}\leq 1/24$.

Next, we turn to a discuss of the bounds of $\gamma_{1,0}$ imposed by the causality and the instabilities. First, it is easy to find that in the limit of large momentum, since $X_{1}=X_{3}=1-4\mathfrak{p}^{2}\hat{\alpha}^{2}\gamma_{1,0}u^{2}$ for only $\gamma_{1,0}$ we have surviving $V_{0\bar{t}}=V_{0y}=f(u)$, which are the dominant terms. Obviously, $V_{0\bar{t}}$ and $V_{0y}$ are independent of the parameter $\gamma_{1,0}$ and satisfy the constraint (40). Meanwhile for the case of the small momentum region, the dominant terms are $V_{1,i}$ ($i=\bar{t},y$), which are shown for representative values of $\gamma_{1,0}$ and $\hat{\alpha}$ in Fig. 16. We can see that there is a negative minimum in $V_{1i}$. So we need to analyze the zero energy bound state of the potentials, which is Myers:2007we

[TABLE]

where $n$ is a positive integer and the potential well in the integral interval $[u_{0},u_{1}]$ is negative. Both $\tilde{n}_{1\bar{t}}$ and $\tilde{n}_{1y}$ as a function $\gamma_{1,0}$ for representative $\hat{\alpha}$ are exhibited in Fig. 17. The detailed analysis indicates that, when $\gamma_{1,0}$ belongs to the region $\gamma_{1,0}<-3/40$ and $1/40<\gamma_{1,0}<1/24$, the $\tilde{n}_{1i}$ are greater than unit and unstable modes develop. Combining the observation from DC conductivity, we can infer that the allowed region for $\gamma_{1,0}$ is

[TABLE]

Also, we have checked that, for finite momentum, no unstable mode appears for the constraint (44). Note that the lower bound of $\gamma_{1,0}$ is consistent with that found in Gouteraux:2016wxj , but the upper bound becomes tighter than that in Gouteraux:2016wxj , which results from the instability of the mode $A_{t}$.

Now, we begin to discuss the bounds on the coupling when both $\gamma_{1,0}$ and $\gamma$ are turned on. We mainly restrict $\gamma$ to the region $-1/12\leq\gamma\leq 1/12$ and explore the constraint of $\gamma_{1,0}$. First, in this case we derive the DC conductivity as

[TABLE]

which we plot as a function of $\hat{\alpha}$ for sample values of $\gamma_{1,0}$ and $\gamma$ in Fig. 18. Since with the increases of $\hat{\alpha}$ the positive $\gamma$ lowers the DC conductivity, it gives a tighter constraint on $\gamma_{1,0}$. Specially, when $\gamma=1/12$, to have positive $\sigma_{0}$, $\gamma_{1,0}\leq 1/36$ should be imposed. It can also be deduced thus: in the limit of $\hat{\alpha}\rightarrow+\infty$, $\sigma_{0}=1-4\gamma-24\gamma_{1,0}$. Second, we examine the potential $V_{0,\bar{t}}$, which is shown in Fig. 19. We see that, for $\gamma_{1,0}=1/36$ and $\gamma=-1/12$, an infinite positive and negative well appears in the limit of $\hat{\alpha}$, which signals an instability. This instability is due to turning on of $\gamma_{1,0}$ for $\gamma=-1/12$. We find that, when we tune $\gamma_{1,0}$ to become smaller, so that $\gamma_{1,0}\leq 7/500$, the infinite well gradually disappears (see Fig. 20). Therefore, if we set $-1/12\leq\gamma\leq 1/12$, then the constraint $-3/40\leq\gamma_{1,0}\leq 7/500$ should be imposed. At the same time, it is easy to see that, for the above range of $\gamma$ and $\gamma_{1,0}$, the potential $V_{0,y}$ satisfies the constraint (40). Third, we analyze the potential $V_{1\bar{t}}$, which is shown in Fig. 21. We see that $V_{1\bar{t}}$ develops a negative minimum. So to determine the range of parameter $\gamma_{1,0}$, we study $\tilde{n}_{1\bar{t}}$ as a function of $\gamma_{1,0}$ for the representative values of $\gamma$ and $\hat{\alpha}$, which are plotted in Fig. 22. A detailed analysis indicates that, when $\gamma_{1,0}$ belongs to the region $-3/40\leq\gamma_{1,0}\leq 1/100$, no unstable mode appears. Therefore, the constraint on $\gamma_{1,0}$ and $\gamma$ is

[TABLE]

A.1.2 Six derivative theory

In this subsection, we study the bounds on the coupling $\gamma_{1,1}$ with other coupling vanishing. First, we derive the DC conductivity:

[TABLE]

We plot it as a function of $\hat{\alpha}$ for sample values of $\gamma_{1,1}$ in FIG.23, in which we see that there are lower and upper bounds set by DC conductivity. By detailed analyzing, we find that $-1/3\leq\gamma_{1,1}\leq 1/24$. Specially, the upper bound can be deduced from that in the limit of $\hat{\alpha}\rightarrow+\infty$, $\sigma_{0}=1-24\gamma_{1,1}$.

Then, we consider the constraint from $V_{0\bar{t}}$, which is shown in Fig.24. We see that, for $\gamma_{1,1}=-1/3$, with the increase of $\hat{\alpha}$, the condition (40) is violated, which indicates that a tighter lower bound should be imposed on $\gamma_{1,1}$. A detailed examination indicates that $-1/50\leq\gamma_{1,1}\leq 1/24$. Also, we examine $V_{0y}$ for this range $\gamma_{1,1}\in[-1/50,1/24]$ and find that it satisfies the condition (40).

Now, we examine the potential $V_{1,\bar{t}}$, which we plot in Fig. 25. As in the previous case, a negative minimum appears in $V_{1,\bar{t}}$. So we further plot $\tilde{n}_{1\bar{t}}$ as a function of $\gamma_{1,1}$ for the representative values of $\hat{\alpha}$, which are shown in FIG.26. We find that, for the range

[TABLE]

no unstable mode appears. A similar analysis also indicates that, for $\gamma_{1,1}$ satisfying the constraint (48), $\tilde{n}_{1y}\leq 1$. In addition, this range of $\gamma_{1,1}$ is also a physically viable region for finite momentum.

A.2 Bounds on the coupling at finite charge density

In this section, we discuss the bounds on the coupling at finite charge density on top of the perturbative black brane geometry in Sect. III. Since the perturbative equations of vectors involve a set of third order differential equations with high nonlinearity, it is hard to decoupling them at finite charge density, like that at zero charge density. Therefore, it is difficult to study the bounds on the coupling at finite charge density by the method of Schrödinger potentials at zero charge density as Appendix A.1 or the quasi-normal modes of vector modes. We hope that these problems can be worked out in the future. Here, we only give the constraints on the coupling parameters at finite charge density from the requirement that the mass of the graviton is real.

It has been demonstrated in Blake:2013owa (also refer to Baggioli:2014roa ; Alberte:2015isw ) that the holographic lattices give the graviton an effective mass. In our present model (1), the effective graviton mass is

[TABLE]

Obviously, $m_{g}^{2}>0$ for the case of zero charge density. In what follows, we shall discuss the bounds on the coupling parameters at finite charge density.

We first turn on $\gamma_{1,0}$. Figure 27 shows $m_{g}^{2}$ as a function of $u$ at zero temperature for representative $\bar{\alpha}$ and $\gamma_{1,0}$. We can see that $m_{g}^{2}>0$ when $\gamma_{1,0}$ satisfies the constraint (44), which is the constraint at zero charge density.

Second, we turn on both $\gamma_{1,0}$ and $\gamma$. Figure 28 shows that if $\gamma_{1,0}$ and $\gamma$ satisfy the constraint (46), set at zero charge density, $m_{g}^{2}>0$ for large $\bar{\alpha}$ (right plot in Fig. 28), but $m_{g}^{2}>0$ is violated for small $\bar{\alpha}$ (left plot in Fig. 28). By detailed analysis, we constrain $\gamma_{1,0}$ and $\gamma$ in the region

[TABLE]

Finally, we analyze the constraint on $\gamma_{1,1}$ ($\gamma_{1,0}$ and $\gamma$ are turned off). The left plot in Fig. 29 shows that, when $\gamma_{1,1}$ reaches its lower bound set at zero charge density, $m_{g}^{2}$ becomes negative for $u$ approaching the horizon. Further analysis indicates that if $\gamma_{1,1}$ satisfies

[TABLE]

$m_{g}^{2}>0$ (see the right plot in Fig. 29).