Holographic entanglement entropy with momentum relaxation

We studied the holographic entanglement entropy for a strip and sharp wedge entangling regions in momentum relaxation systems. In the case of strips, we found analytic and numerical results for the entanglement entropy and examined the effect of the electric field on the entanglement entropy. We also studied the entanglement entropy of wedges and confirmed that there is a linear change due to the electric field. In this case, we showed that the entanglement entropy change is interestingly proportional to the thermoelectric conductivity which is a measurable quantity. Also, we discussed comparable calculations with ABJM theory and we suggested an experiment for our result.


Introduction
Entanglement entropy plays a prominent role in various areas of physics such as quantum field theory, gravity theory and condensed matter physics. However, it is very difficult to evaluate this amount because of non-local property. Fortunately, an effective way to calculate entanglement entropy has recently been developed by Ryu and Takayanagi [1]. Since this method is based on AdS/CFT correspondence [2,3] which relates gravitational systems to quantum field theories, it may shed a light on understanding quantum gravity [4][5][6]. Apart from this, the method is useful to study various interesting cases. For examples, it has been used to distinguish different phases e.g., [7][8][9] and to study thermalization under quantum quenches e.g., [10][11][12][13]. Further studies with the method show us new universal properties of generic CFT e.g., [14][15][16][17][18]. Even for a particular non-conformal field thea e-mail: kimkyungkiu@sejong.ac.kr b e-mail: cyong21@gist.ac.kr c e-mail: junghun.lee@gist.ac.kr d e-mail: bjahn@yonsei.ac.kr ory, it turns out that Ryu and Takayanagi's proposal is still useful by considering a top-down study on the field theory and the corresponding supergravity solution [19][20][21][22][23][24][25].
Entanglement entropy is a very important quantity but, in most cases, it is difficult to calculate. It is also unclear whether this amount can be measured easily. 1 Therefore, it is hard to tell if Ryu and Takayanagi's method is valid in various cases. In this study, we try to find examples where entanglement entropy changes can be represented by other readily measurable physical quantities. Like thermal entropy, the entanglement entropy responses to other external sources. One of such sources which can be produced in a laboratory is the electric field. Thus we will investigate how the entanglement entropy changes under the electric field using holographic methods.
Accordingly, to demonstrate a system in the electric field holographically, we need to consider a gravity dual describing the system. The Ryu-Takayanagi's formula and the gravity dual allow us to study the response of entanglement entropy to the electric field. For the sake of simplicity, we introduce a time-independent electric field, i.e, DC electric field. To realize such a situation, we take the momentum relaxation into account. This is because infinite DC currents can occur without the momentum relaxation. The gravity duals have been developed to study holographic DC conductivities in the context of AdS/CMT e.g.,  and some interesting time-dependent extensions e.g., [47][48][49][50]. Among these backgrounds, we adopted a background geometry of an axion model, where a scalar field breaks translation invariance explicitly [28].
In the absence of the bulk gauge field, an analytic form of holographic entanglement entropy with small momentum relaxation was studied in [51]. In this work, we performed analytic and numerical studies on the strip entangling region in the presence of the bulk electric field, that is dual to the density and the external electric field in the boundary system. We obtained new analytic expressions of the entanglement entropy in the very narrow strip limit. Since our study is about 2+1 dimensional field theory systems, one may compare the result to a field theoretic calculation qualitatively. So, we suggest a comparable situation in the ABJM theory [52]. Also, we provide numerical results without taking any limit. Our numerical results reproduce the earlier work [53] and extend its numerical study. This calculation shows that entanglement entropy becomes larger with more charge density and stronger momentum relaxation. We will discuss the physical implication for the extended parameter space.
Furthermore, we studied how the strip type entanglement entropy changes by turning on the constant electric field. And we found that the minimal surface anchored to the strip is tilted by the electric field and the deformation is proportional to the thermoelectric conductivityᾱ in the conductivity tensor. This is a quite interesting result because the thermoelectric conductivity is a clearly measurable quantity. However, the entanglement entropy changes from the deformation do not appear in the linear order of the electric field. To see the effect, we need to consider the quadratic order gravity dual but it is beyond our present work.
Instead we consider the wedge type entangling region whose symmetric axis is not orthogonal to the electric field . Finally, we took a sharp wedge limit and obtained the entanglement entropy changes at the linear level of the electric field. And we confirmed that the response to the electric field is proportional to the anticipated thermoelectric conductivitȳ α. This result is also remarkable because the thermoelectric conductivity or the Nernst signal could reflect the existence of the quantum critical point even in the normal phase of the superconductor [41,54,55,76]. Thus it is suggested that the entanglement entropy may be associated with the quantum critical point. This paper is organized as follows: In Sect. 2 we introduce the bulk solution dual to a momentum relaxation system with the electric field and summarize how the conductivities can be related to the bulk fields. In Sect. 3, we study on the holographic entanglement entropy for the strip entangling region. In Sect. 4, we consider the holographic entanglement entropy for the wedge region. In Sect. 5. we conclude our work.

