Lifshitz entanglement entropy from holographic cMERA

We study entanglement entropy in free Lifshitz scalar field theories holographically by employing the metrics proposed by Nozaki, Ryu and Takayanagi in \cite{Nozaki:2012zj} obtained from a continuous multi-scale entanglement renormalisation ansatz (cMERA). In these geometries we compute the minimal surface areas governing the entanglement entropy as functions of the dynamical exponent $z$ and we exhibit a transition from an area law to a volume law analytically in the limit of large $z$. We move on to explore the effects of a massive deformation, obtaining results for any $z$ in arbitrary dimension. We then trigger a renormalisation group flow between a Lifshitz theory and a conformal theory and observe a monotonic decrease in entanglement entropy along this flow. We focus on strip regions but also consider a disc in the undeformed theory.


Introduction
Lifshitz theories are characterised by a scaling symmetry under which space and time scale differently: x → λ x, t → λ z t , (1.1) where z is referred to as the dynamical exponent. Such theories govern quantum critical points in many condensed matter systems. As a simple example, consider the following theory of a free massless scalar in d + 1 dimensions: Little is known about entanglement entropy in such theories. For the special case of d = z = 2 one can map the ground state to a Euclidean conformal field theory in 1 + 1 dimensions and compute some subleading universal terms in the entanglement entropy [2,3,4,5,6,7], see also [8,9]. Other work on z = 2 in d = 1 can be found in [10,11]. This will be useful for us to fix the overall normalization factor in the entanglement entropy obtained from cMERA for the case z = 2, as we discuss in Section 3.1.
The discretised theory with d = 1 and arbitrary z, including a mass term, was studied recently in [12,13] and some partial results for d = 2 were obtained in [12]. In [13] an analytical approach was also put forward based on the holographic cMERA technique [1] for Lifshitz scalar fields. In this paper, we follow up and elaborate on this discussion and extend it to higher dimensions, obtaining several new results.
Holography opens up a new way to compute entanglement entropy. The difficult direct field theory calculation is mapped to a geometric extremisation problem in the dual gravity theory via the Ryu-Takayangi (RT) prescription [14,15]. In detail, first one computes the bulk metric that describes the field theory state of interest. The entanglement entropy of a given region is then equal to the area of the minimal area extremal surface, as measured by this metric, that ends on this region at the boundary.
However, the holographic dual of a Lifshitz theory has yet to be universally agreed upon.
Whilst a dual spacetime, termed Lifshitz spacetime, was proposed in [16,17] and has been studied intensively ever since (see [18] for a review), it is unclear whether the RT prescription should be applied to this spacetime. Indeed, a study of various perspectives on the holographic reconstruction of Lifshitz spacetime can be found in [19]. Other recent work suggests that Newton-Cartan geometry may provide a more natural bulk dual for a non-relativistic theory [20]. Regardless, holography typically computes the entanglement entropy for stronglycoupled field theories with large central charges, whereas here we focus on a very different setting.
In this paper we use a method inspired by holography that is applicable to free field theories. In particular, Nozaki-Ryu-Takayanagi (NRT) proposed in [1] that a metric emerges from a continuous version of the multi-scale entanglement renormalisation ansatz (cMERA) [21]. For a given theory, expectation values of the appropriate disentangler operator determine various components of this metric. In some sense, the NRT proposal geometrises the entanglement entropy of free fields. Our goal is to compute entanglement entropies by applying the RT prescription to the cMERA metric for various Lifshitz theories. As we explain later, our results should be viewed as predictions for field theory calculations. We should stress however that, while the NRT proposal is holography inspired, it is not embedded within the AdS/CFT correspondence, since we apply it to a single free scalar field that is neither large N nor at strong coupling. It is well known that free conformal fields do not have a gravity dual, at least not with a weakly coupled gravity sector coupled to matter fields. Nevertheless, the NRT proposal is similar in flavor to the RT prescription at a technical level (extremisation of surface areas using metric geometries), so we can make concrete calculations. The justification of this comes from the MERA approach, and the expectation that a continuum version of it should exist.
We begin in the following section with a brief review of cMERA and the definition of the cMERA metric. As a crucial consistency check we first compute entanglement entropy in the ground state of a relativisitic free massive scalar theory and compare with known results. We then turn to our three main calculations of Lifshitz entanglement entropy in Section 3. We consider the original Lifshitz theory (1.2) and two relevant deformations: a mass term m 2 φ 2 and a relativistic term ( ∇φ) 2 . We conclude with a discussion in Section 4.

