Research ArticlePHYSICS

Holography and criticality in matchgate tensor networks

See allHide authors and affiliations

Science Advances  09 Aug 2019:
Vol. 5, no. 8, eaaw0092
DOI: 10.1126/sciadv.aaw0092


The AdS/CFT correspondence conjectures a holographic duality between gravity in a bulk space and a critical quantum field theory on its boundary. Tensor networks have come to provide toy models to understand these bulk-boundary correspondences, shedding light on connections between geometry and entanglement. We introduce a versatile and efficient framework for studying tensor networks, extending previous tools for Gaussian matchgate tensors in 1 + 1 dimensions. Using regular bulk tilings, we show that the critical Ising theory can be realized on the boundary of both flat and hyperbolic bulk lattices, obtaining highly accurate critical data. Within our framework, we also produce translation-invariant critical states by an efficiently contractible tensor network with the geometry of the multiscale entanglement renormalization ansatz. Furthermore, we establish a link between holographic quantum error–correcting codes and tensor networks. This work is expected to stimulate a more comprehensive study of tensor network models capturing bulk-boundary correspondences.


The notion of holography in the context of bulk-boundary dualities, most famously expressed through the anti–de Sitter space/conformal field theory (AdS/CFT) correspondence (1), has had an enormously stimulating effect on recent developments in theoretical physics. A key feature of these dualities is the relationship between bulk geometry and boundary entanglement entropies (24), prominently elucidated by the Ryu-Takayanagi formula (5). Because of the importance of entanglement in the context of AdS/CFT (6), it was quickly realized that tensor networks are ideally suited for constructing holographic toy models, most notably the multiscale entanglement renormalization ansatz (MERA) (79). The realization that quantum error correction could be realized by a holographic duality (10) further connected to ideas from quantum information theory. Despite the successful construction of several tensor network models that reproduce various aspects of AdS/CFT [see, e.g., (1113)], a general understanding of the features and limits of tensor network holography is still lacking. Particular obstacles are the potentially large parameter spaces of tensor networks and the considerable computational cost of contraction.

In this work, we overcome some of these challenges by applying highly efficient contraction techniques developed for matchgate tensors (14, 15), which replace tensor contraction by a Grassmann-variate integration scheme. These techniques allow us to comprehensively study the interplay of geometry and correlations in Gaussian fermionic tensor networks in a versatile fashion, incorporating toy models for quantum error correction and tensor network approaches for CFT, such as the MERA, into a single framework, highlighting the connections between them. Furthermore, this framework includes highly symmetrical tensor networks based on regular tilings (see Fig. 1, A and B). We are thus in a position to efficiently probe the full space of Gaussian bulk-boundary correspondences from a small set of parameters, including the bulk curvature. We show that matchgate tensor networks with a variety of bulk geometries contain the Ising CFT in their parameter space to remarkably good approximation as a special case, with properties similar to the wavelet MERA model (16, 17). While regular hyperbolic tilings have recently been considered as a MERA alternative (18), we show that flat tilings can lead to very similar boundary states. In our studies, we restrict ourselves to tensor networks that are nonunitary and real, resembling a Euclidean evolution from bulk to boundary. In particular, we do not require the causal constraints of the MERA for efficient contraction, thus providing new approaches in the context of tensor network renormalization (19, 20). While we provide substantial evidence that tensor networks are capable of describing bulk-boundary correspondences beyond known models and introduce a framework for their study, our work is by no means exhaustive. We do hope to provide a starting point for more systematic studies of holography in tensor networks.

Fig. 1 Geometries of tensor networks.

Discretizations of flat (A) and hyperbolic space (B and C) with a triangular tiling (blue edges), into which a tensor network is embedded (black lattice). In the matchgate formalism, joint edges between triangles correspond to an integration over a pair of Grassmann numbers, analogous to tensor network contraction over indices. While (A) and (B) show regular tilings, (C) presents a nonregular MERA-like tiling we call the matchgate MERA (mMERA).


