This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

thanks: [email protected]thanks: [email protected]thanks: [email protected]thanks: [email protected]

Holograhic two-currents model with coupling and its conductivites

Dan Zhang 1    Zhenhua Zhou 2    Guoyang Fu 3    Jian-Pin Wu 3 1 College of Physical Science and Technology, Bohai University, Jinzhou 121013, China 2 School of Physics and Electronic Information, Yunnan Normal University, Kunming, 650500, China 3 Center for Gravitation and Cosmology, College of Physical Science and Technology, Yangzhou University, Yangzhou 225009, China
Abstract

We implement a holographic gravity model of two gauge fields with a coupling between them, which is dual to a two-currents model. An analytical black brane solution is obtained. In particular, we work out the expressions of conductivities with coupling and find that the expressions of conductivities are directly related to the coupling parameter θ\theta. It is the main topic of our present work. As an application, we calculate the conductivities by the scheme outlined here and briefly discuss the properties of the conductivities. An interesting property is that as the coupling θ\theta increases, the dip at low frequency in Re[σA]Re[\sigma_{A}]/Re[σB]Re[\sigma_{B}] becomes deepening and then turns into a hard-gap-like when θ=1\theta=1, which is independent of the doping χ\chi. Some monotonic behaviors of the conductivities are also discussed.

I Introduction

AdS/CFT (Anti-de-Sitter/Conformal field theory) correspondence, also referring to holography, provides us a way to study the dynamics of certain strongly-coupled condensed matter systems Maldacena:1997re ; Gubser:1998bc ; Witten:1998qj ; Aharony:1999ti . In the so-called bottom-up approaches, we can study some universal properties of the dual system by constructing a simple gravitational model. Along this direction, some interesting holographic models, for example, the holographic superconductor Hartnoll:2008vx , holographic metal insulator phase transition (MIT) Donos:2012js ; Ling:2014saa and (non-) Fermi liquid Liu:2009dm , have been implemented. These studies provide some physical insight into the associated mechanisms of the strongly-coupled systems and the universality class of them.

Recently, the holographic two-currents models are increasingly drawn attention, see Kiritsis:2015hoa ; Bigazzi:2011ak ; Huang:2020iyw ; Iqbal:2010eh ; Baggioli:2015dwa ; Rogatko:2017tae ; Rogatko:2020vtz ; Seo:2016vks and references therein. In these models, a pair of U(1)U(1) gauge fields AA and BB in bulk are introduced111Some holographic models with two gauge fields are also studied in Ling:2015exa ; Ling:2020mwm ; Ling:2017naw ; Ling:2016wyr ; Tarrio:2011de ; Alishahiha:2012qu , in which only one of gauge fields is treated as the real Maxwell field and we only concentrate on its transport properties. In Ling:2015exa , another gauge field is treated as an auxiliary field, which is introduced to obtain an insulating phase with a hard gap. On top of this novel holographic insulator constructed in Ling:2015exa , a holographic superconductor is also built Ling:2017naw . Further, Ling:2020mwm also introduces a coupling between these two gauge fields and studies the superconducting instability. In Ling:2016wyr , another gauge field is introduced to induce metal-insulator phase transition over Gubser-Rocha background Gubser:2009qt in the limit of of zero temperature. While in Tarrio:2011de ; Alishahiha:2012qu , the additional gauge field also plays the role of the auxiliary field to implement charged hyperscaling or Lifshitz black hole background. Therefore, these models are still treated as single current models.. Therefore, we have two independently conserved currents, which relate to different kinds of chemical potentials or charged densities in the dual boundary field theory. The mismatch of the two controllable chemical potentials or charge densities induces the unbalance of numbers. In Kiritsis:2015hoa ; Baggioli:2015dwa ; Huang:2020iyw , the ration of the two chemical potentials is proposed to simulate the effect of doping. In Bigazzi:2011ak , they propose the holographic two-currents model has a counterpart of Mott’s two-currents model Mott:1936v1 ; Mott:1936v2 . The chemical potential mismatch is interpreted as a chemical potential for a U(1)BU(1)_{B} “spin” symmetry Bigazzi:2011ak ; Iqbal:2010eh . And then, on top of the two-current model, they construct an unbalanced s-wave superconductor by introducing a charged complex scalar field Bigazzi:2011ak , which is a simple extension of holographic superconductor in Hartnoll:2008vx . Also they study “charge” and “spin” transport Bigazzi:2011ak 222Some related works are also explored, see for example Erdmenger:2011hp ; Dutta:2013osl ; Correa:2019ivh ; Musso:2013ija ; Alsup:2012kr ; Hafshejani:2018svs .. In Rogatko:2017tae ; Rogatko:2020vtz ; Seo:2016vks , the holographic two-currents model is used to describe the nature of graphene.

Most of these works do not contain the coupling between two gauge fields. Recently, there have been a small number of works beginning to concentrate on the effect of the coupling between two gauge fields Baggioli:2015dwa ; Rogatko:2017tae ; Rogatko:2020vtz . The coupling between two gauge fields provides additional degree of freedom in the dual boundary field theory. In Baggioli:2015dwa , they build a holographic superconductor model containing a non-trivial higher derivative therm of axionic field breaking translational symmetry, a complex scalar field breaking U(1)U(1) symmetry and two U(1)U(1) gauge fields. In particular, the complex scalar field provides a non-trivial coupling between two gauge fields. In this way, they implement a superconducting dome-shaped region on the temperature-doping phase diagram. But the study of the transport properties of this model is absent. In Rogatko:2017tae , they introduce a simple coupling term between two gauge fields and explored its thermoelectric transport properties of its holographic dual boundary field theory describing graphene. In Rogatko:2020vtz , they show that there is a bound on the conductivity depending on the coupling between both gauge fields. At the same time, their study also indicates that even strong disorder cannot still induce a MIT in holographic two-currents model as that with single current Grozdanov:2015qia .

In this paper, we shall construct a holographic two-currents model with coupling. We derive the expression of the frequency dependent thermoelectric transport and explore its properties. Especially, we concentrate on the effect of the coupling between two gauge fields. Our paper is organized as what follows. In Sec.II, we describe the holographic framework of the two-currents model with coupling and work out the analytical double charged RN-AdS black brane solution. In Sec.III, by standard holographic renormalized procedure, we obtain the expressions of the holographic conductivities for our two-currents model with coupling. And then, the properties of the conductivities are discussed in Sec.IV. The conclusions and discussions are presented in Sec.V.