cMERA and holographic entanglement entropy
The multi-scale entanglement renormalisation ansatz (MERA) is a variational approach based on the renormalisation group (RG) to construct (ground) states in quantum many body systems and study their entanglement properties [22,23]. A continuum version was developed in [21] for free fields, which we now summarise. We follow the notation of [1].
Choose a quantum field theory in d+1 dimensions and introduce a length cut-off ε. Consider then a one-parameter family of states |Ψ(u) living in the Hilbert space of this theory. The dimensionless parameter u keeps track of the current length scale, with u UV = 0 and u IR = −∞ in the ultraviolet and infrared, respectively. Focus initially on a reference state |Ψ(u IR ) that has no entanglement. Run up the RG scale and generate entanglement by acting with a unitary transformation based on a local operator K(u). Next act with a scale transformation L to introduce new degrees of freedom at shorter length scales. Repeat this process until the UV cut-off is reached. The final state we are interested in can then be written as and the variational principle can then be applied to minimise the energy of this state. The variational parameters are encoded in the coefficients of the interactions in K(u). Note that the states |Ψ(u) are manifestly translationally invariant if these coefficients are independent of x.

Metrics from cMERA
The authors of [1] associate a metric in d + 2 dimensions with a type of cMERA. In particular, they argue that a cMERA yielding a translationally invariant ground state should correspond to a metric of the form The metric in the RG direction parametrised by u is given by the Hilbert-Schmidt distance between cMERA states at nearby scales u. It can be expressed in terms of the variance of the operator K(u) in the state |Ψ(u) 1 : This metric component effectively measures the density of disentanglers at the scale u. The g tt component cannot be determined from the field theory entanglement on a fixed-time slice and therefore plays no role in our discussion. The explicit form of this metric can be calculated for the ground state of a free scalar theory [1], given the choice of disentangler proposed in [21]. It depends purely on the dispersion relation ω(k), with k ≡ | k|, and will be used throughout this paper: The cMERA metric constructed in this way captures how the quantum degrees of freedom in the field theory are entangled with each other at different RG scales. It has the flavour of holography, but the precise connection with AdS/CFT is not understood since we are neither at large N nor at strong coupling and no expression for g tt is provided. For relativistic conformal field theories, the holographic cMERA approach gives the AdS metric [1]. And as we review below, applying the RT prescription to the cMERA metric for massive scalar fields yields the correct answers for the entanglement entropy when the correlation length is small. Based on this, we take a pragmatic approach and apply the holographic cMERA techniques to nonrelativistic theories with Lifshitz scaling and with mass deformations corresponding to small correlation lengths. The method yields predictions for the Lifshitz entanglement entropy for general values of the dynamical exponent z.