We constructed two-dimensional planar tensor networks with fermionic bulk and boundary degrees of freedom. The bulk degrees of freedom are associated with a set V of vertices of a tensor network. At each vertex, vV, a local tensor Tv with kv indices is placed, which can be interpreted as a local fermionic state on kv sites. After contraction over all connected bulk indices, the L remaining open indices are interpreted as boundary sites with the boundary state specified by the full contracted tensor. Because of the planarity of the network, the boundary sites form a loop. The bulk geometry can be flat or negatively curved (a positively curved network closes in on itself after finite distance). We visualized our tensor networks by representing each tensor Tv as a kv-gon whose edges correspond to indices. Thus, the tensor network is represented by a polygon tiling, which determines the bulk geometry. Adjacent edges between two polygons correspond to contracted indices and boundary edges to open ones. See Fig. 1 for examples.

Concretely, each bulk degree of freedom vV is associated with a local tensor Tv : {0,1}×r → ℂ of tensor rank r (equal to the number of edges of the corresponding tile), all of which are contracted to form tensors of higher rank. We denote the tensor component at indices j ∈ {0,1}×r as Tv(j) and the standard computational basis for r boundary spins as jk=1rjk. Each tensor is then equivalent to a state vectorψv=j{0,1}×rTv(j)j(1)

For a broader introduction to tensor networks and their contractions, see (2124).

Instead of explicit tensor contraction along pairs of indices, we used the formalism from (15) using Grassmann integration. Any tensor T can be represented by a Grassmann-variate characteristic functionΦT(θ)=j{0,1}×rT(j)θ1j1θ2j2θrjr(2)where the θk are Grassmann numbers defined by the anticommutation relation θkθk + θkθk = 0. The contraction T1*2 of two tensors T1 and T2 (of rank r1 and r2, respectively) over the last index of T1 and the first index of T2 is given byT1*2(x,y)=z{0,1}T1(x,z)T2(z,y)(3)where x ∈ {0,1}×(r1 − 1), y ∈ {0,1}×(r2 − 1). T1 * 2 has rank r1 + r2 − 2. The characteristic function of the contraction is obtained asΦT1*2(ξ)=dη1dθr1ΦT1(θ)ΦT2(η)exp(θr1η1)(4)where we used ξ = (θ1, …, θr1 − 1, η2, …, ηr2) and ∫ ‍ dη1 ∫ ‍ dθr1 denotes Grassmann integrals, anticommuting multilinear functionals obeying dξjξjzj=δzj,1 [see (15, 2527) for more details]. A self-contained derivation of the equivalence of (4) with tensor contraction, as well as a note on iterated integrals, is given in the Supplementary Materials. Anticommutativity requires an appropriate labeling of all Grassmann variables, but such a labeling can always be found for contractions of planar networks (15). These Grassmann integrations are particularly efficient to compute for the case of matchgate tensors, where their computation scales polynomially in the number of tensor indices.

Consider a rank r tensor T(x) with inputs x ∈ {0,1}×r. One calls T(x) a matchgate if there exists an antisymmetric matrix A ∈ ℂr × r and a z ∈ {0,1}×r so that one can writeT(x)=Pf(AxXORz)T(z)(5)where Pf(A) is the Pfaffian of A and Ax is the principal submatrix of A acting on the subspace supported by x. Furthermore, one calls T(x) an even tensor if T(x) = 0 for any x with odd ∑jxj.

A generic even matchgate has a simple Gaussian characteristic function of the formΦT (θ)=T(0)exp(12j,k=1rAj,kθjθk)(6)where T(0), the tensor component for all-zero input, acts as a normalization factor. Apart from normalization, the full tensor is completely determined by A, which we therefore call the generating matrix. Thus, the rules for contracting matchgate tensors can be written as rules for combining generating matrices. Full derivations of these, including the calculation of physical covariance matrices from the generating matrices, are provided in the Supplementary Materials. With our contraction rules, the computational cost of contracting two tensors is quadratic in the number of indices of the final tensor. Thus, we could bound the total computational cost for contracting an entire network of the type considered here by O(L2N), where L is the number of boundary sites and N is the number of contracted tensors [for similar bounds on matchgate contraction, see (15)].