A gravity dual of momentum relaxation with electric field
In this section we review a gravity dual to a momentum relaxation system in the electric field. The geometry has been discussed to study finite DC conductivity e.g., [28,40]. The momentum relaxation in physical systems is essential to obtain finite DC conductivity. Let us start with a simple model with momentum relaxation: (1) This system admits a black brane solution as follows: with U (r ) = r 2 − β 2 2 − M r + q 2 4r 2 . 2 The massless scalar field, χ I , is associated with momentum relaxation, breaking the translational symmetry and inducing finite DC conductivity. To see this, let us turn on the external electric field and other related fields along the x-direction where E x corresponds to the external field and other fluctuations are dual to the electric current and the heat current in the boundary system. Then, the linearized equations of motion are given as follows : where Eq. (5) comes from the Einstein equation, and Eq. (6) and (7) are the equations of motion of the matter fields. Note that all the equations above are not independent because Eq. (7) can be obtained by rearranging the other equations. As a consequence, only three equations in (5) and (6) determine the profile of the fluctuations for given boundary conditions. Another thing we need to know is that h r x is not a dynamical field and it plays a role of a Lagrange multiplier. In fact, this is a shift vector in an ADM decomposition along the holographic radial direction. We will discuss the gauge fixing later.
The electric current and the heat current using gauge/ gravity duality are given by where In addition, it can be easily checked that J (r ) and Q(r ) remain as constants along the holographic radial direction r by using the equations of motion (5)(6)(7). Therefore, the boundary currents can be computed at the horizon r = r h Near the horizon the fields behave as where the Hawking temperature T is given by T = U (r h ) 4π and the logarithmic term appears due to the ingoing boundary condition. By solving the equations of motion near the horizon, one can obtain Using these results, the linear response theory yields the following electric conductivity (σ = J x /E x ) and thermoelectric conductivity (ᾱ = Q x /E x ): So the constant H r x is given by a combination of physical parameters as H r x = − 1 4πᾱ E x . The following sections will consider minimum surfaces in this geometry. It is convenient to use a formal expansion parameter λ, which will be taken as 1 after the calculation. This trick allows us to write the metric in the following form: In addition to the conductivities, other physical quantities such as temperature, energy density, charge density, and entropy density can be read from the geometry as follows: These quantities characterize the dual system to the unperturbed geometry with λ = 0.
3 Holographic entanglement entropy in electric field: strip entangling region