Consistency checks in relativistic theories
Before embarking on our main calculations, in this section we illustrate the consistency of this approach with a simple example: a free massive scalar field in 1 + 1 dimensions. The action for this theory is (we absorb a factor c 2 / in m such that m has dimension of inverse seconds) and the dispersion relation is ω(k) 2 = c 2 k 2 + m 2 . Computing g uu via (2.4), we should therefore associate the following metric with the ground state of this theory: For m = 0 this reduces to the metric on AdS 3 in Poincaré coordinates with z = εe −u and g tt = −c 2 /z 2 , together with a simple rescaling of x and t. The AdS radius in this normalisation is R AdS = 2, but one can rescale the overall metric to get any radius.
Our task now is to compute entanglement entropy from this metric using the RT prescription. We will focus on an interval in the x-direction of width at t = 0. It is useful to define the dimensionless quantities We want to compute entanglement entropy as a function of J 1 and J 2 . We require J 1 < 1 and J 2 > J 1 to ensure that the cut-off ε is the smallest length scale in the theory. Furthermore, the regime J 2 < 1 means that the correlation length ξ ≡ c/m is larger than the subsystem size , whereas J 2 > 1 means that the correlation length is smaller than the subsystem size.
The entanglement entropy is known to be different in these two regimes [24] -a result we now re-derive. First we change coordinates: The boundaries of the r coordinate are fixed by the limits of the cMERA length scale u: We seek geodesics of this metric parametrised by x(r) that end at x(r UV ) = ± /2. The length of the shortest geodesic is proportional to the entanglement entropy. The steps to find the appropriate geodesics and compute their lengths are identical to those in [14,15]; we adapt them to our setup and only present the results and a few intermediate steps.
Two types of geodesic are relevant for a non-zero mass. See Figure 1 for examples. The first type connects the endpoints and is smooth at the point of deepest extent in r: x(r ) = 0 and x (r ) → ∞. The boundary condition for this type of geodesic relates r to the interval width : At fixed J 1 , the function on the right-hand side is real and non-zero for J 1 (1 + J 2 1 ) −1/2 < a < 1, in accordance with (2.9). It is positive within this range and has a single maximum. Thus, equation (2.10) cannot be satisfied for large J 2 , but may have two solutions for small J 2 . The second type of geodesic exists for all J 2 and consists of two disconnected straight sections x(r) ≡ ± /2 that end at r = r IR .
We must be careful to identify the geodesic with the shortest length for a given J 2 . The lengths of the connected and disconnected geodesics are given respectively by We demonstrate in Figure 2 that the disconnected geodesic is shorter than any connected geodesic above a critical value of J 2 . This value is slightly below that for which a connected geodesic no longer exists.  Figure 2: Ratio of the connected geodesic length to the disconnected geodesic length for J 1 = 1/5 as a function of J 2 . The red, green and blue dots correspond to the curves plotted in Figure 1 for J 2 = 1/2 and the black dot marks the critical value of J 2 for this J 1 . The region to the left of the vertical dashed line is unphysical since J 2 < J 1 therein.
In the massless case (i.e. setting m = 0 from the beginning in (2.6)) the connected geodesic is always the shortest. We can find its length explicitly as a function of /ε: In the last equation, we made an expansion in small ε/ . The first term is the dominant term and leads to the area law, which in 1 + 1 dimensions is logarithmic in . The entanglement entropy is proportional to the length of the geodesic and the proportionality factor is known to be related to the central charge [25]: with c = 1 for a real scalar field.
Given that we have fixed the overall normalisation, we can return to the massive case and perform another consistency check for the case when the correlation length is smaller than the size of the interval. This translates into large J 2 , for which the disconnected bulk curve (2.12) is the shortest. Expanding for small J 1 and using the same normalisation factor c/3 as before, we find the following leading term: where we recall that the correlation length is defined by ξ ≡ c/m. This result for the entanglement entropy matches perfectly with the universal result of Cardy and Calabrese in [24] applied to the case of one interval with two boundary points.
In the opposite limit, i.e. for large correlation lengths compared to the length of the line interval, the cMERA approach seems to not reproduce known results. Namely, it was shown in [26,27] that the entanglement entropy for a free massive scalar field in 1+1 dimensions in the regime ξ contains a term proportional to log(− log(mε/c)) = log(− log J 1 ) in the limit of small J 1 . It is unclear to us how such a double logarithm would appear in the cMERA approach (see also Section 3.3 for the case of general z), and it would be interesting to understand better where this apparent discrepancy comes from. This concludes the consistency checks on the holographic cMERA techniques applied to free relativistic scalar field theories. The generalisation to higher dimensions can also be done, but we include it in the general analysis of Lifshitz theories with arbitrary dynamical exponent z.
When we set z = 1 we recover the relativistic results.

Lifshitz entanglement entropy
We now apply the holographic cMERA technique to compute entanglement entropy in free Lifshitz scalar field theories. This possibility was in fact already pointed out in [1] but not worked out. In fact, in that reference, the case z = 1 treated in the previous section was also not worked out in detail. The case of d = 1 was recently also considered in [13] which we review and extend below.
As suggested in [1], all we need is the dispersion relation and the rest is computational. In this section we perform all the calculations and provide results for arbitrary z and d, treat the massive case as well, and compute the RG flow of the entanglement entropy from z > 1 theories to z = 1.