Using Pauli matrices σα with α ∈ {x, y, z}, one can define Majorana operators γi via the Jordan-Wigner transformationγ2k1=(σz)(k1)σx(12)(rk)(7)γ2k=(σz)(k1)σy(12)(rk)(8)

The computational basis is then equivalent to an occupational basis. In this context, we proved that any fermionic Gaussian state vector in the form of Eq. 1 has coefficients T(j) constituting a matchgate tensor. For details on this proof, refer to the Supplementary Materials. The converse statement is also true, providing a further perspective on the connection to free fermions (28).


The holographic pentagon code

We will now apply our framework to the highly symmetric class of regular bulk tilings, first implementing the holographic error correcting code (HaPPY code) proposed in (12) and then exploiting the versatility of our framework to extend it toward more physical setups. The HaPPY code furnishes a mapping between additional (uncontracted) bulk degrees of freedom on each tensor and the boundary state, realized by a bulk tiling of pentagons. Each pentagon tile encodes one fault-tolerant logical qubit via the encoding isometry of the five-qubit code. This [[5,1,3]] quantum error–correcting code (29) saturates both the quantum Hamming bound (30, 31) and the singleton bound (31) and can be expressed as a stabilizer code (32).

We observe that fixing the bulk degrees of freedom to computational basis states gives rise to a matchgate tensor network, as the logical computational basis states of the holographic pentagon code can be viewed as ground states of a quadratic fermionic Hamiltonian. This can be seen directly by applying Eqs. 7 and 8 onto the stabilizers Sk of the underlying [[5,1,3]] code, thus expressing it in terms of Majorana operators γi and a total parity operator Ptot=(σz)5 asS1=σxσzσzσx12=iγ7γ2S2=12σxσzσzσx=iγ9γ4S3=σx12σxσzσz=iPtotγ6γ1S4=σzσx12σxσz=iPtotγ8γ3S5=σzσzσx12σx=iPtotγ10γ5(9)

As the corresponding stabilizer Hamiltonian is given by H=k=15Sk, we find a doubly degenerate ground state whose degeneracy is lifted by the parity operator Ptot. The resulting two states with parity eigenvalues ±1 correspond to the logical eigenstates 0¯ and 1¯, which are themselves ground states of purely quadratic Hamiltonians with different parity factors. Thus, both computational basis states are pure Gaussian, leading us to the conclusion that for fixed computational input in the bulk, the holographic pentagon code yields a matchgate tensor on the boundary (see Fig. 2). The explicit construction is given in the Supplementary Materials. Using the Schläfli symbol {p, q}, where p is the number of edges per polygon and q is the number of polygons around each corner, we can specify the hyperbolic geometry of the HaPPY model as a regular {5, 4} tiling.

Fig. 2 HaPPY/matchgate equivalence.

The holographic pentagon code of the HaPPY model for fixed computational bulk input (left) is equal to a matchgate tensor network on a hyperbolic pentagon tiling (right).

We find that the correlation structure of this model is best captured in the Majorana picture. Explicitly, consider the pentagon tiling of (12) with all bulk inputs set to the positive-parity eigenvector 0¯. The entries of the Majorana covariance matrix Γj,k=i2ψ[γj,γk]ψ resulting from successive contraction steps are shown in Fig. 3 (A to C). As we can see, both the individual pentagon state and the larger contracted states are characterized by a nonlocal pairing of Majorana fermions. The contractions effectively connect Majorana pairs from each pentagon to a larger chain, so the pairs on the boundary of the contracted network can be seen as end points of a discretized “geodesic” spanning the bulk. While this discontinuous correlation pattern of Γj,k makes the computation of CFT observables difficult, we can estimate the average correlation falloff by counting the relative frequency n(d) of Majorana pairs at distance d = ∣ jk∣ over which they connect points on the boundary. According to the results shown in Fig. 4A, correlation falloff follows a power law n(d) ∝ d−1, as expected of a CFT. Furthermore, we compute the entanglement entropy SA of a subsystem A of size l averaged over all boundary positions, defined asEl(S)=k=1LS[k,k+l](10)The result, shown in Fig. 4B, closely follows the Calabrese-Cardy formula for periodic 1 + 1–dimensional CFTs, given by (33, 34)SA=c3log(LπεsinπlL)c3loglε+O((l/L)2)(11)with a numerical fit yielding c ≈ 4.2 and ε ≈ 1.1 for a cutoff at L = 2605 boundary sites.