II The double charged RN-AdS black brane

The action of the gravity dual for a simple two-currents model with coupling is

S=d4xg(R+6L2+M),\displaystyle S=\int d^{4}x\sqrt{-g}\,\left(R+\frac{6}{L^{2}}+\mathcal{L}_{M}\right)\,,
M=14FμνFμν14GμνGμνθ2FμνGμν,\displaystyle\mathcal{L}_{M}=-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}-\frac{1}{4}G_{\mu\nu}G^{\mu\nu}-\frac{\theta}{2}F_{\mu\nu}G^{\mu\nu}\,, (1)

where we set the AdS radius L=1L=1 in this paper. F=dAF=dA and G=dBG=dB are the field strengths of the two gauge fields AA and BB, respectively. Above we bring in a coupling term between two gauge fields and θ\theta denotes the coupling strength. Notice that when θ2=1\theta^{2}=1, we can combine both gauge fields AA and BB into a new gauge field such that both gauge fields are indistinguishable.

From the action (II), we derive the equations of motion (EOMs) as what follows

Rμν12Rgμν3L2gμν12Tμν(A)12Tμν(B)θTμν(AB)=0,\displaystyle R_{\mu\nu}-\frac{1}{2}Rg_{\mu\nu}-\frac{3}{L^{2}}g_{\mu\nu}-\frac{1}{2}T^{(A)}_{\mu\nu}-\frac{1}{2}T^{(B)}_{\mu\nu}-\theta T^{(AB)}_{\mu\nu}=0\,, (2a)
μFμν=0,\displaystyle\nabla_{\mu}F^{\mu\nu}=0\,, (2b)
μGμν=0.\displaystyle\nabla_{\mu}G^{\mu\nu}=0\,. (2c)

The last three terms in Einstein equation (2a) are defined as

Tμν(A)=FμρFνρ14gμνF2,\displaystyle T^{(A)}_{\mu\nu}=F_{\mu\rho}F_{\nu}^{\ \rho}-\frac{1}{4}g_{\mu\nu}F^{2}\,,\ (3a)
Tμν(B)=GμρGνρ14gμνG2,\displaystyle T^{(B)}_{\mu\nu}=G_{\mu\rho}G_{\nu}^{\ \rho}-\frac{1}{4}g_{\mu\nu}G^{2}\,,\ (3b)
Tμν(AB)=F(μ|ρ|Gν)ρ14gμνFαβGαβ,\displaystyle T^{(AB)}_{\mu\nu}=F_{(\mu|\rho|}G_{\nu)}^{\ \rho}-\frac{1}{4}g_{\mu\nu}F_{\alpha\beta}G^{\alpha\beta}\,, (3c)

where the symmetry bracket means A(μν)=(Aμν+Aνμ)/2A_{(\mu\nu)}=(A_{\mu\nu}+A_{\nu\mu})/2.

One can obtain the following double charged RN-AdS black blane from the theory (II)

ds2=1u2(f(u)dt2+du2f(u)+dx2+dy2),\displaystyle ds^{2}=\frac{1}{u^{2}}\Big{(}-f(u)dt^{2}+\frac{du^{2}}{f(u)}+dx^{2}+dy^{2}\Big{)}\,,\ (4a)
f(u)=1(1+u+4(qA2+qB2+2θqAqB)4)(uu+)3+u+4(qA2+qB2+2θqAqB)4(uu+)4,\displaystyle f(u)=1-\Big{(}1+\frac{u^{4}_{+}(q_{A}^{2}+q_{B}^{2}+2\theta q_{A}q_{B})}{4}\Big{)}\Big{(}\frac{u}{u_{+}}\Big{)}^{3}+\frac{u^{4}_{+}(q_{A}^{2}+q_{B}^{2}+2\theta q_{A}q_{B})}{4}\Big{(}\frac{u}{u_{+}}\Big{)}^{4}\,,\ (4b)
At(u)=μqAu,Bt(u)=δμqBu,\displaystyle A_{t}(u)=\mu-q_{A}\,u\,,~{}~{}~{}~{}~{}~{}B_{t}(u)=\delta\mu-q_{B}\,u\,, (4c)

where the AdS boundary and the horizon are located at u=0u=0 and u=u+u=u_{+}, respectively. μ\mu and δμ\delta\mu are the chemical potentials of the gauge fields AA and BB in the dual boundary field theory. qAq_{A} and qBq_{B} are two integration constants relating to the chemical potentials by qA=μ/u+q_{A}=\mu/u_{+} and qB=δμ/u+q_{B}=\delta\mu/u_{+}, which are determined by the horizon conditions AtA_{t} and BtB_{t} satisfied. Therefore, there are two controllable chemical potentials causing the unbalance of numbers. The ration χ=δμμ\chi=\frac{\delta\mu}{\mu} stands for the strength of the unbalance. In Kiritsis:2015hoa ; Baggioli:2015dwa ; Huang:2020iyw , it was used to simulate the doping. The coupling strength θ\theta characters the charged impurities coupling strength in the dual field theory. When θ=0\theta=0, the black brane solution (4) reduces to that studied in Bigazzi:2011ak .

According to the solution (4), the Hawking temperature can be straightforward calculated as

T=f(u+)4π=14π(3u+u+(μ2+δμ2+2θμδμ)4).\displaystyle T=-\frac{f^{\prime}(u_{+})}{4\pi}=\frac{1}{4\pi}\Big{(}\frac{3}{u_{+}}-\frac{u_{+}(\mu^{2}+\delta\mu^{2}+2\theta\mu\,\delta\mu)}{4}\Big{)}\,. (5)

It corresponds to the temperature of the dual field theory.

To obtain the other thermodynamical quantities, we write down the renormalized action following the strategy in Caldarelli:2016nni

Sren=S+u=εd3xγ(2K4LLR[γ]).\displaystyle S_{ren}=S+\int_{u=\varepsilon}d^{3}x\sqrt{-\gamma}\,\Big{(}2K-\frac{4}{L}-LR[\gamma]\Big{)}\,. (6)