Basic setup
In this subsection, we use the Ryu-Takayanagi formula to obtain entanglement entropy for the strip region. The strip also extends along a direction perpendicular to the electric field. To describe the minimal surface anchored to the strip, we may consider a two dimensional surface with the following map in the target geometry (14): where y coordinate indicates the infinitely extending direction. Then, one can find the action for the minimal surface as follows: where L is the length of the strip along y-direction. This is our basic action for the holographic entanglement entropy. Or, if we define another convenient coordinates z ≡ 1 lr and σ ≡ x/l, we can write down the action in terms of z(σ ) as follows: where f (z) = 1 z 2 −β The Lagrangian in (19) has no explicit dependence of σ . So there is a conserved quantity H given by where z * ≡ z(σ * ) is the location of the tip of the minimal surface. Since we are considering regular surfaces, z (σ ) vanishes at the tip of the surface. From this relation one can find the integrated first order equations of motion as follows : where we took a gauge δg r x = H r x r 2 U (r ) , which is consistent with the boundary conditions. 3 Plugging the above equations into (19), we obtained the regularized area for the minimal surface as follows: where is the length scale corresponding to the UV cutoff.
Here, one can notice that the effect from the electric field is encoded in the location of the tip z * . By using the Ryu-Takayanagi formula, the entanglement entropy for the strip is: where G N is the 4-dimensional Newton constant. In order to show the numerical result, we define the following refined function which is 4 G N l L times the finite part of the entanglement entropy up to the minus sign: 3 In the calculation we took a gauge h r x = Hrx r 2 U (r ) . It is legitimate because h r x plays a role of a Lagrange multiplier at the linear level and the choice satisfies the regularity condition at the horizon and near the boundary of AdS space. Furthermore, our final result in this work depends on the gauge invariant quantities, such as the thermoelectric conductivityᾱ and geometric angle δ. For this reason, the gauge fixing is believed to be adequate to consider this calculation. In addition, we speculate that the gauge fixing would be clearer if we consider the second order effect of the electric field on the geometry. This is because additional constraints must be considered for g r x . Now, let us look at the equation of motion for the minimal surface. The equation of motion can be derived from the action (19) and the solution is expanded in terms of λ as follows : Then the equation of motion for each order is given by Here, the zeroth order solution is given by an even function from the symmetry of (26). Furthermore, suppose z 1 (σ ) is a solution of (27), then one can easily show that −z 1 (−σ ) is also a solution. Therefore the solution z 1 of (27) is an odd function. By using this fact, one can recognize that z * is larger than z 0 (0) because z 0 (0) vanishes but there is a non-vanishing z 1 (0). We found a solution, z 0 (σ ) + λz 1 (σ ), numerically and plot the solution in Fig. 1.
To find the effect on entanglement entropy, let us come back to (22). In order to get the z * , one can consider a shift of the tip along the σ direction. It is denoted by λ , then the regularity condition of the tip is given by This condition gives us = −z 1 (0)/z 0 (0). Therefore, it turns out that z * is From this, we found that the tip change is the order of λ 2 . However, our background (14), is only valid in the linear order. Therefore, one needs a quadratic background metric to see the effect on entanglement entropy for the strip case. We leave this as a future work. Although this result doesn't give us the linear variation of entanglement entropy, it is still worth considering the zeroth order calculation as a leading part of another specific kind of entangling region. We will continue our study to a case with the wedge type entangling region in the Sect. 4. Such a region is regarded as a tail part of the strip case. See the cartoon in Fig. 5. Also, once the very sharp wedge limit is taken Fig. 1 The intersection of a deformed minimal surface by the electric field (Solid curve) withβ = 1,M = 1 andq = 1: The dotted line denotes the minimal surface before applying the electric field. λ = 0.05 is taken for visualization and the right figure is the enlarged picture for the blue square near the tip ( → 0), this limit should give rise to the strip result. Due to these reasons, it is desirable to scrutinize the zeroth order calculation for the strip case numerically and analytically in the following subsections.

Numerical result for holographic entanglement entropy
Now we study the holographic entanglement entropy (22) numerically. As examined in the previous subsection, there is no linear variation of the entanglement entropy to the electric field for the strip entangling region. Therefore, in this subsection, we will focus only on numerical results of the entanglement entropy without the electric field .
First, we define the following function to facilitate numerical computation: where z * 0 ≡ z 0 (0) and˜ ≡ /l. For convenience, we replacẽ M with z h ,β andq using the relation f (z h ) = 0. Then, the useful expression for f (z) is given by By adopting a scaled coordinate u ≡ z/z * 0 , we reparameterize 4 Then the area function (30) becomes 4 Thus the parameters in terms of the original physical quantities are given by ξ = r h Also, the temperature of the system can be rewritten in terms of the dimensionless parameters as In addition, (21) and the regularity determines z * 0 as follows: Therefore, for given ξ,q andβ, the regularized area of minimal surface is determined in terms of˜ . In general, this area has a leading term that is proportional to the inverse of˜ . Subtracting this term and taking the suitable rescaling by the definition in (24), we can get the finite part of the minimal surface area that does not depend on˜ . We show our numerical results in various parameter regions in Fig. 2 and Fig. 3. In [53], the authors studied the system with momentum relaxation and showed their results in some parameter regions. Our results not only reproduce theirs well but also can be expanded to other parameter regions. Our result for the extended parameter space is telling us the role of impurity density, charge density and temperature in the entanglement entropy. As one can see, these parameters make the finite part of the entropy decrease. This finite part actually follows not an area law but a volume law [53].
It is noticeable that the above zeroth order calculation is invariant under the electromagnetic duality (F ↔F) and also the rotation of F andF in the bulk. This property leads to the symmetry exchanging and rotating the charge density and the external magnetic field in the dual field theory [41,59]. Thus the form of entanglement entropy can be extended to the entanglement entropy in the presence of the external magnetic field. The extended entanglement entropy is simply obtained by a transform :q 2 →q 2 +B 2 , whereB = Bl 2 and B is nothing but the magnetic field F xy of the black brane. Therefore our result can cover the dyonic black brane case by a simple parameter change.