Massless Lifshitz scalar
First we consider the simple theory given in (1.2) with d = 1. This has dispersion relation ω(k) = αk z . According to the prescription given in the previous section, we should associate the following metric with its ground state: We wish to compute the entanglement entropy of an interval in this state, and for this we calculate the length of the geodesics in this metric. Remarkably, the result is a simple rescaling of the massless relativistic result (2.13): This can be seen immediately at the level of the metric. Suppose we scale out the factor z 2 in front of the metric. Then up to this overall rescaling, this is the z = 1 metric provided we replace the cutoff ε → zε in (3.1). (Physically, we are not changing the cutoff ε for the Lifshitz theory, it is just a trick to get the answer. Equivalently, one can also rescale x → zx.) It is then clear that the length of geodesics on the t = 0 slice is just a rescaling by z of that obtained from (2.6) at m = 0. This observation was already made in [13], but there the result was only given for large /zε, which does not permit taking large values of z.
Here, we also obtain a sensible large z limit. In a discrete model of a continuum local theory, entanglement entropy is dominated by contributions from nearest neighbors across the boundary of the entangling region. As z is increased, more and more lattice sites are involved and we expect the full interval to contribute at infinite z. This is indeed what we find: i.e. it becomes extensive with the region size, so we recover a volume law in the large z limit. This agrees with the observations and results of [12]. It has to be noted though that this limit is only natural in the continuum limit (where also the cutoff ε is sent to zero), as long as z l/ε.
The entanglement entropy S follows from L by a multiplicative factor. In AdS/CFT, this factor comes from translating Newton's constant into the central charge of the theory, which in 1+1 dimensions gives a multiplicative factor of c/3, leading to S = c/3 log( / ) + c 0 for relativistic CFTs. The cMERA approach, however, does not fix the overall normalization. The best we can do is to replace the central charge c by an overall multiplicative constant c z , independent of but which can depend on the dynamical exponent z. For a relativistic scalar field, we have c z=1 = 1. In general, we call c z the Lifshitz central charge. Then we obtain the following entanglement entropy formula for general values of z: There might still be a non-universal additive constant independent of , just like the coefficient c 0 that is cutoff-dependent. We will leave it out of the discussion here.
The expansion around small values of /zε gives small deviations from the volume law, whereas expanding around the area law /zε 1 gives deviations from the area law, In [13], the assumption was made that c z is independent of z. Recent numerical work, however, shows that this assumption should be relaxed, as the work of [10,11,28] showed that for z = 2, one has c z=2 = 3/4, such that the entanglement entropy for z = 2 starts like S = 1 2 log( ) to leading order. This implies that c z does depend on z. But the fact that we don't know c z does not prevent us from making predictions that could be checked using numerics or other methods.
Indeed for fixed but arbirtray z, we can still check the functional dependence on /z inside the logarithm in (3.4) with numerical methods. Or stated differently, we can take the ratio of the entanglement entropies for large and small values of /z and compare this to lattice results.
This ratio is independent of c z , so is not sensitive to our ignorance of it. Admittedly, this issue needs to be better understood, either from further numerical work, or perhaps analytically, by computing c z from the replica trick or from Lifshitz scale anomaly coefficients. We leave this for further study, and leave the coefficient c z undetermined in this paper. More conveniently, we will focus on just the length or area of the minimal surfaces, which do not involve this factor c z .
We now generalise our treatment to higher dimensions. We make the replacement dx 2 → d x 2 d with x = (x 1 , ..., x d ) and also the change of coordinates e 2u = ε 2 r 2 in (3.1): Our region of interest is now the strip x − 2 ≤ x d ≤ 2 at t = 0. Translational invariance of the metric and strip in the additional spatial directions simplifies the problem dramatically: we seek extremal surfaces of this metric parametrised by x d (r) that fill the other directions and end at x d (ε) = ± /2. The area of the smallest surface is proportional to the entanglement entropy.
We find the following expressions for the width of the strip and the area of the surface extending to r = r , respectively: where we have defined b ≡ r /ε and regulated the infinite volume of the remaining directions with V d−1 . Note that these results are again related to those of the relativistic case with a rescaling by z: we scale ε → zε and multiply A → z d A .
We deduce from these results that the entanglement entropy follows an area law for finite z.
The equation (3.8) relating and r has a unique solution for given d and z. We can invert this asymptotically at large /ε and find the following expansion for the area: (3.10) The first term represents the area term, and the second term is finite and independent of the cutoff. (See also, for example, [15] for the case z = 1 using AdS/CFT and the RT formula.) The entanglement entropy is proportional to A, with a proportionality factor that is not universal for d > 1. In the cMERA approach, this normalisation inherits from that of the disentangler operator K(u) and is usually fixed by comparing with known field theory results for the entanglement entropy. Because of the non-universality, we leave this overall factor undetermined.
In the large z limit, we again obtain a volume law instead of an area law. We can invert (3.8) asymptotically at large z, finding We substitute this into the area expression (3.9) and find which is indeed proportional to the regulated volume V d−1 of the strip.