γ\gamma is the determinant of the boundary induced metric γμν=gμνnμnν\gamma_{\mu\nu}=g_{\mu\nu}-n_{\mu}n_{\nu} and KK is the trace of the extrinsic curvature Kij=12LuuγijK_{ij}=-\frac{1}{2L}u\partial_{u}\gamma_{ij}. They take value at the UV cut-off u=εu=\varepsilon and then is sent to zero following the holographic renormalized procedure. Notice that nμn_{\mu} is the out-point normal vector of the UV cut-off surface.

And then, one obtains the corresponding on-shell renormalized action, which reads

Srenos=d3x[1u+3+μ2+δμ2+2θμδμ4u+].\displaystyle S^{os}_{ren}=\int d^{3}x\Big{[}\frac{1}{u_{+}^{3}}+\frac{\mu^{2}+\delta\mu^{2}+2\theta\mu\,\delta\mu}{4u_{+}}\Big{]}\,. (7)

Immediately, according to holography, the thermal potential Ω\Omega is worked out as

Ω=V2(1u+3+μ2+δμ2+2θμδμ4u+),\displaystyle\Omega=-V_{2}\Big{(}\frac{1}{u_{+}^{3}}+\frac{\mu^{2}+\delta\mu^{2}+2\theta\mu\,\delta\mu}{4u_{+}}\Big{)}\,, (8)

where V2𝑑x𝑑yV_{2}\equiv\int dxdy. Once we have the thermal potential, it is easy to calculate the entropy density ss, charge densities ρA\rho_{A} and ρB\rho_{B}, which are

s=1V2ΩT=4πu+2,ρA=1V2Ωμ=μ+θδμu+,ρB=1V2Ωδμ=θμ+δμu+.\displaystyle s=-\frac{1}{V_{2}}\frac{\partial\Omega}{\partial T}=\frac{4\pi}{u_{+}^{2}}\,,\quad\rho_{A}=-\frac{1}{V_{2}}\frac{\partial\Omega}{\partial\mu}=\frac{\mu+\theta\delta\mu}{u_{+}}\,,\quad\rho_{B}=-\frac{1}{V_{2}}\frac{\partial\Omega}{\partial\,\delta\mu}=\frac{\theta\mu+\delta\mu}{u_{+}}\,. (9)

Also, the press and the energy density of the system can be given by

p=(1u+3+μ2+δμ2+2θμδμ4u+),\displaystyle p=\Big{(}\frac{1}{u_{+}^{3}}+\frac{\mu^{2}+\delta\mu^{2}+2\theta\mu\,\delta\mu}{4u_{+}}\Big{)}\,, (10)
ϵ=2(1u+3+μ2+δμ2+2θμδμ4u+).\displaystyle\epsilon=2\Big{(}\frac{1}{u_{+}^{3}}+\frac{\mu^{2}+\delta\mu^{2}+2\theta\mu\,\delta\mu}{4u_{+}}\Big{)}\,. (11)

Note that the positive definiteness of charge densities ρA\rho_{A} and ρB\rho_{B} requires