Fig. 3 Boundary state correlations.

(A to C) Majorana covariance matrix Γ with color-coded entries for a boundary state of a hyperbolic {5,4} tiling of the HaPPY code with fixed 0¯ input on each tile. Boundary consists of 2L = 10, 40, and 50 Majorana sites, respectively. (D to F) Field correlation matrix 〈ψjψk − ψkψj〉/2 = (Γ2j,2k−1 + Γ2j−1,2k)/4 for boundary states of the {3,6}, {3,7}, and mMERA tiling at criticality with L = 63, 69, and 64 boundary sites, respectively. Matrix entries are normalized to the same color scale. The tiling corresponding to each correlation matrix in (A) to (F) is shown in the lower left corner.

Fig. 4 Critical correlations and entanglement scaling.

(A and B) Boundary state properties of the HaPPY code at 2605 boundary sites. (A) shows average correlations at boundary distance d, computed as the relative frequency n of Majorana pairs. Dashed gray line shows an n(d)~1/d numerical fit. (B) shows the scaling of average entanglement entropy El(S) with subsystem size l. Dashed gray line shows numerical fit using (11). (C) El(S) for regular tilings at the critical values a = 0.580 for a {3,6} tiling (blue) and at a = 0.609 for the {3,7} tiling (yellow) with 348 boundary sites each. The dashed gray line shows the exact c = 1/2 CFT solution.

The peculiar pairwise correlation of boundary Majorana modes, suggesting a connection to Majorana dimer models (35), is more deeply explored in a separate publication (36). However, as the correlation structure breaks the translation and scale invariance expected of CFT ground states, we now consider regular tilings with generic matchgate input.

Regular triangulations

As the boundary states of triangular tilings are necessarily Gaussian (15), we can study their properties comprehensively using matchgate tensors. The simplest such tilings are regular and isotropic, i.e., with each local tensor specified by the same antisymmetric 3 × 3–generating matrix A. Isotropy constrains its components to one parameter a = A1,2 = A1,3 = A2,3. The bulk topology follows from our choice of tiling. For triangular tilings (p = 3), setting q = 6 produces a flat tiling, whereas q > 6 leads to a hyperbolic one (see Fig. 1, A and B). Triangular tilings with q < 6 produce closed polyhedra that are positively curved and lack the notion of an asymptotic boundary. As a convention, we choose the local orientation of the triangles so that the generating matrix for the contracted boundary state satisfies Ai,j>0 for i > j, corresponding to antiperiodic boundary conditions: Covariance matrix entries Γi,j acquire a sign flip when cyclic permutions push either index i or j over the boundary, as relative ordering is reversed.

We now consider the boundary states of {3, k} bulk tilings. The falloff of correlations along the boundary generally depends on k, i.e., the bulk curvature, as shown in Fig. 5 (A and B) for the a = 0.25 case. While correlations between the boundary Majorana fermions of a flat bulk fall off exponentially, a hyperbolic bulk produces a polynomial decay (up to finite-size effects at large distances and rounding errors at very small correlations). In the hyperbolic case, geodesics between boundary points scale logarithmically in boundary distance, so the falloff is still exponential in bulk distance, as we would expect in AdS/CFT (37).

Fig. 5 Boundary correlations for regular tilings.