Analytic result for holographic entanglement entropy
Now we try to obtain analytic expression for the holographic entanglement entropy in the small l limit. Our formulas (32) and (34) are also useful to study the analytic form. In order to study the analytic expression, one may take a small ξ limit on both equations. This limit then allows us to integrate the above equations. We obtained the entanglement entropy in this limit as the following form : where r h is given by 1/(lz h ). In addition, the γ 1 and γ 2 are defined by The physical meaning of this limit(ξ 1) is that the tip distance(z * ) from the boundary is much smaller than that of the horizon of the black hole(z h ). Roughly speaking, this limit is similar to l z h 1, but it is not the equivalent limit in a strict sense. The analytic result (35) is compared with the previous numerical calculations in Fig. 4.
In the present work, we considered the black hole involvling various parameters, temperature, charge density, and momentum relaxation. In general, the black hole horizon is given by a nontrivial function of these parameters. Our result indicates that the finite temperature effect reduces the entanglement entropy by washing out quantum entanglement entropy. In order to see the effects of the charge density and momentum relaxation, we think of their effect with a fixed black hole horizon. In this case, decreasing the charge density and momentum relaxation causes effectively the increase in temperature. This fact implies that similar to the increase in temperature, the decrease of the charge density and momentum relaxation gives the reduction of the entanglement entropy, see Fig. 3. Intuitively, the large momentum relaxation can easily occur in a high-density medium. Thus, we expect that the momentum relaxation gives rise to a similar effect to the density. In addition, as the density decreases, the reduction of the entanglement entropy can be understood as the following way. Noting that the entanglement entropy is linked to the degrees of freedom, we expect that the decrease of the density leads to the reduction of the degrees of freedom which may make quantum entanglement. This is another intuitive understanding for our results.
The above result can also have an extremal limit. At zero temperature, the charge density is given byq = 12 − 2z 2 hβ /z 2 h . So the holographic entanglement entropy at zero temperature is determined to be: This shows the effect of momentum relaxation on the holographic entanglement entropy when the temperature is zero. The above two results quantitatively indicate contribution of impurity to the entanglement entropy. It is very interesting to compare these results with the effects of impurities in weakly coupled field theories. On the other hand, the leading contribution in (37) is proportional to β 2 . It is remarkable that impurities reduce entanglement entropy in a very narrow strip and this result can be compared to known results in condensed matter physics, our system is different from systems studied in condensed matter physics though. See [60] and its references for recent articles. One may read [61] for a recent review.
As a final remark, we may think the strip case with a rotated electric field. This is achieved by replacing the subindex x in (4) by the index i denoting (x, y) with nonvanishing y- WhenT has a small value, the analytic result gives a good approximation components. However, the resulting action is still the same as the previous one up to the redefinition of the coordinate in 17. This fact implies that the rotation of the electric field doesn't give any nontrivial contribution to the action of the minimal surface. In Sect. 4, we will consider the rotated electric field with a wedge entangling region which, even at the linear order, leads to a nontrivial effect on the entanglement entropy.