{θχ,χ1θ1χ,χ1.\displaystyle\begin{cases}\theta\geq-\chi\,,&\chi\leq 1\cr\theta\geq-\frac{1}{\chi}\,,&\chi\geq 1\end{cases}\,. (12)

III Holographic expression of the conductivities

In Bigazzi:2011ak , the holographic expression of the conductivities have been derived for the two-currents model without coupling. In the presence of a coupling between the two gauge fields, the case becomes subtle and we present the detailed derivation for the conductivities in this section.

Thanks to the rotational invariance in xyx-y plane, we only need consider the conductivities along xx-direction. We first describe the conductivity matrix of the holographic two-currents model. We denote the currents, external fields and conductivities associated to the gauge fields AA and BB as (JA,EA,σA)(J_{A},E_{A},\sigma_{A}) and (JB,EB,σB)(J_{B},E_{B},\sigma_{B}), respectively. At the same time, the external field EAE_{A} also leads to the current JBJ_{B}. The associated conductivity is denoted as γ¯AB\bar{\gamma}_{AB}. This process is reciprocal. EBE_{B} also generates the current JAJ_{A} giving the associated conductivity γ¯BA\bar{\gamma}_{BA}. The time-reversal invariance results in γ¯AB=γ¯BAγ¯\bar{\gamma}_{AB}=\bar{\gamma}_{BA}\equiv\bar{\gamma}. In Refs.Mott:1936v1 ; Mott:1936v2 ; Fert:1968 ; Son:1987 ; Johnson:1987 ; Bigazzi:2011ak , σA\sigma_{A} and σB\sigma_{B} are interpreted as electric conductivity and spin-spin conductivity, and correspondingly γ¯\bar{\gamma} is the spin conductivity. Both currents also induce some momentum, acting on the momentum operator TtxT_{tx} as source. Besides, the temperature gradient leads to the heat current Q=TtxμJAδμJBQ=T_{tx}-\mu J_{A}-\delta\mu J_{B}, which induces thermal conductivity κ¯\bar{\kappa}. The another two transport quantities are the thermo-electric and thermo-spin conductivities associated to the transport of heat. We denote them as α\alpha and β\beta. And then, Ohm’s law can be expressed as

(JAQJB)=(σAαTγ¯αTκ¯TβTγ¯βTσB)(EAT/TEB).\displaystyle\left(\begin{array}[]{c}J_{A}\\ Q\\ J_{B}\\ \end{array}\right)=\left(\begin{array}[]{ccc}\sigma_{A}&\alpha T&\bar{\gamma}\\ \alpha T&\bar{\kappa}T&\beta T\\ \bar{\gamma}&\beta T&\sigma_{B}\\ \end{array}\right)\left(\begin{array}[]{c}E_{A}\\ -\nabla T/T\\ E_{B}\\ \end{array}\right)\,. (22)

The conductivity matrix is symmetric because of the time-reversal symmetry.

Now, we turn to derive the expressions of the conductivities in our holographic two-currents model. To this end, we turn on the bulk fluctuations AxA_{x}, BxB_{x} and gtxg_{tx}, which provide the source for the currents JAxJ^{x}_{A} and JBxJ^{x}_{B}, the stress energy tensor component TtxT^{tx} in the dual boundary field theory. Explicitly, we set

ds2=1u2(f(u)dt2+du2f(u)+δijdxidxj+2gtx(u,t)dtdx),\displaystyle ds^{2}=\frac{1}{u^{2}}\Big{(}-f(u)dt^{2}+\frac{du^{2}}{f(u)}+\delta_{ij}dx^{i}dx^{j}+2g_{tx}(u,t)dtdx\Big{)}\,,\ (23a)
A=At(u)dt+Ax(u,t)dx,B=Bt(u)dt+Bx(u,t)dx.\displaystyle A=A_{t}(u)dt+A_{x}(u,t)dx\,,\quad B=B_{t}(u)dt+B_{x}(u,t)dx\,. (23b)

And then, taking a simple time dependence for the fluctuations as

gtx(u,t)1u2htx(u)eiωt,Ax(u,t)Ax(u)eiωt,Bx(u,t)Bx(u)eiωt,\displaystyle g_{tx}(u,t)\equiv\frac{1}{u^{2}}h_{tx}(u)e^{-i\omega t}\,,\quad A_{x}(u,t)\equiv A_{x}(u)e^{-i\omega t}\,,\quad B_{x}(u,t)\equiv B_{x}(u)e^{-i\omega t}\,, (24)

one obtains the linear EOMs for htx,Ax,Bxh_{tx},A_{x},B_{x} as 333One can check that the two non-vanished Einstein equations along txtx and zxzx directions are equivalent, thus we just list one of them here.

htx=u2((qA+θqB)Ax+(qB+θqA)Bx),\displaystyle h^{\prime}_{tx}=u^{2}\Big{(}(q_{A}+\theta q_{B})A_{x}+(q_{B}+\theta q_{A})B_{x}\Big{)}\,, (25a)
(fAx)+ω2AxfqAhtx=0,\displaystyle(fA^{\prime}_{x})^{\prime}+\frac{\omega^{2}A_{x}}{f}-q_{A}h^{\prime}_{tx}=0\,, (25b)
(fBx)+ω2BxfqBhtx=0.\displaystyle(fB^{\prime}_{x})^{\prime}+\frac{\omega^{2}B_{x}}{f}-q_{B}h^{\prime}_{tx}=0\,. (25c)

It is easy to see that among the above three EOMs, only two of them are independent. In the limit of u0u\rightarrow 0, the fields follow

Ax=Ax(0)+Ax(1)u+,\displaystyle A_{x}=A_{x}^{(0)}+A_{x}^{(1)}u+\cdots\,,\ (26a)
Bx=Bx(0)+Bx(1)u+,\displaystyle B_{x}=B_{x}^{(0)}+B_{x}^{(1)}u+\cdots\,,\ (26b)
htx=htx(0)+htx(1)u3+.\displaystyle h_{tx}=h_{tx}^{(0)}+h_{tx}^{(1)}u^{3}+\cdots\,. (26c)

To have a well-defined bulk variational problem, we write down the following renormalized action

Sren=S+u=εd3xγ(2K4).\displaystyle S_{ren}=S+\int_{u=\varepsilon}d^{3}x\sqrt{-\gamma}\,(2K-4)\,. (27)

Making the variation of the on-shell action, one has

δSrenonshell\displaystyle\delta S^{on-shell}_{ren} =u=εd3x|γ|(Kμν+Kγμν2γμν)δγμνu=εd3x|γ|nμ(Fμν+θGμν)δAν\displaystyle=\int_{u=\varepsilon}d^{3}x\sqrt{|\gamma|}\Big{(}-K^{\mu\nu}+K\gamma^{\mu\nu}-2\gamma^{\mu\nu}\Big{)}\delta\gamma_{\mu\nu}-\int_{u=\varepsilon}d^{3}x\sqrt{|\gamma|}n_{\mu}(F^{\mu\nu}+\theta G^{\mu\nu})\delta A_{\nu}
u=εd3x|γ|nμ(Gμν+θFμν)δBν.\displaystyle-\int_{u=\varepsilon}d^{3}x\sqrt{|\gamma|}n_{\mu}(G^{\mu\nu}+\theta F^{\mu\nu})\delta B_{\nu}\,. (28)

Further using the ansatz (23) and the UV expansion (26), we can evaluate the above equation as

δSrenonshell=u=εd3x[2u3(11f)htx12u2htx]eiωtδhtxeiωt\displaystyle\delta S^{on-shell}_{ren}=\int_{u=\varepsilon}d^{3}x\Big{[}\frac{2}{u^{3}}(1-\frac{1}{\sqrt{f}})h_{tx}-\frac{1}{2u^{2}}h^{\prime}_{tx}\Big{]}e^{-i\omega t}\delta h_{tx}e^{-i\omega t}
+u=εd3x[(f(Ax+θBx)(qA+θqB)htx)eiωtδAxeiωt\displaystyle+\int_{u=\varepsilon}d^{3}x\Big{[}\big{(}f(A^{\prime}_{x}+\theta B^{\prime}_{x})-(q_{A}+\theta q_{B})h_{tx}\big{)}e^{-i\omega t}\delta A_{x}e^{-i\omega t}
+(f(θAx+Bx)(θqA+qB)htx)eiωtδBxeiωt]\displaystyle+\big{(}f(\theta A^{\prime}_{x}+B^{\prime}_{x})-(\theta q_{A}+q_{B})h_{tx}\big{)}e^{-i\omega t}\delta B_{x}e^{-i\omega t}\Big{]}
=u=εd3x[2u3(11f)htx(0)12((qA+θqB)Ax(0)+(θqA+qB)Bx(0))]eiωtδhtx(0)eiωt\displaystyle=\int_{u=\varepsilon}d^{3}x\Big{[}\frac{2}{u^{3}}(1-\frac{1}{\sqrt{f}})h^{(0)}_{tx}-\frac{1}{2}\big{(}(q_{A}+\theta q_{B})A^{(0)}_{x}+(\theta q_{A}+q_{B})B^{(0)}_{x}\big{)}\Big{]}e^{-i\omega t}\delta h^{(0)}_{tx}e^{-i\omega t}
+u=εd3x[(Ax(1)+θBx(1)(qA+θqB)htx(0))eiωtδAx(0)eiωt\displaystyle+\int_{u=\varepsilon}d^{3}x\Big{[}\big{(}A^{(1)}_{x}+\theta B^{(1)}_{x}-(q_{A}+\theta q_{B})h^{(0)}_{tx}\big{)}e^{-i\omega t}\delta A^{(0)}_{x}e^{-i\omega t}
+(θAx(1)+Bx(1)(θqA+qB)htx(0))eiωtδBx(0)eiωt].\displaystyle+\big{(}\theta A^{(1)}_{x}+B^{(1)}_{x}-(\theta q_{A}+q_{B})h^{(0)}_{tx}\big{)}e^{-i\omega t}\delta B^{(0)}_{x}e^{-i\omega t}\Big{]}\,. (29)

According to the holographic dictionary

Tμν(t)=2δSrenonshell[gμν(0)]δgμν(0)(t),Jμ(t)=δSrenonshell[Aμ(0)]δAμ(0)(t),\displaystyle\big{\langle}T^{\mu\nu}\big{\rangle}(t)=2\frac{\delta S_{ren}^{on-shell}[g_{\mu\nu}^{(0)}]}{\delta g_{\mu\nu}^{(0)}(t)}\,,\qquad\big{\langle}J^{\mu}\big{\rangle}(t)=\frac{\delta S_{ren}^{on-shell}[A_{\mu}^{(0)}]}{\delta A_{\mu}^{(0)}(t)}\,, (30)

one obtains the expectation values of JAx,JBxJ^{x}_{A},\,J^{x}_{B} and TtxT^{tx} as

Ttx(ω)=4u3(11f)htx(0)((qA+θqB)Ax(0)+(θqA+qB)Bx(0)),\displaystyle\big{\langle}T^{tx}\big{\rangle}(\omega)=\frac{4}{u^{3}}(1-\frac{1}{\sqrt{f}})h^{(0)}_{tx}-\Big{(}(q_{A}+\theta q_{B})A^{(0)}_{x}+(\theta q_{A}+q_{B})B^{(0)}_{x}\Big{)}\,, (31a)
JAx(ω)=Ax(1)+θBx(1)(qA+θqB)htx(0),\displaystyle\big{\langle}J^{x}_{A}\big{\rangle}(\omega)=A^{(1)}_{x}+\theta B^{(1)}_{x}-(q_{A}+\theta q_{B})h^{(0)}_{tx}\,, (31b)
JBx(ω)=θAx(1)+Bx(1)(θqA+qB)htx(0).\displaystyle\big{\langle}J^{x}_{B}\big{\rangle}(\omega)=\theta A^{(1)}_{x}+B^{(1)}_{x}-(\theta q_{A}+q_{B})h^{(0)}_{tx}\,. (31c)

This source-response relation can be wrote in the matrix form,

(JAx(ω)Ttx(ω)JBx(ω))=(δAx(1)δAx(0)+θδBx(1)δAx(0),(qA+θqB),δAx(1)δBx(0)+θδBx(1)δBx(0)(qA+θqB),ϵ,(θqA+qB)θδAx(1)δAx(0)+δBx(1)δAx(0),(θqA+qB),θδAx(1)δBx(0)+δBx(1)δBx(0))(Ax(0)(ω)htx(0)(ω)Bx(0)(ω)),\displaystyle\left(\begin{array}[]{c}J^{x}_{A}(\omega)\\ T^{tx}(\omega)\\ J^{x}_{B}(\omega)\\ \end{array}\right)=\left(\begin{array}[]{ccc}\dfrac{\delta A^{(1)}_{x}}{\delta A^{(0)}_{x}}+\theta\dfrac{\delta B^{(1)}_{x}}{\delta A^{(0)}_{x}},&-(q_{A}+\theta q_{B}),&\dfrac{\delta A^{(1)}_{x}}{\delta B^{(0)}_{x}}+\theta\dfrac{\delta B^{(1)}_{x}}{\delta B^{(0)}_{x}}\\ -(q_{A}+\theta q_{B}),&-\epsilon,&-(\theta q_{A}+q_{B})\\ \theta\dfrac{\delta A^{(1)}_{x}}{\delta A^{(0)}_{x}}+\dfrac{\delta B^{(1)}_{x}}{\delta A^{(0)}_{x}},&-(\theta q_{A}+q_{B}),&\theta\dfrac{\delta A^{(1)}_{x}}{\delta B^{(0)}_{x}}+\dfrac{\delta B^{(1)}_{x}}{\delta B^{(0)}_{x}}\\ \end{array}\right)\left(\begin{array}[]{c}A^{(0)}_{x}(\omega)\\ h_{tx}^{(0)}(\omega)\\ B^{(0)}_{x}(\omega)\\ \end{array}\right)\,, (41)

where the energy density is ϵ=limu04u3(1f1/2)\epsilon=-lim_{u\rightarrow 0}4u^{-3}(1-f^{-1/2}). And then, we have the relation of the heat current QxQ^{x}, electric fields EAx,EBxE_{Ax},E_{Bx}, and temperature gradient xT\nabla_{x}T on the source, which are given by

(JAxQtxJBx)=(100μ1δμ001)(JAxTtxJBx),\displaystyle\left(\begin{array}[]{c}J^{x}_{A}\\ Q^{tx}\\ J^{x}_{B}\\ \end{array}\right)=\left(\begin{array}[]{ccc}1&0&0\\ -\mu&1&-\delta\mu\\ 0&0&1\\ \end{array}\right)\left(\begin{array}[]{c}J^{x}_{A}\\ T^{tx}\\ J^{x}_{B}\\ \end{array}\right)\,,\ (42j)
(Ax(0)htx(0)Bx(0))=1iω(1μ00100δμ1)(EAxxT/TEBx).\displaystyle\left(\begin{array}[]{c}A^{(0)}_{x}\\ h_{tx}^{(0)}\\ B^{(0)}_{x}\\ \end{array}\right)=\frac{1}{i\omega}\left(\begin{array}[]{ccc}1&-\mu&0\\ 0&1&0\\ 0&-\delta\mu&1\\ \end{array}\right)\left(\begin{array}[]{c}E_{Ax}\\ -\nabla_{x}T/T\\ E_{Bx}\\ \end{array}\right)\,. (42t)

Together with (41) and comparing with Ohm’s law (22), one obtains the holographic expressions of alternating current (AC) conductivities of the two-currents model

σA=1iω(δAx(1)δAx(0)+θδBx(1)δAx(0)),\displaystyle\sigma_{A}=\frac{1}{i\omega}\Big{(}\dfrac{\delta A^{(1)}_{x}}{\delta A^{(0)}_{x}}+\theta\dfrac{\delta B^{(1)}_{x}}{\delta A^{(0)}_{x}}\Big{)}\,,\ (43a)
σB=1iω(θδAx(1)δBx(0)+δBx(1)δBx(0)),\displaystyle\sigma_{B}=\frac{1}{i\omega}\Big{(}\theta\dfrac{\delta A^{(1)}_{x}}{\delta B^{(0)}_{x}}+\dfrac{\delta B^{(1)}_{x}}{\delta B^{(0)}_{x}}\Big{)}\,,\ (43b)
γ¯=1iω(δAx(1)δBx(0)+θδBx(1)δBx(0))=1iω(θδAx(1)δAx(0)+δBx(1)δAx(0)),\displaystyle\bar{\gamma}=\frac{1}{i\omega}\Big{(}\dfrac{\delta A^{(1)}_{x}}{\delta B^{(0)}_{x}}+\theta\dfrac{\delta B^{(1)}_{x}}{\delta B^{(0)}_{x}}\Big{)}=\frac{1}{i\omega}\Big{(}\theta\dfrac{\delta A^{(1)}_{x}}{\delta A^{(0)}_{x}}+\dfrac{\delta B^{(1)}_{x}}{\delta A^{(0)}_{x}}\Big{)}\,,\ (43c)
αT=μσAδμγqA+θqBiω,\displaystyle\alpha T=-\mu\sigma_{A}-\delta\mu\,\gamma-\frac{q_{A}+\theta q_{B}}{i\omega}\,,\ (43d)
βT=δμσBμγθqA+qBiω,\displaystyle\beta T=-\delta\mu\sigma_{B}-\mu\,\gamma-\frac{\theta q_{A}+q_{B}}{i\omega}\,,\ (43e)
κ¯T=1iω(ϵp+2μ(qA+θqB)+2δμ(θqA+qB))+μ2σA+(δμ)2σB+2μδμγ.\displaystyle\bar{\kappa}T=\frac{1}{i\omega}\Big{(}-\epsilon-p+2\mu(q_{A}+\theta q_{B})+2\delta\mu(\theta q_{A}+q_{B})\Big{)}+\mu^{2}\sigma_{A}+(\delta\mu)^{2}\sigma_{B}+2\mu\delta\mu\,\gamma\,. (43f)

Since the two gauge fields are directly coupled in the gravity action, the conductivities σA\sigma_{A} and σB\sigma_{B} are also directly related to the coupling strength θ\theta. As a result, the other conductivities γ¯\bar{\gamma}, α\alpha and β\beta are also related to the coupling θ\theta.

Until now, we have worked out the expressions of AC conductivities of the two-currents model in the presence of coupling, which is the main topic of our present paper. Using the expressions, we can explore the transport properties of holographic two-currents model with couple, for example, the superconductivity. As a simple application, here we calculate the AC conductivities of our present model and briefly discuss its properties.

IV Numerical Results

By numerically solving the EOMs (25), we can study the transport properties. In the numerical calculation, by rescaling, the horizon location can be set as unity, i.e., u+=1u_{+}=1. Thanks to the scaling symmetry, we take the chemical potential μ\mu of gauge field AA as scaling unite. Therefore, our theory is specified by the two dimensionless parameters T^Tμ\hat{T}\equiv\frac{T}{\mu} as well as χ\chi, and the coupling parameter θ\theta.

Refer to caption
Refer to caption
Figure 1: The real part of the conductivities σA\sigma_{A} (left) and σB\sigma_{B} (right) at T^=0.05\hat{T}=0.05 for χ=1\chi=1 but for different parameter θ\theta.
Refer to caption
Refer to caption
Figure 2: The real part of the conductivities αT\alpha T (left) and βT\beta T (right) at T^=0.05\hat{T}=0.05 for χ=1\chi=1 but for different parameter θ\theta.

This holographic system possesses the following symmetries

σA(μ,δμ,θ,ω^)=σB(δμ,μ,θ,ω^),\displaystyle\sigma_{A}(\mu,\delta\mu,\theta,\hat{\omega})=\sigma_{B}(\delta\mu,\mu,\theta,\hat{\omega})\,,\ (44a)
αT(μ,δμ,θ,ω^)=βT(δμ,μ,θ,ω^),\displaystyle\alpha T(\mu,\delta\mu,\theta,\hat{\omega})=\beta T(\delta\mu,\mu,\theta,\hat{\omega})\,, (44b)

which can easily deduced from the EOMs (2) and the expressions of the conductivity (43). Numerically, we have also confirmed the above result. It is similar to that without coupling θ\theta Bigazzi:2011ak . Especially, for μ=δμ\mu=\delta\mu, i.e., χ=1\chi=1, one has σA=σB\sigma_{A}=\sigma_{B} and αT=βT\alpha T=\beta T, which are clearly seen in FIG.1 and FIG.2.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 3: The real part of the conductivities σA\sigma_{A} (left) and σB\sigma_{B} (right) at T^=0.05\hat{T}=0.05 for different parameter θ\theta (the plots above for χ=0.5\chi=0.5 and the ones below for χ=1.2\chi=1.2).

Next, we explore the properties of σA\sigma_{A} and σB\sigma_{B}. When θ=0\theta=0, which has been studied in Bigazzi:2011ak , the real part of the conductivities σA\sigma_{A}/σB\sigma_{B} at low frequency exhibits a dip (σB\sigma_{B} is frequency independent when χ=0\chi=0.). With the doping parameter χ\chi increases, the dip of Re[σA]Re[\sigma_{A}] becomes shallow but the one of Re[σB]Re[\sigma_{B}] becomes deepening (see FIG.1 and FIG.3). When χ\chi is fixed, the dip in Re[σA]Re[\sigma_{A}]/Re[σB]Re[\sigma_{B}] becomes more and more deepening as θ\theta increases and finally turns into a hard-gap-like when θ=1\theta=1 is achieved, which is independent of the doping χ\chi. If we further tune θ\theta larger such that it is beyond the unity, the DC conductivities of Re[σA]Re[\sigma_{A}]/Re[σB]Re[\sigma_{B}] will be negative, which violate the positive definiteness of the conductivity. Therefore, the positive definiteness of the conductivity imposes a constraint on the coupling parameter θ\theta. Here, we shall constrain θ\theta in the range of θ1\theta\leq 1. Some higher derivative coupling terms also lead to the violation of the positive definiteness of the conductivity Witczak-Krempa:2013aea ; Fu:2017oqa ; Gouteraux:2016wxj ; Baggioli:2016pia .

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 4: The real part of γ¯\bar{\gamma} and of the thermal conductivity κ¯\bar{\kappa} for sample χ\chi and θ\theta (T^=0.05\hat{T}=0.05).

Further, we report the real part of γ¯\bar{\gamma} and of the thermal conductivity κ¯\bar{\kappa} for sample χ\chi and θ\theta in FIG.4. We find that Re[γ¯]Re[\bar{\gamma}] exhibits a dip at low frequency and approaches a constant at high frequency. At full frequency, Re[γ¯]Re[\bar{\gamma}] increases as θ\theta increases. While Re[κ¯T]Re[\bar{\kappa}T] converges to zero in the limit of ω^0\hat{\omega}\rightarrow 0 independently of the doping χ\chi and the coupling θ\theta, which means that the DC thermal conductivity vanishes. As the frequency increases, the thermal conductivities with different χ\chi and θ\theta increase and separate out, and then approaches different constant value depending on χ\chi and θ\theta. With the increase of χ\chi or θ\theta, this constant value increases.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 5: The real part of the conductivities αT\alpha T (left) and βT\beta T (right) for χ=0.5\chi=0.5 (above) and χ=1.2\chi=1.2 (bottom) at T^=0.05\hat{T}=0.05.

Finally, we also show the real part of the conductivities αT\alpha T and βT\beta T for sample χ\chi and θ\theta in FIG.5 (also see FIG.2 for χ=1\chi=1 and different θ\theta). We see that both αT\alpha T and βT\beta T are negative and at high frequency, they monotonously decreases with the increase of θ\theta.

V Conclusions and discussions

In this paper, we construct a holographic gravity model of two gauge fields with a coupling between them, which corresponds to a two-currents model. An analytical black brane solution is obtained. Also we briefly discuss the thermodynamics. When this coupling is introduced, the expressions of conductivities for holographic two-currents model without coupling studied in Bigazzi:2011ak are no longer applicable. By the standard holographic renormalized procedure, we work out the expressions of conductivities with coupling (see Eqs.(43)). We find that the expressions of conductivities are directly related to the coupling parameter θ\theta. When θ=0\theta=0, they reduce to that without coupling in Bigazzi:2011ak . The expressions of conductivities for holographic two-currents model with coupling are the main topic of our present paper. Our results are also applicable for the holographic two-currents superconductor model with couple or other extension models.

As an application, then we briefly discuss the properties of the conductivities of this holographic two-currents model with coupling. An interesting property is that as the coupling θ\theta increases, the dip at low frequency in Re[σA]Re[\sigma_{A}]/Re[σB]Re[\sigma_{B}] becomes deepening and finally turns into a hard-gap-like when θ=1\theta=1, which is independent of the doping χ\chi. Some monotonic behaviors of the conductivities are also discussed.

Along this direction, there are lots of works deserving further exploration. For example, we can add a charged complex scalar field to study the superconducting instability and the properties of the conductivities. Also it shall be surely interesting to implement the momentum dissipation into our system with coupling and study the properties of the conductivities.

Acknowledgements.
This work is supported by the Natural Science Foundation of China under Grant Nos. 11775036, 11905182 and Fok Ying Tung Education Foundation under Grant No. 171006. Guoyang Fu is supported by the Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX20_2973). J. P. Wu is also supported by Top Talent Support Program from Yangzhou University.