(A and B) Mean value of Majorana covariance Ed(Γ)=k=1LΓk,k+d/L (with Γi,L+j = Γi,j) at boundary distance d. For {3,6} tiling with 150 boundary Majorana fermions (left) and {3,7} tiling with 348 (right). a = 0.25 in both cases. (C and D) Dependence of correlation falloff on a for {3,6} tiling with falloff ~ed (left) and {3,7} tiling with ∝dp (right).

Restricting ourselves to the 0 < a < 1 region, we explore how quickly correlations decay in both settings. At a = 0 and a = 1, the boundary Majorana fermions only have neighboring pair correlations, either pairing within each edge (a = 0) or across the corners (a = 1). Thus, correlation decay becomes infinite in the limits a → 0 and a → 1, independent of bulk geometry. We use numerical fits to study the remaining region 0 < a < 1 (see Fig. 5, C and D). For a hyperbolic bulk geometry, the power law is generic with the slowest decay at a ≈ 0.61, where we see a ∝d−1 falloff over distance d. The exponential decay ∝ed generally produced by a flat bulk geometry, however, slows down to a power law (with correlation length λ diverging) around a ≈ 0.58, where correlations again decay as ∝d−1. At their critical values, the boundary states of both bulk geometries have the same average properties.

Up to finite-size effects, this critical boundary theory turns out to be the Ising CFT, as we confirm by computing a range of critical properties from the covariance matrix, shown in Table 1. The entanglement entropy scaling, shown in Fig. 4C, again matches the expected form (11) irrespective of the choice of tiling. The Ising CFT state that we observe at the critical value of a is the ground state of the HamiltonianH=i(k=1N1γkγk+1+γ1γN)(12)where the sign of the boundary term γ1γN signifies antiperiodic boundary conditions. Triangular tilings also incorporate more generic models: By associating each edge with a bond dimension χ > 2, it is possible to produce boundary theories with central charges c larger than 1/2. In the simplest case, we choose a generating matrix that only couples between sets of fermionic modes, resulting in a boundary theory that consists of multiple copies of the Ising CFT and a corresponding central charge that is a multiple of 1/2 (note that this construction is only possible for {n, k} tilings with even k). Furthermore, by changing the tensor content in a central region of the network, a mass gap can be introduced, highlighting how radii in a hyperbolic bulk correspond to a renormalization scale on the boundary. Details are provided in the Supplementary Materials.

Table 1 Table of conformal data for the regular {3, 6} and {3, 7} bulk tilings as well as the mMERA, compared to the exact results and the wavelet MERA (16).

Listed are the ground-state energy density ϵ0, central charge c, scaling dimensions Δϕ of the fields ϕ=ψ,ψ¯,ϵ,σ, and the structure constant Cσ,σ,ϵ. The nonscaling of the identity 1 is discussed in the Supplementary Materials.

View this table:

Translation invariance and MERA

The regular bulk tilings considered so far have a set of discrete symmetries. When choosing identical tensors on each polygon, the boundary states necessarily inherit these symmetries, breaking translation invariance. To recover it, we consider a tiling with the same geometry as the MERA network. As we restrict ourselves to real generating matrices for the three- and four-leg matchgate tensors in this geometry, our model is not a unitary circuit but a model of Euclidean entanglement renormalization resembling imaginary time evolution, extending ideas from (19, 20). This may provide a more realistic representation of the causal structure of an AdS time slice than the standard MERA. Accordingly, the tensors of our matchgate MERA (mMERA) do not correspond to the usual (norm-preserving) isometries and disentanglers. We can still produce almost perfectly translation-invariant boundary states (Fig. 3F) while optimizing over only three parameters and recover the expected CFT properties (Table 1). In particular, at bond dimension χ = 2, the ground-state energy has a relative error of only 0.02% compared to the exact solution. Note that the optimization process only takes a few minutes on a desktop computer for a network with hundreds of tensors. We also find that the χ = 2 mMERA has a symmetry that allows us to write its four-leg tensors as contractions of simpler three-leg tensors (see Fig. 1C), yielding a nonregular triangular tiling. An interesting question to pursue is whether alternating or quasiperiodic tilings with a larger parameter space than regular tilings can also produce translation-invariant states.