A comparable deformation of ABJM theory
Our results (35) and (37) tell us how the entanglement entropy changes under deforming the CFT vacuum with operators dual to the axion field and the bulk U(1) gauge field in (1). The background solution (2) describes the deformation more explicitly. Since we do not turn on any magnetic field, the operator dual to the gauge field is nothing but the charge density operator Q and the source is given by the chemical potential. In addition, the axion field which consists of two massless real scalar fields are dual to dimension 3 operators O =3 I .
One of well-known examples of gauge/gravity correspondence is the duality between ABJM theory and Ad S 4 ×S 7 /Z k [52]. If we take the planar limit of ABJM theory, then the dual geometry becomes Ad S 4 × CP 3 . Therefore the ABJM theory becomes dual to a 4-dimensional supergravity model that is obtained from the type IIA supergravity compactified on CP 3 . Thus one can find candidates for the fields in (1). The bulk local gauge symmetry appears as a global symmetry in the boundary. The ABJM theory has SU (4) R-symmetry and the U (1) symmetry can be embedded in the R-symmetry. So it is possible to identify Q with the charge operator of the diagonal U (1) of SU (4) R-symmetry, Q R . In addition there are many chiral primary scalar operators with dimension 3, which are candidates for O =3 I . See [66] and table 1 in [67]. Therefore the following deformation of the ABJM theory is dual to (1): where S AB J M is the action of the ABJM theory and S I =3 are source for the scalar operators in the ABJM theory. And it is given by S I =3 = (β x, β y). There could be some difference coming from nonlinear effects of compactification on CP 3 . If, however, we consider small deformations and small size of the entangling region, the correction can be negligible. Then, it is valuable to compare the entanglement entropy with the above deformed action to (35) and (37). We hope to see a comparison of our result and the entanglement entropy obtained by a developing technique of supersymmetric field theory.