References

  • (1) J. M. Maldacena, “The Large N limit of superconformal field theories and supergravity,” Adv. Theor. Math. Phys.  2, 231 (1998) [hep-th/9711200].
  • (2) S. S. Gubser, I. R. Klebanov and A. M. Polyakov, “Gauge theory correlators from noncritical string theory,” Phys. Lett. B 428, 105 (1998) [hep-th/9802109].
  • (3) E. Witten, “Anti-de Sitter space and holography,” Adv. Theor. Math. Phys.  2, 253 (1998) [hep-th/9802150].
  • (4) O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, “Large N field theories, string theory and gravity,” Phys. Rept.  323, 183 (2000) [hep-th/9905111].
  • (5) S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, “Building a Holographic Superconductor,” Phys. Rev. Lett.  101 (2008) 031601. [arXiv:0803.3295 [hep-th]].
  • (6) A. Donos and S. A. Hartnoll, “Interaction-driven localization in holography,” Nature Phys. 9 (2013), 649-655 [arXiv:1212.2998 [hep-th]].
  • (7) Y. Ling, C. Niu, J. Wu, Z. Xian and H. b. Zhang, “Metal-insulator Transition by Holographic Charge Density Waves,” Phys. Rev. Lett. 113, 091602 (2014) [arXiv:1404.0777 [hep-th]].
  • (8) H. Liu, J. McGreevy and D. Vegh, “Non-Fermi liquids from holography,” Phys. Rev. D 83, 065029 (2011) [arXiv:0903.2477 [hep-th]].
  • (9) E. Kiritsis and L. Li, “Holographic Competition of Phases and Superconductivity,” JHEP 1601, 147 (2016) [arXiv:1510.00020 [cond-mat.str-el]].
  • (10) M. Baggioli and M. Goykhman, “Under The Dome: Doped holographic superconductors with broken translational symmetry,” JHEP 01 (2016), 011 [arXiv:1510.06363 [hep-th]].
  • (11) W. Huang, G. Fu, D. Zhang, Z. Zhou and J. P. Wu, “Doped holographic fermionic system,” Eur. Phys. J. C 80 (2020), 608 [arXiv:2002.03343 [hep-th]].
  • (12) F. Bigazzi, A. L. Cotrone, D. Musso, N. Pinzani Fokeeva and D. Seminara, “Unbalanced Holographic Superconductors and Spintronics,” JHEP 1202, 078 (2012) [arXiv:1111.6601 [hep-th]].
  • (13) N. Iqbal, H. Liu, M. Mezei and Q. Si, “Quantum phase transitions in holographic models of magnetism and superconductors,” Phys. Rev. D 82 (2010), 045002 [arXiv:1003.0010 [hep-th]].
  • (14) M. Rogatko and K. I. Wysokinski, “Two interacting current model of holographic Dirac fluid in graphene,” Phys. Rev. D 97 (2018) no.2, 024053 [arXiv:1708.08051 [hep-th]].
  • (15) M. Rogatko and K. I. Wysokinski, “Conductivity bound of the strongly interacting and disordered graphene from gauge/gravity duality,” Phys. Rev. D 101 (2020) no.4, 046019 [arXiv:2002.02177 [hep-th]].
  • (16) Y. Seo, G. Song, P. Kim, S. Sachdev and S. J. Sin, “Holography of the Dirac Fluid in Graphene with two currents,” Phys. Rev. Lett. 118 (2017) no.3, 036601 [arXiv:1609.03582 [hep-th]].
  • (17) N. F. Mott, “The electrical Conductivity of Transition Metals,” Proc. R. Soc. Lond. A 153, 699 (1936).
  • (18) N. F. Mott, “The Resistance and Thermoelectric Properties of the Transition Metals,” Proc. R. Soc. Lond. A 156, 368 (1936).
  • (19) Y. Ling, P. Liu and J. P. Wu, “A novel insulator by holographic Q-lattices,” JHEP 02 (2016), 075 [arXiv:1510.05456 [hep-th]].
  • (20) Y. Ling, P. Liu, J. P. Wu and M. H. Wu, “Holographic superconductor on a novel insulator,” Chin. Phys. C 42 (2018) no.1, 013106 [arXiv:1711.07720 [hep-th]].
  • (21) Y. Ling and M. H. Wu, “The instability of AdSAdS black holes with lattices,” [arXiv:2009.00510 [hep-th]].
  • (22) Y. Ling, P. Liu and J. P. Wu, “Characterization of Quantum Phase Transition using Holographic Entanglement Entropy,” Phys. Rev. D 93 (2016) no.12, 126004 [arXiv:1604.04857 [hep-th]].
  • (23) J. Tarrio and S. Vandoren, “Black holes and black branes in Lifshitz spacetimes,” JHEP 1109, 017 (2011) [arXiv:1105.6335 [hep-th]].
  • (24) M. Alishahiha, E. O Colgain and H. Yavartanoo, “Charged Black Branes with Hyperscaling Violating Factor,” JHEP 11 (2012), 137 [arXiv:1209.3946 [hep-th]].
  • (25) S. S. Gubser, F. D. Rocha, “Peculiar properties of a charged dilatonic black hole in AdS5AdS_{5}”, Phys. Rev. D 81, 046001 (2010), [arXiv:0911.2898 [hep-th]].
  • (26) J. Erdmenger, V. Grass, P. Kerner and T. H. Ngo, “Holographic Superfluidity in Imbalanced Mixtures,” JHEP 08 (2011), 037 [arXiv:1103.4145 [hep-th]].
  • (27) A. Dutta and S. K. Modak, “Holographic entanglement entropy in imbalanced superconductors,” JHEP 01 (2014), 136 [arXiv:1305.6740 [hep-th]].
  • (28) D. Correa, N. Grandi and A. Hernandez, “Doped Holographic Superconductor in an External Magnetic Field,” JHEP 19 (2020), 085 [arXiv:1905.05132 [hep-th]].
  • (29) D. Musso, “Competition/Enhancement of Two Probe Order Parameters in the Unbalanced Holographic Superconductor,” JHEP 06 (2013), 083 [arXiv:1302.7205 [hep-th]].
  • (30) J. Alsup, E. Papantonopoulos and G. Siopsis, “A Novel Mechanism to Generate FFLO States in Holographic Superconductors,” Phys. Lett. B 720 (2013), 379-384 [arXiv:1210.1541 [hep-th]].
  • (31) A. J. Hafshejani and S. A. Hosseini Mansoori, “Unbalanced St ckelberg holographic superconductors with backreaction,” JHEP 01 (2019), 015 Z[arXiv:1808.02628 [hep-th]].
  • (32) S. Grozdanov, A. Lucas, S. Sachdev and K. Schalm, “Absence of disorder-driven metal-insulator transitions in simple holographic models,” Phys. Rev. Lett. 115 (2015) no.22, 221601 [arXiv:1507.00003 [hep-th]].
  • (33) M. M. Caldarelli, A. Christodoulou, I. Papadimitriou and K. Skenderis, “Phases of planar AdS black holes with axionic charge,” JHEP 1704, 001 (2017) [arXiv:1612.07214 [hep-th]].
  • (34) A. Fert and I. A. Campbell, “Two-Current Conduction in Nickel,” Phys. Rev. Lett. 21, 1190 (1968).
  • (35) P. C. van Son, H. van Kempen and P. Wyder, “Boundary Resistance of the Ferromagnetic-Nonferromagnetic Metal Interface,” Phys. Rev. Lett. 58, 2271 (1987).
  • (36) M. Johnson and R. H. Silsbee, “Thermodynamic analysis of interfacial transport and of the thermomagnetoelectric system,” Phys. Rev. B 35, 4959 (1987).
  • (37) W. Witczak-Krempa, “Quantum critical charge response from higher derivatives in holography,” Phys. Rev. B 89 (2014) no.16, 161114 [arXiv:1312.3334 [cond-mat.str-el]].
  • (38) G. Fu, J. P. Wu, B. Xu and J. Liu, “Holographic response from higher derivatives with homogeneous disorder,” Phys. Lett. B 769 (2017), 569-574 [arXiv:1705.06672 [hep-th]].
  • (39) B. Gouteraux, E. Kiritsis and W. J. Li, “Effective holographic theories of momentum relaxation and violation of conductivity bound,” JHEP 04 (2016), 122 [arXiv:1602.01067 [hep-th]].
  • (40) M. Baggioli, B. Gouteraux, E. Kiritsis and W. J. Li, “Higher derivative corrections to incoherent metallic transport in holography,” JHEP 03 (2017), 170 [arXiv:1612.05500 [hep-th]].