In this work, we have studied bulk-boundary correspondences in fermionic Gaussian tensor networks, introducing a versatile framework and a highly efficient contraction method based on matchgate tensors (14, 15) for a wide class of flat and hyperbolic bulk tilings. We showed that our framework includes the holographic pentagon code built from five-qubit stabilizer states for fixed bulk inputs. Its boundary states correspond to a nonlocal bulk pairing of Majorana fermions, opening an avenue to studying the state properties of this holographic model at large sizes. We explicitly computed two-point correlators and entanglement entropies, which were found to exhibit critical scaling. Beyond known models, we showed that critical and gapped Gaussian boundary states can be realized by various bulk tilings. In particular, the average scaling properties of the c=12 Ising CFT (and multiples thereof) can be reproduced using regular one-parameter bulk triangulations with both flat and hyperbolic curvature. This is particularly unexpected for the flat case where boundary theories are typically gapped and raises the question whether this appearance of criticality is retained in strongly interacting models as well. The appearance of equivalent boundary CFT states for flat and hyperbolic bulks resembles the effect of local Weyl transformations in Euclidean path integrals (38). Our reproduction of conformal properties from an isotropic tensor network with only a single parameter further suggests that isotropy is a powerful symmetry for numerical CFT computations. Furthermore, we constructed the mMERA, a Euclidean matchgate tensor network based on the MERA geometry. Beyond the results achievable with regular triangulations, this tiling, which can also be expressed as a triangulation, recovers the Ising CFT with translation invariance while requiring only three free parameters and little computational cost. Within the Gaussian setting, further studies could focus on positively curved bulks, higher-dimensional models, and random tensors. Beyond Gaussianity, one could also explore interacting fermionic tensor networks (3943) by a weak-coupling expansion or under locally restricted interactions. Both of these possible extensions to our framework would require computations scaling only polynomially in the system size, thus still avoiding the prohibitive computational effort of general methods for exact tensor contraction.


Supplementary material for this article is available at

Section S1. Tensor contractions in the Grassmann formalism

Section S2. Matchgates and fermionic Gaussian states

Section S3. Conversion of generating matrices to covariance matrices

Section S4. Contraction rules for generating matrices

Section S5. Explicit generating matrices and numerical results

Fig. S1. Combining tiles of matchgates.

Fig. S2. Tile orientations under contraction.

Fig. S3. Constructing the mMERA.

Fig. S4. Energy convergence of the mMERA.

Fig. S5. Determining scaling dimensions of flat tilings.

Fig. S6. Determining scaling dimensions of hyperbolic tilings.

Fig. S7. Determining scaling dimensions of mMERA.

Fig. S8. Determining structure constants.

Fig. S9. Correlations and entanglement with IR cutoff.

Fig. S10. Construction of triangle states with bond dimension χ = 2, 4, 8.

Table S1. Values of the critical generating matrix parameter a for different {3, k} triangular tilings and ultraviolet cutoffs.

Table S2. Exact conformal scaling dimension of various (quasi-)primary fields ϕ of the Ising CFT.

References (4451)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We would like to thank A. Kubica for pointing out that [[5,1,3]] code states (29) are ground states of quadratic Majorana Hamiltonians. We also thank X.-L. Qi, T. Takayanagi, and P. Caputa for helpful discussions. Funding: We thank the ERC (TAQ), the DFG (CRC 183, EI 519/7-1, EI 519/14-1, and EI 519/15-1), the John Templeton Foundation, the EC (PASQuanS), the Studienstiftung, and the Alexander von Humboldt Foundation for support. Author contributions: F.P. and J.E. conceived the project and provided guidance, while A.J. and M.G. performed analytical calculations. A.J. provided the numerical computations and figure design. All authors discussed the results and wrote the manuscript. Competing interests: The authors declare that they have no competing interests. F.P. is currently employed by Psi Quantum Inc. but provided his main contribution to this work while being employed at Freie Universität Berlin. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. All data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article