Holographic entanglement entropy in electric field: wedge entangling region
In this section, we consider the holographic entanglement entropy with wedge entangling regions. This type of entangling surfaces has been studied in [16,17,68,69] to consider the corner contribution of entanglement entropy. It turns out that this corner contribution contains universal information of conformal field theory, which will be discussed later. See Fig. 5 for a cartoon of a minimal surface anchored to a wedge region. It is useful to take polar coordinates (ρ, θ ) in the x − y plane as the coordinates of the minimal surface. Our main concern is how much the holographic entanglement entropy changes at the linear level under the weak electric field. To see this linear response, the direction of the electric field must be rotated by an angle δ, and one of the Cartesian coordinates x is given by x = −ρ sin(θ − δ). The need for the rotation was inferred by lessons of the previous section.
Considering the above configuration, the background geometry can be written as follows: And we may identify the coordinates of the minimal surface as the following way: Then the action for the minimal surface is given by We are interested only in a variation of the action under the electric field and regard the wedge as a part of a larger minimal surface like the cartoon (b) in Fig. 5. So we consider a part of the minimal surface with the region, ρ < R and we will see how the surface changes for this finite region. By this reason, the integration for the radial coordinate ρ is replaced by the integration between 0 to R. Furthermore, we will take a convenient choice of parameters obtained by a scaling with respect to : Then, the action becomes where f (w) is given by The above action determines the minimal surface by solving the equation of motion for w(ρ, σ ). In general, the equation is not easy to solve. Now let us take a particular limit considering a sharp wedge. Then one can use as another small parameter. The action up to the quadratic order in is given by where Since is considered small, the solution w is closed to the form of ρ times a function of σ . Thus we assume that w can be represented by a polynomial of ρ. Under this assumption, we found a solution of the equation of motion as the following form: +λ H r x 2 g 2,3 (σ )ρ 3 cos δ+ 3 g 3,3 (σ )ρ 3 sin δ . (48) The differential equations for the functions, h p,q and g p,q are given in Appendix A. In addition, the regularized area is given by where α(σ ) + λ β(σ ) is the cut-off line obtained fromz = and the explicit form is given in Appendix A. In this regard, σ max and σ min denote the maximum and minimum values of σ satisfying α(σ ) + λ β(σ ) = R. And they can be expanded as σ max,min ∼ σ (0) max,min + λ σ (1) max,min , respectively. The regularised action to the linear order in λ is Plugging (48) into the above action, we can find the onshell action as follows: where we used σ (0) max . And the first order correction for σ min and σ max is given by This can be easily derived from R = α(σ min,max ) + λ β(σ min,max ) . The zeroth order part in λ is given by the first line in (51). In fact this part shows universal features of underlying CFTs [16,17,70,71], which are summarized as follows. This term has two kinds of divergent terms. One is the usual UV divergent term which is proportional to 1/ that appears commonly in 2+1 dimensional field theory systems. One can see the term also in the strip case, (35). In addition to this, another log-divergent term shows up. Since this term is coming from the singular corner of the wedge, it is called the corner contribution. The general structure of the holographic entanglement entropy is as follows: where L is a typical length scale characterized by the size of the entangling region. Even though the constant B governs the leading term, it crucially depends on UV regulators. C (0) is the finite part of entanglement entropy. The most interesting quantity is the coefficient of log-term, A( ). The characteristic features of A are determined by two limits of : where κ 1 is conjectured to be related to, so called, entropic c-function and σ 1 is given by σ 1 = π 2 24 C T . C T is the central charge appearing in the vacuum two point function as follows: where I μν,αβ is a dimesionless tensor which is completely fixed by conformal symmetry [72]. It is worth noting that the coefficient of the log-divergent term A is UV regulator independent, as opposed to B. In addition these results were generalized to the Reyni entropy. See [17] and references therein.
Even though the above zeroth order part is very interesting, our main interest here is the linear order variation. So we focus on the first order part. If we write down only the first order part using some algebra in Appendix B, the linear action in λ is : where Here the contribution of α(σ ) 2 term to the integration is much smaller than the contribution of R 2 term. See Appendix B for the explicit form of α(σ ) that is order of E. Thus we drop the α(σ ) 2 term. If we drop out the discussed higher orders in and the integrations of some odd functions, 5 the leading change of the regularized minimal surface to the external electric field is given by where h(σ ) = h 0,1 (σ ) and g(σ ) = −g 3,3 (σ ). In addition g and h satisfy the following equations : One can solve the above equations numerically and perform the numerical integration in (59). Finally, we can get the linear response of the holographic entanglement entropy to the electric field as follows.
where N is the integration in (59). The integrand of N is a positive function over − 1 2 ≤ σ ≤ 1 2 , thus N does not vanish and it is about 0.05. Therefore, the linear response of the holographic entanglement entropy to the electric field is proportional to the thermoelectric conductivityᾱ which is a measurable quantity.
Intriguingly, the above result shows a novel relationship between the entanglement entropy and linear response. On the contrary to the present work, the authors in [49,62] investigated the evolution of a quenched system with U(1) current by using the method developed in [63,64]. They showed that in the quenched system, the first law of entanglement entropy satisfies the relation, ∂ t δS E E ∝ σ E 2 . In [49], the authors focus on the fluctuation of a small mass parameter, m ∝ E 2 , which is determined by the regularity condition at the horizon with an electric field [65] However, we showed in the present work that the entanglement entropy of a rotated wedge entangling region can lead to δS E E ∝ᾱ E sin δ. The reason we obtained the entanglement entropy correction proportional to E is that we took into account the rotation (denoted by δ) of the wedge entangling region. For δ = 0, this linear term disappears and the next term proportional to E 2 occurs as a leading correction, which is consistent with the result of [49,62]. In condensed matter community, related physics is a quite active research topic and has produced many interesting results. See [61] for a review. The origin of the linearity of the thermal conductivityᾱ comes from turning on the metric component g r x by the external electric field. Together with a reasonable gauge choice, this component describes the thermoelectric current. This hallmark is universal in various holographic models with momentum relaxation. Therefore, the above linear relation between the entanglement entropy variation and the thermoelectric coefficient is quite noticeable as a universal feature of entanglement in strongly coupled systems. Naturally, one can ask the reason why the leading effect is from the heat current mediated by the thermoelectric coefficient not from the electric current by the electric conductivity. This question is very difficult to answer because such an effect is closely related to the Fermi surface. We leave this investigation for our future work.
The measurement of entanglement is a very intriguing topic in various areas of condensed matter physics. The first measurement of the second Rényi entropy was accomplished in [73]. This work was inspired by [74,75]. Except for this measurement, there are various trials to detect entanglement measure. Although the accomplished detections were seen through few-body correlations, there are many discussions and proposals on extension of the experiments. Among them, a suitable proposal to our result is based on quantum antiferromagnets. See [77]. In the proposal, the entanglement measure can be encoded in magnetization of a system. A circular superconducting plate covers on both sides of an antiferromagnet sheet. By the Meissner effect, the plate plays a role of dividing entangling regions. Assuming that such an experiment is realized, we propose an experiment with a wedge type superconducting plate and applying an rotated electric field along a tangential direction. Then, one may check whether resulting entanglement entropy data is related to the thermoelectric coefficientᾱ or not. Our idea is a bit rough, but worth thinking about more.