Disc geometry
In this section we depart briefly from the strip to consider a region of finite size: the disc. We study the same state of theory (1.2) but simply write the appropriate cMERA metric in different coordinates: The following function describes a surface lying on the constant time slice t = 0 that ends at r = ε on a disc of radius p = R: (3.14) This surface is smooth at p = 0, independent of d, minimises the area functional and has area The area diverges as the cut-off is taken much smaller than the disc radius. For finite z, this divergence is proportional to the area of the disc: The lowest subleading divergence is logarithmic when d is odd; in particular, we recover (3.2) for d = 1 since the interval width satisfies = 2R. For even d, the small x expansion contains a constant term which is universal. All these results are straightforward generalisations and rescalings of the z = 1 case described in, for example, [15].
It is interesting to compare our results with [2], where the entanglement entropy was studied for d = 2 with z = 2. In particular it was found that for the disc subspace, no logarithmic terms were found, see case (a) in Figure 1 in [2]. This is consistent with our general result that for even d, no logarithms appear in the expansion of (3.15). More interesting would be to study the geometries where logarithmic terms do arise, such as the rectangular or half disc subspaces studied in [2] as well. To reproduce the coefficients in front of the log-terms using the cMERA approach, we would need to perform extremisation of surfaces in geometries that end on these rectangles or half discs that have less symmetry. This is much harder in d ≥ 2 but is an interesting problem for separate study. Given that we don't know proportionality coefficient between the area and the entanglement entropy, it is cumbersome to perform such a test at present. In the large z limit we again obtain a volume law: Here B d R is the d-dimensional ball whose boundary is the (d−1)-dimensional sphere of radius R. From now on we focus exclusively on strip regions.

Massive deformation
We now deform the theory (1.2) by adding the mass term m 2 φ 2 . The dispersion relation is modified to ω(k) 2 = α 2 k 2z + m 2 , leading to a cMERA metric The dimensionless parameters that characterise the mass deformation and the region size can be chosen respectively. They reduce to (2.7) in the relativistic limit. We have written them in terms of the correlation length ξ, generalising the one from the relativistic case. We want to compute entanglement entropy as a function of d, z, J 1 and J 2 . Again we require J 1 < 1 and J 2 > J 1 to ensure that the cut-off ε is the smallest length scale in the theory. As in the relativistic case, the correlation length maybe larger or smaller than the interval length or strip width , corresponding to the regimes J 2 < 1 or J 2 > 1 respectively. Just like the relativistic case covered in Section 2.2, there is a competition between two types of extremal surface: connected and disconnected. We analyse this case in a similar fashion. We begin with a change of coordinates: For a connected surface extending to r = r , the width of the strip and the area of the surface can be written

21)
respectively, where it is convenient to parametrise the bulk depth via . (3.23) We can evaluate these integrals exactly for various values of d and z; more generally, they are straightforward to evaluate numerically. Equation (3.21) relates J 2 and r . We find that it may have zero, one or two solutions depending on the value of J 2 , just like the relativistic case. We observe an area law divergence in A C that degenerates to a logarithmic divergence for d = 1.
For example, we have obtained analytical expressions for d = 1, z ∈ N: (3.24) where the γ n satisfy 1 + (−1) z γ z n = 0 =⇒ γ n = e iπ(2n−1+z)/z , n = 1, 2, . . . , z . The disconnected surface exists for all J 2 . Its area is in fact independent of J 2 and can be evaluated in closed form: This result is again proportional to the area of the region. For d = 1 we find a simple expression that diverges logarithmically as J 1 → 0: (3.28) We must now identify which surface has the least area for a given region size J 2 . For fixed d, z and J 1 we find that the disconnected surface has minimal area above a critical value of J 2 .
This value is slightly below that for which a connected surface no longer exists. We observe the same qualitative behaviour as presented in Figure 2 for the relativistic case regardless of d and z. Besides yielding analytical solutions, the point z = d does not appear to be special in this calculation.
For any given J 1 it is straightforward to numerically find the critical value of J 2 at which the two types of surface have equal area. In this way we can construct the phase diagram of a given theory. Our results in d = 1 for various values of z are plotted in Figure 3. Note that the J 2 -intercepts of the phase boundaries follow a roughly linear relationship: J 2 (J 1 = 0) ∼ z/2.
For large values of J 2 the disconnected curve is always the shortest. This is the limit in which the correlation length is smaller than the length of the interval: where we have used the same normalisation factor c z /3 as before that relates the length of the geodesic to the entanglement entropy for a real scalar field. For large correlation lengths ξ > , so small J 2 (and therefore small J 1 ), we get again logarithmically diverging terms − log J 1 . There are no known analytical or numerical results in this case, but we expect similar discrepancies as for z = 1, where the leading diverging term involves a double logarithm, log(− log J 1 ) in the limit of small J 1 .