Discussion
In this work, we have studied the entanglement entropy affected by the external electric field. Though the entanglement entropy has usually been considered as an important and fundamental quantity representing a variety of quantum natures, it still remains a big issue to measure the quantum entanglement entropy in the laboratory. One of the main goals of this work is how we can relate the quantum entanglement entropy to the measurable quantities like the transport coefficients. This would be useful to figure out the relation between the quantum phenomena and various macroscopic quantities. Moreover, the present work may provide new intuitions about how to measure the entanglement entropy in the laboratory.
More specifically, we have taken into account the entanglement entropy of the strip and wedge entangling region with turning on the electric field and momentum relaxation. If there is no momentum relaxation which breaks the translation symmetry, the DC conductivity usually has an infinite value. In order to obtain a finite conductivity in the zero frequency limit, we first considered the geometry with the non-vanishing momentum relaxation which resembles introducing impurities to the dual field theory. This momentum relaxation usually modifies the electric and thermoelectric conductivities. Theoretically, those changes of the transport coefficients can be governed by the response theory and we can easily measure such changes experimentally in the laboratory. In this situation, we can ask how the entanglement entropy is affected by the external electric field and what is the relation between the modified entanglement entropy and transport coefficients. In this paper, we showed how the entanglement entropy modified by the external electric field using the holographic method. It turns out that the entanglement entropy change can be connected to the transport coefficient, especially the thermoelectric conductivity.
For the strip case, we found that turning on the external electric field tilted the minimal surface as expected. In the holographic model, the area of the minimal surface is directly associated with the entanglement entropy. Thus, we may expect the change of the entanglement entropy due to the external electric field. We found that such a change of the entanglement entropy does not occur at least in the linear response theory. This is because the holographic formula of the entanglement entropy in (18) has an invariant form under the parity transformation like x → −x. Despite the tilted minimal surface, the invariance under the parity transformation does not allow the change of the entanglement entropy at least at the linear order. However, the higher order corrections can affect the resulting entanglement entropy. In fact, we explicitly showed that the second order correction caused by the tilted minimal surface can change the entanglement entropy, although we did not regard an additional contribution caused by the background metric deformation. However, we did not present the result in this note because we leave it as part of future works. The additional contribution is related to the response theory at the second order and we hope to report more results on this issue in future works. Anyway, the results in the strip entangling region showed the fact that the change of the entanglement entropy can be represented in terms of transport coefficients.
In order to get more clear and explicit relation at the linear order, we further considered the entanglement entropy in the wedge entangling region which is utilized to extract universal information about the corner contribution. As shown in the strip case, the entanglement entropy invariant under the parity transformation does not give any nontrivial contribution at the linear order, so that we took into account the external electric field which is rotated with an arbitrary angle. Since the external electric field with an arbitrary rotation usually breaks the parity invariance, one can expect the nontrivial contribution to the entanglement entropy even at the linear order. We showed with explicit calculations that the linear order correction to entanglement entropy really occurs in this setup. Intriguingly, we furthermore found that the change of the entanglement entropy is directly related to the thermoelectric conductivity. This is the first example showing how the entanglement entropy can be connected to the transport coefficients. In addition we provided an rough experimental idea for our result of the wedge case and we also described a comparable deformation in the ABJM theory to the strip case. In future works, we hope to report more evidences and understanding of the underlying structure of the connection between the entanglement entropy and transport coefficients.  cos δ, (71) where E = is a scaled cut-off.