RG flow to an IR CFT
In this section we begin with a massless free scalar field in 1 + 1 dimensions and turn on an irrelevant Lifshitz coupling: Our goal is to compute the entanglement entropy along the entire renormalisation group flow from the UV Lifshitz theory to the IR CFT. For sublattice entanglement entropy, such a study was carried out in [13], where it was found that the entanglement entropy decreases along the renormalisation group flow, for any starting value of z > 1. This provides evidence for a generalisation of the c-theorem for entanglement entropy in the relativistic case [29], applicable to flows between two Lifshitz fixed points. Our analysis below will give further support for this. The dispersion relation is given by ω(k) 2 = c 2 k 2 + α 2 k 2z . After a change of coordinates, this results in a cMERA metric of the following form: The function f (r) interpolates between the two limits of the flow as the dimensionless parameter K ≡ α/(c ε z−1 ) is varied: We seek geodesics of the form x(r) extending to r = r that end at x(ε) = ± /2. For example, for z = 2 we find that the interval width and the geodesic length can be written respectively as where we have again defined b ≡ r /ε. Note that the expression for /ε increases monotonically with b (at fixed K) so can be inverted uniquely. It is straightforward to check that these results have the correct limits: (3.2) for large K and (2.13) for small K. These are the UV and IR limits, respectively. We plot the geodesic length as a function of K or /ε in Figure 4. We see that the length decreases monotonically along this RG flow from UV to IR, as expected.
(We have checked that qualitatively similar results can be obtained numerically for other values of z.) L L /ε log 10 K Figure 4: Entanglement entropy for a flow between a CFT in the IR and z = 2 Lifshitz theory in the UV. Left: Length as a function of log 10 K for /ε = 2, 5, 10, 20 (bottom to top). These curves interpolate monotonically between the correct IR (left) and UV (right) limits. Right: Length as a function of /ε for K = 0.5, 5, 50 (bottom to top, solid). The dashed lines correspond to the IR (black) and UV (red) limits. and small N . We have illustrated in this paper that this approach leads to concrete predictions for the entanglement entropy of free scalar fields. In the relativistic case, it reproduces well-known analytical field theory results, as we showed in our analysis. 2 In the Lifshitz case, it leads to new predictions, generalising and extending the observations made in [13]. The overall normalisation for the entanglement entropy is however not fixed from the metrics introduced by [1]; rather, it depends on the choice and normalisation of the disentangler operator. In AdS/CFT this normalisation is fixed by the dictionary that relates Newton's coupling constant to the central charge or number of colours. We have not been able to determine this overall normalization, which we denoted by c z , and which can depend on the value of the dynamical exponent z. It would be interesting to understand if c z is related to Lifshitz scale anomaly coefficients. These anomalies have been studied in [30,31,32], but also in earlier work for z = 2 [33,34] in 2+1 and 3+1 dimensions, and z = 3 in 3+1 dimensions [35]. The results obtained in e.g. [2] show that central charges do appear in certain logarithmic terms entanglement entropy for z = d = 2 (so 2+1 dimensions), but it is not straightforward to see the relation to Lifshitz scale anomaly coefficients. Moreover, for the strip and disc geometries we consider here, these logarithmic terms are absent. It would be important progress to find such a relation if it exists.
In the massless case, our results for free Lifshitz scalars are obtained from a simple rescaling of the relativistic case, but this is not true in the massive case. In the latter case, we generalised the Calabrese-Cardy result to Lifshitz scalar fields with values of the dynamical exponent z > 1.
It would be interesting to repeat our analysis for fermions.
Clearly, it would be important to better understand the validity of the holographic cMERA approach. From this point of view, our results are merely predictions rather than solid results.
It would be nice to reproduce some of our predictions for Lifshitz entanglement based solely on field theoretic techniques. For free Lifshitz scalar fields, one imagines this should not be so difficult by attempting the replica method, but the presence of longer range interactions for large z complicates this. In low dimensions, numerical and lattice methods can also be used, as was illustrated in [12,13]. We leave this topic for future study.